Temperature stress mediates decanalization and dominance of gene expression in Drosophila melanogaster

The modifying effect of mutant LRRK2 on mutant GBA1-associated Parkinson disease

Parkinson disease (PD) is the second most common neurodegenerative disease. While most cases are sporadic, in ~ 5%-10% of PD patients the disease is caused by mutations in several genes, among them GBA1 (glucocerebrosidase beta 1) and LRRK2 (leucine-rich repeat kinase 2), both prevalent among the Ashkenazi Jewish population. LRRK2-associated PD tends to be milder than GBA1-associated PD. Several recent clinical studies have suggested that carriers of both GBA1 and LRRK2 mutations develop milder PD compared to that observed among GBA1 carriers. These findings strongly suggested an interplay between the two genes in the development and progression of PD. In the present study Drosophila was employed as a model to investigate the impact of mutations in the LRRK2 gene on mutant GBA1-associated PD. Our results strongly indicated that flies expressing both mutant genes exhibited milder parkinsonian signs compared to the disease developed in flies expressing only a GBA1 mutation. This was corroborated by a decrease in the ER stress response, increase in the number of dopaminergic cells, elevated levels of tyrosine hydroxylase, reduced neuroinflammation, improved locomotion and extended survival. Furthermore, a significant decrease in the steady-state levels of mutant GBA1-encoded GCase was observed in the presence of mutant LRRK2, strongly implying a role for mutant LRRK2 in degradation of mutant GCase.

Intraspecific gene regulation in cis- and trans

Changes in gene expression underlie much of evolution and occur via either cis-acting mutations, which lie near the affected gene and act in a context-specific manner, or trans-acting mutations, which may be far from the affected gene and act through diffusible molecules such as transcription factors. A commonly held view is that most expression variation within species is controlled in trans- while expression differences between species are largely controlled in cis-. Here, we summarize recent intraspecific gene regulation studies and find, contrary to this widely held view, that many studies in diverse taxa have revealed a large role for cis-acting mutations underlying expression variation within species. A review of the existing literature also shows that preparations using whole organisms rather than individual tissues may be biased toward identifying trans-regulation. Moreover, we note several examples of predominantly cis-acting regulation in recently diverged populations adapted to different environments. We highlight the challenges of drawing general conclusions from comparisons among studies that use different methodologies and we offer suggestions for studies that will address outstanding questions concerning the evolution of gene regulation.

Gene-by-environment Interactions and Adaptive Body Size Variation in Mice From the Americas

The relationship between genotype and phenotype is often mediated by the environment. Moreover, gene-by-environment (GxE) interactions can contribute to variation in phenotypes and fitness. In the last 500 yr, house mice have invaded the Americas. Despite their short residence time, there is evidence of rapid climate adaptation, including shifts in body size and aspects of metabolism with latitude. Previous selection scans have identified candidate genes for metabolic adaptation. However, environmental variation in diet as well as GxE interactions likely impact body mass variation in wild populations. Here, we investigated the role of the environment and GxE interactions in shaping adaptive phenotypic variation. Using new locally adapted inbred strains from North and South America, we evaluated response to a high-fat diet, finding that sex, strain, diet, and the interaction between strain and diet contributed significantly to variation in body size. We also found that the transcriptional response to diet is largely strain-specific, indicating that GxE interactions affecting gene expression are pervasive. Next, we used crosses between strains from contrasting climates to characterize gene expression regulatory divergence on a standard diet and on a high-fat diet. We found that gene regulatory divergence is often condition-specific, particularly for trans-acting changes. Finally, we found evidence for lineage-specific selection on cis-regulatory variation involved in diverse processes, including lipid metabolism. Overlap with scans for selection identified candidate genes for environmental adaptation with diet-specific effects. Together, our results underscore the importance of environmental variation and GxE interactions in shaping adaptive variation in complex traits.

Sex-Specific Dominance of Gene Expression in Seed Beetles

When different alleles are favored in different environments, dominance reversal where alternate alleles are dominant in the environment in which they are favored can generate net balancing selection. The sexes represent two distinct genetic environments and sexually antagonistic (SA) selection can maintain genetic variation, especially when the alleles involved show sex-specific dominance. Sexual dimorphism in gene expression is pervasive and has been suggested to result from SA selection. Yet, whether gene-regulatory variation shows sex-specific dominance is poorly understood. We tested for sex-specific dominance in gene expression using three crosses between homozygous lines derived from a population of a seed beetle, where a previous study documented a signal of dominance reversal for fitness between the sexes. Overall, we found that the dominance effects of variants affecting gene expression were positively correlated between the sexes (r = 0.33 to 0.44). Yet, 586 transcripts showed significant differences in dominance between the sexes. Sex-specific dominance was significantly more common in transcripts with more sex-biased expression, in two of three of our crosses. Among transcripts showing sex-specific dominance, lesser sexual dimorphism in gene expression among heterozygotes was somewhat more common than greater. Gene ontology enrichment analyses showed that functional categories associated with known SA phenotypes in Callosobruchus maculatus were overrepresented among transcripts with sex-specific dominance, including genes involved in metabolic processes and the target-of-rapamycin pathway. Our results support the suggestion that sex-specific dominance of regulatory variants contributes to the maintenance of genetic variation in fitness mediated by SA selection in this species.

Cis-regulatory Variation in Relation to Sex and Sexual Dimorphism in Drosophila melanogaster

Much of sexual dimorphism is likely due to sex-biased gene expression, which results from differential regulation of a genome that is largely shared between males and females. Here, we use allele-specific expression to explore cis-regulatory variation in Drosophila melanogaster in relation to sex. We develop a Bayesian framework to infer the transcriptome-wide joint distribution of cis-regulatory effects across the sexes. We also examine patterns of cis-regulatory variation with respect to two other levels of variation in sexual dimorphism: (i) across genes that vary in their degree of sex-biased expression and (ii) among tissues that vary in their degree of dimorphism (e.g. relatively low dimorphism in heads vs. high dimorphism in gonads). We uncover evidence of widespread cis-regulatory variation in all tissues examined, with female-biased genes being especially enriched for this variation. A sizeable proportion of cis-regulatory variation is inferred to have sex-specific effects, with sex-dependent cis effects being much more frequent in gonads than in heads. Finally, we find some genes where 1 allele contributes to more than 50% of a gene's expression in heterozygous males but <50% of its expression in heterozygous females. Such variants could provide a mechanism for sex-specific dominance reversals, a phenomenon important for sexually antagonistic balancing selection. However, tissue differences in allelic imbalance are approximately as frequent as sex differences, perhaps suggesting that sexual conflict may not be particularly unique in shaping patterns of expression variation.

Functional Analysis of Human GBA1 Missense Mutations in Drosophila: Insights into Gaucher Disease Pathogenesis and Phenotypic Consequences

The human GBA1 gene encodes lysosomal acid β-glucocerebrosidase, whose activity is deficient in Gaucher disease (GD). In Drosophila, there are two GBA1 orthologs, Gba1a and Gba1b, and Gba1b is the bona fide GCase encoding gene. Several fly lines with different deletions in the Gba1b were studied in the past. However, since most GD-associated GBA1 mutations are point mutations, we created missense mutations homologous to the two most common GD mutations: the mild N370S mutation (D415S in Drosophila) and the severe L444P mutation (L494P in Drosophila), using the CRISPR-Cas9 technology. Flies homozygous for the D415S mutation (dubbed D370S hereafter) presented low GCase activity and substrate accumulation, which led to lysosomal defects, activation of the Unfolded Protein Response (UPR), inflammation/neuroinflammation, and neurodegeneration along with earlier death compared to control flies. Surprisingly, the L494P (called L444P hereafter) flies presented higher GCase activity with fewer lysosomal defects and milder disease in comparison to that presented by the D370S homozygous flies. Treatment with ambroxol had a limited effect on all homozygous fly lines tested. Overall, our results underscore the differences between the fly and human GCase enzymes, as evidenced by the distinct phenotypic outcomes of mutations in flies compared to those observed in human GD patients.

