The ribosome-engaged landscape of alternative splicing

Genome-scale exon perturbation screens uncover exons critical for cell fitness

CRISPR-Cas technology has transformed functional genomics, yet understanding of how individual exons differentially shape cellular phenotypes remains limited. Here, we optimized and conducted massively parallel exon deletion and splice-site mutation screens in human cell lines to identify exons that regulate cellular fitness. Fitness-promoting exons are prevalent in essential and highly expressed genes and commonly overlap with protein domains and interaction interfaces. Conversely, fitness-suppressing exons are enriched in nonessential genes, exhibiting lower inclusion levels, and overlap with intrinsically disordered regions and disease-associated mutations. In-depth mechanistic investigation of the screen-hit TAF5 alternative exon-8 revealed that its inclusion is required for assembly of the TFIID general transcription initiation complex, thereby regulating global gene expression output. Collectively, our orthogonal exon perturbation screens established a comprehensive repository of phenotypically important exons and uncovered regulatory mechanisms governing cellular fitness and gene expression.

Uncovering the dynamics and consequences of RNA isoform changes during neuronal differentiation

Static gene expression programs have been extensively characterized in stem cells and mature human cells. However, the dynamics of RNA isoform changes upon cell-state-transitions during cell differentiation, the determinants and functional consequences have largely remained unclear. Here, we established an improved model for human neurogenesis in vitro that is amenable for systems-wide analyses of gene expression. Our multi-omics analysis reveals that the pronounced alterations in cell morphology correlate strongly with widespread changes in RNA isoform expression. Our approach identifies thousands of new RNA isoforms that are expressed at distinct differentiation stages. RNA isoforms mainly arise from exon skipping and the alternative usage of transcription start and polyadenylation sites during human neurogenesis. The transcript isoform changes can remodel the identity and functions of protein isoforms. Finally, our study identifies a set of RNA binding proteins as a potential determinant of differentiation stage-specific global isoform changes. This work supports the view of regulated isoform changes that underlie state-transitions during neurogenesis.

Unique transcriptomes of sensory and non-sensory neurons: insights from Splicing Regulatory States

Alternative Splicing (AS) programs serve as instructive signals of cell type specificity, particularly within the brain, which comprises dozens of molecularly and functionally distinct cell types. Among them, retinal photoreceptors stand out due to their unique transcriptome, making them a particularly well-suited system for studying how AS shapes cell type-specific molecular functions. Here, we use the Splicing Regulatory State (SRS) as a novel framework to discuss the splicing factors governing the unique AS pattern of photoreceptors, and how this pattern may aid in the specification of their highly specialized sensory cilia. In addition, we discuss how other sensory cells with ciliated structures, for which data is much scarcer, also rely on specific SRSs to implement a proteome specialized in the detection of sensory stimuli. By reviewing the general rules of cell type- and tissue-specific AS programs, firstly in the brain and subsequently in specialized sensory neurons, we propose a novel paradigm on how SRSs are established and how they can diversify. Finally, we illustrate how SRSs shape the outcome of mutations in splicing factors to produce cell type-specific phenotypes that can lead to various human diseases.

Quantifying transcriptome diversity: a review

Following the central dogma of molecular biology, gene expression heterogeneity can aid in predicting and explaining the wide variety of protein products, functions and, ultimately, heterogeneity in phenotypes. There is currently overlapping terminology used to describe the types of diversity in gene expression profiles, and overlooking these nuances can misrepresent important biological information. Here, we describe transcriptome diversity as a measure of the heterogeneity in (1) the expression of all genes within a sample or a single gene across samples in a population (gene-level diversity) or (2) the isoform-specific expression of a given gene (isoform-level diversity). We first overview modulators and quantification of transcriptome diversity at the gene level. Then, we discuss the role alternative splicing plays in driving transcript isoform-level diversity and how it can be quantified. Additionally, we overview computational resources for calculating gene-level and isoform-level diversity for high-throughput sequencing data. Finally, we discuss future applications of transcriptome diversity. This review provides a comprehensive overview of how gene expression diversity arises, and how measuring it determines a more complete picture of heterogeneity across proteins, cells, tissues, organisms and species.

Global detection of human variants and isoforms by deep proteome sequencing

An average shotgun proteomics experiment detects approximately 10,000 human proteins from a single sample. However, individual proteins are typically identified by peptide sequences representing a small fraction of their total amino acids. Hence, an average shotgun experiment fails to distinguish different protein variants and isoforms. Deeper proteome sequencing is therefore required for the global discovery of protein isoforms. Using six different human cell lines, six proteases, deep fractionation and three tandem mass spectrometry fragmentation methods, we identify a million unique peptides from 17,717 protein groups, with a median sequence coverage of approximately 80%. Direct comparison with RNA expression data provides evidence for the translation of most nonsynonymous variants. We have also hypothesized that undetected variants likely arise from mutation-induced protein instability. We further observe comparable detection rates for exon-exon junction peptides representing constitutive and alternative splicing events. Our dataset represents a resource for proteoform discovery and provides direct evidence that most frame-preserving alternatively spliced isoforms are translated.

The function of alternative splicing in the proteome: rewiring protein interactomes to put old functions into new contexts

Alternative splicing affects more than 95% of multi-exon genes in the human genome. These changes affect the proteome in a myriad of ways. Here, we review our understanding of the breadth of these changes from their effect on protein structure to their influence on interactions. These changes encompass effects on nucleic acid binding in the nucleus to protein-carbohydrate interactions in the extracellular milieu, altering interactions involving all major classes of biological molecules. Protein isoforms have profound influences on cellular and tissue physiology, for example, by shaping neuronal connections, enhancing insulin secretion by pancreatic beta cells and allowing for alternative viral defense strategies in stem cells. More broadly, alternative splicing enables repurposing proteins from one context to another and thereby contributes to both the evolution of new traits as well as the creation of disease-specific interactomes that drive pathological phenotypes. In this Review, we highlight this universal character of alternative splicing as a central regulator of protein function with implications for almost every biological process.

The status of the human gene catalogue

Scientists have been trying to identify every gene in the human genome since the initial draft was published in 2001. In the years since, much progress has been made in identifying protein-coding genes, currently estimated to number fewer than 20,000, with an ever-expanding number of distinct protein-coding isoforms. Here we review the status of the human gene catalogue and the efforts to complete it in recent years. Beside the ongoing annotation of protein-coding genes, their isoforms and pseudogenes, the invention of high-throughput RNA sequencing and other technological breakthroughs have led to a rapid growth in the number of reported non-coding RNA genes. For most of these non-coding RNAs, the functional relevance is currently unclear; we look at recent advances that offer paths forward to identifying their functions and towards eventually completing the human gene catalogue. Finally, we examine the need for a universal annotation standard that includes all medically significant genes and maintains their relationships with different reference genomes for the use of the human gene catalogue in clinical settings.

