Skip to main content

Advertisement

Gene expression and splicing alterations analyzed by high throughput RNA sequencing of chronic lymphocytic leukemia specimens

BMC Cancer volume 15, Article number: 714 (2015) Cite this article

Article metrics

Abstract

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.

Open Peer Review reports

Background

Chronic lymphocytic leukemia (CLL) is a common leukemia characterized by accumulation of B cells in the blood, marrow and lymphatic tissues. The clinical course is highly variable with biological and genetic heterogeneity in leukemic specimens. A number of genetic alterations have been correlated with prognosis [1–5]; however, the ability to prognosticate outcomes and tailor treatment based on genetic alterations is still limited. To identify genetic alterations in CLL, a number of different methods have been employed including cytogenetic studies [6], and array comparative genomic hybridization CGH [7, 8] and recently whole exome sequencing [9]. The whole exome sequencing of CLL specimens has also resulted in the identification of novel recurring mutations in the MYD88, NOTCH1, KLH6 and SF3B1 genes [10].

To study the complete transcriptome of cells, microarrays have been extensively used, and these studies have identified a number of differentially expressed genes [11–14]. Microarray techniques are, however, subject to a number of limitations including, cross hybridization of transcripts, limitation in coverage, inability to resolve novel transcripts and a falsely higher estimation of low abundance transcripts [15–18]. With the development of massive parallel RNA sequencing (RNA-seq) technology, there have been a growing number of genome-wide studies that have analyzed the complete transcriptome cells in different malignancies [18–22] and non-malignant diseases [23, 24]. Besides analyzing the expression level of genes the RNA-seq technology has the added advantage of analyzing expression at the exon level and provides detailed information about alternative splicing variations, novel transcripts, fusion genes, differential transcription start sites and genomic mutations [25, 26]. As all the RNA transcripts are being directly sequenced, this technology is ideally suited to study altered splicing pattern which is especially relevant in cancer cells as they are known to express unique RNA isoforms with varied biological effects [27, 28].

In this study, we performed RNA-seq analysis on CLL specimens and normal peripheral blood B cells to determine transcriptome differences and splicing variations. The data obtained from the RNA-seq analysis was validated by real time PCR on the RNA-seq cohort and a test cohort of specimens. Besides expression analysis a number of novel differentially spliced genes were also identified and analyzed. These findings will facilitate the identification of novel prognostic markers, therapeutic targets and signaling pathways in CLL.

Methods

Sample isolation and characterization

Primary CLL specimens analyzed in this study were obtained from untreated CLL patients after appropriate human subject approval. The human subject study was approved by the ethics committee of the West Los Angeles VA Medical Center and an informed written consent was obtained from all patients. A peripheral blood draw was performed, and peripheral blood lymphocytes (PBLs) were isolated by ficoll gradient. In all the CLL specimens, more than 90 % of isolated cells were CD19+ by flow cytometry analysis. Total RNA from isolated B cells (five different normal donors, caucasian males) was purchased from ALLCELLS (Alameda, CA). IGVH mutation (Immunoglobulin variable region heavy chain) analysis was performed on the CLL specimens with multiplexed PCR reactions to assess clonality as previously described [29]. Percentage of CLL cells expressing CD38 marker and Zap-70 (intracellular staining) was determined by flow cytometry and specimens with more than 20 % cells expressing Zap-70 were defined as Zap-70 positive. CLL specimens in a separate test cohort (n = 47) were from all clinical stages, chemotherapy naïve, and with more than 90 % CD19+ cells.

RNA-seq and library preparation

For library preparation, the Illumina TruSeq RNA sample Prep Kit v2 (San Diego, CA) was used according to the manufacturer’s protocol. Briefly, 1 μg of total mRNAs from five normal B and ten CLL cells was poly-A purified, fragmented, and first-strand cDNA reverse transcribed using random primers. Following second-strand cDNA synthesis, end repair, addition of a single A base, adaptor ligation, and PCR amplification, the enriched cDNA libraries were sequenced using the Illumina HiSeq 2000 at the UCLA Broad Stem Cell Research Center High Throughput Sequencing Core. The RNA sequencing data is deposited at GEO website, accession number GSE70830.

Primary processing and mapping of RNA-seq reads

50 bp single-end RNA-seq reads were obtained from Illumina HiSeq 2000. Sequence files were generated in FASTQ format (sequence read plus quality information in Phred format). RNA-seq data were analyzed using the UCLA Galaxy server (galaxy.hoffman2.idre.ucla.edu). The quality score of RNA-seq reads was obtained by using the FastQC and the mean quality of each base pair in the samples was 28, indicating a good-quality call in the 50 bp reads [30]. Reads were then processed and aligned to the UCSC H. sapiens reference genome (build hg19) using TopHat v1.3.3 [31–33].

Assembly of transcripts and differential expression

The aligned read BAM files were assembled into transcripts, their abundance estimated and tests for differential expression processed by Cufflinks v2.0.1 [33]. Cufflinks uses the normalized RNA-seq fragment counts to measure the relative abundances of transcripts. The unit of measurement is fragments per Kilobase of exon per Million fragments mapped (FPKM). Confidence intervals for FPKM estimates were calculated using a Bayesian inference method. After assembly with Cufflinks, the output files were sent to Cuffmerge along with a reference annotation file. To normalize multiple samples for differential expression analysis, we utilized a “geometric” method as described in Anders and Huber [34]. For cross-replicate dispersion estimation, a “pooled” method was used in which each replicated condition is used to build a model, and then these models are averaged to provide a single global model for all conditions in the experiment. The expression testing was done at the level of transcripts and genes and pairwise comparisons of expression between normal and CLL samples. Only the comparisons with “q-value” less than 0.05 and expression fold change greater than two fold in the Cuffdiff output were regarded as showing significant differential expression. Downstream analysis for Cuffdiff output was done using CummeRbund [34].

RT-PCR validation of RNA-seq results

The differentially expressed genes were validated by Quantitative Real-time Polymerase Chain Reaction (qRT-PCR) using a StepOnePlus™ Real-Time PCR System (Life technologies). cDNA templates from five normal B cell and ten CLL cells were analyzed for expression of DSP, TRIB2, DUSP1, FOS, JUN, SELPLG, AMICA, MMP9, TYROBP, and LEF1 with taqman probes obtained from Applied Biosystems. The probes selected for these genes provide the best coverage so that the majority of transcripts of the gene are quantified (further information is available on request). To analyze the IGVH subgroups, expression of three genes T, TFEC and IGLL5 was also determined with Taqman probes. Expression of a number of reference genes (Actin, Ribosomal protein large PO, phosphoglycerate kinase, Hypoxanthine phoshoribosyl transferase and Transferrin receptor) was tested for expression in CLL and B cells, and actin was selected as the standard reference gene and the data was analyzed by the method of Pfaffll [35].

Functional annotation of differentially expressed genes

