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

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.

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.

The evolution of mammalian Rem2: unraveling the impact of purifying selection and coevolution on protein function, and implications for human disorders

Rad And Gem-Like GTP-Binding Protein 2 (Rem2), a member of the RGK family of Ras-like GTPases, is implicated in Huntington's disease and Long QT Syndrome and is highly expressed in the brain and endocrine cells. We examine the evolutionary history of Rem2 identified in various mammalian species, focusing on the role of purifying selection and coevolution in shaping its sequence and protein structural constraints. Our analysis of Rem2 sequences across 175 mammalian species found evidence for strong purifying selection in 70% of non-invariant codon sites which is characteristic of essential proteins that play critical roles in biological processes and is consistent with Rem2's role in the regulation of neuronal development and function. We inferred epistatic effects in 50 pairs of codon sites in Rem2, some of which are predicted to have deleterious effects on human health. Additionally, we reconstructed the ancestral evolutionary history of mammalian Rem2 using protein structure prediction of extinct and extant sequences which revealed the dynamics of how substitutions that change the gene sequence of Rem2 can impact protein structure in variable regions while maintaining core functional mechanisms. By understanding the selective pressures, protein- and gene - interactions that have shaped the sequence and structure of the Rem2 protein, we gain a stronger understanding of its biological and functional constraints.

Alternative splicing plays a nonredundant role in greater amberjack (Seriola dumerili) in acclimation to ambient salinity fluctuations

Alternative splicing (AS) is an important post-transcriptional mechanism for adaptation of fish to environmental stress. Here, we performed a genome-wide investigation to AS dynamics in greater amberjack (Seriola dumerili), an economical marine teleost, in response to hypo- (10 ppt) and hyper-salinity (40 ppt) stresses. Totally, 2267-2611 differentially spliced events were identified in gills and kidney upon the exposure to undesired salinity regimes. In gills, genes involved in energy metabolism, stimulus response and epithelial cell differentiation were differentially spliced in response to salinity variation, while sodium ion transport and cellular amide metabolism were enhanced in kidney to combat the adverse impacts of salinity changes. Most of these differentially spliced genes were not differentially expressed, and AS was found to regulate different biological processes from differential gene expression, indicative of the functionally nonredundant role of AS in modulating salinity acclimation in greater amberjack. Together, our study highlights the important contribution of post-transcriptional mechanisms to the adaptation of fish to ambient salinity fluctuations and provides theoretical guidance for the conservation of marine fishery resources against increasingly environmental challenges.

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.

Alternative splicing and environmental adaptation in wild house mice

A major goal of evolutionary genetics is to understand the genetic and molecular mechanisms underlying adaptation. Previous work has established that changes in gene regulation may contribute to adaptive evolution, but most studies have focused on mRNA abundance and only a few studies have investigated the role of post-transcriptional processing. Here, we use a combination of exome sequences and short-read RNA-Seq data from wild house mice (Mus musculus domesticus) collected along a latitudinal transect in eastern North America to identify candidate genes for local adaptation through alternative splicing. First, we identified alternatively spliced transcripts that differ in frequency between mice from the northern-most and southern-most populations in this transect. We then identified the subset of these transcripts that exhibit clinal patterns of variation among all populations in the transect. Finally, we conducted association studies to identify cis-acting splicing quantitative trait loci (cis-sQTL), and we identified cis-sQTL that overlapped with previously ascertained targets of selection from genome scans. Together, these analyses identified a small set of alternatively spliced transcripts that may underlie environmental adaptation in house mice. Many of these genes have known phenotypes associated with body size, a trait that varies clinally in these populations. We observed no overlap between these genes and genes previously identified by changes in mRNA abundance, indicating that alternative splicing and changes in mRNA abundance may provide separate molecular mechanisms of adaptation.

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.

Exploring transcriptional and post-transcriptional epigenetic regulation of crf and 11βhsd2 in rainbow trout brain during chronic social stress

