Ribosome profiling at isoform level reveals evolutionary conserved impacts of differential splicing on the proteome

Long-read sequencing for 29 immune cell subsets reveals disease-linked isoforms

Alternative splicing events are a major causal mechanism for complex traits, but they have been understudied due to the limitation of short-read sequencing. Here, we generate a full-length isoform annotation of human immune cells from an individual by long-read sequencing for 29 cell subsets. This contains a number of unannotated transcripts and isoforms such as a read-through transcript of TOMM40-APOE in the Alzheimer's disease locus. We profile characteristics of isoforms and show that repetitive elements significantly explain the diversity of unannotated isoforms, providing insight into the human genome evolution. In addition, some of the isoforms are expressed in a cell-type specific manner, whose alternative 3'-UTRs usage contributes to their specificity. Further, we identify disease-associated isoforms by isoform switch analysis and by integration of several quantitative trait loci analyses with genome-wide association study data. Our findings will promote the elucidation of the mechanism of complex diseases via alternative splicing.

Specification of neural circuit architecture shaped by context-dependent patterned LAR-RPTP microexons

LAR-RPTPs are evolutionarily conserved presynaptic cell-adhesion molecules that orchestrate multifarious synaptic adhesion pathways. Extensive alternative splicing of LAR-RPTP mRNAs may produce innumerable LAR-RPTP isoforms that act as regulatory "codes" for determining the identity and strength of specific synapse signaling. However, no direct evidence for this hypothesis exists. Here, using targeted RNA sequencing, we detected LAR-RPTP mRNAs in diverse cell types across adult male mouse brain areas. We found pronounced cell-type-specific patterns of two microexons, meA and meB, in Ptprd mRNAs. Moreover, diverse neural circuits targeting the same neuronal populations were dictated by the expression of different Ptprd variants with distinct inclusion patterns of microexons. Furthermore, conditional ablation of Ptprd meA+ variants at presynaptic loci of distinct hippocampal circuits impaired distinct modes of synaptic transmission and objection-location memory. Activity-triggered alterations of the presynaptic Ptprd meA code in subicular neurons mediates NMDA receptor-mediated postsynaptic responses in CA1 neurons and objection-location memory. Our data provide the evidence of cell-type- and/or circuit-specific expression patterns in vivo and physiological functions of LAR-RPTP microexons that are dynamically regulated.

Expression and splicing mediate distinct biological signals

BACKGROUND

Through alternative splicing, most human genes produce multiple isoforms in a cell-, tissue-, and disease-specific manner. Numerous studies show that alternative splicing is essential for development, diseases, and their treatments. Despite these important examples, the extent and biological relevance of splicing are currently unknown.

RESULTS

To solve this problem, we developed pairedGSEA and used it to profile transcriptional changes in 100 representative RNA-seq datasets. Our systematic analysis demonstrates that changes in splicing, on average, contribute to 48.1% of the biological signal in expression analyses. Gene-set enrichment analysis furthermore indicates that expression and splicing both convey shared and distinct biological signals.

CONCLUSIONS

These findings establish alternative splicing as a major regulator of the human condition and suggest that most contemporary RNA-seq studies likely miss out on critical biological insights. We anticipate our results will contribute to the transition from a gene-centric to an isoform-centric research paradigm.

Alternative Splicing Variation: Accessing and Exploiting in Crop Improvement Programs

Alternative splicing (AS) is a gene regulatory mechanism modulating gene expression in multiple ways. AS is prevalent in all eukaryotes including plants. AS generates two or more mRNAs from the precursor mRNA (pre-mRNA) to regulate transcriptome complexity and proteome diversity. Advances in next-generation sequencing, omics technology, bioinformatics tools, and computational methods provide new opportunities to quantify and visualize AS-based quantitative trait variation associated with plant growth, development, reproduction, and stress tolerance. Domestication, polyploidization, and environmental perturbation may evolve novel splicing variants associated with agronomically beneficial traits. To date, pre-mRNAs from many genes are spliced into multiple transcripts that cause phenotypic variation for complex traits, both in model plant Arabidopsis and field crops. Cataloguing and exploiting such variation may provide new paths to enhance climate resilience, resource-use efficiency, productivity, and nutritional quality of staple food crops. This review provides insights into AS variation alongside a gene expression analysis to select for novel phenotypic diversity for use in breeding programs. AS contributes to heterosis, enhances plant symbiosis (mycorrhiza and rhizobium), and provides a mechanistic link between the core clock genes and diverse environmental clues.

