Molecular diagnosis of muscular diseases in outpatient clinics: A Canadian perspective.

OBJECTIVE: To evaluate the diagnostic yield of an 89-gene panel in a large cohort of patients with suspected muscle disorders and to compare the diagnostic yield of gene panel and exome sequencing approaches. METHODS: We tested 1,236 patients from outpatient clinics across Canada using a gene panel and performed exome sequencing for 46 other patients with sequential analysis of 89 genes followed by all mendelian genes. Sequencing and analysis were performed in patients with muscle weakness or symptoms suggestive of a muscle disorder and showing at least 1 supporting clinical laboratory. RESULTS: We identified a molecular diagnosis in 187 (15.1%) of the 1,236 patients tested with the 89-gene panel. Diagnoses were distributed across 40 different genes, but 6 (DMD, RYR1, CAPN3, PYGM, DYSF, and FKRP) explained about half of all cases. Cardiac anomalies, positive family history, age <60 years, and creatine kinase >1,000 IU/L were all associated with increased diagnostic yield. Exome sequencing identified a diagnosis in 10 (21.7%) of the 46 patients tested. Among these, 3 were attributed to genes not included in the 89-gene panel. Despite differences in median coverage, only 1 of the 187 diagnoses that were identified on gene panel in the 1,236 patients could have been potentially missed if exome sequencing had been performed instead. CONCLUSIONS: Our study supports the use of gene panel testing in patients with suspected muscle disorders from outpatient clinics. It also shows that exome sequencing has a low risk of missing diagnoses compared with gene panel, while potentially increasing the diagnostic yield of patients with muscle disorders.

Hypoxia leads to significant changes in alternative splicing and elevated expression of CLK splice factor kinases in PC3 prostate cancer cells.

BACKGROUND: Mounting evidence suggests that one of the ways that cells adapt to hypoxia is through alternative splicing. The aim of this study was firstly to examine the effect of hypoxia on the alternative splicing of cancer associated genes using the prostate cancer cell line PC3 as a model. Secondly, the effect of hypoxia on the expression of several regulators of splicing was examined. METHODS: PC3 cells were grown in 1% oxygen in a hypoxic chamber for 48 h, RNA extracted and sent for high throughput PCR analysis at the RNomics platform at the University of Sherbrooke, Canada. Genes whose exon inclusion rate PSI (ψ) changed significantly were identified, and their altered exon inclusion rates verified by RT-PCR in three cell lines. The expression of splice factors and splice factor kinases in response to hypoxia was examined by qPCR and western blotting. The splice factor kinase CLK1 was inhibited with the benzothiazole TG003. RESULTS: In PC3 cells the exon inclusion rate PSI (ψ) was seen to change by > 25% in 12 cancer-associated genes; MBP, APAF1, PUF60, SYNE2, CDC42BPA, FGFR10P, BTN2A2, UTRN, RAP1GDS1, PTPN13, TTC23 and CASP9 (caspase 9). The expression of the splice factors SRSF1, SRSF2, SRSF3, SAM68, HuR, hnRNPA1, and of the splice factor kinases SRPK1 and CLK1 increased significantly in hypoxia. We also observed that the splice factor kinase CLK3, but not CLK2 and CLK4, was also induced in hypoxic DU145 prostate, HT29 colon and MCF7 breast cancer cell lines. Lastly, we show that the inhibition of CLK1 in PC3 cells with the benzothiazole TG003 increased expression of the anti-apoptotic isoform caspase 9b. CONCLUSIONS: Significant changes in alternative splicing of cancer associated genes occur in prostate cancer cells in hypoxic conditions. The expression of several splice factors and splice factor kinases increases during hypoxia, in particular the Cdc-like splice factor kinases CLK1 and CLK3. We suggest that in hypoxia the elevated expression of these regulators of splicing helps cells adapt through alternative splicing of key cancer-associated genes. We suggest that the CLK splice factor kinases could be targeted in cancers in which hypoxia contributes to resistance to therapy.