The differentially expressed gene lists were submitted to Ingenuity Pathway Analysis (IPA, Ingenuity Systems). The functional annotation identifies the biological functions that are most significant to the data set. A Fisher’s exact test was used to calculate a p-value determining the probability that the association between the genes in the dataset and the functional annotation is explained by chance alone.

Alternative splicing analysis with MATS

The RNA-Seq data of B cells and CLL specimens was analyzed for splicing alterations. To identify such events, MATS 3.0.8 (Multivariate Analysis of Transcript Splicing, ref [36]) was used to determine junctional reads within ENSEMBL human gene annotations. This software implements a Bayesian approach that detects differential AS (alternative splicing) under two conditions by examining whether the difference in the exon-inclusion levels between two samples exceeds a given user-defined threshold. To identify these events, we used the following criteria, 1. Splicing events were labeled significant if the sum of the reads supporting a specific event exceeded 10 reads, 2. P-value was <0.05, and 3. Minimum inclusion level difference as determined by MATS was >0.1 (10 % difference). To validate the splicing alterations RT-PCR analysis was performed by designing primers in the neighboring exons (primer sequences available on request).

Results

Analysis of RNA-seq data

Five normal CD19+ B Cell RNA from different donors (B1 to B5), six IGVH un-mutated primary CLL specimens (CLL6, CLL9, CLL25, CLL28, CLL40, and CLL44) and four IGVH mutated CLL specimens (CLL26, CLL32, CLL37, and CLL39) were subjected to HTS-RNA single-end RNA sequencing (Table 1). The total WBC counts for unmutated IGVH (U-CLL) were higher than mutated IGVH (M-CLL) specimens (Table 1) and the U-CLL specimens were noted to have a higher percentage of leukemic cells expressing CD38 and Zap-70 as described before [4, 5]. The total number of raw reads in B cells (n = 5) and CLL (n = 10) specimens ranged from 31 to 85 million reads, and 37 to 101 million for normal, CLL, respectively (Fig. 1, Additional file 1). To assess the quality of mapping reads to the reference genome hg19, some key metrics were extracted from the TopHat output, and analyzed using the RNA-seq quality control package RseQC [37]. The majority of reads (between 65.5 % and 79.6 %) are uniquely mapped to the reference genome sequences across all samples (Additional file 1). The mean mapping percentage for normal B cells and CLL specimens is 78.3 % and 74.4 % and 5.8 % to 8.8 % of the reads mapped to the known splice junctions respectively.

Table 1 Clinical characteristics of CLL patients and RNA sequencing read count data
Full size table
Fig. 1
figure1

Distribution of sequencing reads in normal B cells and CLL specimens. a The bar diagram represents distribution of uniquely mapped reads to human genome UCSC_hg19 (GRCh37). Each bar depicts the percentage of reads from individual samples (five normal B cell and ten CLL specimens) mapped to coding sequence exon (CDS_exon), 5’ and 3’ untranslated regions (5’ and 3’UTR_Exons), introns and intergenic regions. b Pie charts represent the average percentage of sequencing reads from five normal B cell (left) and ten CLL specimens (right) that map to the above mentioned regions

Full size image

To further examine the read distribution, the uniquely mapped reads were assigned to: exon coding sequence (CDS_Exons), 5’ and 3’ untranslated regions (UTR_Exons), introns and intergenic regions. In Fig. 1a, the distribution of mapped reads is shown across the samples. 41 % to 52 % of reads mapped to exon coding sequence, 2.9 % to 3.8 % mapped to 5’UTR while 18 % to 25 % mapped to 3’UTR. The introns and intergenic regions account for about 15–30 % and 5–9 %, respectively (data for all specimens is in Additional file 1). To compare if there is a difference in read distribution between normal B cell and CLL, mapping data from Fig. 1a was averaged and plotted as a pie chart in Fig. 1b. The exonic reads (CDS_Exons) were higher in CLL specimens as compared to B cells while intronic reads were higher in the B cell specimens (Fig. 1b), 30 % vs. 16 % for normal B cells and CLL specimens. The high number of reads mapping to introns have been reported in other RNA-seq analysis [38] and could be due to genomic DNA contamination, sequencing of pre-mRNA, novel exons, or nascent transcription and co-transcriptional splicing as described in Ameur et al [39].

Defining the transcriptomic profiles of normal B cell, and CLL specimens

To examine the transcriptome profile of normal B cells and primary CLL specimens, transcripts were assembled and their expression values calculated using Cufflinks. Pair-wise comparisons of transcriptomic profiles of normal B cells, CLL specimens as well as disease-subtype as determined by IGVH mutational status (U-CLL, un-mutated IGVH and M-CLL, mutated IGVH), were performed. The transcript abundance was calculated by estimating the fragment per kilobase of exon per million mapped fragments (FPKM). The numbers of assembled transcripts for normal B cell, U-CLL and M-CLL were 10396, 10494, and 10402 and the genes identified for the three sample groups were 10081, 10111, and 10068, respectively (Additional file 2A). Overall, the number of transcripts and genes found in three groups are very similar indicating a uniform sequencing depth in the various groups.

To determine significant differences in the transcriptomic profiles in the three sample groups (B, U-CLL and M-CLL), pair wise scatter plots matrix was generated by CummeRbund [34]. This analysis compares and correlates the FPKM profile of all expressed genes in all three sample groups, and it also shows the density distribution of FPKM for genes expressed. In Additional file 2B, the density plot reveals that the FPKM distributions among three sample groups are similar, and the FPKM of all expressed genes ranged from 0.003 to 3000 (log10FPKM -2.5 to 3.5), with the majority of the genes expressing FPKM range from 1 to 100 (log10FPKM 0 to 2.5). The global profiles of U- and M-CLL show fewer dispersion as compared to plots where normal B cell data is compared to the CLL specimens indicating similar transcriptome profiles of U- and M-CLL specimens.

Analysis of differentially expressed genes

To determine the differentially expressed genes (DEG) between normal B cells and CLL specimens a Cuffdiff analysis was performed. After filtering differential expressed genes with FDR-adjusted (FDR false discovery rate) q value < 0.05 and fold change > 2, there were 2091 DEG genes between CLL specimens and normal B cells (Fig. 2a). Among these genes, 1231 were up-regulated in CLL and 860 genes were down-regulated (complete gene list in Additional file 3), and the top twenty genes in each group are shown in Table 2. The data was also analyzed by segregating CLL specimens based on their IGVH status and comparing them with normal B cells separately. With this analysis 2425 and 1960 DEG genes were identified in U- and M-CLL specimens respectively. Among these genes, 1332 and 1132 were up-regulated and 1093 and 828 were down-regulated in U-CLL and M-CLL (Fig. 2a). In order to find out if there are overlapping genes that are differentially expressed in both U-CLL and M-CLL samples, the gene lists from normal B cells vs. CLLs, normal B cells vs. U-CLL and normal B cells vs. M-CLL were compared to generate a Venn diagram (Fig. 2b). A high number (1382 genes out of 2091) of differentially expressed genes between normal B cells and CLLs were common to the U-CLL and M-CLL specimens, indicating that this subgroup includes a common set of differentially expressed genes.