Pervasive tissue-, genetic background-, and allele-specific gene expression effects in Drosophila melanogaster

The pervasiveness of gene expression variation and its contribution to phenotypic variation and evolution is well known. This gene expression variation is context dependent, with differences in regulatory architecture often associated with intrinsic and environmental factors, and is modulated by regulatory elements that can act in cis (linked) or in trans (unlinked) relative to the genes they affect. So far, little is known about how this genetic variation affects the evolution of regulatory architecture among closely related tissues during population divergence. To address this question, we analyzed gene expression in the midgut, hindgut, and Malpighian tubule as well as microbiome composition in the two gut tissues in four Drosophila melanogaster strains and their F1 hybrids from two divergent populations: one from the derived, European range and one from the ancestral, African range. In both the transcriptome and microbiome data, we detected extensive tissue- and genetic background-specific effects, including effects of genetic background on overall tissue specificity. Tissue-specific effects were typically stronger than genetic background-specific effects, although the two gut tissues were not more similar to each other than to the Malpighian tubules. An examination of allele specific expression revealed that, while both cis and trans effects were more tissue-specific in genes expressed differentially between populations than genes with conserved expression, trans effects were more tissue-specific than cis effects. Despite there being highly variable regulatory architecture, this observation was robust across tissues and genetic backgrounds, suggesting that the expression of trans variation can be spatially fine-tuned as well as or better than cis variation during population divergence and yielding new insights into cis and trans regulatory evolution.

Dominance reversals: the resolution of genetic conflict and maintenance of genetic variation

Beneficial reversals of dominance reduce the costs of genetic trade-offs and can enable selection to maintain genetic variation for fitness. Beneficial dominance reversals are characterized by the beneficial allele for a given context (e.g. habitat, developmental stage, trait or sex) being dominant in that context but recessive where deleterious. This context dependence at least partially mitigates the fitness consequence of heterozygotes carrying one non-beneficial allele for their context and can result in balancing selection that maintains alternative alleles. Dominance reversals are theoretically plausible and are supported by mounting empirical evidence. Here, we highlight the importance of beneficial dominance reversals as a mechanism for the mitigation of genetic conflict and review the theory and empirical evidence for them. We identify some areas in need of further research and development and outline three methods that could facilitate the identification of antagonistic genetic variation (dominance ordination, allele-specific expression and allele-specific ATAC-Seq (assay for transposase-accessible chromatin with sequencing)). There is ample scope for the development of new empirical methods as well as reanalysis of existing data through the lens of dominance reversals. A greater focus on this topic will expand our understanding of the mechanisms that resolve genetic conflict and whether they maintain genetic variation.

Thermal Adaptations in Animals: Genes, Development, and Evolution

Thermal adaptation to environmental temperature is a driving force in animal evolution. This chapter presents thermal adaptation in ectotherms and endotherms from the perspective of developmental biology. In ectotherms, there are known examples of temperature influencing morphological characteristics, such as seasonal color change, melanization, and sex determination. Furthermore, the timing of embryonic development also varies with environmental temperature. This review will introduce the cellular and molecular mechanisms underlying temperature-dependent embryogenesis. The evolution of thermal adaptation in endotherms is also important for survival in cold climates. Recent genome-wide studies have revealed adaptive mutations in the genomes of extant humans as well as extinct species such as woolly mammoths and Neanderthals. These studies have shown that single-nucleotide polymorphisms in physiologically related genes (e.g., CPT1A, LRP5, THATA, PRKG1, and FADS1-3) allow humans to live in cold climates. At the end of this chapter, we present the remaining questions in terms of genetic assimilation, heat shock protein Hsp90, and embryonic development.

Environmentally robust cis-regulatory changes underlie rapid climatic adaptation

Changes in gene expression are thought to play a major role in adaptive evolution. While it is known that gene expression is highly sensitive to the environment, very few studies have determined the influence of genetic and environmental effects on adaptive gene expression differences in natural populations. Here, we utilize allele-specific expression to characterize cis and trans gene regulatory divergence in temperate and tropical house mice in two metabolic tissues under two thermal conditions. First, we show that gene expression divergence is pervasive between populations and across thermal conditions, with roughly 5 to 10% of genes exhibiting genotype-by-environment interactions. Second, we found that most expression divergence was due to cis-regulatory changes that were stable across temperatures. In contrast, patterns of expression plasticity were largely attributable to trans-effects, which showed greater sensitivity to temperature. Nonetheless, we found a small subset of temperature-dependent cis-regulatory changes, thereby identifying loci underlying expression plasticity. Finally, we performed scans for selection in wild house mice to identify genomic signatures of rapid adaptation. Genomic outliers were enriched in genes with evidence for cis-regulatory divergence. Notably, these genes were associated with phenotypes that affected body weight and metabolism, suggesting that cis-regulatory changes are a possible mechanism for adaptive body size evolution between populations. Our results show that gene expression plasticity, largely controlled in trans, may facilitate the colonization of new environments, but that evolved changes in gene expression are largely controlled in cis, illustrating the genetic and nongenetic mechanisms underlying the establishment of populations in new environments.

Sex-specific estimation of cis and trans regulation of gene expression in heads and gonads of Drosophila melanogaster

The regulatory architecture of gene expression is known to differ substantially between sexes in Drosophila, but most studies performed so far used whole-body data and only single crosses, which may have limited their scope to detect patterns that are robust across tissues and biological replicates. Here, we use allele-specific gene expression of parental and reciprocal hybrid crosses between 6 Drosophila melanogaster inbred lines to quantify cis- and trans-regulatory variation in heads and gonads of both sexes separately across 3 replicate crosses. Our results suggest that female and male heads, as well as ovaries, have a similar regulatory architecture. On the other hand, testes display more and substantially different cis-regulatory effects, suggesting that sex differences in the regulatory architecture that have been previously observed may largely derive from testis-specific effects. We also examine the difference in cis-regulatory variation of genes across different levels of sex bias in gonads and heads. Consistent with the idea that intersex correlations constrain expression and can lead to sexual antagonism, we find more cis variation in unbiased and moderately biased genes in heads. In ovaries, reduced cis variation is observed for male-biased genes, suggesting that cis variants acting on these genes in males do not lead to changes in ovary expression. Finally, we examine the dominance patterns of gene expression and find that sex- and tissue-specific patterns of inheritance as well as trans-regulatory variation are highly variable across biological crosses, although these were performed in highly controlled experimental conditions. This highlights the importance of using various genetic backgrounds to infer generalizable patterns.

The effects of mutations on gene expression and alternative splicing