Using dominance hierarchies in juvenile rainbow trout (Oncorhynchus mykiss) as a model of chronic social stress in fish, we explored whether epigenetic transcriptional and post-transcriptional mechanisms are involved in the gene expression of corticotropin-releasing factor (crf) and 11β-hydroxysteroid dehydrogenase (11βhsd2), key factors involved in the regulation of the endocrine stress axis response. In juvenile rainbow trout pairs, subordinate individuals display sustained elevation of circulating cortisol concentrations. Cortisol production is controlled by the hypothalamic-pituitary-interrenal (HPI) axis in fish and initiated by CRF release from the preoptic area (POA). Given that crf is modulated during chronic social stress, and that such stress has been implicated in the epigenetic regulation of crf in other taxa, we probed a role for epigenetic regulation of crf transcript abundance in chronically stressed rainbow trout. We also investigated the regulation of the cortisol-metabolising enzyme 11βhsd2 in the POA, which is upregulated in subordinates. The potential involvement of DNA methylation and microRNAs (miRNAs) in the regulation of crf transcript abundance was investigated during social stress in the POA of fish, as was the potential involvement of miRNAs in 11βhsd2 regulation. Although transcript abundances of crf were elevated in subordinate fish after 4 days, DNA methylation profiles within putative promoter sequences upstream of the crf gene were not significantly affected by chronic stress. An inverse relationship between crf and its predicted posttranscriptional regulator miR-103a-3p in the POA suggests that miRNAs may be involved in mediating the effects of chronic social stress on key components of the endocrine stress axis.

Using individual-based trait frequency distributions to forecast plant-pollinator network responses to environmental change

Determining how and why organisms interact is fundamental to understanding ecosystem responses to future environmental change. To assess the impact on plant-pollinator interactions, recent studies have examined how the effects of environmental change on individual interactions accumulate to generate species-level responses. Here, we review recent developments in using plant-pollinator networks of interacting individuals along with their functional traits, where individuals are nested within species nodes. We highlight how these individual-level, trait-based networks connect intraspecific trait variation (as frequency distributions of multiple traits) with dynamic responses within plant-pollinator communities. This approach can better explain interaction plasticity, and changes to interaction probabilities and network structure over spatiotemporal or other environmental gradients. We argue that only through appreciating such trait-based interaction plasticity can we accurately forecast the potential vulnerability of interactions to future environmental change. We follow this with general guidance on how future studies can collect and analyse high-resolution interaction and trait data, with the hope of improving predictions of future plant-pollinator network responses for targeted and effective conservation.

Alternative splicing of a carboxyl/choline esterase gene enhances the fenpropathrin tolerance of Tetranychus cinnabarinus

Detoxification plays a crucial role in agricultural pests to withstand pesticides, and cytochrome P450s, carboxyl/choline esterases (CCEs), and glutathione-S-transferases are the main proteins responsible for their detoxification ability. The activity of CCEs can be upregulated, downregulated, or modified by mutation. However, few studies have examined the role of alternative splicing in altering the properties of CCEs. We identified 2 variants of TcCCE23 in Tetranychus cinnabarinus: a long version (CCE23-V1) and a short version that is 18 nucleotides shorter than CCE23-V1 (CCE23-V2). Whether splicing affects the activity of TcCCE23 remains unclear. Overexpression of CCE23-V2 in fenpropathrin-resistant T. cinnabarinus revealed that splicing affected the detoxification of fenpropathrin by CCE23-V2. The mortality of mites was significantly higher when the expression of CCE23-V2 was knocked down (43.2% ± 3.3%) via injection of CCE23-dsRNA (double-stranded RNA) compared with the control group injected with green fluorescent protein-dsRNA under fenpropathrin exposure; however, the downregulation of CCE23-V1 (61.3% ± 6.3%) by CCE23-small interfering RNA had no such effect, indicating CCE23-V2 plays a greater role in xenobiotic metabolism than CCE23-V1. The tolerance of flies overexpressing CCE23-V2 to fenpropathrin (50% lethal dose [LD50 ] = 19.47 μg/g) was significantly higher than that of Gal4/UAS-CCE23-V1 transgenic flies (LD50  = 13.11 μg/g). Molecular docking analysis showed that splicing opened a "gate" that enlarges the substrate binding cavity of CCE23-V2, might enhance the ability of CCE23-V2 to harbor fenpropathrin molecules. These findings suggest that splicing might enhance the detoxifying capability of TcCCE23. Generally, our data improve the understanding of the diversity and complexity of the mechanisms underlying the regulation of CCEs.

Sequence variations affect the 5' splice site selection of plant introns

