Correcting for differential transcript coverage reveals a strong relationship between alternative splicing and organism complexity

The role of moonlighting proteins in neurogenesis

The complexity of the mammalian brain must arise from a comparably small number of genes. Proteins with moonlighting functions, i.e. entirely different functions in different compartments or cell types, contribute to multiply functional diversity. Here we review examples of such proteins with moonlighting functions during neurogenesis and in neuronal maturation. These range from cytoskeletal proteins acting as transcriptional regulators or synaptic proteins or exon junction proteins binding to and regulating the cytoskeleton to immediate early gene transcription factors regulating lipid metabolism in the endoplasmic reticulum. We further discuss how proteins with such moonlighting functions contribute to the heterogeneity of organelles shaping cell-type diversity in the brain.

Alternative splicing of Gllac7 regulate lignin degradation in Ganoderma lucidum

Laccases are key enzymes involved in lignin degradation in the white-rot fungus Ganoderma lucidum. Gllac7, which exists as two alternative splice isoforms, Gllac7.1 and Gllac7.2, plays a significant role in the lignin degradation process in G. lucidum. In this study, we examined the functions of the two isoforms and explored the transcriptional regulation mechanism of Gllac7.1. On lignin medium (LM), the expression of Gllac7.1 was 25.11-fold higher than on glucose medium (MM), whereas the expression of Gllac7.2 was 0.62-fold lower than that of MM. On MM, the OE/RNAi_Gllac7.1 and OE/RNAi_Gllac7.2 transformants exhibited no phenotypic differences from the wild-type (WT) strain. However, on LM, the OE_Gllac7.1 and RNAi_Gllac7.2 transformants exhibited accelerated growth and enhanced lignin degradation rates, reaching 1.28 to 2.04 times those of the WT. Moreover, RNAi_Gllac7.2 transformants formed primordia two to three days earlier than that of the WT. A C2H2 zinc finger protein, chr8g0158731, can bind to the promoter of Gllac7.1. Furthermore, RNAi_chr8g0158731 transformants displayed increased expression of Gllac7.1 and demonstrated improved bagasse degradation capabilities. These findings provide valuable insights into the roles of alternative splice isoforms of G. lucidum laccase, offering a foundation for future molecular marker-assisted breeding efforts.

Alternative splicing in a coral during heat stress acclimation and recovery

Climate change has caused drastic declines in corals. As sessile organisms, corals acclimate to environmental shifts through genome-wide changes in gene expression, epigenetic modifications, and alterations in microbiome composition. However, alternative splicing (AS), a conserved mechanism of stress response in many organisms, has been under-explored in corals. Using short-term acute thermal stress assays, we investigated patterns of AS in the scleractinian coral Acropora cervicornis during response to low (33°C), medium (35°C), and high (37°C) heat stress and subsequent overnight recovery. Our findings demonstrate reproducible dynamic shifts in AS of at least 40 percent of all genes during response to heat treatment and the recovery phase. The relative proportion of AS increased in response to heat stress and was primarily dominated by intron retention in specific classes of transcripts, including those related to splicing regulation itself. While AS returned to baseline levels post-exposure to low heat, AS persisted even after reprieve from higher levels of heat stress, which was associated with irreversible loss of photosynthetic efficiency of the symbiont. Our findings demonstrate that, although animals, corals are more plant-like in their likely usage of AS for regulating thermal stress response and recovery.

Organismal complexity strongly correlates with the number of protein families and domains

In the pregenomic era, scientists were puzzled by the observation that haploid genome size (the C-value) did not correlate well with organismal complexity. This phenomenon, called the "C-value paradox," is mostly explained by the fact that protein-coding genes occupy only a small fraction of eukaryotic genomes. When the first genome sequences became available, scientists were even more surprised by the fact that the number of genes (G-value) was also a poor predictor of complexity, which gave rise to the "G-value paradox." The proposed explanations usually invoke mechanisms that increase the information content of each individual gene (protein-protein interactions, intrinsic disorder, posttranslational modifications, alternative splicing, etc.). Less attention has been paid to mechanisms that increase the amount of genetic material but do not increase (or not to the same extent) the amount of information encoded in the genome, such as gene duplication and domain shuffling. Proteins belonging to the same family and/or sharing the same domains often carry out similar or even redundant functions. We thus hypothesized that an organism's number of different protein families and domains should be suitable predictors of organismal complexity. In agreement with our hypothesis, we observed that the number of protein families, clans, domains, and motifs increases from simple to progressively more complex organisms. In addition, these metrics correlate with the number of cell types better than and independently of the number of protein-coding genes and several previously proposed predictors of organismal complexity. Our observations have the potential to represent a resolution to the G-value paradox.

Molecular and Functional Convergences Associated with Complex Multicellularity in Eukarya

A key trait of Eukarya is the independent evolution of complex multicellularity in animals, land plants, fungi, brown algae, and red algae. This phenotype is characterized by the initial exaptation of cell-cell adhesion genes followed by the emergence of mechanisms for cell-cell communication, together with the expansion of transcription factor gene families responsible for cell and tissue identity. The number of cell types is commonly used as a quantitative proxy for biological complexity in comparative genomics studies. While expansions of individual gene families have been associated with variations in the number of cell types within individual complex multicellular lineages, the molecular and functional roles responsible for the independent evolution of complex multicellular across Eukarya remain poorly understood. We employed a phylogeny-aware strategy to conduct a genomic-scale search for associations between the number of cell types and the abundance of genomic components across a phylogenetically diverse set of 81 eukaryotic species, including species from all complex multicellular lineages. Our annotation schemas represent 2 complimentary aspects of genomic information: homology, represented by conserved sequences, and function, represented by Gene Ontology terms. We found many gene families sharing common biological themes that define complex multicellular to be independently expanded in 2 or more complex multicellular lineages, such as components of the extracellular matrix, cell-cell communication mechanisms, and developmental pathways. Additionally, we describe many previously unknown associations of biological themes and biological complexity, such as expansions of genes playing roles in wound response, immunity, cell migration, regulatory processes, and response to natural rhythms. Together, our findings unveil a set of functional and molecular convergences independently expanded in complex multicellular lineages likely due to the common selective pressures in their lifestyles.