Isoform-specific translational control is evolutionarily conserved in primates

Alternative splicing (AS) alters messenger RNA (mRNA) coding capacity, localization, stability, and translation. Here we use comparative transcriptomics to identify cis- acting elements coupling AS to translational control (AS-TC). We sequenced total cytosolic and polyribosome-associated mRNA from human, chimpanzee, and orangutan induced pluripotent stem cells (iPSCs), revealing thousands of transcripts with splicing differences between subcellular fractions. We found both conserved and species-specific polyribosome association patterns for orthologous splicing events. Intriguingly, alternative exons with similar polyribosome profiles between species have stronger sequence conservation than exons with lineage-specific ribosome association. These data suggest that sequence variation underlies differences in the polyribosome association. Accordingly, single nucleotide substitutions in luciferase reporters designed to model exons with divergent polyribosome profiles are sufficient to regulate translational efficiency. We used position specific weight matrices to interpret exons with species-specific polyribosome association profiles, finding that polymorphic sites frequently alter recognition motifs for trans- acting RNA binding proteins. Together, our results show that AS can regulate translation by remodeling the cis- regulatory landscape of mRNA isoforms.

Re-evaluating the impact of alternative RNA splicing on proteomic diversity

Alternative splicing (AS) constitutes a mechanism by which protein-coding genes and long non-coding RNA (lncRNA) genes produce more than a single mature transcript. From plants to humans, AS is a powerful process that increases transcriptome complexity. Importantly, splice variants produced from AS can potentially encode for distinct protein isoforms which can lose or gain specific domains and, hence, differ in their functional properties. Advances in proteomics have shown that the proteome is indeed diverse due to the presence of numerous protein isoforms. For the past decades, with the help of advanced high-throughput technologies, numerous alternatively spliced transcripts have been identified. However, the low detection rate of protein isoforms in proteomic studies raised debatable questions on whether AS contributes to proteomic diversity and on how many AS events are really functional. We propose here to assess and discuss the impact of AS on proteomic complexity in the light of the technological progress, updated genome annotation, and current scientific knowledge.

Prostate cancer resistance leads to a global deregulation of translation factors and unconventional translation

Emerging evidence associates translation factors and regulators to tumorigenesis. However, our understanding of translational changes in cancer resistance is still limited. Here, we generated an enzalutamide-resistant prostate cancer (PCa) model, which recapitulated key features of clinical enzalutamide-resistant PCa. Using this model and poly(ribo)some profiling, we investigated global translation changes that occur during acquisition of PCa resistance. We found that enzalutamide-resistant cells exhibit an overall decrease in mRNA translation with a specific deregulation in the abundance of proteins involved in mitochondrial processes and in translational regulation. However, several mRNAs escape this translational downregulation and are nonetheless bound to heavy polysomes in enzalutamide-resistant cells suggesting active translation. Moreover, expressing these corresponding genes in enzalutamide-sensitive cells promotes resistance to enzalutamide treatment. We also found increased association of long non-coding RNAs (lncRNAs) with heavy polysomes in enzalutamide-resistant cells, suggesting that some lncRNAs are actively translated during enzalutamide resistance. Consistent with these findings, expressing the predicted coding sequences of known lncRNAs JPX, CRNDE and LINC00467 in enzalutamide-sensitive cells drove resistance to enzalutamide. Taken together, this suggests that aberrant translation of specific mRNAs and lncRNAs is a strong indicator of PCa enzalutamide resistance, which points towards novel therapeutic avenues that may target enzalutamide-resistant PCa.

Alternative splicing as a source of phenotypic diversity

