Abstract Background The tumor suppressor gene E-cadherin gene is frequently silenced in chronic lymphocytic leukemia (CLL) cells and results in wnt-pathway activation. We analyzed the role of histone epigenetic modifications in E-cadherin gene silencing. Methods CLL specimens were treated with histone deacetylase inhibitor (HDACi) MS-275 and analyzed for E-cadherin expression with western blot and RT-PCR analysis. The downstream effects of HDACi treated leukemic cells were studied by analyzing the effect on wnt-pathway signaling. HDACi induced alterations in E-cadherin splicing were investigated by transcript specific real time PCR analysis. Results Treatment of CLL specimens with histone deacetylase inhibitors (HDACi) treatment resulted in an increase of the E-cadherin RNA transcript (5 to 119 fold increase, n=10) in eight out of ten CLL specimens indicating that this gene is down regulated by histone hypoacetylation in a majority of CLL specimens. The E-cadherin re-expression in CLL specimens was noted by western blot analysis as well. Besides epigenetic silencing another mechanism of E-cadherin inactivation is aberrant exon 11 splicing resulting in an alternatively spliced transcript that lacks exon 11 and is degraded by the non-sense mediated decay (NMD) pathway. Our chromatin immunoprecipitation experiments show that HDACi increased the acetylation of histones H3 and H4 in the E-cadherin promoter region. This also affected the E-cadherin exon 11 splicing pattern as HDACi treated CLL specimens preferentially expressed the correctly spliced transcript and not the exon 11 skipped aberrant transcript. The re-expressed E- cadherin binds to β-catenin with inhibition of the active wnt-beta-catenin pathway in these cells. This resulted in a down regulation of two wnt target genes, LEF and cyclinD1 and the wnt pathway reporter. Conclusion The E-cadherin gene is epigenetically modified and hypoacetylated in CLL leukemic cells. Treatment of CLL specimens with HDACi MS-275 activates transcription from this silent gene with expression of more correctly spliced E-cadherin transcripts as compared to the aberrant exon11 skipped transcripts that in turn inhibits the wnt signaling pathway. The data highlights the role of epigenetic modifications in altering gene splicing patterns.

Background

To determine differentially expressed and spliced RNA transcripts in chronic lymphocytic leukemia specimens a high throughput RNA-sequencing (HTS RNA-seq) analysis was performed.

Methods

Ten CLL specimens and five normal peripheral blood CD19+ B cells were analyzed by HTS RNA-seq. The library preparation was performed with Illumina TrueSeq RNA kit and analyzed by Illumina HiSeq 2000 sequencing system.

Results

An average of 48.5 million reads for B cells, and 50.6 million reads for CLL specimens were obtained with 10396 and 10448 assembled transcripts for normal B cells and primary CLL specimens respectively. With the Cuffdiff analysis, 2091 differentially expressed genes (DEG) between B cells and CLL specimens based on FPKM (fragments per kilobase of transcript per million reads and false discovery rate, FDR q < 0.05, fold change >2) were identified. Expression of selected DEGs (n = 32) with up regulated and down regulated expression in CLL from RNA-seq data were also analyzed by qRT-PCR in a test cohort of CLL specimens. Even though there was a variation in fold expression of DEG genes between RNA-seq and qRT-PCR; more than 90 % of analyzed genes were validated by qRT-PCR analysis. Analysis of RNA-seq data for splicing alterations in CLL and B cells was performed by Multivariate Analysis of Transcript Splicing (MATS analysis). Skipped exon was the most frequent splicing alteration in CLL specimens with 128 significant events (P-value <0.05, minimum inclusion level difference >0.1).

Conclusion

The RNA-seq analysis of CLL specimens identifies novel DEG and alternatively spliced genes that are potential prognostic markers and therapeutic targets. High level of validation by qRT-PCR for a number of DEG genes supports the accuracy of this analysis. Global comparison of transcriptomes of B cells, IGVH non-mutated CLL (U-CLL) and mutated CLL specimens (M-CLL) with multidimensional scaling analysis was able to segregate CLL and B cell transcriptomes but the M-CLL and U-CLL transcriptomes were indistinguishable. The analysis of HTS RNA-seq data to identify alternative splicing events and other genetic abnormalities specific to CLL is an added advantage of RNA-seq that is not feasible with other genome wide analysis.

Cellular heterogeneity is an integral part of cancer development and progression. Progression can be associated with emergence of cells that exhibit high phenotypic plasticity (including "de-differentiation" to primitive developmental states), and aggressive behavioral properties (including high tumorigenic potentials). We observed that many biomarkers that are used to identify Cancer Stem Cells (CSC) can label cell subsets in an advanced clinical stage of lung cancer (malignant pleural effusions, or MPE). Thus, CSC-biomarkers may be useful for live sorting functionally distinct cell subsets from individual tumors, which may enable investigators to hone in on the molecular basis for functional heterogeneity. We demonstrate that the CD44(hi) (CD44-high) cancer cell subsets display higher clonal, colony forming potential than CD44(lo) cells (n=3) and are also tumorigenic (n=2/2) when transplanted in mouse xenograft model. The CD44(hi) subsets express different levels of embryonal (de-differentiation) markers or chromatin regulators. In archived lung cancer tissues, ALDH markers co-localize more with CD44 in squamous cell carcinoma (n=5/7) than Adeno Carcinoma (n=1/12). MPE cancer cells and a lung cancer cell line (NCI-H-2122) exhibit chromosomal abnormalities and 1p36 deletion (n=3/3). Since miR-34a maps to the 1p36 deletion site, low miR-34a expression levels were detected in these cells. The colony forming efficiency of CD44(hi) cells, characteristic property of CSC, can be inhibited by mir-34a replacement in these samples. In addition the highly tumorigenic CD44(hi) cells are enriched for cells in the G2 phase of cell cycle.