Nanopore sequencing unveils the complexity of the cold-activated murine brown adipose tissue transcriptome

Alternative transcription increases transcriptome complexity by expression of multiple transcripts per gene. Annotation and quantification of transcripts using short-read sequencing is non-trivial. Long-read sequencing aims at overcoming these problems by sequencing full-length transcripts. Activation of brown adipose tissue (BAT) thermogenesis involves major transcriptomic remodeling and positively affects metabolism via increased energy expenditure. We benchmark Oxford Nanopore Technology (ONT) long-read sequencing protocols to Illumina short-read sequencing assessing alignment characteristics, gene and transcript detection and quantification, differential gene and transcript expression, transcriptome reannotation, and differential transcript usage (DTU). We find ONT sequencing is superior to Illumina for transcriptome reassembly, reducing the risk of false-positive events by unambiguously mapping reads to transcripts. We identified novel isoforms of genes undergoing DTU in cold-activated BAT including Cars2, Adtrp, Acsl5, Scp2, Aldoa, and Pde4d, validated by real-time PCR. The reannotated murine BAT transcriptome established here provides a framework for future investigations into the regulation of BAT.

Nuclear OsFKBP20-1b maintains SR34 stability and promotes the splicing of retained introns upon ABA exposure in rice

Alternative splicing (AS) is a critical means by which plants respond to changes in the environment, but few splicing factors contributing to AS have been reported and functionally characterized in rice (Oryza sativa L.). Here, we explored the function and molecular mechanism of the spliceosome-associated protein OsFKBP20-1b during AS. We determined the AS landscape of wild-type and osfkbp20-1b knockout plants upon abscisic acid (ABA) treatment by transcriptome deep sequencing. To capture the dynamics of translating intron-containing mRNAs, we blocked transcription with cordycepin and performed polysome profiling. We also analyzed whether OsFKBP20-1b and the splicing factors OsSR34 and OsSR45 function together in AS using protoplast transfection assays. We show that OsFKBP20-1b interacts with OsSR34 and regulates its stability, suggesting a role as a chaperone-like protein in the spliceosome. OsFKBP20-1b facilitates the splicing of mRNAs with retained introns after ABA treatment; some of these mRNAs are translatable and encode functional transcriptional regulators of stress-responsive genes. In addition, interacting proteins, OsSR34 and OsSR45, regulate the splicing of the same retained introns as OsFKBP20-1b after ABA treatment. Our findings reveal that spliceosome-associated immunophilin functions in alternative RNA splicing in rice by positively regulating the splicing of retained introns to limit ABA response.

Noncoding translation mitigation

Translation is pervasive outside of canonical coding regions, occurring in long noncoding RNAs, canonical untranslated regions and introns1-4, especially in ageing4-6, neurodegeneration5,7 and cancer8-10. Notably, the majority of tumour-specific antigens are results of noncoding translation11-13. Although the resulting polypeptides are often nonfunctional, translation of noncoding regions is nonetheless necessary for the birth of new coding sequences14,15. The mechanisms underlying the surveillance of translation in diverse noncoding regions and how escaped polypeptides evolve new functions remain unclear10,16-19. Functional polypeptides derived from annotated noncoding sequences often localize to membranes20,21. Here we integrate massively parallel analyses of more than 10,000 human genomic sequences and millions of random sequences with genome-wide CRISPR screens, accompanied by in-depth genetic and biochemical characterizations. Our results show that the intrinsic nucleotide bias in the noncoding genome and in the genetic code frequently results in polypeptides with a hydrophobic C-terminal tail, which is captured by the ribosome-associated BAG6 membrane protein triage complex for either proteasomal degradation or membrane targeting. By contrast, canonical proteins have evolved to deplete C-terminal hydrophobic residues. Our results reveal a fail-safe mechanism for the surveillance of unwanted translation from diverse noncoding regions and suggest a possible biochemical route for the preferential membrane localization of newly evolved proteins.

Regulation of pre-mRNA splicing: roles in physiology and disease, and therapeutic prospects

The removal of introns from mRNA precursors and its regulation by alternative splicing are key for eukaryotic gene expression and cellular function, as evidenced by the numerous pathologies induced or modified by splicing alterations. Major recent advances have been made in understanding the structures and functions of the splicing machinery, in the description and classification of physiological and pathological isoforms and in the development of the first therapies for genetic diseases based on modulation of splicing. Here, we review this progress and discuss important remaining challenges, including predicting splice sites from genomic sequences, understanding the variety of molecular mechanisms and logic of splicing regulation, and harnessing this knowledge for probing gene function and disease aetiology and for the design of novel therapeutic approaches.

The physiology of alternative splicing

Alternative splicing is a substantial contributor to the high complexity of transcriptomes of multicellular eukaryotes. In this Review, we discuss the accumulated evidence that most of this complexity is reflected at the protein level and fundamentally shapes the physiology and pathology of organisms. This notion is supported not only by genome-wide analyses but, mainly, by detailed studies showing that global and gene-specific modulations of alternative splicing regulate highly diverse processes such as tissue-specific and species-specific cell differentiation, thermal regulation, neuron self-avoidance, infrared sensing, the Warburg effect, maintenance of telomere length, cancer and autism spectrum disorders (ASD). We also discuss how mastering the control of alternative splicing paved the way to clinically approved therapies for hereditary diseases.

The status of the human gene catalogue

Scientists have been trying to identify all of the genes in the human genome since the initial draft of the genome was published in 2001. Over the intervening years, much progress has been made in identifying protein-coding genes, and the estimated number has shrunk to fewer than 20,000, although the number of distinct protein-coding isoforms has expanded dramatically. The invention of high-throughput RNA sequencing and other technological breakthroughs have led to an explosion in the number of reported non-coding RNA genes, although most of them do not yet have any known function. A combination of recent advances offers a path forward to identifying these functions and towards eventually completing the human gene catalogue. However, much work remains to be done before we have a universal annotation standard that includes all medically significant genes, maintains their relationships with different reference genomes, and describes clinically relevant genetic variants.

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.

COL3A1: Potential prognostic predictor for head and neck cancer based on immune-microenvironment alternative splicing

