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

Decoding the complexity of coding and non-coding RNAs across maize anther development at the isoform level

Anther is a key male reproductive organ that is essential for the plant life cycle, from the sporophyte to the gametophyte generation. To explore isoform-level transcriptional landscape of developing anthers in maize (Zea mays L.), we analyzed Iso-Seq data from anthers collected at 10 developmental stages, together with strand-specific RNA-seq, CAGE-seq, and PAS-seq data. Of the 152,026 high-confidence full-length isoforms identified, 68.8% have not been described; these include 22,365 isoforms that originate from previously unannotated loci and 82,167 novel isoforms that originate from annotated protein-coding genes. Using our newly developed strategy to detect dynamic expression patterns of isoforms, we identified 13,899 differentially variable regions (DVRs); surprisingly, 1,275 genes contain more than two DVRs, revealing highly efficient utilization of limited genic regions. We identified 7,876 long non-coding RNAs (lncRNAs) from 4,098 loci, most of which were preferentially expressed during cell differentiation and meiosis. We also detected 371 long-range interactions involving intergenic lncRNAs (lincRNAs); interestingly, 243 were lincRNA-gene ones, and the interacting genes were highly expressed in anthers, suggesting that many potential lncRNA regulators of key genes are required for anther development. This study provides valuable resources and fundamental information for studying the essential transcripts of key genes during anther development.

Alternative splicing landscape in mouse skeletal muscle and adipose tissue: Effects of intermittent fasting and exercise

Alternative splicing contributes to diversify the cellular protein landscape, but aberrant splicing is implicated in many diseases. To which extent mis-splicing contributes to insulin resistance as the causal defect of type 2 diabetes and whether this can be reversed by lifestyle interventions is largely unknown. Therefore, RNA sequencing data from skeletal muscle and adipose tissue of diabetes-susceptible NZO mice treated with or without intermittent fasting and of healthy C57BL/6J mice subjected to exercise were analyzed for alternative splicing differences using Whippet and rMATS. Diet and exercise interventions triggered comparable levels of splicing changes, although the splicing profile of skeletal muscle appeared to be more flexible than that of adipose tissue, with 72-114 differential splicing events in muscle and less than 25 in adipose tissue. Splicing changes induced by time-restricted feeding, alternate-day fasting and exercise were generally mild, with a maximal percent spliced in (PSI) difference of 67%, indicating that alternative splicing plays a rather minor role in lifestyle-induced adaptations of muscle and adipose tissue in mice. However, intron retention contributed to the regulation of gene expression, influencing genes whose expression was directly linked to phenotypic parameters (e.g. Eno2 and Pan2). Alternate-day fasting promoted skipping of exon 7 in Mlxipl (coding for ChREBP), thereby affecting the glucose sensing module of this carbohydrate-responsive transcription factor. Both intermittent fasting and exercise training led to alternative splicing of known diabetes-related GWAS genes (e.g. Abcc8, Ifnar2, Smarcad1), highlighting the potential metabolic relevance of these changes.

Alternative splicing of transposable elements in human breast cancer

Transposable elements (TEs) drive genome evolution and can affect gene expression through diverse mechanisms. In breast cancer, disrupted regulation of TE sequences may facilitate tumor-specific transcriptomic alterations. We examine 142,514 full-length isoforms derived from long-read RNA sequencing (LR-seq) of 30 breast samples to investigate the effects of TEs on the breast cancer transcriptome. Approximately half of these isoforms contain TE sequences, and these contribute to half of the novel annotated splice junctions. We quantify splicing of these LR-seq derived isoforms in 1,135 breast tumors from The Cancer Genome Atlas (TCGA) and 1,329 healthy tissue samples from the Genotype-Tissue Expression (GTEx), and find 300 TE-overlapping tumor-specific splicing events. Some splicing events are enriched in specific breast cancer subtypes - for example, a TE-driven transcription start site upstream of ERBB2 in HER2 + tumors, and several TE-mediated splicing events are associated with patient survival and poor prognosis. The full-length sequences we capture with LR-seq reveal thousands of isoforms with signatures of RNA editing, including a novel isoform belonging to RHOA; a gene previously implicated in tumor progression. We utilize our full-length isoforms to discover polymorphic TE insertions that alter splicing and validate one of these events in breast cancer cell lines. Together, our results demonstrate the widespread effects of dysregulated TEs on breast cancer transcriptomes and highlight the advantages of long-read isoform sequencing for understanding TE biology. TE-derived isoforms may alter the expression of genes important in cancer and can potentially be used as novel, disease-specific therapeutic targets or biomarkers.One sentence summary: Transposable elements generate alternative isoforms and alter post-transcriptional regulation in human breast cancer.

