The role of alternative splicing in adaptation and evolution

Identification of antennal alternative splicing by combining genome and full-length transcriptome analysis in Bactrocera dorsalis

Alternative splicing is an essential post-transcriptional regulatory mechanism that diversifies gene function by generating multiple protein isoforms from a single gene and act as a crucial role in insect environmental adaptation. Olfaction, a key sense for insect adaptation, relies heavily on the antennae, which are the primary olfactory organs expressing most of the olfactory genes. Despite the extensive annotation of olfactory genes within insect antennal tissues facilitated by high-throughput sequencing technology advancements, systematic analyses of alternative splicing are still relatively less. In this study, we focused on the oriental fruit fly (Bactrocera dorsalis), a significant pest of fruit crops. We performed a detailed analysis of alternative splicing in its antennae by utilizing the full-length transcriptome of its antennal tissue and the insect's genome. The results revealed 8600 non-redundant full-length transcripts identified in the oriental fruit fly antennal full-length transcriptome, spanning 4,145 gene loci. Over 40% of these loci exhibited multiple isoforms. Among these, 161 genes showed sex-biased isoform switching, involving seven different types of alternative splicing. Notably, events involving alternative transcription start sites (ATSS) and alternative transcription termination sites (ATTS) were the most common. Of all the genes undergoing ATSS and ATTS alternative splicing between male and female, 32 genes were alternatively spliced in protein coding regions, potentially affecting protein function. These genes were categorized based on the length of the sex-biased isoforms, with the highest difference in isoform fraction (dIF) associated with the ATSS type, including genes such as BdorABCA13, BdorCAT2, and BdorTSN3. Additionally, transcription factor binding sites for doublesex were identified upstream of both BdorABCA13 and BdorCAT2. Besides being expressed in the antennal tissues, BdorABCA13 and BdorCAT2 are also expressed in the mouthparts, legs, and genitalia of both female and male adults, suggesting their functional diversity. This study reveals alternative splicing events in the antennae of Bactrophora dorsalis from two aspects: odorant receptor genes and other types of genes expressed in the antennae. This study not only provides a research foundation for understanding the regulation of gene function by alternative splicing in the oriental fruit fly but also offers new insights for utilizing olfaction-based behavioral manipulation techniques to manage this pest.

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.

Fused in sarcoma (FUS) inhibits milk production efficiency in mammals

Efficient milk production in mammals confers evolutionary advantages by facilitating the transmission of energy from mother to offspring. However, the regulatory mechanism responsible for the gradual establishment of milk production efficiency in mammals, from marsupials to eutherians, remains elusive. Here, we find that mammary gland of the marsupial sugar glider contained milk components during adolescence, and that mammary gland development is less dynamically cyclic compared to that in placental mammals. Furthermore, fused in sarcoma (FUS) is found to be partially responsible for this establishment of low efficiency. In mouse model, FUS inhibit mammary epithelial cell differentiation through the cyclin-dependent kinase inhibitor p57Kip2, leading to lactation failure and pup starvation. Clinically, FUS levels are negatively correlated with milk production in lactating women. Overall, our results shed light on FUS as a negative regulator of milk production, providing a potential mechanism for the establishment of milk production from marsupial to eutherian mammals.

The landscape of abiotic and biotic stress-responsive splice variants with deep RNA-seq datasets in hot pepper

Alternative splicing (AS) is a widely observed phenomenon in eukaryotes that plays a critical role in development and stress responses. In plants, the large number of RNA-seq datasets in response to different environmental stressors can provide clues for identification of condition-specific and/or common AS variants for preferred agronomic traits. We report RNA-seq datasets (350.7 Gb) from Capsicum annuum inoculated with one of three bacteria, one virus, or one oomycete and obtained additional existing transcriptome datasets. In this study, we investigated the landscape of AS in response to environmental stressors, signaling molecules, and tissues from 425 total samples comprising 841.49 Gb. In addition, we identified genes that undergo AS under specific and shared stress conditions to obtain potential genes that may be involved in enhancing tolerance to stressors. We uncovered 1,642,007 AS events and identified 4,354 differential alternative splicing genes related to environmental stressors, tissues, and signaling molecules. This information and approach provide useful data for basic-research focused on enhancing tolerance to environmental stressors in hot pepper or establishing breeding programs.

Comprehensive deciphering the alternative splicing patterns involved in leaf morphogenesis of Liriodendron chinense

Alternative splicing (AS), a pivotal post-transcriptional regulatory mechanism, profoundly amplifies diversity and complexity of transcriptome and proteome. Liriodendron chinense (Hemsl.) Sarg., an excellent ornamental tree species renowned for its distinctive leaf shape, which resembles the mandarin jacket. Despite the documented potential genes related to leaf development of L. chinense, the underlying post-transcriptional regulatory mechanisms remain veiled. Here, we conducted a comprehensive analysis of the transcriptome to clarify the genome-wide landscape of the AS pattern and the spectrum of spliced isoforms during leaf developmental stages in L. chinense. Our investigation unveiled 50,259 AS events, involving 10,685 genes (32.9%), with intron retention as the most prevalent events. Notably, the initial stage of leaf development witnessed the detection of 804 differentially AS events affiliated with 548 genes. Although both differentially alternative splicing genes (DASGs) and differentially expressed genes (DEGs) were enriched into morphogenetic related pathways during the transition from fishhook (P2) to lobed (P7) leaves, there was only a modest degree of overlap between DASGs and DEGs. Furthermore, we conducted a comprehensively AS analysis on homologous genes involved in leaf morphogenesis, and most of which are subject to post-transcriptional regulation of AS. Among them, the AINTEGUMENTA-LIKE transcript factor LcAIL5 was characterization in detailed, which experiences skipping exon (SE), and two transcripts displayed disparate expression patterns across multiple stages. Overall, these findings yield a comprehensive understanding of leaf development regulation via AS, offering a novel perspective for further deciphering the mechanism of plant leaf morphogenesis.

