Full-length transcriptome reconstruction reveals a large diversity of RNA and protein isoforms in rat hippocampus

Transcriptome profiling and alternative splicing analysis of skeletal muscle in Pseudocaranx dentex: insights into slow-twitch and fast-twitch muscle specialization

Alternative splicing (AS) plays a crucial role in regulating muscle type specialization characteristics in mammals but has been rarely explored in teleost. In this study, we combined Iso-Seq and RNA-Seq technologies to profile the transcriptome and AS landscape of slow-twitch muscle (SM) and fast-twitch muscle (FM) in a migratory teleost, Pseudocaranx dentex. We identified 24,096 full-length transcripts and 14,346 isoforms in SM, and 18,483 full-length transcripts and 10,541 isoforms in FM, revealing extensive transcript and isoform diversity. The 3086 differentially expressed transcripts (DETs) were found to contribute to metabolic and contractile differences between SM and FM. Additionally, we detected 5761 AS events in SM and 4543 in FM, with skipped exons (SE) and intron retention (IR) being the predominant AS types. Furthermore, 325 differentially AS genes (DASGs) were found to regulate differences in metabolic processes, organelles organization, cellular component organization, and microtubule-based processes. Importantly, transcripts of tnni2, capzb, neb, and pdlim5 produced by AS with significant expression differences between SM and FM were determined to associate with sarcomere assembly. This study provides the first comprehensive view of transcriptome complexity and splice variants in teleost skeletal muscle and sheds light on the molecular mechanisms underlying muscle type specialization through post-transcriptional regulation.

Genome-wide profiling and functional characterization of circular RNAs in neural development and injury: insights from a rat model research

Circular RNAs (circRNAs) have re-emerged as promising gene regulators in various physiological and pathological conditions. However, the expression patterns of circRNAs in the developing spinal cord of mammals and the comprehensive distribution of circRNAs across different tissues remain poorly understood. In this study, rats were used as the model organism. We conducted a comprehensive analysis of 15 RNA-Seq datasets comprising 217 rat samples and developed a web-based resource, CiRNat, to facilitate access to these data. We identified 15,251 credible circRNAs and validated them through experimental approaches. Notably, we observed two significant time points for circRNA increase during spinal cord development, approximately at embryonic day 14 (E14d) and postnatal week 4 (P4w). Analysis of circRNA expression in various rat tissues revealed higher expression levels in central nervous system tissues compared to peripheral nervous system tissues and other tissues. Furthermore, some highly abundant circRNAs exhibited tissue- and species-specific expression patterns and differed from their cognate linear RNAs, such as those derived from Gigyf2. Integrating polysome profiling and bioinformatic predictions suggested potential functions of certain circRNAs as miRNA sponges and translational templates. Collectively, this study provides the first comprehensive landscape of circRNAs in the developing spinal cord, offering an important resource and new insights for future exploration of functional circRNAs in central nervous system development and related diseases.

A systematic benchmark of Nanopore long-read RNA sequencing for transcript-level analysis in human cell lines

The human genome contains instructions to transcribe more than 200,000 RNAs. However, many RNA transcripts are generated from the same gene, resulting in alternative isoforms that are highly similar and that remain difficult to quantify. To evaluate the ability to study RNA transcript expression, we profiled seven human cell lines with five different RNA-sequencing protocols, including short-read cDNA, Nanopore long-read direct RNA, amplification-free direct cDNA and PCR-amplified cDNA sequencing, and PacBio IsoSeq, with multiple spike-in controls, and additional transcriptome-wide N6-methyladenosine profiling data. We describe differences in read length, coverage, throughput and transcript expression, reporting that long-read RNA sequencing more robustly identifies major isoforms. We illustrate the value of the SG-NEx data to identify alternative isoforms, novel transcripts, fusion transcripts and N6-methyladenosine RNA modifications. Together, the SG-NEx data provide a comprehensive resource enabling the development and benchmarking of computational methods for profiling complex transcriptional events at isoform-level resolution.

Alternative transcripts recode human genes to express overlapping, frameshifted microproteins

Overlapping genes were thought to be essentially absent from the human genome until the discovery of abundant, frameshifted internal open reading frames (iORFs) nested within annotated protein coding sequences. However, it is currently unclear how many functional human iORFs exist and how they are expressed. We demonstrate that, in hundreds of cases, alternative transcript variants that bypass the start codon of annotated coding sequences (CDSs) can recode a human gene to express the iORF-encoded microprotein. While many human genes generate such non-coding alternative transcripts, they are poorly annotated. Here we develope a new analysis pipeline enabling the assignment of translated human iORFs to alternative transcripts, and provide long-read sequencing and molecular validation of their expression in dozens of cases. Finally, we demonstrate that a conserved DEDD2 iORF switches the function of this gene from pro- to anti-apoptotic. This work thus demonstrates that alternative transcript variants can broadly reprogram human genes to express frameshifted iORFs, revealing new levels of complexity in the human transcriptome and proteome.

Identification of novel cattle (Bos taurus) genes and biological insights of their function in pre-implantation embryo development

BACKGROUND

Appropriate regulation of genes expressed in oocytes and embryos is essential for acquisition of developmental competence in mammals. Here, we hypothesized that several genes expressed in oocytes and pre-implantation embryos remain unknown. Our goal was to reconstruct the transcriptome of oocytes (germinal vesicle and metaphase II) and pre-implantation cattle embryos (blastocysts) using short-read and long-read sequences to identify putative new genes.

RESULTS

