Alternative splicing in colon, bladder, and prostate cancer identified by exon array analysis

A genome-wide aberrant RNA splicing in patients with acute myeloid leukemia identifies novel potential disease markers and therapeutic targets

PURPOSE

Despite new treatments, acute myeloid leukemia (AML) remains an incurable disease. More effective drug design requires an expanded view of the molecular complexity that underlies AML. Alternative splicing of RNA is used by normal cells to generate protein diversity. Growing evidence indicates that aberrant splicing of genes plays a key role in cancer. We investigated genome-wide splicing abnormalities in AML and based on these abnormalities, we aimed to identify novel potential biomarkers and therapeutic targets.

EXPERIMENTAL DESIGN

We used genome-wide alternative splicing screening to investigate alternative splicing abnormalities in two independent AML patient cohorts [Dana-Farber Cancer Institute (DFCI) (Boston, MA) and University Hospital de Nantes (UHN) (Nantes, France)] and normal donors. Selected splicing events were confirmed through cloning and sequencing analysis, and than validated in 193 patients with AML.

RESULTS

Our results show that approximately 29% of expressed genes genome-wide were differentially and recurrently spliced in patients with AML compared with normal donors bone marrow CD34(+) cells. Results were reproducible in two independent AML cohorts. In both cohorts, annotation analyses indicated similar proportions of differentially spliced genes encoding several oncogenes, tumor suppressor proteins, splicing factors, and heterogeneous-nuclear-ribonucleoproteins, proteins involved in apoptosis, cell proliferation, and spliceosome assembly. Our findings are consistent with reports for other malignances and indicate that AML-specific aberrations in splicing mechanisms are a hallmark of AML pathogenesis.

CONCLUSIONS

Overall, our results suggest that aberrant splicing is a common characteristic for AML. Our findings also suggest that splice variant transcripts that are the result of splicing aberrations create novel disease markers and provide potential targets for small molecules or antibody therapeutics for this disease.

Review: Alternative Splicing (AS) of Genes As An Approach for Generating Protein Complexity

Prior to the completion of the human genome project, the human genome was thought to have a greater number of genes as it seemed structurally and functionally more complex than other simpler organisms. This along with the belief of “one gene, one protein”, were demonstrated to be incorrect. The inequality in the ratio of gene to protein formation gave rise to the theory of alternative splicing (AS). AS is a mechanism by which one gene gives rise to multiple protein products. Numerous databases and online bioinformatic tools are available for the detection and analysis of AS. Bioinformatics provides an important approach to study mRNA and protein diversity by various tools such as expressed sequence tag (EST) sequences obtained from completely processed mRNA. Microarrays and deep sequencing approaches also aid in the detection of splicing events. Initially it was postulated that AS occurred only in about 5% of all genes but was later found to be more abundant. Using bioinformatic approaches, the level of AS in human genes was found to be fairly high with 35-59% of genes having at least one AS form. Our ability to determine and predict AS is important as disorders in splicing patterns may lead to abnormal splice variants resulting in genetic diseases. In addition, the diversity of proteins produced by AS poses a challenge for successful drug discovery and therefore a greater understanding of AS would be beneficial.

An integrative framework identifies alternative splicing events in colorectal cancer development

Alternative splicing (AS) is a common mechanism which creates diverse RNA isoforms from a single gene, potentially increasing protein variety. Growing evidence suggests that this mechanism is closely related to cancer progression. In this study, whole transcriptome analysis was performed with GeneChip Human exon 1.0 ST Array from 80 samples comprising 23 normal colon mucosa, 30 primary colorectal cancer and 27 liver metastatic specimens from 46 patients, to identify AS events in colorectal cancer progression. Differentially expressed genes and exons were estimated and AS events were reconstructed by combining exon-level analyses with AltAnalyze algorithms and transcript-level estimations (MMBGX probabilistic method). The number of AS genes in the transition from normal colon mucosa to primary tumor was the most abundant, but fell considerably in the next transition to liver metastasis. 206 genes with probable AS events in colon cancer development and progression were identified, that are involved in processes and pathways relevant to tumor biology, as cell-cell and cell-matrix interactions. Several AS events in VCL, CALD1, B3GNT6 and CTHRC1 genes, differentially expressed during tumor development were validated, at RNA and at protein level. Taken together, these results demonstrate that cancer-specific AS is common in early phases of colorectal cancer natural history.

RNA sequencing of cancer reveals novel splicing alterations

Breast cancer transcriptome acquires a myriad of regulation changes, and splicing is critical for the cell to “tailor-make” specific functional transcripts. We systematically revealed splicing signatures of the three most common types of breast tumors using RNA sequencing: TNBC, non-TNBC and HER2-positive breast cancer. We discovered subtype specific differentially spliced genes and splice isoforms not previously recognized in human transcriptome. Further, we showed that exon skip and intron retention are predominant splice events in breast cancer. In addition, we found that differential expression of primary transcripts and promoter switching are significantly deregulated in breast cancer compared to normal breast. We validated the presence of novel hybrid isoforms of critical molecules like CDK4, LARP1, ADD3, and PHLPP2. Our study provides the first comprehensive portrait of transcriptional and splicing signatures specific to breast cancer sub-types, as well as previously unknown transcripts that prompt the need for complete annotation of tissue and disease specific transcriptome.

Splice variants of zinc finger protein 695 mRNA associated to ovarian cancer

Background

Studies of alternative mRNA splicing (AS) in health and disease have yet to yield the complete picture of protein diversity and its role in physiology and pathology. Some forms of cancer appear to be associated to certain alternative mRNA splice variants, but their role in the cancer development and outcome is unclear.

Methods

We examined AS profiles by means of whole genome exon expression microarrays (Affymetrix GeneChip 1.0) in ovarian tumors and ovarian cancer-derived cell lines, compared to healthy ovarian tissue. Alternatively spliced genes expressed predominantly in ovarian tumors and cell lines were confirmed by RT-PCR.

Results

Among several significantly overexpressed AS genes in malignant ovarian tumors and ovarian cancer cell lines, the most significant one was that of the zinc finger protein ZNF695, with two previously unknown mRNA splice variants identified in ovarian tumors and cell lines. The identity of ZNF695 AS variants was confirmed by cloning and sequencing of the amplicons obtained from ovarian cancer tissue and cell lines.

Conclusions

Alternative ZNF695 mRNA splicing could be a marker of ovarian cancer with possible implications on its pathogenesis.

ToP: a trend-of-disease-progression procedure works well for identifying cancer genes from multi-state cohort gene expression data for human colorectal cancer

Significantly expressed genes extracted from microarray gene expression data have proved very useful for identifying genetic biomarkers of diseases, including cancer. However, deriving a disease related inference from a list of differentially expressed genes has proven less than straightforward. In a systems disease such as cancer, how genes interact with each other should matter just as much as the level of gene expression. Here, in a novel approach, we used the network and disease progression properties of individual genes in state-specific gene-gene interaction networks (GGINs) to select cancer genes for human colorectal cancer (CRC) and obtain a much higher hit rate of known cancer genes when compared with methods not based on network theory. We constructed GGINs by integrating gene expression microarray data from multiple states--healthy control (Nor), adenoma (Ade), inflammatory bowel disease (IBD) and CRC--with protein-protein interaction database and Gene Ontology. We tracked changes in the network degrees and clustering coefficients of individual genes in the GGINs as the disease state changed from one to another. From these we inferred the state sequences Nor-Ade-CRC and Nor-IBD-CRC both exhibited a trend of (disease) progression (ToP) toward CRC, and devised a ToP procedure for selecting cancer genes for CRC. Of the 141 candidates selected using ToP, ∼50% had literature support as cancer genes, compared to hit rates of 20% to 30% for standard methods using only gene expression data. Among the 16 candidate cancer genes that encoded transcription factors, 13 were known to be tumorigenic and three were novel: CDK1, SNRPF, and ILF2. We identified 13 of the 141 predicted cancer genes as candidate markers for early detection of CRC, 11 and 2 at the Ade and IBD states, respectively.

Caldesmon: new insights for diagnosing endometriosis