We aimed to identify a novel prognostic biomarker for head and neck squamous cell carcinoma (HNSCC) based on tumor immunology-related alternative splicing (AS). Data for 502 HNSCC and 44 normal samples were obtained from the TCGA database and used to establish an AS-related risk model through univariate, least absolute shrinkage, and selection operator Cox regression analyses. Fresh HNSCC and normal oral tissues were surgically obtained from 44 HNSCC patients. Western blotting and quantitative reverse transcription-PCR were used to assess gene expression levels. Kaplan-Meier was performed to evaluate patients' overall survival (OS) rate. The CIBERSORT algorithm, single-sample gene set enrichment analysis, and immune checkpoint analyses were performed to compare immune activities between subgroups. The risk model was established using 10 pivotal AS events first. Collagen Type III Alpha 1 Chain (COL3A1) were screened based on |log2FC| ≥ 1 and FDR < 0.05 criteria. COL3A1 expression levels in HNSCC tissues were elevated relative to normal tissues (p < 0.001). Moreover, COL3A1 was a reliable biomarker for HNSCC patients' prognostic prediction in both cohorts (p < 0.001, p = 0.0085, respectively). COL3A1 protein (p = 0.0054) and mRNA (p < 0.0001) levels were correlated with HNSCC differentiation. Furthermore, the T stage was correlated with COL3A1 expression (p = 0.043), and COL3A1 expression was an independent prognostic predictor for HNSCC patients (p = 0.006). Compared with the risk model, COL3A1 was better at evaluating immune cell infiltrations, immune activities, and immune checkpoint gene expressions of HNSCC lesions.

Specific expression of alternatively spliced genes in the turkey (Meleagris gallopavo) reproductive tract revealed their function in spermatogenesis and post-testicular sperm maturation

The tissue-specific profile of alternatively spliced genes (ASGs) and their involvement in reproduction processes characteristic of turkey testis, epididymis, and ductus deferens were investigated for the first time in birds. Deep sequencing of male turkey reproductive tissue RNA samples (n = 6) was performed using Illumina RNA-Seq with 2 independent methods, rMATs and SUPPA2, for differential alternative splicing (DAS) event prediction. The expression of selected ASGs was validated using quantitative real-time reverse transcriptase-polymerase chain reaction. The testis was found to be the site of the highest number of posttranscriptional splicing events within the reproductive tract, and skipping exons were the most frequently occurring class of alternative splicing (AS) among the reproductive tract. Statistical analysis revealed 86, 229, and 6 DAS events in the testis/epididymis, testis/ductus deferens, and epididymis/ductus deferens comparison, respectively. Alternative splicing was found to be a mechanism of gene expression regulation within the turkey reproduction tract. In testis, modification was observed for spermatogenesis specific genes; the changes in 5' UTR could act as regulator of MEIG1 expression (a player during spermatocytes meiosis), and modification of 3' UTR led to diversification of CREM mRNA (modulator of gene expression related to the structuring of mature spermatozoa). Sperm tail formation can be regulated by changes in the 5' UTR of testicular SLC9A3R1 and gene silencing by producing dysfunctional variants of ODF2 in the testis and ATP1B3 in the epididymis. Predicted differentially ASGs in the turkey reproductive tract seem to be involved in the regulation of spermatogenesis, including acrosome formation and sperm tail formation and binding of sperm to the zona pellucida. Several ASGs were classified as cilia by actin and microtubule cytoskeleton organization. Such genes may play a role in the organization of sperm flagellum and post-testicular motility development. To our knowledge, this is the first functional investigation of alternatively spliced genes associated with tissue-specific processes in the turkey reproductive tract.

A post-transcriptional regulatory landscape of aging in the female mouse hippocampus

Aging is associated with substantial physiological changes and constitutes a major risk factor for neurological disorders including dementia. Alterations in gene expression upon aging have been extensively studied; however, an in-depth characterization of post-transcriptional regulatory events remains elusive. Here, we profiled the age-related changes of the transcriptome and translatome in the female mouse hippocampus by RNA sequencing of total RNA and polysome preparations at four ages (3-, 6-, 12-, 20-month-old); and we implemented a variety of bioinformatics approaches to unravel alterations in transcript abundance, alternative splicing, and polyadenylation site selection. We observed mostly well-coordinated transcriptome and translatome expression signatures across age including upregulation of transcripts related to immune system processes and neuroinflammation, though transcripts encoding ribonucleoproteins or associated with mitochondrial functions, calcium signaling and the cell-cycle displayed substantial discordant profiles, suggesting translational control associated with age-related deficits in hippocampal-dependent behavior. By contrast, alternative splicing was less preserved, increased with age and was associated with distinct functionally-related transcripts encoding proteins acting at synapses/dendrites, RNA-binding proteins; thereby predicting regulatory roles for RBM3 and CIRBP. Only minor changes in polyadenylation site selection were identified, indicating pivotal 3'-end selection in young adults compared to older groups. Overall, our study provides a comprehensive resource of age-associated post-transcriptional regulatory events in the mouse hippocampus, enabling further examination of the molecular features underlying age-associated neurological diseases.

Integrated top-down and bottom-up proteomics mass spectrometry for the characterization of endogenous ribosomal protein heterogeneity

Ribosomes are abundant, large RNA-protein complexes that are the sites of all protein synthesis in cells. Defects in ribosomal proteins (RPs), including proteoforms arising from genetic variations, alternative splicing of RNA transcripts, post-translational modifications and alterations of protein expression level, have been linked to a diverse range of diseases, including cancer and aging. Comprehensive characterization of ribosomal proteoforms is challenging but important for the discovery of potential disease biomarkers or protein targets. In the present work, using E. coli 70S RPs as an example, we first developed a top-down proteomics approach on a Waters Synapt G2 Si mass spectrometry (MS) system, and then applied it to the HeLa 80S ribosome. The results were complemented by a bottom-up approach. In total, 50 out of 55 RPs were identified using the top-down approach. Among these, more than 30 RPs were found to have their N-terminal methionine removed. Additional modifications such as methylation, acetylation, and hydroxylation were also observed, and the modification sites were identified by bottom-up MS. In a HeLa 80S ribosomal sample, we identified 98 ribosomal proteoforms, among which multiple truncated 80S ribosomal proteoforms were observed, the type of information which is often overlooked by bottom-up experiments. Although their relevance to diseases is not yet known, the integration of top-down and bottom-up proteomics approaches paves the way for the discovery of proteoform-specific disease biomarkers or targets.

Expression dynamics of periodic transcripts during cancer cell cycle progression and their correlation with anticancer drug sensitivity

BACKGROUND

The cell cycle is at the center of cellular activities and is orchestrated by complex regulatory mechanisms, among which transcriptional regulation is one of the most important components. Alternative splicing dramatically expands the regulatory network by producing transcript isoforms of genes to exquisitely control the cell cycle. However, the patterns of transcript isoform expression in the cell cycle are unclear. Therapies targeting cell cycle checkpoints are commonly used as anticancer therapies, but none of them have been designed or evaluated at the alternative splicing transcript level. The utility of these transcripts as markers of cell cycle-related drug sensitivity is still unknown, and studies on the expression patterns of cell cycle-targeting drug-related transcripts are also rare.

METHODS

