Rapid sex-specific adaptation to high temperature in Drosophila

Reduced Parallel Gene Expression Evolution With Increasing Genetic Divergence-A Hallmark of Polygenic Adaptation

Parallel evolution, the repeated evolution of similar traits in independent lineages, is a topic of considerable interest in evolutionary biology. Although previous studies have focused on the parallelism of phenotypic traits and their underlying genetic basis, the extent of parallelism at the level of gene expression across different levels of genetic divergence is not yet fully understood. This study investigates the evolution of gene expression in replicate Drosophila populations exposed to the same novel environment at three divergence levels: within a population, between populations and between species. We show that adaptive gene expression changes are more heterogeneous with increasing genetic divergence between the compared groups. This finding suggests that the adaptive architecture-comprising factors such as allele frequencies and the effect size of contributing loci-becomes more distinct with increasing divergence. As a result, this leads to a reduction in parallel gene expression evolution. This result implies that redundancy is a crucial factor in both genetic selection responses and gene expression evolution. Hence, our findings are consistent with the omnigenic model, which posits that selection acts on higher-order phenotypes. This work contributes to our understanding of phenotypic evolution and the complex interplay between genomic and molecular responses.

Pleiotropy increases parallel selection signatures during adaptation from standing genetic variation

The phenomenon of parallel evolution, whereby similar genomic and phenotypic changes occur across replicated pairs of populations or species, is widely studied. Nevertheless, the determining factors of parallel evolution remain poorly understood. Theoretical studies have proposed that pleiotropy, the influence of a single gene on multiple traits, is an important factor. In order to gain a deeper insight into the role of pleiotropy for parallel evolution from standing genetic variation, we characterized the interplay between parallelism, polymorphism, and pleiotropy. The present study examined the parallel gene expression evolution in 10 replicated populations of Drosophila simulans, which were adapted from standing variation to the same new temperature regime. The data demonstrate that the parallel evolution of gene expression from standing genetic variation is positively correlated with the strength of pleiotropic effects. The ancestral variation in gene expression is, however, negatively correlated with parallelism. Given that pleiotropy is also negatively correlated with gene expression variation, we conducted a causal analysis to distinguish cause and correlation and evaluate the role of pleiotropy. The causal analysis indicated that both direct (causative) and indirect (correlational) effects of pleiotropy contribute to parallel evolution. The indirect effect is mediated by historic selective constraint in response to pleiotropy. This results in parallel selection responses due to the reduced standing variation of pleiotropic genes. The direct effect of pleiotropy is likely to reflect a genetic correlation among adaptive traits, which in turn gives rise to synergistic effects and higher parallelism.

From whole bodies to single cells: A guide to transcriptomic approaches for ecology and evolutionary biology

RNA sequencing (RNAseq) methodology has experienced a burst of technological developments in the last decade, which has opened up opportunities for studying the mechanisms of adaptation to environmental factors at both the organismal and cellular level. Selecting the most suitable experimental approach for specific research questions and model systems can, however, be a challenge and researchers in ecology and evolution are commonly faced with the choice of whether to study gene expression variation in whole bodies, specific tissues, and/or single cells. A wide range of sometimes polarised opinions exists over which approach is best. Here, we highlight the advantages and disadvantages of each of these approaches to provide a guide to help researchers make informed decisions and maximise the power of their study. Using illustrative examples of various ecological and evolutionary research questions, we guide the readers through the different RNAseq approaches and help them identify the most suitable design for their own projects.

Reproductive isolation arises during laboratory adaptation to a novel hot environment

BACKGROUND

Reproductive isolation can result from adaptive processes (e.g., ecological speciation and mutation-order speciation) or stochastic processes such as "system drift" model. Ecological speciation predicts barriers to gene flow between populations from different environments, but not among replicate populations from the same environment. In contrast, reproductive isolation among populations independently adapted to the same/similar environment can arise from both mutation-order speciation or system drift.

RESULTS