Considerable effort has been invested in searching for less invasive methods of diagnosing endometriosis. Previous studies have indicated altered levels of the CALD1 gene (encoding the protein caldesmon) in endometriosis. The aims of our study were to investigate whether average CALD1 expression and caldesmon protein levels are differentially altered in the endometrium and endometriotic lesions and to evaluate the performance of the CALD1 gene and caldesmon protein as potential biomarkers for endometriosis. Paired biopsies of endometrial tissue (eutopic endometrium) and endometriotic lesions (ectopic endometrium) were obtained from patients with endometriosis to evaluate CALD1 gene expression and caldesmon protein levels by real-time PCR and Western blot analysis, respectively. In addition, immunostaining for caldesmon to determine cellular localization was also performed. Endometrium from women without endometriosis was used as a control. Increased CALD1 expression and caldesmon levels were detected in the endometriotic lesions. The electrophoretic profile of caldesmon by Western blot analysis was clearly different between the control group (endometrium of women without endometriosis) and the group of women with endometriosis (eutopic endometrium and endometriotic lesions). Caldesmon expression as determined by immunostaining showed no variation among the cell types in endometriotic lesions and eutopic endometrium. Stromal cells marked positively in eutopic endometrium from control patients and in the endometriotic lesions. The presence of caldesmon in the endometrium of patients with and without endometriosis permitted diagnoses with 95% sensitivity (specificity 100%) and 100% sensitivity (specificity 100%) for the disease and for minimal to mild endometriosis in the proliferative phase of the menstrual cycle, respectively. In the secretory phase, minimal to mild endometriosis was detected with 90% sensitivity and 93.3% specificity. Caldesmon is a possible predictor of endometrial dysregulation in patients with endometriosis. A potential limitation of our study is the fact that other endometrial diseases were not excluded, and therefore prospective studies are needed to confirm the potential of caldesmon as a biomarker exclusively for endometriosis.

SRp20: an overview of its role in human diseases

Alternative splicing in mRNA maturation has emerged as a major field of study also because of its implications in various diseases. The SR proteins play an important role in the regulation of this process. Evidence indicates that SRp20 (SFSR3), the smallest member of the SR protein family, is involved in numerous biological processes. Here we review the state-of-the-art of knowledge about the SR proteins, in particular SRp20, in terms of its function and misregulation in human diseases including cancer also in view of its potential as a therapeutic target.

RRP1B is a metastasis modifier that regulates the expression of alternative mRNA isoforms through interactions with SRSF1

RRP1B (ribosomal RNA processing 1 homolog B) was first identified as a metastasis susceptibility gene in breast cancer through its ability to modulate gene expression in a manner that can be used to accurately predict prognosis in breast cancer. However, the mechanism(s) by which RRP1B modulates gene expression is currently unclear. Many RRP1B binding candidates are involved in alternative splicing, a mechanism of gene expression regulation that is increasingly recognized to be involved in cancer progression and metastasis. One such target is SRSF1 (SF2/ASF), an essential splicing regulator that also functions as an oncoprotein. Earlier studies demonstrated that splicing and transcription occur concurrently and are coupled processes. Given that RRP1B regulates transcriptional activity, we hypothesized that RRP1B also regulates the expression of alternative mRNA isoforms through its interaction with SRSF1. Interaction between RRP1B and SRSF1 was verified by co-immunoprecipitation and co-immunofluorescence. Treatment of cells with transcriptional inhibitors significantly increased this interaction, demonstrating that the association of these two proteins is transcriptionally regulated. To assess the role of RRP1B in the regulation of alternative isoform expression, RNA-seq data were generated from control and Rrp1b-knockdown cells. Knockdown of Rrp1b induced a significant change in isoform expression in over 600 genes compared to control cell lines. This was verified by qRT-PCR using isoform-specific primers. Pathway enrichment analyses identified cell cycle and checkpoint regulation to be those most affected by Rrp1b knockdown. These data suggest that RRP1B suppresses metastatic progression by altering the transcriptome through its interaction with splicing regulators such as SRSF1.

The impact of splicing on protein domain architecture

Many proteins are composed of protein domains, functional units of common descent. Multidomain forms are common in all eukaryotes making up more than half of the proteome and the evolution of novel domain architecture has been accelerated in metazoans. It is also becoming increasingly clear that alternative splicing is prevalent among vertebrates. Given that protein domains are defined as structurally, functionally and evolutionarily distinct units, one may speculate that some alternative splicing events may lead to clean excisions of protein domains, thus generating a number of different domain architectures from one gene template. However, recent findings indicate that smaller alternative splicing events, in particular in disordered regions, might be more prominent than domain architectural changes. The problem of identifying protein isoforms is, however, still not resolved. Clearly, many splice forms identified through detection of mRNA sequences appear to produce 'nonfunctional' proteins, such as proteins with missing internal secondary structure elements. Here, we review the state of the art methods for identification of functional isoforms and present a summary of what is known, thus far, about alternative splicing with regard to protein domain architectures.

Prognostic significance of aberrantly silenced ANPEP expression in prostate cancer

Background:

Novel biomarkers for prostate cancer (PC) are urgently needed. This study investigates the expression, epigenetic regulation, and prognostic potential of ANPEP in PC.

Methods:

Aminopeptidase N (APN; encoded by ANPEP) expression was analysed by immunohistochemistry using tissue microarrays representing 267 radical prostatectomy (RP) and 111 conservatively treated (CT) PC patients. Clinical end points were recurrence-free survival (RFS) and cancer-specific survival (CSS), respectively. The ANPEP promoter methylation levels were determined by bisulphite sequencing or MethyLight analysis in 278 nonmalignant and PC tissue samples, and in cell lines.

Results:

The APN expression was significantly downregulated in PC compared with nonmalignant prostate tissue samples. Aberrant promoter hypermethylation was frequently observed in PC tissue samples, and 5-aza-2′-deoxycytidine induced ANPEP expression in three hypermethylated prostate cell lines, suggesting epigenetic silencing. Negative APN immunoreactivity was significantly associated with short RFS and short CSS in the RP and CT cohort, respectively, independently of routine clinicopathological predictors. Combining APN with a known angiogenesis marker (vascular endothelial growth factor or microvessel density) improved risk prediction significantly in both cohorts.

Conclusion:

Our results suggest negative APN immunoreactivity as a new independent adverse prognostic factor for patients with clinically localised PC and, furthermore, that epigenetic mechanisms are involved in silencing of ANPEP in PC.

Detection of Alternatively Spliced or Processed RNAs in Cancer Using Oligonucleotide Microarray

Deregulation of gene expression plays a pivotal role in tumorigenesis, so the ability to detect RNA alterations is of great value in cancer diagnosis and management. DNA microarrays have been used to measure changes in mRNA or microRNA level, but less often the change of RNA isoforms. Here we appraise the utilization of microarray in detecting alternatively processed RNAs, which have alternative splice forms, retained introns, or altered 3' untranslated regions. We cover the methodology and focus on cancer studies. Recent development in parallel or deep sequencing used in transcriptome analysis is also discussed.

Up-regulated expression of l-caldesmon associated with malignancy of colorectal cancer

BACKGROUND

Caldesmon (CaD), a major actin-associated protein, is found in smooth muscle and non-muscle cells. Smooth muscle caldesmon, h-CaD, is a multifunctional protein, and non-muscle cell caldesmon, l-CaD, plays a role in cytoskeletal architecture and dynamics. h-CaD is thought to be an useful marker for smooth muscle tumors, but the role(s) of l-CaD has not been examined in tumors.

METHODS

Primary colon cancer and liver metastasis tissues were obtained from colon cancer patients. Prior to chemoradiotherapy (CRT), normal and cancerous tissues were obtained from rectal cancer patients. Whole-tissue protein extracts were analyzed by 2-DE-based proteomics. Expression and phosphorylation level of main cellular signaling proteins were determined by western blot analysis. Cell proliferation after CaD siRNA transfection was monitored by MTT assay.

RESULTS

The expression level of l-CaD was significantly increased in primary colon cancer and liver metastasis tissues compared to the level in the corresponding normal tissues. In cancerous tissues obtained from the patients showing poor response to CRT (Dworak grade 4), the expression of l-CaD was increased compared to that of good response group (Dworak grade 1). In line with, l-CaD positive human colon cancer cell lines were more resistant to 5-fluorouracil (5-FU) and radiation treatment compared to l-CaD negative cell lines. Artificial suppression of l-CaD increased susceptibility of colon cancer cells to 5-FU, and caused an increase of p21 and c-PARP, and a decrease of NF-kB and p-mTOR expression.