Non-concordant epigenetic and transcriptional responses to acute thermal stress in western mosquitofish (Gambusia affinis)

Climate change is intensifying the frequency and severity of extreme temperatures. Understanding the molecular mechanisms underlying the ability to cope with acute thermal stress is key for predicting species' responses to extreme temperature events. While many studies have focused on the individual roles of gene expression, post-transcriptional processes and epigenetic modifications in response to acute thermal stress, the relative contribution of these molecular mechanisms remains unclear. The wide range of thermal limits of western mosquitofish (Gambusia affinis) provides an opportunity to explore this interplay. Here, we quantified changes in gene expression, alternative splicing, DNA methylation and microRNA (miRNA) expression in muscle tissue dissected from mosquitofish immediately after reaching high (CTmax) or low thermal limit (CTmin). Although the numbers of genes showing expression and splicing changes in response to acute temperature stress were small, we found a possibly larger and non-redundant role of splicing compared to gene expression, with more genes being differentially spliced (DSGs) than differentially expressed (DEGs), and little overlap between DSGs and DEGs. We also identified a small proportion of CpGs showing significant methylation change (i.e. differentially methylated cytosines, DMCs) in fish at thermal limits; however, there was no overlap between DEGs and genes annotated with DMCs in both CTmax and CTmin experiments. The weak interplay between epigenetic modifications and gene expression was further supported by our discoveries of no differentially expressed miRNAs. These findings provide novel insights into the relative role of different molecular mechanisms underlying immediate responses to extreme temperatures and demonstrate non-concordant responses of epigenetic and transcriptional mechanisms to acute temperature stress.

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.

Genetic constraints in genes exhibiting splicing plasticity in facultative diapause

Phenotypic plasticity is produced and maintained by processes regulating the transcriptome. While differential gene expression is among the most important of these processes, relatively little is known about other sources of transcriptional variation. Previous work suggests that alternative splicing plays an extensive and functionally unique role in transcriptional plasticity, though plastically spliced genes may be more constrained than the remainder of expressed genes. In this study, we explore the relationship between expression and splicing plasticity, along with the genetic diversity in those genes, in an ecologically consequential polyphenism: facultative diapause. Using 96 samples spread over two tissues and 10 timepoints, we compare the extent of differential splicing and expression between diapausing and direct developing pupae of the butterfly Pieris napi. Splicing differs strongly between diapausing and direct developing trajectories but alters a smaller and functionally unique set of genes compared to differential expression. We further test the hypothesis that among these expressed loci, plastically spliced genes are likely to experience the strongest purifying selection to maintain seasonally plastic phenotypes. Genes with unique transcriptional changes through diapause consistently had the lowest nucleotide diversity, and this effect was consistently stronger among genes that were differentially spliced compared to those with just differential expression through diapause. Further, the strength of negative selection was higher in the population expressing diapause every generation. Our results suggest that maintenance of the molecular mechanisms involved in diapause progression, including post-transcriptional modifications, are highly conserved and likely to experience genetic constraints, especially in northern populations of P. napi.

Gene expression and alternative splicing contribute to adaptive divergence of ecotypes

Regulation of gene expression is a critical link between genotype and phenotype explaining substantial heritable variation within species. However, we are only beginning to understand the ways that specific gene regulatory mechanisms contribute to adaptive divergence of populations. In plants, the post-transcriptional regulatory mechanism of alternative splicing (AS) plays an important role in both development and abiotic stress response, making it a compelling potential target of natural selection. AS allows organisms to generate multiple different transcripts/proteins from a single gene and thus may provide a source of evolutionary novelty. Here, we examine whether variation in alternative splicing and gene expression levels might contribute to adaptation and incipient speciation of dune-adapted prairie sunflowers in Great Sand Dunes National Park, Colorado, USA. We conducted a common garden experiment to assess transcriptomic variation among ecotypes and analyzed differential expression, differential splicing, and gene coexpression. We show that individual genes are strongly differentiated for both transcript level and alternative isoform proportions, even when grown in a common environment, and that gene coexpression networks are disrupted between ecotypes. Furthermore, we examined how genome-wide patterns of sequence divergence correspond to divergence in transcript levels and isoform proportions and find evidence for both cis and trans-regulation. Together, our results emphasize that alternative splicing has been an underappreciated mechanism providing source material for natural selection at short evolutionary time scales.

Differential transcript usage across mammalian oocytes at the germinal vesicle and metaphase II stages