Transcriptome analysis of two isolates of the tomato pathogen Cladosporium fulvum, uncovers genome-wide patterns of alternative splicing during a host infection cycle

Alternative splicing (AS) is a key element of eukaryotic gene expression that increases transcript and proteome diversity in cells, thereby altering their responses to external stimuli and stresses. While AS has been intensively researched in plants and animals, its frequency, conservation, and putative impact on virulence, are relatively still understudied in plant pathogenic fungi. Here, we profiled the AS events occurring in genes of Cladosporium fulvum isolates Race 5 and Race 4, during nearly a complete compatible infection cycle on their tomato host. Our studies revealed extensive heterogeneity in the transcript isoforms assembled from different isolates, infections, and infection timepoints, as over 80% of the transcript isoforms were singletons that were detected in only a single sample. Despite that, nearly 40% of the protein-coding genes in each isolate were predicted to be recurrently AS across the disparate infection timepoints, infections, and the two isolates. Of these, 37.5% were common to both isolates and 59% resulted in multiple protein isoforms, thereby putatively increasing proteome diversity in the pathogen by 31% during infections. An enrichment analysis showed that AS mostly affected genes likely to be involved in the transport of nutrients, regulation of gene expression, and monooxygenase activity, suggesting a role for AS in finetuning adaptation of C. fulvum on its tomato host during infections. Tracing the location of the AS genes on the fungal chromosomes showed that they were mostly located in repeat-rich regions of the core chromosomes, indicating a causal connection between gene location on the genome and propensity to AS. Finally, multiple cases of differential isoform usage in AS genes of C. fulvum were identified, suggesting that modulation of AS at different infection stages may be another way by which pathogens refine infections on their hosts.

Evolutionary genomics of the emergence of brown algae as key components of coastal ecosystems

Brown seaweeds are keystone species of coastal ecosystems, often forming extensive underwater forests, and are under considerable threat from climate change. In this study, analysis of multiple genomes has provided insights across the entire evolutionary history of this lineage, from initial emergence, through later diversification of the brown algal orders, down to microevolutionary events at the genus level. Emergence of the brown algal lineage was associated with a marked gain of new orthologous gene families, enhanced protein domain rearrangement, increased horizontal gene transfer events, and the acquisition of novel signaling molecules and key metabolic pathways, the latter notably related to biosynthesis of the alginate-based extracellular matrix, and halogen and phlorotannin biosynthesis. We show that brown algal genome diversification is tightly linked to phenotypic divergence, including changes in life cycle strategy and zoid flagellar structure. The study also showed that integration of large viral genomes has had a significant impact on brown algal genome content throughout the emergence of the lineage.

Post-transcriptional regulation supports the homeostatic expression of mature RNA

Gene expression regulation is a sophisticated, multi-stage process, and its robustness is critical to normal cell function and the survival of an organism. Previous studies indicate that differential gene expression at the RNA level is typically attenuated at the protein level through translational regulation. However, how post-transcriptional regulation (PTR) influences expression change during the RNA maturation process remains unclear. In this study, we investigated this by quantifying the magnitude of expression change in precursor RNA and mature RNA across a vast range of different biological conditions. We analyzed bulk tissue RNA sequencing data from 4689 samples, including healthy and diseased tissues from human, chimpanzee, rhesus macaque, and murine sources. We demonstrated that PTR tends to support homeostatic expression of mature RNA by amplifying normal tissue-specific expression of precursor RNA, while reducing expression change of precursor RNA in disease contexts. Our study provides insight into the general influence of PTR on gene expression homeostasis. Our analysis also suggests that intronic reads in RNA-seq studies may contain under-utilized information about disease associations. Additionally, our findings may assist in identifying new disease biomarkers and more effective ways of altering gene expression as a therapeutic strategy.

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.

A Procedure for Testing for Tokyo Type 1 Open-Ended Evolution

Tokyo Type 1 open-ended evolution (OEE) is a category of OEE that includes systems exhibiting the ongoing generation of adaptive novelty and ongoing growth in complexity. It can be considered as a necessary foundation for Tokyo Type 2 OEE (ongoing evolution of evolvability) and Tokyo Type 3 OEE (ongoing generation of major transitions). This article brings together five methods of analysis to form a procedure for testing for Tokyo Type 1 OEE. The procedure is presented as simply as possible, isolated from the complexities of any particular evolutionary system, and with a clear rationale for each step. In developing these steps, we also identify five key challenges in OEE. The last of these (achieving a higher order of complexity growth within a system exhibiting indefinitely scalable complexity) can be considered a grand challenge for Tokyo Type 1 OEE. Promising approaches to this grand challenge include also achieving one or both of Tokyo Types 2 and 3 OEE; this can be seen as one answer to why these other types of OEE are important, providing a unified view of OEE.

RNA structure in alternative splicing regulation: from mechanism to therapy

Alternative splicing is a highly intricate process that plays a crucial role in post-transcriptional regulation and significantly expands the functional proteome of a limited number of coding genes in eukaryotes. Its regulation is multifactorial, with RNA structure exerting a significant impact. Aberrant RNA conformations lead to dysregulation of splicing patterns, which directly affects the manifestation of disease symptoms. In this review, the molecular mechanisms of RNA secondary structure-mediated splicing regulation are summarized, with a focus on the complex interplay between aberrant RNA conformations and disease phenotypes resulted from splicing defects. This study also explores additional factors that reshape structural conformations, enriching our understanding of the mechanistic network underlying structure-mediated splicing regulation. In addition, an emphasis has been placed on the clinical role of targeting aberrant splicing corrections in human diseases. The principal mechanisms of action behind this phenomenon are described, followed by a discussion of prospective development strategies and pertinent challenges.

Characterization of alternative splicing during mammalian brain development reveals the extent of isoform diversity and potential effects on protein structural changes

