The importance of alternative splicing in adaptive evolution

Emerging Orchestrator of Ecological Adaptation: m6A Regulation of Post-Transcriptional Mechanisms

Genetic mechanisms have been at the forefront of our exploration into the substrate of adaptive evolution and phenotypic diversification. However, genetic variation only accounts for a fraction of phenotypic variation. In the last decade, the significance of RNA modification mechanisms has become more apparent in the context of organismal adaptation to rapidly changing environments. RNA m6A methylation, the most abundant form of RNA modification, is emerging as a potentially significant player in various biological processes. Despite its fundamental function to regulate other major post-transcriptional mechanisms such as microRNA and alternative splicing, its role in ecology and evolution has been understudied. This review highlights the potential importance of m6A RNA methylation in ecological adaptation, emphasising the need for further research, especially in natural systems. We focus on how m6A not only affects mRNA fate but also influences miRNA-mediated gene regulation and alternative splicing, potentially contributing to organismal adaptation. The aim of this review is to synthesise key background information to enhance our understanding of m6A mechanisms driving species survival in dynamic environments and motivate future research into the dynamics of adaptive RNA methylation.

The roles of different gene expression regulators in acoustic variation in the intermediate horseshoe bat revealed by long-read and short-read RNA sequencing data

Gene expression changes contribute greatly to phenotypic variations in nature. Studying patterns of regulators of gene expression is important to fully understand the molecular mechanism underlying phenotypic variations. In horseshoe bats, the cochleae are finely tuned to echoes of call frequency. Here, using 2 recently diverged subspecies of the intermediate horseshoe bat (Rhinolophus affinis hainanus and R. a. himalayanus) with great acoustic variations as the system, we aim to explore relative roles of different regulators of gene expression (differential gene expression, alternative splicing (AS) and long non-coding RNAs (lncRNAs)) in phenotypic variation with a combination of Illumina short-read and Nanopore long-read RNA-seq data from the cochlea. Compared to R. a. hainanus, R. a. himalayanus exhibited much more upregulated differentially expressed genes (DEGs) and multiple of them may play important roles in the maintenance and damage repair of auditory hair cells. We identified 411 differentially expressed lncRNAs and their target DEGs upregulated in R. a. himalayanus were also mainly involved in a protective mechanism for auditory hair cells. Using 3 different methods of AS analysis, we identified several candidate alternatively spliced genes (ASGs) that expressed different isoforms which may be associated with acoustic divergence of the 2 subspecies. We observed significantly less overlap than expected between DEGs and ASGs, supporting complementary roles of differential gene expression and AS in generating phenotypic variations. Overall, our study highlights the importance of a combination of short-read and long-read RNA-seq data in examining the regulation of gene expression changes responsible for phenotypic variations.

A global survey of bicarbonate stress-induced pre-mRNA alternative splicing in soybean via integrative analysis of Iso-seq and RNA-seq

Alternative splicing (AS) plays important roles in modulating environmental stress responses in plants. However, little is known about the functions of bicarbonate-induced AS in cultivated soybean (Glycine max L. Merr.). In this study, we combined PacBio isoform sequencing (Iso-seq) and Illumina RNA sequencing (RNA-seq) to elucidate the bicarbonate-induced AS events in soybean root and leaf tissues. Compared to RNA-seq, Iso-seq identified more novel genes and transcripts, as well as more AS events, indicating that Iso-seq is more efficient in AS detection. Combining these two technologies, we found that intron retention (IR) is the most frequent AS event type. We identified a total of 913 and 1974 bicarbonate stress-responsive differentially alternative spliced genes (DAGs) in soybean leaves and roots respectively, from our RNA-seq results. Additionally, we determined a transcription factor (GmNTL9) and a splicing factor (GmRSZ22), and validated their roles in bicarbonate stress response by AS. Overall, our study opens an avenue for evaluating plant AS regulatory networks, and the obtained global landscape of alternative splicing provides valuable insights into the AS-mediated bicarbonate-responsive mechanisms in plant species.