We identified 274,342 transcript sequences and 3,033 of those loci do not match a gene present in official annotations and thus are potential new genes. Notably, 63.67% (1,931/3,033) of potential novel genes exhibited coding potential. Also noteworthy, 97.92% of the putative novel genes overlapped annotation with transposable elements. Comparative analysis of transcript abundance identified that 1,840 novel genes (recently added to the annotation) or potential new genes were differentially expressed between developmental stages (FDR < 0.01). We also determined that 522 novel or potential new genes (448 and 34, respectively) were upregulated at eight-cell embryos compared to oocytes (FDR < 0.01). In eight-cell embryos, 102 novel or putative new genes were co-expressed (|r|> 0.85, P < 1 × 10-8) with several genes annotated with gene ontology biological processes related to pluripotency maintenance and embryo development. CRISPR-Cas9 genome editing confirmed that the disruption of one of the novel genes highly expressed in eight-cell embryos reduced blastocyst development (ENSBTAG00000068261, P = 1.55 × 10-7).

CONCLUSIONS

Our results revealed several putative new genes that need careful annotation. Many of the putative new genes have dynamic regulation during pre-implantation development and are important components of gene regulatory networks involved in pluripotency and blastocyst formation.

Thermoregulated transcriptomics: the molecular basis and biological significance of temperature-dependent alternative splicing

Temperature-dependent alternative splicing (AS) is a crucial mechanism for organisms to adapt to varying environmental temperatures. In mammals, even slight fluctuations in body temperature are sufficient to drive significant AS changes in a concerted manner. This dynamic regulation allows organisms to finely tune gene expression and protein isoform diversity in response to temperature cues, ensuring proper cellular function and physiological adaptation. Understanding the molecular mechanisms underlying temperature-dependent AS thus provides valuable insights into the intricate interplay between environmental stimuli and gene expression regulation. In this review, we provide an overview of recent advances in understanding temperature-regulated AS across various biological processes and systems. We will discuss the machinery sensing and translating temperature cues into changed AS patterns, the adaptation of the splicing regulatory machinery to extreme temperatures, the role of temperature-dependent AS in shaping the transcriptome, functional implications and the development of potential therapeutics targeting temperature-sensitive AS pathways.

Poly(A) tale: From A to A; RNA polyadenylation in prokaryotes and eukaryotes

Most eukaryotic mRNAs and different non-coding RNAs undergo a form of 3' end processing known as polyadenylation. Polyadenylation machinery is present in almost all organisms except few species. In bacteria, the machinery has evolved from PNPase, which adds heteropolymeric tails, to a poly(A)-specific polymerase. Differently, a complex machinery for accurate polyadenylation and several non-canonical poly(A) polymerases are developed in eukaryotes. The role of poly(A) tail has also evolved from serving as a degradative signal to a stabilizing modification that also regulates translation. In this review, we discuss poly(A) tail emergence in prokaryotes and its development into a stable, yet dynamic feature at the 3' end of mRNAs in eukaryotes. We also describe how appearance of novel poly(A) polymerases gives cells flexibility to shape poly(A) tail. We explain how poly(A) tail dynamics help regulate cognate RNA metabolism in a context-dependent manner, such as during oocyte maturation. Finally, we describe specific mRNAs in metazoans that bear stem-loops instead of poly(A) tails. We conclude with how recent discoveries about poly(A) tail can be applied to mRNA technology. This article is categorized under: RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution RNA Processing > 3' End Processing RNA Turnover and Surveillance > Regulation of RNA Stability.

Cochlear transcript diversity and its role in auditory functions implied by an otoferlin short isoform

Isoforms of a gene may contribute to diverse biological functions. In the cochlea, the repertoire of alternative isoforms remains unexplored. We integrated single-cell short-read and long-read RNA sequencing techniques and identified 236,012 transcripts, 126,612 of which were unannotated in the GENCODE database. Then we analyzed and verified the unannotated transcripts using RNA-seq, RT-PCR, Sanger sequencing, and MS-based proteomics approaches. To illustrate the importance of identifying spliced isoforms, we investigated otoferlin, a key protein involved in synaptic transmission in inner hair cells (IHCs). Upon deletion of the canonical otoferlin isoform, the identified short isoform is able to support normal hearing thresholds but with reduced sustained exocytosis of IHCs, and further revealed otoferlin functions in endocytic membrane retrieval that was not well-addressed previously. Furthermore, we found that otoferlin isoforms are associated with IHC functions and auditory phenotypes. This work expands our mechanistic understanding of auditory functions at the level of isoform resolution.

Transcriptomic and Functional Analyses Reveal the Different Roles of Vitamins C, E, and K in Regulating Viral Infections in Maize

Maize lethal necrosis (MLN), one of the most important maize viral diseases, is caused by maize chlorotic mottle virus (MCMV) infection in combination with a potyvirid, such as sugarcane mosaic virus (SCMV). However, the resistance mechanism of maize to MLN remains largely unknown. In this study, we obtained isoform expression profiles of maize after SCMV and MCMV single and synergistic infection (S + M) via comparative analysis of SMRT- and Illumina-based RNA sequencing. A total of 15,508, 7567, and 2378 differentially expressed isoforms (DEIs) were identified in S + M, MCMV, and SCMV libraries, which were primarily involved in photosynthesis, reactive oxygen species (ROS) scavenging, and some pathways related to disease resistance. The results of virus-induced gene silencing (VIGS) assays revealed that silencing of a vitamin C biosynthesis-related gene, ZmGalDH or ZmAPX1, promoted viral infections, while silencing ZmTAT or ZmNQO1, the gene involved in vitamin E or K biosynthesis, inhibited MCMV and S + M infections, likely by regulating the expressions of pathogenesis-related (PR) genes. Moreover, the relationship between viral infections and expression of the above four genes in ten maize inbred lines was determined. We further demonstrated that the exogenous application of vitamin C could effectively suppress viral infections, while vitamins E and K promoted MCMV infection. These findings provide novel insights into the gene regulatory networks of maize in response to MLN, and the roles of vitamins C, E, and K in conditioning viral infections in maize.