CONCLUSION

Up-regulated expression of l-CaD may have a role for increasing metastatic property and decreasing CRT susceptibility in colorectal cancer cells.

Placenta-specific novel splice variants of Rho GDP dissociation inhibitor β are highly expressed in cancerous cells

Background

Alternative splicing of pre-mRNA transcripts not only plays a role in normal molecular processes but is also associated with cancer development. While normal transcripts are ubiquitously expressed in normal tissues, splice variants created through abnormal alternative splicing events are often expressed in cancer cells. Although the Rho GDP dissociation inhibitor β (ARHGDIB) gene has been found to be ubiquitously expressed in normal tissues and involved in cancer development, the presence of splice variants of ARHGDIB has not yet been investigated.

Results

Validation analysis for the presence of and exon structures of splice variants of ARHGDIB, performed using reverse-transcriptase polymerase chain reaction and DNA sequencing, successfully identified novel splice variants of ARHGDIB, that is, 6a, 6b, and 6c, in colon, pancreas, stomach, and breast cancer cell lines. Quantitative real-time polymerase chain reaction analysis showed that these variants were also highly expressed in normal placental tissue but not in other types of normal tissue.

Conclusions

Expression of ARHGDIB variants 6a, 6b, and 6c appears to be restricted to cancer cells and normal placental tissue, suggesting that these variants possess cancer-specific functions and, as such, are potential cancer-related biomarkers.

The cancer exome generated by alternative mRNA splicing dilutes predicted HLA class I epitope density

Several studies have shown that cancers actively regulate alternative splicing. Altered splicing mechanisms in cancer lead to cancer-specific transcripts different from the pool of transcripts occurring only in healthy tissue. At the same time, altered presentation of HLA class I epitopes is frequently observed in various types of cancer. Down-regulation of genes related to HLA class I antigen processing has been observed in several cancer types, leading to fewer HLA class I antigens on the cell surface. Here, we use a peptidome wide analysis of predicted alternative splice forms, based on a publicly available database, to show that peptides over-represented in cancer splice variants comprise significantly fewer predicted HLA class I epitopes compared to peptides from normal transcripts. Peptides over-represented in cancer transcripts are in the case of the three most common HLA class I supertype representatives consistently found to contain fewer predicted epitopes compared to normal tissue. We observed a significant difference in amino acid composition between protein sequences associated with normal versus cancer tissue, as transcripts found in cancer are enriched with hydrophilic amino acids. This variation contributes to the observed significant lower likelihood of cancer-specific peptides to be predicted epitopes compared to peptides found in normal tissue.

Analysis of molecular intra-patient variation and delineation of a prognostic 12-gene signature in non-muscle invasive bladder cancer; technology transfer from microarrays to PCR

Background:

Multiple clinical risk factors and genetic profiles have been demonstrated to predict progression of non-muscle invasive bladder cancer; however, no easily clinical applicable gene signature has been developed to predict disease progression independent of disease stage and grade.

Methods:

We measured the intra-patient variation of an 88-gene progression signature using 39 metachronous tumours from 17 patients. For delineation of the optimal quantitative reverse transcriptase PCR panel of markers, we used 115 tumour samples from patients in Denmark, Sweden, UK and Spain.

Results:

Analysis of intra-patient variation of the molecular markers showed 71% similar classification results. A final panel of 12 genes was selected, showing significant correlation with outcome. In multivariate Cox regression analysis, we found that the 12-gene signature was an independent prognostic factor (hazard ratio=7.4 (95% confidence interval: 3.4–15.9), P<0.001) when adjusting for stage, grade and treatment. Independent validation of the 12-gene panel and the determined cut-off values is needed and ongoing.

Conclusion:

Intra-patient marker variation in metachronous tumours is present. Therefore, to increase test sensitivity, it may be necessary to test several metachronous tumours from a patient’s disease course. A PCR-based 12-gene signature significantly predicts disease progression in patients with non-muscle invasive bladder cancer.

Proteomic differential display identifies upregulated vinculin as a possible biomarker of pancreatic cancer

Pancreatic cancer (PC) is characterized by rapid tumor spread, and very few patients with PC survive for more than 5 years. It is imperative to discover additional diagnostic biomarkers or specific therapeutic targets in order to improve the treatment of patients with PC. In search for useful biomarkers, we analyzed ten pairs of non-cancerous and cancer tissues from patients with PC by two-dimensional gel electrophoresis (2-DE). Nineteen protein spots showed differential expression on 2-DE gels between the cancer and non-cancerous tissues. Six upregulated protein spots were identified by liquid chromatography-tandem mass spectrometry (LC-MS/MS) as calreticulin, glutathione synthetase, stathmin, vinculin, α-enolase and glyceraldehyde-3-phosphate dehydrogenase. Western blotting demonstrated that vinculin was predominantly expressed in the pancreatic cancer tissues compared with to non-cancerous tissues. Our findings indicate that vinculin may be a clinically useful biomarker of PC.

Colorectal Neoplasia Differentially Expressed (CRNDE), a Novel Gene with Elevated Expression in Colorectal Adenomas and Adenocarcinomas

An uncharacterized gene locus (Chr16:hCG_1815491), now named colorectal neoplasia differentially expressed (gene symbol CRNDE), is activated early in colorectal neoplasia. The locus is unrelated to any known protein-coding gene. Microarray analysis of 454 tissue specimens (discovery) and 68 previously untested specimens (validation) showed elevated expression of CRNDE in >90% of colorectal adenomas and adenocarcinomas. These findings were confirmed and extended by exon microarray studies and RT-PCR assays. CRNDE transcription start sites were identified in CaCo2 and HCT116 cells by 5'-RACE. The major transcript isoforms in colorectal cancer (CRC) cell lines and colorectal tissue are CRNDE-a, -b, -d, -e, -f, -h, and -j. Except for CRNDE-d, the known CRNDE splice variants are upregulated in neoplastic colorectal tissue; expression levels for CRNDE-h alone demonstrate a sensitivity of 95% and specificity of 96% for adenoma versus normal tissue. A quantitative RT-PCR assay measuring CRNDE-h RNA levels in plasma was (with a threshold of 2(-ΔCt) = 2.8) positive for 13 of 15 CRC patients (87%) but only 1 of 15 healthy individuals (7%). We conclude that individual CRNDE transcripts show promise as tissue and plasma biomarkers, potentially exhibiting high sensitivity and specificity for colorectal adenomas and cancers.

Transcriptome-wide detection of differentially expressed coding and non-coding transcripts and their clinical significance in prostate cancer

Prostate cancer is a clinically and biologically heterogeneous disease. Deregulation of splice variants has been shown to contribute significantly to this complexity. High-throughput technologies such as oligonucleotide microarrays allow for the detection of transcripts that play a role in disease progression in a transcriptome-wide level. In this study, we use a publicly available dataset of normal adjacent, primary tumor, and metastatic prostate cancer samples (GSE21034) to detect differentially expressed coding and non-coding transcripts between these disease states. To achieve this, we focus on transcript-specific probe selection regions, that is, those probe sets that correspond unambiguously to a single transcript. Based on this, we are able to pinpoint at the transcript-specific level transcripts that are differentially expressed throughout prostate cancer progression. We confirm previously reported cases and find novel transcripts for which no prior implication in prostate cancer progression has been made. Furthermore, we show that transcript-specific differential expression has unique prognostic potential and provides a clinically significant source of biomarker signatures for prostate cancer risk stratification. The results presented here serve as a catalog of differentially expressed transcript-specific markers throughout prostate cancer progression that can be used as basis for further development and translation into the clinic.

SplicerEX: a tool for the automated detection and classification of mRNA changes from conventional and splice-sensitive microarray expression data