Understanding the relationship between mutations and their genomic and phenotypic consequences has been a longstanding goal of evolutionary biology. However, few studies have investigated the impact of mutations on gene expression and alternative splicing on the genome-wide scale. In this study, we aim to bridge this knowledge gap by utilizing whole-genome sequencing data and RNA sequencing data from 16 obligately parthenogenetic Daphnia mutant lines to investigate the effects of ethyl methanesulfonate-induced mutations on gene expression and alternative splicing. Using rigorous analyses of mutations, expression changes and alternative splicing, we show that trans-effects are the major contributor to the variance in gene expression and alternative splicing between the wild-type and mutant lines, whereas cis mutations only affected a limited number of genes and do not always alter gene expression. Moreover, we show that there is a significant association between differentially expressed genes and exonic mutations, indicating that exonic mutations are an important driver of altered gene expression.

Fluctuating selection and the determinants of genetic variation

Recent studies of cosmopolitan Drosophila populations have found hundreds to thousands of genetic loci with seasonally fluctuating allele frequencies, bringing temporally fluctuating selection to the forefront of the historical debate surrounding the maintenance of genetic variation in natural populations. Numerous mechanisms have been explored in this longstanding area of research, but these exciting empirical findings have prompted several recent theoretical and experimental studies that seek to better understand the drivers, dynamics, and genome-wide influence of fluctuating selection. In this review, we evaluate the latest evidence for multilocus fluctuating selection in Drosophila and other taxa, highlighting the role of potential genetic and ecological mechanisms in maintaining these loci and their impacts on neutral genetic variation.

A genome-wide RNAi screen for genes important for proliferation of cultured Drosophila cells at low temperature identifies the Ball/VRK protein kinase

A change in ambient temperature is predicted to disrupt cellular homeostasis by affecting all cellular processes in an albeit non-uniform manner. Diffusion is generally less temperature-sensitive than enzymes, for example, and each enzyme has a characteristic individual temperature profile. The actual effects of temperature variation on cells are still poorly understood at the molecular level. Towards an improved understanding, we have performed a genome-wide RNA interference screen with S2R + cells. This Drosophila cell line proliferates over a temperature range comparable to that tolerated by the parental ectothermic organism. Based on effects on cell counts and cell cycle profile after knockdown at 27 and 17 °C, respectively, genes were identified with an apparent greater physiological significance at one or the other temperature. While 27 °C is close to the temperature optimum, the substantially lower 17 °C was chosen to identify genes important at low temperatures, which have received less attention compared to the heat shock response. Among a substantial number of screen hits, we validated a set successfully in cell culture and selected ballchen for further evaluation in the organism. This gene encodes the conserved metazoan VRK protein kinase that is crucial for the release of chromosomes from the nuclear envelope during mitosis. Our analyses in early embryos and larval wing imaginal discs confirmed a higher requirement for ballchen function at temperatures below the optimum. Overall, our experiments validate the genome-wide screen as a basis for future characterizations of genes with increased physiological significance at the lower end of the readily tolerated temperature range.

Genome-wide selection signatures reveal widespread synergistic effects of two different stressors in Drosophila melanogaster

Experimental evolution combined with whole-genome sequencing (evolve and resequence (E&R)) is a powerful approach to study the adaptive architecture of selected traits. Nevertheless, so far the focus has been on the selective response triggered by a single stressor. Building on the highly parallel selection response of founder populations with reduced variation, we evaluated how the presence of a second stressor affects the genomic selection response. After 20 generations of adaptation to laboratory conditions at either 18°C or 29°C, strong genome-wide selection signatures were observed. Only 38% of the selection signatures can be attributed to laboratory adaptation (no difference between temperature regimes). The remaining selection responses are either caused by temperature-specific effects, or reflect the joint effects of temperature and laboratory adaptation (same direction, but the magnitude differs between temperatures). The allele frequency changes resulting from the combined effects of temperature and laboratory adaptation were more extreme in the hot environment for 83% of the affected genomic regions-indicating widespread synergistic effects of the two stressors. We conclude that E&R with reduced genetic variation is a powerful approach to study genome-wide fitness consequences driven by the combined effects of multiple environmental factors.

Deleterious mutations and the rare allele burden on rice gene expression

Deleterious genetic variation is maintained in populations at low frequencies. Under a model of stabilizing selection, rare (and presumably deleterious) genetic variants are associated with increase or decrease in gene expression from some intermediate optimum. We investigate this phenomenon in a population of largely Oryza sativa ssp. indica rice landraces under normal unstressed wet and stressful drought field conditions. We include single nucleotide polymorphisms, insertion/deletion mutations, and structural variants in our analysis and find a stronger association between rare variants and gene expression outliers under the stress condition. We also show an association of the strength of this rare variant effect with linkage, gene expression levels, network connectivity, local recombination rate, and fitness consequence scores, consistent with the stabilizing selection model of gene expression.

Gene expression studies of plastic and evolutionary responses to global warming

Phenotypic plasticity can be a rapid response for coping with global warming, yet may be insufficient to protect species from extinction. Evolutionary adaptation may reinforce adaptive or oppose maladaptive plastic responses. With advances in technology whole transcriptomes can provide us with an unprecedented overview of genes and functional processes underlying the interplay between plasticity and evolution. We advocate that insects provide ideal opportunities to study plasticity in non-adapted and thermally adapted populations to infer reaction norms across the whole transcriptome ('reactionomes'). This can advance our understanding of how the interplay between plasticity and evolution shapes responses to warming. So far, a limited number of studies suggest predominantly maladaptive plastic responses to novel environments that are reduced with time, but much more research is needed to infer general patterns.

Gene Regulatory Evolution in Cold-Adapted Fly Populations Neutralizes Plasticity and May Undermine Genetic Canalization

The relationships between adaptive evolution, phenotypic plasticity, and canalization remain incompletely understood. Theoretical and empirical studies have made conflicting arguments on whether adaptive evolution may enhance or oppose the plastic response. Gene regulatory traits offer excellent potential to study the relationship between plasticity and adaptation, and they can now be studied at the transcriptomic level. Here, we take advantage of three closely related pairs of natural populations of Drosophila melanogaster from contrasting thermal environments that reflect three separate instances of cold tolerance evolution. We measure the transcriptome-wide plasticity in gene expression levels and alternative splicing (intron usage) between warm and cold laboratory environments. We find that suspected adaptive changes in both gene expression and alternative splicing tend to neutralize the ancestral plastic response. Further, we investigate the hypothesis that adaptive evolution can lead to decanalization of selected gene regulatory traits. We find strong evidence that suspected adaptive gene expression (but not splicing) changes in cold-adapted populations are more vulnerable to the genetic perturbation of inbreeding than putatively neutral changes. We find some evidence that these patterns may reflect a loss of genetic canalization accompanying adaptation, although other processes including hitchhiking recessive deleterious variants may contribute as well. Our findings augment our understanding of genetic and environmental effects on gene regulation in the context of adaptive evolution.

Molecular mechanisms underlying plasticity in a thermally varying environment

Adaptation to environmental variability is a prerequisite for species’ persistence in their natural environments. With climate change predicted to increase the frequency and severity of temperature fluctuations, ectothermic organisms may increasingly depend on acclimation capacity to accommodate thermal variability. To elucidate the molecular basis of fluctuating temperature‐induced phenotypic plasticity, we investigated heat tolerance and the mechanisms induced by acclimation to thermal variability as compared to those seen at constant temperature. We ran genome‐wide transcriptomic analysis on Drosophila melanogaster subjected to acclimation at constant (19 ± 0°C) and fluctuating (19 ± 8°C) temperatures and contrasted the induction of molecular mechanisms in adult males, adult females and larvae. We found life stage‐ and sex‐specific dynamics of the acclimation responses to fluctuating temperatures. Adult flies exposed to temperature fluctuations showed a constitutive improvement in heat tolerance while heat tolerance of larvae tracked thermal fluctuations. A constitutive down‐regulation of gene expression was observed for several genes in the larvae exposed to fluctuations. Our results for adult females showed that, for several genes, fluctuating temperature acclimation resulted in canalization of gene expression. Both transcriptional and post‐transcriptional machinery were greatly affected by fluctuations in adult males. Gene ontology analysis showed enrichment of the heat stress response involving several major heat shock proteins in both larvae and adults exposed to fluctuating temperatures, even though fluctuations were in a benign range of temperatures. Finally, molecular mechanisms related to environmental sensing seem to be an important component of insect responses to thermal variability.