Ongoing progress in mRNA-Sequencing technologies has significantly contributed to the refinement of assisted reproductive technologies. However, the prior investigations have predominantly concentrated on alterations in overall gene expression levels, thereby leaving a considerable gap in our understanding of the influence of transcript isoform expression on fundamental cellular mechanisms of oocytes. Given the efficacy of differential transcript usage (DTU) analysis to address such knowledge, we conducted comprehensive DTU analysis utilizing mRNA-Seq datasets of germinal vesicle (GV) and metaphase II (MII) oocytes across six mammalian species from the SRA database, including cow, donkey, horse, human, mouse, and pig. To further illuminate the roles of these genes, we also conducted a rigorous Gene Ontology (GO) term enrichment analysis. While the DTU analysis of each species exhibited several genes with alterations in their transcript isoform usage, referred to as DTU genes, this study focused on only ten cross-species DTU genes sharing among a minimum of five distinct species (FDR≤0.05). These cross-species DTU genes were as follows: ABCF1, CDC6, CFAP36, CNOT10, DNM3, IWS1, NBN, NDEL1, RAD50 and ZCCHC17. GO term enrichment analysis unveiled the alignment of these cross-species DTU gene functions with RNA and cell-cycle control mechanisms across diverse mammalian species, thereby suggesting their vital roles during oocyte maturation. Further exploration of the transcript isoforms of these genes hence bore the potential to uncover novel transcript isoform markers for future reproductive technologies in both human and animal contexts.

SpliceProt 2.0: A Sequence Repository of Human, Mouse, and Rat Proteoforms

SpliceProt 2.0 is a public proteogenomics database that aims to list the sequence of known proteins and potential new proteoforms in human, mouse, and rat proteomes. This updated repository provides an even broader range of computationally translated proteins and serves, for example, to aid with proteomic validation of splice variants absent from the reference UniProtKB/SwissProt database. We demonstrate the value of SpliceProt 2.0 to predict orthologous proteins between humans and murines based on transcript reconstruction, sequence annotation and detection at the transcriptome and proteome levels. In this release, the annotation data used in the reconstruction of transcripts based on the methodology of ternary matrices were acquired from new databases such as Ensembl, UniProt, and APPRIS. Another innovation implemented in the pipeline is the exclusion of transcripts predicted to be susceptible to degradation through the NMD pathway. Taken together, our repository and its applications represent a valuable resource for the proteogenomics community.

PTBP1 as a potential regulator of disease

Polypyrimidine tract-binding protein 1 (PTBP1) is a member of the heterogeneous nuclear ribonucleoprotein (hnRNP) family, which plays a key role in alternative splicing of precursor mRNA and RNA metabolism. PTBP1 is universally expressed in various tissues and binds to multiple downstream transcripts to interfere with physiological and pathological processes such as the tumor growth, body metabolism, cardiovascular homeostasis, and central nervous system damage, showing great prospects in many fields. The function of PTBP1 involves the regulation and interaction of various upstream molecules, including circular RNAs (circRNAs), microRNAs (miRNAs) and long non-coding RNAs (lncRNAs). These regulatory systems are inseparable from the development and treatment of diseases. Here, we review the latest knowledge regarding the structure and molecular functions of PTBP1 and summarize its functions and mechanisms of PTBP1 in various diseases, including controversial studies. Furthermore, we recommend future studies on PTBP1 and discuss the prospects of targeting PTBP1 in new clinical therapeutic approaches.

Identification of cell-type specific alternative transcripts in the multicellular alga Volvox carteri

BACKGROUND

Cell type specialization is a hallmark of complex multicellular organisms and is usually established through implementation of cell-type-specific gene expression programs. The multicellular green alga Volvox carteri has just two cell types, germ and soma, that have previously been shown to have very different transcriptome compositions which match their specialized roles. Here we interrogated another potential mechanism for differentiation in V. carteri, cell type specific alternative transcript isoforms (CTSAI).

METHODS

We used pre-existing predictions of alternative transcripts and de novo transcript assembly with HISAT2 and Ballgown software to compile a list of loci with two or more transcript isoforms, identified a small subset that were candidates for CTSAI, and manually curated this subset of genes to remove false positives. We experimentally verified three candidates using semi-quantitative RT-PCR to assess relative isoform abundance in each cell type.

RESULTS

Of the 1978 loci with two or more predicted transcript isoforms 67 of these also showed cell type isoform expression biases. After curation 15 strong candidates for CTSAI were identified, three of which were experimentally verified, and their predicted gene product functions were evaluated in light of potential cell type specific roles. A comparison of genes with predicted alternative splicing from Chlamydomonas reinhardtii, a unicellular relative of V. carteri, identified little overlap between ortholog pairs with alternative splicing in both species. Finally, we interrogated cell type expression patterns of 126 V. carteri predicted RNA binding protein (RBP) encoding genes and found 40 that showed either somatic or germ cell expression bias. These RBPs are potential mediators of CTSAI in V. carteri and suggest possible pre-adaptation for cell type specific RNA processing and a potential path for generating CTSAI in the early ancestors of metazoans and plants.

CONCLUSIONS

We predicted numerous instances of alternative transcript isoforms in Volvox, only a small subset of which showed cell type specific isoform expression bias. However, the validated examples of CTSAI supported existing hypotheses about cell type specialization in V. carteri, and also suggested new hypotheses about mechanisms of functional specialization for their gene products. Our data imply that CTSAI operates as a minor but important component of V. carteri cellular differentiation and could be used as a model for how alternative isoforms emerge and co-evolve with cell type specialization.

Stochastic modeling of the mRNA life process: A generalized master equation