The key postulate that one gene encodes one protein has been overhauled with the discovery that one gene can generate multiple RNA transcripts through alternative mRNA processing. In this study, we describe SplicerEX, a novel and uniquely motivated algorithm designed for experimental biologists that (1) detects widespread changes in mRNA isoforms from both conventional and splice sensitive microarray data, (2) automatically categorizes mechanistic changes in mRNA processing, and (3) mitigates known technological artifacts of exon array-based detection of alternative splicing resulting from 5' and 3' signal attenuation, background detection limits, and saturation of probe set signal intensity. In this study, we used SplicerEX to compare conventional and exon-based Affymetrix microarray data in a model of EBV transformation of primary human B cells. We demonstrated superior detection of 3'-located changes in mRNA processing by the Affymetrix U133 GeneChip relative to the Human Exon Array. SplicerEX-identified exon-level changes in the EBV infection model were confirmed by RT-PCR and revealed a novel set of EBV-regulated mRNA isoform changes in caspases 6, 7, and 8. Finally, SplicerEX as compared with MiDAS analysis of publicly available microarray data provided more efficiently categorized mRNA isoform changes with a significantly higher proportion of hits supported by previously annotated alternative processing events. Therefore, SplicerEX provides an important tool for the biologist interested in studying changes in mRNA isoform usage from conventional or splice-sensitive microarray platforms, especially considering the expansive amount of archival microarray data generated over the past decade. SplicerEX is freely available upon request.

Downregulation of splicing factor SRSF3 induces p53β, an alternatively spliced isoform of p53 that promotes cellular senescence

Most human pre-mRNA transcripts are alternatively spliced, but the significance and fine-tuning of alternative splicing in different biological processes is only starting to be understood. SRSF3 (SRp20) is a member of a highly conserved family of splicing factors that have critical roles in key biological processes, including tumor progression. Here, we show that SRSF3 regulates cellular senescence, a p53-mediated process to suppress tumorigenesis, through TP53 alternative splicing. Downregulation of SRSF3 was observed in normal human fibroblasts undergoing replicative senescence, and was associated with the upregulation of p53β, an alternatively spliced isoform of p53 that promotes p53-mediated senescence. Knockdown of SRSF3 by short interfering RNA (siRNA) in early-passage fibroblasts induced senescence, which was associated with elevated expression of p53β at mRNA and protein levels. Knockdown of p53 partially rescued SRSF3-knockdown-induced senescence, suggesting that SRSF3 acts on p53-mediated cellular senescence. RNA pulldown assays demonstrated that SRSF3 binds to an alternatively spliced exon uniquely included in p53β mRNA through the consensus SRSF3-binding sequences. RNA crosslinking and immunoprecipitation assays (CLIP) also showed that SRSF3 in vivo binds to endogenous p53 pre-mRNA at the region containing the p53β-unique exon. Splicing assays using a transfected TP53 minigene in combination with siRNA knockdown of SRSF3 showed that SRSF3 functions to inhibit the inclusion of the p53β-unique exon in splicing of p53 pre-mRNA. These data suggest that downregulation of SRSF3 represents an endogenous mechanism for cellular senescence that directly regulates the TP53 alternative splicing to generate p53β. This study uncovers the role for general splicing machinery in tumorigenesis, and suggests that SRSF3 is a direct regulator of p53.

The miR-143/-145 cluster regulates plasminogen activator inhibitor-1 in bladder cancer

BACKGROUND

Upregulation of the proto-oncogene plasminogen activator inhibitor-1 (PAI-1) is a common hallmark of various solid tumours, but the mechanisms controlling its expression are not fully understood.

METHODS

We investigate microRNAs (miRNAs) regulating PAI-1 in a panel of normal bladder urothelial biopsies, superficial Ta bladder tumours and invasive T1-T4 tumours using expression microarrays and qRT-PCR. The prognostic implications of PAI-1 deregulation are established by tissue microarray staining of non-muscle-invasive bladder tumours. MicroRNA repression of PAI-1 is assayed by ectopic miRNA expression, argonaute immunoprecipitation and luciferase assays.

RESULTS

We found that the miR-143/-145 cluster is downregulated in all stages of bladder cancer and inversely correlated with PAI-1 expression. Mature miR-143 and miR-145 are coordinately expressed, and both directly target the PAI-1 3'UTR, leading to reduced PAI-1 mRNA and protein levels. Furthermore, we show that PAI-1 and miR-145 levels may serve as useful prognostic markers for non-muscle-invasive bladder tumours for which accurate progressive outcome is currently difficult to predict.

CONCLUSION

This report provides the first evidence for direct miRNA regulation of PAI-1 in bladder cancer. We also demonstrate mRNA co-targeting by a cluster of non-family miRNAs, and suggest miR-145 and PAI-1 as clinically relevant biomarkers in bladder cancer.

Functional implications of structural predictions for alternative splice proteins expressed in Her2/neu-induced breast cancers

Alternative splicing allows a single gene to generate multiple mRNA transcripts, which can be translated into functionally diverse proteins. However, experimentally determined structures of protein splice isoforms are rare, and homology modeling methods are poor at predicting atomic-level structural differences because of high sequence identity. Here we exploit the state-of-the-art structure prediction method I-TASSER to analyze the structural and functional consequences of alternative splicing of proteins differentially expressed in a breast cancer model. We first successfully benchmarked the I-TASSER pipeline for structure modeling of all seven pairs of protein splice isoforms, which are known to have experimentally solved structures. We then modeled three cancer-related variant pairs reported to have opposite functions. In each pair, we observed structural differences in regions where the presence or absence of a motif can directly influence the distinctive functions of the variants. Finally, we applied the method to five splice variants overexpressed in mouse Her2/neu mammary tumor: anxa6, calu, cdc42, ptbp1, and tax1bp3. Despite >75% sequence identity between the variants, structural differences were observed in biologically important regions of these protein pairs. These results demonstrate the feasibility of integrating proteomic analysis with structure-based conformational predictions of differentially expressed alternative splice variants in cancers and other conditions.

Splicing programs and cancer

Numerous studies report splicing alterations in a multitude of cancers by using gene-by-gene analysis. However, understanding of the role of alternative splicing in cancer is now reaching a new level, thanks to the use of novel technologies allowing the analysis of splicing at a large-scale level. Genome-wide analyses of alternative splicing indicate that splicing alterations can affect the products of gene networks involved in key cellular programs. In addition, many splicing variants identified as being misregulated in cancer are expressed in normal tissues. These observations suggest that splicing programs contribute to specific cellular programs that are altered during cancer initiation and progression. Supporting this model, recent studies have identified splicing factors controlling cancer-associated splicing programs. The characterization of splicing programs and their regulation by splicing factors will allow a better understanding of the genetic mechanisms involved in cancer initiation and progression and the development of new therapeutic targets.

Tumor-specific usage of alternative transcription start sites in colorectal cancer identified by genome-wide exon array analysis

BACKGROUND

Approximately half of all human genes use alternative transcription start sites (TSSs) to control mRNA levels and broaden the transcriptional output in healthy tissues. Aberrant expression patterns promoting carcinogenesis, however, may arise from alternative promoter usage.

RESULTS

By profiling 108 colorectal samples using exon arrays, we identified nine genes (TCF12, OSBPL1A, TRAK1, ANK3, CHEK1, UGP2, LMO7, ACSL5, and SCIN) showing tumor-specific alternative TSS usage in both adenoma and cancer samples relative to normal mucosa. Analysis of independent exon array data sets corroborated these findings. Additionally, we confirmed the observed patterns for selected mRNAs using quantitative real-time reverse-transcription PCR. Interestingly, for some of the genes, the tumor-specific TSS usage was not restricted to colorectal cancer. A comprehensive survey of the nine genes in lung, bladder, liver, prostate, gastric, and brain cancer revealed significantly altered mRNA isoform ratios for CHEK1, OSBPL1A, and TCF12 in a subset of these cancer types.To identify the mechanism responsible for the shift in alternative TSS usage, we antagonized the Wnt-signaling pathway in DLD1 and Ls174T colorectal cancer cell lines, which remarkably led to a shift in the preferred TSS for both OSBPL1A and TRAK1. This indicated a regulatory role of the Wnt pathway in selecting TSS, possibly also involving TP53 and SOX9, as their transcription binding sites were enriched in the promoters of the tumor preferred isoforms together with their mRNA levels being increased in tumor samples. Finally, to evaluate the prognostic impact of the altered TSS usage, immunohistochemistry was used to show deregulation of the total protein levels of both TCF12 and OSBPL1A, corresponding to the mRNA levels observed. Furthermore, the level of nuclear TCF12 had a significant correlation to progression free survival in a cohort of 248 stage II colorectal cancer samples.

CONCLUSIONS