Temperature sensitivity of Notch signaling underlies species-specific developmental plasticity and robustness in amniote brains

Ambient temperature significantly affects developmental timing in animals. The temperature sensitivity of embryogenesis is generally believed to be a consequence of the thermal dependency of cellular metabolism. However, the adaptive molecular mechanisms that respond to variations in temperature remain unclear. Here, we report species-specific thermal sensitivity of Notch signaling in the developing amniote brain. Transient hypothermic conditions increase canonical Notch activity and reduce neurogenesis in chick neural progenitors. Increased biosynthesis of phosphatidylethanolamine, a major glycerophospholipid components of the plasma membrane, mediates hypothermia-induced Notch activation. Furthermore, the species-specific thermal dependency of Notch signaling is associated with developmental robustness to altered Notch signaling. Our results reveal unique regulatory mechanisms for temperature-dependent neurogenic potentials that underlie developmental and evolutionary adaptations to a range of ambient temperatures in amniotes.

Ambient temperature significantly affects embryogenesis, but adaptive molecular mechanisms that respond to temperature remain unclear. Here, the authors identified species-specific thermal sensitivity of Notch signaling in developing amniote brains.

Highly Parallel Genomic Selection Response in Replicated Drosophila melanogaster Populations with Reduced Genetic Variation

Many adaptive traits are polygenic and frequently more loci contributing to the phenotype are segregating than needed to express the phenotypic optimum. Experimental evolution with replicated populations adapting to a new controlled environment provides a powerful approach to study polygenic adaptation. Because genetic redundancy often results in nonparallel selection responses among replicates, we propose a modified evolve and resequence (E&R) design that maximizes the similarity among replicates. Rather than starting from many founders, we only use two inbred Drosophila melanogaster strains and expose them to a very extreme, hot temperature environment (29 °C). After 20 generations, we detect many genomic regions with a strong, highly parallel selection response in 10 evolved replicates. The X chromosome has a more pronounced selection response than the autosomes, which may be attributed to dominance effects. Furthermore, we find that the median selection coefficient for all chromosomes is higher in our two-genotype experiment than in classic E&R studies. Because two random genomes harbor sufficient variation for adaptive responses, we propose that this approach is particularly well-suited for the analysis of polygenic adaptation.

A cis-regulatory element promoting increased transcription at low temperature in cultured ectothermic Drosophila cells

Background

Temperature change affects the myriad of concurrent cellular processes in a non-uniform, disruptive manner. While endothermic organisms minimize the challenge of ambient temperature variation by keeping the core body temperature constant, cells of many ectothermic species maintain homeostatic function within a considerable temperature range. The cellular mechanisms enabling temperature acclimation in ectotherms are still poorly understood. At the transcriptional level, the heat shock response has been analyzed extensively. The opposite, the response to sub-optimal temperature, has received lesser attention in particular in animal species. The tissue specificity of transcriptional responses to cool temperature has not been addressed and it is not clear whether a prominent general response occurs. Cis-regulatory elements (CREs), which mediate increased transcription at cool temperature, and responsible transcription factors are largely unknown.

Results

The ectotherm Drosophila melanogaster with a presumed temperature optimum around 25 °C was used for transcriptomic analyses of effects of temperatures at the lower end of the readily tolerated range (14–29 °C). Comparative analyses with adult flies and cell culture lines indicated a striking degree of cell-type specificity in the transcriptional response to cool. To identify potential cis-regulatory elements (CREs) for transcriptional upregulation at cool temperature, we analyzed temperature effects on DNA accessibility in chromatin of S2R+ cells. Candidate cis-regulatory elements (CREs) were evaluated with a novel reporter assay for accurate assessment of their temperature-dependency. Robust transcriptional upregulation at low temperature could be demonstrated for a fragment from the pastrel gene, which expresses more transcript and protein at reduced temperatures. This CRE is controlled by the JAK/STAT signaling pathway and antagonizing activities of the transcription factors Pointed and Ets97D.

Conclusion

Beyond a rich data resource for future analyses of transcriptional control within the readily tolerated range of an ectothermic animal, a novel reporter assay permitting quantitative characterization of CRE temperature dependence was developed. Our identification and functional dissection of the pst_E1 enhancer demonstrate the utility of resources and assay. The functional characterization of this CoolUp enhancer provides initial mechanistic insights into transcriptional upregulation induced by a shift to temperatures at the lower end of the readily tolerated range.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08057-4.

Differences in temperature-sensitive expression of PcG-regulated genes among natural populations of Drosophila melanogaster

Environmental temperature can affect chromatin-based gene regulation, in particular in ectotherms such as insects. Genes regulated by the Polycomb group (PcG) vary in their transcriptional output in response to changes in temperature. Expression of PcG-regulated genes typically increases with decreasing temperatures. Here, we examined variations in temperature-sensitive expression of PcG target genes in natural populations from different climates of Drosophila melanogaster, and differences thereof across different fly stages and tissues. Temperature-induced expression plasticity was found to be stage- and sex-specific with differences in the specificity between the examined PcG target genes. Some tissues and stages, however, showed a higher number of PcG target genes with temperature-sensitive expression than others. Overall, we found higher levels of temperature-induced expression plasticity in African tropical flies from the ancestral species range than in flies from temperate Europe. We also observed differences between temperate flies, however, with more reduction of expression plasticity in warm-temperate than in cold-temperate populations. Although in general, temperature-sensitive expression appeared to be detrimental in temperate climates, there were also cases in which plasticity was increased in temperate flies, as well as no changes in expression plasticity between flies from different climates.

Phenotypic and Transcriptomic Responses to Stress Differ According to Population Geography in an Invasive Species

Adaptation to rapid environmental changes must occur within a short-time scale. In this context, studies of invasive species may provide insights into the underlying mechanisms of rapid adaptation as these species have repeatedly encountered and adapted to novel environmental conditions. We investigated how invasive and noninvasive genotypes of Drosophila suzukii deal with oxidative stress at the phenotypic and molecular levels. We also studied the impact of transposable element (TE) insertions on the gene expression in response to stress. Our results show that flies from invasive areas (France and the United States) live longer in natural conditions than the ones from native Japanese areas. As expected, lifespan for all genotypes was significantly reduced following exposure to paraquat, but this reduction varied among genotypes (genotype-by-environment interaction) with invasive genotypes appearing more affected by exposure than noninvasive ones. A transcriptomic analysis of genotypes upon paraquat treatment detected many genes differentially expressed (DE). Although a small core set of genes were DE in all genotypes following paraquat exposure, much of the response of each genotype was unique. Moreover, we showed that TEs were not activated after oxidative stress and DE genes were significantly depleted of TEs. In conclusion, it is likely that transcriptomic changes are involved in the rapid adaptation to local environments. We provide new evidence that in the decade since the invasion from Asia, the sampled genotypes in Europe and the United States of D. suzukii diverged from the ones from the native area regarding their phenotypic and genomic response to oxidative stress.

