Global analysis of aberrant pre-mRNA splicing in glioblastoma using exon expression arrays

To RNA-binding and beyond: Emerging facets of the role of Rbfox proteins in development and disease

The RNA-binding Fox-1 homologue (Rbfox) proteins represent an ancient family of splicing factors, conserved through evolution. All members share an RNA recognition motif (RRM), and a particular affinity for the GCAUG signature in target RNA molecules. The role of Rbfox, as a splice factor, deciding the tissue-specific inclusion/exclusion of an exon, depending on its binding position on the flanking introns, is well known. Rbfox often acts in concert with other splicing factors, and forms splicing regulatory networks. Apart from this canonical role, recent studies show that Rbfox can also function as a transcription co-factor, and affects mRNA stability and translation. The repertoire of Rbfox targets is vast, including genes involved in the development of tissue lineages, such as neurogenesis, myogenesis, and erythropoeiesis, and molecular processes, including cytoskeletal dynamics, and calcium handling. A second layer of complexity is added by the fact that Rbfox expression itself is regulated by multiple mechanisms, and, in vertebrates, exhibits tissue-specific expression. The optimum dosage of Rbfox is critical, and its misexpression is etiological to various disease conditions. In this review, we discuss the contextual roles played by Rbfox as a tissue-specific regulator for the expression of many important genes with diverse functions, through the lens of the emerging data which highlights its involvement in many human diseases. Furthermore, we explore the mechanistic details provided by studies in model organisms, with emphasis on the work with Drosophila. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA Turnover and Surveillance > Regulation of RNA Stability RNA Processing > Splicing Regulation/Alternative Splicing.

A Novel Splice Variant of BCAS1 Inhibits β-Arrestin 2 to Promote the Proliferation and Migration of Glioblastoma Cells, and This Effect Was Blocked by Maackiain

Brain-enriched myelin-associated protein 1 (BCAS1) is frequently highly expressed in human cancer, but its detailed function is unclear. Here, we identified a novel splice variant of the BCAS1 gene in glioblastoma multiforme (GBM) named BCAS1-SV1. The expression of BCAS1-SV1 was weak in heathy brain cells but high in GBM cell lines. The overexpression of BCAS1-SV1 significantly increased the proliferation and migration of GBM cells, whereas the RNA-interference-mediated knockdown of BCAS1-SV1 reduced proliferation and migration. Moreover, using a yeast-two hybrid assay, immunoprecipitation, and immunofluorescence staining, we confirmed that β-arrestin 2 is an interaction partner of BCAS1-SV1 but not BCAS1. The downregulation of β-arrestin 2 directly enhanced the malignancy of GBM and abrogated the effects of BCAS1-SV1 on GBM cells. Finally, we used a yeast two-hybrid-based growth assay to identify that maackiain (MK) is a potential inhibitor of the interaction between BCAS1-SV1 and β-arrestin 2. MK treatment lessened the proliferation and migration of GBM cells and prolonged the lifespan of tumor-bearing mice in subcutaneous xenograft and intracranial U87-luc xenograft models. This study provides the first evidence that the gain-of-function BCAS1-SV1 splice variant promotes the development of GBM by suppressing the β-arrestin 2 pathway and opens up a new therapeutic perspective in GBM.

Identifying the Predictive Role of Oxidative Stress Genes in the Prognosis of Glioma Patients

Background

Gliomas are primary aggressive brain tumors with poor prognoses. Oxidative stress plays a crucial role in the tumorigenesis and drug resistance of gliomas. The aim of the present study was to use integrated bioinformatics analyses to evaluate the prognostic value of oxidative stress-related genes (OSRGs) in glioma.

Material/Methods

Disease- and prognosis-associated OSRGs were identified using microarray and clinical data from the Chinese Glioma Genome Atlas database. Functional enrichment, gene-gene interaction, protein-protein interaction, and survival analyses were performed in screened OSRGs. The protein expression was validated by the Human Protein Atlas database. A risk score model was constructed and verified through Cox regression, receiver operating characteristic curve, principal component, and stratified analyses. The Cancer Genome Atlas (TCGA) database was used for external validation. A nomogram was constructed to facilitate the clinical application.

Results

Twenty-one disease-associated and 14 prognosis-associated OSRGs were identified. Enrichment analyses indicated that these signature OSRGs were involved in tumorigenesis and drug resistance of glioma. The risk score model demonstrated a significant difference in overall survival between the high- and low-risk groups. The area under the curve and hazard ratio (1.296) revealed the independent prognostic value of the model. The model exhibited good predictive efficacy in the TCGA cohort. A clinical nomogram was constructed to calculate survival rates in glioma patients at 1, 3, and 5 years.

Conclusions

Our comprehensive study indicated that OSRGs were valuable for prognosis prediction in glioma, which provides a novel insight into the relationship between oxidative stress and glioma and a potential therapeutic strategy for glioma patients.

TPD52L2 Is a Prognostic Biomarker and Correlated With Immune Infiltration in Lung Adenocarcinoma

Tumor protein D52-like 2 (TPD52L2) belongs to the members of the TPD52 family. TPD52L2 was reported to regulate proliferation and apoptosis in cancer cells. However, its role in lung adenocarcinoma (LUAD) was uncertain. We evaluated the expression, methylation, copy number alteration, and prognostic significance of TPD52L2 using RNA-seq data from The Cancer Genome Atlas (TCGA). Enrichment analysis of TPD52L2 was conducted using the R package “clusterProfiler.” We further assessed the association between TPD52L2 and immune cell infiltration level, immunosuppressive genes, and tumor mutational burden (TMB). The difference of gene mutant frequency in high- and low-TPD52L2 groups was also analyzed. The results showed that TPD52L2 was over-expressed and predicted worse survival status in LUAD. We also found that TPD52L2 expression was positively associated with the infiltration levels of immunosuppressive cells, such as regulatory T cells (Tregs) and tumor-associated macrophages (TAMs), and negatively correlated with immune killer cells, such as CD8⁺ T and NK cells in pan-cancer, including LUAD. In addition, TPD52L2 expression was associated with immunosuppressive genes and TMB. High expression of TPD52L2 was with more mutant frequency of TP53. In summary, our results show that TPD52L2 is an oncogene and a potential prognostic biomarker in LUAD. High TPD52L2 expression is a possible indicator of immune infiltration and associated with tumor immunosuppressive status in LUAD.

Genome-wide profiling of alternative splicing in glioblastoma and their clinical value

Background

Alternative splicing (AS), one of the main post-transcriptional biological regulation mechanisms, plays a key role in the progression of glioblastoma (GBM). Systematic AS profiling in GBM is limited and urgently needed.

Methods