Alternative TSS usage in colorectal adenoma and cancer samples has been shown for nine genes, and OSBPL1A and TRAK1 were found to be regulated in vitro by Wnt signaling. TCF12 protein expression was upregulated in cancer samples and correlated with progression free survival.

Low ANXA10 expression is associated with disease aggressiveness in bladder cancer

Background:

Markers for outcome prediction in bladder cancer are urgently needed. We have previously identified a molecular signature for predicting progression in non-muscle-invasive bladder cancer. ANXA10 was one of the markers included in the signature and we now validated the prognostic relevance of ANXA10 at the protein level.

Methods:

We investigated ANXA10 expression by immunohistochemistry using a tissue microarray with 249 Ta and T1 urothelial carcinomas. The expression of ANXA10 was also investigated in an additional set of 97 more advanced tumours. The functional role of ANXA10 in cell lines was investigated by siRNA-mediated ANXA10 knockdown using wound-healing assays, proliferation assays, and ingenuity pathway analysis.

Results:

Low expression of ANXA10 correlated with shorter progression-free survival in patients with stage Ta and T1 tumours (P<0.00001). Furthermore, patients with more advanced tumours and low ANXA10 expression had an unfavourable prognosis (P<0.00001). We found that ANXA10 siRNA transfected cells grew significantly faster compared with control siRNA transfected cells. Furthermore, a wound-healing assay showed that ANXA10 siRNA transfected cells spread along wound edges faster than control transfected cells.

Conclusion:

We conclude that ANXA10 may be a clinical relevant marker for predicting outcome in both early and advanced stages of bladder cancer.

Alternative splicing in oncogenic kinases: from physiological functions to cancer

Among the 518 protein kinases encoded by the human kinome, several of them act as oncoproteins in human cancers. Like other eukaryotic genes, oncogenes encoding protein kinases are frequently subjected to alternative splicing in coding as well as noncoding sequences. In the present paper, we will illustrate how alternative splicing can significantly impact on the physiological functions of oncogenic protein kinases, as demonstrated by mouse genetic model studies. This includes examples of membrane-bound tyrosine kinases receptors (FGFR2, Ret, TrkB, ErbB4, and VEGFR) as well as cytosolic protein kinases (B-Raf). We will further discuss how regular alternative splicing events of these kinases are in some instances implicated in oncogenic processes during tumor progression (FGFR, TrkB, ErbB2, Abl, and AuroraA). Finally, we will present typical examples of aberrant splicing responsible for the deregulation of oncogenic kinases activity in cancers (AuroraB, Jak2, Kit, Met, and Ron).

Attenuation of the beta-catenin/TCF4 complex in colorectal cancer cells induces several growth-suppressive microRNAs that target cancer promoting genes

Aberrant activation of the Wnt signaling pathway is causally involved in the formation of most colorectal cancers (CRCs). Although detailed knowledge exists regarding Wnt-regulated protein-coding genes, much less is known about the possible involvement of non-coding RNAs. Here we used TaqMan Array MicroRNA Cards, capable of detecting 664 unique human microRNAs (miRNAs), to describe changes of the miRNA transcriptome following disruption of beta-catenin/TCF4 activity in DLD1 CRC cells. Most miRNAs appeared to respond independent of host gene regulation and proximal TCF4 chromatin occupancy as inferred from expression microarray and ChIP-chip data. A module of miRNAs induced by abrogated Wnt signaling in vitro was downregulated in two independent series of human primary CRCs (n=76) relative to normal adjacent mucosa (n=34). Several of these miRNAs (miR-145, miR-126, miR-30e-3p and miR-139-5p) markedly inhibited CRC cell growth in vitro when ectopically expressed. By using an integrative approach of proteomics and expression microarrays, we found numerous mRNAs and proteins to be affected by ectopic miR-30e-3p levels. This included HELZ and PIK3C2A that were directly repressed by several miRNA binding sites as confirmed by luciferase reporter assays in combination with mutational analyses. Finally, small interfering RNA-mediated downregulation of PIK3C2A, but not HELZ, was sufficient on its own to restrict CRC cell growth. Collectively, our study demonstrates that multiple miRNAs are upregulated as a consequence of forced attenuation of Wnt signaling in CRC cells, and some of these miRNAs inhibit cell growth with concomitant suppression of several growth-stimulatory cancer-related genes.

Tumor-specific expression and alternative splicing of the COL6A3 gene in pancreatic cancer

BACKGROUND

Pancreatic ductal adenocarcinoma (PDA) is a highly lethal disease; a prominent desmoplastic reaction is a defining characteristic. Fibrillar collagens, such as collagen I and to a lesser extent, collagens III and V, comprise the majority of this stromal fibrosis. Type VI collagen (COL6) forms a microfibrillar network associated with type I collagen fibrils. The expression of COL6 has been linked with inflammation and survival. Importantly, tumor-specific alternative splicing in COL6A3 has been identified in several cancers by genome exon arrays. We evaluated the expression and localization of COL6A3 in PDA and premalignant lesions and explored the presence of alternative splicing events.

METHODS

We analyzed paired PDA-normal (n = 18), intraductal papillary mucinous neoplasms (IPMN; n = 5), pancreatic cystadenoma (n = 5), and 8 PDA cell lines with reverse transcriptase polymerase chain reaction, using unique primers that identify total COL6A3 gene and alternative splicing sites in several of its exons. Western blot analysis and immunohistochemistry were used to analyze the expression levels and localization of COL6A3 protein in the different lesions, and in 2 animal models of PDA.

RESULTS

COL6A3 protein levels were significantly upregulated in 77% of the paired PDA-adjacent tissue examined. COL6A3 was mainly present in the desmoplastic stroma of PDA, with high deposition around the malignant ducts and in between the sites of stromal fatty infiltration. Analysis of the COL6A3 splice variants showed tumor-specific consistent inclusion of exons 3 and 6 in 17 of the 18 (94%) paired PDA-adjacent tissues. Inclusion of exon 4 was exclusively tumor specific, with barely detectable expression in the adjacent tissues. IPMN and pancreatic cystadenomas showed no expression of any of the examined exons. Total COL6A3 mRNA and exon 6 were identified in 6 PDA cell lines, but only 2 cell lines (MIA PACA-2 and ASPC-1) expressed exons 3 and 4. In both the xenograft and transgenic models of PDA, COL6A3 immunoreactivity was present in the stroma and some PDA cells.

CONCLUSION

We have described, for the first time, a dynamic process of tumor-specific alternative splicing in several exons of stromal COL6A3. Alternatively spliced proteins may contribute to the etiology or progression of cancer and may serve as markers for cancer diagnosis. Identification of COL6A3 isoforms as PDA-specific provides the basis for future studies to explore the oncogenic and diagnostic potential of these alternative splicing events.

Repression of KIAA1199 attenuates Wnt-signalling and decreases the proliferation of colon cancer cells

BACKGROUND

The KIAA1199 transcript is upregulated in colon adenomas and downregulated upon β-catenin knockdown.

METHODS

Transcript profiling was performed on >500 colon biopsies, methylation profiling data were compared with transcript data. Immunohistochemistry assessed KIAA1199 protein expression in 270 stage II/III tumours (>3 years follow-up). The effects of stable KIAA1199 knockdown in SW480 cells (three different constructs) were studied using transcriptional profiling, proliferation and protein analysis.

RESULTS

The KIAA1199 transcript was strongly upregulated in 95% of adenocarcinomas. Absent expression in normal mucosa correlated with KIAA1199 promotor methylation. Nuclear and cytoplasmic KIAA1199 protein expression was identified in colon adenocarcinomas and other types of cancers. A subpopulation of patients with tumours strongly expressing KIAA1199 in the nucleus showed a better outcome with regard to recurrence as lung or liver metastases. The KIAA1199 knockdown affected the cell cycle and the Wnt-signalling pathway. Reduced cellular proliferation and decreased KI67, phosphorylated retinoblastoma, β-catenin and ASCL2 protein expression supported these findings. Eighteen Wnt-signalling genes differentially expressed upon KIAA1199 knockdown correlated with the KIAA1199 expression profile in clinical specimens.

CONCLUSION

The KIAA1199 knockdown attenuates the effects of the Wnt/β-catenin signalling and it may thus be regarded as a regulatory part of this pathway.

Alternative splicing and tumor progression