The mRNA life cycle is a complex biochemical process, involving transcription initiation, elongation, termination, splicing, and degradation. Each of these molecular events is multistep and can create a memory. The effect of this molecular memory on gene expression is not clear, although there are many related yet scattered experimental reports. To address this important issue, we develop a general theoretical framework formulated as a master equation in the sense of queue theory, which can reduce to multiple previously studied gene models in limiting cases. This framework allows us to interpret experimental observations, extract kinetic parameters from experimental data, and identify how the mRNA kinetics vary under regulatory influences. Notably, it allows us to evaluate the influences of elongation processes on mature RNA distribution; e.g., we find that the non-exponential elongation time can induce the bimodal mRNA expression and there is an optimal elongation noise intensity such that the mature RNA noise achieves the lowest level. In a word, our framework can not only provide insight into complex mRNA life processes but also bridge a dialogue between theoretical studies and experimental data.

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.

Genome sequencing-based transcriptomic analysis reveals novel genes in Peucedanum praeruptorum

BACKGROUND

Peucedanum praeruptorum Dunn, a traditional Chinese herbal medicine, contains coumarin and volatile oil components that have clinical application value. However, early bolting often occurs in the medicinal materials of Apiaceae plants. The rhizomes of the medicinal parts are gradually lignified after bolting, resulting in a sharp decrease in the content of coumarins. At present, the link between coumarin biosynthesis and early bolting in P. praeruptorum has not been elucidated.

RESULTS

Combining the genome sequencing and the previous transcriptome sequencing results, we reanalyzed the differential transcripts of P. praeruptorum before and after bolting. A total of 62,088 new transcripts were identified, of which 31,500 were unknown transcripts. Functional classification and annotation showed that many genes were involved in the regulation of transcription, defense response, and carbohydrate metabolic processes. The main domains are the pentatricopeptide repeat, protein kinase, RNA recognition motif, leucine-rich repeat, and ankyrin repeat domains, indicating their pivotal roles in protein modification and signal transduction. Gene structure analysis showed that skipped exon (SE) was the most dominant alternative splicing, followed by the alternative 3' splice site (A3SS) and the alternative 5' splice site (A5SS). Functional enrichment of differentially expressed genes showed that these differentially expressed genes mainly include transmembrane transporters, channel proteins, DNA-binding proteins, polysaccharide-binding proteins, etc. In addition, genes involved in peroxisome, hexose phosphate pathway, phosphatidylinositol signaling system, and inositol phosphate metabolism pathway were greatly enriched. A protein-protein interaction network analysis discoverd 1,457 pairs of proteins that interact with each other. The expression levels of six UbiA genes, three UGT genes, and four OMT genes were higher during the bolting stage. This observation suggests their potential involvement in the catalytic processes of prenylation, glycosylation, and methylation of coumarins, respectively. A total of 100 peroxidase (PRX) genes were identified being involved in lignin polymerization, but only nine PRX genes were highly expressed at the bolting stage. It is worth noting that 73 autophagy-related genes (ATGs) were first identified from the KEGG pathway-enriched genes. Some ATGs, such as BHQH00009837, BHQH00013830, and novel8944, had higher expression levels after bolting.

CONCLUSIONS

Comparative transcriptome analysis and large-scale genome screening provide guidance and new opinions for the identification of bolting-related genes in P. praeruptorum.

Bioinformatics Analysis of MSH1 Genes of Green Plants: Multiple Parallel Length Expansions, Intron Gains and Losses, Partial Gene Duplications, and Alternative Splicing

MutS homolog 1 (MSH1) is involved in the recombining and repairing of organelle genomes and is essential for maintaining their stability. Previous studies indicated that the length of the gene varied greatly among species and detected species-specific partial gene duplications in Physcomitrella patens. However, there are critical gaps in the understanding of the gene size expansion, and the extent of the partial gene duplication of MSH1 remains unclear. Here, we screened MSH1 genes in 85 selected species with genome sequences representing the main clades of green plants (Viridiplantae). We identified the MSH1 gene in all lineages of green plants, except for nine incomplete species, for bioinformatics analysis. The gene is a singleton gene in most of the selected species with conserved amino acids and protein domains. Gene length varies greatly among the species, ranging from 3234 bp in Ostreococcus tauri to 805,861 bp in Cycas panzhihuaensis. The expansion of MSH1 repeatedly occurred in multiple clades, especially in Gymnosperms, Orchidaceae, and Chloranthus spicatus. MSH1 has exceptionally long introns in certain species due to the gene length expansion, and the longest intron even reaches 101,025 bp. And the gene length is positively correlated with the proportion of the transposable elements (TEs) in the introns. In addition, gene structure analysis indicated that the MSH1 of green plants had undergone parallel intron gains and losses in all major lineages. However, the intron number of seed plants (gymnosperm and angiosperm) is relatively stable. All the selected gymnosperms contain 22 introns except for Gnetum montanum and Welwitschia mirabilis, while all the selected angiosperm species preserve 21 introns except for the ANA grade. Notably, the coding region of MSH1 in algae presents an exceptionally high GC content (47.7% to 75.5%). Moreover, over one-third of the selected species contain species-specific partial gene duplications of MSH1, except for the conserved mosses-specific partial gene duplication. Additionally, we found conserved alternatively spliced MSH1 transcripts in five species. The study of MSH1 sheds light on the evolution of the long genes of green plants.

NF-YA isoforms with alternative splicing of exon-5 in Aves