Regulation of gene expression is critical for fate commitment of stem and progenitor cells during tissue formation. In the context of mammalian brain development, a plethora of studies have described how changes in the expression of individual genes characterize cell types across ontogeny and phylogeny. However, little attention has been paid to the fact that different transcripts can arise from any given gene through alternative splicing (AS). Considered a key mechanism expanding transcriptome diversity during evolution, assessing the full potential of AS on isoform diversity and protein function has been notoriously difficult. Here, we capitalize on the use of a validated reporter mouse line to isolate neural stem cells, neurogenic progenitors and neurons during corticogenesis and combine the use of short- and long-read sequencing to reconstruct the full transcriptome diversity characterizing neurogenic commitment. Extending available transcriptional profiles of the mammalian brain by nearly 50,000 new isoforms, we found that neurogenic commitment is characterized by a progressive increase in exon inclusion resulting in the profound remodeling of the transcriptional profile of specific cortical cell types. Most importantly, we computationally infer the biological significance of AS on protein structure by using AlphaFold2, revealing how radical protein conformational changes can arise from subtle changes in isoforms sequence. Together, our study reveals that AS has a greater potential to impact protein diversity and function than previously thought, independently from changes in gene expression.

GTDrift: a resource for exploring the interplay between genetic drift, genomic and transcriptomic characteristics in eukaryotes

We present GTDrift, a comprehensive data resource that enables explorations of genomic and transcriptomic characteristics alongside proxies of the intensity of genetic drift in individual species. This resource encompasses data for 1506 eukaryotic species, including 1413 animals and 93 green plants, and is organized in three components. The first two components contain approximations of the effective population size, which serve as indicators of the extent of random genetic drift within each species. In the first component, we meticulously investigated public databases to assemble data on life history traits such as longevity, adult body length and body mass for a set of 979 species. The second component includes estimations of the ratio between the rate of non-synonymous substitutions and the rate of synonymous substitutions (dN/dS) in protein-coding sequences for 1324 species. This ratio provides an estimate of the efficiency of natural selection in purging deleterious substitutions. Additionally, we present polymorphism-derived N e estimates for 66 species. The third component encompasses various genomic and transcriptomic characteristics. With this component, we aim to facilitate comparative transcriptomics analyses across species, by providing easy-to-use processed data for more than 16 000 RNA-seq samples across 491 species. These data include intron-centered alternative splicing frequencies, gene expression levels and sequencing depth statistics for each species, obtained with a homogeneous analysis protocol. To enable cross-species comparisons, we provide orthology predictions for conserved single-copy genes based on BUSCO gene sets. To illustrate the possible uses of this database, we identify the most frequently used introns for each gene and we assess how the sequencing depth available for each species affects our power to identify major and minor splice variants.

Heterobifunctional small molecules to modulate RNA function

RNA has diverse cellular functionality, including regulating gene expression, protein translation, and cellular response to stimuli, due to its intricate structures. Over the past decade, small molecules have been discovered that target functional structures within cellular RNAs and modulate their function. Simple binding, however, is often insufficient, resulting in low or even no biological activity. To overcome this challenge, heterobifunctional compounds have been developed that can covalently bind to the RNA target, alter RNA sequence, or induce its cleavage. Herein, we review the recent progress in the field of RNA-targeted heterobifunctional compounds using representative case studies. We identify critical gaps and limitations and propose a strategic pathway for future developments of RNA-targeted molecules with augmented functionalities.

Nucleotide-level distance metrics to quantify alternative splicing implemented in TranD

Advances in affordable transcriptome sequencing combined with better exon and gene prediction has motivated many to compare transcription across the tree of life. We develop a mathematical framework to calculate complexity and compare transcript models. Structural features, i.e. intron retention (IR), donor/acceptor site variation, alternative exon cassettes, alternative 5'/3' UTRs, are compared and the distance between transcript models is calculated with nucleotide level precision. All metrics are implemented in a PyPi package, TranD and output can be used to summarize splicing patterns for a transcriptome (1GTF) and between transcriptomes (2GTF). TranD output enables quantitative comparisons between: annotations augmented by empirical RNA-seq data and the original transcript models; transcript model prediction tools for longread RNA-seq (e.g. FLAIR versus Isoseq3); alternate annotations for a species (e.g. RefSeq vs Ensembl); and between closely related species. In C. elegans, Z. mays, D. melanogaster, D. simulans and H. sapiens, alternative exons were observed more frequently in combination with an alternative donor/acceptor than alone. Transcript models in RefSeq and Ensembl are linked and both have unique transcript models with empirical support. D. melanogaster and D. simulans, share many transcript models and long-read RNAseq data suggests that both species are under-annotated. We recommend combined references.

Random genetic drift sets an upper limit on mRNA splicing accuracy in metazoans

Most eukaryotic genes undergo alternative splicing (AS), but the overall functional significance of this process remains a controversial issue. It has been noticed that the complexity of organisms (assayed by the number of distinct cell types) correlates positively with their genome-wide AS rate. This has been interpreted as evidence that AS plays an important role in adaptive evolution by increasing the functional repertoires of genomes. However, this observation also fits with a totally opposite interpretation: given that 'complex' organisms tend to have small effective population sizes (Ne), they are expected to be more affected by genetic drift, and hence more prone to accumulate deleterious mutations that decrease splicing accuracy. Thus, according to this 'drift barrier' theory, the elevated AS rate in complex organisms might simply result from a higher splicing error rate. To test this hypothesis, we analyzed 3496 transcriptome sequencing samples to quantify AS in 53 metazoan species spanning a wide range of Ne values. Our results show a negative correlation between Ne proxies and the genome-wide AS rates among species, consistent with the drift barrier hypothesis. This pattern is dominated by low abundance isoforms, which represent the vast majority of the splice variant repertoire. We show that these low abundance isoforms are depleted in functional AS events, and most likely correspond to errors. Conversely, the AS rate of abundant isoforms, which are relatively enriched in functional AS events, tends to be lower in more complex species. All these observations are consistent with the hypothesis that variation in AS rates across metazoans reflects the limits set by drift on the capacity of selection to prevent gene expression errors.

Alternative Splicing Variation: Accessing and Exploiting in Crop Improvement Programs

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

Interplays between cis- and trans-Acting Factors for Alternative Splicing in Response to Environmental Changes during Biological Invasions of Ascidians