Alternative splicing is a key molecular mechanism for increasing the functional diversity of the eukaryotic proteomes. A large body of experimental data implicates aberrant splicing in various human diseases, including cancer. Both mutations in cis-acting splicing elements and alterations in the expression and/or activity of splicing regulatory factors drastically affect the splicing profile of many cancer-associated genes. In addition, the splicing profile of several cancer-associated genes is altered in particular types of cancer arguing for a direct role of specific splicing isoforms in tumor progression. Deciphering the mechanisms underlying aberrant splicing in cancer may prove crucial to understand how splicing machinery is controlled and integrated with other cellular processes, in particular transcription and signaling pathways. Moreover, the characterization of splicing deregulation in cancer will lead to a better comprehension of malignant transformation. Cancer-associated alternative splicing variants may be new tools for the diagnosis and classification of cancers and could be the targets for innovative therapeutical interventions based on highly selective splicing correction approaches.

Downregulation of zinc finger protein 132 in prostate cancer is associated with aberrant promoter hypermethylation and poor prognosis

This study investigates the expression and biomarker potential of zinc finger protein 132 (ZNF132) in prostate cancer (PC) by transcriptional profiling and immunohistochemical analysis of tissue microarrays, including tumor specimens from 615 radical prostatectomy (RP) patients and 199 conservatively treated patients. Primary clinical endpoints were time to PSA recurrence and cancer-specific death, respectively. Compared to normal prostate epithelial cells from men without PC, ZNF132 transcript levels were significantly reduced in PC cells from patients with localized PC and further downregulated in metastatic PC. Likewise, ZNF132 protein expression was significantly lower in primary tumors from patients with metastatic compared to localized PC and further reduced in castrate-refractory PC, indicating that ZNF132 downregulation correlates with disease progression. Reduced ZNF132 immunoreactivity was significantly associated with high Gleason score and advanced T stage in both PC patient cohorts. By univariate analysis, no/weak ZNF132 staining was a significant adverse predictor of PSA recurrence after RP (p = 0.024) and cancer-specific death following conservative treatment (p = 0.009). In multivariate models, however, ZNF132 did not add significant independent value to established prognostic factors. Finally, bisulfite sequencing revealed frequent promoter hypermethylation of ZNF132 in both PC cell lines and PC tissue samples, indicating that ZNF132 is epigenetically silenced in PC. In summary, our results show that downregulation of ZNF132 is associated with aggressive PC and furthermore identify ZNF132 as a new candidate methylation marker for PC.

Exon array analysis reveals neuroblastoma tumors have distinct alternative splicing patterns according to stage and MYCN amplification status

BACKGROUND

Neuroblastoma (NB) tumors are well known for their pronounced clinical and molecular heterogeneity. The global gene expression and DNA copy number alterations have been shown to have profound differences in tumors of low or high stage and those with or without MYCN amplification. RNA splicing is an important regulatory mechanism of gene expression, and differential RNA splicing may be associated with the clinical behavior of a tumor.

METHODS

In this study, we used exon array profiling to investigate global alternative splicing pattern of 47 neuroblastoma samples in stage 1 and stage 4 with normal or amplified MYCN copy number (stage 1-, 4- and 4+). The ratio of exon-level expression to gene-level expression was used to detect alternative splicing events, while the gene-level expression was applied to characterize whole gene expression change.

RESULTS

Principal component analysis (PCA) demonstrated distinct splicing pattern in three groups of samples. Pairwise comparison identified genes with splicing changes and/or whole gene expression changes in high stage tumors. In stage 4- compared with stage 1- tumors, alternatively spliced candidate genes had little overlap with genes showing whole gene expression changes, and most of them were involved in different biological processes. In contrast, a larger number of genes exhibited either exon-level splicing, gene-level expression or both changes in stage 4+ versus stage 1- tumors. Those biological processes involved in stage 4- tumors were disrupted to a greater extent by both splicing and transcription regulations in stage 4+ tumors.

CONCLUSIONS

Our results demonstrated a significant role of alternative splicing in high stage neuroblastoma, and suggested a MYCN-associated splicing regulation pathway in stage 4+ tumors. The identification of differentially spliced genes and pathways in neuroblastoma tumors of different stages and molecular subtypes may be important to the understanding of cancer biology and the discovery of diagnostic markers or therapeutic targets in neuroblastoma.

Melanoma spheroids grown under neural crest cell conditions are highly plastic migratory/invasive tumor cells endowed with immunomodulator function

Background

The aggressiveness of melanoma tumors is likely to rely on their well-recognized heterogeneity and plasticity. Melanoma comprises multi-subpopulations of cancer cells some of which may possess stem cell-like properties. Although useful, the sphere-formation assay to identify stem cell-like or tumor initiating cell subpopulations in melanoma has been challenged, and it is unclear if this model can predict a functional phenotype associated with aggressive tumor cells.

Methodology/Principal Findings

We analyzed the molecular and functional phenotypes of melanoma spheroids formed in neural crest cell medium. Whether from metastatic or advanced primary tumors, spheroid cells expressed melanoma-associated markers. They displayed higher capacity to differentiate along mesenchymal lineages and enhanced expression of SOX2, NANOG, KLF4, and/or OCT4 transcription factors, but not enhanced self-renewal or tumorigenicity when compared to their adherent counterparts. Gene expression profiling attributed a neural crest cell signature to these spheroids and indicated that a migratory/invasive and immune-function modulating program could be associated with these cells. In vitro assays confirmed that spheroids display enhanced migratory/invasive capacities. In immune activation assays, spheroid cells elicited a poorer allogenic response from immune cells and inhibited mitogen-dependent T cells activation and proliferation more efficiently than their adherent counterparts. Our findings reveal a novel immune-modulator function of melanoma spheroids and suggest specific roles for spheroids in invasion and in evasion of antitumor immunity.

Conclusion/Significance

The association of a more plastic, invasive and evasive, thus a more aggressive tumor phenotype with melanoma spheroids reveals a previously unrecognized aspect of tumor cells expanded as spheroid cultures. While of limited efficiency for melanoma initiating cell identification, our melanoma spheroid model predicted aggressive phenotype and suggested that aggressiveness and heterogeneity of melanoma tumors can be supported by subpopulations other than cancer stem cells. Therefore, it could be constructive to investigate melanoma aggressiveness, relevant to patients and clinical transferability.

Characterization of alternative spliceoforms and the RNA splicing machinery in pancreatic cancer

OBJECTIVES AND METHODS

Alternative splicing provides proteomic diversity that can have profound effects. The extent, pattern, and roles of alternative splicing in pancreatic cancer have not been systematically investigated. We have utilized a spliceoform-specific microarray and polymerase chain reaction to evaluate all known splice variants in human pancreatic cancer cell lines representing a spectrum of differentiation, from near-normal HPDE6 to Capan-1 and poorly differentiated MiaPaCa2 cells. Validation of altered spliceoforms was verified in primary cancer specimens and normal pancreatic ductal cells. In addition, expression of 92 spliceosomal genes was examined to better understand the mechanism for observed differences in mRNA splicing.

RESULTS

A statistically significant reduction in alternative splicing was found in the pancreatic cancer cell lines compared with HPDE6 cells. Many splice variants identified in Capan-1 and MiaPaCa2 cells were observed in grades 3 and 4 tumors. Analysis of genes encoding spliceosomal proteins revealed that 28 of 92 genes had significantly decreased expression in cancer compared with normal pancreas.

CONCLUSIONS

Pancreatic cancer has reduced alternative splicing diversity compared with normal pancreas. This is demonstrated in both cell lines and primary tumors, with the loss in splicing diversity correlated with relative reduction in expression of spliceosomal genes.

Alternative splicing in bone following mechanical loading