In experimentally evolved populations adapting to a hot environment for over 100 generations, we find evidence for pre- and postmating reproductive isolation. On one hand, an altered lipid metabolism and cuticular hydrocarbon composition pointed to possible premating barriers between the ancestral and replicate evolved populations. On the other hand, the pronounced gene expression differences in male reproductive genes may underlie the postmating isolation among replicate evolved populations adapting to the same environment with the same standing genetic variation.

CONCLUSION

Our study confirms that replicated evolution experiments provide valuable insights into the mechanisms of speciation. The rapid emergence of the premating reproductive isolation during temperature adaptation showcases incipient ecological speciation. The potential evidence of postmating reproductive isolation among replicates gave rise to two hypotheses: (1) mutation-order speciation through a common selection on early fecundity leading to an inherent inter-locus sexual conflict; (2) system drift with genetic drift along the neutral ridges.

The evolution of sexual dimorphism in gene expression in response to a manipulation of mate competition

Many genes are differentially expressed between males and females and patterns of sex-biased gene expression (SBGE) vary among species. Some of this variation is thought to have evolved in response to differences in mate competition among species that cause varying patterns of sex-specific selection. We used experimental evolution to test this by quantifying SBGE and sex-specific splicing in 15 Drosophila melanogaster populations that evolved for 104 generations in mating treatments that removed mate competition via enforced monogamy, or allowed mate competition in either small, simple, or larger, structurally more complex mating environments. Consistent with sex-specific selection affecting SBGE, initially sex-biased genes diverged in expression more among treatments than unbiased genes, and there was greater expression divergence for male- than female-biased genes. It has been suggested the transcriptome should be "feminized" under monogamy because of the removal of sexual selection on males; we did not observe this, likely because selection differs in additional ways between monogamy vs. polygamy. Significant divergence in average expression dimorphism between treatments was observed and, in some treatment comparisons, the direction of the divergence differed across different sex-bias categories. There was not a generalized reduction in expression dimorphism under enforced monogamy.

Insecticide resistant mosquitoes remain thermal stress resistant, without loss of thermal plasticity

A major component of mosquito's climate change response is their heat tolerance, and any ability to rapidly adjust to extreme environmental conditions through phenotypic plasticity. The excessive use of insecticides for the control of major mosquito species leads to resistant populations, however it is largely unclear if this concurrently impacts thermal stress resistance and their potential to adjust tolerance via phenotypic plasticity. Culex pipiens pipiens, Culex pipiens molestus and Aedes albopictus populations obtained from the same region were subjected for 12 generations to selection trials to larvicides Diflubenzuron (DFB) and Bacillus thuringiensis subsp. israelensis (Bti) to develop insecticide resistance. Adults emerging from the selected populations were acclimated at different temperatures and the upper and lower critical thermal limits (CTmax and CTmin) were estimated using dynamic thermal assays. In addition, the supercooling points (SCPs) of non-acclimated adults of resistant and control populations were determined. Our results revealed marked differences in thermal response among the three species, the different acclimation regimes and sexes. Aedes albopictus was more resistant in high than low temperatures compared to both Culex pipiens biotypes. Culex forms responded similarly to heat but differently to cold stress. In both forms, females responded better than males to all thermal stressors. Acclimation at higher and lower temperatures improves CTmax and CTmin values, respectively in both insecticide resistant and control populations of all three species. Overall, selection to insecticides did not affect the thermal performance of adults. Hence, insecticide-resistant mosquito populations perform similarly to untreated ones and are capable of readily adapting to new environmental changes rising concerns regarding their geographic range expansion and disease transmission globally.

Positive selection and relaxed purifying selection contribute to rapid evolution of male-biased genes in a dioecious flowering plant