NF-YA, the regulatory subunit of the trimeric CCAAT-binding transcription factor NF-Y, is present in vertebrates in two major alternative spliced isoforms: NF-YAl and NF-YAs, differing for the presence of exon-3. NF-YAx, a third isoform without exon-3/-5, was reported only in human neuronal cells and tumors. These events affect the Trans-Activation Domain. We provide here evidence for the expression of NF-YAx and for the existence of a new isoform, NF-YAg, skipping only exon-5. These isoforms are abundant in Aves, but not in reptiles, and are the prevalent transcripts in the initial phases of embryo development in chicken. Finally, we analyzed NF-YAg and NF-YAx amino acid sequence using AlphaFold: absence of exon-5 denotes a global reduction of β-stranded elements, while removal of the disordered exon-3 sequence has limited effects on TAD architecture. These data identify an expanded program of NF-YA isoforms within the TAD in Aves, implying a role during early development.

Transposons contribute to the functional diversification of the head, gut, and ovary transcriptomes across Drosophila natural strains

Transcriptomes are dynamic, with cells, tissues, and body parts expressing particular sets of transcripts. Transposable elements (TEs) are a known source of transcriptome diversity; however, studies often focus on a particular type of chimeric transcript, analyze single body parts or cell types, or are based on incomplete TE annotations from a single reference genome. In this work, we have implemented a method based on de novo transcriptome assembly that minimizes the potential sources of errors while identifying a comprehensive set of gene-TE chimeras. We applied this method to the head, gut, and ovary dissected from five Drosophila melanogaster natural strains, with individual reference genomes available. We found that ∼19% of body part-specific transcripts are gene-TE chimeras. Overall, chimeric transcripts contribute a mean of 43% to the total gene expression, and they provide protein domains for DNA binding, catalytic activity, and DNA polymerase activity. Our comprehensive data set is a rich resource for follow-up analysis. Moreover, because TEs are present in virtually all species sequenced to date, their role in spatially restricted transcript expression is likely not exclusive to the species analyzed in this work.

Multifaceted framework for defining conservation units: An example from Atlantic salmon (Salmo salar) in Canada

Conservation units represent important components of intraspecific diversity that can aid in prioritizing and protecting at-risk populations, while also safeguarding unique diversity that can contribute to species resilience. In Canada, identification and assessments of conservation units is done by the Committee on the Status of Endangered Wildlife in Canada (COSEWIC). COSEWIC can recognize conservation units below the species level (termed "designatable units"; DUs) if the unit has attributes that make it both discrete and evolutionarily significant. There are various ways in which a DU can meet criteria of discreteness and significance, and increasing access to "big data" is providing unprecedented information that can directly inform both criteria. Specifically, the incorporation of genomic data for an increasing number of non-model species is informing more COSEWIC assessments; thus, a repeatable, robust framework is needed for integrating these data into DU characterization. Here, we develop a framework that uses a multifaceted, weight of evidence approach to incorporate multiple data types, including genetic and genomic data, to inform COSEWIC DUs. We apply this framework to delineate DUs of Atlantic salmon (Salmo salar, L.), an economically, culturally, and ecologically significant species, that is also characterized by complex hierarchical population structure. Specifically, we focus on an in-depth example of how our approach was applied to a previously data limited region of northern Canada that was defined by a single large DU. Application of our framework with newly available genetic and genomic data led to subdividing this DU into three new DUs. Although our approach was developed to meet criteria of COSEWIC, it is widely applicable given similarities in the definitions of a conservation unit.

A serpin gene from a parasitoid wasp disrupts host immunity and exhibits adaptive alternative splicing

Alternative splicing (AS) is a major source of protein diversity in eukaryotes, but less is known about its evolution compared to gene duplication (GD). How AS and GD interact is also largely understudied. By constructing the evolutionary trajectory of the serpin gene PpSerpin-1 (Pteromalus puparum serpin 1) in parasitoids and other insects, we found that both AS and GD jointly contribute to serpin protein diversity. These two processes are negatively correlated and show divergent features in both protein and regulatory sequences. Parasitoid wasps exhibit higher numbers of serpin protein/domains than nonparasitoids, resulting from more GD but less AS in parasitoids. The potential roles of AS and GD in the evolution of parasitoid host-effector genes are discussed. Furthermore, we find that PpSerpin-1 shows an exon expansion of AS compared to other parasitoids, and that several isoforms are involved in the wasp immune response, have been recruited to both wasp venom and larval saliva, and suppress host immunity. Overall, our study provides an example of how a parasitoid serpin gene adapts to parasitism through AS, and sheds light on the differential features of AS and GD in the evolution of insect serpins and their associations with the parasitic life strategy.

A brain specific alternatively spliced isoform of nonmuscle myosin IIA lacks its mechanoenzymatic activities

Recent genomic studies reported that 90 to 95% of human genes can undergo alternative splicing, by which multiple isoforms of proteins are synthesized. However, the functional consequences of most of the isoforms are largely unknown. Here, we report a novel alternatively spliced isoform of nonmuscle myosin IIA (NM IIA), called NM IIA2, which is generated by the inclusion of 21 amino acids near the actin-binding region (loop 2) of the head domain of heavy chains. Expression of NM IIA2 is found exclusively in the brain tissue, where it reaches a maximum level at 24 h during the circadian rhythm. The actin-dependent Mg2+-ATPase activity and in vitro motility assays reveal that NM IIA2 lacks its motor activities but localizes with actin filaments in cells. Interestingly, NM IIA2 can also make heterofilaments with NM IIA0 (noninserted isoform of NM IIA) and can retard the in vitro motility of NM IIA, when the two are mixed. Altogether, our findings provide the functional importance of a previously unknown alternatively spliced isoform, NM IIA2, and its potential physiological role in regulating NM IIA activity in the brain.

