Genome Sequencing and RNA-Motif Analysis Reveal Novel Damaging Noncoding Mutations in Human Tumors

Breaking paradigms: Long non-coding RNAs forming gene fusions with potential implications in cancer

Long non-coding RNAs (lncRNAs) are non-coding RNAs longer than 200 nucleotides with dynamic regulatory functions. They interact with a wide range of molecules such as DNA, RNA, and proteins to modulate diverse cellular functions through several mechanisms and, if deregulated, they can lead to cancer development and progression. Recently, it has been described that lncRNAs are susceptible to form gene fusions with mRNAs or other lncRNAs, breaking the paradigm of gene fusions consisting mainly of protein-coding genes. However, their biological significance in the tumor phenotype is still uncertain. Therefore, their recent identification opens a new line of research to study their biological role in tumorigenesis, and their potential as biomarkers with clinical relevance or as therapeutic targets. The present study aimed to review the lncRNA fusions identified so far and to know which of them have been associated with a potential function. We address the current challenges to deepen their study as well as the reasons why they represent a future therapeutic window in cancer.

Mutation density analyses on long noncoding RNA reveal comparable patterns to protein-coding RNA and prognostic value

Mutations and gene expression are the two most studied genomic features in cancer research. In the last decade, the combined advances in genomic technology and computational algorithms have broadened mutation research with the concept of mutation density and expanded the traditional scope of protein-coding RNA to noncoding RNAs. However, mutation density analysis had yet to be integrated with non-coding RNAs. In this study, we examined long non-coding RNA (lncRNA) mutation density patterns of 57 unique cancer types using 80 cancer cohorts. Our analysis revealed that lncRNAs exhibit mutation density patterns reminiscent to those of protein-coding mRNAs. These patterns include mutation peak and dip around transcription start sites of lncRNA. In many cohorts, these patterns justified statistically significant transcription strand bias, and the transcription strand bias was shared between lncRNAs and mRNAs. We further quantified transcription strand biases with a Log Odds Ratio metric and showed that some of these biases are associated with patient prognosis. The prognostic effect may be exerted due to strong Transcription-coupled repair mechanisms associated with the individual patient. For the first time, our study combined mutational density patterns with lncRNA mutations, and the results demonstrated remarkably comparable patterns between protein-coding mRNA and lncRNA, further illustrating lncRNA's potential roles in cancer research.

A convergent malignant phenotype in B-cell acute lymphoblastic leukemia involving the splicing factor SRRM1

A significant proportion of infant B-cell acute lymphoblastic leukemia (B-ALL) patients remains with a dismal prognosis due to yet undetermined mechanisms. We performed a comprehensive multicohort analysis of gene expression, gene fusions, and RNA splicing alterations to uncover molecular signatures potentially linked to the observed poor outcome. We identified 87 fusions with significant allele frequency across patients and shared functional impacts, suggesting common mechanisms across fusions. We further identified a gene expression signature that predicts high risk independently of the gene fusion background and includes the upregulation of the splicing factor SRRM1. Experiments in B-ALL cell lines provided further evidence for the role of SRRM1 on cell survival, proliferation, and invasion. Supplementary analysis revealed that SRRM1 potentially modulates splicing events associated with poor outcomes through protein-protein interactions with other splicing factors. Our findings reveal a potential convergent mechanism of aberrant RNA processing that sustains a malignant phenotype independently of the underlying gene fusion and that could potentially complement current clinical strategies in infant B-ALL.

Inferring Potential Cancer Driving Synonymous Variants

Synonymous single nucleotide variants (sSNVs) are often considered functionally silent, but a few cases of cancer-causing sSNVs have been reported. From available databases, we collected four categories of sSNVs: germline, somatic in normal tissues, somatic in cancerous tissues, and putative cancer drivers. We found that screening sSNVs for recurrence among patients, conservation of the affected genomic position, and synVep prediction (synVep is a machine learning-based sSNV effect predictor) recovers cancer driver variants (termed proposed drivers) and previously unknown putative cancer genes. Of the 2.9 million somatic sSNVs found in the COSMIC database, we identified 2111 proposed cancer driver sSNVs. Of these, 326 sSNVs could be further tagged for possible RNA splicing effects, RNA structural changes, and affected RBP motifs. This list of proposed cancer driver sSNVs provides computational guidance in prioritizing the experimental evaluation of synonymous mutations found in cancers. Furthermore, our list of novel potential cancer genes, galvanized by synonymous mutations, may highlight yet unexplored cancer mechanisms.