Sex-biased genes offer insights into the evolution of sexual dimorphism. Sex-biased genes, especially those with male bias, show elevated evolutionary rates of protein sequences driven by positive selection and relaxed purifying selection in animals. Although rapid sequence evolution of sex-biased genes and evolutionary forces have been investigated in animals and brown algae, less is known about evolutionary forces in dioecious angiosperms. In this study, we separately compared the expression of sex-biased genes between female and male floral buds and between female and male flowers at anthesis in dioecious Trichosanthes pilosa (Cucurbitaceae). In floral buds, sex-biased gene expression was pervasive, and had significantly different roles in sexual dimorphism such as physiology. We observed higher rates of sequence evolution for male-biased genes in floral buds compared to female-biased and unbiased genes. Male-biased genes under positive selection were mainly associated with functions to abiotic stress and immune responses, suggesting that high evolutionary rates are driven by adaptive evolution. Additionally, relaxed purifying selection may contribute to accelerated evolution in male-biased genes generated by gene duplication. Our findings, for the first time in angiosperms, suggest evident rapid evolution of male-biased genes, advance our understanding of the patterns and forces driving the evolution of sexual dimorphism in dioecious plants.

Impact of Heat Stress on Transposable Element Expression and Derived Small RNAs in Drosophila subobscura

Global warming is forcing insect populations to move and adapt, triggering adaptive genetic responses. Thermal stress is known to alter gene expression, repressing the transcription of active genes, and inducing others, such as those encoding heat shock proteins. It has also been related to the activation of some specific transposable element (TE) families. However, the actual magnitude of this stress on the whole genome and the factors involved in these genomic changes are still unclear. We studied mRNAs and small RNAs in gonads of two Drosophila subobscura populations, considered a good model to study adaptation to temperature changes. In control conditions, we found that a few genes and TE families were differentially expressed between populations, pointing out their putative involvement in the adaptation of populations to their different environments. Under heat stress, sex-specific changes in gene expression together with a trend toward overexpression, mainly of heat shock response-related genes, were observed. We did not observe large changes of TE expression nor small RNA production due to stress. Only population and sex-specific expression changes of some TE families (mainly retrotransposons), or the amounts of siRNAs and piRNAs, derived from specific TE families were observed, as well as the piRNA production from some piRNA clusters. Changes in small RNA amounts and TE expression could not be clearly correlated, indicating that other factors as chromatin modulation could also be involved. This work provides the first whole transcriptomic study including genes, TEs, and small RNAs after a heat stress in D. subobscura.

How predictable is adaptation from standing genetic variation? Experimental evolution in Drosophila highlights the central role of redundancy and linkage disequilibrium

Experimental evolution is well-suited to test the predictability of evolution without the confounding effects of inaccurate forecasts about future environments. Most of the literature about parallel (and thus predictable) evolution has been carried out in asexual microorganisms, which adapt by de novo mutations. Nevertheless, parallel evolution has also been studied in sexual species at the genomic level. Here, I review the evidence for parallel evolution in Drosophila, the best-studied obligatory outcrossing model for adaptation from standing genetic variation in the laboratory. Similar to asexual microorganisms, evidence for parallel evolution varies between the focal hierarchical levels. Selected phenotypes consistently respond in a very predicable way, but the underlying allele frequency changes are much less predictable. The most important insight is that the predictability of the genomic selection response for polygenic traits depends highly on the founder population and to a much lesser extent on the selection regime. This implies that predicting adaptive genomic response is challenging and requires a good understanding of the adaptive architecture (including linkage disequilibrium) in the ancestral populations. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

Thermal biology of invasive Aedes mosquitoes in the context of climate change

The increasing incidence of arboviral diseases in tropical endemic areas and their emergence in new temperate countries is one of the most important challenges that Public Health agencies are currently facing. Because mosquitoes are poikilotherms, shifts in temperature influence physiological functions besides egg viability. These traits impact not only vector density, but also their interaction with their hosts and arboviruses. As such the relationship among mosquitoes, arboviral diseases and temperature is complex. Here, we summarize current knowledge on the thermal biology of Aedes invasive mosquitoes, highlighting differences among species. We also emphasize the need to expand knowledge on the variability in thermal sensitivity across populations within a species, especially in light of climate change that encompasses increase not only in mean environmental temperature but also in the frequency of hot and cold snaps. Finally, we suggest a novel experimental approach to investigate the molecular architecture of thermal adaptation in mosquitoes.

Natural variation in Drosophila shows weak pleiotropic effects

Background