TCGA SpliceSeq data and the corresponding clinical data were downloaded from the TCGA data portal. Survival-related AS events were identified through Kaplan–Meier survival analysis and univariate Cox analysis. Then, splicing correlation network was constructed based on these AS events and associated splicing factors. LASSO regression followed by multivariate Cox analysis was performed to validate independent AS biomarkers and to construct a risk prediction model. Enrichment analysis was subsequently conducted to explore potential signaling pathways of these AS events.

Results

A total of 132 TCGA GBM samples and 45,610 AS events were included in our study, among which 416 survival-related AS events were identified. An AS correlation network, including 54 AS events and 94 splicing factors, was constructed, and further functional enrichment was performed. Moreover, the novel risk prediction model we constructed displayed moderate performance (the area under the curves were > 0.7) at both one, two and three years.

Conclusions

Survival-related AS events may be vital factors of both biological function and prognosis. Our findings in this study can deepen the understanding of the complicated mechanisms of AS in GBM and provide novel insights for further study. Moreover, our risk prediction model is ready for preliminary clinical applications. Further verification is required.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08681-z.

Tumor protein D52 is upregulated in oral squamous carcinoma cells under hypoxia in a hypoxia-inducible-factor-independent manner and is involved in cell death resistance

Background

Tumor protein D52 (TPD52) reportedly plays an important role in the proliferation and metastasis of various cancer cells, including oral squamous cell carcinoma (OSCC) cells, and is expressed strongly at the center of the tumor, where the microenvironment is hypoxic. Thus, the present study investigated the roles of TPD52 in the survival and death of OSCC cells under hypoxia, and the relationship with hypoxia-inducible factor (HIF). We examined the expression of TPD52 in OSCC cells under hypoxic conditions and analyzed the effects of HIF on the modulation of TPD52 expression. Finally, the combinational effects of TPD52 knockdown and HIF inhibition were investigated both in vitro and in vivo.

Results

The mRNA and protein levels of TPD52 increased in OSCC cells under hypoxia. However, the increase was independent of HIF transcription. Importantly, the observation was due to upregulation of mRNA stability by binding of mRNA to T-cell intercellular antigen (TIA) 1 and TIA-related protein (TIAR). Simultaneous knockdown of TPD52 and inhibition of HIF significantly reduced cell viability. In addition, the in vivo tumor-xenograft experiments showed that TPD52 acts as an autophagy inhibitor caused by a decrease in p62.

Conclusions

This study showed that the expression of TPD52 increases in OSCC cells under hypoxia in a HIF-independent manner and plays an important role in the proliferation and survival of the cells in concordance with HIF, suggesting that novel cancer therapeutics might be led by TPD52 suppression.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13578-021-00634-0.

RBFOX1 Regulates the Permeability of the Blood-Tumor Barrier via the LINC00673/MAFF Pathway

The blood-tumor barrier limits the delivery of therapeutic drugs to brain tumor tissues. Selectively opening the blood-tumor barrier is considered crucial for effective chemotherapy of glioma. RNA-binding proteins have emerged as crucial regulators in various biologic processes. This study found that RNA-binding Fox-1 homolog 1 (RBFOX1) was downregulated in glioma vascular endothelial cells derived from glioma tissues, and in glioma endothelial cells obtained by co-culturing endothelial cells with glioma cells. Overexpression of RBFOX1 impaired the integrity of the blood-tumor barrier and increased its permeability. Additionally, RBFOX1 overexpression decreased the expression of tight junction proteins ZO-1, occludin, and claudin-5. Subsequent analysis of the mechanism indicated that the overexpression of RBFOX1 increased musculoaponeurotic fibrosarcoma protein basic leucine zipper [bZIP] transcription factor F (MAFF) expression by downregulating LINC00673, which stabilized MAFF messenger RNA (mRNA) through Staufen1-mediated mRNA decay. Moreover, MAFF could bind to the promoter region and inhibit the promoter activities of ZO-1, occludin, and claudin-5, which reduced its expression. The combination of RBFOX1 upregulation and LINC00673 downregulation promoted doxorubicin delivery across the blood-tumor barrier, resulting in apoptosis of glioma cells. In conclusion, this study indicated that overexpression of RBFOX1 increased blood-tumor barrier permeability through the LINC00673/MAFF pathway, which might provide a new useful target for future enhancement of blood-tumor barrier permeability.

Splicing Dysregulation as Oncogenic Driver and Passenger Factor in Brain Tumors

Brain tumors are a heterogeneous group of neoplasms ranging from almost benign to highly aggressive phenotypes. The malignancy of these tumors mostly relies on gene expression reprogramming, which is frequently accompanied by the aberrant regulation of RNA processing mechanisms. In brain tumors, defects in alternative splicing result either from the dysregulation of expression and activity of splicing factors, or from mutations in the genes encoding splicing machinery components. Aberrant splicing regulation can generate dysfunctional proteins that lead to modification of fundamental physiological cellular processes, thus contributing to the development or progression of brain tumors. Herein, we summarize the current knowledge on splicing abnormalities in brain tumors and how these alterations contribute to the disease by sustaining proliferative signaling, escaping growth suppressors, or establishing a tumor microenvironment that fosters angiogenesis and intercellular communications. Lastly, we review recent efforts aimed at developing novel splicing-targeted cancer therapies, which employ oligonucleotide-based approaches or chemical modulators of alternative splicing that elicit an impact on brain tumor biology.

The aberrant splicing of BAF45d links splicing regulation and transcription in glioblastoma

Background

Glioblastoma, the most aggressive primary brain tumor, is genetically heterogeneous. Alternative splicing (AS) plays a key role in numerous pathologies, including cancer. The objectives of our study were to determine whether aberrant AS could play a role in the malignant phenotype of glioma and to understand the mechanism underlying its aberrant regulation.

Methods

We obtained surgical samples from patients with glioblastoma who underwent 5-aminolevulinic fluorescence-guided surgery. Biopsies were taken from the tumor center as well as from adjacent normal-appearing tissue. We used a global splicing array to identify candidate genes aberrantly spliced in these glioblastoma samples. Mechanistic and functional studies were performed to elucidate the role of our top candidate splice variant, BAF45d, in glioblastoma.

Results

BAF45d is part of the switch/sucrose nonfermentable complex and plays a key role in the development of the CNS. The BAF45d/6A isoform is present in 85% of over 200 glioma samples that have been analyzed and contributes to the malignant glioma phenotype through the maintenance of an undifferentiated cellular state. We demonstrate that BAF45d splicing is mediated by polypyrimidine tract-binding protein 1 (PTBP1) and that BAF45d regulates PTBP1, uncovering a reciprocal interplay between RNA splicing regulation and transcription.

Conclusions

Our data indicate that AS is a mechanism that contributes to the malignant phenotype of glioblastoma. Understanding the consequences of this biological process will uncover new therapeutic targets for this devastating disease.