Temperature-dependent effects of house fly proto-Y chromosomes on gene expression could be responsible for fitness differences that maintain polygenic sex determination

Sex determination, the developmental process by which sexually dimorphic phenotypes are established, evolves fast. Evolutionary turnover in a sex determination pathway may occur via selection on alleles that are genetically linked to a new master sex determining locus on a newly formed proto-sex chromosome. Species with polygenic sex determination, in which master regulatory genes are found on multiple different proto-sex chromosomes, are informative models to study the evolution of sex determination and sex chromosomes. House flies are such a model system, with male determining loci possible on all six chromosomes and a female-determiner on one of the chromosomes as well. The two most common male-determining proto-Y chromosomes form latitudinal clines on multiple continents, suggesting that temperature variation is an important selection pressure responsible for maintaining polygenic sex determination in this species. Temperature-dependent fitness effects could be manifested through temperature-dependent gene expression differences across proto-Y chromosome genotypes. These gene expression differences may be the result of cis regulatory variants that affect the expression of genes on the proto-sex chromosomes, or trans effects of the proto-Y chromosomes on genes elswhere in the genome. We used RNA-seq to identify genes whose expression depends on proto-Y chromosome genotype and temperature in adult male house flies. We found no evidence for ecologically meaningful temperature-dependent expression differences of sex determining genes between male genotypes, but we were probably not sampling an appropriate developmental time-point to identify such effects. In contrast, we identified many other genes whose expression depends on the interaction between proto-Y chromosome genotype and temperature, including genes that encode proteins involved in reproduction, metabolism, lifespan, stress response, and immunity. Notably, genes with genotype-by-temperature interactions on expression were not enriched on the proto-sex chromosomes. Moreover, there was no evidence that temperature-dependent expression is driven by chromosome-wide cis-regulatory divergence between the proto-Y and proto-X alleles. Therefore, if temperature-dependent gene expression is responsible for differences in phenotypes and fitness of proto-Y genotypes across house fly populations, these effects are driven by a small number of temperature-dependent alleles on the proto-Y chromosomes that may have trans effects on the expression of genes on other chromosomes.

Parallel and Population-specific Gene Regulatory Evolution in Cold-Adapted Fly Populations

Changes in gene regulation at multiple levels may comprise an important share of the molecular changes underlying adaptive evolution in nature. However, few studies have assayed within- and between-population variation in gene regulatory traits at a transcriptomic scale, and therefore inferences about the characteristics of adaptive regulatory changes have been elusive. Here, we assess quantitative trait differentiation in gene expression levels and alternative splicing (intron usage) between three closely related pairs of natural populations of Drosophila melanogaster from contrasting thermal environments that reflect three separate instances of cold tolerance evolution. The cold-adapted populations were known to show population genetic evidence for parallel evolution at the SNP level, and here we find evidence for parallel expression evolution between them, with stronger parallelism at larval and adult stages than for pupae. We also implement a flexible method to estimate cis- vs trans-encoded contributions to expression or splicing differences at the adult stage. The apparent contributions of cis- vs trans-regulation to adaptive evolution vary substantially among population pairs. While two of three population pairs show a greater enrichment of cis-regulatory differences among adaptation candidates, trans-regulatory differences are more likely to be implicated in parallel expression changes between population pairs. Genes with significant cis-effects are enriched for signals of elevated genetic differentiation between cold- and warm-adapted populations, suggesting that they are potential targets of local adaptation. These findings expand our knowledge of adaptive gene regulatory evolution and our ability to make inferences about this important and widespread process.

Testcrosses are an efficient strategy for identifying cis-regulatory variation: Bayesian analysis of allele-specific expression (BayesASE)