A multi-omics view of neuronal subcellular protein synthesis

While it has long been known that protein synthesis is necessary for long-term memory in the brain, the logistics of neuronal protein synthesis is complicated by the extensive subcellular compartmentalization of the neuron. Local protein synthesis solves many of the logistic problems posed by the extreme complexity of dendritic and axonal arbors and the huge number of synapses. Here we review recent multi-omic and quantitative studies that elaborate a systems view of decentralized neuronal protein synthesis. We highlight recent insights from the transcriptomic, translatomic, and proteomic levels, discuss the nuanced logic of local protein synthesis for different protein features, and list the missing information needed to build a comprehensive logistic model for neuronal protein supply.

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.

Integrative analysis of Iso-Seq and RNA-seq reveals dynamic changes of alternative promoter, alternative splicing and alternative polyadenylation during Angiotensin II-induced senescence in rat primary aortic endothelial cells

In eukaryotes, alternative promoter (AP), alternative splicing (AS), and alternative polyadenylation (APA) are three crucial regulatory mechanisms that modulate message RNA (mRNA) diversity. Although AP, AS and APA are involved in diverse biological processess, whether they have dynamic changes in Angiotensin II (Ang II) induced senescence in rat primary aortic endothelial cells (RAECs), an important cellular model for studying cardiovascular disease, remains unclear. Here we integrated both PacBio single-molecule long-read isoform sequencing (Iso-Seq) and Illumina short-read RNA sequencing (RNA-seq) to analyze the changes of AP, AS and APA in Ang II-induced senescent RAECs. Iso-Seq generated 36,278 isoforms from 10,145 gene loci and 65.81% of these isoforms are novel, which were further cross-validated by public data obtained by other techonologies such as CAGE, PolyA-Seq and 3'READS. APA contributed most to novel isoforms, followed by AS and AP. Further investigation showed that AP, AS and APA could all contribute to the regulation of isoform, but AS has more dynamic changes compared to AP and APA upon Ang II stimulation. Genes undergoing AP, AS and APA in Ang II-treated cells are enriched in various pathways related to aging or senescence, suggesting that these molecular changes are involved in functional alterations during Ang II-induced senescence. Together, the present study largely improved the annotation of rat genome and revealed gene expression changes at isoform level, extending the understanding of the complexity of gene regulation in Ang II-treated RAECs, and also provided novel clues for discovering the regulatory mechanism undelying Ang II caused vascular senescence and diseases.

Genome assembly and isoform analysis of a highly heterozygous New Zealand fisheries species, the tarakihi (Nemadactylus macropterus)

Although being some of the most valuable and heavily exploited wild organisms, few fisheries species have been studied at the whole-genome level. This is especially the case in New Zealand, where genomics resources are urgently needed to assist fisheries management. Here, we generated 55 Gb of short Illumina reads (92× coverage) and 73 Gb of long Nanopore reads (122×) to produce the first genome assembly of the marine teleost tarakihi [Nemadactylus macropterus (Forster, 1801)], a highly valuable fisheries species in New Zealand. An additional 300 Mb of Iso-Seq reads were obtained to assist in gene annotation. The final genome assembly was 568 Mb long with an N50 of 3.37 Mb. The genome completeness was high, with 97.8% of complete Actinopterygii Benchmarking Universal Single-Copy Orthologs. Heterozygosity values estimated through k-mer counting (1.00%) and bi-allelic SNPs (0.64%) were high compared with the same values reported for other fishes. Iso-Seq analysis recovered 91,313 unique transcripts from 15,515 genes (mean ratio of 5.89 transcripts per gene), and the most common alternative splicing event was intron retention. This highly contiguous genome assembly and the isoform-resolved transcriptome will provide a useful resource to assist the study of population genomics and comparative eco-evolutionary studies in teleosts and related organisms.

Context-specific regulation and function of mRNA alternative polyadenylation

Alternative cleavage and polyadenylation (APA) is a widespread mechanism to generate mRNA isoforms with alternative 3' untranslated regions (UTRs). The expression of alternative 3' UTR isoforms is highly cell type specific and is further controlled in a gene-specific manner by environmental cues. In this Review, we discuss how the dynamic, fine-grained regulation of APA is accomplished by several mechanisms, including cis-regulatory elements in RNA and DNA and factors that control transcription, pre-mRNA cleavage and post-transcriptional processes. Furthermore, signalling pathways modulate the activity of these factors and integrate APA into gene regulatory programmes. Dysregulation of APA can reprogramme the outcome of signalling pathways and thus can control cellular responses to environmental changes. In addition to the regulation of protein abundance, APA has emerged as a major regulator of mRNA localization and the spatial organization of protein synthesis. This role enables the regulation of protein function through the addition of post-translational modifications or the formation of protein-protein interactions. We further discuss recent transformative advances in single-cell RNA sequencing and CRISPR-Cas technologies, which enable the mapping and functional characterization of alternative 3' UTRs in any biological context. Finally, we discuss new APA-based RNA therapeutics, including compounds that target APA in cancer and therapeutic genome editing of degenerative diseases.