Specific inhibition of splicing factor activity by decoy RNA oligonucleotides

Alternative splicing, a fundamental step in gene expression, is deregulated in many diseases. Splicing factors (SFs), which regulate this process, are up- or down regulated or mutated in several diseases including cancer. To date, there are no inhibitors that directly inhibit the activity of SFs. We designed decoy oligonucleotides, composed of several repeats of a RNA motif, which is recognized by a single SF. Here we show that decoy oligonucleotides targeting splicing factors RBFOX1/2, SRSF1 and PTBP1, can specifically bind to their respective SFs and inhibit their splicing and biological activities both in vitro and in vivo. These decoy oligonucleotides present an approach to specifically downregulate SF activity in conditions where SFs are either up-regulated or hyperactive.

Alternative splicing, critical for gene expression, is deregulated in many diseases. Here the authors develop decoy oligonucleotides to specifically downregulate splicing factors activity.

Identification of a Novel Splice Variant of Neural Cell Adhesion Molecule in Glioblastoma Through Proteogenomics Analysis

Splice variants are known to be important in the pathophysiology of tumors, including the brain cancers. We applied a proteogenomics pipeline to identify splice variants in glioblastoma (GBM, grade IV glioma), a highly malignant brain tumor, using in-house generated mass spectrometric proteomic data and public domain RNASeq dataset. Our analysis led to the identification of a novel exon that maps to the long isoform of Neural cell adhesion molecule 1 (NCAM1), expressed on the surface of glial cells and neurons, important for cell adhesion and cell signaling. The presence of the novel exon is supported with the identification of five peptides spanning it. Additional peptides were also detected in sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) gel separated proteins from GBM patient tissue, underscoring the presence of the novel peptides in the intact brain protein. The novel exon was detected in the RNASeq dataset in 18 of 25 GBM samples and separately validated in additional 10 GBM tumor tissues using quantitative real-time-polymerase chain reaction (qRT-PCR). Both transcriptomic and proteomic data indicate downregulation of NCAM1, including the novel variant, in GBM. Domain analysis of the novel NCAM1 sequence indicates that the insertion of the novel exon contributes extra low-complexity region in the protein that may be important for protein-protein interactions and hence for cell signaling associated with tumor development. Taken together, the novel NCAM1 variant reported in this study exemplifies the importance of future multiomics research and systems biology applications in GBM.

Elucidating the microRNA-203 specific biological processes in glioblastoma cells from comprehensive RNA-sequencing transcriptome profiling

Glioblastoma (GBM) is the most common primary malignant intracranial adult brain tumor. Allelic deletion on chromosome 14q plays an essential role in GBM pathogenesis, and this chromosome 14q site was thought to harbor multiple tumor suppressor genes associated with GBM, a region that also encodes microRNA-203 (miR-203). This study was conducted to identify whole transcriptome profile changes associated with miR-203 expression by high-throughput RNA sequencing. Enrichment analyses for gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that miR-203 expression had a strong, negative effect on a number of fundamental and interconnected biological processes involved in cell growth and proliferation. The biological processes mostly influenced were p53 signaling pathway, FoxO signaling pathway, DNA replication, cell cycle, MAPK signaling pathway, and apoptosis. In total, 847 upregulated and 345 downregulated differentially expressed genes were identified in control versus miR-203 expressing glioma cells. After GO enrichment, the downregulated differentially expressed genes such as BCL2, SPARC were found to be mainly enriched in cell cycle regulation and apoptosis processes, whereas the upregulated differentially expressed genes such as CCND1, E2F1 were involved in the DNA replication and cell cycle regulation. We also performed miR-203 target analysis and found BCL2, AKT, SPARC, ROBO1, c-JUN, PDGFA, and CREB were predicted target of miR-203 and miR-203 expression suppressed the protein and mRNA levels of these target genes by western blotting and qRT-PCR analysis. Moreover, co-transfection experiments using a luciferase-based reporter assay demonstrated that miR-203 directly regulated BCL-2 expression and BCL-2 overexpression suppressed miR-203 mediated glioma cell apoptosis. These results indicate that overexpression of miR-203 coordinately regulates several oncogenic pathways in GBM.

Probing Glioblastoma Tissue Heterogeneity with Laser Capture Microdissection

Among various methods now available to isolate distinct cell populations or even single cells for DNA/RNA and proteomic analysis, laser capture microdissection (LCM) offers a unique opportunity to study cells in their topological contexts. This chapter focuses on the preparation of LCM membrane slides, tissue staining and laser microdissection of cells of interest from frozen or formalin-fixed, paraffin-embedded glioblastoma tissue.

Revealing the Determinants of Widespread Alternative Splicing Perturbation in Cancer

It is increasingly appreciated that alternative splicing plays a key role in generating functional specificity and diversity in cancer. However, the mechanisms by which cancer mutations perturb splicing remain unknown. Here, we developed a network-based strategy, DrAS-Net, to investigate more than 2.5 million variants across cancer types and link somatic mutations with cancer-specific splicing events. We identified more than 40,000 driver variant candidates and their 80,000 putative splicing targets deregulated in 33 cancer types and inferred their functional impact. Strikingly, tumors with splicing perturbations show reduced expression of immune system-related genes and increased expression of cell proliferation markers. Tumors harboring different mutations in the same gene often exhibit distinct splicing perturbations. Further stratification of 10,000 patients based on their mutation-splicing relationships identifies subtypes with distinct clinical features, including survival rates. Our work reveals how single-nucleotide changes can alter the repertoires of splicing isoforms, providing insights into oncogenic mechanisms for precision medicine.

Transcriptional signature of lymphoblastoid cell lines of BRCA1, BRCA2 and non-BRCA1/2 high risk breast cancer families

Approximately 25% of hereditary breast cancer cases are associated with a strong familial history which can be explained by mutations in BRCA1 or BRCA2 and other lower penetrance genes. The remaining high-risk families could be classified as BRCAX (non-BRCA1/2) families.

Gene expression involving alternative splicing represents a well-known mechanism regulating the expression of multiple transcripts, which could be involved in cancer development. Thus using RNA-seq methodology, the analysis of transcriptome was undertaken to potentially reveal transcripts implicated in breast cancer susceptibility and development.

RNA was extracted from immortalized lymphoblastoid cell lines of 117 women (affected and unaffected) coming from BRCA1, BRCA2 and BRCAX families. Anova analysis revealed a total of 95 transcripts corresponding to 85 different genes differentially expressed (Bonferroni corrected p-value <0.01) between those groups. Hierarchical clustering allowed distinctive subgrouping of BRCA1/2 subgroups from BRCAX individuals. We found 67 transcripts, which could discriminate BRCAX from BRCA1/BRCA2 individuals while 28 transcripts discriminate affected from unaffected BRCAX individuals.