Pleiotropy describes the phenomenon in which a gene affects multiple phenotypes. The extent of pleiotropy is still disputed, mainly because of issues of inadequate power of analyses. A further challenge is that empirical tests of pleiotropy are restricted to a small subset of all possible phenotypes. To overcome these limitations, we propose a new measurement of pleiotropy that integrates across many phenotypes and multiple generations to improve power.

Results

We infer pleiotropy from the fitness cost imposed by frequency changes of pleiotropic loci. Mixing Drosophila simulans populations, which adapted independently to the same new environment using different sets of genes, we show that the adaptive frequency changes have been accompanied by measurable fitness costs.

Conclusions

Unlike previous studies characterizing the molecular basis of pleiotropy, we show that many loci, each of weak effect, contribute to genome-wide pleiotropy. We propose that the costs of pleiotropy are reduced by the modular architecture of gene expression, which facilitates adaptive gene expression changes with low impact on other functions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02680-4.

Evolution of phenotypic variance in response to a novel hot environment

Shifts in trait means are widely considered as evidence for adaptive responses, but the impact on phenotypic variance remains largely unexplored. Classic quantitative genetics provides a theoretical framework to predict how selection on phenotypic mean affects the variance. In addition to this indirect effect, it is also possible that the variance of the trait is the direct target of selection, but experimentally characterized cases are rare. Here, we studied gene expression variance of Drosophila simulans males before and after 100 generations of adaptation to a novel hot laboratory environment. In each of the two independently evolved populations, the variance of 125 and 97 genes was significantly reduced. We propose that the drastic loss in environmental complexity from nature to the laboratory may have triggered selection for reduced variance. Our observation that selection could drive changes in the variance of gene expression could have important implications for studies of adaptation processes in natural and experimental populations.

Sexual conflict in a changing environment

Sexual conflict has extremely important consequences for various evolutionary processes including its effect on local adaptation and extinction probability during environmental change. The awareness that the intensity and dynamics of sexual conflict is highly dependent on the ecological setting of a population has grown in recent years, but much work is yet to be done. Here, we review progress in our understanding of the ecology of sexual conflict and how the environmental sensitivity of such conflict feeds back into population adaptivity and demography, which, in turn, determine a population's chances of surviving a sudden environmental change. We link two possible forms of sexual conflict - intralocus and interlocus sexual conflict - in an environmental context and identify major gaps in our knowledge. These include sexual conflict responses to fluctuating and oscillating environmental changes and its influence on the interplay between interlocus and intralocus sexual conflict, among others. We also highlight the need to move our investigations into more natural settings and to investigate sexual conflict dynamics in wild populations.

Neuroendocrine Regulation of Reproductive Dormancy in the Fruit Fly Drosophila melanogaster: A Review of Juvenile Hormone-Dependent Regulation

Animals can adjust their physiology, helping them survive and reproduce under a wide range of environmental conditions. One of the strategies to endure unfavorable environmental conditions such as low temperature and limited food supplies is dormancy. In some insect species, this may manifest as reproductive dormancy, which causes their reproductive organs to be severely depleted under conditions unsuitable for reproduction. Reproductive dormancy in insects is induced by a reduction in juvenile hormones synthesized in the corpus allatum (pl. corpora allata; CA) in response to winter-specific environmental cues, such as low temperatures and short-day length. In recent years, significant progress has been made in the study of dormancy-inducing conditions dependent on CA control mechanisms in Drosophila melanogaster. This review summarizes dormancy control mechanisms in D. melanogaster and discusses the implications for future studies of insect dormancy, particularly focusing on juvenile hormone-dependent regulation.

Evolution of sex-specific heat stress tolerance and larval Hsp70 expression in populations of Drosophila melanogaster adapted to larval crowding