SF3B1: from core splicing factor to oncogenic driver

Highly recurrent somatic mutations in the gene encoding the core splicing factor SF3B1 are drivers of multiple cancer types. SF3B1 is a scaffold protein that orchestrates multivalent protein-protein interactions within the spliceosome that are essential for recognizing the branchsite (BS) and selecting the 3' splice site during the earliest stage of pre-mRNA splicing. In this review, we first describe the molecular mechanism by which multiple oncogenic SF3B1 mutations disrupt splicing. This involves perturbation of an early spliceosomal trimeric protein complex necessary for accurate BS recognition in a subset of introns, which leads to activation of upstream branchpoints and selection of cryptic 3' splice sites. We next discuss how specific transcripts affected by aberrant splicing in SF3B1-mutant cells contribute to the initiation and progression of cancer. Finally, we highlight the prognostic value and disease phenotypes of different cancer-associated SF3B1 mutations, which is critical for developing new targeted therapeutics against SF3B1-mutant cancers still lacking in the clinic.

Splicing dysregulation in glioblastoma alters the function of cell migration-related genes

Glioblastoma (GBM) has a poor prognosis with a high recurrence and low survival rate. Previous RNA-seq analyses have revealed that alternative splicing (AS) plays a role in GBM progression. Here, we present a novel AS analysis method (Semi-Q) and describe its use to identify GBM-specific AS events. We analyzed RNA-seq data from normal brain (NB), normal human astrocytes (NHAs) and GBM samples, and found that comparison between NHA and GBM was especially informative. Importantly, this analysis revealed that genes encoding cell migration-related proteins, including filamins (FLNs) and actinins (ACTNs), were among those most affected by differential AS. Functional assays revealed that dysregulated AS of FLNA, B and C transcripts produced protein isoforms that not only altered transcription of cell proliferation-related genes but also led to enhanced cell migration, resistance to cell death and/or mitochondrial respiratory function, while a dysregulated AS isoform of ACTN4 enhanced cell migration. Together, our results indicate that cell migration and actin cytoskeleton-related genes are differentially regulated by AS in GBM, supporting a role for AS in facilitating tumor growth and invasiveness.

Neuronal cell adhesion molecule (NRCAM) variant defined by microexon skipping is an essential, antigenically distinct, and targetable proteoform in high-grade glioma

To overcome the paucity of known tumor-specific surface antigens in pediatric high-grade glioma (pHGG), we contrasted splicing patterns in pHGGs and normal brain samples. Among alternative splicing events affecting extracellular protein domains, the most pervasive alteration was the skipping of ≤30 nucleotide-long microexons. Several of these skipped microexons mapped to L1-IgCAM family members, such as NRCAM . Bulk and single-nuclei short- and long-read RNA-seq revealed uniform skipping of NRCAM microexons 5 and 19 in virtually every pHGG sample. Importantly, the Δex5Δex19 (but not the full-length) NRCAM proteoform was essential for pHGG cell migration and invasion in vitro and tumor growth in vivo. We developed a monoclonal antibody selective for Δex5Δex19 NRCAM and demonstrated that "painting" of pHGG cells with this antibody enables killing by T cells armed with an FcRI-based universal immune receptor. Thus, pHGG-specific NRCAM and possibly other L1-IgCAM proteoforms are promising and highly selective targets for adoptive immunotherapies.

Statement of significance

Existing targets for chimeric antigen receptors (CAR)-armed T cells are often shared by CNS tumors and normal tissues, creating the potential for on-target/off-tumor toxicities. Here we demonstrate that in CNS tumors of glial origin, cell adhesion molecules have alternatively spliced proteoforms, which could be targeted by highly selective therapeutic antibodies.

TranscriptDB: a transcript-centric database to study eukaryotic transcript conservation and evolution

Eukaryotic genes can encode multiple distinct transcripts through the alternative splicing (AS) of genes. Interest in the AS mechanism and its evolution across different species has stimulated numerous studies, leading to several databases that provide information on AS and transcriptome data across multiple eukaryotic species. However, existing resources do not offer information on transcript conservation and evolution between genes of multiple species. Similarly to genes, identifying conserved transcripts-those from homologous genes that have retained a similar exon composition-is useful for determining transcript homology relationships, studying transcript functions and reconstructing transcript phylogenies. To address this gap, we have developed TranscriptDB, a database dedicated to studying the conservation and evolution of transcripts within gene families. TranscriptDB offers an extensive catalog of conserved transcripts and phylogenies for 317 annotated eukaryotic species, sourced from Ensembl database version 111. It serves multiple purposes, including the exploration of gene and transcript evolution. Users can access TranscriptDB through various browsing and querying tools, including a user-friendly web interface. The incorporated web servers enable users to retrieve information on transcript evolution using their own data as input. Additionally, a REST application programming interface is available for programmatic data retrieval. A data directory is also available for bulk downloads. TranscriptDB and its resources are freely accessible at https://transcriptdb.cobius.usherbrooke.ca.