SRSF10: an atypical splicing regulator with critical roles in stress response, organ development, and viral replication

Serine/arginine splicing factor 10 (SRSF10) is a member of the family of mammalian splicing regulators known as SR proteins. Like several of its SR siblings, the SRSF10 protein is composed of an RNA binding domain (RRM) and of arginine and serine-rich auxiliary domains (RS) that guide interactions with other proteins. The phosphorylation status of SRSF10 is of paramount importance for its activity and is subjected to changes during mitosis, heat-shock, and DNA damage. SRSF10 overexpression has functional consequences in a growing list of cancers. By controlling the alternative splicing of specific transcripts, SRSF10 has also been implicated in glucose, fat, and cholesterol metabolism, in the development of the embryonic heart, and in neurological processes. SRSF10 is also important for the proper expression and processing of HIV-1 and other viral transcripts. We discuss how SRSF10 could become a potentially appealing therapeutic target to combat cancer and viral infections.

ISOTOPE: ISOform-guided prediction of epiTOPEs in cancer

Immunotherapies provide effective treatments for previously untreatable tumors and identifying tumor-specific epitopes can help elucidate the molecular determinants of therapy response. Here, we describe a pipeline, ISOTOPE (ISOform-guided prediction of epiTOPEs In Cancer), for the comprehensive identification of tumor-specific splicing-derived epitopes. Using RNA sequencing and mass spectrometry for MHC-I associated proteins, ISOTOPE identified neoepitopes from tumor-specific splicing events that are potentially presented by MHC-I complexes. Analysis of multiple samples indicates that splicing alterations may affect the production of self-epitopes and generate more candidate neoepitopes than somatic mutations. Although there was no difference in the number of splicing-derived neoepitopes between responders and non-responders to immune therapy, higher MHC-I binding affinity was associated with a positive response. Our analyses highlight the diversity of the immunogenic impacts of tumor-specific splicing alterations and the importance of studying splicing alterations to fully characterize tumors in the context of immunotherapies. ISOTOPE is available at https://github.com/comprna/ISOTOPE.

Exome sequencing in genetic disease: recent advances and considerations

Over the past decade, exome sequencing (ES) has allowed significant advancements to the field of disease research. By targeting the protein-coding regions of the genome, ES combines the depth of knowledge on protein-altering variants with high-throughput data generation and ease of analysis. New discoveries continue to be made using ES, and medical science has benefitted both theoretically and clinically from its continued use. In this review, we describe recent advances and successes of ES in disease research. Through selected examples of recent publications, we explore how ES continues to be a valuable tool to find variants that might explain disease etiology or provide insight into the biology underlying the disease. We then discuss shortcomings of ES in terms of variant discoveries made by other sequencing technologies that would be missed because of the scope and techniques of ES. We conclude with a brief outlook on the future of ES, suggesting that although newer and more thorough sequencing methods will soon supplant ES, its results will continue to be useful for disease research.

Roles and mechanisms of alternative splicing in cancer - implications for care

Removal of introns from messenger RNA precursors (pre-mRNA splicing) is an essential step for the expression of most eukaryotic genes. Alternative splicing enables the regulated generation of multiple mRNA and protein products from a single gene. Cancer cells have general as well as cancer type-specific and subtype-specific alterations in the splicing process that can have prognostic value and contribute to every hallmark of cancer progression, including cancer immune responses. These splicing alterations are often linked to the occurrence of cancer driver mutations in genes encoding either core components or regulators of the splicing machinery. Of therapeutic relevance, the transcriptomic landscape of cancer cells makes them particularly vulnerable to pharmacological inhibition of splicing. Small-molecule splicing modulators are currently in clinical trials and, in addition to splice site-switching antisense oligonucleotides, offer the promise of novel and personalized approaches to cancer treatment.

Prevalence and architecture of posttranscriptionally impaired synonymous mutations in 8,320 genomes across 22 cancer types