The ability to tolerate temperature stress is an important component of adult fitness. In holometabolous insects like Drosophila melanogaster, adult stress resistance can be affected by growth conditions experienced during the larval stages. Although evolution under crowded larval conditions is known to lead to the correlated evolution of many adult traits, its consequences on adult heat stress tolerance have not been investigated. Therefore, in the present study, we assessed the adult heat stress tolerance in populations of D. melanogaster adapted to a stressful larval crowding environment. We used replicate populations of D. melanogaster, selected for adaptation to larval crowding stress (MCUs), for more than 230 generations, and their respective controls (MBs). Larvae from selected and control populations were grown under crowded and uncrowded conditions, and their adult heat shock resistance at two different temperatures was measured. Further, we compared Hsp70 expression in crowded and uncrowded larvae of both populations and also measured the Hsp70 expression after a mild heat treatment in adults of selected and control populations. Our results showed that adaptation to larval crowding leads to the evolution of Hsp70 gene expression in larval stages and improves adult heat stress tolerance ability in males, but not in females.

Evolutionary Responses to Warming

Climate change is predicted to dramatically alter biological diversity and distributions, driving extirpations, extinctions, and extensive range shifts across the globe. Warming can also, however, lead to phenotypic or behavioural plasticity, as species adapt to new conditions. Recent genomic research indicates that some species are capable of rapid evolution as selection favours adaptive responses to environmental change and altered or novel niche spaces. New advances are providing mechanistic insights into how temperature might accelerate evolution in the Anthropocene. These discoveries highlight intriguing new research directions - such as using geothermal and polar systems combined with powerful genomic tools - that will help us to understand the processes underpinning adaptive evolution and better project how ecosystems will change in a warming world.

Sexual Dimorphism through the Lens of Genome Manipulation, Forward Genetics, and Spatiotemporal Sequencing

Sexual reproduction often leads to selection that favors the evolution of sex-limited traits or sex-specific variation for shared traits. These sexual dimorphisms manifest due to sex-specific genetic architectures and sex-biased gene expression across development, yet the molecular mechanisms underlying these patterns are largely unknown. The first step is to understand how sexual dimorphisms arise across the genotype-phenotype-fitness map. The emergence of "4D genome technologies" allows for efficient, high-throughput, and cost-effective manipulation and observations of this process. Studies of sexual dimorphism will benefit from combining these technological advances (e.g., precision genome editing, inducible transgenic systems, and single-cell RNA sequencing) with clever experiments inspired by classic designs (e.g., bulked segregant analysis, experimental evolution, and pedigree tracing). This perspective poses a synthetic view of how manipulative approaches coupled with cutting-edge observational methods and evolutionary theory are poised to uncover the molecular genetic basis of sexual dimorphism with unprecedented resolution. We outline hypothesis-driven experimental paradigms for identifying genetic mechanisms of sexual dimorphism among tissues, across development, and over evolutionary time.

Parallel gene expression evolution in natural and laboratory evolved populations

Ecological adaptation is frequently inferred by the comparison of natural populations from different environments. Nevertheless, inference of the selective forces suffers the challenge that many environmental factors covary. With well-controlled environmental conditions, experimental evolution provides a powerful approach to complement the analysis of natural populations. On the other hand, it is apparent that laboratory conditions differ in many ways from natural environments, which raises the question as to what extent selection responses in experimental evolution studies can inform us about adaptation processes in the wild. In this study, we compared the expression profiles of replicated Drosophila melanogaster populations which have been exposed to two distinct temperature regimes (18/28 and 10/20°C) in the laboratory for more than 80 generations. Using gene-wise differential expression analysis and co-expression network analysis, we identified 541 genes and three coregulated gene modules that evolved in the same direction in both temperature regimes, and most of these changes probably reflect an adaptation to the space constraint or diurnal temperature fluctuation that is common in both selection regimes. In total, 203 genes and seven modules evolved temperature-specific expression changes. Remarkably, we detected a significant overlap of these temperature-adaptive genes/modules from experimental evolution with temperature-adaptive genes inferred from natural Drosophila populations covering two different temperature clines. We conclude that well-designed experimental evolution studies are a powerful tool to dissect evolutionary responses.

The Evolution of Phenotypic Plasticity in Response to Temperature Stress

