Transcriptome instability in colorectal cancer identified by exon microarray analyses: Associations with splicing factor expression levels and patient survival
- Anita Sveen†1, 2,
- Trude H Ågesen†1, 2,
- Arild Nesbakken2, 3,
- Torleiv O Rognum4,
- Ragnhild A Lothe1, 2 and
- Rolf I Skotheim1, 2Email author
© Sveen et al.; licensee BioMed Central Ltd. 2011
Received: 14 April 2011
Accepted: 27 May 2011
Published: 27 May 2011
Colorectal cancer (CRC) is a heterogeneous disease that, on the molecular level, can be characterized by inherent genomic instabilities; chromosome instability and microsatellite instability. In the present study we analyze genome-wide disruption of pre-mRNA splicing, and propose transcriptome instability as a characteristic that is analogous to genomic instability on the transcriptome level.
Exon microarray profiles from two independent series including a total of 160 CRCs were investigated for their relative amounts of exon usage differences. Each exon in each sample was assigned an alternative splicing score calculated by the FIRMA algorithm. Amounts of deviating exon usage per sample were derived from exons with extreme splicing scores.
There was great heterogeneity within both series in terms of sample-wise amounts of deviating exon usage. This was strongly associated with the expression levels of approximately half of 280 splicing factors (54% and 48% of splicing factors were significantly correlated to deviating exon usage amounts in the two series). Samples with high or low amounts of deviating exon usage, associated with overall transcriptome instability, were almost completely separated into their respective groups by hierarchical clustering analysis of splicing factor expression levels in both sample series. Samples showing a preferential tendency towards deviating exon skipping or inclusion were associated with skewed transcriptome instability. There were significant associations between transcriptome instability and reduced patient survival in both sample series. In the test series, patients with skewed transcriptome instability showed the strongest prognostic association (P = 0.001), while a combination of the two characteristics showed the strongest association with poor survival in the validation series (P = 0.03).
We have described transcriptome instability as a characteristic of CRC. This transcriptome instability has associations with splicing factor expression levels and poor patient survival.
Colorectal cancer (CRC) is a prevalent disease with a world-wide incidence of more than one million new cases each year, making it the third most commonly diagnosed cancer among men and women . Colorectal tumors are heterogeneous and evolve through multiple pathways. Malignant transformation is dependent on the accumulation of numerous genetic changes over years. Such genetic instability provides a way of classifying tumors into different molecular subtypes . CRCs with the microsatellite instability (MSI) phenotype have a defective mismatch repair system, which results in a high proportion of mutations in nucleotide repeats (microsatellites) throughout the genome. Mismatch mutations of microsatellites located in coding-regions may exert tumorigenic effects - for example, by inactivating tumor suppressor genes [3, 4]. Cancers with MSI generally have few numeric changes or rearrangements at the chromosomal level. In contrast, CRCs with chromosomal instability (CIN) exhibit frequent chromosomal rearrangements and aneuploidy . MSI and CIN tumors represent two different types of virtually mutually exclusive genetic instabilities, and also two separate mechanistic mutational pathways for CRC development . A third subgroup, comprising both MSI and CIN tumors, is denoted as having a CpG island methylator phenotype (CIMP). Although these tumors often harbor BRAF mutations, CIMP is considered an epigenetically rather than genetically acting phenotype, and is characterized by extensive DNA methylation at promoter regions and associated gene silencing [6, 7]. Despite the fact that the above described refinements to CRC classification take into account molecular characteristics, it is evident that the total genetic definition of this heterogeneous disease is yet to be described .
Tumor staging remains the most important factor for determining the prognosis of CRC patients . Colorectal tumors are divided into four stages by the tumor-node-metastasis (TNM) system, taking into account depth of infiltration of the tumors, the extent of lymph node involvement, and the presence of distant metastases . During the past few years, molecular markers with potential prognostic value have been identified, several with relationships to the genetic and epigenetic types of instabilities mentioned above . However, no such markers are currently recommended for clinical use, and the need for discovery of novel biomarkers to aid the determination of prognosis in CRC patients remains great.
Alternative pre-mRNA splicing is an important source of functional diversity in the majority of mammalian genes . Nearly all multi-exon genes are expressed in multiple isoforms . Alternative splicing may occur as cassette alternative exons, mutually exclusive exons, intron retentions, or alternative 3' or 5' splice sites. Additionally, transcript variation may be caused by alternative promoter usage, resulting in alternative first exons, or alternative polyadenylation sites, presented as alternative terminal exons . Three consensus sequence elements are required for splicing to occur, two at each border of the intron, as well as the intronic branch site. Additionally, splicing is regulated by other cis-acting elements (exonic and intronic splicing silencers and enhancers) as well as the numerous trans-acting factors constituting the splicing machinery [14, 15]. The integrity of this wide range of elements is crucial for splicing fidelity and the maintenance of a stable and properly functioning transcriptome.
Aberrant splicing patterns have been associated with cancer . One example is the shift in expression towards the anti-apoptotic isoform (BCL-X L ) of the apoptosis regulator BCL-X . Cancer-specific transcript variation has also been found in CRC [18, 19]. Recently, it has been shown that global splicing patterns are likely to be distorted in several cancers , and preferential utilization of alternative splice sites is found to be a property of tumors [21, 22]. Moreover, differential expression of splicing factors is reported in cancer tissues compared to normal tissue . For some cancers - for example, CRC - this is, for the most part, seen as up-regulated expression . Despite previous publications on general changes in expression levels of splicing factors in cancer, it is yet to be determined what implications this imposes on the cancer phenotype.
In this study, we describe the disruption of alternative splicing as a global event in CRC tissue using exon microarray analysis. We show a great extent of variation in the amount of deviating exon usage among the samples. This transcriptome instability is found to be associated with prognosis in CRC patients, and also has clear associations with the expression levels of approximately half of all splicing factors.
Clinicopathological and molecular characteristics of the two independent colorectal cancer study populations
Test series (n= 83)
Validation series (n= 77)
Age at diagnosis (mean ± SD)
66.0 ± 11.7
72.7 ± 13.5
Sex (male; female)
Stage (II; III)
Location (right; left; rectum)
26; 25; 32
46; 20; 11
Mean follow-up, years (minimum; maximum)
6.7 (0.7; 10.0)
3.5 (0.2; 5.0)
Number of events (deaths from CRC)
Either sTIN, oTIN, or both
Both TIN phenotypes
Exon microarray analysis
RNA (1 μg) from each sample was individually amplified, reverse transcribed, fragmented, and labeled using the Affymetrix GeneChip® Whole Transcript (WT) Sense Target Labeling Assay . Labeled sense strand DNA was hybridized onto the Affymetrix GeneChip Human Exon 1.0 ST Array for 16 to 18 hours . Each array contains 1.4 million probe sets, of which 289,961 target well annotated full-length human mRNAs ('core' probe sets), and the remaining probe sets are derived from annotations of lower confidence levels, as well as computer predictions . A probe set corresponds approximately to one exon, and will be referred to as such herein. The arrays were finally washed, stained and scanned according to the manufacturer's protocol.
Scanning of the microarrays and preprocessing of raw image intensity data were controlled by the Affymetrix GeneChip Command Console software (version 1.0). For each microarray, the software generated cell intensity (CEL) files storing probe-level intensity data calculated from scanned image files containing pixel intensity values. CEL data files were used as input for preprocessing and alternative splicing detection with the Finding Isoforms using Robust Multichip Analysis (FIRMA) method  (Additional file 1). As part of the FIRMA method, the first two preprocessing steps of the microarrays were performed according to the robust multi-array average (RMA) approach, involving background correction of perfect match probes and inter-chip quantile normalization . The summarization step estimating gene expression levels was slightly modified from standard RMA, not taking into account the chip-exon effect, that is, ruling out the relative change for the sample in a particular exon. For this purpose, a custom made chip definition file containing 284,258 probe sets targeting exons belonging to the 'core' set of well annotated exons was downloaded from aroma.affymetrix . Applying this annotation file, the collective set of exons made up 18,708 transcript clusters, or genes. Alternative splicing scores, FIRMA scores, were calculated for each individual exon in each individual sample to represent a measure for whether differential exon usage has occurred. These scores were calculated as exon-level intensities deviating from the corresponding gene level, assessed indirectly as the residual after fitting the gene-level model to the actual data. Large residuals indicated differential expression of the particular exon compared to the corresponding gene level . The FIRMA scores were log-2 transformed. The microarray data can be accessed from NCBI's Gene Expression Omnibus (GEO) with the accession number [GEO:GSE24551].
To provide a global estimate of the relative amount of differential exon usage per sample, we counted the sample-wise numbers of probe set level FIRMA scores belonging to the upper and lower 1st percentiles of all FIRMA scores in the data sets.
For further statistical analyses, the software SPSS 15.0 (SPSS Inc., Chicago, IL, USA) was used. This includes t-statistics, multinomial logistic regression, generation of Kaplan-Meier plots, Cox regression analyses for calculation of hazard ratios (HR) and corresponding 95% confidence intervals (CI), Fisher's exact test, and Mantel-Cox test for equality of survival distributions. P-values < 0.05 were considered significant. Hierarchical clustering analysis was done using J-Express 2011 (MolMine AS, Bergen, Norway).
A list of 280 human splicing factors (Table S1 in Additional file 2) was created by combining results from the Gene Ontology project  and Swiss-Prot at the UniProt Knowledgebase  in July 2009. Using the AmiGO web application , the Gene Ontology database was searched for the terms 'nuclear mRNA splicing, via spliceosome' (GO:0000398) and 'spliceosomal complex' (GO:0005681). The ExPASy proteomics server  was used to search Swiss-Prot for human proteins with the terms 'splicing' and 'spliceosome'. Thirty-one additional genes were added to the list based on their splicing-related descriptions, as found using the GeneCards Human Gene Database . Gene level expression data for these splicing factor genes were independently obtained from the CEL files of the CRC samples in the two series. The expression data were summarized on background-corrected and quantile-normalized data using the RMA algorithm implemented in the Affymetrix Expression Console 1.1 software.
For comparison, 100 gene sets with 280 genes each were created by random sampling using the R statistical software (Additional file 1). Expression levels for these genes were obtained from the test series of CRC samples, in the same manner as for the splicing factor genes.
Variation in the amounts of aberrant alternative exon usage among colorectal cancer tissue samples
We consider tumors with large amounts of deviating exon usage as tumors with transcriptome instability (TIN). Tumors showing preferential exon skipping or inclusion - that is, having a skewed distribution in the amounts of deviating skipping and inclusion events (difference in relative amounts of deviating skipping and inclusion greater than ± 0.7; n = 12) - are considered to have a skewed TIN (sTIN) subtype. Tumors where the overall amounts of deviating exon usage differ from the average (total relative amounts of deviating exon usage greater than ± 1.0; n = 14) are considered to have an overall TIN (oTIN) subtype. The patients showed no significant associations between sTIN or oTIN and MSI status, tumor stage, patient age and gender, or tumor location (by multinomial logistic regression).
To investigate whether this large sample-to-sample variation in deviating exon usage amounts was cancer specific, comparisons were made between 13 paired CRC and normal colonic mucosa samples from the validation series. This comparison was conducted in a manner that explored differential exon usage in each sample pair individually, normalized against the background level of differential exon usage occurring in normal colonic mucosa (Additional file 1). The total relative amounts of deviating exon usage were significantly higher in the cancer samples compared to their normal counterparts (P = 0.003 by paired samples t-test; Figure S2A in Additional file 2). To ensure that the scores given by the FIRMA algorithm truly are sensitive to differential exon usage between CRC and normal colonic mucosa, we investigated the alternative splicing scores of splicing events that have previously been predicted and validated in exon microarray analyses [18, 19]. Of 13 exons known to be alternatively spliced between CRC and normal colonic mucosa (indicated with exon array probe set IDs in Table S2 in Additional file 2), 11 showed a mean difference in alternative splicing scores between the paired samples according to expectation (Figure S2B in Additional file 2). Moreover, across the 13 CRC samples, 36% of the probe sets targeting these alternatively spliced exons were assigned a score exceeding the upper or lower 1st percentile thresholds, and have accordingly been accounted for in the total relative amounts of deviating exon usage in CRC compared to normal colonic mucosa.
Correlation between oTIN and expression levels of splicing factors
Associations between TIN and poor patient survival
In a multivariate model including tumor stage, MSI status, patient age and gender, as well as tumor location, sTIN was the strongest independent prognostic predictor (HR = 3.5; 95% CI, 1.5 to 8.3; P = 0.004). Among the other variables, only tumor stage showed significant associations with patient prognosis in this model.
Patients whose tumors were characterized by oTIN had a slightly poorer survival rate than patients with average amounts of deviating exon usage, although this finding was not significant (Table S5A in Additional file 2). However, patients with either sTIN or oTIN cancers (n = 24) had significantly poorer survival than patients with TIN-negative cancers (HR = 2.1; 95% CI 1.1 to 3.9; P = 0.02). This difference was significant also in a multivariate model including tumor stage, MSI status, patient age and gender, as well as tumor location (HR = 3.2; 95% CI, 1.5 to 6.7; P = 0.002). Two patients had tumors with overlapping subtypes of TIN, that is, characterized by both sTIN and oTIN. These patients died from metastatic disease 0.8 and 4.3 years after surgical removal of their primary tumor.
Validation of TIN in an independent series of stage II and III colorectal cancers
As in the test series, there were no significant associations between either of the two TIN subtypes and MSI status, tumor stage, patient age or gender, or tumor location (by multinomial logistic regression).
Also in the validation series, there was a strong association between oTIN and expression levels of splicing factors. Here, the sample-wise expression levels of 48% of splicing factors (133 of 280) were significantly correlated to the total relative amounts of deviating exon usage (compared to 54% in the test sample series; Pearson correlation, P < 0.05; Figure S4A in Additional file 2). Furthermore, the significant shift towards stronger negative correlation among splicing factor genes compared to random gene sets was also indicated in the validation series. In fact, 129 of 133 (97%) significantly correlated splicing factor genes were negatively correlated.
As in the test series, samples in the validation series were separated according to the oTIN subtype by hierarchical clustering analysis of expression levels of splicing factor genes (Euclidean distance metrics, complete linkage; Figure S4B in Additional file 2). Although the amount of oTIN samples was higher than in the test sample series (n = 30 compared to n = 14), the oTIN samples were almost completely separated into the two respective groups (Figure S4C in Additional file 2). Also in the validation series, the clustering was independent of tumor stage and MSI status.
The strong associations between TIN and disease-specific survival among the patients in the test series were not as clearly indicated in the validation series (Table S5B in Additional file 2). Notably, the mean follow-up period of the patients was considerably shorter (3.5 years) than in the test series (6.7 years). However, using the same stratification thresholds as in the test series, patients whose tumors were characterized with both of the TIN subtypes (n = 11) had a significantly poorer 5-year survival rate (64%) than patients whose tumors were characterized with neither or only one of the TIN subtypes (n = 66, 91% survival rate, P = 0.03; Figure 5B; HR = 3.7; 95% CI, 1.1 to 13.6). In a multivariate model including tumor stage, MSI status, patient age and gender, as well as tumor location, TIN was close to being an independent predictor of poor prognosis (multivariate HR = 3.1; 95% CI 0.9 to 11.2; P = 0.08).
In this study we have described TIN, characterized by deviating mRNA splice variant patterns, in CRC. Two main subtypes of this characteristic are described, oTIN and sTIN. The oTIN subtype was demonstrated through great variation in the total amounts of deviating exon usage among CRC tissue samples. This characteristic was found to be associated with the expression levels of approximately half of all splicing factors in two independent sample series. The sTIN subtype separated the samples by the preferred usage of alternative exon skipping or inclusion. TIN was shown to be associated with poor patient survival in two independent sample series, although prognostic stratification was achieved by different TIN subtypes in the two series. In the test series, sTIN was shown to have the strongest prognostic value, while a combination of the two subtypes was the most strongly associated with poor patient survival in the validation series.
Genomic instability is an important classifier of colorectal tumors because of its associations with certain risk factors and clinical features . Such genomic instability includes the virtually non-redundant differentiation between CIN and MSI. More recently, a classifier of the epigenetic state of the genome, CIMP, has also been introduced as an important phenotype describing the molecular nature of CRC [37, 38]. In analogy to these molecular classifiers of the CRC genome, TIN is here proposed as a potential classifier of the CRC transcriptome.
Transcriptome instability has potential clinical value. From our analyses it is not evident which of the two TIN subtypes is more closely related to patient survival, as the results differed between the test and validation series. It is not clear whether this may be due to the different clinicopathological constitution of the two patient series. Notably, the mean follow-up period of the patients was considerably shorter in the validation series (3.5 years, compared to 6.7 years in the test series). These patients received adjuvant chemotherapy according to the current standard, whereas patients in the historical test series received no other treatment than surgery. The mean age at diagnosis was also different between the two patient series (66 versus 73 years). These factors may have contributed to the observed different associations between TIN and patient survival. Although we found no associations between MSI and TIN, the different amounts of patients with MSI tumors in the two series (16% versus 31%) may also have influence in this regard. However, both series showed significantly reduced survival for patients stratified according to TIN, suggesting that high amounts of inter-tumor deviations in exon usage patterns may indeed be associated with poor patient survival. We would welcome additional effort to verify the prognostic value of this molecular characteristic.
The amounts of samples assigned to either or both of the TIN subtypes varied between the two sample series. In the test series, 29% (24 of 83) of the tumors were characterized by TIN, compared to 56% in the validation series (43 of 77). Although resulting in an increased amount of TIN-positive samples, we chose to keep a fixed scoring threshold for TIN to avoid introducing subjective bias into the validation. Despite the different frequencies, the reproducibility of the strong associations between oTIN and expression levels of splicing factor genes provides strong evidence for a biological explanation of the observed variability in deviating exon usage amounts. However, it is not certain to what extent this has influenced the associations between TIN and patient survival in the two series.
The amounts of deviating exon usage were significantly higher in the cancer samples compared to paired normal colonic mucosa. A possible bias may have been introduced to these comparisons by the use of different RNA extraction protocols for the two sample groups. However, when analyzing known splicing events, 11 of 13 exons were found to have a mean difference in alternative splicing scores between the paired samples according to expectation. Also, 36% of the probe sets across the 13 CRC samples included in this analysis of known splicing events were designated as differentially spliced (exceeding the upper and lower 1st percentiles of alternative splicing scores). Accordingly, it seems evident that alternative exon usage is indeed reflected in the analytical approach, where individual exons in individual samples are scored according to the likelihood of differential splicing. For the majority of the known splicing events tested, the recurrence rate in CRC tissue is uncertain, due to limited sample numbers used for discovery . The cancer-specific splicing event of SLC39A14, however, is reported to have high cancer sensitivity . In accordance with this, exon 4A in this gene was designated as differentially excluded in all CRC samples relative to normal mucosa (<1st percentile of alternative splicing scores). In contrast, in the inter-tumor comparisons that are the main focus here, the majority of the cancer samples showed no signs of differential splicing for this exon, reflecting the nature of the current study, investigating exon usage variation among tumor samples and not between tumor and normal samples. Since the amounts of aberrant splicing were found to be higher in the cancer samples than their normal counterparts, we find it likely that TIN is most relevant to studies of cancer tissues. It remains uncertain whether this is specific for CRC or may be a common characteristic for cancers in other tissues as well.
It is striking that the expression levels of the majority of splicing factors negatively correlated to the numbers of deviating exon usage. Irrespective of the role of the splicing factor, whether it is predominantly a splicing enhancer or silencer, and irrespective of the type of exon usage event (exon skipping or inclusion), this was true for more than 95% of the significantly correlated splicing factor genes in two independent sample series. In fact, the correlations among splicing factor genes were significantly stronger in the direction of negative correlation than among random gene sets of equal size. This strong correlation is further supported by unsupervised hierarchical clustering analyses, which were based on the expression levels of the total set of splicing factors and separated the samples according to the oTIN phenotype in both series. Therefore, the association between low expression levels of splicing factors and increased variability in exon usage seems to be indicative of a critical role of splicing factor activity for the maintenance of a stable transcriptome. This suggests a biological rationale for the differences in exon usage observed in this study.
Consistent with this great variation in exon usage among CRC samples, other recent studies have also shown that the extent of alternative splicing far exceeds previous estimates. On average, individual multi-exon genes are suggested to undergo at least seven alternative splicing events across various human tissues and cell lines [12, 39]. This provides the possibility for tremendous variation in transcriptome composition. Several studies have reported genome-wide distortion of alternative splicing in cancer tissues [21, 22]. It has also been shown that alternative splicing in cancer exhibits both tissue dependency and dependency upon the type of splicing event considered . Intron retention and cassette alternative exons have been suggested to be more prevalent in normal tissues, and alternative 3' and 5' splice sites to occur more often in cancers. Various splicing events are also believed to occur at different levels in different cancer types. This indicates a complex distortion of exon splicing in cancer. Due to the composition of the exon microarrays and the nature of the splicing detection algorithm, the current study describes primarily the alternative skipping and inclusion of individual exons. In terms of splicing events, this essentially represents intron retention and cassette alternative exons. Patterns of mutual exclusion among exons are not considered. Transcript variation not attributed to alternative splicing, but rather to alternative promoter usage or polyadenylation sites, is also represented in the analysis, as all flanking exons are treated similarly to exons contained internally within transcripts. A detailed analysis of alternative usage of promoters and polyadenylation sites requires a detailed description of transcript structure and abundance, which is beyond the scope of this study.
Although the approach taken here to analyze sample-wise amounts of deviating exon usage was shown to correctly classify known splicing events, the main purpose of the study was to describe splicing on a genome-wide scale and does not allow for detailed analysis of individual transcript structures and individual splicing events. Furthermore, the current analyses did not intend to provide insights into the functional consequences of individual splicing events, that is, whether the predicted splicing events yield functionally different protein isoforms. There is increasing evidence that a great amount of expressed transcripts result from splicing noise [11, 40, 41]. This is true especially for the large number of non-specific, non-conserved, and non-abundant transcripts that are frequently subjected to degradation by the regulatory mechanism nonsense-mediated decay .
This study provides evidence for a high degree of variation among CRC samples with regards to amounts of differential exon usage. Based on this, we suggest TIN as a characteristic of CRC, which can be further dissected into the oTIN and sTIN subclasses. The oTIN subtype, reflecting sample-wise total relative amounts of deviating exon usage, is negatively correlated to the expression level of the majority of splicing factor encoding genes. Furthermore, analyses of corresponding clinical data demonstrate that TIN is associated with poor patient survival.
CpG island methylator phenotype
Finding Isoforms using Robust Multichip Analysis
overall transcriptome instability
robust multi-array average
skewed transcriptome instability
AS has a PhD grant from the Research Council at Rikshospitalet-Radiumhospitalet Health Enterprise (RAL). THÅ has a PhD grant from the Norwegian Cancer Society. The study has been financed by grants from the Norwegian Cancer Society (PR-2006-0442; RAL, PR-2007-0166; RIS).
- Ferlay J, Shin HR, Bray F, Forman D, Mathers C, Parkin DM: Estimates of worldwide burden of cancer in 2008: GLOBOCAN 2008. Int J Cancer. 2010, 127: 2893-2917. 10.1002/ijc.25516.View ArticlePubMedGoogle Scholar
- Jass JR: Classification of colorectal cancer based on correlation of clinical, morphological and molecular features. Histopathology. 2007, 50: 113-130. 10.1111/j.1365-2559.2006.02549.x.View ArticlePubMedGoogle Scholar
- Ionov Y, Peinado MA, Malkhosyan S, Shibata D, Perucho M: Ubiquitous somatic mutations in simple repeated sequences reveal a new mechanism for colonic carcinogenesis. Nature. 1993, 363: 558-561. 10.1038/363558a0.View ArticlePubMedGoogle Scholar
- Markowitz S, Wang J, Myeroff L, Parsons R, Sun L, Lutterbaugh J, Fan RS, Zborowska E, Kinzler KW, Vogelstein B: Inactivation of the type II TGF-beta receptor in colon cancer cells with microsatellite instability. Science. 1995, 268: 1336-1338. 10.1126/science.7761852.View ArticlePubMedGoogle Scholar
- Kinzler KW, Vogelstein B: Lessons from hereditary colorectal cancer. Cell. 1996, 87: 159-170. 10.1016/S0092-8674(00)81333-1.View ArticlePubMedGoogle Scholar
- Samowitz WS: Genetic and epigenetic changes in colon cancer. Exp Mol Pathol. 2008, 85: 64-67. 10.1016/j.yexmp.2008.03.008.View ArticlePubMedGoogle Scholar
- Hawkins N, Norrie M, Cheong K, Mokany E, Ku SL, Meagher A, O'Connor T, Ward R: CpG island methylation in sporadic colorectal cancers and its relationship to microsatellite instability. Gastroenterology. 2002, 122: 1376-1387. 10.1053/gast.2002.32997.View ArticlePubMedGoogle Scholar
- Wolpin BM, Mayer RJ: Systemic treatment of colorectal cancer. Gastroenterology. 2008, 134: 1296-1310. 10.1053/j.gastro.2008.02.098.View ArticlePubMed CentralPubMedGoogle Scholar
- American Joint Committee on Cancer: Colon and rectum. AJCC Cancer Staging Manual. Edited by: Greene FL, Page DL, Fleming ID, Fritz AG, Balch CM, Haller DG, Morrow M. 2002, New York: Springer-Verlag, 113-124.View ArticleGoogle Scholar
- Huerta S: Recent advances in the molecular diagnosis and prognosis of colorectal cancer. Expert Rev Mol Diagn. 2008, 8: 277-288. 10.1586/1473718.104.22.1687.View ArticlePubMedGoogle Scholar
- Graveley BR: Alternative splicing: increasing diversity in the proteomic world. Trends Genet. 2001, 17: 100-107. 10.1016/S0168-9525(00)02176-4.View ArticlePubMedGoogle Scholar
- Pan Q, Shai O, Lee LJ, Frey BJ, Blencowe BJ: Deep surveying of alternative splicing complexity in the human transcriptome by high-throughput sequencing. Nat Genet. 2008, 40: 1413-1415. 10.1038/ng.259.View ArticlePubMedGoogle Scholar
- Blencowe BJ: Alternative splicing: new insights from global analyses. Cell. 2006, 126: 37-47. 10.1016/j.cell.2006.06.023.View ArticlePubMedGoogle Scholar
- Chen M, Manley JL: Mechanisms of alternative splicing regulation: insights from molecular and genomics approaches. Nat Rev Mol Cell Biol. 2009, 10: 741-754.PubMed CentralPubMedGoogle Scholar
- Query CC, Konarska MM: Splicing fidelity revisited. Nat Struct Mol Biol. 2006, 13: 472-474. 10.1038/nsmb0606-472.View ArticlePubMedGoogle Scholar
- Venables JP: Aberrant and alternative splicing in cancer. Cancer Res. 2004, 64: 7647-7654. 10.1158/0008-5472.CAN-04-1910.View ArticlePubMedGoogle Scholar
- Mercatante DR, Bortner CD, Cidlowski JA, Kole R: Modification of alternative splicing of Bcl-x pre-mRNA in prostate and breast cancer cells. Analysis of apoptosis and cell death. J Biol Chem. 2001, 276: 16411-16417. 10.1074/jbc.M009256200.View ArticlePubMedGoogle Scholar
- Gardina PJ, Clark TA, Shimada B, Staples MK, Yang Q, Veitch J, Schweitzer A, Awad T, Sugnet C, Dee S, Davies C, Williams A, Turpaz Y: Alternative splicing and differential gene expression in colon cancer detected by a whole genome exon array. BMC Genomics. 2006, 7: 325-10.1186/1471-2164-7-325.View ArticlePubMed CentralPubMedGoogle Scholar
- Thorsen K, Mansilla F, Schepeler T, Øster B, Rasmussen MH, Dyrskjøt L, Karni R, Akerman M, Krainer AR, Laurberg S, Andersen CL, Ørntoft TF: Alternative splicing of SLC39A14 in colorectal cancer is regulated by the Wnt pathway. Mol Cell Proteomics. 2011, 10: M110.002998-10.1074/mcp.M110.002998.View ArticlePubMed CentralPubMedGoogle Scholar
- Melamud E, Moult J: Stochastic noise in splicing machinery. Nucleic Acids Res. 2009, 37: 4873-4886. 10.1093/nar/gkp471.View ArticlePubMed CentralPubMedGoogle Scholar
- He C, Zhou F, Zuo Z, Cheng H, Zhou R: A global view of cancer-specific transcript variants by subtractive transcriptome-wide analysis. PLoS ONE. 2009, 4: e4732-10.1371/journal.pone.0004732.View ArticlePubMed CentralPubMedGoogle Scholar
- Ritchie W, Granjeaud S, Puthier D, Gautheret D: Entropy measures quantify global splicing disorders in cancer. PLOS Comput Biol. 2008, 4: e1000011-10.1371/journal.pcbi.1000011.View ArticlePubMed CentralPubMedGoogle Scholar
- Skotheim RI, Nees M: Alternative splicing in cancer: Noise, functional, or systematic?. Int J Biochem Cell Biol. 2007, 39: 1432-1449. 10.1016/j.biocel.2007.02.016.View ArticlePubMedGoogle Scholar
- Kirschbaum-Slager N, Lopes GMP, Galante PAF, Riggins GJ: Splicing factors are differentially expressed in tumors. Genet Mol Res. 2004, 3: 512-520.PubMedGoogle Scholar
- Thorstensen L, Diep CB, Meling GI, Aagesen TH, Ahrens CH, Rognum TO, Lothe RA: WNT1 inducible signaling pathway protein 3, WISP-3, a novel target gene in colorectal carcinomas with microsatellite instability. Gastroenterology. 2001, 121: 1275-1280. 10.1053/gast.2001.29570.View ArticlePubMedGoogle Scholar
- GeneChip Whole Transcript (WT) Sense Target Labeling Assay Manual. [http://www.affymetrix.com/support/downloads/manuals/wt_sensetarget_label_manual.pdf]
- GeneChip Exon Array Design. [http://www.affymetrix.com/support/technical/technotes/exon_array_design_technote.pdf]
- Exon Probeset Annotations and Transcript Cluster Groupings. [http://www.affymetrix.com/support/technical/whitepapers/exon_probeset_trans_clust_whitepaper.pdf]
- Purdom E, Simpson KM, Robinson MD, Conboy JG, Lapuk AV, Speed TP: FIRMA: a method for detection of alternative splicing from exon array data. Bioinformatics. 2008, 24: 1707-1714. 10.1093/bioinformatics/btn284.View ArticlePubMed CentralPubMedGoogle Scholar
- Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003, 4: 249-264. 10.1093/biostatistics/4.2.249.View ArticlePubMedGoogle Scholar
- HuEx-1_0-st-v2,coreR3,A20071112,EP.cdf. [http://bcgc.lbl.gov/cdfFiles/HuEx-1_0-st-v2,A20071112,EP/HuEx-1_0-st-v2,coreR3,A20071112,EP.cdf]
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25: 25-29. 10.1038/75556.View ArticlePubMed CentralPubMedGoogle Scholar
- The UniProt Consortium: The Universal Protein Resource (UniProt) 2009. Nucleic Acids Res. 2008, 37: D169-D174.View ArticlePubMed CentralGoogle Scholar
- Carbon S, Ireland A, Mungall CJ, Shu SQ, Marshall B, Lewis S, the AmiGO Hub, the Web Presence Working Group: AmiGO: online access to ontology and annotation data. Bioinformatics. 2009, 25: 288-289. 10.1093/bioinformatics/btn615.View ArticlePubMed CentralPubMedGoogle Scholar
- Gasteiger E, Gattiker A, Hoogland C, Ivanyi I, Appel RD, Bairoch A: ExPASy: The proteomics server for in-depth protein knowledge and analysis. Nucleic Acids Res. 2003, 31: 3784-3788. 10.1093/nar/gkg563.View ArticlePubMed CentralPubMedGoogle Scholar
- Rebhan M, Chalifa-Caspi V, Prilusky J, Lancet D: GeneCards: a novel functional genomics compendium with automated data mining and query reformulation support. Bioinformatics. 1998, 14: 656-664. 10.1093/bioinformatics/14.8.656.View ArticlePubMedGoogle Scholar
- Toyota M, Ahuja N, Ohe-Toyota M, Herman JG, Baylin SB, Issa JP: CpG island methylator phenotype in colorectal cancer. Proc Natl Acad Sci USA. 1999, 96: 8681-8686. 10.1073/pnas.96.15.8681.View ArticlePubMed CentralPubMedGoogle Scholar
- Toyota M, Ohe-Toyota M, Ahuja N, Issa JP: Distinct genetic profiles in colorectal tumors with or without the CpG island methylator phenotype. Proc Natl Acad Sci USA. 2000, 97: 710-715. 10.1073/pnas.97.2.710.View ArticlePubMed CentralPubMedGoogle Scholar
- Wang ET, Sandberg R, Luo S, Khrebtukova I, Zhang L, Mayr C, Kingsmore SF, Schroth GP, Burge CB: Alternative isoform regulation in human tissue transcriptomes. Nature. 2008, 456: 470-476. 10.1038/nature07509.View ArticlePubMed CentralPubMedGoogle Scholar
- Hiller M, Szafranski K, Backofen R, Platzer M: Alternative splicing at NAGNAG acceptors: Simply noise or noise and more?. PLoS Genet. 2006, 2: 1944-View ArticleGoogle Scholar
- Sugnet CW, Kent WJ, Ares M, Haussler D: Transcriptome and genome conservation of alternative splicing events in humans and mice. Pac Symp Biocomput. 2004, 66-77.Google Scholar
- Bengtsson H, Simpson K, Bullard J, Hansen K: aroma.affymetrix: A generic framework in R for analyzing small to very large Affymetrix data sets in bounded memory. Technical Report 745. 2008, Department of Statistics, University of California, Berkeley, [http://statistics.berkeley.edu/tech-reports/745.pdf]Google Scholar
- aroma.affymetrix. [http://www.aroma-project.org/]
- Langfelder P, Horvath S: WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008, 9: 559-10.1186/1471-2105-9-559.View ArticlePubMed CentralPubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.