Alternative splicing (AS), a pivotal biological process contributing to phenotypic plasticity, creates a bridge linking genotypes with phenotypes. Despite its importance, the AS mechanisms underlying environmental response and adaptation have not been well studied, and more importantly, the cis- and trans-acting factors influencing AS variation remain unclear. Using the model invasive congeneric ascidians, Ciona robusta, and Ciona savignyi, we compared their AS responses to environmental changes and explored the potential determinants. Our findings unveiled swift and dynamic AS changes in response to environmental challenges, and differentially alternative spliced genes (DASGs) were functionally enriched in transmembrane transport processes. Interestingly, both the prevalence and level of AS in C. robusta were lower than those observed in C. savignyi. Furthermore, these two indices were higher under temperature stresses compared to salinity stresses in C. savignyi. All the observed patterns underscore the species-specific and environmental context-dependent AS responses to environmental challenges. The dissimilarities in genomic structure and exon/intron size distributions between these two species likely contributed to the observed AS variation. Moreover, we identified a total of 11 and 9 serine/arginine-rich splicing factors (SRSFs) with conserved domains and gene structures in the genomes of C. robusta and C. savignyi, respectively. Intriguingly, our analysis revealed that all detected SRSFs did not exhibit prevalent AS regulations. Instead, we observed AS control over a set of genes related to splicing factors and spliceosome components. Altogether, our results elucidate species-specific and environmental challenge-dependent AS response patterns in closely related invasive ascidians. The identified splicing factors and spliceosome components under AS control offer promising candidates for further investigations into AS-mediated rapid responses to environmental challenges complementary to SRSFs.

Inter-species association mapping links splice site evolution to METTL16 and SNRNP27K

Eukaryotic genes are interrupted by introns that are removed from transcribed RNAs by splicing. Patterns of splicing complexity differ between species, but it is unclear how these differences arise. We used inter-species association mapping with Saccharomycotina species to correlate splicing signal phenotypes with the presence or absence of splicing factors. Here, we show that variation in 5' splice site sequence preferences correlate with the presence of the U6 snRNA N6-methyladenosine methyltransferase METTL16 and the splicing factor SNRNP27K. The greatest variation in 5' splice site sequence occurred at the +4 position and involved a preference switch between adenosine and uridine. Loss of METTL16 and SNRNP27K orthologs, or a single SNRNP27K methionine residue, was associated with a preference for +4 U. These findings are consistent with splicing analyses of mutants defective in either METTL16 or SNRNP27K orthologs and models derived from spliceosome structures, demonstrating that inter-species association mapping is a powerful orthogonal approach to molecular studies. We identified variation between species in the occurrence of two major classes of 5' splice sites, defined by distinct interaction potentials with U5 and U6 snRNAs, that correlates with intron number. We conclude that variation in concerted processes of 5' splice site selection by U6 snRNA is associated with evolutionary changes in splicing signal phenotypes.

MicroRNA-Mediated Obstruction of Stem-loop Alternative Splicing (MIMOSAS): a global mechanism for the regulation of alternative splicing

While RNA secondary structures are critical to regulate alternative splicing of long-range pre-mRNA, the factors that modulate RNA structure and interfere with the recognition of the splice sites are largely unknown. Previously, we identified a small, non-coding microRNA that sufficiently affects stable stem structure formation of Nmnat pre-mRNA to regulate the outcomes of alternative splicing. However, the fundamental question remains whether such microRNA-mediated interference with RNA secondary structures is a global molecular mechanism for regulating mRNA splicing. We designed and refined a bioinformatic pipeline to predict candidate microRNAs that potentially interfere with pre-mRNA stem-loop structures, and experimentally verified splicing predictions of three different long-range pre-mRNAs in the Drosophila model system. Specifically, we observed that microRNAs can either disrupt or stabilize stem-loop structures to influence splicing outcomes. Our study suggests that MicroRNA-Mediated Obstruction of Stem-loop Alternative Splicing (MIMOSAS) is a novel regulatory mechanism for the transcriptome-wide regulation of alternative splicing, increases the repertoire of microRNA function and further indicates cellular complexity of post-transcriptional regulation.

HCS-Splice: A High-Content Screening Method to Advance the Discovery of RNA Splicing-Modulating Therapeutics

Nucleic acid therapeutics have demonstrated an impressive acceleration in recent years. They work through multiple mechanisms of action, including the downregulation of gene expression and the modulation of RNA splicing. While several drugs based on the former mechanism have been approved, few target the latter, despite the promise of RNA splicing modulation. To improve our ability to discover novel RNA splicing-modulating therapies, we developed HCS-Splice, a robust cell-based High-Content Screening (HCS) assay. By implementing the use of a two-colour (GFP/RFP) fluorescent splicing reporter plasmid, we developed a versatile, effective, rapid, and robust high-throughput strategy for the identification of potent splicing-modulating molecules. The HCS-Splice strategy can also be used to functionally confirm splicing mutations in human genetic disorders or to screen drug candidates. As a proof-of-concept, we introduced a dementia-related splice-switching mutation in the Microtubule-Associated Protein Tau (MAPT) exon 10 splicing reporter. We applied HCS-Splice to the wild-type and mutant reporters and measured the functional change in exon 10 inclusion. To demonstrate the applicability of the method in cell-based drug discovery, HCS-Splice was used to evaluate the efficacy of an exon 10-targeting siRNA, which was able to restore the correct alternative splicing balance.

MicroRNA-Mediated Obstruction of Stem-loop Alternative Splicing (MIMOSAS): a global mechanism for the regulation of alternative splicing

While RNA secondary structures are critical to regulate alternative splicing of long-range pre-mRNA, the factors that modulate RNA structure and interfere with the recognition of the splice sites are largely unknown. Previously, we identified a small, non-coding microRNA that sufficiently affects stable stem structure formation of Nmnat pre-mRNA to regulate the outcomes of alternative splicing. However, the fundamental question remains whether such microRNA-mediated interference with RNA secondary structures is a global molecular mechanism for regulating mRNA splicing. We designed and refined a bioinformatic pipeline to predict candidate microRNAs that potentially interfere with pre-mRNA stem-loop structures, and experimentally verified splicing predictions of three different long-range pre-mRNAs in the Drosophila model system. Specifically, we observed that microRNAs can either disrupt or stabilize stem-loop structures to influence splicing outcomes. Our study suggests that MicroRNA-Mediated Obstruction of Stem-loop Alternative Splicing (MIMOSAS) is a novel regulatory mechanism for the transcriptome-wide regulation of alternative splicing, increases the repertoire of microRNA function and further indicates cellular complexity of post-transcriptional regulation.

