Different levels of alternative splicing among eukaryotes

RNA splicing controls organ-wide maturation of postnatal heart in mice

Postnatal cardiac development requires the orchestrated maturation of diverse cellular components for which unifying control mechanisms are still lacking. Using full-length sequencing, we examined the transcriptomic landscape of the maturating mouse heart (E18.5-P28) at single-cell and transcript isoform resolution. We identified dynamically changing intercellular networks as a molecular basis of the maturing heart and alternative splicing (AS) as a common mechanism that distinguished developmental age. Manipulation of RNA-binding proteins (RBPs) remodeled the AS landscape, cardiac cell maturation, and intercellular communication through direct binding of splice targets, which were enriched for functions related to general, as well as cell-type-specific, maturation. Overexpression of an RBP nuclear cap-binding protein subunit 2 (NCBP2) in neonatal hearts repressed cardiac maturation. Together, our data suggest AS regulation by RBPs as an organ-level control mechanism in mammalian postnatal cardiac development and provide insight into the possibility of manipulating RBPs for therapeutic purposes.

Phylogenetic Analysis of 590 Species Reveals Distinct Evolutionary Patterns of Intron-Exon Gene Structures Across Eukaryotic Lineages

Introns are highly prevalent in most eukaryotic genomes. Despite the accumulating evidence for benefits conferred by the possession of introns, their specific roles and functions, as well as the processes shaping their evolution, are still only partially understood. Here, we explore the evolution of the eukaryotic intron-exon gene structure by focusing on several key features such as the intron length, the number of introns, and the intron-to-exon length ratio in protein-coding genes. We utilize whole-genome data from 590 species covering the main eukaryotic taxonomic groups and analyze them within a statistical phylogenetic framework. We found that the basic gene structure differs markedly among the main eukaryotic groups, with animals, and particularly chordates, displaying intron-rich genes, compared with plants and fungi. Reconstruction of gene structure evolution suggests that these differences evolved prior to the divergence of the main phyla and have remained mostly conserved within groups. We revisit the previously reported association between the genome size and the mean intron length and report that this association differs considerably among phyla. Analyzing a large and diverse dataset of species with whole-genome information while applying advanced modeling techniques allowed us to obtain a global evolutionary perspective. Our findings may indicate that introns play different molecular and evolutionary roles in different organisms.

Augmenting Transcriptome Annotations through the Lens of Splicing Evolution

Alternative splicing (AS) is a ubiquitous mechanism in eukaryotes. It is estimated that 90% of human genes are alternatively spliced. Despite enormous efforts, transcriptome annotations remain, nevertheless, incomplete. Conventional means of annotation were largely driven by experimental data such as RNA-seq and protein sequences, while little insight was shed on understanding transcriptomes and alternative splicings from the perspective of evolution. This study addresses this critical gap by presenting TENNIS (Transcript EvolutioN for New Isoform Splicing), an evolution-based model to predict unannotated isoforms and refine existing annotations without requiring additional data. The model of TENNIS is based on two minimal premises-AS isoforms evolve sequentially from existing isoforms, and each evolutionary step involves a single AS event. We formulate the identification of missing transcripts as an optimization problem and parsimoniously find the minimal number of novel transcripts. Our analysis showed approximately 80% of multi-transcript groups from six transcriptome annotations satisfy our evolutionary model. At a high confidence level, 40% of isoforms predicted by TENNIS were validated by deep long-read RNA-seq. In a simulated incomplete annotation scenario, TENNIS dramatically outperforms two randomized baseline approaches by a 2.25-3 fold-change in precision or a 3.5-3.9 fold-change in recall, after controlling the same level of recall or precision of the baseline methods. These results demonstrate that TENNIS effectively identifies missing transcripts by complying with minimal propositions, offering a powerful approach for transcriptome augmentations through the lens of alternative splicing evolutions. TENNIS is freely available at https://github.com/Shao-Group/tennis .

Functional analysis of the GPAT4 gene mutation predicted to affect splicing

The GPAT4 gene is considered as a potential functional candidate for single nucleotide polymorphism (SNP) studies in dairy cattle breeding due to its association with dairy performance in cattle by encoding an enzyme responsible for the presence of diacylglycerols and triacylglycerols in milk. Using the example of the GPAT4 gene, we applied the minigene splicing assay to analyze the functional consequences of its variant that was predicted to affect normal splicing. The results of functional analysis revealed the sequence variations (rs442541537), transfection experiments in a wild type and mutant cell line model system demonstrated that the investigated mutation in the second intron of the GPAT4 gene was responsible for the presence of a second exon in mature messenger RNA (mRNA). The cases of its absence in the spliced mature mRNA transcript resulted in a truncated dysfunctional protein due to the appearance of a stop codon. Thus, the discovered SNP led to alternative splicing in pre-mRNA by the 'cassette exon' ('exon skipping') mechanism. The studied mutation can potentially be a molecular genetic marker for alternative splicing for the GPAT4 gene and, therefore contributes to economic benefits in cattle breeding programs.

GWAS of multiple neuropathology endophenotypes identifies new risk loci and provides insights into the genetic risk of dementia

Genome-wide association studies (GWAS) have identified >80 Alzheimer's disease and related dementias (ADRD)-associated genetic loci. However, the clinical outcomes used in most previous studies belie the complex nature of underlying neuropathologies. Here we performed GWAS on 11 ADRD-related neuropathology endophenotypes with participants drawn from the following three sources: the National Alzheimer's Coordinating Center, the Religious Orders Study and Rush Memory and Aging Project, and the Adult Changes in Thought study (n = 7,804 total autopsied participants). We identified eight independent significantly associated loci, of which four were new (COL4A1, PIK3R5, LZTS1 and APOC2). Separately testing known ADRD loci, 19 loci were significantly associated with at least one neuropathology after false-discovery rate adjustment. Genetic colocalization analyses identified pleiotropic effects and quantitative trait loci. Methylation in the cerebral cortex at two sites near APOC2 was associated with cerebral amyloid angiopathy. Studies that include neuropathology endophenotypes are an important step in understanding the mechanisms underlying genetic ADRD risk.

RNA-Binding Protein-Mediated Alternative Splicing Regulates Abiotic Stress Responses in Plants

The alternative splicing of pre-mRNA generates distinct mRNA variants from a pre-mRNA, thereby modulating a gene's function. The splicing of pre-mRNA depends on splice sites and regulatory elements in pre-mRNA, as well as the snRNA and proteins that recognize these sequences. Among these, RNA-binding proteins (RBPs) are the primary regulators of pre-mRNA splicing and play a critical role in the regulation of alternative splicing by recognizing the elements in pre-mRNA. However, little is known about the function of RBPs in stress response in plants. Here, we summarized the RBPs involved in the alternative splicing of pre-mRNA and their recognizing elements in pre-mRNA, and the recent advance in the role of RBP-mediated alternative splicing in response to abiotic stresses in plants. This review proposes that the regulation of pre-mRNA alternative splicing by RBPs is an important way for plants to adapt to abiotic stresses, and the regulation of alternative splicing by RBPs is a promising direction for crop breeding.

Full-length transcriptome sequencing of pepper fruit during development and construction of a transcript variation database