To explore alternative splicing patterns during cell cycle progression, we performed sequential transcriptomic assays following cell cycle synchronization in colon cancer HCT116 and breast cancer MDA-MB-231 cell lines, using flow cytometry and reference cell cycle transcripts to confirm the cell cycle phases of samples, and we developed a new algorithm to describe the periodic patterns of transcripts fluctuating during the cell cycle. Genomics of Drug Sensitivity in Cancer (GDSC) drug sensitivity datasets and Cancer Cell Line Encyclopedia (CCLE) transcript datasets were used to assess the correlation of genes and their transcript isoforms with drug sensitivity. We identified transcripts associated with typical drugs targeting cell cycle by determining correlation coefficients. Cytotoxicity assays were used to confirm the effect of ENST00000257904 against cyclin dependent kinase 4/6 (CDK4/6) inhibitors. Finally, alternative splicing transcripts associated with mitotic (M) phase arrest were analyzed using an RNA synthesis inhibition assay and transcriptome analysis.

RESULTS

We established high-resolution transcriptome datasets of synchronized cell cycle samples from colon cancer HCT116 and breast cancer MDA-MB-231 cells. The results of the cell cycle assessment showed that 43,326, 41,578 and 29,244 transcripts were found to be periodically expressed in HeLa, HCT116 and MDA-MB-231 cells, respectively, among which 1280 transcripts showed this expression pattern in all three cancer cell lines. Drug sensitivity assessments showed that a large number of these transcripts displayed a higher correlation with drug sensitivity than their corresponding genes. Cell cycle-related drug screening showed that the level of the CDK4 transcript ENST00000547281 was more significantly associated with the resistance of cells to CDK4/6 inhibitors than the level of the CDK4 reference transcript ENST00000257904. The transcriptional inhibition assay following M phase arrest further confirmed the M-phase-specific expression of the splicing transcripts. Combined with the cell cycle-related drug screening, the results also showed that a set of periodic transcripts, for example, ENST00000314392 (a dolichyl-phosphate mannosyltransferase polypeptide 2 isoform transcript), was more associated with drug sensitivity than the levels of their corresponding gene transcripts.

CONCLUSIONS

In summary, we identified a panel of cell cycle-related periodic transcripts and found that the levels of transcripts of drug target genes showed different values for predicting drug sensitivity, providing novel insights into alternative splicing-related drug development and evaluation.

Where protein structure and cell diversity meet

Protein-protein interaction networks - interactomes - are charted with the hope to understand how phenotypes emerge and how they are altered in disease states. Early efforts to map interactomes have focused on the assembly of context agnostic, reference networks. However, recent studies have mapped interactomes across different cell lines and tissues, finding highly variable interactomes due to the rewiring of protein-protein interactions in different contexts. Increasing evidence points to significant links between protein structure and interactome diversity seen across cell types and tissues. We discuss how recent findings support the key role of alternative splicing and phosphorylation, two well-established regulators of protein structural and functional diversity, in defining cell type- and tissue-specific interactomes. Moreover, we show that intrinsically disordered protein regions are most favorably equipped to support interactome rewiring by acting as hubs of protein structure and function regulation.

Patterns and Crucial Regulation of Alternative Splicing During Early Development in Zebrafish

Many vertebrate genes generate multiple transcript variants that may encode functionally distinct protein isoforms, but the transcriptomes of various developmental stages are poorly defined. Identifying the transcriptome and its regulation during the normal developmental process is the key to deciphering the developmental stage-specific functions of genes. Here we presented a systematic assessment of the temporal alternative splicing (AS) events during the critical development stages to capture the dynamic gene expression changes and AS in zebrafish. An unexpected transcriptome complexity generated by AS was observed during zebrafish development. The patterns of AS events varied substantially among developmental stages despite the similarities in the total proportion of AS genes. We further found that AS afforded substantial functional diversification of genes through the generation of stage-specific AS events from broadly protein-coding genes as an essential developmental regulatory mechanism. Skipped exon (SE) showed the strongest signals among developmental AS (devAS), suggesting that devAS events generated by SE may be necessary for the normal development of zebrafish. Most developmental genes regulated by AS mechanisms were not modulated in terms of their overall expression levels, indicating that AS shaped the transcriptome independently from transcriptional regulation during development. 128-cell stage was a critical stage for gene transcription during embryonic development. Splicing factors as an essential developmental regulator underwent AS in the potential autoregulatory feedback loop and expressed multiple isoforms. Thus, zebrafish development was shaped by an interplay of programs controlling gene expression levels and AS. Overall, we provided a global view of developmental patterns of AS during zebrafish development and revealed that AS transitions were the crucial regulatory component of zebrafish embryonic development.

RNA splicing regulators play critical roles in neurogenesis

Alternative RNA splicing increases transcript diversity in different cell types and under varying conditions. It is executed with the help of RNA splicing regulators (RSRs), which are operationally defined as RNA-binding proteins (RBPs) that regulate alternative splicing, but not directly catalyzing the chemical reactions of splicing. By systematically searching for RBPs and manually identifying those that regulate splicing, we curated 305 RSRs in the human genome. Surprisingly, most of the RSRs are involved in neurogenesis. Among these RSRs, we focus on nine families (PTBP, NOVA, RBFOX, ELAVL, CELF, DBHS, MSI, PCBP, and MBNL) that play essential roles in the neurogenic pathway. A better understanding of their functions will provide novel insights into the role of splicing in brain development, health, and disease. This comprehensive review serves as a stepping-stone to explore the diverse and complex set of RSRs as fundamental regulators of neural development. This article is categorized under: RNA-Based Catalysis &gt; RNA Catalysis in Splicing and Translation RNA Interactions with Proteins and Other Molecules &gt; Protein-RNA Interactions: Functional Implications RNA Processing &gt; Splicing Regulation/Alternative Splicing.

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.

Regulating the regulator: a survey of mechanisms from transcription to translation controlling expression of mammalian cell cycle kinase Aurora A

Aurora Kinase A (AURKA) is a positive regulator of mitosis with a strict cell cycle-dependent expression pattern. Recently, novel oncogenic roles of AURKA have been uncovered that are independent of the kinase activity and act within multiple signalling pathways, including cell proliferation, survival and cancer stem cell phenotypes. For this, cellular abundance of AURKA protein is per se crucial and must be tightly fine-tuned. Indeed, AURKA is found overexpressed in different cancers, typically as a result of gene amplification or enhanced transcription. It has however become clear that impaired processing, decay and translation of AURKA mRNA can also offer the basis for altered AURKA levels. Accordingly, the involvement of gene expression mechanisms controlling AURKA expression in human diseases is increasingly recognized and calls for much more research. Here, we explore and create an integrated view of the molecular processes regulating AURKA expression at the level of transcription, post-transcription and translation, intercalating discussion on how impaired regulation underlies disease. Given that targeting AURKA levels might affect more functions compared to inhibiting the kinase activity, deeper understanding of its gene expression may aid the design of alternative and therapeutically more successful ways of suppressing the AURKA oncogene.

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.

A transcriptomic atlas underlying developmental plasticity of seasonal forms of Bicyclus anynana butterflies