Splicing arrays reveal novel RBM10 targets, including SMN2 pre-mRNA.

BACKGROUND: RBM10 is an RNA binding protein involved in message stabilization and alternative splicing regulation. The objective of the research described herein was to identify novel targets of RBM10-regulated splicing. To accomplish this, we downregulated RBM10 in human cell lines, using small interfering RNAs, then monitored alternative splicing, using a reverse transcription-PCR screening platform. RESULTS: RBM10 knockdown (KD) provoked alterations in splicing events in 10-20% of the pre-mRNAs, most of which had not been previously identified as RBM10 targets. Hierarchical clustering of the genes affected by RBM10 KD revealed good conservation of alternative exon inclusion or exclusion across cell lines. Pathway annotation showed RAS signaling to be most affected by RBM10 KD. Of particular interest was the finding that splicing of SMN pre-mRNA, encoding the survival of motor neuron (SMN) protein, was influenced by RBM10 KD. Inhibition of RBM10 resulted in preferential expression of the full-length, exon 7 retaining, SMN transcript in four cancer cell lines and one normal skin fibroblast cell line. SMN protein is expressed from two genes, SMN1 and SMN2, but the SMN1 gene is homozygously disrupted in people with spinal muscular atrophy; as a consequence, all of the SMN that is expressed in people with this disease is from the SMN2 gene. Expression analyses using primary fibroblasts from control, carrier and spinal muscle atrophy donors demonstrated that RBM10 KD resulted in preferential expression of the full-length, exon 7 retaining, SMN2 transcript. At the protein level, upregulation of the full-length SMN2 was also observed. Re-expression of RBM10, in a stable RBM10 KD cancer cell line, correlated with a reversion of the KD effect, demonstrating specificity. CONCLUSION: Our work has not only expanded the number of pre-mRNA targets for RBM10, but identified RBM10 as a novel regulator of SMN2 alternative inclusion.

Transcriptome-wide analysis of alternative RNA splicing events in Epstein-Barr virus-associated gastric carcinomas.

Multiple human diseases including cancer have been associated with a dysregulation in RNA splicing patterns. In the current study, modifications to the global RNA splicing landscape of cellular genes were investigated in the context of Epstein-Barr virus-associated gastric cancer. Global alterations to the RNA splicing landscape of cellular genes was examined in a large-scale screen from 295 primary gastric adenocarcinomas using high-throughput RNA sequencing data. RT-PCR analysis, mass spectrometry, and co-immunoprecipitation studies were also used to experimentally validate and investigate the differential alternative splicing (AS) events that were observed through RNA-seq studies. Our study identifies alterations in the AS patterns of approximately 900 genes such as tumor suppressor genes, transcription factors, splicing factors, and kinases. These findings allowed the identification of unique gene signatures for which AS is misregulated in both Epstein-Barr virus-associated gastric cancer and EBV-negative gastric cancer. Moreover, we show that the expression of Epstein-Barr nuclear antigen 1 (EBNA1) leads to modifications in the AS profile of cellular genes and that the EBNA1 protein interacts with cellular splicing factors. These findings provide insights into the molecular differences between various types of gastric cancer and suggest a role for the EBNA1 protein in the dysregulation of cellular AS.

Global Profiling of the Cellular Alternative RNA Splicing Landscape during Virus-Host Interactions.

Alternative splicing (AS) is a central mechanism of genetic regulation which modifies the sequence of RNA transcripts in higher eukaryotes. AS has been shown to increase both the variability and diversity of the cellular proteome by changing the composition of resulting proteins through differential choice of exons to be included in mature mRNAs. In the present study, alterations to the global RNA splicing landscape of cellular genes upon viral infection were investigated using mammalian reovirus as a model. Our study provides the first comprehensive portrait of global changes in the RNA splicing signatures that occur in eukaryotic cells following infection with a human virus. We identify 240 modified alternative splicing events upon infection which belong to transcripts frequently involved in the regulation of gene expression and RNA metabolism. Using mass spectrometry, we also confirm modifications to transcript-specific peptides resulting from AS in virus-infected cells. These findings provide additional insights into the complexity of virus-host interactions as these splice variants expand proteome diversity and function during viral infection.