One-Sentence Summary

MicroRNA-Mediated Obstruction of Stem-loop Alternative Splicing (MIMOSAS) is a novel regulatory mechanism for the transcriptome-wide regulation of alternative splicing.

Alternative splicing of the SUMO1/2/3 transcripts affects cellular SUMOylation and produces functionally distinct SUMO protein isoforms

Substantial increases in the conjugation of the main human SUMO paralogs, SUMO1, SUMO2, and SUMO3, are observed upon exposure to different cellular stressors, and such increases are considered important to facilitate cell survival to stress. Despite their critical cellular role, little is known about how the levels of the SUMO modifiers are regulated in the cell, particularly as it relates to the changes observed upon stress. Here we characterize the contribution of alternative splicing towards regulating the expression of the main human SUMO paralogs under normalcy and three different stress conditions, heat-shock, cold-shock, and Influenza A Virus infection. Our data reveal that the normally spliced transcript variants are the predominant mature mRNAs produced from the SUMO genes and that the transcript coding for SUMO2 is by far the most abundant of all. We also provide evidence that alternatively spliced transcripts coding for protein isoforms of the prototypical SUMO proteins, which we refer to as the SUMO alphas, are also produced, and that their abundance and nuclear export are affected by stress in a stress- and cell-specific manner. Additionally, we provide evidence that the SUMO alphas are actively synthesized in the cell as their coding mRNAs are found associated with translating ribosomes. Finally, we provide evidence that the SUMO alphas are functionally different from their prototypical counterparts, with SUMO1α and SUMO2α being non-conjugatable to protein targets, SUMO3α being conjugatable but targeting a seemingly different subset of protein from those targeted by SUMO3, and all three SUMO alphas displaying different cellular distributions from those of the prototypical SUMOs. Thus, alternative splicing appears to be an important contributor to the regulation of the expression of the SUMO proteins and the cellular functions of the SUMOylation system.

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.

Gradistics: An underappreciated dimension in evolutionary space

The growth of complexity is an unsolved and underappreciated problem. We consider possible causes of this growth, hypotheses testing, molecular mechanisms, complexity measures, cases of simplification, and significance for biomedicine. We focus on a general ability of regulation, which is based on the growing information storage and processing capacities, as the main proxy of complexity. Natural selection is indifferent to complexity. However, complexification can be inferred from the same first principle, on which natural selection is founded. Natural selection depends on potentially unlimited reproduction under limited environmental conditions. Because of the demographic pressure, the simple ecological niches become fulfilled and diversified (due to species splitting and divergence). Diversification increases complexity of biocenoses. After the filling and diversification of simple niches, the more complex niches can arise. This is the 'atomic orbitals' (AO) model. Complexity has many shortcomings but it has an advantage. This advantage is ability to regulatory adaptation, including behavioral, formed in the evolution by means of genetic adaptation. Regulatory adaptation is much faster than genetic one because it is based on the information previously accumulated via genetic adaptation and learning. Regulatory adaptation further increases complexity of biocenoses. This is the 'regulatory advantage' (RA) model. The comparison of both models allows testable predictions. We focus on the animal evolution because of the appearance of higher regulatory level (nervous system), which is absent in other lineages, and relevance to humans (including biomedical aspects).

Exon Elongation Added Intrinsically Disordered Regions to the Encoded Proteins and Facilitated the Emergence of the Last Eukaryotic Common Ancestor

Most prokaryotic proteins consist of a single structural domain (SD) with little intrinsically disordered regions (IDRs) that by themselves do not adopt stable structures, whereas the typical eukaryotic protein comprises multiple SDs and IDRs. How eukaryotic proteins evolved to differ from prokaryotic proteins has not been fully elucidated. Here, we found that the longer the internal exons are, the more frequently they encode IDRs in eight eukaryotes including vertebrates, invertebrates, a fungus, and plants. Based on this observation, we propose the "small bang" model from the proteomic viewpoint: the protoeukaryotic genes had no introns and mostly encoded one SD each, but a majority of them were subsequently divided into multiple exons (step 1). Many exons unconstrained by SDs elongated to encode IDRs (step 2). The elongated exons encoding IDRs frequently facilitated the acquisition of multiple SDs to make the last common ancestor of eukaryotes (step 3). One prediction of the model is that long internal exons are mostly unconstrained exons. Analytical results of the eight eukaryotes are consistent with this prediction. In support of the model, we identified cases of internal exons that elongated after the rat-mouse divergence and discovered that the expanded sections are mostly in unconstrained exons and preferentially encode IDRs. The model also predicts that SDs followed by long internal exons tend to have other SDs downstream. This prediction was also verified in all the eukaryotic species analyzed. Our model accounts for the dichotomy between prokaryotic and eukaryotic proteins and proposes a selective advantage conferred by IDRs.

The Theory of Carcino-Evo-Devo and Its Non-Trivial Predictions

To explain the sources of additional cell masses in the evolution of multicellular organisms, the theory of carcino-evo-devo, or evolution by tumor neofunctionalization, has been developed. The important demand for a new theory in experimental science is the capability to formulate non-trivial predictions which can be experimentally confirmed. Several non-trivial predictions were formulated using carcino-evo-devo theory, four of which are discussed in the present paper: (1) The number of cellular oncogenes should correspond to the number of cell types in the organism. The evolution of oncogenes, tumor suppressor and differentiation gene classes should proceed concurrently. (2) Evolutionarily new and evolving genes should be specifically expressed in tumors (TSEEN genes). (3) Human orthologs of fish TSEEN genes should acquire progressive functions connected with new cell types, tissues and organs. (4) Selection of tumors for new functions in the organism is possible. Evolutionarily novel organs should recapitulate tumor features in their development. As shown in this paper, these predictions have been confirmed by the laboratory of the author. Thus, we have shown that carcino-evo-devo theory has predictive power, fulfilling a fundamental requirement for a new theory.

Where protein structure and cell diversity meet

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