Comparative transcriptome analysis reveals regional specialization of gene expression in larval silkworm (Bombyx mori) midgut

Insect midgut plays a central role in food digestion and nutrition absorption. Larval silkworm midgut could be divided into 3 distinct regions based on their morphological colors. However, it remains rudimentary of regional gene expression and physiological function in larval silkworm midgut. Through transcriptome sequencing of 3 midgut compartments, a comprehensive analysis of gene expression atlas along the anterior-posterior axis was conducted. Posterior midgut was found transcriptionally divergent from anterior and middle midgut. Differentially expressed gene analysis revealed the regional specialization of digestive enzyme production, transmembrane transport, chitin metabolism, and hormone regulation in different midgut regions. In addition, gene subsets of pan-midgut and region-specific transcription factors (TFs) along the length of midgut were also identified. The results suggested that homeobox TFs might play an essential role in transcriptional variations across the midgut. Altogether, our study provided the first fundamental resource to investigate physiological function and regulation mechanism in larval midgut compartmentalization.

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.

Targeted proteoform mapping uncovers specific Neurexin-3 variants required for dendritic inhibition

The diversification of cell adhesion molecules by alternative splicing is proposed to underlie molecular codes for neuronal wiring. Transcriptomic approaches mapped detailed cell-type-specific mRNA splicing programs. However, it has been hard to probe the synapse-specific localization and function of the resulting protein splice isoforms, or “proteoforms,” in vivo. We here apply a proteoform-centric workflow in mice to test the synapse-specific functions of the splice isoforms of the synaptic adhesion molecule Neurexin-3 (NRXN3). We uncover a major proteoform, NRXN3 AS5, that is highly expressed in GABAergic interneurons and at dendrite-targeting GABAergic terminals. NRXN3 AS5 abundance significantly diverges from Nrxn3 mRNA distribution and is gated by translation-repressive elements. Nrxn3 AS5 isoform deletion results in a selective impairment of dendrite-targeting interneuron synapses in the dentate gyrus without affecting somatic inhibition or glutamatergic perforant-path synapses. This work establishes cell- and synapse-specific functions of a specific neurexin proteoform and highlights the importance of alternative splicing regulation for synapse specification.

•

Translational regulation guides alternative Neurexin proteoform expression

•

NRXN3 AS5 proteoforms are concentrated at dendrite-targeting interneuron synapses

•

A proteome-centric workflow uncovers NRXN3 AS5 interactors in vivo

•

Loss of NRXN3 AS5 leads to selective impairments in dendritic inhibition

Generation and network analysis of an RNA-seq transcriptional atlas for the rat

The laboratory rat is an important model for biomedical research. To generate a comprehensive rat transcriptomic atlas, we curated and downloaded 7700 rat RNA-seq datasets from public repositories, downsampled them to a common depth and quantified expression. Data from 585 rat tissues and cells, averaged from each BioProject, can be visualized and queried at http://biogps.org/ratatlas. Gene co-expression network (GCN) analysis revealed clusters of transcripts that were tissue or cell type restricted and contained transcription factors implicated in lineage determination. Other clusters were enriched for transcripts associated with biological processes. Many of these clusters overlap with previous data from analysis of other species, while some (e.g. expressed specifically in immune cells, retina/pineal gland, pituitary and germ cells) are unique to these data. GCN analysis on large subsets of the data related specifically to liver, nervous system, kidney, musculoskeletal system and cardiovascular system enabled deconvolution of cell type-specific signatures. The approach is extensible and the dataset can be used as a point of reference from which to analyse the transcriptomes of cell types and tissues that have not yet been sampled. Sets of strictly co-expressed transcripts provide a resource for critical interpretation of single-cell RNA-seq data.

Long read sequencing reveals novel isoforms and insights into splicing regulation during cell state changes

BACKGROUND

Alternative splicing is a key mechanism underlying cellular differentiation and a driver of complexity in mammalian neuronal tissues. However, understanding of which isoforms are differentially used or expressed and how this affects cellular differentiation remains unclear. Long read sequencing allows full-length transcript recovery and quantification, enabling transcript-level analysis of alternative splicing processes and how these change with cell state. Here, we utilise Oxford Nanopore Technologies sequencing to produce a custom annotation of a well-studied human neuroblastoma cell line SH-SY5Y, and to characterise isoform expression and usage across differentiation.

RESULTS

We identify many previously unannotated features, including a novel transcript of the voltage-gated calcium channel subunit gene, CACNA2D2. We show differential expression and usage of transcripts during differentiation identifying candidates for future research into state change regulation.

CONCLUSIONS

Our work highlights the potential of long read sequencing to uncover previously unknown transcript diversity and mechanisms influencing alternative splicing.

Improved Reference Genome Annotation of Brassica rapa by Pacific Biosciences RNA Sequencing

The species Brassica rapa includes several important vegetable crops. The draft reference genome of B. rapa ssp. pekinensis was completed in 2011, and it has since been updated twice. The pangenome with structural variations of 18 B. rapa accessions was published in 2021. Although extensive genomic analysis has been conducted on B. rapa, a comprehensive genome annotation including gene structure, alternative splicing (AS) events, and non-coding genes is still lacking. Therefore, we used the Pacific Biosciences (PacBio) single-molecular long-read technology to improve gene models and produced the annotated genome version 3.5. In total, we obtained 753,041 full-length non-chimeric (FLNC) reads and collapsed these into 92,810 non-redundant consensus isoforms, capturing 48% of the genes annotated in the B. rapa reference genome annotation v3.1. Based on the isoform data, we identified 830 novel protein-coding genes that were missed in previous genome annotations, defined the untranslated regions (UTRs) of 20,340 annotated genes and corrected 886 wrongly spliced genes. We also identified 28,564 AS events and 1,480 long non-coding RNAs (lncRNAs). We produced a relatively complete and high-quality reference transcriptome for B. rapa that can facilitate further functional genomic research.