Global profiling of alternative RNA splicing events provides insights into molecular differences between various types of hepatocellular carcinoma.

BACKGROUND: Dysregulations in alternative splicing (AS) patterns have been associated with many human diseases including cancer. In the present study, alterations to the global RNA splicing landscape of cellular genes were investigated in a large-scale screen from 377 liver tissue samples using high-throughput RNA sequencing data. RESULTS: Our study identifies modifications in the AS patterns of transcripts encoded by more than 2500 genes such as tumor suppressor genes, transcription factors, and kinases. These findings provide insights into the molecular differences between various types of hepatocellular carcinoma (HCC). Our analysis allowed the identification of 761 unique transcripts for which AS is misregulated in HBV-associated HCC, while 68 are unique to HCV-associated HCC, 54 to HBV&HCV-associated HCC, and 299 to virus-free HCC. Moreover, we demonstrate that the expression pattern of the RNA splicing factor hnRNPC in HCC tissues significantly correlates with patient survival. We also show that the expression of the HBx protein from HBV leads to modifications in the AS profiles of cellular genes. Finally, using RNA interference and a reverse transcription-PCR screening platform, we examined the implications of cellular proteins involved in the splicing of transcripts involved in apoptosis and demonstrate the potential contribution of these proteins in AS control. CONCLUSIONS: This study provides the first comprehensive portrait of global changes in the RNA splicing signatures that occur in hepatocellular carcinoma. Moreover, these data allowed us to identify unique signatures of genes for which AS is misregulated in the different types of HCC.

RBFOX1 cooperates with MBNL1 to control splicing in muscle, including events altered in myotonic dystrophy type 1.

With the goal of identifying splicing alterations in myotonic dystrophy 1 (DM1) tissues that may yield insights into targets or mechanisms, we have surveyed mis-splicing events in three systems using a RT-PCR screening and validation platform. First, a transgenic mouse model expressing CUG-repeats identified splicing alterations shared with other mouse models of DM1. Second, using cell cultures from human embryonic muscle, we noted that DM1-associated splicing alterations were significantly enriched in cytoskeleton (e.g. SORBS1, TACC2, TTN, ACTN1 and DMD) and channel (e.g. KCND3 and TRPM4) genes. Third, of the splicing alterations occurring in adult DM1 tissues, one produced a dominant negative variant of the splicing regulator RBFOX1. Notably, half of the splicing events controlled by MBNL1 were co-regulated by RBFOX1, and several events in this category were mis-spliced in DM1 tissues. Our results suggest that reduced RBFOX1 activity in DM1 tissues may amplify several of the splicing alterations caused by the deficiency in MBNL1.

Tumor microenvironment-associated modifications of alternative splicing.

Pre-mRNA alternative splicing is modified in cancer, but the origin and specificity of these changes remain unclear. Here, we probed ovarian tumors to identify cancer-associated splicing isoforms and define the mechanism by which splicing is modified in cancer cells. Using high-throughput quantitative PCR, we monitored the expression of splice variants in laser-dissected tissues from ovarian tumors. Surprisingly, changes in alternative splicing were not limited to the tumor tissues but were also found in the tumor microenvironment. Changes in the tumor-associated splicing events were found to be regulated by splicing factors that are differentially expressed in cancer tissues. Overall, ∼20% of the alternative splicing events affected by the down-regulation of the splicing factors QKI and RBFOX2 were altered in the microenvironment of ovarian tumors. Together, our results indicate that the tumor microenvironment undergoes specific changes in alternative splicing orchestrated by a limited number of splicing factors.

Redirecting splicing with bifunctional oligonucleotides.