MicroRNAs are deeply linked to the emergence of the complex octopus brain

Soft-bodied cephalopods such as octopuses are exceptionally intelligent invertebrates with a highly complex nervous system that evolved independently from vertebrates. Because of elevated RNA editing in their nervous tissues, we hypothesized that RNA regulation may play a major role in the cognitive success of this group. We thus profiled messenger RNAs and small RNAs in three cephalopod species including 18 tissues of the Octopus vulgaris. We show that the major RNA innovation of soft-bodied cephalopods is an expansion of the microRNA (miRNA) gene repertoire. These evolutionarily novel miRNAs were primarily expressed in adult neuronal tissues and during the development and had conserved and thus likely functional target sites. The only comparable miRNA expansions happened, notably, in vertebrates. Thus, we propose that miRNAs are intimately linked to the evolution of complex animal brains.

m6A modification of U6 snRNA modulates usage of two major classes of pre-mRNA 5' splice site

Alternative splicing of messenger RNAs is associated with the evolution of developmentally complex eukaryotes. Splicing is mediated by the spliceosome, and docking of the pre-mRNA 5' splice site into the spliceosome active site depends upon pairing with the conserved ACAGA sequence of U6 snRNA. In some species, including humans, the central adenosine of the ACAGA box is modified by N6 methylation, but the role of this m6A modification is poorly understood. Here, we show that m6A modified U6 snRNA determines the accuracy and efficiency of splicing. We reveal that the conserved methyltransferase, FIONA1, is required for Arabidopsis U6 snRNA m6A modification. Arabidopsis fio1 mutants show disrupted patterns of splicing that can be explained by the sequence composition of 5' splice sites and cooperative roles for U5 and U6 snRNA in splice site selection. U6 snRNA m6A influences 3' splice site usage. We generalise these findings to reveal two major classes of 5' splice site in diverse eukaryotes, which display anti-correlated interaction potential with U5 snRNA loop 1 and the U6 snRNA ACAGA box. We conclude that U6 snRNA m6A modification contributes to the selection of degenerate 5' splice sites crucial to alternative splicing.

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.

An Arabidopsis pre-RNA processing8a (prp8a) missense allele restores splicing of a subset of mis-spliced mRNAs

Eukaryotic precursor mRNAs often harbor noncoding introns that must be removed prior to translation. Accurate splicing of precursor messenger RNA depends on placement and assembly of small nuclear ribonucleoprotein (snRNP) sub-complexes of the spliceosome. Yeast (Saccharomyces cerevisiae) studies established a role in splice-site selection for PRE-RNA PROCESSING8 (PRP8), a conserved spliceosome scaffolding protein of the U5 snRNP. However, analogous splice-site selection studies in multicellular eukaryotes are lacking. Such studies are crucial for a comprehensive understanding of alternative splicing, which is extensive in plants and animals but limited in yeast. In this work, we describe an Arabidopsis (Arabidopsis thaliana) prp8a mutant that modulates splice-site selection. We isolated prp8a-14 from a screen for suppressors of pex14-6, which carries a splice-site mutation in the PEROXIN14 (PEX14) peroxisome biogenesis gene. To elucidate Arabidopsis PRP8A function in spliceosome fidelity, we combined prp8a-14 with various pex14 splice-site mutations and monitored the double mutants for physiological and molecular consequences of dysfunctional and functional peroxisomes that correspond to impaired and recovered splicing, respectively. prp8a-14 restored splicing and PEX14 function to alleles with mutations in the exonic guanine of the 5'-splice site but did not restore splicing or function to alleles with mutations in the intronic guanine of 5'- or 3'-splice sites. We used RNA-seq to reveal the systemic impact of prp8a-14 and found hundreds of differentially spliced transcripts and thousands of transcripts with significantly altered levels. Among differentially spliced transcripts, prp8a-14 significantly altered 5'- and 3'-splice-site utilization to favor sites resulting in shorter introns. This study provides a genetic platform for probing splicing in plants and hints at a role for plant PRP8 in splice-site selection.

Alternative splicing as a source of phenotypic diversity

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

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

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

LPS-pretreated MSC-conditioned medium optimized with 10-kDa filter attenuates the injury of H9c2 cardiomyocytes in a model of hypoxia/reoxygenation

Mesenchymal stem cell-derived conditioned medium (MSC-CM) improves cardiac function, which is partly attributed to released paracrine factors. Since such cardioprotection is moderate and transient, it's essential to optimize MSC-CM effective components to alleviate myocardial injury. To optimize MSC-CM, MSCs were treated with or without lipopolysaccharides (LPSs) for 48 h (serum-free), and the supernatant was collected. Then, LPS-CM (MSC stimulated by LPS) was further treated with LPS remover (LPS Re-CM) or was concentrated with a 10-kDa cutoff filter (10 kDa-CM). ELISA showed that all pretreatments increased levels of VEGF, HGF, and IGF except LPS remover; 10 kDa-CM was superior to other-CM. CCK-8 displayed that viability of injured H9c2 cells enhanced with the increase of MSC-CM concentration. We also found 10 kDa-CM significantly alleviated H9c2 hypoxia/reoxygenation (H/R) injury, as evidenced by increased Bcl-2/Bax ratio, decreased the levels of LDH and cTn. TEM, TUNEL, and H&E staining confirmed 10 kDa-CM inhibited H/R-induced H9c2 morphological changes. Proteomic analysis identified 41 differentially expressed proteins in 10 kDa-CM, among which anti-inflammation, pro-angiogenesis, and anti-apoptosis were related to cardiac protection. This study indicates that 10 kDa-CM protects H9c2 cardiomyocytes from H/R injury by preserving most of the protective factors, such as VEGF, HGF, and IGF, in MSC-CM.

In rice splice variants that restore the reading frame after frameshifting indel introduction are common, often induced by the indels and sometimes lead to organism-level rescue