A chromosome-level reference genome assembly and a full-length transcriptome assembly of the giant freshwater prawn (Macrobrachium rosenbergii)

The giant freshwater prawn (Macrobrachium rosenbergii) is a key species in the aquaculture industry in several Asian, African, and South American countries. Despite a considerable growth in its production worldwide, the genetic complexities of M. rosenbergii various morphotypes pose challenges in cultivation. This study reports the first chromosome-scale reference genome and a high-quality full-length transcriptome assembly for M. rosenbergii. We employed the PacBio High Fidelity (HiFi) sequencing to obtain an initial draft assembly and further scaffolded it with the chromatin contact mapping (Hi-C) technique to achieve a final assembly of 3.73-Gb with an N50 scaffold length of 33.6 Mb. Repetitive elements constituted nearly 60% of the genome assembly, with simple sequence repeats and retrotransposons being the most abundant. The availability of both the chromosome-scale assembly and the full-length transcriptome assembly enabled us to thoroughly probe alternative splicing events in M. rosenbergii. Among the 2,041 events investigated, exon skipping represented the most prevalent class, followed by intron retention. Interestingly, specific isoforms were observed across multiple tissues. Additionally, within a single tissue type, transcripts could undergo alternative splicing, yielding multiple isoforms. We believe that the availability of a chromosome-level reference genome for M. rosenbergii, along with its full-length transcriptome, will be instrumental in advancing our understanding of the giant freshwater prawn biology and enhancing its molecular breeding programs, paving the way for the development of M. rosenbergii with valuable traits in commercial aquaculture.

New insights into the dynamics of m6A epitranscriptome: hybrid-seq identifies novel mRNAs of the m6A writers METTL3/14

Background: N6-methyladenosine (m6A), a prevalent mRNA modification, is dynamically regulated by methyltransferases, including METTL3 and METTL14.Materials & methods: In the current study, we employed a custom hybrid-seq method to identify novel METTL3/14 transcripts, explore their protein-coding capacities and predict the putative role of the METTL isoforms.Results: Demultiplexing of the hybrid-seq barcoded datasets unraveled the expression patterns of the newly identified mRNAs in major malignancies as well as in non-malignant cells, providing a deeper understanding of the methylation pathways. Open reading frame query revealed novel METTL3/14 isoforms, broadening our perspective for the structural diversity within METTL family.Conclusion: Our findings offer significant insights into the intricate transcriptional landscape of METTL3/14, shedding light on the regulatory mechanisms underlying methylation in mRNAs.

When splicing meets transgene DNA methylation to uncover elite alleles for plant yield

This article is a Commentary on Zhao et al. (2024), 243: 1724–1741.

Alternative splicing in parallel evolution and the evolutionary potential in sticklebacks

Repeatability of adaptation to similar environments provides opportunity to evaluate the predictability of natural selection. While many studies have investigated gene expression differences between populations adapted to contrasting environments, the role of post-transcriptional processes such as alternative splicing has rarely been evaluated in the context of parallel adaptation. To address the aforementioned knowledge gap, we reanalysed transcriptomic data from three pairs of threespine stickleback (Gasterosteus aculeatus) ecotypes adapted to marine or freshwater environment. First, we identified genes with repeated expression or splicing divergence across ecotype pairs, and compared the genetic architecture and biological processes between parallelly expressed and parallelly spliced loci. Second, we analysed the extent to which parallel adaptation was reflected at gene expression and alternative splicing levels. Finally, we tested how the two axes of transcriptional variation differed in their potential for evolutionary change. Although both repeated differential splicing and differential expression across ecotype pairs showed tendency for parallel divergence, the degree of parallelism was lower for splicing than expression. Furthermore, parallel divergences in splicing and expression were likely to be associated with distinct cis-regulatory genetic variants and functionally unique set of genes. Finally, we found that parallelly spliced genes showed higher nucleotide diversity than parallelly expressed genes, indicating splicing is less susceptible to genetic variation erosion during parallel adaptation. Our results provide novel insight into the role of splicing in parallel adaptation, and underscore the contribution of splicing to the evolutionary potential of wild populations under environmental change.