It is estimated that more than 90% of human genes express multiple mRNA transcripts due to alternative splicing. Consequently, the proteins produced by different splice variants will likely have different functions and expression levels. Several genes with splice variants are known in bone, with functions that affect osteoblast function and bone formation. The primary goal of this study was to evaluate the extent of alternative splicing in a bone subjected to mechanical loading and subsequent bone formation. We used the rat forelimb loading model, in which the right forelimb was loaded axially for 3 min, while the left forearm served as a non-loaded control. Animals were subjected to loading sessions every day, with 24 h between sessions. Ulnae were sampled at 11 time points, from 4 h to 32days after beginning loading. RNA was isolated and mRNA abundance was measured at each time point using Affymetrix exon arrays (GeneChip® Rat Exon 1.0 ST Arrays). An ANOVA model was used to identify potential alternatively spliced genes across the time course, and five alternatively spliced genes were validated with qPCR: Akap12, Fn1, Pcolce, Sfrp4, and Tpm1. The number of alternatively spliced genes varied with time, ranging from a low of 68 at 12h to a high of 992 at 16d. We identified genes across the time course that encoded proteins with known functions in bone formation, including collagens, matrix proteins, and components of the Wnt/β-catenin and TGF-β signaling pathways. We also identified alternatively spliced genes encoding cytokines, ion channels, muscle-related genes, and solute carriers that do not have a known function in bone formation and represent potentially novel findings. In addition, a functional characterization was performed to categorize the global functions of the alternatively spliced genes in our data set. In conclusion, mechanical loading induces alternative splicing in bone, which may play an important role in the response of bone to mechanical loading.

New insights in gene-derived therapy: the example of Duchenne muscular dystrophy

The two therapeutic approaches currently most advanced in clinical trials for Duchenne muscular dystrophy are antisense-mediated exon skipping and forced read-through of premature stop codons. Interestingly, these approaches target the gene product rather than the gene itself. This review will explain the rationale and current state of affairs of these approaches and will then discuss how these gene-derived therapies might also be applicable to other diseases.

Detection of clinically relevant exonic copy-number changes by array CGH

Array comparative genomic hybridization (aCGH) is a powerful tool for the molecular elucidation and diagnosis of disorders resulting from genomic copy-number variation (CNV). However, intragenic deletions or duplications--those including genomic intervals of a size smaller than a gene--have remained beyond the detection limit of most clinical aCGH analyses. Increasing array probe number improves genomic resolution, although higher cost may limit implementation, and enhanced detection of benign CNV can confound clinical interpretation. We designed an array with exonic coverage of selected disease and candidate genes and used it clinically to identify losses or gains throughout the genome involving at least one exon and as small as several hundred base pairs in size. In some patients, the detected copy-number change occurs within a gene known to be causative of the observed clinical phenotype, demonstrating the ability of this array to detect clinically relevant CNVs with subkilobase resolution. In summary, we demonstrate the utility of a custom-designed, exon-targeted oligonucleotide array to detect intragenic copy-number changes in patients with various clinical phenotypes.

Exon array analysis using re-defined probe sets results in reliable identification of alternatively spliced genes in non-small cell lung cancer

Background

Treatment of non-small cell lung cancer with novel targeted therapies is a major unmet clinical need. Alternative splicing is a mechanism which generates diverse protein products and is of functional relevance in cancer.

Results

In this study, a genome-wide analysis of the alteration of splicing patterns between lung cancer and normal lung tissue was performed. We generated an exon array data set derived from matched pairs of lung cancer and normal lung tissue including both the adenocarcinoma and the squamous cell carcinoma subtypes. An enhanced workflow was developed to reliably detect differential splicing in an exon array data set. In total, 330 genes were found to be differentially spliced in non-small cell lung cancer compared to normal lung tissue. Microarray findings were validated with independent laboratory methods for CLSTN1, FN1, KIAA1217, MYO18A, NCOR2, NUMB, SLK, SYNE2, TPM1, (in total, 10 events) and ADD3, which was analysed in depth. We achieved a high validation rate of 69%. Evidence was found that the activity of FOX2, the splicing factor shown to cause cancer-specific splicing patterns in breast and ovarian cancer, is not altered at the transcript level in several cancer types including lung cancer.

Conclusions

This study demonstrates how alternatively spliced genes can reliably be identified in a cancer data set. Our findings underline that key processes of cancer progression in NSCLC are affected by alternative splicing, which can be exploited in the search for novel targeted therapies.

Specific variants in the MLH1 gene region may drive DNA methylation, loss of protein expression, and MSI-H colorectal cancer

BACKGROUND

We previously identified an association between a mismatch repair gene, MLH1, promoter SNP (rs1800734) and microsatellite unstable (MSI-H) colorectal cancers (CRCs) in two samples. The current study expanded on this finding as we explored the genetic basis of DNA methylation in this region of chromosome 3. We hypothesized that specific polymorphisms in the MLH1 gene region predispose it to DNA methylation, resulting in the loss of MLH1 gene expression, mismatch-repair function, and consequently to genome-wide microsatellite instability.

METHODOLOGY/PRINCIPAL FINDINGS

We first tested our hypothesis in one sample from Ontario (901 cases, 1,097 controls) and replicated major findings in two additional samples from Newfoundland and Labrador (479 cases, 336 controls) and from Seattle (591 cases, 629 controls). Logistic regression was used to test for association between SNPs in the region of MLH1 and CRC, MSI-H CRC, MLH1 gene expression in CRC, and DNA methylation in CRC. The association between rs1800734 and MSI-H CRCs, previously reported in Ontario and Newfoundland, was replicated in the Seattle sample. Two additional SNPs, in strong linkage disequilibrium with rs1800734, showed strong associations with MLH1 promoter methylation, loss of MLH1 protein, and MSI-H CRC in all three samples. The logistic regression model of MSI-H CRC that included MLH1-promoter-methylation status and MLH1 immunohistochemistry status fit most parsimoniously in all three samples combined. When rs1800734 was added to this model, its effect was not statistically significant (P-value  = 0.72 vs. 2.3×10(-4) when the SNP was examined alone).

CONCLUSIONS/SIGNIFICANCE

The observed association of rs1800734 with MSI-H CRC occurs through its effect on the MLH1 promoter methylation, MLH1 IHC deficiency, or both.

Alternative splicing of SLC39A14 in colorectal cancer is regulated by the Wnt pathway

Alternative splicing is a crucial step in the generation of protein diversity and its misregulation is observed in many human cancer types. By analyzing 143 colorectal samples using exon arrays, SLC39A14, a divalent cation transporter, was identified as being aberrantly spliced in tumor samples. SLC39A14 contains two mutually exclusive exons 4A and 4B and the exon 4A/4B ratio was significantly altered in adenomas (p = 3.6 × 10(-10)) and cancers (p = 9.4 × 10(-11)), independent of microsatellite stability status. The findings were validated in independent exon array data sets and by quantitative real-time reverse-transcription PCR (qRT-PCR). Aberrant Wnt signaling is a hallmark of colorectal tumorigenesis and is characterized by nuclear β-catenin. Experimental inactivation of Wnt signaling in DLD1 and Ls174T cells by knockdown of β-catenin or overexpression of dominant negative TCFs (TCF1 and TCF4) altered the 4A/4B ratio, indicating that SLC39A14 splicing is regulated by the Wnt pathway. An altered 4A/4B ratio was also observed in gastric and lung cancer where Wnt signaling is also known to be aberrantly activated. The splicing factor SRSF1 and its regulator, the kinase SRPK1, were found to be deregulated upon Wnt inactivation in colorectal carcinoma cells. SRPK1 was also found up-regulated in both adenoma samples (p = 1.5 × 10(-5)) and cancer samples (p = 5 × 10(-4)). In silico splicing factor binding analysis predicted SRSF1 to bind predominantly to the cancer associated exon 4B, hence, it was hypothesized that SRPK1 activates SRSF1 through phosphorylation, followed by SRSF1 binding to exon 4B and regulation of SLC39A14 splicing. Indeed, siRNA-mediated knockdown of SRPK1 and SRSF1 in DLD1 and SW480 colorectal cancer cells led to a change in the 4A/4B isoform ratio, supporting a role of these factors in the regulation of SLC39A14 splicing. In conclusion, alternative splicing of SLC39A14 was identified in colorectal tumors and found to be regulated by the Wnt pathway, most likely through regulation of SRPK1 and SRSF1.

Alternative splicing of SLC39A14 in colorectal cancer is regulated by the Wnt pathway