Fig. 2
figure2

Transcriptomic expression profiles and validation. a The number of statistically significant differentially (Up and Down regulated) expressed genes identified from Cuffdiff analysis in various groups relative to B cells are shown in a table format. The differentially expressed genes (DEG, FDR-adjusted q-value < 0.05, Fold change > 2) in all CLL specimens (n = 10), U-CLL (n = 6) and M-CLL (n = 4) was compared to normal B cells. b Venn diagram illustrates the overlapped DEG between the three groups in panel A

Full size image
Table 2 Top twenty up (positive fold change) and down-regulated (negative fold change) genes in CLL versus Normal B cells
Full size table

To validate the RNA-seq data, a number of differentially expressed genes with potential biological relevance to CLL were selected from this analysis, and their FPKM data was compared to the expression level by real time RT- PCR (qRT-PCR). In an initial experiment the expression level of a number of reference genes in normal B cells and CLL specimens was determined to identify the appropriate reference gene. Actin, Ribosomal protein large PO, Phosphoglycerate kinase, Hypoxanthine phoshoribosyl transferase and Transferrin receptor expression was analyzed with Taqman probes and the expression of actin was the most abundant in all the CLL (n = 3) and B cell specimens (n = 3) and was selected as the standard reference gene. FOS (# 111), JUN (#152), DSP (desmoplakin #2), TRIB2 (Tribbles homolog 2, #66) and DUSP1 (dual specificity phosphatase 1 # 49) were selected from the set of genes that have a lower expression in CLL specimens than B cells (Table 2, Additional file 3, # represent the position of the gene based on the FPKM data, lower number corresponds to higher down regulation). AMICA1 (#48), MMP9 (#2), TYROBP (#49), SELPG (# 604), LEF1 (#64) were selected as candidate genes to compare the fold over-expression by the two methodologies (# represents fold over-expression relative to B cell based on FPKM data, smaller number indicates higher fold overexpression).

The RNA from the identical 10 CLL specimens and five normal B cells (control) was used to perform Taqman probe based qRT-PCR assays. Probes selected for expression analysis provide the best coverage for a particular gene. Figure 3a, b shows three sets of data for each gene expression (n = 10), expression based of FPKM values in RNA-seq cohort (n = 10), qRT-PCR data from the identical specimens (n = 10) as RNA-seq cohort (relative to actin) and qRT-PCR data of a test cohort (n = 47, relative to actin) of CLL specimens. Figure 3a bar diagram depicts average ∆∆cT values in the three cohorts with the table below showing the p values of data in Fig. 3a. For down regulated genes in CLL, only DSP and TRIB2 expression is significantly lower (p < 0.5) as compared to B cells, while in the set of up-regulated genes, the expression of SELPG, AMICA, TYROBP and LEF1 is significantly higher (p < 0.5) in the test cohort. MMP9 expression though significantly higher in the smaller RNA-seq cohort is not significantly higher in the test cohort.

Fig. 3
figure3

Validation analysis of selected differentially expressed genes. a qRT-PCR of selected genes on B cells (n = 5), CLL specimens, RNA seq cohort (n = 5) and CLL specimens, test cohort (n = 47). Data shown is the delta delta cT relative to actin. (Mean and standard deviation). Table below panel A shows the P-values of the qRT-PCR data for the comparison of B cells and CLL RNA seq cohort (n = 10), and B cells and CLL test cohort. (t-test). b Fold expression of selected genes in the larger CLL cohort (n = 47) based on qRT-PCR analysis. * PTPRK expression was not detected in normal B cells therefore fold change could not be calculated

Full size image

Figure 3b Table compares the fold expression obtained by these analysis, as an example in the case of DSP the difference in cT values between actin and DSP RNA is around 10 cycles while the expression in CLL cohorts is around 8 cycles lower, i.e. 256 fold fold down regulation of DSP expression in CLL specimens as compared to B cells. This lower DSP expression in CLL specimens is similar to the results obtained from FPKM analysis (179 fold lower expression in CLL). FOS and JUN expression based on RNA-seq FPKM data is 7.9 and 6.2 fold less than B cells while based on the qRT-PCR analysis their expression is 4.6 and 4.2 fold less than B cells. However in the test cohort (n = 47) the lower expression of JUN cannot be confirmed and for FOS the fold lower expression is less than the RNA-seq data (1.7 vs 7.9). Similar variability in expression is observed for MMP9 and AMICA1 expression as the fold expression vary 20 to 35 fold when analyzed by RNA-seq FPKM and qRT-PCR. The analysis shows that genes identified as differentially expressed by RNA-seq can be confirmed by qRT-PCR analysis, however the fold expression obtained by the two analysis are variable. Also confirmation with qRT-PCR in additional primary CLL specimens is required as there is significant variability of expression in leukemic cells.

Based on this analysis, additional DEG genes were selected to further compare the two methodologies for RNA expression (Table 3). FPKM and qRT-PCR fold expression levels were compared in the RNA-seq cohort and a test cohort of CLL specimens (n = 22). Nine downregulated genes from the RNA-seq data were randomly selected and their expression compared to qRT-PCR analysis. In the case of PTPRK, expression in normal B cells was not identified by qRT-PCR and therefore the RNA-seq data could not be validated. In the case of CCD69, the expression by RNA-seq and qRT-PCR is similar but this lower expression is not observed in our test cohort. Besides these two examples, qRT-PCR confirms a lower expression of these genes in CLL specimens as compared to normal B cells. Twelve genes with a range of over-expression were randomly selected from the list of over-expressed genes from the RNA-seq analysis (Additional file 3) and analyzed by qRT-PCR (Table 3). All the genes were found to be over-expressed based on qRT-PCR in the RNA-seq and test cohort however the expression was variable. Difference in fold expression was observed when the identical specimen was tested by both methodologies as well. Possible explanations for this discrepancy are the normalization of RNA-seq data and the use FPKM for calculation while qRT-PCR analysis is relative expression to a housekeeping gene and the Taqman probe may not provide coverage for all the transcript variants.

Table 3 Validation of twenty one differentially expressed genes in CLL. Data from RNA seq analysis (n = 10), qRT-PCR of identical specimens (n = 10) and qRT-PCR from a test cohort (n = 22) of CLL specimens
Full size table

Functional pathway analysis

The functional analysis tool was used to categorize genes that were differentially expressed in CLL specimens. Genes from Additional file 3 were analyzed by IPA analysis. The output of the functional annotation is shown in Additional file 4 and the list of genes in each pathway are in Additional file 5. The highest number of DEG genes are in the cell death and survival group correlating well with the unique biological characteristic of CLL, namely resistance to apoptosis. Other significant clustering of genes is observed in cellular movement, cellular development, cellular growth and proliferation and cancer pathways.

Comparison the CLL IGVH mutated and non-mutated transcriptomes

Based on the Cuffdiff analysis in Fig. 2a and b a number of genes are differentially expressed in the two CLL subsets, M and U-CLL. A total of 679 genes were more than 2 fold up or down regulated when the average FPKM data of all the genes was compared in the two subsets (Fig. 4a, Additional file 4). To determine whether global transcriptome analysis could segregate the CLL specimens based on IGVH status, a multi-dimensional scaling (MDS) plot (Fig. 4b) was constructed based on their complete transcriptomes. This analysis visualizes the level of similarity of individual samples within a group. MDS analysis was able to segregate the five normal B cells (B1-B5) as they cluster together away from the ten CLL specimens. The CLL specimens, U-CLL (closed boxes) and M-CLL (closed triangles) appear to be separate from each other but there is overlap of CLL specimens #25, #39, #37. Lack of clear separation of specimens on this plot indicates that based on the transcriptome data, M- and U-CLL specimens are not fully distinguishable.

Fig. 4
figure4

Transcriptomic comparison of IgVH mutated (M-CLL) and non-mutated (U-CLL). a Table with Cuffdiff data showing significant differentially expressed genes between M- and U-CLL specimens. b MDS plot (Multi-Dimensional Scaling) shows the clustering of the transcriptomic expression profiles of normal B cells (B1-B5), U-CLL and M-CLL samples (numbered as in Table 1). Axes in the MDS plot (M1 and M2) are arbitrary, and the values on the axes are distance units. c, d, e qRT-PCR data from the RNA-seq cohort of CLL specimens (n = 10) for three selected genes (T, IGLL5 and TFEC) (relative to actin, log scale). These panels show the scatter-plot qRT-PCR data in a separate cohort of CLL specimens and compare the expression of the three selected genes in M and U-CLL specimens. The dotted line separates the M- and U-CLL specimens

Full size image

In Additional file 6, the list of differentially expressed between the two groups (U-CLL and M-CLL) is shown that was obtained by dividing the mean FPKMs of the two sub-groups. From this list, we identified three genes for further analysis on an additional cohort of CLL specimens to determine if the expression of these genes is different in these two subgroups. The expression of IGLL5 (immunoglobulin lambda-like polypeptide 5, immunoglobulin lambda-like polypeptide) and T (brachyury homolog, embryonic nuclear transcription factor), was higher in the U-CLL group (top twenty most over-expressed genes in U-CLL as compared to M-CLL) and the expression of TFEC (transcription factor EC) was similarly higher in M-CLL group. The expression of these genes was determined by qRT-PCR in a separate cohort of 21 CLL specimens (relative to actin) and is shown in Fig. 4d, e and f scatter plots. The dotted line divides the U- and M-CLL specimens and the expression of these genes in the two sub-groups of CLL specimens is not significantly different when additional CLL specimens are analyzed. The expression of these three specific genes and the transcriptome as a whole for the U- and M-CLL specimens are similar.

RNA splicing alterations in CLL specimens

Besides accurately identifying the expression of the genes, the RNA–seq data is also useful in characterizing alternative splicing events. Splicing alterations in CLL specimens can alter the type of transcripts and thereby function of a large number of cellular proteins that may provide the cell with survival advantage. To define the splicing alterations in CLL specimens, the available RNA-seq data was analyzed by MATS (Material and Methods). Fig. 5a is a schematic of the various alternative splicing (AS) events that were analyzed and the number of events identified are listed in Fig. 5b. The analysis identifies AS events both in normal B cells and CLL specimens. Skipped exon (SE) is the most common splicing abnormality with 40974 events of which 128 events passing the threshold for significance. The complete lists of all the significant splice events in Fig. 5b table is in Additional file 7. The significant events in Fig. 5b table are divided into two columns, B and CLL which indicate whether the splicing event led to a higher inclusion of the exon in B or CLL specimens, e.g. 78 SE events resulted in a higher inclusion of exon in B cells and in 50 events, the inclusion of the exon was higher in CLL specimens.

Fig. 5
figure5

Alternative splicing events in B and CLL specimens. a Schematic showing alternative splicing (AS) events (from MATS analysis website). b Table with MATS analysis data with different AS events, total events and significant events are shown. B and CLL columns indicate the events out of all the significant events that had higher inclusion levels in either B or CLL specimens

Full size image

As the SE events, were by far the predominant events they were analyzed by RT-PCR analysis. Sixteen genes (listed in Additional file 8) were selected for initial analysis and primers were designed in the neighboring exons. To confirm DNA amplification of alternatively spliced exons, RT-PCR analysis was performed on B and CLL specimen (Additional file 8). Out of sixteen SE events in two genes there was no PCR amplification and in three genes only one DNA fragment corresponding to a single transcript was amplified (gels in Additional file 8). Low level of transcripts that are not amplified by the PCR is a likely reason that SE events could not be confirmed in these three genes. From the remaining eleven SE events TRIP11, TP53, MBNL2, ARGLU1, PER1, and PTPRC genes were randomly selected for further analysis and RT-PCR analysis was performed on B cell (n = 5) and CLL specimens (n = 9, Fig. 6). For each SE event Fig. 6 describes the exon of the gene that is alternatively spliced, expected base-pair size of the transcript with and without skipped exon along with average Inc. level (inclusion level, based on DNA band densitometry). TRIP11 (thyroid hormone receptor interactor 11 protein), tumor suppressor p53 (TP53), ARGLU1 (arginine and glutamate rich-1), Per1 (period 1) and PTPRC (CD45) demonstrate at least a two-fold difference in inclusion level of the SE exons in their transcripts. The analysis for MBNL2 (Muscleblind-like splicing regulator 2) did not show any inclusion level difference between normal B and CLL specimens.

Fig. 6
figure6

Validation of alternative splicing events. RT-PCR analysis of six AS events. For each gene, five B cell specimens and nine CLL specimen was analyzed. Expected bp (base pair) of the DNA fragments, with schematic of the skipped exon and mean Inc level (inclusion level, based on gel densitometry) are shown

Full size image

Discussion

Accurate transcriptome analysis is crucial for determining the expression of genes and thereby activity of signaling pathways that result in growth and survival of leukemic cells. The data from HTS RNA sequencing is an improvement over previous methodologies to effectively and efficiently evaluate the entire transcriptome. The RNA-seq data allows additional analysis of splicing alterations, transcriptional start sites, identification of novel signaling pathways and molecular categorization of specimens that is not feasible with prior genome analytic techniques. With improvement in HTS-sequencing technologies and reduction in the cost of sequencing it is now feasible to compare clinical and biological characteristics of CLL specimens with their global transcriptome profile.

In this study 20 % of genes are identified as differentially expressed (FDR q value < 0.05 and fold change > 2) in CLL specimens as compared to primary B cells. Recently, Ferreira et al has reported RNA-seq and transcriptome analysis of a large cohort of CLL specimens [40]. They report 1089 differentially expressed genes (DEG) between normal B cells and CLL specimens (FDR < 0.01 and median fold change of more than 3). This compares well to our analysis of 2091 DEG genes with a slightly less stringent FDR of < 0.05 and fold change more than 2. A number of DEG genes identified by this analysis were also reported by Ferreira et al [40], e.g. FOS, JUN, CYBRD1, GZMB, FMOD, CTLA-4, etc. In this study, data from the RNA-seq analysis was additionally also validated by qRT-PCR in a separate cohort of CLL specimens. Even though a similar expression trend in expression is observed when the genes from the RNA-seq analysis are validated by qRT-PCR, in some instances there can be wide variation based on FPKM and qRT-PCR. There are a number of reasons for the observed differences, e.g. the library preparation for RNA-seq analysis uses mRNA as starting material while total RNA was used for qRT-PCR. Another reason is that the RNA-seq analysis uses the FPKM method for normalization while actin (reference gene) was used as a control with qRT-PCR. Even though the Taqman probes with maximum coverage were used for the assay, but it is possible that some transcript variants were not analyzed by qRT-PCR as we observe that qRT-PCR under-estimates the level of expression as compared to the RNA-seq data in some genes.

Some of the DEG genes identified in this study have been reported earlier in microarray studies, e.g. MMP-9 and FMOD (fibromodulin) over-expression in CLL specimens has been described [41, 42]. MMP-9, matrix metallopeptidase 9, functions by degrading a number of matrix proteins such as type IV collagen, the major component of basement membranes. This gene was found to be highly expressed by CLL cells present in the bone marrow and lymph nodes, and contribute to B-CLL progression by facilitating cell migration and tissue invasion. FMOD, fibromodulin, is a member of a family of small interstitial proteoglycans and a component of the extracellular matrix that may also regulate TGF-beta activities by sequestering TGF-beta into the extracellular matrix [43]. SiRNA knock-down of this gene results in apoptosis of CLL, indicating its role in CLL survival [42]. SEPLG, TYROBP, LEF1 and AMICA1 genes were significantly up regulated in a number of CLL specimens. SELPLG (CD162) is a cell adhesion molecule that is the counter-receptor for selectins and plays a role in lymphocyte trafficking. High expression of SELPLG could potentially aid leukemic cells in trans-endothelial migration by interacting with the selectins on the endothelium cells [44]. TYROBP or Dap-12 is transmembrane protein that contains ITAM motifs (immunoreceptor tyrosine based activation motif) that are also present in the B-cell receptor (BCR) signaling components [45]. ITAM motifs are central to BCR signaling as a number of signaling molecules and adapter proteins assemble at these motifs. LEF1 (lymphoid enhancer-binding factor-1) gene encodes a transcription factor that participates in wnt signaling pathway that is active in CLL specimens. LEF1 is also involved in the transcriptional activation of Myc and CyclinD1, and both these genes are also up regulated in CLL leukemic cells [46]. AMICA1 expression was marginally higher in the larger cohort of CLL specimens (p = 0.032) and is a membrane protein that interacts with CXADR antigen expressed on epithelial and endothelial cells [47]. Table 3 lists additional DEG genes that were confirmed by qRT-PCR analysis. Pim1 kinase over-expression has been reported and this phosphorylates CXCR4 receptor that in turn mediates microenvironment signaling [48]. Similarly PDE4 transcripts in CLL specimens have been described [49] and Lck is associated with B-cell receptor signaling and blocking LCk function results in apoptosis [50].

40–50 % of the DEG genes demonstrate loss of expression in CLL specimens as compared to B cells. FOS and JUN down regulation in CLL specimens has been reported but in this study, the expression was not significantly lower when this was studied in a larger cohort of specimens [51–53]. One of the mechanisms of FOS down regulation is by its interaction with TCL1 oncogene that is a potential mechanism of resistance to apoptosis observed in CLL cells. DSP (desmoplakin) and TRIB2 expression were significantly lower in CLL specimens. DSP is a key component of the desmosomes that form intercellular junctions and loss of its expression is associated with more invasive behavior of cells [54]. It can also potentially function as a tumor suppressor gene by inhibiting the Wnt/β-catenin signaling pathway [55]. TRIB2 is a member of the Tribbles family of proteins that are similar to serine-threonine kinases but lack catalytic function. These proteins are highly conserved and modulate a number of signaling pathways [56]. DUSP1, is a phosphatase that controls cell proliferation and its expression was not significantly lower in the larger cohort of CLL specimens in this analysis [57] as compared to B cells. Additional genes that are downregulated or silenced in CLL specimens are listed in Table 3. UACA (Uveal autoantigen with coiled-coil domains and Ankyrin repeats) that regulates expression of an apoptotic regulator APAF1 [58]. JUP (junctional plakoglobin or gamma-catenin) associates with cytoplasmic domains of cadherins and has tumor and metastasis suppressor activity [59]. Based on their reported functions both these genes are potential tumor suppressor genes in this leukemia as well.

Based on the MDS (multi-dimensional scaling) analysis of the RNA-seq data, the normal B cells and CLL specimens could be segregated on a two-dimensional plot scaling plot. However, the transcriptome data does not clearly distinguish the U- and M-CLL transcriptomes as there is overlap on the scaling plot. The two sub-groups have important biological and clinical differences [60], but their transcriptomes are not significantly different. Expression analysis of selected genes (T, IGLL5 and TFEC, ref [61–63]) in the two IGVH sub-groups gave a similar result with no significant difference of expression. The study by Ferreira et al [40] reached an identical conclusion as their analysis could not detect significant transcriptome differences in these two groups. This has also been the observation of other groups that have performed microarray analysis of CLL specimens and have reached a similar conclusion [13, 14].

Alternative splicing events add another layer of complexity besides genes expression as they can alter the structure and function of cellular proteins. Skipped exons are the most common alternative spliced events in CLL specimens in this study and cancer cells in general [27, 28]. Splicing abnormalities have increasingly become more relevant in CLL with the identification of mutations in SF3B1, a splicing factor in a small subset of CLL patients [9] that confer poorer prognosis and can alter RNA splicing patterns. In our analysis, we focused on differentially exon skipping (SE) events as they were by far the most frequent events as compared to alternative 5’, alternative 3’, mutually exclusive exons, and retained intron. Confirmation of different inclusion levels in CD45 (PTPRC, a phosphotyrosine phosphatase), TP53, ARGLU1, PER1, and Trip11 genes by RT-PCR indicates that RNA-seq data can be analyzed for splicing alterations. Alternative exon usage in some of these genes such as PTPRC and p53 is well described in previous studies [64–68]. PTPRC is a member of the protein tyrosine phosphatase family that has a role in antigen receptor signaling, B cell development and may modulate signaling via integrins and cytokine receptors [66]. A number of studies have characterized expression of PTPRC (CD45) isoforms in CLL leukemic cells due to splicing in exons 4, 5 and 6 that alter the extra-cellular domain of the protein. It is however, not well understood whether exon skipping and expression of a particular isoform changes the function of this phosphatase. Period1 is a gene expressed in a circadian pattern with probable tumor suppressor function and its alternative splice forms though reported have not been characterized [69]. Alternative exon usage in less well characterized genes Trip11 and ARGLU1 was identified by this RNA-seq analysis and confirmed by RT-PCR. However their role in CLL biology is currently not clear and this will require additional studies to sequence novel transcripts in leukemic cells and to determine whether alternative exon usage alters the function of the expressed protein.

Conclusion

The main strength of RNA sequencing data is that besides providing expression analysis it can be further mined for a number of other genetic abnormalities, including splicing alterations, fusion transcripts, alternate transcription start sites, point mutations, novel transcripts, fusion genes etc. that will provide novel insights in this leukemia. As there is variability of expression in primary leukemic specimens and occasionally between RNA-seq data and qRT-PCR, further confirmation of RNA-seq data is required to obtain accurate information. Novel DEG and spliced transcripts were identified that potentially have biological significance in this leukemia and are valuable leads for discovery of novel biomarkers and therapeutic targets in this disease.

Abbreviations

AS:

Alternative splicing

CGH:

Comparative genomic hybridization

CLL:

Chronic lymphocytic leukemia

DEG:

Differentially expressed gene

FPKM:

Fragments per Kilobase of exon per million fragments mapped

IGVH:

Immunoglobulin variable region heavy chain

M-CLL:

IGVH mutated CLL

MATS:

Multivariate Analysis of Transcript Splicing

SE:

Skipped exon

U-CLL:

IGVH non-mutated CLL

References

  1. 1.

    Damle RN, Wasil T, Fais F, Ghiotto F, Valetto A, Allen SL, et al. Ig V gene mutation status and CD38 expression as novel prognostic indicators in chronic lymphocytic leukemia. Blood. 1999;94(6):1840–7.

  2. 2.

    Hamblin TJ, Davis Z, Gardiner A, Oscier DG, Stevenson FK. Unmutated Ig V(H) genes are associated with a more aggressive form of chronic lymphocytic leukemia. Blood. 1999;94(6):1848–54.

  3. 3.

    Ibrahim S, Keating M, Do KA, O'Brien S, Huh YO, Jilani I, et al. CD38 expression as an important prognostic factor in B-cell chronic lymphocytic leukemia. Blood. 2001;98(1):181–6.

  4. 4.

    Crespo M, Bosch F, Villamor N, Bellosillo B, Colomer D, Rozman M, et al. ZAP-70 expression as a surrogate for immunoglobulin-variable-region mutations in chronic lymphocytic leukemia. N Engl J Med. 2003;348(18):1764–75.

  5. 5.

    Chen L, Widhopf G, Huynh L, Rassenti L, Rai KR, Weiss A, et al. Expression of ZAP-70 is associated with increased B-cell receptor signaling in chronic lymphocytic leukemia. Blood. 2002;100(13):4609–14.

  6. 6.

    Dohner H, Stilgenbauer S, Dohner K, Bentz M, Lichter P. Chromosome aberrations in B-cell chronic lymphocytic leukemia: reassessment based on molecular cytogenetic analysis. J Mol Med. 1999;77(2):266–81.

  7. 7.

    Gunn SR, Mohammed MS, Gorre ME, Cotter PD, Kim J, Bahler DW, et al. Whole-genome scanning by array comparative genomic hybridization as a clinical tool for risk assessment in chronic lymphocytic leukemia. J Mol Diagn. 2008;10(5):442–51.

  8. 8.

    Higgins RA, Gunn SR, Robetorye RS. Clinical application of array-based comparative genomic hybridization for the identification of prognostically important genetic alterations in chronic lymphocytic leukemia. Mol Diagn Ther. 2008;12(5):271–80.

  9. 9.

    Quesada V, Conde L, Villamor N, Ordonez GR, Jares P, Bassaganyas L, et al. Exome sequencing identifies recurrent mutations of the splicing factor SF3B1 gene in chronic lymphocytic leukemia. Nat Genet. 2011;44(1):47–52.

  10. 10.

    Puente XS, Pinyol M, Quesada V, Conde L, Ordonez GR, Villamor N, et al. Whole-genome sequencing identifies recurrent mutations in chronic lymphocytic leukaemia. Nature. 2011;475(7354):101–5.

  11. 11.

    Klein U, Tu Y, Stolovitzky GA, Mattioli M, Cattoretti G, Husson H, et al. Gene expression profiling of B cell chronic lymphocytic leukemia reveals a homogeneous phenotype related to memory B cells. J Exp Med. 2001;194(11):1625–38.

  12. 12.

    Falt S, Merup M, Gahrton G, Lambert B, Wennborg A. Identification of progression markers in B-CLL by gene expression profiling. Exp Hematol. 2005;33(8):883–93.

  13. 13.

    Haslinger C, Schweifer N, Stilgenbauer S, Dohner H, Lichter P, Kraut N, et al. Microarray gene expression profiling of B-cell chronic lymphocytic leukemia subgroups defined by genomic aberrations and VH mutation status. J Clin Oncol. 2004;22(19):3937–49.

  14. 14.

    Rosenwald A, Alizadeh AA, Widhopf G, Simon R, Davis RE, Yu X, et al. Relation of gene expression phenotype to immunoglobulin mutation genotype in B cell chronic lymphocytic leukemia. J Exp Med. 2001;194(11):1639–47.

  15. 15.

    Pawitan Y, Michiels S, Koscielny S, Gusnanto A, Ploner A. False discovery rate, sensitivity and sample size for microarray studies. Bioinformatics. 2005;21(13):3017–24.

  16. 16.

    Marioni JC, Mason CE, Mane SM, Stephens M, Gilad Y. RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays. Genome Res. 2008;18(9):1509–17.

  17. 17.

    Liu S, Lin L, Jiang P, Wang D, Xing Y. A comparison of RNA-Seq and high-density exon array for detecting differential gene expression between closely related species. Nucleic Acids Res. 2010;39(2):578–88.

  18. 18.

    Eswaran J, Cyanam D, Mudvari P, Reddy SD, Pakala SB, Nair SS, et al. Transcriptomic landscape of breast cancers through mRNA sequencing. Sci Rep. 2012;2:264.

  19. 19.

    Huang Q, Lin B, Liu H, Ma X, Mo F, Yu W, et al. RNA-Seq analyses generate comprehensive transcriptomic landscape and reveal complex transcript patterns in hepatocellular carcinoma. PLoS One. 2011;6(10), e26168.

  20. 20.

    Ma S, Bao JY, Kwan PS, Chan YP, Tong CM, Fu L, et al. Identification of PTK6, via RNA sequencing analysis, as a suppressor of esophageal squamous cell carcinoma. Gastroenterology. 2012;143(3):675-686–e671-612.

  21. 21.

    Ren S, Peng Z, Mao JH, Yu Y, Yin C, Gao X, et al. RNA-seq analysis of prostate cancer in the Chinese population identifies recurrent gene fusions, cancer-associated long noncoding RNAs and aberrant alternative splicings. Cell Res. 2012;22(5):806–21.

  22. 22.

    Shah SP, Roth A, Goya R, Oloumi A, Ha G, Zhao Y, et al. The clonal and mutational evolution spectrum of primary triple-negative breast cancers. Nature. 2012;486(7403):395–9.

  23. 23.

    Twine NA, Janitz K, Wilkins MR, Janitz M. Whole transcriptome sequencing reveals gene expression and splicing differences in brain regions affected by Alzheimer's disease. PLoS One. 2011;6(1), e16266.

  24. 24.

    Zhang LQ, Cheranova D, Gibson M, Ding S, Heruth DP, Fang D, et al. RNA-seq reveals novel transcriptome of genes and their isoforms in human pulmonary microvascular endothelial cells treated with thrombin. PLoS One. 2012;7(2), e31229.

  25. 25.

    Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10(1):57–63.

  26. 26.

    Sultan M, Schulz MH, Richard H, Magen A, Klingenhoff A, Scherf M, et al. A global view of gene activity and alternative splicing by deep sequencing of the human transcriptome. Science. 2008;321(5891):956–60.

  27. 27.

    Venables JP. Unbalanced alternative splicing and its significance in cancer. Bioessays. 2006;28(4):378–86.

  28. 28.

    David CJ, Manley JL. Alternative pre-mRNA splicing regulation in cancer: pathways and programs unhinged. Genes Dev. 2010;24(21):2343–64.

  29. 29.

    van Dongen JJ, Langerak AW, Bruggemann M, et al. Design and standardization of PCR primers and protocols for detection of clonal immunoglobulin and T-cell receptor gene recombinations in suspect lymphoproliferations: report of the BIOMED-2 Concerted Action BMH4-CT98-3936. Leukemia. 2003;17:2257–317.

  30. 30.

    FastQC web site Babraham Bioinformatics. A quality control tool for high throughput sequence data available at http://www.bioinformatics.bbsrc.ac.uk/projects/fastqc.

  31. 31.

    Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25(9):1105–11.

  32. 32.

    Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012;7(3):562–78.

  33. 33.

    Jiang H, Wong WH. Statistical inferences for isoform expression in RNA-Seq. Bioinformatics. 2009;25(8):1026–32.

  34. 34.

    Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):R106.

  35. 35.

    Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29(9), e45.

  36. 36.

    Shen S, Park JW, Huang J, Dittmar KA, Lu ZX, Zhou Q, et al. MATS: a Bayesian framework for flexible detection of differential alternative splicing from RNA-Seq data. Nucleic Acids Res. 2012;40(8), e61.

  37. 37.

    Wang L, Wang S, Li W. RSeQC: quality control of RNA-seq experiments. Bioinformatics. 2012;28(16):2184–5.

  38. 38.

    Kapranov P, St Laurent G, Raz T, Ozsolak F, Reynolds CP, Sorensen PH, et al. The majority of total nuclear-encoded non-ribosomal RNA in a human cell is 'dark matter' un-annotated RNA. BMC Biol. 2010;8:149.

  39. 39.

    Ameur A, Zaghlool A, Halvardson J, Wetterbom A, Gyllensten U, Cavelier L, et al. Total RNA sequencing reveals nascent transcription and widespread co-transcriptional splicing in the human brain. Nat Struct Mol Biol. 2010;18(12):1435–40.

  40. 40.

    Ferreira PG, Jares P, Rico D, Gomez-Lopez G, Martinez-Trillos A, Villamor N, et al. Transcriptome characterization by RNA sequencing identifies a major molecular and clinical subdivision in chronic lymphocytic leukemia. Genome Res. 2014;24(2):212–26.

  41. 41.

    Kamiguti AS, Lee ES, Till KJ, Harris RJ, Glenn MA, Lin K, et al. The role of matrix metalloproteinase 9 in the pathogenesis of chronic lymphocytic leukaemia. Br J Haematol. 2004;125(2):128–40.

  42. 42.

    Choudhury A, Derkow K, Daneshmanesh AH, Mikaelsson E, Kiaii S, Kokhaei P, et al. Silencing of ROR1 and FMOD with siRNA results in apoptosis of CLL cells. Br J Haematol. 2010;151(4):327–35.

  43. 43.

    Soo C, Hu FY, Zhang X, Wang Y, Beanes SR, Lorenz HP, et al. Differential expression of fibromodulin, a transforming growth factor-beta modulator, in fetal skin development and scarless repair. Am J Pathol. 2000;157(2):423–33.

  44. 44.

    Hidalgo A, Peired AJ, Wild MK, Vestweber D, Frenette PS. Complete identification of E-selectin ligands on neutrophils reveals distinct functions of PSGL-1, ESL-1, and CD44. Immunity. 2007;26(4):477–89.

  45. 45.

    Ormsby T, Schlecker E, Ferdin J, Tessarz AS, Angelisova P, Koprulu AD, et al. Btk is a positive regulator in the TREM-1/DAP12 signaling pathway. Blood. 2011;118(4):936–45.

  46. 46.

    Gandhirajan RK, Poll-Wolbeck SJ, Gehrke I, Kreuzer KA. Wnt/beta-catenin/LEF-1 signaling in chronic lymphocytic leukemia (CLL): a target for current and potential therapeutic options. Curr Cancer Drug Targets. 2010;10(7):716–27.

  47. 47.

    Zen K, Liu Y, McCall IC, Wu T, Lee W, Babbin BA, et al. Neutrophil migration across tight junctions is mediated by adhesive interactions between epithelial coxsackie and adenovirus receptor and a junctional adhesion molecule-like protein on neutrophils. Mol Biol Cell. 2005;16(6):2694–703.

  48. 48.

    Decker S, Finter J, Forde AJ, Kissel S, Schwaller J, Mack TS, et al. PIM kinases are essential for chronic lymphocytic leukemia cell survival (PIM2/3) and CXCR4-mediated microenvironmental interactions (PIM1). Mol Cancer Ther. 2014;13(5):1231–45.

  49. 49.

    Moon E, Lee R, Near R, Weintraub L, Wolda S, Lerner A. Inhibition of PDE3B augments PDE4 inhibitor-induced apoptosis in a subset of patients with chronic lymphocytic leukemia. Clin Cancer Res. 2002;8(2):589–95.

  50. 50.

    Talab F, Allen JC, Thompson V, Lin K, Slupsky JR. LCK is an important mediator of B-cell receptor signaling in chronic lymphocytic leukemia cells. Mol Cancer Res. 2013;11(5):541–54.

  51. 51.

    Pekarsky Y, Palamarchuk A, Maximov V, Efanov A, Nazaryan N, Santanam U, et al. Tcl1 functions as a transcriptional regulator and is directly involved in the pathogenesis of CLL. Proc Natl Acad Sci U S A. 2008;105(50):19643–8.

  52. 52.

    Inada K, Okada S, Phuchareon J, Hatano M, Sugimoto T, Moriya H, et al. c-Fos induces apoptosis in germinal center B cells. J Immunol. 1998;161(8):3853–61.

  53. 53.

    Colotta F, Polentarutti N, Sironi M, Mantovani A. Expression and involvement of c-fos and c-jun protooncogenes in programmed cell death induced by growth factor deprivation in lymphoid cell lines. J Biol Chem. 1992;267(26):18278–83.

  54. 54.

    Shinohara M, Hiraki A, Ikebe T, Nakamura S, Kurahara S, Shirasuna K, et al. Immunohistochemical study of desmosomes in oral squamous cell carcinoma: correlation with cytokeratin and E-cadherin staining, and with tumour behaviour. J Pathol. 1998;184(4):369–81.

  55. 55.

    Yang L, Chen Y, Cui T, Knosel T, Zhang Q, Albring KF, et al. Desmoplakin acts as a tumor suppressor by inhibition of the Wnt/beta-catenin signaling pathway in human lung cancer. Carcinogenesis. 2012;33(10):1863–70.

  56. 56.

    Gilby DC, Sung HY, Winship PR, Goodeve AC, Reilly JT, Kiss-Toth E. Tribbles-1 and -2 are tumour suppressors, down-regulated in human acute myeloid leukaemia. Immunol Lett. 2010;130(1-2):115–24.

  57. 57.

    Bermudez O, Pages G, Gimond C. The dual-specificity MAP kinase phosphatases: critical roles in development and cancer. Am J Physiol Cell Physiol. 2010;299(2):C189–202.

  58. 58.

    Chiorazzi N, Efremov DG. Chronic lymphocytic leukemia: a tale of one or two signals? Cell Res. 2013;23(2):182–5.

  59. 59.

    Fernando RI, Litzinger M, Trono P, Hamilton DH, Schlom J, Palena C. The T-box transcription factor Brachyury promotes epithelial-mesenchymal transition in human tumor cells. J Clin Invest. 2010;120(2):533–44.

  60. 60.

    Guglielmi P, Davi F. Expression of a novel type of immunoglobulin C lambda transcripts in human mature B lymphocytes producing kappa light chains. Eur J Immunol. 1991;21(2):501–8.

  61. 61.

    Moravcikova E, Krepela E, Prochazka J, Rousalova I, Cermak J, Benkova K. Down-regulated expression of apoptosis-associated genes APIP and UACA in non-small cell lung carcinoma. Int J Oncol. 2012;40(6):2111–21.

  62. 62.

    Aktary Z, Pasdar M. Plakoglobin: Role in Tumorigenesis and Metastasis. Int J Cell Biol. 2012;2012.

  63. 63.

    Rehli M, Sulzbacher S, Pape S, Ravasi T, Wells CA, Heinz S, et al. Transcription factor Tfec contributes to the IL-4-inducible expression of a small group of genes in mouse macrophages including the granulocyte colony-stimulating factor receptor. J Immunol. 2005;174(11):7111–22.

  64. 64.

    Yu Y, Rabinowitz R, Polliack A, Ben-Bassat H, Schlesinger M. B-lymphocytes in CLL and NHL differ in the mRNA splicing pattern of the CD45 molecule. Eur J Haematol. 2000;64(6):376–84.

  65. 65.

    Vilpo J, Tobin G, Hulkkonen J, Hurme M, Thunberg U, Sundstrom C, et al. Surface antigen expression and correlation with variable heavy-chain gene mutation status in chronic lymphocytic leukemia. Eur J Haematol. 2003;70(1):53–9.

  66. 66.

    Hermiston ML, Xu Z, Weiss A. CD45: a critical regulator of signaling thresholds in immune cells. Annu Rev Immunol. 2003;21:107–37.

  67. 67.

    Russell SM, Sparrow RL, McKenzie IF, Purcell DF. Tissue-specific and allelic expression of the complement regulator CD46 is controlled by alternative splicing. Eur J Immunol. 1992;22(6):1513–8.

  68. 68.

    Courtois S, Caron de Fromentel C, Hainaut P. p53 protein variants: structural and functional similarities with p63 and p73 isoforms. Oncogene. 2004;23(3):631–8.

  69. 69.

    Kelleher FC, Rao A, Maguire A. Circadian molecular clocks and cancer. Cancer Lett. 2014;342(1):9–18.