To our knowledge, this represents the first study identifying transcripts differentially expressed in lymphoblastoid cell lines from major classes of mutation-related breast cancer subgroups, namely BRCA1, BRCA2 and BRCAX. Moreover, some transcripts could discriminate affected from unaffected BRCAX individuals, which could represent potential therapeutic targets for breast cancer treatment.

RNA processing as an alternative route to attack glioblastoma

Genomic analyses have become an important tool to identify new avenues for therapy. This is especially true for cancer types with extremely poor outcomes, since our lack of effective therapies offers no tangible clinical starting point to build upon. The highly malignant brain tumor glioblastoma (GBM) exemplifies such a refractory cancer, with only 15 month average patient survival. Analyses of several hundred GBM samples compiled by the TCGA (The Cancer Genome Atlas) have produced an extensive transcriptomic map, identified prevalent chromosomal alterations, and defined important driver mutations. Unfortunately, clinical trials based on these results have not yet delivered an improvement on outcome. It is, therefore, necessary to characterize other regulatory routes known for playing a role in tumor relapse and response to treatment. Alternative splicing affects more than 90% of the human coding genes and it is an important source for transcript variation and gene regulation. Mutations and alterations in splicing factors are highly prevalent in multiple cancers, demonstrating the potential for splicing to act as a tumor driver. As a result, numerous genes are expressed as cancer-specific splicing isoforms that are functionally distinct from the canonical isoforms found in normal tissue. These include genes that regulate cancer-critical pathways such as apoptosis, DNA repair, cell proliferation, and migration. Splicing defects can even induce genomic instability, a common characteristic of cancer, and a driver of tumor evolution. Importantly, components of the splicing machinery are targetable; multiple drugs can inhibit splicing factors or promote changes in splicing which could be exploited to begin improving clinical outcomes. Here, we review the current literature and present a case for exploring RNA processing as therapeutic route for the treatment of GBM.

Splice variants as novel targets in pancreatic ductal adenocarcinoma

Despite a wealth of genomic information, a comprehensive alternative splicing (AS) analysis of pancreatic ductal adenocarcinoma (PDAC) has not been performed yet. In the present study, we assessed whole exome-based transcriptome and AS profiles of 43 pancreas tissues using Affymetrix exon array. The AS analysis of PDAC indicated on average two AS probe-sets (ranging from 1-28) in 1,354 significantly identified protein-coding genes, with skipped exon and alternative first exon being the most frequently utilised. In addition to overrepresented extracellular matrix (ECM)-receptor interaction and focal adhesion that were also seen in transcriptome differential expression (DE) analysis, Fc gamma receptor-mediated phagocytosis and axon guidance AS genes were also highly represented. Of note, the highest numbers of AS probe-sets were found in collagen genes, which encode the characteristically abundant stroma seen in PDAC. We also describe a set of 37 'hypersensitive' genes which were frequently targeted by somatic mutations, copy number alterations, DE and AS, indicating their propensity for multidimensional regulation. We provide the most comprehensive overview of the AS landscape in PDAC with underlying changes in the spliceosomal machinery. We also collate a set of AS and DE genes encoding cell surface proteins, which present promising diagnostic and therapeutic targets in PDAC.

Tumor protein D52 expression is post-transcriptionally regulated by T-cell intercellular antigen (TIA) 1 and TIA-related protein via mRNA stability

Although tumor protein D52 (TPD52) family proteins were first identified nearly 20 years ago, their molecular regulatory mechanisms remain unclear. Therefore, we investigated the post-transcriptional regulation of TPD52 family genes. An RNA immunoprecipitation (RIP) assay showed the potential binding ability of TPD52 family mRNAs to several RNA-binding proteins, and an RNA degradation assay revealed that TPD52 is subject to more prominent post-transcriptional regulation than are TPD53 and TPD54. We subsequently focused on the 3'-untranslated region (3'-UTR) of TPD52 as a cis-acting element in post-transcriptional gene regulation. Several deletion mutants of the 3'-UTR of TPD52 mRNA were constructed and ligated to the 3'-end of a reporter green fluorescence protein gene. An RNA degradation assay revealed that a minimal cis-acting region, located in the 78-280 region of the 5'-proximal region of the 3'-UTR, stabilized the reporter mRNA. Biotin pull-down and RIP assays revealed specific binding of the region to T-cell intracellular antigen 1 (TIA-1) and TIA-1-related protein (TIAR). Knockdown of TIA-1/TIAR decreased not only the expression, but also the stability of TPD52 mRNA; it also decreased the expression and stability of the reporter gene ligated to the 3'-end of the 78-280 fragment. Stimulation of transforming growth factor-β and epidermal growth factor decreased the binding ability of these factors, resulting in decreased mRNA stability. These results indicate that the 78-280 fragment and TIA-1/TIAR concordantly contribute to mRNA stability as a cis-acting element and trans-acting factor(s), respectively. Thus, we here report the specific interactions between these elements in the post-transcriptional regulation of the TPD52 gene.

Evaluation of autoantibody signatures in meningioma patients using human proteome arrays

Meningiomas are one of the most common tumors of the Central nervous system (CNS). This study aims to identify the autoantibody biomarkers in meningiomas using high-density human proteome arrays (~17,000 full-length recombinant human proteins). Screening of sera from 15 unaffected healthy individuals, 10 individuals with meningioma grade I and 5 with meningioma grade II was performed. This comprehensive proteomics based investigation revealed the dysregulation of 489 and 104 proteins in grades I and II of meningioma, respectively, along with the enrichment of several signalling pathways, which might play a crucial role in the manifestation of the disease. Autoantibody targets like IGHG4, CRYM, EFCAB2, STAT6, HDAC7A and CCNB1 were significantly dysregulated across both the grades. Further, we compared this to the tissue proteome and gene expression profile from GEO database. Previously reported upregulated proteins from meningioma tissue-based proteomics obtained from high-resolution mass spectrometry demonstrated an aggravated autoimmune response, emphasizing the clinical relevance of these targets. Some of these targets like SELENBP1 were tested for their presence in tumor tissue using immunoblotting. In the light of highly invasive diagnostic modalities employed to diagnose CNS tumors like meningioma, these autoantibody markers offer a minimally invasive diagnostic platform which could be pursued further for clinical translation.

MYCN controls an alternative RNA splicing program in high-risk metastatic neuroblastoma