Ectopic modulators of alternative splicing are important tools to study the function of splice variants and for correcting mis-splicing events that cause human diseases. Such modulators can be bifunctional oligonucleotides made of an antisense portion that determines target specificity, and a non-hybridizing tail that recruits proteins or RNA/protein complexes that affect splice site selection (TOSS and TOES, respectively, for targeted oligonucleotide silencer of splicing and targeted oligonucleotide enhancer of splicing). The use of TOSS and TOES has been restricted to a handful of targets. To generalize the applicability and demonstrate the robustness of TOSS, we have tested this approach on more than 50 alternative splicing events. Moreover, we have developed an algorithm that can design active TOSS with a success rate of 80%. To produce bifunctional oligonucleotides capable of stimulating splicing, we built on the observation that binding sites for TDP-43 can stimulate splicing and improve U1 snRNP binding when inserted downstream from 5' splice sites. A TOES designed to recruit TDP-43 improved exon 7 inclusion in SMN2. Overall, our study shows that bifunctional oligonucleotides can redirect splicing on a variety of genes, justifying their inclusion in the molecular arsenal that aims to alter the production of splice variants.

RBFOX2 is an important regulator of mesenchymal tissue-specific splicing in both normal and cancer tissues.

Alternative splicing provides a critical and flexible layer of regulation intervening in many biological processes to regulate the diversity of proteins and impact cell phenotype. To identify alternative splicing differences that distinguish epithelial from mesenchymal tissues, we have investigated hundreds of cassette exons using a high-throughput reverse transcription-PCR (RT-PCR) platform. Extensive changes in splicing were noted between epithelial and mesenchymal tissues in both human colon and ovarian tissues, with many changes from mostly one splice variant to predominantly the other. Remarkably, many of the splicing differences that distinguish normal mesenchymal from normal epithelial tissues matched those that differentiate normal ovarian tissues from ovarian cancer. Furthermore, because splicing profiling could classify cancer cell lines according to their epithelial/mesenchymal characteristics, we used these cancer cell lines to identify regulators for these specific splicing signatures. By knocking down 78 potential splicing factors in five cell lines, we provide an extensive view of the complex regulatory landscape associated with the epithelial and mesenchymal states, thus revealing that RBFOX2 is an important driver of mesenchymal tissue-specific splicing.

Proteins associated with the exon junction complex also control the alternative splicing of apoptotic regulators.

Several apoptotic regulators, including Bcl-x, are alternatively spliced to produce isoforms with opposite functions. We have used an RNA interference strategy to map the regulatory landscape controlling the expression of the Bcl-x splice variants in human cells. Depleting proteins known as core (Y14 and eIF4A3) or auxiliary (RNPS1, Acinus, and SAP18) components of the exon junction complex (EJC) improved the production of the proapoptotic Bcl-x(S) splice variant. This effect was not seen when we depleted EJC proteins that typically participate in mRNA export (UAP56, Aly/Ref, and TAP) or that associate with the EJC to enforce nonsense-mediated RNA decay (MNL51, Upf1, Upf2, and Upf3b). Core and auxiliary EJC components modulated Bcl-x splicing through different cis-acting elements, further suggesting that this activity is distinct from the established EJC function. In support of a direct role in splicing control, recombinant eIF4A3, Y14, and Magoh proteins associated preferentially with the endogenous Bcl-x pre-mRNA, interacted with a model Bcl-x pre-mRNA in early splicing complexes, and specifically shifted Bcl-x alternative splicing in nuclear extracts. Finally, the depletion of Y14, eIF4A3, RNPS1, SAP18, and Acinus also encouraged the production of other proapoptotic splice variants, suggesting that EJC-associated components are important regulators of apoptosis acting at the alternative splicing level.

Alternative splicing of SYK regulates mitosis and cell survival.