Download references

Acknowledgements

SS is supported by a grant from Flight Attendants Medical Research Institute (FAMRI), and Veterans Administration Merit Research award. We thank the Broad Stem Cell Research Institute at UCLA for their help in high throughput RNA sequencing.

Author information

Affiliations

  1. Division of Hematology-Oncology, UCLA-VA Greater Los Angeles Healthcare System, Los Angeles, CA, USA

    • Wei Liao
    • , Gwen Jordaan
    •  & Sanjai Sharma

  2. Department of Pathology, VA Greater Los Angeles Healthcare System, Los Angeles, CA, USA

    • Phillipp Nham
    •  & Ryan T. Phan

  3. Department of Molecular, Cell and Developmental Biology, UCLA, Los Angeles, CA, USA

    • Matteo Pelegrini

  4. UCLA West Los Angeles VA Medical Center, 11301 Wilshire Blvd, Bldg 304, Rm E1-115, Los Angeles, CA, 90073, USA

    • Sanjai Sharma

Authors

  1. Search for Wei Liao in:

  2. Search for Gwen Jordaan in:

  3. Search for Phillipp Nham in:

  4. Search for Ryan T. Phan in:

  5. Search for Matteo Pelegrini in:

  6. Search for Sanjai Sharma in:

Corresponding author