Introns are noncoding sequences spliced out of pre-mRNAs by the spliceosome to produce mature mRNAs. The 5' ends of introns mostly begin with GU and have a conserved sequence motif of AG/GUAAGU that could base-pair with the core sequence of U1 snRNA of the spliceosome. Intriguingly, ∼ 1% of introns in various eukaryotic species begin with GC. This occurrence could cause misannotation of genes; however, the underlying splicing mechanism is unclear. We analyzed the sequences around the intron 5' splice site (ss) in Arabidopsis (Arabidopsis thaliana) and found sequences at the GC intron ss are much more stringent than those of GT introns. Mutational analysis at various positions of the intron 5' ss revealed that although mutations impair base pairing, different mutations at the same site can have different effects, suggesting that steric hindrance also affects splicing. Moreover, mutations of 5' ss often activate a hidden ss nearby. Our data suggest that the 5' ss is selected via a competition between the major ss and the nearby minor ss. This work not only provides insights into the splicing mechanism of intron 5' ss but also improves the accuracy of gene annotation and the study of the evolution of intron 5' ss.

Comparative transcriptomes reveal geographic differences in the ability of the liver of plateau zokors (Eospalax baileyi) to respond and adapt to toxic plants

BACKGROUND

Environmental changes are expected to intensify in the future. The invasion of toxic plants under environmental changes may change herbivore feeding environments. Herbivores living long-term in toxic plant-feeding environments will inevitably ingest plant secondary metabolites (PSMs), and under different feeding environments are likely to have unique protection mechanisms that support improved adaptation to PSMs in their habitat. We aimed to compare different subterranean herbivore population responses and adaptations to toxic plants to unveil their feeding challenges.

RESULTS

Here, we investigated the adaptive capacity of the liver in two geographically separated populations of plateau zokors (Eospalax baileyi) before and after exposure to the toxic plant Stellera chamaejasme (SC), at the organ, biochemical, and transcriptomic levels. The results showed no significant liver granules or inflammatory reactions in the Tianzhu (TZ) population after the SC treatment. The transaminase level in the TZ population was significantly lower than that in the Luqu population. Transcriptome analysis revealed that the TZ population exhibited interactions with other detoxification metabolic pathways by oxytocin pathway-associated genes, including diacylglycerol lipase alpha (Dagla), calcium/calmodulin dependent protein kinase II Alpha (Camk2a), and CD38 molecule (Cd38). The phase II process of liver drug metabolism increased to promote the rate of metabolism. We found that alternative splicing (AS) and the expression of the cyclin D (Ccnd1) gene interact-a TZ population hallmark-reduced liver inflammatory responses.

CONCLUSION

Our study supports the detoxification limitation hypothesis that differences in liver detoxification metabolism gene expression and AS are potential factors in herbivore adaptation to PSMs and may be a strategy of different herbivore populations to improve toxic plant adaptability.

Evolutionary Shortcuts via Multinucleotide Substitutions and Their Impact on Natural Selection Analyses

Inference and interpretation of evolutionary processes, in particular of the types and targets of natural selection affecting coding sequences, are critically influenced by the assumptions built into statistical models and tests. If certain aspects of the substitution process (even when they are not of direct interest) are presumed absent or are modeled with too crude of a simplification, estimates of key model parameters can become biased, often systematically, and lead to poor statistical performance. Previous work established that failing to accommodate multinucleotide (or multihit, MH) substitutions strongly biases dN/dS-based inference towards false-positive inferences of diversifying episodic selection, as does failing to model variation in the rate of synonymous substitution (SRV) among sites. Here, we develop an integrated analytical framework and software tools to simultaneously incorporate these sources of evolutionary complexity into selection analyses. We found that both MH and SRV are ubiquitous in empirical alignments, and incorporating them has a strong effect on whether or not positive selection is detected (1.4-fold reduction) and on the distributions of inferred evolutionary rates. With simulation studies, we show that this effect is not attributable to reduced statistical power caused by using a more complex model. After a detailed examination of 21 benchmark alignments and a new high-resolution analysis showing which parts of the alignment provide support for positive selection, we show that MH substitutions occurring along shorter branches in the tree explain a significant fraction of discrepant results in selection detection. Our results add to the growing body of literature which examines decades-old modeling assumptions (including MH) and finds them to be problematic for comparative genomic data analysis. Because multinucleotide substitutions have a significant impact on natural selection detection even at the level of an entire gene, we recommend that selection analyses of this type consider their inclusion as a matter of routine. To facilitate this procedure, we developed, implemented, and benchmarked a simple and well-performing model testing selection detection framework able to screen an alignment for positive selection with two biologically important confounding processes: site-to-site synonymous rate variation, and multinucleotide instantaneous substitutions.

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.