UNAGI: Yeast Transcriptome Reconstruction and Gene Discovery Using Nanopore Sequencing

Computational approaches are the main approaches used in genome annotation. However, accuracy is low. Untranslated regions are not identified, complex isoforms are not predicted correctly and discovery rate of noncoding RNA is low. RNA-seq has revolutionized transcriptome reconstruction over the last decade. However, fragmentation included in cDNA sequencing leads to information loss, requiring transcripts to be assembled and reconstructed, thus affecting the accuracy of reconstructed transcriptome. Recently, long-read sequencing has been introduced with technologies such as Oxford Nanopore sequencing. cDNA is sequenced directly without fragmentation producing long reads that don't need to be assembled keeping the transcript structure intact and increasing the accuracy of transcriptome reconstruction.Here we present a protocol and a pipeline to reconstruct the transcriptome of compact genomes including yeasts. It involves generating full-length cDNA and using Oxford Nanopore ligation-based sequencing kit to sequence multiple samples in the same run. The pipeline (1) strands the generated long reads, (2) corrects the reads by mapping them to the reference genome, (3) identifies transcripts including 5'UTR and 3'UTR, (4) profiles the isoforms, filtering out artifacts resulting from low accuracy in sequencing, and (5) improves accuracy of provided annotations. Using long reads improves the accuracy of transcriptome reconstruction and helps in discovering a significant number of novel RNAs.

Foster thy young: enhanced prediction of orphan genes in assembled genomes

Proteins encoded by newly-emerged genes ('orphan genes') share no sequence similarity with proteins in any other species. They provide organisms with a reservoir of genetic elements to quickly respond to changing selection pressures. Here, we systematically assess the ability of five gene prediction pipelines to accurately predict genes in genomes according to phylostratal origin. BRAKER and MAKER are existing, popular ab initio tools that infer gene structures by machine learning. Direct Inference is an evidence-based pipeline we developed to predict gene structures from alignments of RNA-Seq data. The BIND pipeline integrates ab initio predictions of BRAKER and Direct inference; MIND combines Direct Inference and MAKER predictions. We use highly-curated Arabidopsis and yeast annotations as gold-standard benchmarks, and cross-validate in rice. Each pipeline under-predicts orphan genes (as few as 11 percent, under one prediction scenario). Increasing RNA-Seq diversity greatly improves prediction efficacy. The combined methods (BIND and MIND) yield best predictions overall, BIND identifying 68% of annotated orphan genes, 99% of ancient genes, and give the highest sensitivity score regardless dataset in Arabidopsis. We provide a light weight, flexible, reproducible, and well-documented solution to improve gene prediction.

Full-Length Transcriptome Sequencing: An Insight Into the Dog Model of Heart Failure

Heart failure (HF) leads to a progressive increase in morbidity and mortality rates. This study aimed to explore the transcriptional landscape during HF and identify differentially expressed transcripts (DETs) and alternative splicing events associated with HF. We generated a dog model of HF (n = 3) using right ventricular pacemaker implantation. We performed full-length transcriptome sequencing (based on nanopore platform) on the myocardial tissues and analyzed the transcripts using differential expression analysis and functional annotation methods [Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses]. Additionally, we estimated the expression of the selected genes by quantitative real-time PCR (qRT-PCR) and detected the proportion of immune cells using flow cytometry. We found that increased B-type natriuretic peptide reduced ejection fraction, and apparent clinical signs were observed in the dog model of HF. We identified 67,458 transcripts using full-length transcriptome sequencing. A total of 785 DETs were obtained from the HF and control groups. These DETs were mainly enriched in the immune responses, especially Th1, Th2, and Th17 cell differentiation processes. Furthermore, flow cytometry results revealed that the proportion of Th1 and Th17 cells increased in patients with HF compared to controls, while the proportion of Th2 cells decreased. Differentially expressed genes in the HF and control groups associated with Th1, Th2, and Th17 cell differentiation were quantified using qRT-PCR. We also identified variable splicing events of sarcomere genes (e.g., MYBPC3, TNNT2, TTN, FLNC, and TTNI3). In addition, we detected 4,892 transcription factors and 406 lncRNAs associated with HF. Our analysis based on full-length transcript sequencing provided an analysis perspective in a dog model of HF, which is valuable for molecular research in an increasingly relevant large animal model of HF.

Probing Isoform Switching Events in Various Cancer Types: Lessons From Pan-Cancer Studies

Alternative splicing is an essential regulatory mechanism for gene expression in mammalian cells contributing to protein, cellular, and species diversity. In cancer, alternative splicing is frequently disturbed, leading to changes in the expression of alternatively spliced protein isoforms. Advances in sequencing technologies and analysis methods led to new insights into the extent and functional impact of disturbed alternative splicing events. In this review, we give a brief overview of the molecular mechanisms driving alternative splicing, highlight the function of alternative splicing in healthy tissues and describe how alternative splicing is disrupted in cancer. We summarize current available computational tools for analyzing differential transcript usage, isoform switching events, and the pathogenic impact of cancer-specific splicing events. Finally, the strategies of three recent pan-cancer studies on isoform switching events are compared. Their methodological similarities and discrepancies are highlighted and lessons learned from the comparison are listed. We hope that our assessment will lead to new and more robust methods for cancer-specific transcript detection and help to produce more accurate functional impact predictions of isoform switching events.