bZIP Transcription Factors: Structure, Modification, Abiotic Stress Responses and Application in Plant Improvement

Plant growth, yield, and distribution are significantly impacted by abiotic stresses, affecting global ecosystems and forestry practices. However, plants have evolved complex adaptation mechanisms governed by numerous genes and transcription factors (TFs) to manage these stresses. Among these, bZIP (basic leucine zipper) is a crucial regulator orchestrating morphological adaptations. This review aims to elucidate the multifaceted roles of bZIP TFs in plant species. We discuss the morphological changes induced by stress stimuli and the pivotal functions of bZIP TFs in mediating these responses. While several publications have explored the mechanisms of bZIP TFs in response to abiotic stresses, this review delves into the intricate regulatory networks, summarizing alternative splicing and post-translational modifications, signaling networks interacting with bZIP TFs, and genetic engineering of bZIP TFs. By synthesizing current research, this review provides an updated discussion on bZIP interactions with other proteins to regulate stresses such as cold, heat, drought, and salt. Additionally, it offers avenues for future research and applications of bZIP TFs to improve abiotic stress resilience in plants through genetic engineering.

On the resolution of sexual conflict over shared traits

Anisogamy, different-sized male and female gametes, sits at the heart of sexual selection and conflict between the sexes. Sperm producers (males) and egg producers (females) of the same species generally share most, if not all, of the same genome, but selection frequently favours different trait values in each sex for traits common to both. The extent to which this conflict might be resolved, and the potential mechanisms by which this can occur, have been widely debated. Here, we summarize recent findings and emphasize that once the sexes evolve, sexual selection is ongoing, and therefore new conflict is always possible. In addition, sexual conflict is largely a multivariate problem, involving trait combinations underpinned by networks of interconnected genes. Although these complexities can hinder conflict resolution, they also provide multiple possible routes to decouple male and female phenotypes and permit sex-specific evolution. Finally, we highlight difficulty in the study of sexual conflict over shared traits and promising directions for future research.

Understanding the role of temperature in seasonal timing: Effects on behavioural, physiological and molecular phenotypes

Organisms adapt to daily and seasonal environmental changes to maximise their metabolic and reproductive fitness. For seasonally breeding animals, photoperiod is considered the most robust cue to drive these changes. It, however, does not explain the interannual variations in different seasonal phenotypes. Several studies have repeatedly shown the influence of ambient temperature on the timing of different seasonal physiologies including the timing of migration, reproduction and its associated behaviours, etc. In the present review, we have discussed the effects of changes in ambient temperature on different seasonal events in endotherms with a focus on migratory birds as they have evolved to draw benefits from distinct but largely predictable seasonal patterns of natural resources. We have further discussed the physiological and molecular mechanisms by which temperature affects seasonal timings. The primary brain area involved in detecting temperature changes is the hypothalamic preoptic area. This area receives thermal inputs via sensory neurons in the peripheral ganglia that measure changes in thermoregulatory tissues such as the skin and spinal cord. For the input signals, several thermal sensory TRP (transient receptor potential ion channels) channels have been identified across different classes of vertebrates. These channels are activated at specific thermal ranges. Once perceived, this information should activate an effector function. However, the link between temperature sensation and the effector pathways is not properly understood yet. Here, we have summarised the available information that may help us understand how temperature information is translated into seasonal timing.

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.

Genome-Wide Identification of Glutathione S-Transferase Genes in Eggplant (Solanum melongena L.) Reveals Their Potential Role in Anthocyanin Accumulation on the Fruit Peel