Organisms residing in regions with alternating seasons often develop different phenotypes, or forms, in each season. These forms are often adaptations to each season and result from an altered developmental response to specific environmental cues such as temperature. While multiple studies have examined form-specific gene expression profiles in a diversity of species, little is known about how environments and developmental transitions, cued by hormone pulses, alter post-transcriptional patterns. In this study, we examine how gene expression, alternative splicing, and miRNA-mediated gene silencing in Bicyclus anynana butterfly hindwing tissue, varies across two rearing temperatures at four developmental timepoints. These timepoints flank two temperature-sensitive periods that coincide with two pulses of the insect hormone 20E. Our results suggest that developmental transitions, coincident with 20E pulses, elicit a greater impact on all these transcriptomic patterns than rearing temperatures per se. More similar transcriptomic patterns are observed pre-20E pulses than those observed post-20E pulses. We also found functionally distinct sets of differentially expressed and differentially spliced genes in the seasonal forms. Furthermore, around 10% of differentially expressed genes are predicted to be direct targets of, and regulated by, differentially expressed miRNAs between the seasonal forms. Many differentially expressed genes, miRNAs, or differentially spliced genes potentially regulate eyespot size plasticity, and we validated the differential splicing pattern of one such gene, daughterless. We present a comprehensive and interactive transcriptomic atlas of the hindwing tissue of both seasonal forms of B. anynana throughout development, a model organism of seasonal plasticity.

Mechanisms of non-coding RNA-modulated alternative splicing in cancer

Alternative splicing (AS) is a common and pivotal process for eukaryotic gene expression regulation, which enables a precursor RNA to produce multiple transcript variants with diverse cellular functions. Aberrant AS represents a hallmark of cancer, engaged in all stages of tumorigenesis from initiation to metastasis. Accumulating pieces of evidence have revealed the involvement of non-coding RNAs (ncRNAs) in regulating AS in human cancers. In this review, we overview the underlying mechanisms of non-coding RNAs, including microRNAs (miRNAs), long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) modulated AS at diverse levels in human cancers, and summarize their regulatory functions in tumorigenesis.

Enhanced protein isoform characterization through long-read proteogenomics

BACKGROUND

The detection of physiologically relevant protein isoforms encoded by the human genome is critical to biomedicine. Mass spectrometry (MS)-based proteomics is the preeminent method for protein detection, but isoform-resolved proteomic analysis relies on accurate reference databases that match the sample; neither a subset nor a superset database is ideal. Long-read RNA sequencing (e.g., PacBio or Oxford Nanopore) provides full-length transcripts which can be used to predict full-length protein isoforms.

RESULTS

We describe here a long-read proteogenomics approach for integrating sample-matched long-read RNA-seq and MS-based proteomics data to enhance isoform characterization. We introduce a classification scheme for protein isoforms, discover novel protein isoforms, and present the first protein inference algorithm for the direct incorporation of long-read transcriptome data to enable detection of protein isoforms previously intractable to MS-based detection. We have released an open-source Nextflow pipeline that integrates long-read sequencing in a proteomic workflow for isoform-resolved analysis.

CONCLUSIONS

Our work suggests that the incorporation of long-read sequencing and proteomic data can facilitate improved characterization of human protein isoform diversity. Our first-generation pipeline provides a strong foundation for future development of long-read proteogenomics and its adoption for both basic and translational research.

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

Computational Analysis of Alternative Splicing Using VAST-TOOLS and the VastDB Framework

Alternative splicing (AS) can vastly expand animal transcriptomes and proteomes. Two main open questions in the field are how AS is regulated across cell/tissue types and disease, and what roles different AS events play. To facilitate AS research, we have created the computational VastDB framework, which comprises a series of complementary software and resources that we describe in this chapter. The VastDB framework is especially designed to aid biomedical researchers without a strong computational background. It offers tools and resources to: (a) quantify AS and identify differentially spliced AS events using RNA-seq data (vast-tools), (b) perform multiple genomic and sequence analyses for investigating AS events (Matt), (c) identify AS events with genomic and regulatory conservation among species (ExOrthist), and (d) help with the biological interpretation of the results, and, ultimately, with the identification of interesting AS events to design wet-lab experiments (VastDB and PastDB).

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.

An Optimized Protocol for the Mapping of Cell Type-Specific Ribosome-Associated Transcript Isoforms from Small Mouse Brain Regions

Over the past years, technological advances in transcriptomics provided deep insights into gene expression programs and their role in tissue organization and cellular functions. The isolation of ribosome-associated transcripts is a powerful approach for deep profiling of cell type-specific transcripts, and particularly well-suited for quantitative analysis of transcript isoforms. This method employs conditional ribosome epitope-tagging in genetically defined cell types, followed by affinity-isolation of ribosome-associated mRNAs. Advantages of this approach are twofold: first, the method enables rapid retrieval of mRNAs without tissue dissociation and cell sorting steps. Second, capturing of ribosome-associated mRNAs, enriches for transcripts recruited for active translation, therefore providing an approximation to the cellular translatome. Here, we describe one application of this method for the identification of the transcriptome of excitatory neuronal cells (mitral and tufted cells) of the mouse olfactory bulb, through RiboTag isolation from the vGlut2-IRES-cre mouse line as genetic driver of endogenously tagged ribosome expression.

Temporal Modulation of Differential Alternative Splicing in HaCaT Human Keratinocyte Cell Line Chronically Exposed to Arsenic for up to 28 Wk

BACKGROUND

Chronic arsenic exposure via drinking water is associated with an increased risk of developing cancer and noncancer chronic diseases. Pre-mRNAs are often subject to alternative splicing, generating mRNA isoforms encoding functionally distinct protein isoforms. The resulting imbalance in isoform species can result in pathogenic changes in critical signaling pathways. Alternative splicing as a mechanism of arsenic-induced toxicity and carcinogenicity is understudied.

OBJECTIVE

This study aimed to accurately profile differential alternative splicing events in human keratinocytes induced by chronic arsenic exposure that might play a role in carcinogenesis.

METHODS

Independent quadruplicate cultures of immortalized human keratinocytes (HaCaT) were maintained continuously for 28 wk with 0 or 100 nM sodium arsenite. RNA-sequencing (RNA-Seq) was performed with poly(A) RNA isolated from cells harvested at 7, 19, and 28 wk with subsequent replicate multivariate analysis of transcript splicing (rMATS) analysis to detect and quantify differential alternative splicing events. Reverse transcriptase-polymerase chain reaction (RT-PCR) for selected alternative splicing events was performed to validate RNA-Seq predictions. Functional enrichment was performed by gene ontology (GO) analysis of the differential alternative splicing event data set at each time point.

RESULTS