A major goal of evolutionary genetics is to understand the genetic processes that give rise to phenotypic diversity in multicellular organisms. Alternative splicing generates multiple transcripts from a single gene, enriching the diversity of proteins and phenotypic traits. It is well established that alternative splicing contributes to key innovations over long evolutionary timescales, such as brain development in bilaterians. However, recent developments in long-read sequencing and the generation of high-quality genome assemblies for diverse organisms has facilitated comparisons of splicing profiles between closely related species, providing insights into how alternative splicing evolves over shorter timescales. Although most splicing variants are probably non-functional, alternative splicing is nonetheless emerging as a dynamic, evolutionarily labile process that can facilitate adaptation and contribute to species divergence.

Native RNA sequencing in fission yeast reveals frequent alternative splicing isoforms

The unicellular yeast Schizosaccharomyces pombe (fission yeast) retains many of the splicing features observed in humans and is thus an excellent model to study the basic mechanisms of splicing. Nearly half the genes contain introns, but the impact of alternative splicing in gene regulation and proteome diversification remains largely unexplored. Here we leverage Oxford Nanopore Technologies native RNA sequencing (dRNA), as well as ribosome profiling data, to uncover the full range of polyadenylated transcripts and translated open reading frames. We identify 332 alternative isoforms affecting the coding sequences of 262 different genes, 97 of which occur at frequencies >20%, indicating that functional alternative splicing in S. pombe is more prevalent than previously suspected. Intron retention events make ∼80% of the cases; these events may be involved in the regulation of gene expression and, in some cases, generate novel protein isoforms, as supported by ribosome profiling data in 18 of the intron retention isoforms. One example is the rpl22 gene, in which intron retention is associated with the translation of a protein of only 13 amino acids. We also find that lowly expressed transcripts tend to have longer poly(A) tails than highly expressed transcripts, highlighting an interdependence between poly(A) tail length and transcript expression level. Finally, we discover 214 novel transcripts that are not annotated, including 158 antisense transcripts, some of which also show translation evidence. The methodologies described in this work open new opportunities to study the regulation of splicing in a simple eukaryotic model.

A comprehensive long-read isoform analysis platform and sequencing resource for breast cancer

Tumors display widespread transcriptome alterations, but the full repertoire of isoform-level alternative splicing in cancer is unknown. We developed a long-read (LR) RNA sequencing and analytical platform that identifies and annotates full-length isoforms and infers tumor-specific splicing events. Application of this platform to breast cancer samples identifies thousands of previously unannotated isoforms; ~30% affect protein coding exons and are predicted to alter protein localization and function. We performed extensive cross-validation with -omics datasets to support transcription and translation of novel isoforms. We identified 3059 breast tumor–specific splicing events, including 35 that are significantly associated with patient survival. Of these, 21 are absent from GENCODE and 10 are enriched in specific breast cancer subtypes. Together, our results demonstrate the complexity, cancer subtype specificity, and clinical relevance of previously unidentified isoforms and splicing events in breast cancer that are only annotatable by LR-seq and provide a rich resource of immuno-oncology therapeutic targets.

RiboChat: a chat-style web interface for analysis and annotation of ribosome profiling data