The introduction of frameshifting non-3n indels enables the identification of gene-trait associations. However, it has been hypothesised that recovery of the original reading frame owing to usage of non-canonical splice forms could cause rescue. To date there is very little evidence for organism-level rescue by such a mechanism and it is unknown how commonly indels induce, or are otherwise associated with, frame-restoring splice forms. We perform CRISPR/Cas9 editing of randomly selected loci in rice to investigate these issues. We find that the majority of loci have a frame-restoring isoform. Importantly, three quarters of these isoforms are not seen in the absence of the indels, consistent with indels commonly inducing novel isoforms. This is supported by analysis in the context of NMD knockdowns. We consider in detail the two top rescue candidates, in wax deficient anther 1 (wda1) and brittle culm (bc10), finding that organismal-level rescue in both cases is strong but owing to different splice modification routes. More generally, however, as frame-restoring isoforms are low abundance and possibly too disruptive, such rescue we suggest to be the rare exception, not the rule. Nonetheless, assuming that indels commonly induce frame-restoring isoforms, these results emphasize the need to examine RNA level effects of non-3n indels and suggest that multiple non-3n indels in any given gene are advisable to probe a gene’s trait associations.

As protein coding genes are read in units of three (codons), insertions or deletions (indels) that are not a multiple of three long (non 3n) are expected to be especially harmful. Whether they are is important both for interpreting the results of non-3n indel experiments to probe a gene’s functional importance and for diagnostics. Particularly enigmatic are incidences where some non-3n changes in a gene compromise phenotypes while other seemingly comparable ones do not. One explanation for the latter is that a non-3n indel might be rescued via a frame-restoring splice form. Here we examine this hypothesis by inducing non-3n indels in many genes in rice and find that many non-3n indels are associated with a splice form that restores the reading frame. In the majority of these cases the indel appears to induce the potential rescuing splice form. We examine two top hit cases in detail and show functional rescue by splice modification. More generally, the frame-restoring forms are, however, low abundance and probably result in compromised proteins. We conclude then that splice mediated rescue is possible, but probably uncommon. Nonetheless it should not be overlooked in experimental design and interpretation.

Alternative splicing in seasonal plasticity and the potential for adaptation to environmental change

Seasonal plasticity is accomplished via tightly regulated developmental cascades that translate environmental cues into trait changes. Little is known about how alternative splicing and other posttranscriptional molecular mechanisms contribute to plasticity or how these mechanisms impact how plasticity evolves. Here, we use transcriptomic and genomic data from the butterfly Bicyclus anynana, a model system for seasonal plasticity, to compare the extent of differential expression and splicing and test how these axes of transcriptional plasticity differ in their potential for evolutionary change. Between seasonal morphs, we find that differential splicing affects a smaller but functionally unique set of genes compared to differential expression. Further, we find strong support for the novel hypothesis that spliced genes are more susceptible than differentially expressed genes to erosion of genetic variation due to selection on seasonal plasticity. Our results suggest that splicing plasticity is especially likely to experience genetic constraints that could affect the potential of wild populations to respond to rapidly changing environments.

Identification and Characterization of a New Splicing Variant of Prokineticin 2

Prokineticin 2 (PROK2) is a secreted bioactive peptide that regulates a variety of biological responses via two GPCRs, the prokineticin receptors (PROKRs). The aim of this study was to characterize a new alternatively spliced product of the prok2 gene consisting of four exons. The 40-amino acid peptide, designated PROK2C, is encoded by exon 1 and exon 4, and its expression was detected in the hippocampus and spinal cord of mice. PROK2C was expressed in a heterologous system, Pichia pastoris, and its binding specificity to the amino-terminal regions of PROKR1 and PROKR2 was investigated by GST pull-down experiments. In addition, the introduction of the unnatural amino acid p-benzoyl-L-phenylalanine using amber codon suppression technology demonstrated the role of tryptophan at position 212 of PROKR2 for PROK2C binding by photoactivatable cross-linking. The functional significance of this new isoform was determined in vivo by nociceptive experiments, which showed that PROK2C elicits strong sensitization of peripheral nociceptors to painful stimuli. In order to analyze the induction of PROK2C signal transduction, STAT3 and ERK phosphorylation levels were determined in mammalian CHO cells expressing PROKR1 and PROKR2. Our data show by in vivo and in vitro experiments that PROK2C can bind and activate both prokineticin receptors.

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.

The role of alternative splicing in adaptation and evolution

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

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.

Growth of Biological Complexity from Prokaryotes to Hominids Reflected in the Human Genome

The growth of complexity in evolution is a most intriguing phenomenon. Using gene phylostratigraphy, we showed this growth (as reflected in regulatory mechanisms) in the human genome, tracing the path from prokaryotes to hominids. Generally, the different regulatory gene families expanded at different times, yet only up to the Euteleostomi (bony vertebrates). The only exception was the expansion of transcription factors (TF) in placentals; however, we argue that this was not related to increase in general complexity. Surprisingly, although TF originated in the Prokaryota while chromatin appeared only in the Eukaryota, the expansion of epigenetic factors predated the expansion of TF. Signaling receptors, tumor suppressors, oncogenes, and aging- and disease-associated genes (indicating vulnerabilities in terms of complex organization and strongly enrichment in regulatory genes) also expanded only up to the Euteleostomi. The complexity-related gene properties (protein size, number of alternative splicing mRNA, length of untranslated mRNA, number of biological processes per gene, number of disordered regions in a protein, and density of TF–TF interactions) rose in multicellular organisms and declined after the Euteleostomi, and possibly earlier. At the same time, the speed of protein sequence evolution sharply increased in the genes that originated after the Euteleostomi. Thus, several lines of evidence indicate that molecular mechanisms of complexity growth were changing with time, and in the phyletic lineage leading to humans, the most salient shift occurred after the basic vertebrate body plan was fixed with bony skeleton. The obtained results can be useful for evolutionary medicine.

Differential transcript usage analysis of bulk and single-cell RNA-seq data with DTUrtle

Motivation

Each year, the number of published bulk and single-cell RNA-seq datasets is growing exponentially. Studies analyzing such data are commonly looking at gene-level differences, while the collected RNA-seq data inherently represents reads of transcript isoform sequences. Utilizing transcriptomic quantifiers, RNA-seq reads can be attributed to specific isoforms, allowing for analysis of transcript-level differences. A differential transcript usage (DTU) analysis is testing for proportional differences in a gene’s transcript composition, and has been of rising interest for many research questions, such as analysis of differential splicing or cell-type identification.

Results