The molecular mechanisms underlying the aggressive behavior of MYCN driven neuroblastoma (NBL) is under intense investigation; however, little is known about the impact of this family of transcription factors on the splicing program. Here we used high-throughput RNA sequencing to systematically study the expression of RNA isoforms in stage 4 MYCN-amplified NBL, an aggressive subtype of metastatic NBL. We show that MYCN-amplified NBL tumors display a distinct gene splicing pattern affecting multiple cancer hallmark functions. Six splicing factors displayed unique differential expression patterns in MYCN-amplified tumors and cell lines, and the binding motifs for some of these splicing factors are significantly enriched in differentially-spliced genes. Direct binding of MYCN to promoter regions of the splicing factors PTBP1 and HNRNPA1 detected by ChIP-seq demonstrates that MYCN controls the splicing pattern by direct regulation of the expression of these key splicing factors. Furthermore, high expression of PTBP1 and HNRNPA1 was significantly associated with poor overall survival of stage4 NBL patients (p ≤ 0.05). Knocking down PTBP1, HNRNPA1 and their downstream target PKM2, an isoform of pro-tumor-growth, result in repressed growth of NBL cells. Therefore, our study reveals a novel role of MYCN in controlling global splicing program through regulation of splicing factors in addition to its well-known role in the transcription program. These findings suggest a therapeutically potential to target the key splicing factors or gene isoforms in high-risk NBL with MYCN-amplification.

Identifying the genes regulated by IDH1 via gene-chip in glioma cell U87

Glioma is the most common form of primary brain tumor. Increasing evidence show that IDH1 gene mutation is implicated in glioma. However, the mechanism involved in the progression of glioma remains unclear until now. In the study reported here, we used gene chip to identifying the genes regulated with IDH mutanted at R132. The results showed that IDH1-mutant leads to 1255 up-regulated genes and 1862 down-regulated genes in U87 cell lines. Meanwhile, GO and gene-network was performed and shown IDH1-mutant mainly affect small molecule metabolic process, mitotic cell cycle and apoptosis. This result will lay a foundation for further study of IDH1 gene function in the future.

Membrane proteome analysis of glioblastoma cell invasion

Glioblastoma multiforme (GBM) tumor invasion is facilitated by cell migration and degradation of the extracellular matrix. Invadopodia are actin-rich structures that protrude from the plasma membrane in direct contact with the extracellular matrix and are proposed to participate in epithelial-mesenchymal transition. We characterized the invasiveness of 9 established GBM cell lines using an invadopodia assay and performed quantitative mass spectrometry-based proteomic analyses on enriched membrane fractions. All GBM cells produced invadopodia, with a 65% difference between the most invasive cell line (U87MG) and the least invasive cell line (LN229) (p = 0.0001). Overall, 1,141 proteins were identified in the GBM membrane proteome; the levels of 49 proteins correlated with cell invasiveness. Ingenuity Pathway Analysis predicted activation "cell movement" (z-score = 2.608, p = 3.94E(-04)) in more invasive cells and generated a network of invasion-associated proteins with direct links to key regulators of invadopodia formation. Gene expression data relating to the invasion-associated proteins ITGA5 (integrin α5), CD97, and ANXA1 (annexin A1) showed prognostic significance in independent GBM cohorts. Fluorescence microscopy demonstrated ITGA5, CD97, and ANXA1 localization in invadopodia assays, and small interfering RNA knockdown of ITGA5 reduced invadopodia formation in U87MG cells. Thus, invasion-associated proteins, including ITGA5, may prove to be useful anti-invasive targets; volociximab, a therapeutic antibody against integrin α5β1, may be useful for treatment of patients with GBM.

Identification of recurrent regulated alternative splicing events across human solid tumors

Cancer is a complex disease that involves aberrant gene expression regulation. Discriminating the modified expression patterns driving tumor biology from the many that have no or little contribution is important for understanding cancer molecular basis. Recurrent deregulation patterns observed in multiple cancer types are enriched for such driver events. Here, we studied splicing alterations in hundreds of matched tumor and normal RNA-seq samples of eight solid cancer types. We found hundreds of cassette exons for which splicing was altered in multiple cancer types and identified a set of highly frequent altered splicing events. Specific splicing regulators, including RBFOX2, MBNL1/2 and QKI, appear to account for many splicing alteration events in multiple cancer types. Together, our results provide a first global analysis of regulated splicing alterations in cancer and identify common events with a potential causative role in solid tumor development.

Lineage-specific splicing of a brain-enriched alternative exon promotes glioblastoma progression

Tissue-specific alternative splicing is critical for the emergence of tissue identity during development, yet the role of this process in malignant transformation is undefined. Tissue-specific splicing involves evolutionarily conserved, alternative exons that represent only a minority of the total alternative exons identified. Many of these conserved exons have functional features that influence signaling pathways to profound biological effect. Here, we determined that lineage-specific splicing of a brain-enriched cassette exon in the membrane-binding tumor suppressor annexin A7 (ANXA7) diminishes endosomal targeting of the EGFR oncoprotein, consequently enhancing EGFR signaling during brain tumor progression. ANXA7 exon splicing was mediated by the ribonucleoprotein PTBP1, which is normally repressed during neuronal development. PTBP1 was highly expressed in glioblastomas due to loss of a brain-enriched microRNA (miR-124) and to PTBP1 amplification. The alternative ANXA7 splicing trait was present in precursor cells, suggesting that glioblastoma cells inherit the trait from a potential tumor-initiating ancestor and that these cells exploit this trait through accumulation of mutations that enhance EGFR signaling. Our data illustrate that lineage-specific splicing of a tissue-regulated alternative exon in a constituent of an oncogenic pathway eliminates tumor suppressor functions and promotes glioblastoma progression. This paradigm may offer a general model as to how tissue-specific regulatory mechanisms can reprogram normal developmental processes into oncogenic ones.

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.

Most expression and splicing changes in myotonic dystrophy type 1 and type 2 skeletal muscle are shared with other muscular dystrophies

The prevailing pathomechanistic paradigm for myotonic dystrophy (DM) is that aberrant expression of embryonic/fetal mRNA/protein isoforms accounts for most aspects of the pleiotropic phenotype. To identify aberrant isoforms in skeletal muscle of DM1 and DM2 patients, we performed exon-array profiling and RT-PCR validation on the largest DM sample set to date, including Duchenne, Becker and tibial muscular dystrophy (NMD) patients as disease controls, and non-disease controls. Strikingly, most expression and splicing changes in DM patients were shared with NMD controls. Comparison between DM and NMD identified almost no significant differences. We conclude that DM1 and DM2 are essentially identical for dysregulation of gene expression, and DM expression changes represent a subset of broader spectrum dystrophic changes. We found no evidence for qualitative splicing differences between DM1 and DM2. While some DM-specific splicing differences exist, most of the DM splicing differences were also seen in NMD controls. SSBP3 exon 6 missplicing was observed in all diseased muscle and led to reduced protein. We conclude there is no widespread DM-specific spliceopathy in skeletal muscle and suggest that missplicing in DM (and NMD) may not be the driving mechanism for the muscle pathology, since the same pathways show expression changes unrelated to splicing.