At least 600 differential alternative splicing events were detected at each time point tested, comprising all the five main types of alternative splicing and occurring in both open reading frames (ORFs) and untranslated regions (UTRs). Based on functional relevance ELK4, SHC1, and XRRA1 were selected for validation of predicted alternative splicing events at 7 wk by RT-PCR. Densitometric analysis of RT-PCR data corroborated the rMATS predicted alternative splicing for all three events. Protein expression validation of the selected alternative splicing events was challenging given that very few isoform-specific antibodies are available. GO analysis demonstrated that the enriched terms in differential alternatively spliced mRNAs changed dynamically with the time of exposure. Notably, RNA metabolism and splicing regulation pathways were enriched at the 7-wk time point, when the greatest number of differentially alternatively spliced mRNAs are detected. Our preliminary proteomic analysis demonstrated that the expression of the canonical isoforms of the splice regulators DDX42, RMB25, and SRRM2 were induced upon chronic arsenic exposure, corroborating the splicing predictions.

DISCUSSION

These results using cultures of HaCaT cells suggest that arsenic exposure disrupted an alternative splice factor network and induced time-dependent genome-wide differential alternative splicing that likely contributed to the changing proteomic landscape in arsenic-induced carcinogenesis. However, significant challenges remain in corroborating alternative splicing data at the proteomic level. https://doi.org/10.1289/EHP9676.

Transcriptome-Wide Analyses Identify Dominant as the Predominantly Non-Conservative Alternative Splicing Inheritance Patterns in F1 Chickens

Transcriptome analysis has been used to investigate many economically traits in chickens; however, alternative splicing still lacks a systematic method of study that is able to promote proteome diversity, and fine-tune expression dynamics. Hybridization has been widely utilized in chicken breeding due to the resulting heterosis, but the dynamic changes in alternative splicing during this process are significant yet unclear. In this study, we performed a reciprocal crossing experiment involving the White Leghorn and Cornish Game chicken breeds which exhibit major differences in body size and reproductive traits, and conducted RNA sequencing of the brain, muscle, and liver tissues to identify the inheritance patterns. A total of 40 515 and 42 612 events were respectively detected in the brain and muscle tissues, with 39 843 observed in the liver; 2807, 4242, and 4538 events significantly different between two breeds were identified in the brain, muscle, and liver tissues, respectively. The hierarchical cluster of tissues from different tissues from all crosses, based on the alternative splicing profiles, suggests high tissue and strain specificity. Furthermore, a comparison between parental strains and hybrid crosses indicated that over one third of alternative splicing genes showed conserved patterns in all three tissues, while the second prevalent pattern was non-additive, which included both dominant and transgressive patterns; this meant that the dominant pattern plays a more important role than suppression. Our study provides an overview of the inheritance patterns of alternative splicing in layer and broiler chickens, to better understand post-transcriptional regulation during hybridization.

Differential fates of introns in gene expression due to global alternative splicing

The discovery of introns over four decades ago revealed a new vision of genes and their interrupted arrangement. Throughout the years, it has appeared that introns play essential roles in the regulation of gene expression. Unique processing of excised introns through the formation of lariats suggests a widespread role for these molecules in the structure and function of cells. In addition to rapid destruction, these lariats may linger on in the nucleus or may even be exported to the cytoplasm, where they remain stable circular RNAs (circRNAs). Alternative splicing (AS) is a source of diversity in mature transcripts harboring retained introns (RI-mRNAs). Such RNAs may contain one or more entire retained intron(s) (RIs), but they may also have intron fragments resulting from sequential excision of smaller subfragments via recursive splicing (RS), which is characteristic of long introns. There are many potential fates of RI-mRNAs, including their downregulation via nuclear and cytoplasmic surveillance systems and the generation of new protein isoforms with potentially different functions. Various reports have linked the presence of such unprocessed transcripts in mammals to important roles in normal development and in disease-related conditions. In certain human neurological-neuromuscular disorders, including myotonic dystrophy type 2 (DM2), frontotemporal dementia/amyotrophic lateral sclerosis (FTD/ALS) and Duchenne muscular dystrophy (DMD), peculiar processing of long introns has been identified and is associated with their pathogenic effects. In this review, we discuss different mechanisms involved in the processing of introns during AS and the functions of these large sections of the genome in our biology.

Comparative ribosome profiling reveals distinct translational landscapes of salt-sensitive and -tolerant rice

BACKGROUND

Soil salinization represents a serious threat to global rice production. Although significant research has been conducted to understand salt stress at the genomic, transcriptomic and proteomic levels, few studies have focused on the translatomic responses to this stress. Recent studies have suggested that transcriptional and translational responses to salt stress can often operate independently.

RESULTS

We sequenced RNA and ribosome-protected fragments (RPFs) from the salt-sensitive rice (O. sativa L.) cultivar 'Nipponbare' (NB) and the salt-tolerant cultivar 'Sea Rice 86' (SR86) under normal and salt stress conditions. A large discordance between salt-induced transcriptomic and translatomic alterations was found in both cultivars, with more translationally regulated genes being observed in SR86 in comparison to NB. A biased ribosome occupancy, wherein RPF depth gradually increased from the 5' ends to the 3' ends of coding regions, was revealed in NB and SR86. This pattern was strengthened by salt stress, particularly in SR86. On the contrary, the strength of ribosome stalling was accelerated in salt-stressed NB but decreased in SR86.

CONCLUSIONS

This study revealed that translational reprogramming represents an important layer of salt stress responses in rice, and the salt-tolerant cultivar SR86 adopts a more flexible translationally adaptive strategy to cope with salt stress compared to the salt susceptible cultivar NB. The differences in translational dynamics between NB and SR86 may derive from their differing levels of ribosome stalling under salt stress.

The three as: Alternative splicing, alternative polyadenylation and their impact on apoptosis in immune function

The latest advances in next-generation sequencing studies and transcriptomic profiling over the past decade have highlighted a surprising frequency of genes regulated by RNA processing mechanisms in the immune system. In particular, two control steps in mRNA maturation, namely alternative splicing and alternative polyadenylation, are now recognized to occur in the vast majority of human genes. Both have the potential to alter the identity of the encoded protein, as well as control protein abundance or even protein localization or association with other factors. In this review, we will provide a summary of the general mechanisms by which alternative splicing (AS) and alternative polyadenylation (APA) occur, their regulation within cells of the immune system, and their impact on immunobiology. In particular, we will focus on how control of apoptosis by AS and APA is used to tune cell fate during an immune response.

Isoform Age - Splice Isoform Profiling Using Long-Read Technologies