Chili pepper is an important spice and a model plant for fruit development studies. Large-scale omics information on chili pepper plant development continues to be gathered for understanding development as well as capsaicin biosynthesis. In this study, a full-spectrum transcriptome data of eight chili pepper tissues at five growth stages using the Oxford Nanopore long-read sequencing approach was generated. Of the 485 351 transcripts, 35 336 were recorded as reference transcripts (genes), while 450 015 were novel including coding, lnc, and other non-coding RNAs. These novel transcripts belonged to unknown/intergenic (347703), those retained introns (26336), and had multi-exons with at least one junction match (20333). In terms of alternative splicing, retained intron had the highest proportion (14795). The number of tissue-specific expressed transcripts ranged from 22 925 (stem) to 40 289 (flower). The expression changes during fruit and placenta development are discussed in detail. Integration of gene expression and capsaicin content quantification throughout the placental development clarifies that capsaicin biosynthesis in pepper is mainly derived from valine, leucin, and isoleucine degradation as well as citrate cycle and/or pyrimidine metabolism pathways. Most importantly, a user-friendly Pepper Full-Length Transcriptome Variation Database (PFTVD 1.0) (http://pepper-database.cn/) has been developed. PFTVD 1.0 provides transcriptomics and genomics information and allows users to analyse the data using various tools implemented. This work highlights the potential of long-read sequencing to discover novel genes and transcripts and their diversity in plant developmental biology.

The Functional Relationship Between RNA Splicing and the Chromatin Landscape

Chromatin is a highly regulated and dynamic structure that has been shown to play an essential role in transcriptional and co-transcriptional regulation. In the context of RNA splicing, early evidence suggested a loose connection between the chromatin landscape and splicing. More recently, it has been shown that splicing occurs in a co-transcriptional manner, meaning that the splicing process occurs in the context of chromatin. Experimental and computational evidence have also shown that chromatin dynamics can influence the splicing process and vice versa. However, much of this evidence provides mainly correlative relationships between chromatin and splicing with just a few concrete examples providing defined molecular mechanisms by which these two processes are functionally related. Nevertheless, it is clear that chromatin and RNA splicing are tightly interconnected to one another. In this review, we highlight the current state of knowledge of the relationship between chromatin and splicing.

SimSpliceEvol2: alternative splicing-aware simulation of biological sequence evolution and transcript phylogenies

BACKGROUND

SimSpliceEvol is a tool for simulating the evolution of eukaryotic gene sequences that integrates exon-intron structure evolution as well as the evolution of the sets of transcripts produced from genes. It takes a guide gene tree as input and generates a gene sequence with its transcripts for each node of the tree, from the root to the leaves. However, the sets of transcripts simulated at different nodes of the guide gene tree lack evolutionary connections. Consequently, SimSpliceEvol is not suitable for evaluating methods for transcript phylogeny inference or gene phylogeny inference that rely on transcript conservation.

RESULTS

Here, we introduce SimSpliceEvol2, which, compared to the first version, incorporates an explicit model of transcript evolution for simulating alternative transcripts along the branches of a guide gene tree, as well as the transcript phylogenies inferred. We offer a comprehensive software with a graphical user interface and an updated version of the web server, ensuring easy and user-friendly access to the tool.

CONCLUSION

SimSpliceEvol2 generates synthetic datasets that are useful for evaluating methods and tools for spliced RNA sequence analysis, such as spliced alignment methods, methods for identifying conserved transcripts, and transcript phylogeny reconstruction methods. The web server is accessible at https://simspliceevol.cobius.usherbrooke.ca , where you can also download the standalone software. Comprehensive documentation for the software is available at the same address. For developers interested in the source code, which requires the installation of all prerequisites to run, it is provided at  https://github.com/UdeS-CoBIUS/SimSpliceEvol .

Surveying the landscape of RNA isoform diversity and expression across 9 GTEx tissues using long-read sequencing data

Even though alternative RNA splicing was discovered nearly 50 years ago (1977), we still understand very little about most isoforms arising from a single gene, including in which tissues they are expressed and if their functions differ. Human gene annotations suggest remarkable transcriptional complexity, with approximately 252,798 distinct RNA isoform annotations from 62,710 gene bodies (Ensembl v109; 2023), emphasizing the need to understand their biological effects. For example, 256 gene bodies have ≥50 annotated isoforms and 30 have ≥100, where one protein-coding gene (MAPK10) even has 192 distinct RNA isoform annotations. Whether such isoform diversity results from biological redundancy or spurious alternative splicing (i.e., noise), or whether individual isoforms have specialized functions (even if subtle) remains a mystery for most genes. Recent studies by Aguzzoli-Heberle et al., Leung et al., and Glinos et al. demonstrated long-read RNAseq enables improved RNA isoform quantification for essentially any tissue, cell type, or biological condition (e.g., disease, development, aging, etc.), making it possible to better assess individual isoform expression and function. While each study provided important discoveries related to RNA isoform diversity, deeper exploration is needed. We sought to quantify and characterize real isoform usage across tissues (compared to annotations). We used long-read RNAseq data from 58 GTEx samples across nine tissues (three brain, two heart, muscle, lung, liver, and cultured fibroblasts) generated by Glinos et al. and found considerable isoform diversity within and across tissues. Cerebellar hemisphere was the most transcriptionally complex tissue (22,522 distinct isoforms; 3,726 unique); liver was least diverse (12,435 distinct isoforms; 1,039 unique). We highlight gene clusters exhibiting high tissue-specific isoform diversity per tissue (e.g., TPM1 expresses 19 in heart's atrial appendage). We also validated 447 of the 700 new isoforms discovered by Aguzzoli-Heberle et al. and found that 88 were expressed in all nine tissues, while 58 were specific to a single tissue. This study represents a broad survey of the RNA isoform landscape, demonstrating isoform diversity across nine tissues and emphasizes the need to better understand how individual isoforms from a single gene body contribute to human health and disease.

Unveiling a Microexon Switch: Novel Regulation of the Activities of Sugar Assimilation and Plant-Cell-Wall-Degrading Xylanases and Cellulases by Xlr2 in Trichoderma virens

Functional microexons have not previously been described in filamentous fungi. Here, we describe a novel mechanism of transcriptional regulation in Trichoderma requiring the inclusion of a microexon from the Xlr2 gene. In low-glucose environments, a long mRNA including the microexon encodes a protein with a GAL4-like DNA-binding domain (Xlr2-α), whereas in high-glucose environments, a short mRNA that is produced encodes a protein lacking this DNA-binding domain (Xlr2-β). Interestingly, the protein isoforms differ in their impact on cellulase and xylanase activity. Deleting the Xlr2 gene reduced both xylanase and cellulase activity and growth on different carbon sources, such as carboxymethylcellulose, xylan, glucose, and arabinose. The overexpression of either Xlr2-α or Xlr2-β in T. virens showed that the short isoform (Xlr2-β) caused higher xylanase activity than the wild types or the long isoform (Xlr2-α). Conversely, cellulase activity did not increase when overexpressing Xlr2-β but was increased with the overexpression of Xlr2-α. This is the first report of a novel transcriptional regulation mechanism of plant-cell-wall-degrading enzyme activity in T. virens. This involves the differential expression of a microexon from a gene encoding a transcriptional regulator.

Variant- and vaccination-specific alternative splicing profiles in SARS-CoV-2 infections

The COVID-19 pandemic, driven by the SARS-CoV-2 virus and its variants, highlights the important role of understanding host-viral molecular interactions influencing infection outcomes. Alternative splicing post-infection can impact both host responses and viral replication. We analyzed RNA splicing patterns in immune cells across various SARS-CoV-2 variants, considering immunization status. Using a dataset of 190 RNA-seq samples from our prior studies, we observed a substantial deactivation of alternative splicing and RNA splicing-related genes in COVID-19 patients. The alterations varied significantly depending on the infecting variant and immunization history. Notably, Alpha or Beta-infected patients differed from controls, while Omicron-infected patients displayed a splicing profile closer to controls. Particularly, vaccinated Omicron-infected individuals showed a distinct dynamic in alternative splicing patterns not widely shared among other groups. Our findings underscore the intricate interplay between SARS-CoV-2 variants, vaccination-induced immunity, and alternative splicing, emphasizing the need for further investigations to deepen understanding and guide therapeutic development.

Exitrons: offering new roles to retained introns-the novel regulators of protein diversity and utility

Exitrons are exonic introns. This subclass of intron retention alternative splicing does not contain a Pre-Terminating stop Codon. Therefore, when retained, they are always a part of a protein. Intron retention is a frequent phenomenon predominantly found in plants, which results in either the degradation of the transcripts or can serve as a stable intermediate to be processed upon induction by specific signals or the cell status. Interestingly, exitrons have coding ability and may confer additional attributes to the proteins that retain them. Therefore, exitron-containing and exitron-spliced isoforms will be a driving force for creating protein diversity in the proteome of an organism. This review establishes a basic understanding of exitron, discussing its genesis, key features, identification methods and functions. We also try to depict its other potential roles. The present review also aims to provide a fundamental background to those who found such exitronic sequences in their gene(s) and to speculate the future course of studies.

Bulk and single-cell alternative splicing analyses reveal roles of TRA2B in myogenic differentiation

Alternative splicing (AS) disruption has been linked to disorders of muscle development, as well as muscular atrophy. However, the precise changes in AS patterns that occur during myogenesis are not well understood. Here, we employed isoform long-reads RNA-seq (Iso-seq) and single-cell RNA-seq (scRNA-seq) to investigate the AS landscape during myogenesis. Our Iso-seq data identified 61,146 full-length isoforms representing 11,682 expressed genes, of which over 52% were novel. We identified 38,022 AS events, with most of these events altering coding sequences and exhibiting stage-specific splicing patterns. We identified AS dynamics in different types of muscle cells through scRNA-seq analysis, revealing genes essential for the contractile muscle system and cytoskeleton that undergo differential splicing across cell types. Gene-splicing analysis demonstrated that AS acts as a regulator, independent of changes in overall gene expression. Two isoforms of splicing factor TRA2B play distinct roles in myogenic differentiation by triggering AS of TGFBR2 to regulate canonical TGF-β signalling cascades differently. Our study provides a valuable transcriptome resource for myogenesis and reveals the complexity of AS and its regulation during myogenesis.

Age-Related Alternative Splicing: Driver or Passenger in the Aging Process?

Alternative splicing changes are closely linked to aging, though it remains unclear if they are drivers or effects. As organisms age, splicing patterns change, varying gene isoform levels and functions. These changes may contribute to aging alterations rather than just reflect declining RNA quality control. Three main splicing types-intron retention, cassette exons, and cryptic exons-play key roles in age-related complexity. These events modify protein domains and increase nonsense-mediated decay, shifting protein isoform levels and functions. This may potentially drive aging or serve as a biomarker. Fluctuations in splicing factor expression also occur with aging. Somatic mutations in splicing genes can also promote aging and age-related disease. The interplay between splicing and aging has major implications for aging biology, though differentiating correlation and causation remains challenging. Declaring a splicing factor or event as a driver requires comprehensive evaluation of the associated molecular and physiological changes. A greater understanding of how RNA splicing machinery and downstream targets are impacted by aging is essential to conclusively establish the role of splicing in driving aging, representing a promising area with key implications for understanding aging, developing novel therapeutical options, and ultimately leading to an increase in the healthy human lifespan.

The subgenome Saccharum spontaneum contributes to sugar accumulation in sugarcane as revealed by full-length transcriptomic analysis

INTRODUCTION

Modern sugarcane cultivars (Saccharum spp. hybrids) derived from crosses between S. officinarum and S. spontaneum, with high-sugar traits and excellent stress tolerance inherited respectively. However, the contribution of the S. spontaneum subgenome to sucrose accumulation is still unclear.

OBJECTIVE

To compensate for the absence of a high-quality reference genome, a transcriptome analysis method is needed to analyze the molecular basis of differential sucrose accumulation in sugarcane hybrids and to find clues to the contribution of the S. spontaneum subgenome to sucrose accumulation.

METHODS

PacBio full-length sequencing was used to complement genome annotation, followed by the identification of differential genes between the high and low sugar groups using differential alternative splicing analysis and differential expression analysis. At the subgenomic level, the factors responsible for differential sucrose accumulation were investigated from the perspective of transcriptional and post-transcriptional regulation.

RESULTS

A full-length transcriptome annotated at the subgenomic level was provided, complemented by 263,378 allele-defined transcript isoforms and 139,405 alternative splicing (AS) events. Differential alternative splicing (DA) analysis and differential expression (DE) analysis identified differential genes between high and low sugar groups and explained differential sucrose accumulation factors by the KEGG pathways. In some gene models, different or even opposite expression patterns of alleles from the same gene were observed, reflecting the potential evolution of these alleles toward novel functions in polyploid sugarcane. Among DA and DE genes in the sucrose source-sink complex pathway, we found some alleles encoding sucrose accumulation-related enzymes derived from the S. spontaneum subgenome were differentially expressed or had DA events between the two contrasting sugarcane hybrids.

CONCLUSION

Full-length transcriptomes annotated at the subgenomic level could better characterize sugarcane hybrids, and the S. spontaneum subgenome was found to contribute to sucrose accumulation.

Variant- and Vaccination-Specific Alternative Splicing Profiles in SARS-CoV-2 Infections

The COVID-19 pandemic, caused by the coronavirus SARS-CoV-2, and its subsequent variants has underscored the importance of understanding the host-viral molecular interactions to devise effective therapeutic strategies. A significant aspect of these interactions is the role of alternative splicing in modulating host responses and viral replication mechanisms. Our study sought to delineate the patterns of alternative splicing of RNAs from immune cells across different SARS-CoV-2 variants and vaccination statuses, utilizing a robust dataset of 190 RNA-seq samples from our previous studies, encompassing an average of 212 million reads per sample. We identified a dynamic alteration in alternative splicing and genes related to RNA splicing were highly deactivated in COVID-19 patients and showed variant- and vaccination-specific expression profiles. Overall, Omicron-infected patients exhibited a gene expression profile akin to healthy controls, unlike the Alpha or Beta variants. However, significantly, we found identified a subset of infected individuals, most pronounced in vaccinated patients infected with Omicron variant, that exhibited a specific dynamic in their alternative splicing patterns that was not widely shared amongst the other groups. Our findings underscore the complex interplay between SARS-CoV-2 variants, vaccination-induced immune responses, and alternative splicing, emphasizing the necessity for further investigations into these molecular cross-talks to foster deeper understanding and guide strategic therapeutic development.

A Case-Control Study of the Relationship Between Genetic Polymorphism and Cretinism in Xinjiang

Background

Cretinism is a subtype of congenital hypothyroidism, an endocrine disorder resulting from inadequate thyroid hormone production or receptor deficiency. Genetic abnormalities play a major role in the development of thyroid dysfunction.

Methods

We recruited 183 participants with cretinism and 119 healthy participants from the Xinjiang Uyghur Autonomous Region and randomly selected 29 tag single nucleotide polymorphisms (tSNPs) in TSHB, PAX8, TPO, NKX2-5, and TSHR in all participants. We compared genotype and allele frequencies between cases and controls utilizing the chi-squared test, logistic regression analysis, and haplotype analysis.

Results

Using the chi-squared test, a single SNP was found to be associated with cretinism (recessive model: rs3754363, OR = 0.46, 95% CI = 0.27-0.80, P = 0.00519; genotype model: P = 0.01677). We stratified neurological, myxedematous, and mixed type and determined that another SNP was associated with a higher risk when comparing myxedematous type to the neurological type (rs2277923).

Conclusion

rs3754363 has a statistically significant protective effect on people with cretinism, while rs2277923 may play a greater role in promoting the development of neurocretinism.

Recursive splicing discovery using lariats in total RNA sequencing

Recursive splicing is a non-canonical splicing mechanism in which an intron is removed in segments via multiple splicing reactions. Relatively few recursive splice sites have been identified with high confidence in human introns, and more comprehensive analyses are needed to better characterize where recursive splicing happens and whether or not it has a regulatory function. In this study, we use an unbiased approach using intron lariats to search for recursive splice sites in constitutive introns and alternative exons in the human transcriptome. We find evidence for recursive splicing in a broader range of intron sizes than previously reported and detail a new location for recursive splicing at the distal ends of cassette exons. In addition, we identify evidence for the conservation of these recursive splice sites among higher vertebrates and the use of these sites to influence alternative exon exclusion. Together, our data demonstrate the prevalence of recursive splicing and its potential influence on gene expression through alternatively spliced isoforms.

Design, execution, and interpretation of plant RNA-seq analyses

Genomics has transformed our understanding of the genetic architecture of traits and the genetic variation present in plants. Here, we present a review of how RNA-seq can be performed to tackle research challenges addressed by plant sciences. We discuss the importance of experimental design in RNA-seq, including considerations for sampling and replication, to avoid pitfalls and wasted resources. Approaches for processing RNA-seq data include quality control and counting features, and we describe common approaches and variations. Though differential gene expression analysis is the most common analysis of RNA-seq data, we review multiple methods for assessing gene expression, including detecting allele-specific gene expression and building co-expression networks. With the production of more RNA-seq data, strategies for integrating these data into genetic mapping pipelines is of increased interest. Finally, special considerations for RNA-seq analysis and interpretation in plants are needed, due to the high genome complexity common across plants. By incorporating informed decisions throughout an RNA-seq experiment, we can increase the knowledge gained.

ATTIC is an integrated approach for predicting A-to-I RNA editing sites in three species

A-to-I editing is the most prevalent RNA editing event, which refers to the change of adenosine (A) bases to inosine (I) bases in double-stranded RNAs. Several studies have revealed that A-to-I editing can regulate cellular processes and is associated with various human diseases. Therefore, accurate identification of A-to-I editing sites is crucial for understanding RNA-level (i.e. transcriptional) modifications and their potential roles in molecular functions. To date, various computational approaches for A-to-I editing site identification have been developed; however, their performance is still unsatisfactory and needs further improvement. In this study, we developed a novel stacked-ensemble learning model, ATTIC (A-To-I ediTing predICtor), to accurately identify A-to-I editing sites across three species, including Homo sapiens, Mus musculus and Drosophila melanogaster. We first comprehensively evaluated 37 RNA sequence-derived features combined with 14 popular machine learning algorithms. Then, we selected the optimal base models to build a series of stacked ensemble models. The final ATTIC framework was developed based on the optimal models improved by the feature selection strategy for specific species. Extensive cross-validation and independent tests illustrate that ATTIC outperforms state-of-the-art tools for predicting A-to-I editing sites. We also developed a web server for ATTIC, which is publicly available at http://web.unimelb-bioinfortools.cloud.edu.au/ATTIC/. We anticipate that ATTIC can be utilized as a useful tool to accelerate the identification of A-to-I RNA editing events and help characterize their roles in post-transcriptional regulation.

Decoding and Recoding of mRNA Sequences by the Ribosome

Faithful translation of messenger RNA (mRNA) into protein is essential to maintain protein homeostasis in the cell. Spontaneous translation errors are very rare due to stringent selection of cognate aminoacyl transfer RNAs (tRNAs) and the tight control of the mRNA reading frame by the ribosome. Recoding events, such as stop codon readthrough, frameshifting, and translational bypassing, reprogram the ribosome to make intentional mistakes and produce alternative proteins from the same mRNA. The hallmark of recoding is the change of ribosome dynamics. The signals for recoding are built into the mRNA, but their reading depends on the genetic makeup of the cell, resulting in cell-specific changes in expression programs. In this review, I discuss the mechanisms of canonical decoding and tRNA-mRNA translocation; describe alternative pathways leading to recoding; and identify the links among mRNA signals, ribosome dynamics, and recoding.

Recent expansion of metabolic versatility in Diplonema papillatum, the model species of a highly speciose group of marine eukaryotes

BACKGROUND

Diplonemid flagellates are among the most abundant and species-rich of known marine microeukaryotes, colonizing all habitats, depths, and geographic regions of the world ocean. However, little is known about their genomes, biology, and ecological role.

RESULTS

We present the first nuclear genome sequence from a diplonemid, the type species Diplonema papillatum. The ~ 280-Mb genome assembly contains about 32,000 protein-coding genes, likely co-transcribed in groups of up to 100. Gene clusters are separated by long repetitive regions that include numerous transposable elements, which also reside within introns. Analysis of gene-family evolution reveals that the last common diplonemid ancestor underwent considerable metabolic expansion. D. papillatum-specific gains of carbohydrate-degradation capability were apparently acquired via horizontal gene transfer. The predicted breakdown of polysaccharides including pectin and xylan is at odds with reports of peptides being the predominant carbon source of this organism. Secretome analysis together with feeding experiments suggest that D. papillatum is predatory, able to degrade cell walls of live microeukaryotes, macroalgae, and water plants, not only for protoplast feeding but also for metabolizing cell-wall carbohydrates as an energy source. The analysis of environmental barcode samples shows that D. papillatum is confined to temperate coastal waters, presumably acting in bioremediation of eutrophication.

CONCLUSIONS

Nuclear genome information will allow systematic functional and cell-biology studies in D. papillatum. It will also serve as a reference for the highly diverse diplonemids and provide a point of comparison for studying gene complement evolution in the sister group of Kinetoplastida, including human-pathogenic taxa.

Insights into established and emerging roles of SR protein family in plants and animals

Splicing of pre-mRNA is an essential part of eukaryotic gene expression. Serine-/arginine-rich (SR) proteins are highly conserved RNA-binding proteins present in all metazoans and plants. SR proteins are involved in constitutive and alternative splicing, thereby regulating the transcriptome and proteome diversity in the organism. In addition to their role in splicing, SR proteins are also involved in mRNA export, nonsense-mediated mRNA decay, mRNA stability, and translation. Due to their pivotal roles in mRNA metabolism, SR proteins play essential roles in normal growth and development. Hence, any misregulation of this set of proteins causes developmental defects in both plants and animals. SR proteins from the animal kingdom are extensively studied for their canonical and noncanonical functions. Compared with the animal kingdom, plant genomes harbor more SR protein-encoding genes and greater diversity of SR proteins, which are probably evolved for plant-specific functions. Evidence from both plants and animals confirms the essential role of SR proteins as regulators of gene expression influencing cellular processes, developmental stages, and disease conditions. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Processing > Splicing Regulation/Alternative Splicing.

Sex differences in gene expression and alternative splicing in the Chinese horseshoe bat

Sexually dimorphic traits are common in sexually reproducing organisms and can be encoded by differential gene regulation between males and females. Although alternative splicing is common mechanism in generating transcriptional diversity, its role in generating sex differences relative to differential gene expression is less clear. Here, we investigate the relative roles of differential gene expression and alternative splicing between male and female the horseshoe bat species, Rhinolophus sinicus. Horseshoe bats are an excellent model to study acoustic differences between sexes. Using RNA-seq analyses of two somatic tissues (brain and liver) from males and females of the same population, we identified 3,471 and 2,208 differentially expressed genes between the sexes (DEGs) in the brain and liver, respectively. DEGs were enriched with functional categories associated with physiological difference of the sexes (e.g.,gamete generation and energy production for reproduction in females). In addition, we also detected many differentially spliced genes between the sexes (DSGs, 2,231 and 1,027 in the brain and liver, respectively) which were mainly involved in regulation of RNA splicing and mRNA metabolic process. Interestingly, we found a significant enrichment of DEGs on the X chromosome, but not for DSGs. As for the extent of overlap between the two sets of genes, more than expected overlap of DEGs and DSGs was observed in the brain but not in the liver. This suggests that more complex tissues, such as the brain, may require the intricate and simultaneous interplay of both differential gene expression and splicing of genes to govern sex-specific functions. Overall, our results support that variation in gene expression and alternative splicing are important and complementary mechanisms governing sex differences.

A2TEA: Identifying trait-specific evolutionary adaptations

Background: Plants differ in their ability to cope with external stresses (e.g., drought tolerance). Genome duplications are an important mechanism to enable plant adaptation. This leads to characteristic footprints in the genome, such as protein family expansion. We explore genetic diversity and uncover evolutionary adaptation to stresses by exploiting genome comparisons between stress tolerant and sensitive species and RNA-Seq data sets from stress experiments. Expanded gene families that are stress-responsive based on differential expression analysis could hint at species or clade-specific adaptation, making these gene families exciting candidates for follow-up tolerance studies and crop improvement. Software: Integration of such cross-species omics data is a challenging task, requiring various steps of transformation and filtering. Ultimately, visualization is crucial for quality control and interpretation. To address this, we developed A2TEA: Automated Assessment of Trait-specific Evolutionary Adaptations, a Snakemake workflow for detecting adaptation footprints in silico. It functions as a one-stop processing pipeline, integrating protein family, phylogeny, expression, and protein function analyses. The pipeline is accompanied by an R Shiny web application that allows exploring, highlighting, and exporting the results interactively. This allows the user to formulate hypotheses regarding the genomic adaptations of one or a subset of the investigated species to a given stress. Conclusions: While our research focus is on crops, the pipeline is entirely independent of the underlying species and can be used with any set of species. We demonstrate pipeline efficiency on real-world datasets and discuss the implementation and limits of our analysis workflow as well as planned extensions to its current state. The A2TEA workflow and web application are publicly available at: https://github.com/tgstoecker/A2TEA.Workflow and https://github.com/tgstoecker/A2TEA.WebApp, respectively.

The implications of alternative pre-mRNA splicing in cell signal transduction

Cells produce multiple mRNAs through alternative splicing, which ensures proteome diversity. Because most human genes undergo alternative splicing, key components of signal transduction pathways are no exception. Cells regulate various signal transduction pathways, including those associated with cell proliferation, development, differentiation, migration, and apoptosis. Since proteins produced through alternative splicing can exhibit diverse biological functions, splicing regulatory mechanisms affect all signal transduction pathways. Studies have demonstrated that proteins generated by the selective combination of exons encoding important domains can enhance or attenuate signal transduction and can stably and precisely regulate various signal transduction pathways. However, aberrant splicing regulation via genetic mutation or abnormal expression of splicing factors negatively affects signal transduction pathways and is associated with the onset and progression of various diseases, including cancer. In this review, we describe the effects of alternative splicing regulation on major signal transduction pathways and highlight the significance of alternative splicing.

Branch point strength controls species-specific CAMK2B alternative splicing and regulates LTP

Regulation and functionality of species-specific alternative splicing has remained enigmatic to the present date. Calcium/calmodulin-dependent protein kinase IIβ (CaMKIIβ) is expressed in several splice variants and plays a key role in learning and memory. Here, we identify and characterize several primate-specific CAMK2B splice isoforms, which show altered kinetic properties and changes in substrate specificity. Furthermore, we demonstrate that primate-specific CAMK2B alternative splicing is achieved through branch point weakening during evolution. We show that reducing branch point and splice site strengths during evolution globally renders constitutive exons alternative, thus providing novel mechanistic insight into cis-directed species-specific alternative splicing regulation. Using CRISPR/Cas9, we introduce a weaker, human branch point sequence into the mouse genome, resulting in strongly altered Camk2b splicing in the brains of mutant mice. We observe a strong impairment of long-term potentiation in CA3-CA1 synapses of mutant mice, thus connecting branch point-controlled CAMK2B alternative splicing with a fundamental function in learning and memory.

Identification of sex-specific splicing via comparative transcriptome analysis of gonads from sea cucumber Apostichopus japonicus

Alternative splicing (AS) is an essential post-transcriptional regulation mechanism for sex differentiation and gonadal development, which has rarely been reported in marine invertebrates. Sea cucumber (Apostichopus japonicus) is an economically important marine benthic echinoderm with a potential XX/XY sex determination mechanism, whose molecular mechanism in the gonadal differentiation has not been clearly understood. In this study, we analyzed available RNA-seq datasets of male and female gonads to explore if AS mechanism exerts an essential function in sex differentiation and gonadal development of A. japonicus. In our results, a total of 20,338 AS events from 7219 alternatively spliced genes, and 189 sexually differential alternative splicing (DAS) events from 156 genes were identified in gonadal transcriptome of sea cucumber. Gene Ontology analysis indicated that these DAS genes were significantly enriched in spermatogenesis-related GO terms. Maximal Clique Centrality (MCC) was then applied for protein-protein interaction (PPI) analysis to search for protein interactions and hub DAS gene. Among all DAS genes, we identified 10 DAS genes closely related to spermatogenesis and (or) sperm motility and a hub gene dnah1. Thus, this study revealed that alternative isoforms were generated from certain genes in female and male gonads through alternative splicing, which may provide direct evidence that alternative splicing mechanisms participate in female and male gonads. These results suggested a novel perspective for explaining the molecular mechanisms underlying gonadal differentiation between male and female sea cucumbers.

Origins and Evolution of Human Tandem Duplicated Exon Substitution Events

The mutually exclusive splicing of tandem duplicated exons produces protein isoforms that are identical save for a homologous region that allows for the fine tuning of protein function. Tandem duplicated exon substitution events are rare, yet highly important alternative splicing events. Most events are ancient, their isoforms are highly expressed, and they have significantly more pathogenic mutations than other splice events. Here, we analyzed the physicochemical properties and functional roles of the homologous polypeptide regions produced by the 236 tandem duplicated exon substitutions annotated in the human gene set. We find that the most important structural and functional residues in these homologous regions are maintained, and that most changes are conservative rather than drastic. Three quarters of the isoforms produced from tandem duplicated exon substitution events are tissue-specific, particularly in nervous and cardiac tissues, and tandem duplicated exon substitution events are enriched in functional terms related to structures in the brain and skeletal muscle. We find considerable evidence for the convergent evolution of tandem duplicated exon substitution events in vertebrates, arthropods, and nematodes. Twelve human gene families have orthologues with tandem duplicated exon substitution events in both Drosophila melanogaster and Caenorhabditis elegans. Six of these gene families are ion transporters, suggesting that tandem exon duplication in genes that control the flow of ions into the cell has an adaptive benefit. The ancient origins, the strong indications of tissue-specific functions, and the evidence of convergent evolution suggest that these events may have played important roles in the evolution of animal tissues and organs.

Can tandem alternative splicing and evasion of premature termination codon surveillance contribute to attenuated Peutz-Jeghers syndrome?

Alternative use of short distance tandem sites such as NAGN_n AG are a common mechanism of alternative splicing; however, single nucleotide variants are rarely reported as likely to generate or to disrupt tandem splice sites. We identify a pathogenic intron 5 STK11 variant (NM_000455.4:c.[735-6A>G];[=]) segregating with the mucocutaneous features but not the hamartomatous polyps of Peutz-Jeghers syndrome in two individuals. By RNAseq analysis of peripheral blood mRNA, this variant was shown to generate a novel and preferentially used tandem proximal splice acceptor (AAGTGAAG). The variant transcript (NM_000455.4:c.734_734 + 1insTGAAG), which encodes a frameshift (p.[Tyr246Glufs*43]) constituted 36%-43% of STK11 transcripts suggesting partial escape from nonsense mediated mRNA decay and translation of a truncated protein. A review of the ClinVar database identified other similar variants. We suggest that nucleotide changes creating or disrupting tandem alternative splice sites are a pertinent disease mechanism and require contextualization for clinical reporting. Additionally, we hypothesize that some pathogenic STK11 variants cause an attenuated phenotype.

Single-molecule Real-time (SMRT) Sequencing Facilitates Transcriptome Research and Genome Annotation of the Fish Sillago sinica

As a newly described Sillaginidae species, Chinese sillago (Sillago sinica) needs a better understanding of gene annotation information. In this study, we reported the first full-length transcriptome data of S. sinica using the PacBio isoform sequencing Iso-seq and a description of transcriptome structure analysis. A total of 454,979 high-quality full-length transcripts were obtained by single-molecule real-time (SMRT) sequencing, which was corrected by Illumina sequencing data. After that, 66,948 non-redundant full-length transcripts were generated after mapping to the reference genome of S. sinica, including 49 fusion isoforms and 9,250 novel isoforms. 63,459 isoforms were successfully annotated by one of the Nr, Nt, SwissProt, Pfam, KOG, GO, and KEGG databases. Additionally, 30,987 alternative polyadenylation (APA) sites, 451,867 alternative splicing (AS) events, 21,928 long non-coding RNAs (lncRNAs) and 12,911 transcription factors (TFs) were identified. The full-length transcripts of S. sinica would provide a precious resource for characterizing the transcriptome of S. sinica and for the further study of gene function and regulatory mechanism of this species.

Integrated omics analysis identified genes and their splice variants involved in fruit development and metabolites production in Capsicum species

To date, several transcriptomic studies during fruit development have been reported; however, no comprehensive integrated study on expression diversity, alternative splicing, and metabolomic profiling was reported in Capsicum. This study analyzed RNA-seq data and untargeted metabolomic profiling from early green (EG), mature green (MG), and breaker (Br) fruit stages from two Capsicum species, i.e., C. annuum (Cann) and C. frutescens (Cfrut) from Northeast India. A total of 117,416 and 96,802 alternatively spliced events (AltSpli-events) were identified from Cann and Cfrut, respectively. Among AltSpli-events, intron retention (IR; 32.2% Cann and 25.75% Cfrut) followed by alternative acceptor (AA; 15.4% Cann and 18.9% Cfrut) were the most abundant in Capsicum. Around 7600 genes expressed in at least one fruit stage of Cann and Cfrut were AltSpli. The study identified spliced variants of genes including transcription factors (TFs) potentially involved in fruit development/ripening (Aux/IAA 16-like, ETR, SGR1, ARF, CaGLK2, ETR, CaAGL1, MADS-RIN, FUL1, SEPALLATA1), carotenoid (PDS, CA1, CCD4, NCED3, xanthoxin dehydrogenase, CaERF82, CabHLH100, CaMYB3R-1, SGR1, CaWRKY28, CaWRKY48, CaWRKY54), and capsaicinoids or flavonoid biosynthesis (CaMYB48, CaWRKY51), which were significantly differentially spliced (DS) between consecutive Capsicum fruit stages. Also, this study observed that differentially expressed isoforms (DEiso) from 38 genes with differentially spliced events (DSE) were significantly enriched in various metabolic pathways such as starch and sucrose metabolism, amino acid metabolism, cysteine cutin suberin and wax biosynthesis, and carotenoid biosynthesis. Furthermore, the metabolomic profiling revealed that metabolites from aforementioned pathways such as carbohydrates (mainly sugars such as D-fructose, D-galactose, maltose, and sucrose), organic acids (carboxylic acids), and peptide groups significantly altered during fruit development. Taken together, our findings could help in alternative splicing-based targeted studies of candidate genes involved in fruit development and ripening in Capsicum crop.

The nucleolus is the site for inflammatory RNA decay during infection

Inflammatory cytokines are key signaling molecules that can promote an immune response, thus their RNA turnover must be tightly controlled during infection. Most studies investigate the RNA decay pathways in the cytosol or nucleoplasm but never focused on the nucleolus. Although this organelle has well-studied roles in ribosome biogenesis and cellular stress sensing, the mechanism of RNA decay within the nucleolus is not completely understood. Here, we report that the nucleolus is an essential site of inflammatory pre-mRNA instability during infection. RNA-sequencing analysis reveals that not only do inflammatory genes have higher intronic read densities compared with non-inflammatory genes, but their pre-mRNAs are highly enriched in nucleoli during infection. Notably, nucleolin (NCL) acts as a guide factor for recruiting cytosine or uracil (C/U)-rich sequence-containing inflammatory pre-mRNAs and the Rrp6-exosome complex to the nucleolus through a physical interaction, thereby enabling targeted RNA delivery to Rrp6-exosomes and subsequent degradation. Consequently, Ncl depletion causes aberrant hyperinflammation, resulting in a severe lethality in response to LPS. Importantly, the dynamics of NCL post-translational modifications determine its functional activity in phases of LPS. This process represents a nucleolus-dependent pathway for maintaining inflammatory gene expression integrity and immunological homeostasis during infection.

The nucleolus is the traditional site for ribosomal RNA biogenesis. Here, the authors find that the nucleolus is a site of inflammatory pre-mRNA turnover and elucidated how immune homeostasis can be maintained by controlling inflammatory gene expression.

Full-Length Transcriptome Sequencing Reveals Alternative Splicing and lncRNA Regulation during Nodule Development in Glycine max

Alternative splicing (AS) is a ubiquitous phenomenon among eukaryotic intron-containing genes, which greatly contributes to transcriptome and proteome diversity. Here we performed the isoform sequencing (Iso-Seq) of soybean underground tissues inoculated and uninoculated with Rhizobium and obtained 200,681 full-length transcripts covering 26,183 gene loci. It was found that 80.78% of the multi-exon loci produced more than one splicing variant. Comprehensive analysis of these identified 7874 differentially splicing events with highly diverse splicing patterns during nodule development, especially in defense and transport-related processes. We further profiled genes with differential isoform usage and revealed that 2008 multi-isoform loci underwent stage-specific or simultaneous major isoform switches after Rhizobium inoculation, indicating that AS is a vital way to regulate nodule development. Moreover, we took the lead in identifying 1563 high-confidence long non-coding RNAs (lncRNAs) in soybean, and 157 of them are differentially expressed during nodule development. Therefore, our study uncovers the landscape of AS during the soybean-Rhizobium interaction and provides systematic transcriptomic data for future study of multiple novel directions in soybean.

Meta-Analysis Suggests That Intron Retention Can Affect Quantification of Transposable Elements from RNA-Seq Data

Simple Summary

Transposable elements (TEs) are repetitive sequences comprising more than one third of the human genome with the original ability to change their location within the genome. Owing to their repetitive nature, the quantification of TEs results often challenging. RNA-seq is a useful tool for genome-wide TEs quantification, nevertheless it also presents technical issues, including low reads mappability and erroneous quantification derived from the transcription of TEs fragments embedded in canonical transcripts. Fragments derived from TEs are found within the introns of most genes, which led to the hypothesis that intron retention (IR) can affect the unbiased quantification of TEs expression. Performing meta-analysis of public RNA-seq datasets, here we observe that IR can indeed impact the quantification of TEs by increasing the number of reads mapped on intronic TE copies. Our work highlights a correlation between IR and TEs expression measurement by RNA-seq that should be taken into account to achieve reliable TEs quantification, especially in samples characterized by extensive IR, because differential IR might be confused with differential TEs expression.

Abstract

Transposable elements (TEs), also known as “jumping genes”, are repetitive sequences with the capability of changing their location within the genome. They are key players in many different biological processes in health and disease. Therefore, a reliable quantification of their expression as transcriptional units is crucial to distinguish between their independent expression and the transcription of their sequences as part of canonical transcripts. TEs quantification faces difficulties of different types, the most important one being low reads mappability due to their repetitive nature preventing an unambiguous mapping of reads originating from their sequences. A large fraction of TEs fragments localizes within introns, which led to the hypothesis that intron retention (IR) can be an additional source of bias, potentially affecting accurate TEs quantification. IR occurs when introns, normally removed from the mature transcript by the splicing machinery, are maintained in mature transcripts. IR is a widespread mechanism affecting many different genes with cell type-specific patterns. We hypothesized that, in an RNA-seq experiment, reads derived from retained introns can introduce a bias in the detection of overlapping, independent TEs RNA expression. In this study we performed meta-analysis using public RNA-seq data from lymphoblastoid cell lines and show that IR can impact TEs quantification using established tools with default parameters. Reads mapped on intronic TEs were indeed associated to the expression of TEs and influence their correct quantification as independent transcriptional units. We confirmed these results using additional independent datasets, demonstrating that this bias does not appear in samples where IR is not present and that differential TEs expression does not impact on IR quantification. We concluded that IR causes the over-quantification of intronic TEs and differential IR might be confused with differential TEs expression. Our results should be taken into account for a correct quantification of TEs expression from RNA-seq data, especially in samples in which IR is abundant.

Intron retention coupled with nonsense-mediated decay is involved in cellulase biosynthesis in cellulolytic fungi

BACKGROUND

Knowledge on regulatory networks associated with cellulase biosynthesis is prerequisite for exploitation of such regulatory systems in enhancing cellulase production with low cost. The biological functions of intron retention (IR) and nonsense-mediated mRNA decay (NMD) in filamentous fungi is lack of study, let alone their roles in cellulase biosynthesis.

RESULTS

We found that major cellulase genes (cel7a, cel7b, and cel3a) exhibited concomitant decrease in IR rates and increase in their gene expression in T. reesei under cellulase-producing condition (cellulose and lactose) that was accompanied with a more active NMD pathway, as compared to cellulase non-producing condition (glucose). In the presence of the NMD pathway inhibitor that successfully repressed the NMD pathway, the mRNA levels of cellulase genes were sharply down-regulated, but the rates of IR in these genes were significantly up-regulated. Consistently, the cellulase activities were severely inhibited. In addition, the NMD pathway inhibitor caused the downregulated mRNA levels of two important genes of the target of rapamycin (TOR) pathway, trfkbp12 and trTOR1. The absence of gene trfkbp12 made the cellulase production in T. reesei more sensitive to the NMD pathway inhibitor.

CONCLUSIONS

All these findings suggest that the IR of cellulase genes regulates their own gene expression by coupling with the NMD pathway, which might involve the TOR pathway. Our results provide better understanding on intron retention, the NMD pathway, and cellulase production mechanism in filamentous fungi.

Identification and characterization of a new germline-specific marker vasa gene and its promoter in the giant freshwater prawn Macrobrachium rosenbergii

Vasa gene encodes a protein member of DEAD-box superfamily of ATP-dependent RNA helicases, which plays a key role in germline development in metazoans. In present study, we identified a new germline-specific marker Mrvasa in the prawn Macrobrachium rosenbergii, whose genomic DNA sequence consists of 14 exons and 13 introns. A 2516 bp of full-length Mrvasa cDNA encodes a protein of 603 amino acids. It contains nine conserved motifs, a zinc-finger motif, and RGG repeats. RT-PCR indicated that Mrvasa mRNA was specifically expressed in gonads. QPCR analysis further revealed that the expression of Mrvasa mRNA is much higher in testis than in ovary. In testis, the relative expression level of Mrvasa mRNA in late developing stage is significantly higher than that in early-middle developing stage. During ovarian development, no significant difference in expression was found. In situ hybridization demonstrated that Mrvasa mRNA was localized in germline cells including spermatogonia, spermatocytes, and spermatozoa in testes, and previtellogenic and vitellogenic oocytes in ovary. We then isolated the Mrvasa promoter and determined the transcription core region of this promoter. This is the first report on identification of vasa core promoter in crustaceans. Our results will provide a useful germline-specific marker Mrvasa for tracing germline cell formation and development in M. rosenbergii.

The alternative splicing landscape of a coral reef fish during a marine heatwave

Alternative splicing is a molecular mechanism that enables a single gene to encode multiple transcripts and proteins by post-transcriptional modification of pre-RNA molecules. Changes in the splicing scheme of genes can lead to modifications of the transcriptome and the proteome. This mechanism can enable organisms to respond to environmental fluctuations. In this study, we investigated patterns of alternative splicing in the liver of the coral reef fish Acanthochromis polyacanthus in response to the 2016 marine heatwave on the Great Barrier Reef. The differentially spliced (DS; n = 40) genes during the onset of the heatwave (i.e., 29.49°C or +1°C from average) were related to essential cellular functions such as the MAPK signaling system, Ca(2+) binding, and homeostasis. With the persistence of the heatwave for a period of one month (February to March), 21 DS genes were detected, suggesting that acute warming during the onset of the heatwave is more influential on alternative splicing than the continued exposure to elevated temperatures. After the heatwave, the water temperature cooled to ~24.96°C, and fish showed differential splicing of genes related to cyto-protection and post-damage recovery (n = 26). Two-thirds of the DS genes detected across the heatwave were also differentially expressed, revealing that the two molecular mechanisms act together in A. polyacanthus to cope with the acute thermal change. This study exemplifies how splicing patterns of a coral reef fish can be modified by marine heatwaves. Alternative splicing could therefore be a potential mechanism to adjust cellular physiological states under thermal stress and aid coral reef fishes in their response to more frequent acute thermal fluctuations in upcoming decades.

Cross-examination of engineered nanomaterials in crop production: Application and related implications

The review presents the current knowledge on the development and implementation of nanotechnology in crop production, giving particular attention to potential opportunities and challenges of the use of nano-sensors, nano-pesticides, and nano-fertilizers. Due to the size-dependent properties, e.g. high reactivity, targeted and controlled delivery of active ingredients, engineered nanomaterials (ENMs) are expected to be more efficient agrochemicals than conventional agents. Growing production and usage of ENMs result in the spread of ENMs in the environment. Because plants constitute an important component of the agri-ecosystem, they are subjected to the ENMs activity. A number of studies have confirmed the uptake and translocation of ENMs by plants as well as their positive/negative effects on plants. Here, these endpoints are briefly summarized to show the diversity of plant responses to ENMs. The review includes a detailed molecular analysis of ENMs-plant interactions. The transcriptomics, proteomics and metabolomics tools have been very recently employed to explore ENMs-induced effects in planta. The omics approach allows a comprehensive understanding of the specific machinery of ENMs occurring at the molecular level. The summary of data will be valuable in defining future studies on the ENMs-plant system, which is crucial for developing a suitable strategy for the ENMs usage.

Principles and correction of 5'-splice site selection

In Eukarya, immature mRNA transcripts (pre-mRNA) often contain coding sequences, or exons, interleaved by non-coding sequences, or introns. Introns are removed upon splicing, and further regulation of the retained exons leads to alternatively spliced mRNA. The splicing reaction requires the stepwise assembly of the spliceosome, a macromolecular machine composed of small nuclear ribonucleoproteins (snRNPs). This review focuses on the early stage of spliceosome assembly, when U1 snRNP defines each intron 5'-splice site (5'ss) in the pre-mRNA. We first introduce the splicing reaction and the impact of alternative splicing on gene expression regulation. Thereafter, we extensively discuss splicing descriptors that influence the 5'ss selection by U1 snRNP, such as sequence determinants, and interactions mediated by U1-specific proteins or U1 small nuclear RNA (U1 snRNA). We also include examples of diseases that affect the 5'ss selection by U1 snRNP, and discuss recent therapeutic advances that manipulate U1 snRNP 5'ss selectivity with antisense oligonucleotides and small-molecule splicing switches.

OTC intron 4 variations mediate pathogenic splicing patterns caused by the c.386G>A mutation in humans and spfash mice, and govern susceptibility to RNA-based therapies

Background

Aberrant splicing is a common outcome in the presence of exonic or intronic variants that might hamper the intricate network of interactions defining an exon in a specific gene context. Therefore, the evaluation of the functional, and potentially pathological, role of nucleotide changes remains one of the major challenges in the modern genomic era. This aspect has also to be taken into account during the pre-clinical evaluation of innovative therapeutic approaches in animal models of human diseases. This is of particular relevance when developing therapeutics acting on splicing, an intriguing and expanding research area for several disorders. Here, we addressed species-specific splicing mechanisms triggered by the OTC c.386G>A mutation, relatively frequent in humans, leading to Ornithine TransCarbamylase Deficiency (OTCD) in patients and spf^ash mice, and its differential susceptibility to RNA therapeutics based on engineered U1snRNA.

Methods

Creation and co-expression of engineered U1snRNAs with human and mouse minigenes, either wild-type or harbouring different nucleotide changes, in human (HepG2) and mouse (Hepa1-6) hepatoma cells followed by analysis of splicing pattern. RNA pulldown studies to evaluate binding of specific splicing factors.

Results

Comparative nucleotide analysis suggested a role for the intronic +10-11 nucleotides, and pull-down assays showed that they confer preferential binding to the TIA1 splicing factor in the mouse context, where TIA1 overexpression further increases correct splicing. Consistently, the splicing profile of the human minigene with mouse +10-11 nucleotides overlapped that of mouse minigene, and restored responsiveness to TIA1 overexpression and to compensatory U1snRNA. Swapping the human +10-11 nucleotides into the mouse context had opposite effects.

Moreover, the interplay between the authentic and the adjacent cryptic 5′ss in the human OTC dictates pathogenic mechanisms of several OTCD-causing 5′ss mutations, and only the c.386+5G>A change, abrogating the cryptic 5′ss, was rescuable by engineered U1snRNA.

Conclusions

Subtle intronic variations explain species-specific OTC splicing patterns driven by the c.386G>A mutation, and the responsiveness to engineered U1snRNAs, which suggests careful elucidation of molecular mechanisms before proposing translation of tailored therapeutics from animal models to humans.

Supplementary Information

The online version contains supplementary material available at 10.1186/s10020-021-00418-9.

Tissue-specific regulation of translational readthrough tunes functions of the traffic jam transcription factor

Translational readthrough (TR) occurs when the ribosome decodes a stop codon as a sense codon, resulting in two protein isoforms synthesized from the same mRNA. TR has been identified in several eukaryotic organisms; however, its biological significance and mechanism remain unclear. Here, we quantify TR of several candidate genes in Drosophila melanogaster and characterize the regulation of TR in the large Maf transcription factor Traffic jam (Tj). Using CRISPR/Cas9-generated mutant flies, we show that the TR-generated Tj isoform is expressed in a subset of neural cells of the central nervous system and is excluded from the somatic cells of gonads. Control of TR in Tj is critical for preservation of neuronal integrity and maintenance of reproductive health. The tissue-specific distribution of a release factor splice variant, eRF1H, plays a critical role in modulating differential TR of leaky stop codon contexts. Fine-tuning of gene regulatory functions of transcription factors by TR provides a potential mechanism for cell-specific regulation of gene expression.

Genome-wide Analysis of Alternative Gene Splicing Associated with Virulence in the Brown Planthopper Nilaparvata lugens (Hemiptera: Delphacidae)

Alternative splicing of protein coding genes plays a profound role in phenotypic variation for many eukaryotic organisms. The development of high-throughput sequencing and bioinformatics algorithms provides the possibility of genome-wide identification of alternative splicing events in eukaryotes. However, for the brown planthopper Nilaparvata lugens, a destructive pest of rice crops, whole-genome distribution of alternative splicing events and the role of alternative splicing in the phenotypic plasticity of virulence have not previously been estimated. Here, we developed an analysis pipeline to identify alternative splicing events in the genome of N. lugens. Differential expression analysis and functional annotation were performed on datasets related to different virulence phenotypes. In total, 27,880 alternative splicing events corresponding to 9,787 multi-exon genes were detected in N. lugens. Among them, specifically expressed alternative splicing transcripts in the virulent Mudgo population were enriched in metabolic process categories, while transcripts in the avirulent TN1 population were enriched in regulator activity categories. In addition, genes encoding odorant receptor, secreted saliva protein and xenobiotic metabolic P450 monooxygenase showed different splicing patterns between Mudgo population and TN1 population. Host change experiment also revealed that an isoform of a P450 gene could be specially induced by the stimulation of resistant rice variety Mudgo. This research pioneered a genome-wide study of alternative gene splicing in the rice brown planthopper. Differences in alternative splicing between virulent and avirulent populations indicated that alternative splicing might play an important role in the formation of virulence phenotypes in N. lugens.

Alternative splicing level related to intron size and organism complexity

Background

Alternative splicing is the process of selecting different combinations of splice sites to produce variably spliced mRNAs. However, the relationships between alternative splicing prevalence and level (ASP/L) and variations of intron size and organism complexity (OC) remain vague. Here, we developed a robust protocol to analyze the relationships between ASP/L and variations of intron size and OC. Approximately 8 Tb raw RNA-Seq data from 37 eumetazoan species were divided into three sets of species based on variations in intron size and OC.

Results

We found a strong positive correlation between ASP/L and OC, but no correlation between ASP/L and intron size across species. Surprisingly, ASP/L displayed a positive correlation with mean intron size of genes within individual genomes. Moreover, our results revealed that four ASP/L-related pathways contributed to the differences in ASP/L that were associated with OC. In particular, the spliceosome pathway displayed distinct genomic features, such as the highest gene expression level, conservation level, and fraction of disordered regions. Interestingly, lower or no obvious correlations were observed among these genomic features.

Conclusions

The positive correlation between ASP/L and OC ubiquitously exists in eukaryotes, and this correlation is not affected by the mean intron size of these species. ASP/L-related splicing factors may play an important role in the evolution of OC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08172-2.

Evolution of the Early Spliceosomal Complex-From Constitutive to Regulated Splicing

Pre-mRNA splicing is a major process in the regulated expression of genes in eukaryotes, and alternative splicing is used to generate different proteins from the same coding gene. Splicing is a catalytic process that removes introns and ligates exons to create the RNA sequence that codifies the final protein. While this is achieved in an autocatalytic process in ancestral group II introns in prokaryotes, the spliceosome has evolved during eukaryogenesis to assist in this process and to finally provide the opportunity for intron-specific splicing. In the early stage of splicing, the RNA 5' and 3' splice sites must be brought within proximity to correctly assemble the active spliceosome and perform the excision and ligation reactions. The assembly of this first complex, termed E-complex, is currently the least understood process. We focused in this review on the formation of the E-complex and compared its composition and function in three different organisms. We highlight the common ancestral mechanisms in S. cerevisiae, S. pombe, and mammals and conclude with a unifying model for intron definition in constitutive and regulated co-transcriptional splicing.

miRNA-Based Regulation of Alternative RNA Splicing in Metazoans

Alternative RNA splicing is an important regulatory process used by genes to increase their diversity. This process is mainly executed by specific classes of RNA binding proteins that act in a dosage-dependent manner to include or exclude selected exons in the final transcripts. While these processes are tightly regulated in cells and tissues, little is known on how the dosage of these factors is achieved and maintained. Several recent studies have suggested that alternative RNA splicing may be in part modulated by microRNAs (miRNAs), which are short, non-coding RNAs (~22 nt in length) that inhibit translation of specific mRNA transcripts. As evidenced in tissues and in diseases, such as cancer and neurological disorders, the dysregulation of miRNA pathways disrupts downstream alternative RNA splicing events by altering the dosage of splicing factors involved in RNA splicing. This attractive model suggests that miRNAs can not only influence the dosage of gene expression at the post-transcriptional level but also indirectly interfere in pre-mRNA splicing at the co-transcriptional level. The purpose of this review is to compile and analyze recent studies on miRNAs modulating alternative RNA splicing factors, and how these events contribute to transcript rearrangements in tissue development and disease.

Comparative Analysis of Alternative Splicing Profiles in Th Cell Subsets Reveals Extensive Cell Type-Specific Effects Modulated by a Network of Transcription Factors and RNA-Binding Proteins

Alternative splicing (AS) plays an important role in the development of many cell types; however, its contribution to Th subsets has been clearly defined. In this study, we compare mice naive CD4+ Th cells with Th1, Th2, Th17, and T regulatory cells and observed that the majority of AS events were retained intron, followed by skipped-exon events, with at least 1200 genes across cell types affected by AS events. A significant fraction of the AS events, especially retained intron events from the 72-h time point, were no longer observed 2 wk postdifferentiation, suggesting a role for AS in early activation and differentiation via preferential expression of specific isoforms required during T cell activation, but not for differentiation or effector function. Examining the protein consequence of the exon-skipping events revealed an abundance of structural proteins encoding for intrinsically unstructured peptide regions, followed by transmembrane helices, β strands, and polypeptide turn motifs. Analyses of expression profiles of RNA-binding proteins (RBPs) and their cognate binding sites flanking the discovered AS events revealed an enrichment for specific RBP recognition sites in each of the Th subsets. Integration with publicly available chromatin immunoprecipitation sequencing datasets for transcription factors support a model wherein lineage-determining transcription factors impact the RBP profile within the differentiating cells, and this differential expression contributes to AS of the transcriptome via a cascade of cell type-specific posttranscriptional rewiring events.

CircRNAs and their regulatory roles in cancers

Circular RNAs (circRNAs), a novel type of non-coding RNAs (ncRNAs), have a covalently closed circular structure resulting from pre-mRNA back splicing via spliceosome and ribozymes. They can be classified differently in accordance with different criteria. As circRNAs are abundant, conserved, and stable, they can be used as diagnostic markers in various diseases and targets to develop new therapies. There are various functions of circRNAs, including sponge for miR/proteins, role of scaffolds, templates for translation, and regulators of mRNA translation and stability. Without m7G cap and poly-A tail, circRNAs can still be degraded in several ways, including RNase L, Ago-dependent, and Ago-independent degradation. Increasing evidence indicates that circRNAs can be modified by N-6 methylation (m6A) in many aspects such as biogenesis, nuclear export, translation, and degradation. In addition, they have been proved to play a regulatory role in the progression of various cancers. Recently, methods of detecting circRNAs with high sensitivity and specificity have also been reported. This review presents a detailed overview of circRNAs regarding biogenesis, biomarker, functions, degradation, and dynamic modification as well as their regulatory roles in various cancers. It’s particularly summarized in detail in the biogenesis of circRNAs, regulation of circRNAs by m6A modification and mechanisms by which circRNAs affect tumor progression respectively. Moreover, existing circRNA detection methods and their characteristics are also mentioned.