Aberrant splicing in neurological diseases

Splicing of precursor messenger RNA (pre-mRNA) removes the intervening sequences (introns) and joins the expressed regions (exons) in the nucleus, before an intron-containing eukaryotic mRNA transcript can be exported and translated into proteins in the cytoplasm. While some sequences are always included or excluded (constitutive splicing), others can be selectively used (alternative splicing) in this process. Particularly by alternative splicing, up to tens of thousands of variant transcripts can be produced from a single gene, which contributes greatly to the proteomic diversity for such complex cellular functions as 'wiring' neurons in the nervous system. Disruption of this process leads to aberrant splicing, which accounts for the defects of up to 50% of mutations that cause certain human genetic diseases. In this review, we describe the different mechanisms of aberrant splicing that cause or have been associated with neurological diseases.

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.

Identification and characterization of alternative exon usage linked glioblastoma multiforme survival

Background

Alternative exon usage (AEU) is an important component of gene regulation. Exon expression platforms allow the detection of associations between AEU and phenotypes such as cancer. Numerous studies have identified associations between gene expression and the brain cancer glioblastoma multiforme (GBM). The few consistent gene expression biomarkers of GBM that have been reported may be due to the limited consideration of AEU and the analytical approaches used. The objectives of this study were to develop a model that accounts for the variations in expression present between the exons within a gene and to identify AEU biomarkers of GBM survival.

Methods

The expression of exons corresponding to 25,403 genes was related to the survival of 250 individuals diagnosed with GBM in a training data set. Genes exhibiting AEU in the training data set were confirmed in an independent validation data set of 78 patients. A hierarchical mixed model that allows the consideration of covariation between exons within a gene and of the effect of the epidemiological characteristics of the patients was developed to identify associations between exon expression and patient survival. This general model describes all three possible scenarios: multi-exon genes with and without AEU, and single-exon genes.

Results

AEU associated with GBM survival was identified on 2477 genes (P-value < 5.0E-04 or FDR-adjusted P-value < 0.05). G-protein coupled receptor 98 (Gpr98) and epidermal growth factor (Egf) were among the genes exhibiting AEU with 30 and 9 exons associated with GBM survival, respectively. Pathways enriched among the AEU genes included focal adhesion, ECM-receptor interaction, ABC transporters and pathways in cancer. In addition, 24 multi-exon genes without AEU and 8 single-exon genes were associated with GBM survival (FDR-adjusted P-value < 0.05).

Conclusions

The inferred patterns of AEU were consistent with in silico AS models. The hierarchical model used offered a flexible and simple way to interpret and identify associations between survival that accommodates multi-exon genes with or without AEU and single exon genes. Our results indicate that differential expression of AEU could be used as biomarker for GBM and potentially other cancers.

PTBP1-dependent regulation of USP5 alternative RNA splicing plays a role in glioblastoma tumorigenesis

Aberrant RNA splicing is thought to play a key role in tumorigenesis. The assessment of its specific contributions is limited by the complexity of information derived from genome-wide array-based approaches. We describe how performing splicing factor-specific comparisons using both tumor and cell line data sets may more readily identify physiologically relevant tumor-specific splicing events. Affymetrix exon array data derived from glioblastoma (GBM) tumor samples with defined polypyrimidine tract-binding protein 1 (PTBP1) levels were compared with data from U251 GBM cells with and without PTBP1 knockdown. This comparison yielded overlapping gene sets that comprised only a minor fraction of each data set. The identification of a novel GBM-specific splicing event involving the USP5 gene led us to further examine its role in tumorigenesis. In GBM, USP5 generates a shorter isoform 2 through recognition of a 5' splice site within exon 15. Production of the USP5 isoform 2 was strongly correlated with PTBP1 expression in GBM tumor samples and cell lines. Splicing regulation was consistent with the presence of an intronic PTBP1 binding site and could be modulated through antisense targeting of the isoform 2 splice site to force expression of isoform 1 in GBM cells. The forced expression of USP5 isoform 1 in two GBM cell lines inhibited cell growth and migration, implying an important role for USP5 splicing in gliomagenesis. These results support a role for aberrant RNA splicing in tumorigenesis and suggest that changes in relatively few genes may be sufficient to drive the process.

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.

Integrated analysis of whole genome exon array and array-comparative genomic hybridization in gastric and colorectal cancer cells

Whole genome-scale integrated analyses of exon array and array-comparative genomic hybridization are expected to enable the identification of unknown genetic features of cancer cells. Here, we evaluated this approach in 22 gastric and colorectal cancer cell lines, focusing on protein kinase genes and genes belonging to the cadherin-catenin family. Regarding alternative splicing patterns, several cancer cell lines predominantly expressed isoform 1 of protein kinase A catalytic subunit beta (PRKACB). Paired gastric cancer specimens demonstrated that isoform 1 of PRKACB was a novel cancer-related variant transcript in gastric cancers. In addition, the exon array analysis clearly identified exon 3 or exon 3-4 skipping in catenin beta 1, a short intron insertion with exon 9 skipping in CDH1, and a deletional transcript of CDH13. These abnormal transcripts were shown to have arisen from small genomic deletions. Meanwhile, an integrated analysis of 11 gastric cancer cell lines revealed that four cell lines amplified fibroblast growth factor receptor 2, with truncated forms observed in two of the cell lines. Gene amplification, and not the truncated form, was found to determine the sensitivity to a fibroblast growth factor receptor inhibitor, indicating that our cell line panel might be useful for cell-based evaluations of specific inhibitors. Using an integrated analysis, we identified several abnormal transcripts and genomic alterations in gastric and colorectal cancer cells. Our approach might enable genetic changes to be identified more efficiently, and the present results warrant further investigation using clinical samples and integrated analyses.

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.

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).

Splicing factor hnRNPH drives an oncogenic splicing switch in gliomas

This study reveals two alternative splicing events that contribute to the development of glioma. HnRNPH is shown to control production of a pro-survival splice variant of the death-domain adaptor protein IG20-MADD and the motility-enhancing isoform of the RON receptor tyrosine kinase.

In tumours, aberrant splicing generates variants that contribute to multiple aspects of tumour establishment, progression and maintenance. We show that in glioblastoma multiforme (GBM) specimens, death-domain adaptor protein Insuloma-Glucagonoma protein 20 (IG20) is consistently aberrantly spliced to generate an antagonist, anti-apoptotic isoform (MAP-kinase activating death domain protein, MADD), which effectively redirects TNF-α/TRAIL-induced death signalling to promote survival and proliferation instead of triggering apoptosis. Splicing factor hnRNPH, which is upregulated in gliomas, controls this splicing event and similarly mediates switching to a ligand-independent, constitutively active Recepteur d′Origine Nantais (RON) tyrosine kinase receptor variant that promotes migration and invasion. The increased cell death and the reduced invasiveness caused by hnRNPH ablation can be rescued by the targeted downregulation of IG20/MADD exon 16- or RON exon 11-containing variants, respectively, using isoform-specific knockdown or splicing redirection approaches. Thus, hnRNPH activity appears to be involved in the pathogenesis and progression of malignant gliomas as the centre of a splicing oncogenic switch, which might reflect reactivation of stem cell patterns and mediates multiple key aspects of aggressive tumour behaviour, including evasion from apoptosis and invasiveness.