The increasing volume of ribosome profiling (Ribo-seq) data, computational complexity of its data processing and operational handicap of related analytical procedures present a daunting set of informatics challenges. These impose a substantial barrier to researchers particularly with no or limited bioinformatics expertise in analyzing and decoding translation information from Ribo-seq data, thus driving the need for a new research paradigm for data computation and information extraction. In this knowledge base, we herein present a novel interactive web platform, RiboChat (https://db.cngb.org/ribobench/chat.html), for direct analyzing and annotating Ribo-seq data in the form of a chat conversation. It consists of a user-friendly web interface and a backend cloud-computing service. When typing a data analysis question into the chat window, the object-text detection module will be run to recognize relevant keywords from the input text. Based on the features identified in the input, individual analytics modules are then scored to find the perfect-matching candidate. The corresponding analytics module will be further executed after checking the completion status of the uploading of datasets and configured parameters. Overall, RiboChat represents an important step forward in the emerging direction of next-generation data analytics and will enable the broad research community to conveniently decipher translation information embedded within Ribo-seq data.

APPRIS: selecting functionally important isoforms

APPRIS (https://appris.bioinfo.cnio.es) is a well-established database housing annotations for protein isoforms for a range of species. APPRIS selects principal isoforms based on protein structure and function features and on cross-species conservation. Most coding genes produce a single main protein isoform and the principal isoforms chosen by the APPRIS database best represent this main cellular isoform. Human genetic data, experimental protein evidence and the distribution of clinical variants all support the relevance of APPRIS principal isoforms. APPRIS annotations and principal isoforms have now been expanded to 10 model organisms. In this paper we highlight the most recent updates to the database. APPRIS annotations have been generated for two new species, cow and chicken, the protein structural information has been augmented with reliable models from the EMBL-EBI AlphaFold database, and we have substantially expanded the confirmatory proteomics evidence available for the human genome. The most significant change in APPRIS has been the implementation of TRIFID functional isoform scores. TRIFID functional scores are assigned to all splice isoforms, and APPRIS uses the TRIFID functional scores and proteomics evidence to determine principal isoforms when core methods cannot.

Uncovering the impacts of alternative splicing on the proteome with current omics techniques

The high‐throughput sequencing of cellular RNAs has underscored a broad effect of isoform diversification through alternative splicing on the transcriptome. Moreover, the differential production of transcript isoforms from gene loci has been recognized as a critical mechanism in cell differentiation, organismal development, and disease. Yet, the extent of the impact of alternative splicing on protein production and cellular function remains a matter of debate. Multiple experimental and computational approaches have been developed in recent years to address this question. These studies have unveiled how molecular changes at different steps in the RNA processing pathway can lead to differences in protein production and have functional effects. New and emerging experimental technologies open exciting new opportunities to develop new methods to fully establish the connection between messenger RNA expression and protein production and to further investigate how RNA variation impacts the proteome and cell function.

This article is categorized under:

RNA Processing > Splicing Regulation/Alternative Splicing

Translation > Regulation

RNA Evolution and Genomics > Computational Analyses of RNA

How alternative splicing changes the properties of plant proteins

Most plant primary transcripts undergo alternative splicing (AS), and its impact on protein diversity is a subject of intensive investigation. Several studies have uncovered various mechanisms of how particular protein splice isoforms operate. However, the common principles behind the AS effects on protein function in plants have rarely been surveyed. Here, on the selected examples, we highlight diverse tissue expression patterns, subcellular localization, enzymatic activities, abilities to bind other molecules and other relevant features. We describe how the protein isoforms mutually interact to underline their intriguing roles in altering the functionality of protein complexes. Moreover, we also discuss the known cases when these interactions have been placed inside the autoregulatory loops. This review is particularly intended for plant cell and developmental biologists who would like to gain inspiration on how the splice variants encoded by their genes of interest may coordinately work.

A high-quality genome and comparison of short- versus long-read transcriptome of the palaearctic duck Aythya fuligula (tufted duck)

BACKGROUND

The tufted duck is a non-model organism that experiences high mortality in highly pathogenic avian influenza outbreaks. It belongs to the same bird family (Anatidae) as the mallard, one of the best-studied natural hosts of low-pathogenic avian influenza viruses. Studies in non-model bird species are crucial to disentangle the role of the host response in avian influenza virus infection in the natural reservoir. Such endeavour requires a high-quality genome assembly and transcriptome.

FINDINGS

This study presents the first high-quality, chromosome-level reference genome assembly of the tufted duck using the Vertebrate Genomes Project pipeline. We sequenced RNA (complementary DNA) from brain, ileum, lung, ovary, spleen, and testis using Illumina short-read and Pacific Biosciences long-read sequencing platforms, which were used for annotation. We found 34 autosomes plus Z and W sex chromosomes in the curated genome assembly, with 99.6% of the sequence assigned to chromosomes. Functional annotation revealed 14,099 protein-coding genes that generate 111,934 transcripts, which implies a mean of 7.9 isoforms per gene. We also identified 246 small RNA families.

CONCLUSIONS

This annotated genome contributes to continuing research into the host response in avian influenza virus infections in a natural reservoir. Our findings from a comparison between short-read and long-read reference transcriptomics contribute to a deeper understanding of these competing options. In this study, both technologies complemented each other. We expect this annotation to be a foundation for further comparative and evolutionary genomic studies, including many waterfowl relatives with differing susceptibilities to avian influenza viruses.

The role of alternative splicing in adaptation and evolution

Regulation of gene expression plays a central role in adaptive divergence and evolution. Although the role of gene regulation in microevolutionary processes is gaining wide acceptance, most studies have only investigated the evolution of transcript levels, ignoring the potentially significant role of transcript structures. We argue that variation in alternative splicing plays an important and widely unexplored role in adaptation (e.g., by increasing transcriptome and/or proteome diversity, or buffering potentially deleterious genetic variation). New studies increasingly highlight the potential for independent evolution in alternative splicing and transcript level, providing alternative paths for selection to act upon. We propose that alternative splicing and transcript levels can provide contrasting, nonredundant mechanisms of equal importance for adaptive diversification of gene function and regulation.

Translation complex stabilization on messenger RNA and footprint profiling to study the RNA responses and dynamics of protein biosynthesis in the cells

During protein biosynthesis, ribosomes bind to messenger (m)RNA, locate its protein-coding information, and translate the nucleotide triplets sequentially as codons into the corresponding sequence of amino acids, forming proteins. Non-coding mRNA features, such as 5' and 3' untranslated regions (UTRs), start sites or stop codons of different efficiency, stretches of slower or faster code and nascent polypeptide interactions can alter the translation rates transcript-wise. Most of the homeostatic and signal response pathways of the cells converge on individual mRNA control, as well as alter the global translation output. Among the multitude of approaches to study translational control, one of the most powerful is to infer the locations of translational complexes on mRNA based on the mRNA fragments protected by these complexes from endonucleolytic hydrolysis, or footprints. Translation complex profiling by high-throughput sequencing of the footprints allows to quantify the transcript-wise, as well as global, alterations of translation, and uncover the underlying control mechanisms by attributing footprint locations and sizes to different configurations of the translational complexes. The accuracy of all footprint profiling approaches critically depends on the fidelity of footprint generation and many methods have emerged to preserve certain or multiple configurations of the translational complexes, often in challenging biological material. In this review, a systematic summary of approaches to stabilize translational complexes on mRNA for footprinting is presented and major findings are discussed. Future directions of translation footprint profiling are outlined, focusing on the fidelity and accuracy of inference of the native in vivo translation complex distribution on mRNA.

The dedicated chaperones of eL43, Puf6 and Loc1, can also bind RPL43 mRNA and regulate the production of this ribosomal protein

The level of ribosome biogenesis is highly associated with cell growth rate. Because many ribosomal proteins have extraribosomal functions, overexpression or insufficient supply of these proteins may impair cellular growth. Therefore, the supply of ribosomal proteins is tightly controlled in response to rRNA syntheses and environmental stimuli. In our previous study, 2 RNA-binding proteins, Puf6 and Loc1, were identified as dedicated chaperones of the ribosomal protein eL43, with which they associate to maintain its protein level and proper loading. In this study, we demonstrate that Puf6 and Loc1 interact with RPL43 mRNA. Notably, Puf6 and Loc1 usually function as a dimeric complex to bind other mRNAs; however, in this instance, the individual proteins, but not the complex form, can bind RPL43 mRNA. Thus, Puf6 or Loc1 could bind RPL43 mRNA in loc1Δ or puf6Δ, respectively. The binding of Puf6 or Loc1 caused negative effects for eL43 production: decreased RNA stability and translation of RPL43A/B mRNA. The present results suggest that these dedicated chaperones control the protein levels of eL43 from the standpoint of stability and through regulating its production.

Spatiotemporal Regulation of Transcript Isoform Expression in the Hippocampus

Proper development and plasticity of hippocampal neurons require specific RNA isoforms to be expressed in the right place at the right time. Precise spatiotemporal transcript regulation requires the incorporation of essential regulatory RNA sequences into expressed isoforms. In this review, we describe several RNA processing strategies utilized by hippocampal neurons to regulate the spatiotemporal expression of genes critical to development and plasticity. The works described here demonstrate how the hippocampus is an ideal investigative model for uncovering alternate isoform-specific mechanisms that restrict the expression of transcripts in space and time.

RNA Dysregulation: An Expanding Source of Cancer Immunotherapy Targets

Cancer transcriptomes frequently exhibit RNA dysregulation. As the resulting aberrant transcripts may be translated into cancer-specific proteins, there is growing interest in exploiting RNA dysregulation as a source of tumor antigens (TAs) and thus novel immunotherapy targets. Recent advances in high-throughput technologies and rapid accumulation of multiomic cancer profiling data in public repositories have provided opportunities to systematically characterize RNA dysregulation in cancer and identify antigen targets for immunotherapy. However, given the complexity of cancer transcriptomes and proteomes, important conceptual and technological challenges exist. Here, we highlight the expanding repertoire of TAs arising from RNA dysregulation and introduce multiomic and big data strategies for identifying optimal immunotherapy targets. We discuss extant barriers for translating these targets into effective therapies as well as the implications for future research.

Silencing of SRRM4 suppresses microexon inclusion and promotes tumor growth across cancers

RNA splicing is widely dysregulated in cancer, frequently due to altered expression or activity of splicing factors (SFs). Microexons are extremely small exons (3–27 nucleotides long) that are highly evolutionarily conserved and play critical roles in promoting neuronal differentiation and development. Inclusion of microexons in mRNA transcripts is mediated by the SF Serine/Arginine Repetitive Matrix 4 (SRRM4), whose expression is largely restricted to neural tissues. However, microexons have been largely overlooked in prior analyses of splicing in cancer, as their small size necessitates specialized computational approaches for their detection. Here, we demonstrate that despite having low expression in normal nonneural tissues, SRRM4 is further silenced in tumors, resulting in the suppression of normal microexon inclusion. Remarkably, SRRM4 is the most consistently silenced SF across all tumor types analyzed, implying a general advantage of microexon down-regulation in cancer independent of its tissue of origin. We show that this silencing is favorable for tumor growth, as decreased SRRM4 expression in tumors is correlated with an increase in mitotic gene expression, and up-regulation of SRRM4 in cancer cell lines dose-dependently inhibits proliferation in vitro and in a mouse xenograft model. Further, this proliferation inhibition is accompanied by induction of neural-like expression and splicing patterns in cancer cells, suggesting that SRRM4 expression shifts the cell state away from proliferation and toward differentiation. We therefore conclude that SRRM4 acts as a proliferation brake, and tumors gain a selective advantage by cutting off this brake.

Using data from The Cancer Genome Atlas, this study shows that the splicing factor SRRM4 and its program of differentiation-promoting microexons are downregulated across tumor types with remarkable consistency, providing tumors with a proliferative advantage.

OpenProt 2021: deeper functional annotation of the coding potential of eukaryotic genomes

OpenProt (www.openprot.org) is the first proteogenomic resource supporting a polycistronic annotation model for eukaryotic genomes. It provides a deeper annotation of open reading frames (ORFs) while mining experimental data for supporting evidence using cutting-edge algorithms. This update presents the major improvements since the initial release of OpenProt. All species support recent NCBI RefSeq and Ensembl annotations, with changes in annotations being reported in OpenProt. Using the 131 ribosome profiling datasets re-analysed by OpenProt to date, non-AUG initiation starts are reported alongside a confidence score of the initiating codon. From the 177 mass spectrometry datasets re-analysed by OpenProt to date, the unicity of the detected peptides is controlled at each implementation. Furthermore, to guide the users, detectability statistics and protein relationships (isoforms) are now reported for each protein. Finally, to foster access to deeper ORF annotation independently of one's bioinformatics skills or computational resources, OpenProt now offers a data analysis platform. Users can submit their dataset for analysis and receive the results from the analysis by OpenProt. All data on OpenProt are freely available and downloadable for each species, the release-based format ensuring a continuous access to the data. Thus, OpenProt enables a more comprehensive annotation of eukaryotic genomes and fosters functional proteomic discoveries.