The Responses of Alternative Splicing during Heat Stress in the Pacific White Shrimp Litopenaeus vannamei

Heat tolerance is increasingly becoming a crucial trait for aquaculture species in the face of rapidly changing climate conditions. Alternative splicing (AS) is a vital mechanism within cells that modulates gene abundance and functional diversity, enabling organisms to effectively respond to diverse stressful conditions, including thermal stress. However, it is still uncertain whether AS contributes to heat tolerance in shrimp. In this study, we conducted an extensive transcriptome analysis on the Pacific white shrimp, Litopenaeus vannamei, revealing a total of 1267, 987, and 130 differential AS events (DAS) in the gill, hepatopancreas, and muscle, respectively, following exposure to heat stress. Among all of the DAS events, exon skipping (ES) was the predominant form of splicing modification observed. Interestingly, a minor portion of DAS genes exhibited overlap across the three tissues, implying that heat stress exerts unique effects on various tissue types. Moreover, the functional enrichment analysis demonstrated that commonly identified DAS genes were primarily associated with the "spliceosome" pathway, indicating that the AS of splicing-related genes played a crucial role in the response to heat stress. Our findings also revealed that heat stress tended to induce longer mRNA isoforms through differential alternative 3' splice site (A3SS) events. Notably, A3SS events exhibited the highest proportion of maintained open reading frames (ORFs) under heat stress. Interestingly, we observed a limited overlap between the genes exhibiting DAS and those showing differential gene expression (DEG), indicating that AS may function as a distinct regulatory mechanism independent of transcriptional regulation in response to heat stress. This is the first comprehensive study on AS in crustacea species under heat stress, which broadens our understanding of the regulatory mechanisms governing the crustaceans' response to environmental stress, providing valuable insights for the aquaculture breeding of shrimp and other aquatic animals.

A functional link between lariat debranching enzyme and the intron-binding complex is defective in non-photosensitive trichothiodystrophy

The pre-mRNA life cycle requires intron processing; yet, how intron-processing defects influence splicing and gene expression is unclear. Here, we find that TTDN1/MPLKIP, which is encoded by a gene implicated in non-photosensitive trichothiodystrophy (NP-TTD), functionally links intron lariat processing to spliceosomal function. The conserved TTDN1 C-terminal region directly binds lariat debranching enzyme DBR1, whereas its N-terminal intrinsically disordered region (IDR) binds the intron-binding complex (IBC). TTDN1 loss, or a mutated IDR, causes significant intron lariat accumulation, as well as splicing and gene expression defects, mirroring phenotypes observed in NP-TTD patient cells. A Ttdn1-deficient mouse model recapitulates intron-processing defects and certain neurodevelopmental phenotypes seen in NP-TTD. Fusing DBR1 to the TTDN1 IDR is sufficient to recruit DBR1 to the IBC and circumvents the functional requirement for TTDN1. Collectively, our findings link RNA lariat processing with splicing outcomes by revealing the molecular function of TTDN1.

Differential impact of BRAFV600E isoforms on tumorigenesis in a zebrafish model of melanoma

BRAFV600E comes as two main splicing variants. The well-studied ref isoform and the recently discovered X1 isoform are co-expressed in cancer cells and differ in terms of 3'UTR length and sequence, as well as C-term protein sequence. Here, we use a melanoma model in zebrafish to study the role played by each isoform in larval pigmentation, nevi formation, and their progression into melanoma tumours. We show that both BRAFV600E-ref and BRAFV600E-X1 proteins promote larval pigmentation and nevi formation, while melanoma-free survival curves performed in adult fish indicate that BRAFV600E-ref protein is a much stronger melanoma driver that BRAFV600E-X1 protein. Crucially, we also show that the presence of the 3'UTR suppresses the effect of ref protein. Our data highlight the necessity to undertake a systematic study of BRAFV600E isoforms, in order to uncover the full spectrum of their kinase-(in)dependent and coding-(in)dependent functions, hence to develop more informed strategies for therapeutic targeting.

The roles of sexual selection and sexual conflict in shaping patterns of genome and transcriptome variation

Sexual dimorphism is one of the most prevalent, and often the most extreme, examples of phenotypic variation within species, and arises primarily from genomic variation that is shared between females and males. Many sexual dimorphisms arise through sex differences in gene expression, and sex-biased expression is one way that a single, shared genome can generate multiple, distinct phenotypes. Although many sexual dimorphisms are expected to result from sexual selection, and many studies have invoked the possible role of sexual selection to explain sex-specific traits, the role of sexual selection in the evolution of sexually dimorphic gene expression remains difficult to differentiate from other forms of sex-specific selection. In this Review, we propose a holistic framework for the study of sex-specific selection and transcriptome evolution. We advocate for a comparative approach, across tissues, developmental stages and species, which incorporates an understanding of the molecular mechanisms, including genomic variation and structure, governing gene expression. Such an approach is expected to yield substantial insights into the evolution of genetic variation and have important applications in a variety of fields, including ecology, evolution and behaviour.

Resurrecting the alternative splicing landscape of archaic hominins using machine learning