Most human genes produce multiple mRNA isoforms through alternative splicing. However, the biological relevance of most splice variants remains unclear. In this study, we evaluated the functional impact of alternative splicing in cancer cells. We modulated the splicing pattern of 41 cancer-associated splicing events and scored the effects on cell growth, viability and apoptosis, identifying three isoforms essential for cell survival. Specifically, changing the splicing pattern of the spleen tyrosine kinase gene (SYK) impaired cell-cycle progression and anchorage-independent growth. Notably, exposure of cancer cells to epithelial growth factor modulated the SYK splicing pattern to promote the pro-survival isoform that is associated with cancer tissues in vivo. The data suggest that splicing of selected genes is specifically modified during tumor development to allow the expression of isoforms that promote cancer cell survival.

High-throughput quantification of splicing isoforms.

Most human messenger RNAs (mRNAs) are alternatively spliced and many exhibit disease-specific splicing patterns. However, the contribution of most splicing events to the development and maintenance of human diseases remains unclear. As the contribution of alternative splicing events to diagnosis and prognosis is becoming increasingly recognized, it becomes important to develop precise methods to quantify the abundance of these isoforms in clinical samples. Here we present a pipeline for real-time PCR annotation of splicing events (RASE) that allows accurate identification of a large number of splicing isoforms in human tissues. The RASE automatically designed specific primer pairs for 81% of all alternative splicing events in the NCBI build 36 database. Experimentally, the majority of the RASE designed primers resulted in isoform-specific amplification suitable for quantification in human cell lines or in formalin-fixed, paraffin-embedded (FFPE) RNA extract. Using this pipeline it is now possible to rapidly identify splicing isoform signatures in different types of human tissues or to validate complete sets of data generated by microarray expression profiling and deep sequencing techniques.

Cancer-associated regulation of alternative splicing.

Alternative splicing of pre-mRNA increases the diversity of protein functions. Here we show that about half of all active alternative splicing events in ovarian and breast tissues are changed in tumors, and many seem to be regulated by a single factor; sequence analysis revealed binding sites for the RNA binding protein FOX2 downstream of one-third of the exons skipped in cancer. High-resolution analysis of FOX2 binding sites defined the precise positions relative to alternative exons at which the protein may function as either a silencer or an enhancer. Most of the identified targets were shifted in the same direction by FOX2 depletion in cell lines as they were in breast and ovarian cancer tissues. Notably, we found expression of FOX2 itself is downregulated in ovarian cancer and its splicing is altered in breast cancer samples. These results suggest that the decreased expression of FOX2 in cancer tissues modulates splicing and controls proliferation.

Identification of alternative splicing markers for breast cancer.

Breast cancer is the most common cause of cancer death among women under age 50 years, so it is imperative to identify molecular markers to improve diagnosis and prognosis of this disease. Here, we present a new approach for the identification of breast cancer markers that does not measure gene expression but instead uses the ratio of alternatively spliced mRNAs as its indicator. Using a high-throughput reverse transcription-PCR-based system for splicing annotation, we monitored the alternative splicing profiles of 600 cancer-associated genes in a panel of 21 normal and 26 cancerous breast tissues. We validated 41 alternative splicing events that significantly differed in breast tumors relative to normal breast tissues. Most cancer-specific changes in splicing that disrupt known protein domains support an increase in cell proliferation or survival consistent with a functional role for alternative splicing in cancer. In a blind screen, a classifier based on the 12 best cancer-associated splicing events correctly identified cancer tissues with 96% accuracy. Moreover, a subset of these alternative splicing events could order tissues according to histopathologic grade, and 5 markers were validated in a further blind set of 19 grade 1 and 19 grade 3 tumor samples. These results provide a simple alternative for the classification of normal and cancerous breast tumor tissues and underscore the putative role of alternative splicing in the biology of cancer.

Multiple and specific mRNA processing targets for the major human hnRNP proteins.