Allelic imbalance (AI) occurs when alleles in a diploid individual are differentially expressed and indicates cis acting regulatory variation. What is the distribution of allelic effects in a natural population? Are all alleles the same? Are all alleles distinct? The approach described applies to any technology generating allele-specific sequence counts, for example for chromatin accessibility and can be applied generally including to comparisons between tissues or environments for the same genotype. Tests of allelic effect are generally performed by crossing individuals and comparing expression between alleles directly in the F1. However, a crossing scheme that compares alleles pairwise is a prohibitive cost for more than a handful of alleles as the number of crosses is at least (n2-n)/2 where n is the number of alleles. We show here that a testcross design followed by a hypothesis test of AI between testcrosses can be used to infer differences between nontester alleles, allowing n alleles to be compared with n crosses. Using a mouse data set where both testcrosses and direct comparisons have been performed, we show that the predicted differences between nontester alleles are validated at levels of over 90% when a parent-of-origin effect is present and of 60%-80% overall. Power considerations for a testcross, are similar to those in a reciprocal cross. In all applications, the testing for AI involves several complex bioinformatics steps. BayesASE is a complete bioinformatics pipeline that incorporates state-of-the-art error reduction techniques and a flexible Bayesian approach to estimating AI and formally comparing levels of AI between conditions. The modular structure of BayesASE has been packaged in Galaxy, made available in Nextflow and as a collection of scripts for the SLURM workload manager on github (https://github.com/McIntyre-Lab/BayesASE).

Models for Eco-Evolutionary Extinction Vortices under Balancing Selection

AbstractThe smaller a population is, the faster it loses genetic diversity as a result of genetic drift. Loss of genetic diversity can reduce population growth rate, making populations even smaller and more vulnerable to loss of genetic diversity. Ultimately, the population can be driven to extinction by this "eco-evolutionary extinction vortex." While there are already quantitative models for extinction vortices resulting from inbreeding depression and mutation accumulation, to date extinction vortices resulting from loss of genetic diversity at loci under various forms of balancing selection have been mainly described verbally. To understand better when such extinction vortices arise and to develop methods for detecting them, we propose quantitative eco-evolutionary models, both stochastic individual-based simulations and deterministic approximations, linking loss of genetic diversity and population decline. Using mathematical analysis and simulations, we identify parameter combinations that exhibit strong interactions between population size and genetic diversity and match our definition of an eco-evolutionary vortex (i.e., per capita population decline rates and per-locus fixation rates increase with decreasing population size and number of polymorphic loci). We further highlight cues that may be exhibited by such populations but find that classical early-warning signals are of limited use in detecting populations undergoing an eco-evolutionary extinction vortex.

Allele-Specific Expression and Evolution of Gene Regulation Underlying Acute Heat Stress Response and Local Adaptation in the Copepod Tigriopus californicus

Geographic variation in environmental temperature can select for local adaptation among conspecific populations. Divergence in gene expression across the transcriptome is a key mechanism for evolution of local thermal adaptation in many systems, yet the genetic mechanisms underlying this regulatory evolution remain poorly understood. Here we examine gene expression in 2 locally adapted Tigriopus californicus populations (heat tolerant San Diego, SD, and less tolerant Santa Cruz, SC) and their F1 hybrids during acute heat stress response. Allele-specific expression (ASE) in F1 hybrids was used to determine cis-regulatory divergence. We found that the number of genes showing significant allelic imbalance increased under heat stress compared to unstressed controls. This suggests that there is significant population divergence in cis-regulatory elements underlying heat stress response. Specifically, the number of genes showing an excess of transcripts from the more thermal tolerant (SD) population increased with heat stress while that number of genes with an SC excess was similar in both treatments. Inheritance patterns of gene expression also revealed that genes displaying SD-dominant expression phenotypes increase in number in response to heat stress; that is, across loci, gene expression in F1's following heat stress showed more similarity to SD than SC, a pattern that was absent in the control treatment. The observed patterns of ASE and inheritance of gene expression provide insight into the complex processes underlying local adaptation and thermal stress response.

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.

Molecular and evolutionary processes generating variation in gene expression

Heritable variation in gene expression is common within and between species. This variation arises from mutations that alter the form or function of molecular gene regulatory networks that are then filtered by natural selection. High-throughput methods for introducing mutations and characterizing their cis- and trans-regulatory effects on gene expression (particularly, transcription) are revealing how different molecular mechanisms generate regulatory variation, and studies comparing these mutational effects with variation seen in the wild are teasing apart the role of neutral and non-neutral evolutionary processes. This integration of molecular and evolutionary biology allows us to understand how the variation in gene expression we see today came to be and to predict how it is most likely to evolve in the future.

Ocean acidification induces distinct transcriptomic responses across life history stages of the sea urchin Heliocidaris erythrogramma

Ocean acidification (OA) from seawater uptake of rising carbon dioxide emissions impairs development in marine invertebrates, particularly in calcifying species. Plasticity in gene expression is thought to mediate many of these physiological effects, but how these responses change across life history stages remains unclear. The abbreviated lecithotrophic development of the sea urchin Heliocidaris erythrogramma provides a valuable opportunity to analyse gene expression responses across a wide range of life history stages, including the benthic, post-metamorphic juvenile. We measured the transcriptional response to OA in H. erythrogramma at three stages of the life cycle (embryo, larva, and juvenile) in a controlled breeding design. The results reveal a broad range of strikingly stage-specific impacts of OA on transcription, including changes in the number and identity of affected genes; the magnitude, sign, and variance of their expression response; and the developmental trajectory of expression. The impact of OA on transcription was notably modest in relation to gene expression changes during unperturbed development and much smaller than genetic contributions from parentage. The latter result suggests that natural populations may provide an extensive genetic reservoir of resilience to OA. Taken together, these results highlight the complexity of the molecular response to OA, its substantial life history stage specificity, and the importance of contextualizing the transcriptional response to pH stress in light of normal development and standing genetic variation to better understand the capacity for marine invertebrates to adapt to OA.

Shadow enhancers can suppress input transcription factor noise through distinct regulatory logic

Shadow enhancers, groups of seemingly redundant enhancers, are found in a wide range of organisms and are critical for robust developmental patterning. However, their mechanism of action is unknown. We hypothesized that shadow enhancers drive consistent expression levels by buffering upstream noise through a separation of transcription factor (TF) inputs at the individual enhancers. By measuring the transcriptional dynamics of several Kruppel shadow enhancer configurations in live Drosophila embryos, we showed that individual member enhancers act largely independently. We found that TF fluctuations are an appreciable source of noise that the shadow enhancer pair can better buffer than duplicated enhancers. The shadow enhancer pair is also uniquely able to maintain low levels of expression noise across a wide range of temperatures. A stochastic model demonstrated the separation of TF inputs is sufficient to explain these findings. Our results suggest the widespread use of shadow enhancers is partially due to their noise suppressing ability.

Genetic Analysis of the Transition from Wild to Domesticated Cotton (Gossypium hirsutum L.)

The evolution and domestication of cotton is of great interest from both economic and evolutionary standpoints. Although many genetic and genomic resources have been generated for cotton, the genetic underpinnings of the transition from wild to domesticated cotton remain poorly known. Here we generated an intraspecific QTL mapping population specifically targeting domesticated cotton phenotypes. We used 466 F₂ individuals derived from an intraspecific cross between the wild Gossypium hirsutum var. yucatanense (TX2094) and the elite cultivar G. hirsutum cv. Acala Maxxa, in two environments, to identify 120 QTL associated with phenotypic changes under domestication. While the number of QTL recovered in each subpopulation was similar, only 22 QTL were considered coincident (i.e., shared) between the two locations, eight of which shared peak markers. Although approximately half of QTL were located in the A-subgenome, many key fiber QTL were detected in the D-subgenome, which was derived from a species with unspinnable fiber. We found that many QTL are environment-specific, with few shared between the two environments, indicating that QTL associated with G. hirsutum domestication are genomically clustered but environmentally labile. Possible candidate genes were recovered and are discussed in the context of the phenotype. We conclude that the evolutionary forces that shape intraspecific divergence and domestication in cotton are complex, and that phenotypic transformations likely involved multiple interacting and environmentally responsive factors.

Constitutive and Plastic Gene Expression Variation Associated with Desiccation Resistance Differences in the Drosophila americana Species Group

Stress response mechanisms are ubiquitous and important for adaptation to heterogenous environments and could be based on constitutive or plastic responses to environmental stressors. Here we quantify constitutive and plastic gene expression differences under ambient and desiccation stress treatments, in males and females of three species of Drosophila known to differ in desiccation resistance. Drosophila novamexicana survives desiccation trials significantly longer than the two subspecies of Drosophila americana, consistent with its natural species range in the desert southwest USA. We found that desiccation stress reduces global expression differences between species—likely because many general stress response mechanisms are shared among species—but that all species showed plastic expression changes at hundreds of loci during desiccation. Nonetheless, D. novamexicana had the fewest genes with significant plastic expression changes, despite having the highest desiccation resistance. Of the genes that were significantly differentially expressed between species—either within each treatment (>200 loci), constitutively regardless of treatment (36 loci), or with different species-specific plasticity (26 loci)—GO analysis did not find significant enrichment of any major gene pathways or broader functions associated with desiccation stress. Taken together, these data indicate that if gene expression changes contribute to differential desiccation resistance between species, these differences are likely shaped by a relatively small set of influential genes rather than broad genome-wide differentiation in stress response mechanisms. Finally, among the set of genes with the greatest between-species plasticity, we identified an interesting set of immune-response genes with consistent but opposing reaction norms between sexes, whose potential functional role in sex-specific mechanisms of desiccation resistance remains to be determined.

Compensatory trans-regulatory alleles minimizing variation in TDH3 expression are common within Saccharomyces cerevisiae

Heritable variation in gene expression is common within species. Much of this variation is due to genetic differences outside of the gene with altered expression and is trans-acting. This trans-regulatory variation is often polygenic, with individual variants typically having small effects, making the genetic architecture and evolution of trans-regulatory variation challenging to study. Consequently, key questions about trans-regulatory variation remain, including the variability of trans-regulatory variation within a species, how selection affects trans-regulatory variation, and how trans-regulatory variants are distributed throughout the genome and within a species. To address these questions, we isolated and measured trans-regulatory differences affecting TDH3 promoter activity among 56 strains of Saccharomyces cerevisiae, finding that trans-regulatory backgrounds varied approximately twofold in their effects on TDH3 promoter activity. Comparing this variation to neutral models of trans-regulatory evolution based on empirical measures of mutational effects revealed that despite this variability in the effects of trans-regulatory backgrounds, stabilizing selection has constrained trans-regulatory differences within this species. Using a powerful quantitative trait locus mapping method, we identified ∼100 trans-acting expression quantitative trait locus in each of three crosses to a common reference strain, indicating that regulatory variation is more polygenic than previous studies have suggested. Loci altering expression were located throughout the genome, and many loci were strain specific. This distribution and prevalence of alleles is consistent with recent theories about the genetic architecture of complex traits. In all mapping experiments, the nonreference strain alleles increased and decreased TDH3 promoter activity with similar frequencies, suggesting that stabilizing selection maintained many trans-acting variants with opposing effects. This variation may provide the raw material for compensatory evolution and larger scale regulatory rewiring observed in developmental systems drift among species.

Novel approach to quantitative spatial gene expression uncovers genetic stochasticity in the developing Drosophila eye

Robustness in development allows for the accumulation of genetically based variation in expression. However, this variation is usually examined in response to large perturbations, and examination of this variation has been limited to being spatial, or quantitative, but because of technical restrictions not both. Here we bridge these gaps by investigating replicated quantitative spatial gene expression using rigorous statistical models, in different genotypes, sexes, and species (Drosophila melanogaster and D. simulans). Using this type of quantitative approach with molecular developmental data allows for comparison among conditions, such as different genetic backgrounds. We apply this approach to the morphogenetic furrow, a wave of differentiation that patterns the developing eye disc. Within the morphogenetic furrow, we focus on four genes, hairy, atonal, hedgehog, and Delta. Hybridization chain reaction quantitatively measures spatial gene expression, co-staining for all four genes simultaneously. We find considerable variation in the spatial expression pattern of these genes in the eye between species, genotypes, and sexes. We also find that there has been evolution of the regulatory relationship between these genes, and that their spatial interrelationships have evolved between species. This variation has no phenotypic effect, and could be buffered by network thresholds or compensation from other genes. Both of these mechanisms could potentially be contributing to long term developmental systems drift.

Dominance reversals and the maintenance of genetic variation for fitness

Antagonistic selection between different fitness components (e.g., survival versus fertility) or different types of individuals in a population (e.g., females versus males) can potentially maintain genetic diversity and thereby account for the high levels of fitness variation observed in natural populations. However, the degree to which antagonistic selection can maintain genetic variation critically depends on the dominance relations between antagonistically selected alleles in diploid individuals. Conditions for stable polymorphism of antagonistically selected alleles are narrow, particularly when selection is weak, unless the alleles exhibit "dominance reversals"-in which each allele is partially or completely dominant in selective contexts in which it is favored and recessive in contexts in which it is harmful. Although theory predicts that dominance reversals should emerge under biologically plausible conditions, evidence for dominance reversals is sparse. In this primer, we review theoretical arguments and data supporting a role for dominance reversals in the maintenance of genetic variation. We then highlight an illuminating new study by Grieshop and Arnqvist, which reports a genome-wide signal of dominance reversals between male and female fitness in seed beetles.

A 24 h Age Difference Causes Twice as Much Gene Expression Divergence as 100 Generations of Adaptation to a Novel Environment

Gene expression profiling is one of the most reliable high-throughput phenotyping methods, allowing researchers to quantify the transcript abundance of expressed genes. Because many biotic and abiotic factors influence gene expression, it is recommended to control them as tightly as possible. Here, we show that a 24 h age difference of Drosophilasimulans females that were subjected to RNA sequencing (RNA-Seq) five and six days after eclosure resulted in more than 2000 differentially expressed genes. This is twice the number of genes that changed expression during 100 generations of evolution in a novel hot laboratory environment. Importantly, most of the genes differing in expression due to age introduce false positives or negatives if an adaptive gene expression analysis is not controlled for age. Our results indicate that tightly controlled experimental conditions, including precise developmental staging, are needed for reliable gene expression analyses, in particular in an evolutionary framework.

Systemic signalling and local effectors in developmental stability, body symmetry, and size

Symmetric growth and the origins of fluctuating asymmetry are unresolved phenomena of biology. Small, and sometimes noticeable, deviations from perfect bilateral symmetry reflect the vulnerability of development to perturbations. The degree of asymmetry is related to the magnitude of the perturbations and the ability of an individual to cope with them. As the left and right sides of an individual were presumed to be genetically identical, deviations of symmetry were traditionally attributed to non-genetic effects such as environmental and developmental noise. In this review, we draw attention to other possible sources of variability, especially to somatic mutations and transposons. Mutations are a major source of phenotypic variability and recent genomic data have highlighted somatic mutations as ubiquitous, even in phenotypically normal individuals. We discuss the importance of factors that are responsible for buffering and stabilizing the genome and for maintaining size robustness and quality through elimination of less-fit or damaged cells. However, the important question that arises from these studies is whether this self-correcting capacity and intrinsic organ size controls are sufficient to explain how symmetric structures can reach an identical size and shape. Indeed, recent discoveries in the fruit fly have uncovered a conserved hormone of the insulin/IGF/relaxin family, Dilp8, that is responsible for stabilizing body size and symmetry in the face of growth perturbations. Dilp8 alarm signals periphery growth status to the brain, where it acts on its receptor Lgr3. Loss of Dilp8-Lgr3 signaling renders flies incapable of detecting growth perturbations and thus maintaining a stable size and symmetry. These findings help to understand how size and symmetry of somatic tissues remain undeterred in noisy environments, after injury or illnesses, and in the presence of accumulated somatic mutations.

Prenatal programming of postnatal plasticity revisited-And extended

Two sets of evidence reviewed herein, one indicating that prenatal stress is associated with elevated behavioral and physiological dysregulation and the other that such phenotypic functioning is itself associated with heightened susceptibility to positive and negative environmental influences postnatally, raises the intriguing hypothesis first advanced by Pluess and Belsky (2011) that prenatal stress fosters, promotes, or "programs" postnatal developmental plasticity. Here we review further evidence consistent with this proposition, including new experimental research systematically manipulating both prenatal stress and postnatal rearing. Collectively this work would seem to explain why prenatal stress has so consistently been linked to problematic development: stresses encountered prenatally are likely to continue postnatally, thereby adversely affecting the development of children programmed (by prenatal stress) to be especially susceptible to environmental effects. Less investigated are the potential benefits prenatal stress may promote, due to increased plasticity, when the postnatal environment proves to be favorable. Future directions of research pertaining to potential mechanisms instantiating postnatal plasticity and moderators of such prenatal-programming effects are outlined.

Transgenerational Effects of Extended Dauer Diapause on Starvation Survival and Gene Expression Plasticity in Caenorhabditis elegans

Phenotypic plasticity is facilitated by epigenetic regulation, and remnants of such regulation may persist after plasticity-inducing cues are gone. However, the relationship between plasticity and transgenerational epigenetic memory is not understood. Dauer diapause in Caenorhabditis elegans provides an opportunity to determine how a plastic response to the early-life environment affects traits later in life and in subsequent generations. We report that, after extended diapause, postdauer worms initially exhibit reduced reproductive success and greater interindividual variation. In contrast, F3 progeny of postdauers display increased starvation resistance and lifespan, revealing potentially adaptive transgenerational effects. Transgenerational effects are dependent on the duration of diapause, indicating an effect of extended starvation. In agreement, RNA-seq demonstrates a transgenerational effect on nutrient-responsive genes. Further, postdauer F3 progeny exhibit reduced gene expression plasticity, suggesting a trade-off between plasticity and epigenetic memory. This work reveals complex effects of nutrient stress over different time scales in an animal that evolved to thrive in feast and famine.

Gene Regulatory Variation in Drosophila melanogaster Renal Tissue

Genetic variation influencing levels of gene expression is abundant in natural populations, and may exert its effects through complex mechanisms that depend on an organism's genetic background and the tissue in which expression is measured. We investigated natural variation in gene expression in the Malpighian tubules of three inbred Drosophila melanogaster strains and their F₁ hybrids. One of the strains was from a population in the species' ancestral range (Zambia), while the other two were from a more recently derived population (Sweden). Although closely related, the two Swedish strains differed greatly in terms of their expression inheritance when hybridized with the Zambian strain, with one Swedish strain showing a large excess of genes with recessive expression inheritance, as well as a large number of genes with overdominant inheritance. Although most expression variation could be attributed to trans-regulation, there were ∼200 genes that showed allele-specific expression differences in each of the between-population hybrids, indicating that cis-regulation contributes as well. The cis-regulated genes were enriched with cytochrome P450 genes, and the upstream regions of six of these genes were incorporated into transgenic reporter gene constructs to test their effects on expression. Differential expression was observed for five of the six reporter genes in the Malpighian tubule, suggesting that a large proportion of cis-regulatory variation lies directly upstream of the affected gene. In most cases, the differential expression was specific to the Malpighian tubule or greater in this tissue than in the rest of the body, highlighting the importance of single-tissue studies of gene expression variation.

The Evolution of Gene Expression in cis and trans

There is abundant variation in gene expression between individuals, populations, and species. The evolution of gene regulation and expression within and between species is thought to frequently contribute to adaptation. Yet considerable evidence suggests that the primary evolutionary force acting on variation in gene expression is stabilizing selection. We review here the results of recent studies characterizing the evolution of gene expression occurring in cis (via linked polymorphisms) or in trans (through diffusible products of other genes) and their contribution to adaptation and response to the environment. We review the evidence for buffering of variation in gene expression at the level of both transcription and translation, and the possible mechanisms for this buffering. Lastly, we summarize unresolved questions about the evolution of gene regulation.

Molecular dissection of a natural transposable element invasion

The first tracking of the dynamics of a natural invasion by a transposable element (TE) provides unprecedented details on the establishment of host defense mechanisms against TEs. We captured a D. simulans population at an early stage of a P-element invasion and studied the spread of the TE in replicated experimentally evolving populations kept under hot and cold conditions. We analyzed the factors controlling the invasion by NGS, RNA-FISH, and gonadal dysgenesis assays. Under hot conditions, the P-element spread rapidly for 20 generations, but no further spread was noted later on. This plateauing of the invasion was mediated by the rapid emergence of P-element-specific piRNAs. Under cold conditions, we observed a lower expression of the P-element and a slower emergence of the piRNA defense, resulting in a three times slower invasion that continued beyond 40 generations. We conclude that the environment is a major factor determining the evolution of TEs in their host.

Strong phenotypic plasticity limits potential for evolutionary responses to climate change

Phenotypic plasticity, the expression of multiple phenotypes from one genome, is a widespread adaptation to short-term environmental fluctuations, but whether it facilitates evolutionary adaptation to climate change remains contentious. Here, we investigate seasonal plasticity and adaptive potential in an Afrotropical butterfly expressing distinct phenotypes in dry and wet seasons. We assess the transcriptional architecture of plasticity in a full-factorial analysis of heritable and environmental effects across 72 individuals, and reveal pervasive gene expression differences between the seasonal phenotypes. Strikingly, intra-population genetic variation for plasticity is largely absent, consistent with specialisation to a particular environmental cue reliably predicting seasonal transitions. Under climate change, deteriorating accuracy of predictive cues will likely aggravate maladaptive phenotype-environment mismatches and increase selective pressures on reaction norms. However, the observed paucity of genetic variation for plasticity limits evolutionary responses, potentially weakening prospects for population persistence. Thus, seasonally plastic species may be especially vulnerable to climate change.

Direct Testing for Allele-Specific Expression Differences Between Conditions

Allelic imbalance (AI) indicates the presence of functional variation in cis regulatory regions. Detecting cis regulatory differences using AI is widespread, yet there is no formal statistical methodology that tests whether AI differs between conditions. Here, we present a novel model and formally test differences in AI across conditions using Bayesian credible intervals. The approach tests AI by environment (G×E) interactions, and can be used to test AI between environments, genotypes, sex, and any other condition. We incorporate bias into the modeling process. Bias is allowed to vary between conditions, making the formulation of the model general. As gene expression affects power for detection of AI, and, as expression may vary between conditions, the model explicitly takes coverage into account. The proposed model has low type I and II error under several scenarios, and is robust to large differences in coverage between conditions. We reanalyze RNA-seq data from a Drosophila melanogaster population panel, with F1 genotypes, to compare levels of AI between mated and virgin female flies, and we show that AI × genotype interactions can also be tested. To demonstrate the use of the model to test genetic differences and interactions, a formal test between two F1s was performed, showing the expected 20% difference in AI. The proposed model allows a formal test of G×E and G×G, and reaffirms a previous finding that cis regulation is robust between environments.

Seasonally fluctuating selection can maintain polymorphism at many loci via segregation lift

Most natural populations are affected by seasonal changes in temperature, rainfall, or resource availability. Seasonally fluctuating selection could potentially make a large contribution to maintaining genetic polymorphism in populations. However, previous theory suggests that the conditions for multilocus polymorphism are restrictive. Here, we explore a more general class of models with multilocus seasonally fluctuating selection in diploids. In these models, the multilocus genotype is mapped to fitness in two steps. The first mapping is additive across loci and accounts for the relative contributions of heterozygous and homozygous loci-that is, dominance. The second step uses a nonlinear fitness function to account for the strength of selection and epistasis. Using mathematical analysis and individual-based simulations, we show that stable polymorphism at many loci is possible if currently favored alleles are sufficiently dominant. This general mechanism, which we call "segregation lift," requires seasonal changes in dominance, a phenomenon that may arise naturally in situations with antagonistic pleiotropy and seasonal changes in the relative importance of traits for fitness. Segregation lift works best under diminishing-returns epistasis, is not affected by problems of genetic load, and is robust to differences in parameters across loci and seasons. Under segregation lift, loci can exhibit conspicuous seasonal allele-frequency fluctuations, but often fluctuations may be small and hard to detect. An important direction for future work is to formally test for segregation lift in empirical data and to quantify its contribution to maintaining genetic variation in natural populations.

Unexpected high genetic diversity in small populations suggests maintenance by associative overdominance

The effective population size (N_e ) is a central factor in determining maintenance of genetic variation. The neutral theory predicts that loss of variation depends on N_e , with less genetic drift in larger populations. We monitored genetic drift in 42 Drosophila melanogaster populations of different adult census population sizes (10, 50 or 500) using pooled RAD sequencing. In small populations, variation was lost at a substantially lower rate than expected. This observation was consistent across two ecological relevant thermal regimes, one stable and one with a stressful increase in temperature across generations. Estimated ratios between N_e and adult census size were consistently higher in small than in larger populations. The finding provides evidence for a slower than expected loss of genetic diversity and consequently a higher than expected long-term evolutionary potential in small fragmented populations. More genetic diversity was retained in areas of low recombination, suggesting that associative overdominance, driven by disfavoured homozygosity of recessive deleterious alleles, is responsible for the maintenance of genetic diversity in smaller populations. Consistent with this hypothesis, the X-chromosome, which is largely free of recessive deleterious alleles due to hemizygosity in males, fits neutral expectations even in small populations. Our experiments provide experimental answers to a range of unexpected patterns in natural populations, ranging from variable diversity on X-chromosomes and autosomes to surprisingly high levels of nucleotide diversity in small populations.