We present the R package DTUrtle, the first DTU analysis workflow for both bulk and single-cell RNA-seq datasets, and the first package to conduct a ‘classical’ DTU analysis in a single-cell context. DTUrtle extends established statistical frameworks, offers various result aggregation and visualization options and a novel detection probability score for tagged-end data. It has been successfully applied to bulk and single-cell RNA-seq data of human and mouse, confirming and extending key results. In addition, we present novel potential DTU applications like the identification of cell-type specific transcript isoforms as biomarkers.

Availability and implementation

The R package DTUrtle is available at https://github.com/TobiTekath/DTUrtle with extensive vignettes and documentation at https://tobitekath.github.io/DTUrtle/.

Supplementary information

Supplementary data are available at Bioinformatics online.

A Participatory Approach for Balancing Accuracy and Complexity in Modeling Resilience and Robustness

Robustness and resilience are widely used in the biological sciences and related disciplines to describe how systems respond to change. Robustness is the ability to tolerate change without adapting or moving to another state. Resilience refers to the ability for a system to sustain a perturbation and maintain critical functions. Robustness and resilience transcend levels of biological organization, though they do not scale directly across levels. We live in an era of novel stressors and unprecedented change, including climate change, emerging environmental contaminants, and changes to earth's biogeochemical and hydrological cycles. We envision a common framework for developing models to predict the robustness and resilience of biological functions associated with complex systems that can transcend disciplinary boundaries. Conceptual and quantitative models of robustness and resilience must consider cross-scale interactions of potentially infinite complexity, but it is impossible to capture everything within a single model. Here, we discuss the need to balance accuracy and complexity when designing models, data collection, and downstream analyses to study robustness and resilience. We also consider the difficulties in defining the spatiotemporal domain when studying robustness and resilience as an emergent property of a complex system. We suggest a framework for implementing transdisciplinary research on robustness and resilience of biological systems that draws on participatory stakeholder engagement methods from the fields of conservation and natural resources management. Further, we suggest that a common, simplified model development framework for describing complex biological systems will provide new, broadly relevant educational tools. Efficient interdisciplinary collaboration to accurately develop a model of robustness and resilience would enable rapid, context-specific assessment of complex biological systems with benefits for a broad range of societally relevant problems.

Direct Nanopore Sequencing of mRNA Reveals Landscape of Transcript Isoforms in Apicomplexan Parasites

Alternative splicing is a widespread phenomenon in metazoans by which single genes are able to produce multiple isoforms of the gene product. However, this has been poorly characterized in apicomplexans, a major phylum of some of the most important global parasites. Efforts have been hampered by atypical transcriptomic features, such as the high AU content of Plasmodium RNA, but also the limitations of short-read sequencing in deciphering complex splicing events. In this study, we utilized the long read direct RNA sequencing platform developed by Oxford Nanopore Technologies to survey the alternative splicing landscape of Toxoplasma gondii and Plasmodium falciparum. We find that while native RNA sequencing has a reduced throughput, it allows us to obtain full-length or nearly full-length transcripts with comparable quantification to Illumina sequencing. By comparing these data with available gene models, we find widespread alternative splicing, particularly intron retention, in these parasites. Most of these transcripts contain premature stop codons, suggesting that in these parasites, alternative splicing represents a pathway to transcriptomic diversity, rather than expanding proteomic diversity. Moreover, alternative splicing rates are comparable between parasites, suggesting a shared splicing machinery, despite notable transcriptomic differences between the parasites. This study highlights a strategy in using long-read sequencing to understand splicing events at the whole-transcript level and has implications in the future interpretation of transcriptome sequencing studies.

IMPORTANCE We have used a novel nanopore sequencing technology to directly analyze parasite transcriptomes. The very long reads of this technology reveal the full-length genes of the parasites that cause malaria and toxoplasmosis. Gene transcripts must be processed in a process called splicing before they can be translated to protein. Our analysis reveals that these parasites very frequently only partially process their gene products, in a manner that departs dramatically from their human hosts.

Modulation of Yorkie activity by alternative splicing is required for developmental stability

The Hippo signaling pathway is a major regulator of organ growth, which controls the activity of the transcription coactivator Yorkie (Yki) in Drosophila and its homolog YAP in mammals. Both Yki and YAP proteins exist as alternatively spliced isoforms containing either one or two WW domains. The biological importance of this conserved alternative splicing event is unknown. Here, we identify the splicing factor B52 as a regulator of yki alternative splicing in Drosophila and show that B52 modulates growth in part through modulation of yki alternative splicing. Yki isoforms differ by their transcriptional activity as well as their ability to bind and bridge PPxY motifs-containing partners, and can compete in vivo. Strikingly, flies in which yki alternative splicing has been abrogated, thus expressing only Yki2 isoform, exhibit fluctuating wing asymmetry, a signal of developmental instability. Our results identify yki alternative splicing as a new level of modulation of the Hippo pathway, that is required for growth equilibration during development. This study provides the first demonstration that the process of alternative splicing contributes to developmental robustness.

MeDAS: a Metazoan Developmental Alternative Splicing database

Alternative splicing is widespread throughout eukaryotic genomes and greatly increases transcriptomic diversity. Many alternative isoforms have functional roles in developmental processes and are precisely temporally regulated. To facilitate the study of alternative splicing in a developmental context, we created MeDAS, a Metazoan Developmental Alternative Splicing database. MeDAS is an added-value resource that re-analyses publicly archived RNA-seq libraries to provide quantitative data on alternative splicing events as they vary across the time course of development. It has broad temporal and taxonomic scope and is intended to assist the user in identifying trends in alternative splicing throughout development. To create MeDAS, we re-analysed a curated set of 2232 Illumina polyA+ RNA-seq libraries that chart detailed time courses of embryonic and post-natal development across 18 species with a taxonomic range spanning the major metazoan lineages from Caenorhabditis elegans to human. MeDAS is freely available at https://das.chenlulab.com both as raw data tables and as an interactive browser allowing searches by species, tissue, or genomic feature (gene, transcript or exon ID and sequence). Results will provide details on alternative splicing events identified for the queried feature and can be visualised at the gene-, transcript- and exon-level as time courses of expression and inclusion levels, respectively.