Alternative splicing is a key mechanism regulating gene expression, and it is often used to produce antagonistic activities particularly in apoptotic genes. Heterogeneous nuclear ribonucleoparticle (hnRNP) proteins form a family of RNA-binding proteins that coat nascent pre-mRNAs. Many but not all major hnRNP proteins have been shown to participate in splicing control. The range and specificity of hnRNP protein action remain poorly documented, even for those affecting splice site selection. We used RNA interference and a reverse transcription-PCR screening platform to examine the implications of 14 of the major hnRNP proteins in the splicing of 56 alternative splicing events in apoptotic genes. Out of this total of 784 alternative splicing reactions tested in three human cell lines, 31 responded similarly to a knockdown in at least two different cell lines. On the other hand, the impact of other hnRNP knockdowns was cell line specific. The broadest effects were obtained with hnRNP K and C, two proteins whose role in alternative splicing had not previously been firmly established. Different hnRNP proteins affected distinct sets of targets with little overlap even between closely related hnRNP proteins. Overall, our study highlights the potential contribution of all of these major hnRNP proteins in alternative splicing control and shows that the targets for individual hnRNP proteins can vary in different cellular contexts.

Rheumatoid arthritis specific anti-Sa antibodies target citrullinated vimentin.

Antibodies directed to the Sa antigen are highly specific for rheumatoid arthritis (RA) and can be detected in approximately 40% of RA sera. The antigen, a doublet of protein bands of about 50 kDa, is present in placenta and in RA synovial tissue. Although it has been stated that the Sa antigen is citrullinated vimentin, experimental proof for this claim has never been published. In this study, we investigated the precise nature of the antigen. Peptide sequences that were obtained from highly purified Sa antigen were unique to vimentin. Recombinant vimentin, however, was not recognized by anti-Sa reference sera. In vivo, vimentin is subjected to various post-translational modifications, including citrullination. Since antibodies to citrullinated proteins are known to be highly specific for RA, we investigated whether Sa is citrullinated and found that Sa indeed is citrullinated vimentin. Anti-Sa antibodies thus belong to the family of anticitrullinated protein/peptide antibodies. The presence of the Sa antigen in RA synovial tissue, and the recent observation that vimentin is citrullinated in dying human macrophages, make citrullinated vimentin an interesting candidate autoantigen in RA and may provide new insights into the potential role of citrullinated synovial antigens and the antibodies directed to them in the pathophysiology of RA.

Small interfering RNA-mediated reduction in heterogeneous nuclear ribonucleoparticule A1/A2 proteins induces apoptosis in human cancer cells but not in normal mortal cell lines.

To prevent their recognition as DNA breaks, the ends of linear chromosomes are organized into telomeres, which are made of proteins bound to telomere-specific, double-stranded repeats and to single-stranded DNA extensions, the G-tails. The mammalian heterogeneous nuclear ribonucleoparticule A1 and A2 proteins can bind with high affinity to such G-tails. Moreover, previous work established that in certain mouse cells a severe reduction in the level of A1 is associated with shortened telomeric repeat tracts, and restoring A1 expression increases telomere length. Here, we document that the expression of A1/A2 proteins is elevated in a variety of human cancers, whereas A1/A2 expression is lower or absent in normal tissues. To determine whether the status of A1/A2 proteins could be improved from cancer markers to cancer targets, we used small interfering RNA-mediated RNA interference to elicit a reduction in A1/A2 proteins in a variety of human cell lines. We show that this treatment provoked specific and rapid cell death by apoptosis in cell lines derived from cervical, colon, breast, ovarian, and brain cancers. Cancer cell lines that lack p53 or express a defective p53 protein were equally sensitive to a small interfering RNA-mediated decrease in A1/A2 expression. The reduction in A1/A2 levels in HeLa cells was associated with a change in the distribution of the lengths of G-tails, an event not observed when apoptosis was induced with staurosporine. Remarkably, comparable decreases in the expression of A1/A2 in several mortal human fibroblastic and epithelial cell lines did not promote cell death. Thus, manipulating the level and activity of A1/A2 proteins may constitute a potent and specific approach in the treatment of human cancers of various origins.