Steering research on mRNA splicing in cancer towards clinical translation

Splicing factors are affected by recurrent somatic mutations and copy number variations in several types of haematologic and solid malignancies, which is often seen as prima facie evidence that splicing aberrations can drive cancer initiation and progression. However, numerous spliceosome components also 'moonlight' in DNA repair and other cellular processes, making their precise role in cancer difficult to pinpoint. Still, few would deny that dysregulated mRNA splicing is a pervasive feature of most cancers. Correctly interpreting these molecular fingerprints can reveal novel tumour vulnerabilities and untapped therapeutic opportunities. Yet multiple technological challenges, lingering misconceptions, and outstanding questions hinder clinical translation. To start with, the general landscape of splicing aberrations in cancer is not well defined, due to limitations of short-read RNA sequencing not adept at resolving complete mRNA isoforms, as well as the shallow read depth inherent in long-read RNA-sequencing, especially at single-cell level. Although individual cancer-associated isoforms are known to contribute to cancer progression, widespread splicing alterations could be an equally important and, perhaps, more readily actionable feature of human cancers. This is to say that in addition to 'repairing' mis-spliced transcripts, possible therapeutic avenues include exacerbating splicing aberration with small-molecule spliceosome inhibitors, targeting recurrent splicing aberrations with synthetic lethal approaches, and training the immune system to recognize splicing-derived neoantigens.

Long-read Ribo-STAMP simultaneously measures transcription and translation with isoform resolution

Transcription and translation are intertwined processes in which mRNA isoforms are crucial intermediaries. However, methodological limitations in analyzing translation at the mRNA isoform level have left gaps in our understanding of critical biological processes. To address these gaps, we developed an integrated computational and experimental framework called long-read Ribo-STAMP (LR-Ribo-STAMP) that capitalizes on advancements in long-read sequencing and RNA-base editing-mediated technologies to simultaneously profile translation and transcription at both the gene and mRNA isoform levels. We also developed the EditsC metric to quantify editing and leverage the single-molecule, full-length transcript information provided by long-read sequencing. Here, we report concordance between gene-level translation profiles obtained with long-read and short-read Ribo-STAMP. We show that LR-Ribo-STAMP successfully profiles translation of mRNA isoforms and links regulatory features, such as upstream open reading frames (uORFs), to translation measurements. We apply LR-Ribo-STAMP to discovering translational differences at both the gene and isoform levels in a triple-negative breast cancer cell line under normoxia and hypoxia and find that LR-Ribo-STAMP effectively delineates orthogonal transcriptional and translation shifts between conditions. We also discover regulatory elements that distinguish translational differences at the isoform level. We highlight GRK6, in which hypoxia is observed to increase expression and translation of a shorter mRNA isoform, giving rise to a truncated protein without the AGC Kinase domain. Overall, LR-Ribo-STAMP is an important advance in our repertoire of methods that measures mRNA translation with isoform sensitivity.

Functional and dynamic profiling of transcript isoforms reveals essential roles of alternative splicing in interferon response

Type I interferon (IFN-I) plays an important role in the innate immune response through inducing IFN-I-stimulated genes (ISGs). However, how alternative splicing (AS) events, especially over time, affect their function remains poorly understood. We generated an annotation (113,843 transcripts) for IFN-I-stimulated human B cells called isoISG using high-accuracy long-read sequencing data from PacBio Sequel II/IIe. Transcript isoform profiling using isoISG revealed that isoform switching occurred in the early response to IFN-I so that ISGs would gain functional domains (e.g., C4B) or higher protein production (e.g., IRF3). Conversely, isoforms lacking functional domains increased during the late phase of IFN-I response, mainly due to intron retention events. This suggests that isoform switching both triggers and terminates IFN-I responses at the translation and protein levels. Furthermore, genetic variants influencing the isoform ratio of ISGs were associated with immunological and infectious diseases. AS has essential roles in regulating innate immune response and associated diseases.

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.

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.

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.

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.