Tumor profiling: development of prognostic and predictive factors to guide brain tumor treatment

Primary brain tumors are a heterogeneous group of malignancies with highly variable outcomes, and diagnosis is largely based on the histological appearance of the tumors. However, the diversity of primary brain tumors has made prognostic determinations based purely on clinicopathologic variables difficult. There is an increasing body of data suggesting a significant amount of molecular diversity accounts for the heterogeneity of clinical observations, such as response to treatment and time to progression. The last decade has witnessed an explosive advance in our knowledge of the molecular genetics of brain tumors, due in large part to the availability of high-throughput profiling techniques and to the completion of the human genome sequencing project. The large amount of data generated by these efforts has enabled the identification of prognostic and predictive factors and helping to identify pathways which are driving tumor growth. Identification of biomarkers will enable better patient stratification and individualization of treatment.

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.

The use of global profiling in biomarker development for gliomas

The diffuse gliomas are a heterogeneous group of malignancies with highly variable outcomes and diagnosis is largely based on the histological appearance of the tumors. Tumor classification according to cello type and grade provides some prognostic information. However, the diversity of gliomas, within tumor type and grade categories, has made prognostic determinations based purely on clinicopathologic variables difficult. There is an increasing body of data suggesting a significant amount of molecular diversity accounts for the heterogeneity of clinical observations, such as response to treatment and time to progression. The last decade has witnessed an explosive advance in our knowledge of the molecular genetics of brain tumors, due in large part to the availability of high-throughput profiling techniques, including new sequencing methodologies as well as multidimensional profiling by the Cancer Genome Atlas project. The large amount of data generated by these efforts has enabled the identification of prognostic and predictive factors and helping to identify pathways that are driving tumor growth. Identification of biomarkers, especially when coupled to clinical trials of newer targeted therapies, will enable better patient stratification and individualization of treatment.

LATENT RANK CHANGE DETECTION FOR ANALYSIS OF SPLICE-JUNCTION MICROARRAYS WITH NONLINEAR EFFECTS

Alternative splicing of gene transcripts greatly expands the functional capacity of the genome, and certain splice isoforms may indicate specific disease states such as cancer. Splice junction microarrays interrogate thousands of splice junctions, but data analysis is difficult and error prone because of the increased complexity compared to differential gene expression analysis. We present Rank Change Detection (RCD) as a method to identify differential splicing events based upon a straightforward probabilistic model comparing the over- or underrepresentation of two or more competing isoforms. RCD has advantages over commonly used methods because it is robust to false positive errors due to nonlinear trends in microarray measurements. Further, RCD does not depend on prior knowledge of splice isoforms, yet it takes advantage of the inherent structure of mutually exclusive junctions, and it is conceptually generalizable to other types of splicing arrays or RNA-Seq. RCD specifically identifies the biologically important cases when a splice junction becomes more or less prevalent compared to other mutually exclusive junctions. The example data is from different cell lines of glioblastoma tumors assayed with Agilent microarrays.

Altered MEF2 isoforms in myotonic dystrophy and other neuromuscular disorders

Because of their central role in muscle development and maintenance, MEF2 family members represent excellent candidate effectors of the muscle pathology in myotonic dystrophy (DM). We investigated the expression and alternative splicing of all four MEF2 genes in muscle from neuromuscular disorder (NMD) patients, including DM1 and DM2. We observed MEF2A and MEF2C overexpression in all NMD muscle, including 12 MEF2-interacting genes. Exon 4 and 5 usage in MEF2A and MEF2C was different between DM and normal muscle, with DM showing the embryonic isoform. Similar splicing differences were observed in other NMD muscle. For MEF2C, missplicing was more pronounced in DM than in other dystrophies. Our data confirm dysregulation of MEF2A and MEF2C expression and splicing in several NMD, including DM. Our findings demonstrate that aberrant splicing in NMD is independent from expression of mutant repeats, and suggests that some aberrant splicing, even in DM, may be compensatory rather than primary.

Drosophila ataxin 2-binding protein 1 marks an intermediate step in the molecular differentiation of female germline cysts

In the Drosophila ovary, extrinsic signaling from the niche and intrinsic translational control machinery regulate the balance between germline stem cell maintenance and the differentiation of their daughters. However, the molecules that promote the continued stepwise development of ovarian germ cells after their exit from the niche remain largely unknown. Here, we report that the early development of germline cysts depends on the Drosophila homolog of the human ataxin 2-binding protein 1 (A2BP1) gene. Drosophila A2BP1 protein expression is first observed in the cytoplasm of 4-, 8- and 16-cell cysts, bridging the expression of the early differentiation factor Bam with late markers such as Orb, Rbp9 and Bruno encoded by arrest. The expression of A2BP1 is lost in bam, sans-fille (snf) and mei-P26 mutants, but is still present in other mutants such as rbp9 and arrest. A2BP1 alleles of varying strength produce mutant phenotypes that include germline counting defects and cystic tumors. Phenotypic analysis reveals that strong A2BP1 alleles disrupt the transition from mitosis to meiosis. These mutant cells continue to express high levels of mitotic cyclins and fail to express markers of terminal differentiation. Biochemical analysis reveals that A2BP1 isoforms bind to each other and associate with Bruno, a known translational repressor protein. These data show that A2BP1 promotes the molecular differentiation of ovarian germline cysts.

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.

Inhibition of Y-box binding protein-1 slows the growth of glioblastoma multiforme and sensitizes to temozolomide independent O6-methylguanine-DNA methyltransferase

Glioblastoma multiforme (GBM) is an aggressive type of brain tumor where <3% of newly diagnosed cases in the patients will survive >5 years. In adults, GBM is the most common type of brain tumor. It is rarer in children, where it constitutes approximately 15% of all brain tumors diagnosed. These tumors are often invasive, making surgical resection difficult. Further, they can be refractory to current therapies such as temozolomide. The current dogma is that temozolomide resistance rests on the expression of O6-methylguanine-DNA methyltransferase (MGMT) because it cleaves methylated DNA adducts formed by the drug. Our laboratory recently reported that another drug resistance gene known as the Y-box binding protein-1 (YB-1) is highly expressed in primary GBM but not in normal brain tissues based on the evaluation of primary tumors. We therefore questioned whether GBM depend on YB-1 for growth and/or response to temozolomide. Herein, we report that YB-1 inhibition reduced tumor cell invasion and growth in monolayer as well as in soft agar. Moreover, blocking this protein ultimately delayed tumor onset in mice. Importantly, inhibiting YB-1 enhanced temozolomide sensitivity in a manner that was independent of MGMT in models of adult and pediatric GBM. In conclusion, inhibiting YB-1 may be a novel way to improve the treatment of GBM.