Anthocyanins are ubiquitous pigments derived from the phenylpropanoid compound conferring red, purple and blue pigmentations to various organs of horticultural crops. The metabolism of flavonoids in the cytoplasm leads to the biosynthesis of anthocyanin, which is then conveyed to the vacuoles for storage by plant glutathione S-transferases (GST). Although GST is important for transporting anthocyanin in plants, its identification and characterization in eggplant (Solanum melongena L.) remains obscure. In this study, a total of 40 GST genes were obtained in the eggplant genome and classified into seven distinct chief groups based on the evolutionary relationship with Arabidopsis thaliana GST genes. The seven subgroups of eggplant GST genes (SmGST) comprise: dehydroascorbate reductase (DHAR), elongation factor 1Bγ (EF1Bγ), Zeta (Z), Theta(T), Phi(F), Tau(U) and tetra-chlorohydroquinone dehalogenase TCHQD. The 40 GST genes were unevenly distributed throughout the 10 eggplant chromosomes and were predominantly located in the cytoplasm. Structural gene analysis showed similarity in exons and introns within a GST subgroup. Six pairs of both tandem and segmental duplications have been identified, making them the primary factors contributing to the evolution of the SmGST. Light-related cis-regulatory elements were dominant, followed by stress-related and hormone-responsive elements. The syntenic analysis of orthologous genes indicated that eggplant, Arabidopsis and tomato (Solanum lycopersicum L.) counterpart genes seemed to be derived from a common ancestry. RNA-seq data analyses showed high expression of 13 SmGST genes with SmGSTF1 being glaringly upregulated on the peel of purple eggplant but showed no or low expression on eggplant varieties with green or white peel. Subsequently, SmGSTF1 had a strong positive correlation with anthocyanin content and with anthocyanin structural genes like SmUFGT (r = 0.9), SmANS (r = 0.85), SmF3H (r = 0.82) and SmCHI2 (r = 0.7). The suppression of SmGSTF1 through virus-induced gene silencing (VIGs) resulted in a decrease in anthocyanin on the infiltrated fruit surface. In a nutshell, results from this study established that SmGSTF1 has the potential of anthocyanin accumulation in eggplant peel and offers viable candidate genes for the improvement of purple eggplant. The comprehensive studies of the SmGST family genes provide the foundation for deciphering molecular investigations into the functional analysis of SmGST genes in eggplant.

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.

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.

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.

Alternative splicing variants of stimulator of interferon genes (STING) from Asian seabass (Lates calcarifer) and their immune response against red spotted grouper nervous necrosis virus (RGNNV)

The Stimulator of Interferon Genes (STING, also known as MITA/ERYS/MPYS) is an adaptor molecule that plays a crucial role in the RLR pathway and responds to DNA and RNA viruses. In the present study, we have identified two novel isoforms of STING (the canonical form named as LcSTINGa and its alternative splicing isoform named as LcSTINGb) from teleost Lates calcarifer. LcSTINGa has an ORF of 1230 bp, encoding a 409 amino acid protein, while its alternative splicing variant, LcSTINGb, features an ORF of 987 bp, encoding 328 amino acids. LcSTINGa is predicted to contain four transmembrane helices, whereas LcSTINGb has only two. The Lates STING protein showed about 86.85% identity with Perca flavescens, 86.45% with Seriola and 39.51% with Homo sapiens. The tissue distribution studies revealed that the STING variants were constitutively expressed in all the tissues examined, with the highest expression in blood. In-vivo upregulation of LcSTINGa and LcSTINGb mRNA following immune challenge with poly (I:C), Red-spotted grouper nervous necrosis virus (RGNNV) and zymosan A suggests its significance in the immune response.

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.

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.

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.

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.

Differential gene expression in response to fungal pathogen exposure in the aquatic invertebrate, Daphnia dentifera

While vertebrate immune systems are appreciated for their complexity and adaptability, invertebrate immunity is often considered to be less complex. However, immune responses in many invertebrates likely involve sophisticated processes. Interactions between the crustacean host Daphnia dentifera and its fungal pathogen Metschnikowia bicuspidata provide an excellent model for exploring the mechanisms underlying crustacean immunity. To explore the genomic basis of immunity in Daphnia, we used RNA-sequencing technology to quantify differential gene expression between individuals of a single host genotype exposed or unexposed to M. bicuspidata over 24 h. Transcriptomic analyses showed that the number of differentially expressed genes between the control (unexposed) and experimental (exposed) groups increased over time. Gene ontology enrichment analysis revealed that differentially expressed genes were enriched for immune-related molecules and processes, such as cuticle development, prostaglandin, and defense response processes. Our findings provide a suite of immunologically relevant genes and suggest the presence of a rapidly upregulated immune response involving the cuticle in Daphnia. Studies involving gene expression responses to pathogen exposure shine a light on the processes occurring during the course of infection. By leveraging knowledge on the genetic basis for immunity, immune mechanisms can be more thoroughly understood to refine our understanding of disease spread within invertebrate populations.