Somatic synonymous mutations are one of the most frequent genetic variants occurring in the coding region of cancer genomes, while their contributions to cancer development remain largely unknown. To assess whether synonymous mutations involved in post-transcriptional regulation contribute to the genetic etiology of cancers, we collected whole exome data from 8,320 patients across 22 cancer types. By employing our developed algorithm, PIVar, we identified a total of 22,948 posttranscriptionally impaired synonymous SNVs (pisSNVs) spanning 2,042 genes. In addition, 35 RNA binding proteins impacted by these identified pisSNVs were significantly enriched. Remarkably, we discovered markedly elevated ratio of somatic pisSNVs across all 22 cancer types, and a high pisSNV ratio was associated with worse patient survival in five cancer types. Intriguing, several well-established cancer genes, including PTEN, RB1 and PIK3CA, appeared to contribute to tumorigenesis at both protein function and posttranscriptional regulation levels, whereas some pisSNV-hosted genes, including UBR4, EP400 and INTS1, exerted their function during carcinogenesis mainly via posttranscriptional mechanisms. Moreover, we predicted three drugs associated with two pisSNVs, and numerous compounds associated with expression signature of pisSNV-hosted genes. Our study reveals the prevalence and clinical relevance of pisSNVs in cancers, and emphasizes the importance of considering posttranscriptional impaired synonymous mutations in cancer biology.

Alternative mRNA splicing in cancer immunotherapy

Immunotherapies are yielding effective treatments for several previously untreatable cancers. Still, the identification of suitable antigens specific to the tumour that can be targets for cancer vaccines and T cell therapies is a challenge. Alternative processing of mRNA, a phenomenon that has been shown to alter the proteomic diversity of many cancers, may offer the potential of a broadened target space. Here, we discuss the promise of analysing mRNA processing events in cancer cells, with an emphasis on mRNA splicing, for the identification of potential new targets for cancer immunotherapy. Further, we highlight the challenges that must be overcome for this new avenue to have clinical applicability.

Gain-of-Function Mutations: An Emerging Advantage for Cancer Biology

Advances in next-generation sequencing have identified thousands of genomic variants that perturb the normal functions of proteins, further contributing to diverse phenotypic consequences in cancer. Elucidating the functional pathways altered by loss-of-function (LOF) or gain-of-function (GOF) mutations will be crucial for prioritizing cancer-causing variants and their resultant therapeutic liabilities. In this review, we highlight the fundamental function of GOF mutations and discuss the potential mechanistic effects in the context of signaling networks. We also summarize advances in experimental and computational resources, which will dramatically help with studies on the functional and phenotypic consequences of mutations. Together, systematic investigations of the function of GOF mutations will provide an important missing piece for cancer biology and precision therapy.

Discovery of Allele-Specific Protein-RNA Interactions in Human Transcriptomes

Gene expression is tightly regulated at the post-transcriptional level through splicing, transport, translation, and decay. RNA-binding proteins (RBPs) play key roles in post-transcriptional gene regulation, and genetic variants that alter RBP-RNA interactions can affect gene products and functions. We developed a computational method ASPRIN (Allele-Specific Protein-RNA Interaction) that uses a joint analysis of CLIP-seq (cross-linking and immunoprecipitation followed by high-throughput sequencing) and RNA-seq data to identify genetic variants that alter RBP-RNA interactions by directly observing the allelic preference of RBP from CLIP-seq experiments as compared to RNA-seq. We used ASPRIN to systematically analyze CLIP-seq and RNA-seq data for 166 RBPs in two ENCODE (Encyclopedia of DNA Elements) cell lines. ASPRIN identified genetic variants that alter RBP-RNA interactions by modifying RBP binding motifs within RNA. Moreover, through an integrative ASPRIN analysis with population-scale RNA-seq data, we showed that ASPRIN can help reveal potential causal variants that affect alternative splicing via allele-specific protein-RNA interactions.

SNPs and Somatic Mutation on Long Non-Coding RNA: New Frontier in the Cancer Studies?