Alternative splicing (AS) of RNA is a key mechanism that results in the expression of multiple transcript isoforms from single genes and leads to an increase in the complexity of both the transcriptome and proteome. Regulation of AS is critical for the correct functioning of many biological pathways, while disruption of AS can be directly pathogenic in diseases such as cancer or cause risk for complex disorders. Current short-read sequencing technologies achieve high read depth but are limited in their ability to resolve complex isoforms. In this review we examine how long-read sequencing (LRS) technologies can address this challenge by covering the entire RNA sequence in a single read and thereby distinguish isoform changes that could impact RNA regulation or protein function. Coupling LRS with technologies such as single cell sequencing, targeted sequencing and spatial transcriptomics is producing a rapidly expanding suite of technological approaches to profile alternative splicing at the isoform level with unprecedented detail. In addition, integrating LRS with genotype now allows the impact of genetic variation on isoform expression to be determined. Recent results demonstrate the potential of these techniques to elucidate the landscape of splicing, including in tissues such as the brain where AS is particularly prevalent. Finally, we also discuss how AS can impact protein function, potentially leading to novel therapeutic targets for a range of diseases.

Assessing conservation of alternative splicing with evolutionary splicing graphs

Understanding how protein function has evolved and diversified is of great importance for human genetics and medicine. Here, we tackle the problem of describing the whole transcript variability observed in several species by generalizing the definition of splicing graph. We provide a practical solution to construct parsimonious evolutionary splicing graphs where each node is a minimal transcript building block defined across species. We show a clear link between the functional relevance, tissue regulation, and conservation of alternative transcripts on a set of 50 genes. By scaling up to the whole human protein-coding genome, we identify a few thousand genes where alternative splicing modulates the number and composition of pseudorepeats. We have implemented our approach in ThorAxe, an efficient, versatile, robust, and freely available computational tool.

The clinical importance of tandem exon duplication-derived substitutions

Most coding genes in the human genome are annotated with multiple alternative transcripts. However, clear evidence for the functional relevance of the protein isoforms produced by these alternative transcripts is often hard to find. Alternative isoforms generated from tandem exon duplication-derived substitutions are an exception. These splice events are rare, but have important functional consequences. Here, we have catalogued the 236 tandem exon duplication-derived substitutions annotated in the GENCODE human reference set. We find that more than 90% of the events have a last common ancestor in teleost fish, so are at least 425 million years old, and twenty-one can be traced back to the Bilateria clade. Alternative isoforms generated from tandem exon duplication-derived substitutions also have significantly more clinical impact than other alternative isoforms. Tandem exon duplication-derived substitutions have >25 times as many pathogenic and likely pathogenic mutations as other alternative events. Tandem exon duplication-derived substitutions appear to have vital functional roles in the cell and may have played a prominent part in metazoan evolution.

A Dynamic Splicing Program Ensures Proper Synaptic Connections in the Developing Cerebellum

Tight coordination of gene expression in the developing cerebellum is crucial for establishment of neuronal circuits governing motor and cognitive function. However, transcriptional changes alone do not explain all of the switches underlying neuronal differentiation. Here we unveiled a widespread and highly dynamic splicing program that affects synaptic genes in cerebellar neurons. The motifs enriched in modulated exons implicated the splicing factor Sam68 as a regulator of this program. Sam68 controls splicing of exons with weak branchpoints by directly binding near the 3' splice site and competing with U2AF recruitment. Ablation of Sam68 disrupts splicing regulation of synaptic genes associated with neurodevelopmental diseases and impairs synaptic connections and firing of Purkinje cells, resulting in motor coordination defects, ataxia, and abnormal social behavior. These findings uncover an unexpectedly dynamic splicing regulatory network that shapes the synapse in early life and establishes motor and cognitive circuitry in the developing cerebellum.

An extended catalogue of tandem alternative splice sites in human tissue transcriptomes

Tandem alternative splice sites (TASS) is a special class of alternative splicing events that are characterized by a close tandem arrangement of splice sites. Most TASS lack functional characterization and are believed to arise from splicing noise. Based on the RNA-seq data from the Genotype Tissue Expression project, we present an extended catalogue of TASS in healthy human tissues and analyze their tissue-specific expression. The expression of TASS is usually dominated by one major splice site (maSS), while the expression of minor splice sites (miSS) is at least an order of magnitude lower. Among 46k miSS with sufficient read support, 9k (20%) are significantly expressed above the expected noise level, and among them 2.5k are expressed tissue-specifically. We found significant correlations between tissue-specific expression of RNA-binding proteins (RBP), tissue-specific expression of miSS, and miSS response to RBP inactivation by shRNA. In combination with RBP profiling by eCLIP, this allowed prediction of novel cases of tissue-specific splicing regulation including a miSS in QKI mRNA that is likely regulated by PTBP1. The analysis of human primary cell transcriptomes suggested that both tissue-specific and cell-type-specific factors contribute to the regulation of miSS expression. More than 20% of tissue-specific miSS affect structured protein regions and may adjust protein-protein interactions or modify the stability of the protein core. The significantly expressed miSS evolve under the same selection pressure as maSS, while other miSS lack signatures of evolutionary selection and conservation. Using mixture models, we estimated that not more than 15% of maSS and not more than 54% of tissue-specific miSS are noisy, while the proportion of noisy splice sites among non-significantly expressed miSS is above 63%.

Pre-mRNA splicing is an important step in the processing of the genomic information during gene expression. During splicing, introns are excised from a gene transcript, and the remaining exons are ligated. Our work concerns one its particular subtype, which involves the so-called tandem alternative splice sites, a group of closely located exon borders that are used alternatively. We analyzed RNA-seq measurements of gene expression provided by the Genotype-Tissue Expression (GTEx) project, the largest to-date collection of such measurements in healthy human tissues, and constructed a detailed catalogue of tandem alternative splice sites. Within this catalogue, we characterized patterns of tissue-specific expression, regulation, impact on protein structure, and evolutionary selection acting on tandem alternative splice sites. In a number of genes, we predicted regulatory mechanisms that could be responsible for choosing one of many tandem alternative splice sites. The results of this study provide an invaluable resource for molecular biologists studying alternative splicing.

Hepatocyte nuclear factor 1 beta: A perspective in cancer

Hepatocyte nuclear factor 1 beta (HNF1 β/B) exists as a homeobox transcription factor having a vital role in the embryonic development of organs mainly liver, kidney and pancreas. Initially described as a gene causing maturity‐onset diabetes of the young (MODY), HNF1β expression deregulation and single nucleotide polymorphisms in HNF1β have now been associated with several tumours including endometrial, prostate, ovarian, hepatocellular, renal and colorectal cancers. Its function has been studied either as homodimer or heterodimer with HNF1α. In this review, the role of HNF1B in different cancers will be discussed along with the role of its splice variants, and its emerging role as a potential biomarker in cancer.

Isoform-resolved correlation analysis between mRNA abundance regulation and protein level degradation