Alternative Splicing during Fiber Development in G. hirsutum

Cotton is a valuable cash crop in many countries. Cotton fiber is a trichome that develops from a single epidermal cell and serves as an excellent model for understanding cell differentiation and other life processes. Alternative splicing (AS) of genes is a common post-transcriptional regulatory process in plants that is essential for plant growth and development. The process of AS during cotton fiber formation, on the other hand, is mainly unknown. A substantial number of multi-exon genes were discovered to be alternatively spliced during cotton fiber formation in this study, accounting for 23.31% of the total number of genes in Gossypium hirsutum. Retention intron (RI) is not necessarily the most common AS type, indicating that AS genes and processes during fiber development are very temporal and tissue-specific. When compared to fiber samples, AS is more prevalent at the fiber initiation stages and in the ovule, indicating that development stages and tissues use different AS strategies. Genes involved in fiber development have gone through stage-specific AS, demonstrating that AS regulates cotton fiber development. Furthermore, AS can be regulated by trans-regulation elements such as splicing factor and cis-regulation elements such as gene length, exon numbers, and GC content, particularly at exon-intron junction sites. Our findings also suggest that increased DNA methylation may aid in the efficiency of AS, and that gene body methylation is key in AS control. Finally, our research will provide useful information about the roles of AS during the cotton fiber development process.

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.

A first annotated genome sequence for Haliotis midae with genomic insights into abalone evolution and traits of economic importance

Haliotis midae or "perlemoen" is one of five abalone species endemic to South Africa, and being palatable, the only commercially important abalone species with a high international demand. The higher demand for this abalone species has resulted in the decrease of natural stocks due to overexploitation by capture fisheries and poaching. Facilitating aquaculture production of H. midae should assist in minimising the pressure on the wild populations. Here, the draft genome of H. midae has been sequenced, assembled, and annotated. The draft assembly resulted in a total length of 1.5 Gb, contig N50 of 0.238 Mb, scaffold N50 of 0. 238 Mb and GC level of 40%. Gene annotation, combining ab initio and evidence-based pipelines identified 52,280 genes with protein coding potential. The genes identified were used to predict orthologous genes shared among the four other abalone species (H. laevigata, H. rubra, H. discus hannai and H. rufescens) and 4702 orthologous genes were shared across the five species. Among the orthologous genes in abalones, single copy genes were further analysed for signatures of selection and several molecular regulatory proteins involved in developmental functions were found to be under positive selection in specific abalone lineages. Furthermore, whole genome SNP-based phylogenomic assessment was performed to confirm the evolutionary relationship among the considered abalone species with draft genomes, reaffirming that H. midae is closely related to the Australian Greenlip (H. laevigata) and Blacklip (H. rubra). The study assists in the understanding of genes related to various biological systems underscoring the evolution and development of abalones, with potential applications for genetic improvement of commercial stocks.

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.

PlaASDB: a comprehensive database of plant alternative splicing events in response to stress

BACKGROUND

Alternative splicing (AS) is a co-transcriptional regulatory mechanism of plants in response to environmental stress. However, the role of AS in biotic and abiotic stress responses remains largely unknown. To speed up our understanding of plant AS patterns under different stress responses, development of informative and comprehensive plant AS databases is highly demanded.

DESCRIPTION

In this study, we first collected 3,255 RNA-seq data under biotic and abiotic stresses from two important model plants (Arabidopsis and rice). Then, we conducted AS event detection and gene expression analysis, and established a user-friendly plant AS database termed PlaASDB. By using representative samples from this highly integrated database resource, we compared AS patterns between Arabidopsis and rice under abiotic and biotic stresses, and further investigated the corresponding difference between AS and gene expression. Specifically, we found that differentially spliced genes (DSGs) and differentially expressed genes (DEG) share very limited overlapping under all kinds of stresses, suggesting that gene expression regulation and AS seemed to play independent roles in response to stresses. Compared with gene expression, Arabidopsis and rice were more inclined to have conserved AS patterns under stress conditions.