Alternative splicing is a crucial step in the generation of protein diversity and its misregulation is observed in many human cancer types. By analyzing 143 colorectal samples using exon arrays, SLC39A14, a divalent cation transporter, was identified as being aberrantly spliced in tumor samples. SLC39A14 contains two mutually exclusive exons 4A and 4B and the exon 4A/4B ratio was significantly altered in adenomas (p = 3.6 × 10(-10)) and cancers (p = 9.4 × 10(-11)), independent of microsatellite stability status. The findings were validated in independent exon array data sets and by quantitative real-time reverse-transcription PCR (qRT-PCR). Aberrant Wnt signaling is a hallmark of colorectal tumorigenesis and is characterized by nuclear β-catenin. Experimental inactivation of Wnt signaling in DLD1 and Ls174T cells by knockdown of β-catenin or overexpression of dominant negative TCFs (TCF1 and TCF4) altered the 4A/4B ratio, indicating that SLC39A14 splicing is regulated by the Wnt pathway. An altered 4A/4B ratio was also observed in gastric and lung cancer where Wnt signaling is also known to be aberrantly activated. The splicing factor SRSF1 and its regulator, the kinase SRPK1, were found to be deregulated upon Wnt inactivation in colorectal carcinoma cells. SRPK1 was also found up-regulated in both adenoma samples (p = 1.5 × 10(-5)) and cancer samples (p = 5 × 10(-4)). In silico splicing factor binding analysis predicted SRSF1 to bind predominantly to the cancer associated exon 4B, hence, it was hypothesized that SRPK1 activates SRSF1 through phosphorylation, followed by SRSF1 binding to exon 4B and regulation of SLC39A14 splicing. Indeed, siRNA-mediated knockdown of SRPK1 and SRSF1 in DLD1 and SW480 colorectal cancer cells led to a change in the 4A/4B isoform ratio, supporting a role of these factors in the regulation of SLC39A14 splicing. In conclusion, alternative splicing of SLC39A14 was identified in colorectal tumors and found to be regulated by the Wnt pathway, most likely through regulation of SRPK1 and SRSF1.

Identification of tumor-associated cassette exons in human cancer through EST-based computational prediction and experimental validation

BACKGROUND

Many evidences report that alternative splicing, the mechanism which produces mRNAs and proteins with different structures and functions from the same gene, is altered in cancer cells. Thus, the identification and characterization of cancer-specific splice variants may give large impulse to the discovery of novel diagnostic and prognostic tumour biomarkers, as well as of new targets for more selective and effective therapies.

RESULTS

We present here a genome-wide analysis of the alternative splicing pattern of human genes through a computational analysis of normal and cancer-specific ESTs from seventeen anatomical groups, using data available in AspicDB, a database resource for the analysis of alternative splicing in human. By using a statistical methodology, normal and cancer-specific genes, splice sites and cassette exons were predicted in silico. The condition association of some of the novel normal/tumoral cassette exons was experimentally verified by RT-qPCR assays in the same anatomical system where they were predicted. Remarkably, the presence in vivo of the predicted alternative transcripts, specific for the nervous system, was confirmed in patients affected by glioblastoma.

CONCLUSION

This study presents a novel computational methodology for the identification of tumor-associated transcript variants to be used as cancer molecular biomarkers, provides its experimental validation, and reports specific biomarkers for glioblastoma.

Position-dependent alternative splicing activity revealed by global profiling of alternative splicing events regulated by PTB

To gain global insights into the role of the well-known repressive splicing regulator PTB, we analyzed the consequences of PTB knockdown in HeLa cells using high-density oligonucleotide splice-sensitive microarrays. The major class of identified PTB-regulated splicing event was PTB-repressed cassette exons, but there was also a substantial number of PTB-activated splicing events. PTB-repressed and PTB-activated exons showed a distinct arrangement of motifs with pyrimidine-rich motif enrichment within and upstream of repressed exons but downstream of activated exons. The N-terminal half of PTB was sufficient to activate splicing when recruited downstream of a PTB-activated exon. Moreover, insertion of an upstream pyrimidine tract was sufficient to convert a PTB-activated exon to a PTB-repressed exon. Our results show that PTB, an archetypal splicing repressor, has variable splicing activity that predictably depends upon its binding location with respect to target exons.

Exon-level transcriptome profiling in murine breast cancer reveals splicing changes specific to tumors with different metastatic abilities

BACKGROUND

Breast cancer is the second most frequent type of cancer affecting women. We are increasingly aware that changes in mRNA splicing are associated with various characteristics of cancer. The most deadly aspect of cancer is metastasis, the process by which cancer spreads from the primary tumor to distant organs. However, little is known specifically about the involvement of alternative splicing in the formation of macroscopic metastases. Our study investigates transcript isoform changes that characterize tumors of different abilities to form growing metastases.

METHODS AND FINDINGS

To identify alternative splicing events (ASEs) that are associated with the fully metastatic phenotype in breast cancer, we used Affymetrix Exon Microarrays to profile mRNA isoform variations genome-wide in weakly metastatic (168FARN and 4T07) and highly metastatic (4T1) mammary carcinomas. Statistical analysis identified significant expression changes in 7606 out of 155,994 (4%) exons and in 1725 out of 189,460 (1%) intronic regions, which affect 2623 out of 16,654 (16%) genes. These changes correspond to putative alternative isoforms-several of which are novel-that are differentially expressed between tumors of varying metastatic phenotypes. Gene pathway analysis showed that 1224 of genes expressing alternative isoforms were involved in cell growth, cell interactions, cell proliferation, cell migration and cell death and have been previously linked to cancers and genetic disorders. We chose ten predicted splice variants for RT-PCR validation, eight of which were successfully confirmed (MED24, MFI2, SRRT, CD44, CLK1 and HNRNPH1). These include three novel intron retentions in CD44, a gene in which isoform variations have been previously associated with the metastasis of several cancers.

CONCLUSION

Our findings reveal that various genes are differently spliced and/or expressed in association with the metastatic phenotype of tumor cells. Identification of metastasis-specific isoforms may contribute to the development of improved breast cancer stage identification and targeted therapies.

The dietary isothiocyanate sulforaphane modulates gene expression and alternative gene splicing in a PTEN null preclinical murine model of prostate cancer

Background

Dietary or therapeutic interventions to counteract the loss of PTEN expression could contribute to the prevention of prostate carcinogenesis or reduce the rate of cancer progression. In this study, we investigate the interaction between sulforaphane, a dietary isothiocyanate derived from broccoli, PTEN expression and gene expression in pre malignant prostate tissue.

Results

We initially describe heterogeneity in expression of PTEN in non-malignant prostate tissue of men deemed to be at risk of prostate cancer. We subsequently use the mouse prostate-specific PTEN deletion model, to show that sulforaphane suppresses transcriptional changes induced by PTEN deletion and induces additional changes in gene expression associated with cell cycle arrest and apoptosis in PTEN null tissue, but has no effect on transcription in wild type tissue. Comparative analyses of changes in gene expression in mouse and human prostate tissue indicate that similar changes can be induced in humans with a broccoli-rich diet. Global analyses of exon expression demonstrated that sulforaphane interacts with PTEN deletion to modulate alternative gene splicing, illustrated through a more detailed analysis of DMBT1 splicing.

Conclusion

To our knowledge, this is the first report of how diet may perturb changes in transcription induced by PTEN deletion, and the effects of diet on global patterns of alternative gene splicing. The study exemplifies the complex interaction between diet, genotype and gene expression, and the multiple modes of action of small bioactive dietary components.

Development of a novel splice array platform and its application in the identification of alternative splice variants in lung cancer

Background

Microarrays strategies, which allow for the characterization of thousands of alternative splice forms in a single test, can be applied to identify differential alternative splicing events. In this study, a novel splice array approach was developed, including the design of a high-density oligonucleotide array, a labeling procedure, and an algorithm to identify splice events.

Results

The array consisted of exon probes and thermodynamically balanced junction probes. Suboptimal probes were tagged and considered in the final analysis. An unbiased labeling protocol was developed using random primers. The algorithm used to distinguish changes in expression from changes in splicing was calibrated using internal non-spliced control sequences. The performance of this splice array was validated with artificial constructs for CDC6, VEGF, and PCBP4 isoforms. The platform was then applied to the analysis of differential splice forms in lung cancer samples compared to matched normal lung tissue. Overexpression of splice isoforms was identified for genes encoding CEACAM1, FHL-1, MLPH, and SUSD2. None of these splicing isoforms had been previously associated with lung cancer.

Conclusions

This methodology enables the detection of alternative splicing events in complex biological samples, providing a powerful tool to identify novel diagnostic and prognostic biomarkers for cancer and other pathologies.