Correspondence to Sanjai Sharma.

Additional information

Competing interests

The authors declare no competing interests as defined by this journal or any other conflicting interests.

Authors’ contributions

WL performed experiments, analyzed data and wrote the manuscript, GJ, PP, and RP performed experiments, MP, RP and SS designed experiments, supervised the study and wrote the manuscript. All authors have read and approved the manuscript.

Additional files

Additional file 1:

Alignment statistic summary of all 15 samples signal-end reads mapped to UCSC H. sapiens reference genome (build hg19) using Tophat alignment program. (DOCX 19 kb)

Additional file 2:

Number of transcripts and genes in B cells, U-CLL and M-CLL. Pair wise scatter plot matrix. (DOCX 148 kb)

Additional file 3:

Differentially expressed genes in CLL relative to B cells based on FPKM analysis. (XLSX 214 kb)

Additional file 4:

IPA functional annotation of differentially expressed genes in CLL specimens. (DOCX 19 kb)

Additional file 5:

List of genes in each functional pathway from IPA analysis. (XLSX 27 kb)

Additional file 6:

Differentially expressed genes in M- and U-CLL specimens. (XLSX 78 kb)

Additional file 7:

List of alternative spliced genes from MATS analysis, comparing CLL and B cell data. (XLSX 45 kb)

Additional file 8:

List of skipped exon events tested in CLL specimens with PCR analysis. (PPTX 119 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Liao, W., Jordaan, G., Nham, P. et al. Gene expression and splicing alterations analyzed by high throughput RNA sequencing of chronic lymphocytic leukemia specimens. BMC Cancer 15, 714 (2015) doi:10.1186/s12885-015-1708-9

Download citation

Keywords

  • CLL
  • RNA-sequencing
  • Differential gene expression
  • Alternative splicing

BMC Cancer

ISSN: 1471-2407

Contact us