Pervasive 3'-UTR Isoform Switches During Mouse Oocyte Maturation

Oocyte maturation is the foundation for developing healthy individuals of mammals. Upon germinal vesicle breakdown, oocyte meiosis resumes and the synthesis of new transcripts ceases. To quantitatively profile the transcriptomic dynamics after meiotic resumption throughout the oocyte maturation, we generated transcriptome sequencing data with individual mouse oocytes at three main developmental stages: germinal vesicle (GV), metaphase I (MI), and metaphase II (MII). When clustering the sequenced oocytes, results showed that isoform-level expression analysis outperformed gene-level analysis, indicating isoform expression provided extra information that was useful in distinguishing oocyte stages. Comparing transcriptomes of the oocytes at the GV stage and the MII stage, in addition to identification of differentially expressed genes (DEGs), we detected many differentially expressed transcripts (DETs), some of which came from genes that were not identified as DEGs. When breaking down the isoform-level changes into alternative RNA processing events, we found the main source of isoform composition changes was the alternative usage of polyadenylation sites. With detailed analysis focusing on the alternative usage of 3′-UTR isoforms, we identified, out of 3,810 tested genes, 512 (13.7%) exhibiting significant switches of 3′-UTR isoforms during the process of moues oocyte maturation. Altogether, our data and analyses suggest the importance of examining isoform abundance changes during oocyte maturation, and further investigation of the pervasive 3′-UTR isoform switches in the transition may deepen our understanding on the molecular mechanisms underlying mammalian early development.

Hybrid Sequencing in Different Types of Goat Skeletal Muscles Reveals Genes Regulating Muscle Development and Meat Quality

Domestic goats are commonly reared for meat and milk production in several regions of the world. However, the genetic mechanism underlying muscle development and meat quality of goats is limited. Therefore, the aim of this study was to identify known and novel genes regulating muscle development and meat quality of goats using second- and third-generation sequencing technologies. To achieve this, the meat quality and transcriptomes of longissimus dorsi (LD) and biceps femoris (BF) muscle tissues of Lingqiu Greyback goats were examined and compared. Differentially expressed genes (DEGs) and isoforms (DEIs) were functionally annotated. Results showed that 45,574 full-length transcripts covering 18,491 loci were characterized, and 12,566 genes were co-expressed in all samples. Differential expression analysis identified 231 DEGs, including 45 novel genes in the LD and BF muscles of the goats. Additionally, 1173 DEIs were found, in which 642 novel isoforms were identified in this study. Functional annotation and pathway analysis of the DEGs and DEIs revealed that some of them were associated with muscle growth and lipid metabolism. Overall, the findings of this study contribute to the understanding of the transcriptomic diversity underlying meat quality and muscle development of goat.

Linking transport and translation of mRNAs with endosomes and mitochondria

In eukaryotic cells, proteins are targeted to their final subcellular locations with precise timing. A key underlying mechanism is the active transport of cognate mRNAs, which in many systems can be linked intimately to membrane trafficking. A prominent example is the long-distance endosomal transport of mRNAs and their local translation. Here, we describe current highlights of fundamental mechanisms of the underlying transport process as well as of biological functions ranging from endosperm development in plants to fungal pathogenicity and neuronal processes. Translation of endosome-associated mRNAs often occurs at the cytoplasmic surface of endosomes, a process that is needed for membrane-assisted formation of heteromeric protein complexes and for accurate subcellular targeting of proteins. Importantly, endosome-coupled translation of mRNAs encoding mitochondrial proteins, for example, seems to be particularly important for efficient organelle import and for regulating subcellular mitochondrial activity. In essence, these findings reveal a new mechanism of loading newly synthesised proteins onto endocytic membranes enabling intimate crosstalk between organelles. The novel link between endosomes and mitochondria adds an inspiring new level of complexity to trafficking and organelle biology.

A Functional Dissection of the mRNA and Locally Synthesized Protein Population in Neuronal Dendrites and Axons

Neurons are characterized by a complex morphology that enables the generation of subcellular compartments with unique biochemical and biophysical properties, such as dendrites, axons, and synapses. To sustain these different compartments and carry a wide array of elaborate operations, neurons express a diverse repertoire of gene products. Extensive regulation at both the messenger RNA (mRNA) and protein levels allows for the differentiation of subcellular compartments as well as numerous forms of plasticity in response to variable stimuli. Among the multiple mechanisms that control cellular functions, mRNA translation is manipulated by neurons to regulate where and when a protein emerges. Interestingly, transcriptomic and translatomic profiles of both dendrites and axons have revealed that the mRNA population only partially predicts the local protein population and that this relation significantly varies between different gene groups. Here, we describe the space that local translation occupies within the large molecular and regulatory complexity of neurons, in contrast to other modes of regulation. We then discuss the specialized organization of mRNAs within different neuronal compartments, as revealed by profiles of the local transcriptome. Finally, we discuss the features and functional implications of both locally correlated-and anticorrelated-mRNA-protein relations both under baseline conditions and during synaptic plasticity. Expected final online publication date for the Annual Review of Genetics, Volume 55 is November 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Transposable element sequence fragments incorporated into coding and noncoding transcripts modulate the transcriptome of human pluripotent stem cells