Phenotypic plasticity is the ability of a single genotype to produce different phenotypes in response to environmental variation. The importance of phenotypic plasticity in natural populations and its contribution to phenotypic evolution during rapid environmental change is widely debated. Here, we show that thermal plasticity of gene expression in natural populations is a key component of its adaptation: evolution to novel thermal environments increases ancestral plasticity rather than mean genetic expression. We determined the evolution of plasticity in gene expression by conducting laboratory natural selection on a Drosophila simulans population in hot and cold environments. After more than 60 generations in the hot environment, 325 genes evolved a change in plasticity relative to the natural ancestral population. Plasticity increased in 75% of these genes, which were strongly enriched for several well-defined functional categories (e.g., chitin metabolism, glycolysis, and oxidative phosphorylation). Furthermore, we show that plasticity in gene expression of populations exposed to different temperatures is rather similar across species. We conclude that most of the ancestral plasticity can evolve further in more extreme environments.

Unique genetic signatures of local adaptation over space and time for diapause, an ecologically relevant complex trait, in Drosophila melanogaster

Organisms living in seasonally variable environments utilize cues such as light and temperature to induce plastic responses, enabling them to exploit favorable seasons and avoid unfavorable ones. Local adapation can result in variation in seasonal responses, but the genetic basis and evolutionary history of this variation remains elusive. Many insects, including Drosophila melanogaster, are able to undergo an arrest of reproductive development (diapause) in response to unfavorable conditions. In D. melanogaster, the ability to diapause is more common in high latitude populations, where flies endure harsher winters, and in the spring, reflecting differential survivorship of overwintering populations. Using a novel hybrid swarm-based genome wide association study, we examined the genetic basis and evolutionary history of ovarian diapause. We exposed outbred females to different temperatures and day lengths, characterized ovarian development for over 2800 flies, and reconstructed their complete, phased genomes. We found that diapause, scored at two different developmental cutoffs, has modest heritability, and we identified hundreds of SNPs associated with each of the two phenotypes. Alleles associated with one of the diapause phenotypes tend to be more common at higher latitudes, but these alleles do not show predictable seasonal variation. The collective signal of many small-effect, clinally varying SNPs can plausibly explain latitudinal variation in diapause seen in North America. Alleles associated with diapause are segregating in Zambia, suggesting that variation in diapause relies on ancestral polymorphisms, and both pro- and anti-diapause alleles have experienced selection in North America. Finally, we utilized outdoor mesocosms to track diapause under natural conditions. We found that hybrid swarms reared outdoors evolved increased propensity for diapause in late fall, whereas indoor control populations experienced no such change. Our results indicate that diapause is a complex, quantitative trait with different evolutionary patterns across time and space.

Early response to heat stress in Chinese tongue sole (Cynoglossus semilaevis): performance of different sexes, candidate genes and networks

BACKGROUND

Temperature is known to affect living organisms and alter the expression of responsive genes, which affects a series of life processes, such as development, reproduction and metabolism. Several genes and gene families have been involved in high temperature responses, such as heat shock protein (hsp) family, Jumonji family and genes related to cortisol synthesis. Gonad is a vital organ related to the existence of a species. However, the comprehensive understanding of gonadal responses to environmental temperature is limited.

RESULTS

To explore the effects of environmental temperature on genes and gene networks in gonads, we performed acute heat treatment (48 h) on Chinese tongue sole (Cynoglossus semilaevis). Gonadal transcriptome analysis was conducted on females, pseudomales and males exposed to high (28 °C) and normal (22 °C) temperatures. A total of 1226.24 million clean reads were obtained from 18 libraries. Principal component analysis (PCA) and differentially expressed gene (DEG) analysis revealed different performance of sex responses to heat stress. There were 4565, 790 and 1117 specific genes altered their expression level in females, pseudomales and males, respectively. Of these, genes related to hsp gene family, cortisol synthesis and metabolism and epigenetic regulation were involved in early heat response. Furthermore, a total of 1048 DEGs were shared among females, pesudomales and males, which may represent the inherent difference between high and normal temperatures. Genes, such as eef1akmt3, eef1akmt4, pnmt and hsp family members, were found.

CONCLUSIONS

Our results depicted for the first time the gonadal gene expression under acute high temperature treatment in Chinese tongue sole. The findings may provide a clue for understanding the responses of genes and networks to environmental temperature.