SplicerAV: a tool for mining microarray expression data for changes in RNA processing

Background

Over the past two decades more than fifty thousand unique clinical and biological samples have been assayed using the Affymetrix HG-U133 and HG-U95 GeneChip microarray platforms. This substantial repository has been used extensively to characterize changes in gene expression between biological samples, but has not been previously mined en masse for changes in mRNA processing. We explored the possibility of using HG-U133 microarray data to identify changes in alternative mRNA processing in several available archival datasets.

Results

Data from these and other gene expression microarrays can now be mined for changes in transcript isoform abundance using a program described here, SplicerAV. Using in vivo and in vitro breast cancer microarray datasets, SplicerAV was able to perform both gene and isoform specific expression profiling within the same microarray dataset. Our reanalysis of Affymetrix U133 plus 2.0 data generated by in vitro over-expression of HRAS, E2F3, beta-catenin (CTNNB1), SRC, and MYC identified several hundred oncogene-induced mRNA isoform changes, one of which recognized a previously unknown mechanism of EGFR family activation. Using clinical data, SplicerAV predicted 241 isoform changes between low and high grade breast tumors; with changes enriched among genes coding for guanyl-nucleotide exchange factors, metalloprotease inhibitors, and mRNA processing factors. Isoform changes in 15 genes were associated with aggressive cancer across the three breast cancer datasets.

Conclusions

Using SplicerAV, we identified several hundred previously uncharacterized isoform changes induced by in vitro oncogene over-expression and revealed a previously unknown mechanism of EGFR activation in human mammary epithelial cells. We analyzed Affymetrix GeneChip data from over 400 human breast tumors in three independent studies, making this the largest clinical dataset analyzed for en masse changes in alternative mRNA processing. The capacity to detect RNA isoform changes in archival microarray data using SplicerAV allowed us to carry out the first analysis of isoform specific mRNA changes directly associated with cancer survival.

Genome-wide analysis of immune activation in human T and B cells reveals distinct classes of alternatively spliced genes

Alternative splicing of pre-mRNA is a mechanism that increases the protein diversity of a single gene by differential exon inclusion/exclusion during post-transcriptional processing. While alternative splicing is established to occur during lymphocyte activation, little is known about the role it plays during the immune response. Our study is among the first reports of a systematic genome-wide analysis of activated human T and B lymphocytes using whole exon DNA microarrays integrating alternative splicing and differential gene expression. Purified human CD2⁺ T or CD19⁺ B cells were activated using protocols to model the early events in post-transplant allograft immunity and sampled as a function of time during the process of immune activation. Here we show that 3 distinct classes of alternatively spliced and/or differentially expressed genes change in an ordered manner as a function of immune activation. We mapped our results to function-based canonical pathways and demonstrated that some are populated by only one class of genes, like integrin signaling, while other pathways, such as purine metabolism and T cell receptor signaling, are populated by all three classes of genes. Our studies augment the current view of T and B cell activation in immunity that has been based exclusively upon differential gene expression by providing evidence for a large number of molecular networks populated as a function of time and activation by alternatively spliced genes, many of which are constitutively expressed.

Comparison of Affymetrix Gene Array with the Exon Array shows potential application for detection of transcript isoform variation

BACKGROUND

The emergence of isoform-sensitive microarrays has helped fuel in-depth studies of the human transcriptome. The Affymetrix GeneChip Human Exon 1.0 ST Array (Exon Array) has been previously shown to be effective in profiling gene expression at the isoform level. More recently, the Affymetrix GeneChip Human Gene 1.0 ST Array (Gene Array) has been released for measuring gene expression and interestingly contains a large subset of probes from the Exon Array. Here, we explore the potential of using Gene Array probes to assess expression variation at the sub-transcript level. Utilizing datasets of the high quality Microarray Quality Control (MAQC) RNA samples previously assayed on the Exon Array and Gene Array, we compare the expression measurements of the two platforms to determine the performance of the Gene Array in detecting isoform variations.

RESULTS

Overall, we show that the Gene Array is comparable to the Exon Array in making gene expression calls. Moreover, to examine expression of different isoforms, we modify the Gene Array probe set definition file to enable summarization of probe intensity values at the exon level and show that the expression profiles between the two platforms are also highly correlated. Next, expression calls of previously known differentially spliced genes were compared and also show concordant results. Splicing index analysis, representing estimates of exon inclusion levels, shows a lower but good correlation between platforms. As the Gene Array contains a significant subset of probes from the Exon Array, we note that, in comparison, the Gene Array overlaps with fewer but still a high proportion of splicing events annotated in the Known Alt Events UCSC track, with abundant coverage of cassette exons. We discuss the ability of the Gene Array to detect alternative splicing and isoform variation and address its limitations.

CONCLUSION

The Gene Array is an effective expression profiling tool at gene and exon expression level, the latter made possible by probe set annotation modifications. We demonstrate that the Gene Array is capable of detecting alternative splicing and isoform variation. As expected, in comparison to the Exon Array, it is limited by reduced gene content coverage and is not able to detect as wide a range of alternative splicing events. However, for the events that can be monitored by both platforms, we estimate that the selectivity and sensitivity levels are comparable. We hope our findings will shed light on the potential extension of the Gene Array to detect alternative splicing. It should be particularly suitable for researchers primarily interested in gene expression analysis, but who may be willing to look for splicing and isoform differences within their dataset. However, we do not suggest it to be an equivalent substitute to the more comprehensive Exon Array.

Singular Value Decomposition-based Alternative Splicing Detection

Altered alternative splicing has been identified as an important factor in tumorigenesis. The Affymetrix exon tiling array is designed for detecting alternative splicing events in a transcriptome-wide fashion; however, there are currently few analysis tools that are well studied for effective detection of alternative splicing events. We propose a new screening procedure based on singular value decomposition (SVD) of the residual matrix from a robust additive model fit to probe selection region (PSR) data. With this approach, we analyze the exon tiling array data from a brain cancer study conducted at the M. D. Anderson Cancer Center, and show that the proposed SVD-based approach is able to better accommodate outlying measures and capitalize on the multidimensional group-by-PSR gene expression profiles for more effective detection of group-specific alternative splicing events as well as the PSRs that are most likely associated with the alternative splicing. Lab validation confirmed some of our findings, but the list of candidates detected with our proposed method may provide a better signpost to guide further investigations.