Transposable elements (TEs) occupy nearly 40% of mammalian genomes and, whilst most are fragmentary and no longer capable of transposition, they can nevertheless contribute to cell function. TEs within genes transcribed by RNA polymerase II can be copied as parts of primary transcripts; however, their full contribution to mature transcript sequences remains unresolved. Here, using long and short read (LR and SR) RNA sequencing data, we show that 26% of coding and 65% of noncoding transcripts in human pluripotent stem cells (hPSCs) contain TE-derived sequences. Different TE families are incorporated into RNAs in unique patterns, with consequences to transcript structure and function. The presence of TE sequences within a transcript is correlated with TE-type specific changes in its subcellular distribution, alterations in steady-state levels and half-life, and differential association with RNA Binding Proteins (RBPs). We identify hPSC-specific incorporation of endogenous retroviruses (ERVs) and LINE:L1 into protein-coding mRNAs, which generate TE sequence-derived peptides. Finally, single cell RNA-seq reveals that hPSCs express ERV-containing transcripts, whilst differentiating subpopulations lack ERVs and express SINE and LINE-containing transcripts. Overall, our comprehensive analysis demonstrates that the incorporation of TE sequences into the RNAs of hPSCs is more widespread and has a greater impact than previously appreciated.

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.

Assessing the functional relevance of splice isoforms

Alternative splicing of messenger RNA can generate an array of mature transcripts, but it is not clear how many go on to produce functionally relevant protein isoforms. There is only limited evidence for alternative proteins in proteomics analyses and data from population genetic variation studies indicate that most alternative exons are evolving neutrally. Determining which transcripts produce biologically important isoforms is key to understanding isoform function and to interpreting the real impact of somatic mutations and germline variations. Here we have developed a method, TRIFID, to classify the functional importance of splice isoforms. TRIFID was trained on isoforms detected in large-scale proteomics analyses and distinguishes these biologically important splice isoforms with high confidence. Isoforms predicted as functionally important by the algorithm had measurable cross species conservation and significantly fewer broken functional domains. Additionally, exons that code for these functionally important protein isoforms are under purifying selection, while exons from low scoring transcripts largely appear to be evolving neutrally. TRIFID has been developed for the human genome, but it could in principle be applied to other well-annotated species. We believe that this method will generate valuable insights into the cellular importance of alternative splicing.

Phenotypic Differences between the Alzheimer's Disease-Related hAPP-J20 Model and Heterozygous Zbtb20 Knock-Out Mice

Diverse gene products contribute to the pathogenesis of Alzheimer’s disease (AD). Experimental models have helped elucidate their mechanisms and impact on brain functions. Human amyloid precursor protein (hAPP) transgenic mice from line J20 (hAPP-J20 mice) are widely used to simulate key aspects of AD. However, they also carry an insertional mutation in noncoding sequence of one Zbtb20 allele, a gene involved in neural development. We demonstrate that heterozygous hAPP-J20 mice have reduced Zbtb20 expression in some AD-relevant brain regions, but not others, and that Zbtb20 levels are higher in hAPP-J20 mice than heterozygous Zbtb20 knock-out (Zbtb20^+/–) mice. Whereas hAPP-J20 mice have premature mortality, severe deficits in learning and memory, other behavioral alterations, and prominent nonconvulsive epileptiform activity, Zbtb20^+/– mice do not. Thus, the insertional mutation in hAPP-J20 mice does not ablate the affected Zbtb20 allele and is unlikely to account for the AD-like phenotype of this model.

The genome of the European estuarine calanoid copepod Eurytemora affinis: Potential use in molecular ecotoxicology

In this study, we sequenced and assembled the genome of a European estuarine calanoid copepod using Oxford Nanopore PromethION and Illumina HiSeq 2500 platforms. The length of the assembled genome was 776.1 Mb with N50 = 474.9 kb (BUSCO 85.9%), and the genome consisted of 2473 contigs. A total of 18,014 genes were annotated and orthologous gene clusters were analyzed in comparison to other copepods. In addition, genome-wide identification of cytochrome P450s, glutathione S-transferases, and ATP-binding cassette transporters in E. affinis was performed to determine gene repertoire of these detoxification-related gene families. Results revealed the presence of species-specific gene inventories, indicating that these gene families have evolved through species-specific gene loss/expansion processes, possibly due to adaptation to different environmental stressors. Our study provides a new inventory of the European estuarine calanoid copepod E. affinis genome with emphasis on phase I, II, and III detoxification systems.

A Novel, Universally Active C-terminal Protein Degradation Signal Generated by Alternative Splicing

Proteome integrity is crucial for cellular homeostasis and adaptation to stress conditions such as hypoxia. One mechanism for rapid adaptation of the proteome in response to changing environmental signals is alternative splicing. In addition to generating different protein isoforms, alternative splicing is also capable of controlling total protein levels by the regulated synthesis of non-productive mRNA isoforms. The hypoxia-induced isoform E of the tumor suppressor MAX is produced by retention and translation of the last intron. This leads to an alternative C-terminus that harbors a potent C-degron, the isoE degron. Strikingly, the isoE degron represents a universal protein degradation signal that is not only functional in mammalian cells, but also in yeast and even in bacteria. Essential for efficient protein decay is a conserved (F/W)xxW motif. Degradation of isoE tagged proteins is mediated by the proteasome in eukaryotes and Lon protease in bacteria. Thus, the isoE degron is a broadly applicable and highly efficient tool in protein analyses.

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.