Profiling of biological relationships between different molecular layers dissects regulatory mechanisms that ultimately determine cellular function. To thoroughly assess the role of protein post‐translational turnover, we devised a strategy combining pulse stable isotope‐labeled amino acids in cells (pSILAC), data‐independent acquisition mass spectrometry (DIA‐MS), and a novel data analysis framework that resolves protein degradation rate on the level of mRNA alternative splicing isoforms and isoform groups. We demonstrated our approach by the genome‐wide correlation analysis between mRNA amounts and protein degradation across different strains of HeLa cells that harbor a high grade of gene dosage variation. The dataset revealed that specific biological processes, cellular organelles, spatial compartments of organelles, and individual protein isoforms of the same genes could have distinctive degradation rate. The protein degradation diversity thus dissects the corresponding buffering or concerting protein turnover control across cancer cell lines. The data further indicate that specific mRNA splicing events such as intron retention significantly impact the protein abundance levels. Our findings support the tight association between transcriptome variability and proteostasis and provide a methodological foundation for studying functional protein degradation.

Sex-Specific Selection Drives the Evolution of Alternative Splicing in Birds

Males and females of the same species share the majority of their genomes, yet they are frequently exposed to conflicting selection pressures. Gene regulation is widely assumed to resolve these conflicting sex-specific selection pressures, and although there has been considerable focus on elucidating the role of gene expression level in sex-specific adaptation, other regulatory mechanisms have been overlooked. Alternative splicing enables different transcripts to be generated from the same gene, meaning that exons which have sex-specific beneficial effects can in theory be retained in the gene product, whereas exons with detrimental effects can be skipped. However, at present, little is known about how sex-specific selection acts on broad patterns of alternative splicing. Here, we investigate alternative splicing across males and females of multiple bird species. We identify hundreds of genes that have sex-specific patterns of splicing and establish that sex differences in splicing are correlated with phenotypic sex differences. Additionally, we find that alternatively spliced genes have evolved rapidly as a result of sex-specific selection and suggest that sex differences in splicing offer another route to sex-specific adaptation when gene expression level changes are limited by functional constraints. Overall, our results shed light on how a diverse transcriptional framework can give rise to the evolution of phenotypic sexual dimorphism.

Differential contribution of transcriptomic regulatory layers in the definition of neuronal identity

Previous transcriptomic profiling studies have typically focused on separately analyzing mRNA expression, alternative splicing and alternative polyadenylation differences between cell and tissue types. However, the relative contribution of these three transcriptomic regulatory layers to cell type specification is poorly understood. This question is particularly relevant to neurons, given their extensive heterogeneity associated with brain location, morphology and function. In the present study, we generated profiles for the three regulatory layers from developmentally and regionally distinct subpopulations of neurons from the mouse hippocampus and broader nervous system. Multi-omics factor analyses revealed differing contributions of each transcriptomic layer in the discrimination of neurons based on their stage of development, region, and function. Importantly, profiles of differential alternative splicing and polyadenylation better discriminated specific neuronal subtype populations than gene expression patterns. These results provide evidence for differential relative contributions of coordinated gene regulatory layers in the specification of neuronal subtypes.

Insights Into Translatomics in the Nervous System

Most neurological disorders are caused by abnormal gene translation. Generally, dysregulation of elements involved in the translational process disrupts homeostasis in neurons and neuroglia. Better understanding of how the gene translation process occurs requires detailed analysis of transcriptomic and proteomic profile data. However, a lack of strictly direct correlations between mRNA and protein levels limits translational investigation by combining transcriptomic and proteomic profiling. The much better correlation between proteins and translated mRNAs than total mRNAs in abundance and insufficiently sensitive proteomics approach promote the requirement of advances in translatomics technology. Translatomics which capture and sequence the mRNAs associated with ribosomes has been effective in identifying translational changes by genetics or projections, ribosome stalling, local translation, and transcript isoforms in the nervous system. Here, we place emphasis on the main three translatomics methods currently used to profile mRNAs attached to ribosome-nascent chain complex (RNC-mRNA). Their prominent applications in neurological diseases including glioma, neuropathic pain, depression, fragile X syndrome (FXS), neurodegenerative disorders are outlined. The content reviewed here expands our understanding on the contributions of aberrant translation to neurological disease development.

RNA in cancer

While the processing of mRNA is essential for gene expression, recent findings have highlighted that RNA processing is systematically altered in cancer. Mutations in RNA splicing factor genes and the shortening of 3' untranslated regions are widely observed. Moreover, evidence is accumulating that other types of RNAs, including circular RNAs, can contribute to tumorigenesis. In this Review, we highlight how altered processing or activity of coding and non-coding RNAs contributes to cancer. We introduce the regulation of gene expression by coding and non-coding RNA and discuss both established roles (microRNAs and long non-coding RNAs) and emerging roles (selective mRNA processing and circular RNAs) for RNAs, highlighting the potential mechanisms by which these RNA subtypes contribute to cancer. The widespread alteration of coding and non-coding RNA demonstrates that altered RNA biogenesis contributes to multiple hallmarks of cancer.

An analysis of tissue-specific alternative splicing at the protein level

The role of alternative splicing is one of the great unanswered questions in cellular biology. There is strong evidence for alternative splicing at the transcript level, and transcriptomics experiments show that many splice events are tissue specific. It has been suggested that alternative splicing evolved in order to remodel tissue-specific protein-protein networks. Here we investigated the evidence for tissue-specific splicing among splice isoforms detected in a large-scale proteomics analysis. Although the data supporting alternative splicing is limited at the protein level, clear patterns emerged among the small numbers of alternative splice events that we could detect in the proteomics data. More than a third of these splice events were tissue-specific and most were ancient: over 95% of splice events that were tissue-specific in both proteomics and RNAseq analyses evolved prior to the ancestors of lobe-finned fish, at least 400 million years ago. By way of contrast, three in four alternative exons in the human gene set arose in the primate lineage, so our results cannot be extrapolated to the whole genome. Tissue-specific alternative protein forms in the proteomics analysis were particularly abundant in nervous and muscle tissues and their genes had roles related to the cytoskeleton and either the structure of muscle fibres or cell-cell connections. Our results suggest that this conserved tissue-specific alternative splicing may have played a role in the development of the vertebrate brain and heart.

We manually curated a set of 255 splice events detected in a large-scale tissue-based proteomics experiment and found that more than a third had evidence of significant tissue-specific differences. Events that were significantly tissue-specific at the protein level were highly conserved; almost 75% evolved over 400 million years ago. The tissues in which we found most evidence for tissue-specific splicing were nervous tissues and cardiac tissues. Genes with tissue-specific events in these two tissues had functions related to important cellular structures in brain and heart tissues. These splice events may have been essential for the development of vertebrate heart and muscle. However, our data set may not be representative of alternative exons as a whole. We found that most tissue specific splicing was strongly conserved, but just 5% of annotated alternative exons in the human gene set are ancient. More than three quarters of alternative exons are primate-derived. Although the analysis does not provide a definitive answer to the question of the functional role of alternative splicing, our results do indicate that alternative splice variants may have played a significant part in the evolution of brain and heart tissues in vertebrates.