Alternative splicing contributes to adaptation and divergence in many species. However, it has not been possible to directly compare splicing between modern and archaic hominins. Here, we unmask the recent evolution of this previously unobservable regulatory mechanism by applying SpliceAI, a machine-learning algorithm that identifies splice-altering variants (SAVs), to high-coverage genomes from three Neanderthals and a Denisovan. We discover 5,950 putative archaic SAVs, of which 2,186 are archaic-specific and 3,607 also occur in modern humans via introgression (244) or shared ancestry (3,520). Archaic-specific SAVs are enriched in genes that contribute to traits potentially relevant to hominin phenotypic divergence, such as the epidermis, respiration and spinal rigidity. Compared to shared SAVs, archaic-specific SAVs occur in sites under weaker selection and are more common in genes with tissue-specific expression. Further underscoring the importance of negative selection on SAVs, Neanderthal lineages with low effective population sizes are enriched for SAVs compared to Denisovan and shared SAVs. Finally, we find that nearly all introgressed SAVs in humans were shared across the three Neanderthals, suggesting that older SAVs were more tolerated in human genomes. Our results reveal the splicing landscape of archaic hominins and identify potential contributions of splicing to phenotypic differences among hominins.

Multidimensional plasticity jointly contributes to rapid acclimation to environmental challenges during biological invasions

Rapid plastic response to environmental changes, which involves extremely complex underlying mechanisms, is crucial for organismal survival during many ecological and evolutionary processes such as those in global change and biological invasions. Gene expression is among the most studied molecular plasticity, while co- or posttranscriptional mechanisms are still largely unexplored. Using a model invasive ascidian Ciona savignyi, we studied multidimensional short-term plasticity in response to hyper- and hyposalinity stresses, covering the physiological adjustment, gene expression, alternative splicing (AS), and alternative polyadenylation (APA) regulations. Our results demonstrated that rapid plastic response varied with environmental context, timescales, and molecular regulatory levels. Gene expression, AS, and APA regulations independently acted on different gene sets and corresponding biological functions, highlighting their nonredundant roles in rapid environmental adaptation. Stress-induced gene expression changes illustrated the use of a strategy of accumulating free amino acids under high salinity and losing/reducing them during low salinity to maintain the osmotic homoeostasis. Genes with more exons were inclined to use AS regulations, and isoform switches in functional genes such as SLC2a5 and Cyb5r3 resulted in enhanced transporting activities by up-regulating the isoforms with more transmembrane regions. The extensive 3'-untranslated region (3'UTR) shortening through APA was induced by both salinity stresses, and APA regulation predominated transcriptomic changes at some stages of stress response. The findings here provide evidence for complex plastic mechanisms to environmental changes, and thereby highlight the importance of systemically integrating different levels of regulatory mechanisms in studying initial plasticity in evolutionary trajectories.

Alternative splicing in ABA signaling during seed germination

Seed germination is an essential step in a plant's life cycle. It is controlled by complex physiological, biochemical, and molecular mechanisms and external factors. Alternative splicing (AS) is a co-transcriptional mechanism that regulates gene expression and produces multiple mRNA variants from a single gene to modulate transcriptome diversity. However, little is known about the effect of AS on the function of generated protein isoforms. The latest reports indicate that alternative splicing (AS), the relevant mechanism controlling gene expression, plays a significant role in abscisic acid (ABA) signaling. In this review, we present the current state of the art about the identified AS regulators and the ABA-related changes in AS during seed germination. We show how they are connected with the ABA signaling and the seed germination process. We also discuss changes in the structure of the generated AS isoforms and their impact on the functionality of the generated proteins. Also, we point out that the advances in sequencing technology allow for a better explanation of the role of AS in gene regulation by more accurate detection of AS events and identification of full-length splicing isoforms.

HBV Enhances Sorafenib Resistance in Hepatocellular Carcinoma by Reducing Ferroptosis via SRSF2-Mediated Abnormal PCLAF Splicing

Hepatocellular carcinoma (HCC) is one of the most lethal human cancers. Hepatitis B virus (HBV) infection accounts for nearly 50% of HCC cases. Recent studies indicate that HBV infection induces resistance to sorafenib, the first-line systemic treatment for advanced HCC for more than a decade, from 2007 to 2020. Our previous research shows that variant 1 (tv1) of proliferating cell nuclear antigen clamp-associated factor (PCLAF), overexpressed in HCC, protects against doxorubicin-induced apoptosis. However, there are no reports on the relevance of PCLAF in sorafenib resistance in HBV-related HCC. In this article, we found that PCLAF levels were higher in HBV-related HCC than in non-virus-related HCC using bioinformatics analysis. Immunohistochemistry (IHC) staining of clinical samples and the splicing reporter minigene assay using HCC cells revealed that PCLAF tv1 was elevated by HBV. Furthermore, HBV promoted the splicing of PCLAF tv1 by downregulating serine/arginine-rich splicing factor 2 (SRSF2), which hindered the inclusion of PCLAF exon 3 through a putative cis-element (116-123), "GATTCCTG". The CCK-8 assay showed that HBV decreased cell susceptibility to sorafenib through SRSF2/PCLAF tv1. HBV reduced ferroptosis by decreasing intracellular Fe²⁺ levels and activating GPX4 expression via the SRSF2/PCLAF tv1 axis, according to a mechanism study. Suppressed ferroptosis, on the other hand, contributed to HBV-mediated sorafenib resistance through SRSF2/PCLAF tv1. These data suggested that HBV regulated PCLAF abnormal alternative splicing by suppressing SRSF2. HBV caused sorafenib resistance by reducing ferroptosis via the SRSF2/PCLAF tv1 axis. As a result, the SRSF2/PCLAF tv1 axis may be a prospective molecular therapeutic target in HBV-related HCC, as well as a predictor of sorafenib resistance. The inhibition of the SRSF2/PCLAF tv1 axis may be crucial in the emergence of systemic chemotherapy resistance in HBV-associated HCC.