Subcellular sequencing of single neurons reveals the dendritic transcriptome of GABAergic interneurons

Although mRNAs are localized in the processes of excitatory neurons, it is still unclear whether interneurons also localize a large population of mRNAs. In addition, the variability in the localized mRNA population within and between cell types is unknown. Here we describe the unbiased transcriptomic characterization of the subcellular compartments of hundreds of single neurons. We separately profiled the dendritic and somatic transcriptomes of individual rat hippocampal neurons and investigated mRNA abundances in the soma and dendrites of single glutamatergic and GABAergic neurons. We found that, like their excitatory counterparts, interneurons contain a rich repertoire of ~4000 mRNAs. We observed more cell type-specific features among somatic transcriptomes than their associated dendritic transcriptomes. Finally, using celltype-specific metabolic labeling of isolated neurites, we demonstrated that the processes of glutamatergic and, notably, GABAergic neurons were capable of local translation, suggesting mRNA localization and local translation are general properties of neurons.

A Picture Worth a Thousand Molecules-Integrative Technologies for Mapping Subcellular Molecular Organization and Plasticity in Developing Circuits

A key challenge in developmental neuroscience is identifying the local regulatory mechanisms that control neurite and synaptic refinement over large brain volumes. Innovative molecular techniques and high-resolution imaging tools are beginning to reshape our view of how local protein translation in subcellular compartments drives axonal, dendritic, and synaptic development and plasticity. Here we review recent progress in three areas of neurite and synaptic study in situ-compartment-specific transcriptomics/translatomics, targeted proteomics, and super-resolution imaging analysis of synaptic organization and development. We discuss synergies between sequencing and imaging techniques for the discovery and validation of local molecular signaling mechanisms regulating synaptic development, plasticity, and maintenance in circuits.

The PacBio Full-Length Transcriptome of the Tea Aphid as a Reference Resource

The tea aphid, Aphis aurantii, has become one of the destructive pests in tea plantations in the tropics and subtropics. Very few functional studies have so far focused on the developmental and reproductive biology at a molecular level, because of the lack of comprehensive genetic information. Full-length transcriptomes represent a very highly efficient approach to obtain reference gene sequences in non-model insects. In the present study, the transcriptome of A. aurantii was comprehensively sequenced using PacBio Iso-Seq technology. A total of 46.8 Gb nucleotides and 15,938 non-redundant full-length transcripts were obtained, 13,498 (84.69%) of which were annotated into seven databases. Of these transcripts, 2,029 alternative splicing events and 15,223 simple sequence repeats were detected. Among these transcripts, 4,571 (28.68%) and 11,367 (71.32%) were long non-coding RNAs (lncRNAs) and protein-coding genes, respectively. Five hundred and ninety transcription factors were detected. The first full-length transcriptome represents a significant increase in the known genetic information of A. aurantii. It will assist the future functional study of genes involved in its development and reproduction.

Improving the diversity of captured full-length isoforms using a normalized single-molecule RNA-sequencing method

Human genes form a large variety of isoforms after transcription, encoding distinct transcripts to exert different functions. Single-molecule RNA sequencing facilitates accurate identification of the isoforms by extending nucleotide read length significantly. However, the gene or isoform diversity is lowly represented by the mRNA molecules captured by single-molecule RNA sequencing. Here, we show that a cDNA normalization procedure before the library preparation for PacBio RS II sequencing captures 3.2–6.0 fold more full-length high-quality isoform species for different human samples, as compared to the non-normalized capture procedure. Many lowly expressed, functionally important isoforms can be detected. In addition, normalized PacBio RNA sequencing also resolves more allele-specific haplotype transcripts. Finally, we apply the cDNA normalization based long-read RNA sequencing method to profile the transcriptome of human gastric signet-ring cell carcinomas, identify new cancer-specific transcriptome signatures, and thus, bring out the utility of the improved protocols in gene expression studies.

Hu et al. combine cDNA normalization before library preparation with a software tool to increase the capture of RNA isoform species in single-molecule RNA sequencing. They demonstrate that this approach can detect previously unknown transcripts in gastric signet-ring cell carcinomas that are not present in non-malignant tissue.

Splicing-accessible coding 3'UTRs control protein stability and interaction networks

BACKGROUND

3'-Untranslated regions (3'UTRs) play crucial roles in mRNA metabolism, such as by controlling mRNA stability, translation efficiency, and localization. Intriguingly, in some genes the 3'UTR is longer than their coding regions, pointing to additional, unknown functions. Here, we describe a protein-coding function of 3'UTRs upon frameshift-inducing alternative splicing in more than 10% of human and mouse protein-coding genes.

RESULTS

3'UTR-encoded amino acid sequences show an enrichment of PxxP motifs and lead to interactome rewiring. Furthermore, an elevated proline content increases protein disorder and reduces protein stability, thus allowing splicing-controlled regulation of protein half-life. This could also act as a surveillance mechanism for erroneous skipping of penultimate exons resulting in transcripts that escape nonsense mediated decay. The impact of frameshift-inducing alternative splicing on disease development is emphasized by a retinitis pigmentosa-causing mutation leading to translation of a 3'UTR-encoded, proline-rich, destabilized frameshift-protein with altered protein-protein interactions.

CONCLUSIONS

We describe a widespread, evolutionarily conserved mechanism that enriches the mammalian proteome, controls protein expression and protein-protein interactions, and has important implications for the discovery of novel, potentially disease-relevant protein variants.