CONCLUSION

PlaASDB is a comprehensive plant-specific AS database that mainly integrates the AS and gene expression data of Arabidopsis and rice in stress response. Through large-scale comparative analyses, the global landscape of AS events in Arabidopsis and rice was observed. We believe that PlaASDB could help researchers understand the regulatory mechanisms of AS in plants under stresses more conveniently. PlaASDB is freely accessible at http://zzdlab.com/PlaASDB/ASDB/index.html .

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

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

Genomic basis of Y-linked dwarfism in cichlids pursuing alternative reproductive tactics

Sexually antagonistic selection, which favours different optima in males and females, is predicted to play an important role in the evolution of sex chromosomes. Body size is a sexually antagonistic trait in the shell-brooding cichlid fish Lamprologous callipterus, as "bourgeois" males must be large enough to carry empty snail shells to build nests whereas females must be small enough to fit into shells for breeding. In this species, there is also a second male morph: smaller "dwarf" males employ an alternative reproductive strategy by wriggling past spawning females into shells to fertilize eggs. L. callipterus male morphology is passed strictly from father to son, suggesting Y-linkage. However, sex chromosomes had not been previously identified in this species, and the genomic basis of size dimorphism was unknown. Here we used whole-genome sequencing to identify a 2.4-Mb sex-linked region on scaffold_23 with reduced coverage and single nucleotide polymorphism density in both male morphs compared to females. Within this sex region, distinct Y-haplotypes delineate the two male morphs, and candidate genes for body size (GHRHR, a known dwarfism gene) and sex determination (ADCYAP1R1) are in high linkage disequilibrium. Because differences in body size between females and males are under strong selection in L. callipterus, we hypothesize that sexual antagonism over body size initiated early events in sex chromosome evolution, followed by Y divergence to give rise to bourgeois and dwarf male reproductive strategies. Our results are consistent with the hypothesis that sexually antagonistic traits should be linked to young sex chromosomes.

Temporal Analysis of Gene Expression and Isoform Switching in Brown Bears (Ursus arctos)

Hibernation in brown bears is an annual process involving multiple physiologically distinct seasons-hibernation, active, and hyperphagia. While recent studies have characterized broad patterns of differential gene regulation and isoform usage between hibernation and active seasons, patterns of gene and isoform expression during hyperphagia remain relatively poorly understood. The hyperphagia stage occurs between active and hibernation seasons and involves the accumulation of large fat reserves in preparation for hibernation. Here, we use time-series analyses of gene expression and isoform usage to interrogate transcriptomic regulation associated with all three seasons. We identify a large number of genes with significant differential isoform usage (DIU) across seasons and show that these patterns of isoform usage are largely tissue-specific. We also show that DIU and differential gene-level expression responses are generally non-overlapping, with only a small subset of multi-isoform genes showing evidence of both gene-level expression changes and changes in isoform usage across seasons. Additionally, we investigate nuanced regulation of candidate genes involved in the insulin signaling pathway and find evidence of hyperphagia-specific gene expression and isoform regulation that may enhance fat accumulation during hyperphagia. Our findings highlight the value of using temporal analyses of both gene- and isoform-level gene expression when interrogating complex physiological phenotypes and provide new insight into the mechanisms underlying seasonal changes in bear physiology.

Relocation to avoid costs: A hypothesis on red carotenoid-based signals based on recent CYP2J19 gene expression data