The landscape of expression and alternative splicing variation across human traits

Understanding the consequences of individual transcriptome variation is fundamental to deciphering human biology and disease. We implement a statistical framework to quantify the contributions of 21 individual traits as drivers of gene expression and alternative splicing variation across 46 human tissues and 781 individuals from the Genotype-Tissue Expression project. We demonstrate that ancestry, sex, age, and BMI make additive and tissue-specific contributions to expression variability, whereas interactions are rare. Variation in splicing is dominated by ancestry and is under genetic control in most tissues, with ribosomal proteins showing a strong enrichment of tissue-shared splicing events. Our analyses reveal a systemic contribution of types 1 and 2 diabetes to tissue transcriptome variation with the strongest signal in the nerve, where histopathology image analysis identifies novel genes related to diabetic neuropathy. Our multi-tissue and multi-trait approach provides an extensive characterization of the main drivers of human transcriptome variation in health and disease.

FishExp: A comprehensive database and analysis platform for gene expression and alternative splicing of fish species

The publicly archived RNA-seq data has grown exponentially, while its valuable information has not yet been fully discovered and utilized, such as alternative splicing and its integration with gene expression. This is especially true for fish species which play important roles in ecology, research and the food industry. Furthermore, there is a lack of online platform to analyze users’ new data individually and jointly with existing data for the comprehensive analysis of alternative splicing and gene expression. Here, we present FishExp, a web-based data platform covering gene expression and alternative splicing in 26,081 RNA-seq experiments from 44 fishes. It allows users to query the data in a variety of ways, including gene identifier/symbol, functional term, and BLAST alignment. Moreover, users can customize experiments and tools to perform differential/specific expression and alternative splicing analysis, co-expression and cross-species analysis. In addition, functional enrichment is provided to confer biological significance. Notably, users are allowed to submit their own data and perform various analyses using the new data alone or alongside existing data in FishExp. Results of retrieval and analysis can be visualized on the gene-, transcript- and splicing event-level webpage in a highly interactive and intuitive manner. All data in FishExp can be downloaded for more in-depth analysis. The manually curated sample information, uniform data processing and various tools make it efficient for users to gain new insights from these large data sets, facilitating scientific hypothesis generation. FishExp is freely accessible at https://bioinfo.njau.edu.cn/fishExp.

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.

Genomic reaction norms inform predictions of plastic and adaptive responses to climate change

Genomic reaction norms represent the range of gene expression phenotypes (usually mRNA transcript levels) expressed by a genotype along an environmental gradient. Reaction norms derived from common‐garden experiments are powerful approaches for disentangling plastic and adaptive responses to environmental change in natural populations. By treating gene expression as a phenotype in itself, genomic reaction norms represent invaluable tools for exploring causal mechanisms underlying organismal responses to climate change across multiple levels of biodiversity.

Our goal is to provide the context, framework and motivation for applying genomic reaction norms to study the responses of natural populations to climate change.

Here, we describe the utility of integrating genomics with common‐garden‐gradient experiments under a reaction norm analytical framework to answer fundamental questions about phenotypic plasticity, local adaptation, their interaction (i.e. genetic variation in plasticity) and future adaptive potential.

An experimental and analytical framework for constructing and analysing genomic reaction norms is presented within the context of polygenic climate change responses of structured populations with gene flow. Intended for a broad eco‐evo readership, we first briefly review adaptation with gene flow and the importance of understanding the genomic basis and spatial scale of adaptation for conservation and management of structured populations under anthropogenic change. Then, within a high‐dimensional reaction norm framework, we illustrate how to distinguish plastic, differentially expressed (difference in reaction norm intercepts) and differentially plastic (difference in reaction norm slopes) genes, highlighting the areas of opportunity for applying these concepts.

We conclude by discussing how genomic reaction norms can be incorporated into a holistic framework to understand the eco‐evolutionary dynamics of climate change responses from molecules to ecosystems. We aim to inspire researchers to integrate gene expression measurements into common‐garden experimental designs to investigate the genomics of climate change responses as sequencing costs become increasingly accessible.

The Importance of a Genome-Wide Association Analysis in the Study of Alternative Splicing Mutations in Plants with a Special Focus on Maize

Alternative splicing is an important mechanism for regulating gene expressions at the post-transcriptional level. In eukaryotes, the genes are transcribed in the nucleus to produce pre-mRNAs and alternative splicing can splice a pre-mRNA to eventually form multiple different mature mRNAs, greatly increasing the number of genes and protein diversity. Alternative splicing is involved in the regulation of various plant life activities, especially the response of plants to abiotic stresses and is also an important process of plant growth and development. This review aims to clarify the usefulness of a genome-wide association analysis in the study of alternatively spliced variants by summarizing the application of alternative splicing, genome-wide association analyses and genome-wide association analyses in alternative splicing, as well as summarizing the related research progress.