In the last decade, it has been demonstrated that long non-coding RNAs (lncRNAs) are involved in cancer development. The great majority of studies on lncRNAs report alterations, principally on their expression profiles, in several tumor types with respect to the normal tissues of origin. Conversely, since lncRNAs constitute a relatively novel class of RNAs compared to protein-coding transcripts (mRNAs), the landscape of their mutations and variations has not yet been extensively studied. However, in recent years an ever-increasing number of articles have described mutations of lncRNAs. Single-nucleotide polymorphisms (SNPs) that occur within the lncRNA transcripts can affect the structure and function of these RNA molecules, while the presence of a SNP in the promoter region of a lncRNA could alter its expression level. Also, somatic mutations that occur within lncRNAs have been shown to exert important effects in cancer and preliminary data are promising. Overall, the evidence suggests that SNPs and somatic mutation on lncRNAs may play a role in the pathogenesis of cancer, and indicates strong potential for further development of lncRNAs as biomarkers.

Alternative splicing links histone modifications to stem cell fate decision

BACKGROUND

Understanding the embryonic stem cell (ESC) fate decision between self-renewal and proper differentiation is important for developmental biology and regenerative medicine. Attention has focused on mechanisms involving histone modifications, alternative pre-messenger RNA splicing, and cell-cycle progression. However, their intricate interrelations and joint contributions to ESC fate decision remain unclear.

RESULTS

We analyze the transcriptomes and epigenomes of human ESC and five types of differentiated cells. We identify thousands of alternatively spliced exons and reveal their development and lineage-dependent characterizations. Several histone modifications show dynamic changes in alternatively spliced exons and three are strongly associated with 52.8% of alternative splicing events upon hESC differentiation. The histone modification-associated alternatively spliced genes predominantly function in G2/M phases and ATM/ATR-mediated DNA damage response pathway for cell differentiation, whereas other alternatively spliced genes are enriched in the G1 phase and pathways for self-renewal. These results imply a potential epigenetic mechanism by which some histone modifications contribute to ESC fate decision through the regulation of alternative splicing in specific pathways and cell-cycle genes. Supported by experimental validations and extended datasets from Roadmap/ENCODE projects, we exemplify this mechanism by a cell-cycle-related transcription factor, PBX1, which regulates the pluripotency regulatory network by binding to NANOG. We suggest that the isoform switch from PBX1a to PBX1b links H3K36me3 to hESC fate determination through the PSIP1/SRSF1 adaptor, which results in the exon skipping of PBX1.

CONCLUSION

We reveal the mechanism by which alternative splicing links histone modifications to stem cell fate decision.

MBNL1 alternative splicing isoforms play opposing roles in cancer

MBNL1 proteins lacking exon 7 (−ex7) are antisurvival factors with tumor suppressive role that cancer cells tend to down-regulate in favor of MBNL +ex7 isoforms.

The extent of and the oncogenic role played by alternative splicing (AS) in cancer are well documented. Nonetheless, only few studies have attempted to dissect individual gene function at an isoform level. Here, we focus on the AS of splicing factors during prostate cancer progression, as these factors are known to undergo extensive AS and have the potential to affect hundreds of downstream genes. We identified exon 7 (ex7) in the MBNL1 (Muscleblind-like 1) transcript as being the most differentially included exon in cancer, both in cell lines and in patients' samples. In contrast, MBNL1 overall expression was down-regulated, consistently with its described role as a tumor suppressor. This observation holds true in the majority of cancer types analyzed. We first identified components associated to the U2 splicing complex (SF3B1, SF3A1, and PHF5A) as required for efficient ex7 inclusion and we confirmed that this exon is fundamental for MBNL1 protein homodimerization. We next used splice-switching antisense oligonucleotides (AONs) or siRNAs to compare the effect of MBNL1 splicing isoform switching with knockdown. We report that whereas the absence of MBNL1 is tolerated in cancer cells, the expression of isoforms lacking ex7 (MBNL1 Δex7) induces DNA damage and inhibits cell viability and migration, acting as dominant negative proteins. Our data demonstrate the importance of studying gene function at the level of alternative spliced isoforms and support our conclusion that MBNL1 Δex7 proteins are antisurvival factors with a defined tumor suppressive role that cancer cells tend to down-regulate in favor of MBNL +ex7 isoforms.