In many vertebrates, the enzymatic oxidation of dietary yellow carotenoids generates red keto-carotenoids giving color to ornaments. The oxidase CYP2J19 is here a key effector. Its purported intracellular location suggests a shared biochemical pathway between trait expression and cell functioning. This might guarantee the reliability of red colorations as individual quality signals independent of production costs. We hypothesize that the ornament type (feathers vs. bare parts) and production costs (probably CYP2J19 activity compromising vital functions) could have promoted tissue-specific gene relocation. We review current avian tissue-specific CYP2J19 expression data. Among the ten red-billed species showing CYP2J19 bill expression, only one showed strong hepatic expression. Moreover, a phylogenetically-controlled analysis of 25 red-colored species shows that those producing red bare parts are less likely to have strong hepatic CYP2J19 expression than species with only red plumages. Thus, both production costs and shared pathways might have contributed to the evolution of red signals.

High-quality reference transcriptome construction improves RNA-seq quantification in Oryza sativa indica

The Reference Transcriptomic Dataset (RTD) is an accurate and comprehensive collection of transcripts originating from a given organism. It holds the key to precise transcript quantification and downstream analysis of differential expressions and regulations. Currently, transcriptome annotations for most crop plants are far from complete. For example, Oryza sativa indica (O. sativa indica) is reported to have 40,759 transcripts in the Ensembl database without alternative transcript isoforms and alternative splicing (AS) events. To generate a high-quality RTD, we conducted RNA sequencing of rice leaf samples collected at various time points during Rhizoctonia solani infection. The obtained reads were analyzed by adopting the recently developed computational analysis pipeline to assemble the RTD with increased transcript and AS diversity for O. sativa indica (IndicaRTD). After stringent quality filtering, the newly constructed transcriptome annotation was comprised of 122,968 non-redundant transcripts from 53,695 genes. This study identified many novel transcripts compared to Ensembl deposited data that are important for regulating molecular and physiological processes in the plant system. Currently, the assembled IndicaRTD must allow fast quantification of transcript and gene expression with high precision.

Detecting signatures of selection on gene expression

A substantial amount of phenotypic diversity results from changes in gene expression levels and patterns. Understanding how the transcriptome evolves is therefore a key priority in identifying mechanisms of adaptive change. However, in contrast to powerful models of sequence evolution, we lack a consensus model of gene expression evolution. Furthermore, recent work has shown that many of the comparative approaches used to study gene expression are subject to biases that can lead to false signatures of selection. Here we first outline the main approaches for describing expression evolution and their inherent biases. Next, we bridge the gap between the fields of phylogenetic comparative methods and transcriptomics to reinforce the main pitfalls of inferring selection on expression patterns and use simulation studies to show that shifts in tissue composition can heavily bias inferences of selection. We close by highlighting the multi-dimensional nature of transcriptional variation and identifying major unanswered questions in disentangling how selection acts on the transcriptome.

Acute thermal stress elicits interactions between gene expression and alternative splicing in a fish of conservation concern

Transcriptomic research provides a mechanistic understanding of an organism's response to environmental challenges such as increasing temperatures, which can provide key insights into the threats posed by thermal challenges associated with urbanization and climate change. Differential gene expression and alternative splicing are two elements of the transcriptomic stress response that may work in tandem, but relatively few studies have investigated these interactions in fishes of conservation concern. We studied the imperilled redside dace (Clinostomus elongatus) as thermal stress is hypothesised to be an important cause of population declines. We tested the hypothesis that gene expression-splicing interactions contribute to the thermal stress response. Wild fish exposed to acute thermal stress were compared with both handling controls and fish sampled directly from a river. Liver tissue was sampled to study the transcriptomic stress response. With a gene set enrichment analysis, we found that thermally stressed fish showed a transcriptional response related to transcription regulation and responses to unfolded proteins, and alternatively spliced genes related to gene expression regulation and metabolism. One splicing factor, prpf38b, was upregulated in the thermally stressed group compared to the other treatments. This splicing factor may have a role in the Jun/AP-1 cellular stress response, a pathway with wide-ranging and context-dependent effects. Given large gene interaction networks and the context-dependent nature of transcriptional responses, our results highlight the importance of understanding interactions between gene expression and splicing for understanding transcriptomic responses to thermal stress. Our results also reveal transcriptional pathways that can inform conservation breeding, translocation, and reintroduction programs for redside dace and other imperilled species by identifying appropriate source populations.