Phenotypic plasticity of the Drosophila transcriptome

Impact of tRNA-induced proline-to-serine mistranslation on the transcriptome of Drosophila melanogaster

Mistranslation is the misincorporation of an amino acid into a polypeptide. Mistranslation has diverse effects on multicellular eukaryotes and is implicated in several human diseases. In Drosophila melanogaster, a serine transfer RNA (tRNA) that misincorporates serine at proline codons (P→S) affects male and female flies differently. The mechanisms behind this discrepancy are currently unknown. Here, we compare the transcriptional response of male and female flies to P→S mistranslation to identify genes and cellular processes that underlie sex-specific differences. Both males and females downregulate genes associated with various metabolic processes in response to P→S mistranslation. Males downregulate genes associated with extracellular matrix organization and response to negative stimuli such as wounding, whereas females downregulate aerobic respiration and ATP synthesis genes. Both sexes upregulate genes associated with gametogenesis, but females also upregulate cell cycle and DNA repair genes. These observed differences in the transcriptional response of male and female flies to P→S mistranslation have important implications for the sex-specific impact of mistranslation on disease and tRNA therapeutics.

Transcriptomics of differences in thermal plasticity associated with selection for an exaggerated male sexual trait

The information about the magnitude of differences in thermal plasticity both between and within populations, as well as identification of the underlying molecular mechanisms are key to understanding the evolution of thermal plasticity. In particular, genes underlying variation in the physiological response to temperature can provide raw material for selection acting on plastic traits. Using RNAseq, we investigate the transcriptional response to temperature in males and females from bulb mite populations selected for the increased frequency of one of two discrete male morphs (fighter- and scrambler-selected populations) that differ in relative fitness depending on temperature. We show that different mechanisms underlie the divergence in thermal response between fighter- and scrambler-selected populations at decreased vs. increased temperature. Temperature decrease to 18 °C was associated with higher transcriptomic plasticity of males with more elaborate armaments, as indicated by a significant selection-by-temperature interaction effect on the expression of 40 genes, 38 of which were upregulated in fighter-selected populations in response to temperature decrease. In response to 28 °C, no selection-by-temperature interaction in gene expression was detected. Hence, differences in phenotypic response to temperature increase likely depended on genes associated with their distinct morph-specific thermal tolerance. Selection of males also drove gene expression patterns in females. These patterns could be associated with temperature-dependent fitness differences between females from fighter- vs. scrambler-selected populations reported in previous studies. Our study shows that selection for divergent male sexually selected morphologies and behaviors has a potential to drive divergence in metabolic pathways underlying plastic response to temperature in both sexes.

Omics data analysis reveals the system-level constraint on cellular amino acid composition

Proteins play a pivotal role in coordinating the functions of organisms, essentially governing their traits, as the dynamic arrangement of diverse amino acids leads to a multitude of folded configurations within peptide chains. Despite dynamic changes in amino acid composition of an individual protein (referred to as AAP) and great variance in protein expression levels under different conditions, our study, utilizing transcriptomics data from four model organisms uncovers surprising stability in the overall amino acid composition of the total cellular proteins (referred to as AACell). Although this value may vary between different species, we observed no significant differences among distinct strains of the same species. This indicates that organisms enforce system-level constraints to maintain a consistent AACell, even amid fluctuations in AAP and protein expression. Further exploration of this phenomenon promises insights into the intricate mechanisms orchestrating cellular protein expression and adaptation to varying environmental challenges.

Impact of tRNA-induced proline-to-serine mistranslation on the transcriptome of Drosophila melanogaster

Mistranslation is the misincorporation of an amino acid into a polypeptide. Mistranslation has diverse effects on multicellular eukaryotes and is implicated in several human diseases. In Drosophila melanogaster, a serine transfer RNA (tRNA) that misincorporates serine at proline codons (P→S) affects male and female flies differently. The mechanisms behind this discrepancy are currently unknown. Here, we compare the transcriptional response of male and female flies to P→S mistranslation to identify genes and cellular processes that underlie sex-specific differences. Both males and females downregulate genes associated with various metabolic processes in response to P→S mistranslation. Males downregulate genes associated with extracellular matrix organization and response to negative stimuli such as wounding, whereas females downregulate aerobic respiration and ATP synthesis genes. Both sexes upregulate genes associated with gametogenesis, but females also upregulate cell cycle and DNA repair genes. These observed differences in the transcriptional response of male and female flies to P→S mistranslation have important implications for the sex-specific impact of mistranslation on disease and tRNA therapeutics.

Genetic, Environmental, and Stochastic Components of Lifespan Variability: The Drosophila Paradigm

Lifespan is a complex quantitative trait involving genetic and non-genetic factors as well as the peculiarities of ontogenesis. As with all quantitative traits, lifespan shows considerable variation within populations and between individuals. Drosophila, a favourite object of geneticists, has greatly advanced our understanding of how different forms of variability affect lifespan. This review considers the role of heritable genetic variability, phenotypic plasticity and stochastic variability in controlling lifespan in Drosophila melanogaster. We discuss the major historical milestones in the development of the genetic approach to study lifespan, the breeding of long-lived lines, advances in lifespan QTL mapping, the environmental factors that have the greatest influence on lifespan in laboratory maintained flies, and the mechanisms, by which individual development affects longevity. The interplay between approaches to study ageing and lifespan limitation will also be discussed. Particular attention will be paid to the interaction of different types of variability in the control of lifespan.

Perspectives on the Drosophila melanogaster Model for Advances in Toxicological Science

The use of Drosophila melanogaster for studies of toxicology has grown considerably in the last decade. The Drosophila model has long been appreciated as a versatile and powerful model for developmental biology and genetics because of its ease of handling, short life cycle, low cost of maintenance, molecular genetic accessibility, and availability of a wide range of publicly available strains and data resources. These features, together with recent unique developments in genomics and metabolomics, make the fly model especially relevant and timely for the development of new approach methodologies and movements toward precision toxicology. Here, we offer a perspective on how flies can be leveraged to identify risk factors relevant to environmental exposures and human health. First, we review and discuss fundamental toxicologic principles for experimental design with Drosophila. Next, we describe quantitative and systems genetics approaches to resolve the genetic architecture and candidate pathways controlling susceptibility to toxicants. Finally, we summarize the current state and future promise of the emerging field of Drosophila metabolomics for elaborating toxic mechanisms. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC.

Dissecting the effect of soil on plant phenology and berry transcriptional plasticity in two Italian grapevine varieties (Vitis vinifera L.)

Grapevine embodies a fascinating species as regards phenotypic plasticity and genotype-per-environment interactions. The terroir, namely the set of agri-environmental factors to which a variety is subjected, can influence the phenotype at the physiological, molecular, and biochemical level, representing an important phenomenon connected to the typicality of productions. We investigated the determinants of plasticity by conducting a field-experiment where all terroir variables, except soil, were kept as constant as possible. We isolated the effect of soils collected from different areas, on phenology, physiology, and transcriptional responses of skin and flesh of a red and a white variety of great economic value: Corvina and Glera. Molecular results, together with physio-phenological parameters, suggest a specific effect of soil on grapevine plastic response, highlighting a higher transcriptional plasticity of Glera in respect to Corvina and a marked response of skin compared to flesh. Using a novel statistical approach, we identified clusters of plastic genes subjected to the specific influence of soil. These findings could represent an issue of applicative value, posing the basis for targeted agricultural practices to enhance the desired characteristics for any soil/cultivar combination, to improve vineyards management for a better resource usage and to valorize vineyards uniqueness maximizing the terroir-effect.

The genetic architecture of behavioral canalization

Behaviors are components of fitness and contribute to adaptive evolution. Behaviors represent the interactions of an organism with its environment, yet innate behaviors display robustness in the face of environmental change, which we refer to as 'behavioral canalization'. We hypothesize that positive selection of hub genes of genetic networks stabilizes the genetic architecture for innate behaviors by reducing variation in the expression of interconnected network genes. Robustness of these stabilized networks would be protected from deleterious mutations by purifying selection or suppressing epistasis. We propose that, together with newly emerging favorable mutations, epistatically suppressed mutations can generate a reservoir of cryptic genetic variation that could give rise to decanalization when genetic backgrounds or environmental conditions change to allow behavioral adaptation.

Differential expression of Phlebotomus tobbi Adler, Theodor & Lourie, 1930 (Diptera: Psychodidae) genes under different environmental conditions

Phlebotomus tobbi is a widely distributed sand fly species in Turkey and is the proven vector of Leishmania infantum and several Phleboviruses. Information regarding the genetic basis of phenotypic plasticity is crucial for managing vector-borne diseases, as the changing environmental conditions have consequences for the survival of arthropods and the disease agents they transmit. However, limited data is available on the impacts of environmental conditions on the traits associated with sand fly survival, reproduction, and vectorial competence. The present study aimed to reveal the changes in the expression levels of three selected P. tobbi genes using laboratory-reared and wild-caught populations. A nervous system protein, Cacophony (PtCac), related to the life history traits of sand flies, and two sand fly salivary protein genes, PtSP32 and PtSP38, influence the infection of the vertebrate hosts, were assessed. Sand flies were maintained at 23 °C and 27 °C in the laboratory to evaluate the relationship between temperature and the expressed phenotypes. Field collections were carried out in three climatically distinct regions of Turkey to establish the regional differences in the gene expression levels of natural P. tobbi populations. In the laboratory, PtCac expression increased with the temperature. However, PtCac expression was negatively correlated with local temperature and humidity conditions. No differences were detected in the PtSP32 gene expression levels of both laboratory-reared and wild-caught females, but a negative correlation was observed with relative humidity in natural populations. Although the expression levels of PtSP38 did not differ among the females collected from distinct regions, a positive correlation was detected in the laboratory-reared colony. The findings indicated that changes in environmental conditions could drive the expression levels of P. tobbi genes, which influence population dynamics and the transmission risk of the disease.

Pleiotropic fitness effects of a Drosophila odorant-binding protein

Insect odorant-binding proteins (OBPs) are members of a rapidly evolving multigene family traditionally thought to facilitate chemosensation. However, studies on Drosophila have shown that members of this family have evolved functions beyond chemosensation, as evident from their expression in reproductive tissues and the brain. Previous studies implicated diverse functions of Obp56h, a member of the largest gene cluster of the D. melanogaster Obp repertoire. Here, we examined the effect of CRISPR/Cas9-mediated deletion of Obp56h on 2 fitness phenotypes, on resistance to starvation stress and heat stress, and on locomotion and sleep phenotypes. Obp56h-/- mutants show a strong sexually dimorphic effect on starvation stress survival, with females being more resistant to starvation stress than the control. In contrast, Obp56h-/- females, but not males, are highly sensitive to heat stress. Both sexes show changes in locomotion and sleep patterns. Transcriptional profiling of RNA from heads of Obp56h-/- flies and the wildtype control reveals differentially expressed genes, including gene products associated with antimicrobial immune responses and members of the Turandot family of stress-induced secreted peptides. In addition, differentially expressed genes of unknown function were identified in both sexes. Genes encoding components of the mitochondrial electron transport chain, cuticular proteins, gene products associated with regulation of feeding behavior (Lst and CCHa2), ribosomal proteins, lncRNAs, snoRNAs, tRNAs, and snRNAs show changes in transcript abundances in Obp56h-/- females. These differentially expressed genes are likely to contribute to Obp56h-mediated effects on the diverse phenotypes that arise upon deletion of this OBP.

Signaling through the dystrophin glycoprotein complex affects the stress-dependent transcriptome in Drosophila

Deficiencies in the human dystrophin glycoprotein complex (DGC), which links the extracellular matrix with the intracellular cytoskeleton, cause muscular dystrophies, a group of incurable disorders associated with heterogeneous muscle, brain and eye abnormalities. Stresses such as nutrient deprivation and aging cause muscle wasting, which can be exacerbated by reduced levels of the DGC in membranes, the integrity of which is vital for muscle health and function. Moreover, the DGC operates in multiple signaling pathways, demonstrating an important function in gene expression regulation. To advance disease diagnostics and treatment strategies, we strive to understand the genetic pathways that are perturbed by DGC mutations. Here, we utilized a Drosophila model to investigate the transcriptomic changes in mutants of four DGC components under temperature and metabolic stress. We identified DGC-dependent genes, stress-dependent genes and genes dependent on the DGC for a proper stress response, confirming a novel function of the DGC in stress-response signaling. This perspective yields new insights into the etiology of muscular dystrophy symptoms, possible treatment directions and a better understanding of DGC signaling and regulation under normal and stress conditions.

Distinct signals of clinal and seasonal allele frequency change at eQTLs in Drosophila melanogaster

Populations of short-lived organisms can respond to spatial and temporal environmental heterogeneity through local adaptation. Local adaptation can be reflected on both phenotypic and genetic levels, and it has been documented in many organisms. Although complex fitness-related phenotypes have been shown to vary across latitudinal clines and seasons in similar ways in Drosophila melanogaster populations, the comparative signals of local adaptation across space and time remain poorly understood. Here, we examined patterns of allele frequency change across a latitudinal cline and between seasons at previously reported expression quantitative trait loci (eQTLs). We divided eQTLs into groups by using differential expression profiles of fly populations collected across latitudinal clines or exposed to different environmental conditions. In general, we find that eQTLs are enriched for clinally varying polymorphisms, and that these eQTLs change in frequency in concordant ways across the cline and in response to starvation and chill-coma. The enrichment of eQTLs among seasonally varying polymorphisms is more subtle, and the direction of allele frequency change at eQTLs appears to be somewhat idiosyncratic. Taken together, we suggest that clinal adaptation at eQTLs is at least partially distinct from seasonal adaptation.

Transcriptomic analysis provides insights into molecular mechanisms of thermal physiology

Physiological trait variation underlies health, responses to global climate change, and ecological performance. Yet, most physiological traits are complex, and we have little understanding of the genes and genomic architectures that define their variation. To provide insight into the genetic architecture of physiological processes, we related physiological traits to heart and brain mRNA expression using a weighted gene co-expression network analysis. mRNA expression was used to explain variation in six physiological traits (whole animal metabolism (WAM), critical thermal maximum (CT_max), and four substrate specific cardiac metabolic rates (CaM)) under 12 °C and 28 °C acclimation conditions. Notably, the physiological trait variations among the three geographically close (within 15 km) and genetically similar F. heteroclitus populations are similar to those found among 77 aquatic species spanning 15–20° of latitude (~ 2,000 km). These large physiological trait variations among genetically similar individuals provide a powerful approach to determine the relationship between mRNA expression and heritable fitness related traits unconfounded by interspecific differences. Expression patterns explained up to 82% of metabolic trait variation and were enriched for multiple signaling pathways known to impact metabolic and thermal tolerance (e.g., AMPK, PPAR, mTOR, FoxO, and MAPK) but also contained several unexpected pathways (e.g., apoptosis, cellular senescence), suggesting that physiological trait variation is affected by many diverse genes.

Using phenotypic plasticity to understand the structure and evolution of the genotype-phenotype map

Deciphering the genotype-phenotype map necessitates relating variation at the genetic level to variation at the phenotypic level. This endeavour is inherently limited by the availability of standing genetic variation, the rate of spontaneous mutation to novo genetic variants, and possible biases associated with induced mutagenesis. An interesting alternative is to instead rely on the environment as a source of variation. Many phenotypic traits change plastically in response to the environment, and these changes are generally underlain by changes in gene expression. Relating gene expression plasticity to the phenotypic plasticity of more integrated organismal traits thus provides useful information about which genes influence the development and expression of which traits, even in the absence of genetic variation. We here appraise the prospects and limits of such an environment-for-gene substitution for investigating the genotype-phenotype map. We review models of gene regulatory networks, and discuss the different ways in which they can incorporate the environment to mechanistically model phenotypic plasticity and its evolution. We suggest that substantial progress can be made in deciphering this genotype-environment-phenotype map, by connecting theory on gene regulatory network to empirical patterns of gene co-expression, and by more explicitly relating gene expression to the expression and development of phenotypes, both theoretically and empirically.

Lack of genotype-by-environment interaction suggests limited potential for evolutionary changes in plasticity in the eastern oyster, Crassostrea virginica

Eastern oysters in the northern Gulf of Mexico are facing rapid environmental changes and can respond to this change via plasticity or evolution. Plasticity can act as an immediate buffer against environmental change, but this buffering could impact the organism's ability to evolve in subsequent generations. While plasticity and evolution are not mutually exclusive, the relative contribution and interaction between them remains unclear. In this study, we investigate the roles of plastic and evolved responses to environmental variation and Perkinsus marinus infection in Crassostrea virginica by using a common garden experiment with 80 oysters from six families outplanted at two field sites naturally differing in salinity. We use growth data, P. marinus infection intensities, 3' RNA sequencing (TagSeq) and low-coverage whole-genome sequencing to identify the effect of genotype, environment and genotype-by-environment interaction on the oyster's response to site. As one of first studies to characterize the joint effects of genotype and environment on transcriptomic and morphological profiles in a natural setting, we demonstrate that C. virginica has a highly plastic response to environment and that this response is parallel among genotypes. We also find that genes responding to genotype have distinct and opposing profiles compared to genes responding to environment with regard to expression levels, Ka/Ks ratios and nucleotide diversity. Our findings suggest that C. virginica may be able to buffer the immediate impacts of future environmental changes by altering gene expression and physiology, but the lack of genetic variation in plasticity suggests limited capacity for evolved responses.

Comparative Transcriptome and DNA Methylation Analysis of Phenotypic Plasticity in the Pacific Abalone (Haliotis discus hannai)

Phenotypic plasticity is an adaptive mechanism used by organisms to cope with environmental fluctuations. Pacific abalone (Haliotis discus hannai) are large-scale farmed in the temperate area of northern China and in the warmer waters of southern China. RNA-seq and comparative transcriptomic analysis here were performed to determine if the northern and southern populations have evolved divergent plasticity and if functional differences are associated with protein synthesis and growth-related biological progress. The DNA methylation (5mC) landscape of H. discus hannai from the two populations using whole genomic bisulfite sequencing (WGBS), exhibited different epigenetic patterns. The southern population had significant genomic hypo-methylation that may have resulted from long-term acclimation to heat stress. Combining 790 differentially expressed genes (DEGs) and 7635 differentially methylated genes (DMGs), we found that methylation within the gene body might be important in predicting abalone gene expression. Genes related to growth, development, transduction, and apoptosis may be regulated by methylation and could explain the phenotypic divergence of H. discus hannai. Our findings not only emphasize the significant roles of adaptive plasticity in the acclimation of H. discus hannai to high temperatures but also provide a new understanding of the epigenetic mechanism underlying the phenotypic plasticity in adaptation to climate change for marine organisms.

False and true positives in arthropod thermal adaptation candidate gene lists

Genome-wide studies are prone to false positives due to inherently low priors and statistical power. One approach to ameliorate this problem is to seek validation of reported candidate genes across independent studies: genes with repeatedly discovered effects are less likely to be false positives. Inversely, genes reported only as many times as expected by chance alone, while possibly representing novel discoveries, are also more likely to be false positives. We show that, across over 30 genome-wide studies that reported Drosophila and Daphnia genes with possible roles in thermal adaptation, the combined lists of candidate genes and orthologous groups are rapidly approaching the total number of genes and orthologous groups in the respective genomes. This is consistent with the expectation of high frequency of false positives. The majority of these spurious candidates have been identified by one or a few studies, as expected by chance alone. In contrast, a noticeable minority of genes have been identified by numerous studies with the probabilities of such discoveries occurring by chance alone being exceedingly small. For this subset of genes, different studies are in agreement with each other despite differences in the ecological settings, genomic tools and methodology, and reporting thresholds. We provide a reference set of presumed true positives among Drosophila candidate genes and orthologous groups involved in response to changes in temperature, suitable for cross-validation purposes. Despite this approach being prone to false negatives, this list of presumed true positives includes several hundred genes, consistent with the "omnigenic" concept of genetic architecture of complex traits.

Functional Diversification, Redundancy, and Epistasis among Paralogs of the Drosophila melanogaster Obp50a-d Gene Cluster

Large multigene families, such as the insect odorant-binding proteins (OBPs), are thought to arise through functional diversification after repeated gene duplications. Whereas many OBPs function in chemoreception, members of this family are also expressed in tissues outside chemosensory organs. Paralogs of the Obp50 gene cluster are expressed in metabolic and male reproductive tissues, but their functions and interrelationships remain unknown. Here, we report the genetic dissection of four members of the Obp50 cluster, which are in close physical proximity without intervening genes. We used CRISPR technology to excise the entire cluster while introducing a PhiC31 reintegration site to reinsert constructs in which different combinations of the constituent Obp genes were either intact or rendered inactive. We performed whole transcriptome sequencing and assessed sexually dimorphic changes in transcript abundances (transcriptional niches) associated with each gene-edited genotype. Using this approach, we were able to estimate redundancy, additivity, diversification, and epistasis among Obp50 paralogs. We analyzed the effects of gene editing of this cluster on organismal phenotypes and found a significant skewing of sex ratios attributable to Obp50a, and sex-specific effects on starvation stress resistance attributable to Obp50d. Thus, there is functional diversification within the Obp50 cluster with Obp50a contributing to development and Obp50d to stress resistance. The deletion-reinsertion approach we applied to the Obp50 cluster provides a general paradigm for the genetic dissection of paralogs of multigene families.

Diversifying maize genomic selection models

Genomic selection (GS) is one of the most powerful tools available for maize breeding. Its use of genome-wide marker data to estimate breeding values translates to increased genetic gains with fewer breeding cycles. In this review, we cover the history of GS and highlight particular milestones during its adaptation to maize breeding. We discuss how GS can be applied to developing superior maize inbreds and hybrids. Additionally, we characterize refinements in GS models that could enable the encapsulation of non-additive genetic effects, genotype by environment interactions, and multiple levels of the biological hierarchy, all of which could ultimately result in more accurate predictions of breeding values. Finally, we suggest the stages in a maize breeding program where it would be beneficial to apply GS. Given the current sophistication of high-throughput phenotypic, genotypic, and other -omic level data currently available to the maize community, now is the time to explore the implications of their incorporation into GS models and thus ensure that genetic gains are being achieved as quickly and efficiently as possible.

Dynamic global analysis of transcription reveals the role of miRNAs in synergistic stabilization of gene expression

Buffering exogenous perturbation is crucial to maintain transcriptional homeostasis during development. While miRNAs have been speculated to play a role in stability maintenance, previous studies seeking to check this conjecture focused on measurements of transcript levels at steady state or involved individual miRNA targets. We measured whole-genome expression dynamics by introducing a transient perturbation and establishing a perturbation and recovery system in Drosophila larvae. We inhibited all transcription and assayed transcriptomes at several time points during recovery from inhibition. We performed these experiments in the wild type and miRNA-deficient genetic backgrounds. Consistent with theories about miRNAs' function in stabilizing the transcriptome, we find that attenuating miRNA expression leads to weak impairment in degradation of targets but strong destabilization of target genes when transcription is re-activated. We further fitted a model that captures the essential aspects of transcription dynamics in our experiments and found that the miRNA target transcripts uniformly overshoot the original steady state as they recover from a general inhibition of transcription if global miRNA levels are reduced. Collectively, our results provide experimental evidence for the idea that miRNAs act cumulatively to stabilize the transcriptional regulatory network. We therefore found a promising approach to assess the effect of these molecules on transcription dynamics.

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.

Evolution and Plasticity of the Transcriptome Under Temperature Fluctuations in the Fungal Plant Pathogen Zymoseptoria tritici

Most species live in a variable environment in nature. Yet understanding the evolutionary processes underlying molecular adaptation to fluctuations remains a challenge. In this study we investigate the transcriptome of the fungal wheat pathogen Zymoseptoria tritici after experimental evolution under stable or fluctuating temperature, by comparing ancestral and evolved populations simultaneously. We found that temperature regimes could have a large and pervasive effect on the transcriptome evolution, with as much as 38% of the genes being differentially expressed between selection regimes. Although evolved lineages showed different changes of gene expression based on ancestral genotypes, we identified a set of genes responding specifically to fluctuation. We found that transcriptome evolution in fluctuating conditions was repeatable between parallel lineages initiated from the same genotype for about 60% of the differentially expressed genes. Further, we detected several hotspots of significantly differentially expressed genes in the genome, in regions known to be enriched in repetitive elements, including accessory chromosomes. Our findings also evidenced gene expression evolution toward a gain of robustness (loss of phenotypic plasticity) associated with the fluctuating regime, suggesting robustness is adaptive in changing environment. This work provides valuable insight into the role of transcriptional rewiring for rapid adaptation to abiotic changes in filamentous plant pathogens.

[Phenotypic plasticity in insects]

Insects represent 85% of the animals. They have adapted to many environments and play a major role in ecosystems. Many insect species exhibit phenotypic plasticity. We here report on the mechanisms involved in phenotypic plasticity of different insects (aphids, migratory locust, map butterfly, honeybee) and also on the nutritional size plasticity in Drosophila and the plasticity of the wing eye-spots of the butterfly Bicyclus anynana. We also describe in more detail our work concerning the thermal plasticity of pigmentation in Drosophila. We have shown that the expression of the tan, yellow and Ddc genes, encoding enzymes of the melanin synthesis pathway, is modulated by temperature and that it is a consequence, at least in part, of the temperature-sensitive expression of the bab locus genes that repress them.

[Phenotypic plasticity: a brief introduction]

Phenotypic plasticity describes the ability of a given genotype to produce different phenotypes in response to distinct environmental conditions. It has major implications in agronomy, animal husbandry and medicine and is also thought to facilitate evolution. Phenotypic plasticity is widely observed in the wild. It is only relatively recently that the mechanisms involved in phenotypic plasticity have been analysed. Thanks to laboratory experiments we understand better how environmental conditions are involved in phenotypic variations. This article introduces major concepts from the phenotypic plasticity field, presents briefly mechanisms involved in phenotypic plasticity and discusses the links between phenotypic plasticity and evolution.

Integrating GWAS and Transcriptomics to Identify the Molecular Underpinnings of Thermal Stress Responses in Drosophila melanogaster

Thermal tolerance of an organism depends on both the ability to dynamically adjust to a thermal stress and preparatory developmental processes that enhance thermal resistance. However, the extent to which standing genetic variation in thermal tolerance alleles influence dynamic stress responses vs. preparatory processes is unknown. Here, using the model species Drosophila melanogaster, we used a combination of Genome Wide Association mapping (GWAS) and transcriptomic profiling to characterize whether genes associated with thermal tolerance are primarily involved in dynamic stress responses or preparatory processes that influence physiological condition at the time of thermal stress. To test our hypotheses, we measured the critical thermal minimum (CTmin) and critical thermal maximum (CTmax) of 100 lines of the Drosophila Genetic Reference Panel (DGRP) and used GWAS to identify loci that explain variation in thermal limits. We observed greater variation in lower thermal limits, with CTmin ranging from 1.81 to 8.60°C, while CTmax ranged from 38.74 to 40.64°C. We identified 151 and 99 distinct genes associated with CTmin and CTmax, respectively, and there was strong support that these genes are involved in both dynamic responses to thermal stress and preparatory processes that increase thermal resistance. Many of the genes identified by GWAS were involved in the direct transcriptional response to thermal stress (72/151 for cold; 59/99 for heat), and overall GWAS candidates were more likely to be differentially expressed than other genes. Further, several GWAS candidates were regulatory genes that may participate in the regulation of stress responses, and gene ontologies related to development and morphogenesis were enriched, suggesting many of these genes influence thermal tolerance through effects on development and physiological status. Overall, our results suggest that thermal tolerance alleles can influence both dynamic plastic responses to thermal stress and preparatory processes that improve thermal resistance. These results also have utility for directly comparing GWAS and transcriptomic approaches for identifying candidate genes associated with thermal tolerance.

Different genetic basis for alcohol dehydrogenase activity and plasticity in a novel alcohol environment for Drosophila melanogaster

Phenotypic plasticity is known to enhance population persistence, facilitate adaptive evolution and initiate novel phenotypes in novel environments. How plasticity can contribute or hinder adaptation to different environments hinges on its genetic architecture. Even though plasticity in many traits is genetically controlled, whether and how plasticity's genetic architecture might change in novel environments is still unclear. Because much of gene expression can be environmentally influenced, each environment may trigger different sets of genes that influence a trait. Using a quantitative trait loci (QTL) approach, we investigated the genetic basis of plasticity in a classic functional trait, alcohol dehydrogenase (ADH) activity in D. melanogaster, across both historical and novel alcohol environments. Previous research in D. melanogaster has also demonstrated that ADH activity is plastic in response to alcohol concentration in substrates used by both adult flies and larvae. We found that across all environments tested, ADH activity was largely influenced by a single QTL encompassing the Adh-coding gene and its known regulatory locus, delta-1. After controlling for the allelic variation of the Adh and delta-1 loci, we found additional but different minor QTLs in the 0 and 14% alcohol environments. In contrast, we discovered no major QTL for plasticity itself, including the Adh locus, regardless of the environmental gradients. This suggests that plasticity in ADH activity is likely influenced by many loci with small effects, and that the Adh locus is not environmentally sensitive to dietary alcohol.

Life-History Evolution and the Genetics of Fitness Components in Drosophila melanogaster

Life-history traits or "fitness components"-such as age and size at maturity, fecundity and fertility, age-specific rates of survival, and life span-are the major phenotypic determinants of Darwinian fitness. Analyzing the evolution and genetics of these phenotypic targets of selection is central to our understanding of adaptation. Due to its simple and rapid life cycle, cosmopolitan distribution, ease of maintenance in the laboratory, well-understood evolutionary genetics, and its versatile genetic toolbox, the "vinegar fly" Drosophila melanogaster is one of the most powerful, experimentally tractable model systems for studying "life-history evolution." Here, I review what has been learned about the evolution and genetics of life-history variation in D. melanogaster by drawing on numerous sources spanning population and quantitative genetics, genomics, experimental evolution, evolutionary ecology, and physiology. This body of work has contributed greatly to our knowledge of several fundamental problems in evolutionary biology, including the amount and maintenance of genetic variation, the evolution of body size, clines and climate adaptation, the evolution of senescence, phenotypic plasticity, the nature of life-history trade-offs, and so forth. While major progress has been made, important facets of these and other questions remain open, and the D. melanogaster system will undoubtedly continue to deliver key insights into central issues of life-history evolution and the genetics of adaptation.

Evolution of Epistatic Networks and the Genetic Basis of Innate Behaviors

Instinctive behaviors are genetically programmed behaviors that occur independent of experience. How genetic programs that give rise to the manifestation of such behaviors evolve remains an unresolved question. I propose that evolution of species-specific innate behaviors is accomplished through progressive modifications of pre-existing genetic networks composed of allelic variants. I hypothesize that changes in frequencies of one or more constituent allelic variants within the network leads to changes in gene network connectivity and the emergence of a reorganized network that can support the emergence of a novel behavioral phenotype and becomes stabilized when key allelic variants are driven to fixation.

Genomics of Developmental Plasticity in Animals

Developmental plasticity refers to the property by which the same genotype produces distinct phenotypes depending on the environmental conditions under which development takes place. By allowing organisms to produce phenotypes adjusted to the conditions that adults will experience, developmental plasticity can provide the means to cope with environmental heterogeneity. Developmental plasticity can be adaptive and its evolution can be shaped by natural selection. It has also been suggested that developmental plasticity can facilitate adaptation and promote diversification. Here, we summarize current knowledge on the evolution of plasticity and on the impact of plasticity on adaptive evolution, and we identify recent advances and important open questions about the genomics of developmental plasticity in animals. We give special attention to studies using transcriptomics to identify genes whose expression changes across developmental environments and studies using genetic mapping to identify loci that contribute to variation in plasticity and can fuel its evolution.

Sex allocation plasticity on a transcriptome scale: Socially sensitive gene expression in a simultaneous hermaphrodite

Phenotypic plasticity can enable organisms to produce optimal phenotypes in multiple environments. A crucial life history trait that is often highly plastic is sex allocation, which in simultaneous hermaphrodites describes the relative investment into the male versus female sex functions. Theory predicts-and morphological evidence supports-that greater investment into the male function is favoured with increasing group size, due to the increasing importance of sperm competition for male reproductive success. Here, we performed a genome-wide gene expression assay to test for such sex allocation plasticity in a model simultaneous hermaphrodite, the free-living flatworm Macrostomum lignano. Based on RNA-Seq data from 16 biological replicates spanning four different group size treatments, we demonstrate that at least 10% of the >75,000 investigated transcripts in M. lignano are differentially expressed according to the social environment, rising to >30% of putative gonad-specific transcripts (spermatogenesis and oogenesis candidates) and tail-specific transcripts (seminal fluid candidates). This transcriptional response closely corresponds to the expected shift away from female and towards male reproductive investment with increasing sperm competition level. Using whole-mount in situ hybridization, we then confirm that many plastic transcripts exhibit the expected organ-specific expression, and RNA interference of selected testis- and ovary-specific candidates establishes that these indeed function in gametogenesis pathways. We conclude that a large proportion of sex-specific transcripts in M. lignano are differentially expressed according to the prevailing ecological conditions and that these are functionally relevant to key reproductive phenotypes. Our study thus begins to bridge organismal and molecular perspectives on sex allocation plasticity.

Effects of adult temperature on gene expression in a butterfly: identifying pathways associated with thermal acclimation

BACKGROUND

Phenotypic plasticity is a pervasive property of all organisms and considered to be of key importance for dealing with environmental variation. Plastic responses to temperature, which is one of the most important ecological factors, have received much attention over recent decades. A recurrent pattern of temperature-induced adaptive plasticity includes increased heat tolerance after exposure to warmer temperatures and increased cold tolerance after exposure to cooler temperatures. However, the mechanisms underlying these plastic responses are hitherto not well understood. Therefore, we here investigate effects of adult acclimation on gene expression in the tropical butterfly Bicyclus anynana, using an RNAseq approach.

RESULTS

We show that several antioxidant markers (e.g. peroxidase, cytochrome P450) were up-regulated at a higher temperature compared with a lower adult temperature, which might play an important role in the acclamatory responses subsequently providing increased heat tolerance. Furthermore, several metabolic pathways were up-regulated at the higher temperature, likely reflecting increased metabolic rates. In contrast, we found no evidence for a decisive role of the heat shock response.

CONCLUSIONS

Although the important role of antioxidant defence mechanisms in alleviating detrimental effects of oxidative stress is firmly established, we speculate that its potentially important role in mediating heat tolerance and survival under stress has been underestimated thus far and thus deserves more attention.

Symbiont type and environmental factors affect transcriptome-wide gene expression in the coral Montipora capitata

Reef-building corals may harbor genetically distinct lineages of endosymbiotic dinoflagellates in the genus Symbiodinium, which have been shown to affect important colony properties, including growth rates and resilience against environmental stress. However, the molecular processes underlying these differences are not well understood. In this study, we used whole transcriptome sequencing (RNA-seq) to assess gene expression differences between 27 samples of the coral Montipora capitata predominantly hosting two different Symbiodinium types in clades C and D. The samples were further characterized by their origin from two field sites on Hawai'i Island with contrasting environmental conditions. We found that transcriptome-wide gene expression profiles clearly separated by field site first, and symbiont clade second. With 273 differentially expressed genes (DEGs, 1.3% of all host transcripts), symbiont clade had a measurable effect on host gene expression, but the effect of field site proved almost an order of magnitude higher (1,957 DEGs, 9.6%). According to SNP analysis, we found moderate evidence for host genetic differentiation between field sites (F ST = 0.046) and among corals harboring alternative symbiont clades (F ST = 0.036), suggesting that site-related gene expression differences are likely due to a combination of local adaptation and acclimatization to environmental factors. The correlation between host gene expression and symbiont clade may be due to several factors, including host genotype or microhabitat selecting for alternative clades, host physiology responding to different symbionts, or direct modulation of host gene expression by Symbiodinium. However, the magnitude of these effects at the level of transcription was unexpectedly small considering the contribution of symbiont type to holobiont phenotype.

Effects of mutation and selection on plasticity of a promoter activity in Saccharomyces cerevisiae

Phenotypic plasticity is an evolvable property of biological systems that can arise from environment-specific regulation of gene expression. To better understand the evolutionary and molecular mechanisms that give rise to plasticity in gene expression, we quantified the effects of 235 single-nucleotide mutations in the Saccharomyces cerevisiae TDH3 promoter (P_TDH3 ) on the activity of this promoter in media containing glucose, galactose, or glycerol as a carbon source. We found that the distributions of mutational effects differed among environments because many mutations altered the plastic response exhibited by the wild-type allele. Comparing the effects of these mutations with the effects of 30 P_TDH3 polymorphisms on expression plasticity in the same environments provided evidence of natural selection acting to prevent the plastic response in P_TDH3 activity between glucose and galactose from becoming larger. The largest changes in expression plasticity were observed between fermentable (glucose or galactose) and nonfermentable (glycerol) carbon sources and were caused by mutations located in the RAP1 and GCR1 transcription factor binding sites. Mutations altered expression plasticity most frequently between the two fermentable environments, with mutations causing significant changes in plasticity between glucose and galactose distributed throughout the promoter, suggesting they might affect chromatin structure. Taken together, these results provide insight into the molecular mechanisms underlying gene-by-environment interactions affecting gene expression as well as the evolutionary dynamics affecting natural variation in plasticity of gene expression.

Mating Changes Sexually Dimorphic Gene Expression in the Seed Beetle Callosobruchus maculatus

Sexually dimorphic phenotypes arise largely from sex-specific gene expression, which has mainly been characterized in sexually naïve adults. However, we expect sexual dimorphism in transcription to be dynamic and dependent on factors such as reproductive status. Mating induces many behavioral and physiological changes distinct to each sex and is therefore expected to activate regulatory changes in many sex-biased genes. Here, we first characterized sexual dimorphism in gene expression in Callosobruchus maculatus seed beetles. We then examined how females and males respond to mating and how it affects sex-biased expression, both in sex-limited (abdomen) and sex-shared (head and thorax) tissues. Mating responses were largely sex-specific and, as expected, females showed more genes responding compared with males (∼2,000 vs. ∼300 genes in the abdomen, ∼500 vs. ∼400 in the head and thorax, respectively). Of the sex-biased genes present in virgins, 16% (1,041 genes) in the abdomen and 17% (243 genes) in the head and thorax altered their relative expression between the sexes as a result of mating. Sex-bias status changed in 2% of the genes in the abdomen and 4% in the head and thorax following mating. Mating responses involved de-feminization of females and, to a lesser extent, de-masculinization of males relative to their virgin state: mating decreased rather than increased dimorphic expression of sex-biased genes. The fact that regulatory changes of both types of sex-biased genes occurred in both sexes suggests that male- and female-specific selection is not restricted to male- and female-biased genes, respectively, as is sometimes assumed.

Molecular Ecological Basis of Grasshopper (Oedaleus asiaticus) Phenotypic Plasticity under Environmental Selection

While ecological adaptation in insects can be reflected by plasticity of phenotype, determining the causes and molecular mechanisms for phenotypic plasticity (PP) remains a crucial and still difficult question in ecology, especially where control of insect pests is involved. Oedaleus asiaticus is one of the most dominant pests in the Inner Mongolia steppe and represents an excellent system to study phenotypic plasticity. To better understand ecological factors affecting grasshopper phenotypic plasticity and its molecular control, we conducted a full transcriptional screening of O. asiaticus grasshoppers reared in four different grassland patches in Inner Mongolia. Grasshoppers showed different degrees of PP associated with unique gene expressions and different habitat plant community compositions. Grasshopper performance variables were susceptible to habitat environment conditions and closely associated with plant architectures. Intriguingly, eco-transcriptome analysis revealed five potential candidate genes playing important roles in grasshopper performance, with gene expression closely relating to PP and plant community factors. By linking the grasshopper performances to gene profiles and ecological factors using canonical regression, we first demonstrated the eco-transcriptomic architecture (ETA) of grasshopper phenotypic traits (ETAGPTs). ETAGPTs revealed plant food type, plant density, coverage, and height were the main ecological factors influencing PP, while insect cuticle protein (ICP), negative elongation factor A (NELFA), and lactase-phlorizin hydrolase (LCT) were the key genes associated with PP. Our study gives a clear picture of gene-environment interaction in the formation and maintenance of PP and enriches our understanding of the transcriptional events underlying molecular control of rapid phenotypic plasticity associated with environmental variability. The findings of this study may also provide new targets for pest control and highlight the significance of ecological management practice on grassland conservation.

Incorporation of Biological Knowledge Into the Study of Gene-Environment Interactions

A growing knowledge base of genetic and environmental information has greatly enabled the study of disease risk factors. However, the computational complexity and statistical burden of testing all variants by all environments has required novel study designs and hypothesis-driven approaches. We discuss how incorporating biological knowledge from model organisms, functional genomics, and integrative approaches can empower the discovery of novel gene-environment interactions and discuss specific methodological considerations with each approach. We consider specific examples where the application of these approaches has uncovered effects of gene-environment interactions relevant to drug response and immunity, and we highlight how such improvements enable a greater understanding of the pathogenesis of disease and the realization of precision medicine.

The quantitative genetic basis of clinal divergence in phenotypic plasticity

Phenotypic plasticity is thought to be an important mechanism for adapting to environmental heterogeneity. Nonetheless, the genetic basis of plasticity is still not well understood. In Drosophila melanogaster and D. simulans, body size and thermal stress resistance show clinal patterns along the east coast of Australia, and exhibit plastic responses to different developmental temperatures. The genetic basis of thermal plasticity, and whether the genetic effects underlying clinal variation in traits and their plasticity are similar, remains unknown. Here, we use line-cross analyses between a tropical and temperate population of Drosophila melanogaster and D. simulans developed at three constant temperatures (18°C, 25°C, and 29°C) to investigate the quantitative genetic basis of clinal divergence in mean thermal response (elevation) and plasticity (slope and curvature) for thermal stress and body size traits. Generally, the genetic effects underlying divergence in mean response and plasticity differed, suggesting that different genetic models may be required to understand the evolution of trait means and plasticity. Furthermore, our results suggest that nonadditive genetic effects, in particular epistasis, may commonly underlie plastic responses, indicating that current models that ignore epistasis may be insufficient to understand and predict evolutionary responses to environmental change.

Genomic responses to the socio-sexual environment in male Drosophila melanogaster exposed to conspecific rivals

Socio-sexual environments have profound effects on fitness. Local sex ratios can alter the threat of sexual competition, to which males respond via plasticity in reproductive behaviors and ejaculate composition. In Drosophila melanogaster, males detect the presence of conspecific, same-sex mating rivals prior to mating using multiple, redundant sensory cues. Males that respond to rivals gain significant fitness benefits by altering mating duration and ejaculate composition. Here we investigated the underlying genome-wide changes involved. We used RNA-seq to analyze male transcriptomic responses 2, 26, and 50 h after exposure to rivals, a time period that was previously identified as encompassing the major facets of male responses to rivals. The results showed a strong early activation of multiple sensory genes in the head-thorax (HT), prior to the expression of any phenotypic differences. This gene expression response was reduced by 26 h, at the time of maximum phenotypic change, and shut off by 50 h. In the abdomen (A), fewer genes changed in expression and gene expression responses appeared to increase over time. The results also suggested that different sets of functionally equivalent genes might be activated in different replicates. This could represent a mechanism by which robustness is conferred upon highly plastic traits. Overall, our study reveals that mRNA-seq can identify subtle genomic signatures characteristic of flexible behavioral phenotypes.

Response of Gene Expression and Alternative Splicing to Distinct Growth Environments in Tomato

Phenotypic plasticity is the phenomenon that one particular genotype produces different phenotypes under different environmental conditions, but its underlying molecular and genetic mechanisms are poorly understood. Plastic traits may be under the control of genes whose expression is modulated by environmental cues. In this study, we investigated phenotypic plasticity in tomato (Solanum lycopersicum) and its ancestral species S. pimpinellifolium by comparing the global gene expression of young seedlings grown under two distinct growth conditions. Our results show that more than 7000 genes exhibited differential expression in response to environmental changes from phytotron to a plastic greenhouse, and 98 environmentally sensitive genes displayed the same patterns of expression response across the two tomato species. We also found that growth conditions had a remarkable impact on transcriptome complexity, attributable to alternative splicing (AS), in which 665 splice variants showed differential expression in response to the environmental changes. Moreover, more splice variants and AS events per gene were detected in plastic greenhouse-grown seedlings than their phytotron counterparts, and these seedlings also had higher percentages of intron retention events. The identification of the conserved environmentally-sensitive genes and the splice variants in this study will be useful for further analysis of gene regulation of environmental response in tomato and other crops.

Modulation of yellow expression contributes to thermal plasticity of female abdominal pigmentation in Drosophila melanogaster

Phenotypic plasticity describes the ability of a given genotype to produce distinct phenotypes in different environments. We use the temperature sensitivity of abdominal pigmentation in Drosophila melanogaster females as a model to analyse the effect of the environment on development. We reported previously that thermal plasticity of abdominal pigmentation in females involves the pigmentation gene tan (t). However, the expression of the pigmentation gene yellow (y) was also modulated by temperature in the abdominal epidermis of pharate females. We investigate here the contribution of y to female abdominal pigmentation plasticity. First, we show that y is required for the production of black Dopamine-melanin. Then, using in situ hybridization, we show that the expression of y is strongly modulated by temperature in the abdominal epidermis of pharate females but not in bristles. Interestingly, these two expression patterns are known to be controlled by distinct enhancers. However, the activity of the y-wing-body epidermal enhancer only partially mediates the effect of temperature suggesting that additional regulatory sequences are involved. In addition, we show that y and t co-expression is needed to induce strong black pigmentation indicating that y contributes to female abdominal pigmentation plasticity.

Effects of temperature on transcriptome and cuticular hydrocarbon expression in ecologically differentiated populations of desert Drosophila

We assessed the effects of temperature differences on gene expression using whole-transcriptome microarrays and cuticular hydrocarbon variation in populations of cactophilic Drosophila mojavensis. Four populations from Baja California and mainland Mexico and Arizona were each reared on two different host cacti, reared to sexual maturity on laboratory media, and adults were exposed for 12 hr to 15, 25, or 35°C. Temperature differences influenced the expression of 3,294 genes, while population differences and host plants affected >2,400 each in adult flies. Enriched, functionally related groups of genes whose expression changed at high temperatures included heat response genes, as well as genes affecting chromatin structure. Gene expression differences between mainland and peninsular populations included genes involved in metabolism of secondary compounds, mitochondrial activity, and tRNA synthases. Flies reared on the ancestral host plant, pitaya agria cactus, showed upregulation of genes involved in metabolism, while flies reared on organ pipe cactus had higher expression of DNA repair and chromatin remodeling genes. Population × environment (G × E) interactions had widespread effects on the transcriptome where population × temperature interactions affected the expression of >5,000 orthologs, and there were >4,000 orthologs that showed temperature × host plant interactions. Adults exposed to 35°C had lower amounts of most cuticular hydrocarbons than those exposed to 15 or 25°C, including abundant unsaturated alkadienes. For insects adapted to different host plants and climatic regimes, our results suggest that temperature shifts associated with climate change have large and significant effects on transcriptomes of genetically differentiated natural populations.

Microenvironmental Gene Expression Plasticity Among Individual Drosophila melanogaster

Differences in phenotype among genetically identical individuals exposed to the same environmental condition are often noted in genetic studies. Despite this commonplace observation, little is known about the causes of this variability, which has been termed microenvironmental plasticity. One possibility is that stochastic or technical sources of variance produce these differences. A second possibility is that this variation has a genetic component. We have explored gene expression robustness in the transcriptomes of 730 individual Drosophila melanogaster of 16 fixed genotypes, nine of which are infected with Wolbachia. Three replicates of flies were grown, controlling for food, day/night cycles, humidity, temperature, sex, mating status, social exposure, and circadian timing of RNA extraction. Despite the use of inbred genotypes, and carefully controlled experimental conditions, thousands of genes were differentially expressed, revealing a unique and dynamic transcriptional signature for each individual fly. We found that 23% of the transcriptome was differentially expressed among individuals, and that the variability in gene expression among individuals is influenced by genotype. This transcriptional variation originated from specific gene pathways, suggesting a plastic response to the microenvironment; but there was also evidence of gene expression differences due to stochastic fluctuations. These observations reveal previously unappreciated genetic sources of variability in gene expression among individuals, which has implications for complex trait genetics and precision medicine.

Thermal plasticity in Drosophila melanogaster populations from eastern Australia: quantitative traits to transcripts

The flexibility afforded to genotypes in different environments by phenotypic plasticity is of interest to biologists studying thermal adaptation because of the thermal lability of many traits. Differences in thermal performance and reaction norms can provide insight into the evolution of thermal adaptation to explore broader questions such as species distributions and persistence under climate change. One approach is to study the effects of temperature on fitness, morphological and more recently gene expression traits in populations from different climatic origins. The diverse climatic conditions experienced by Drosophila melanogaster along the eastern Australian temperate-tropical gradient are ideal given the high degree of continuous trait differentiation, but reaction norm variation has not been well studied in this system. Here, we reared a tropical and temperate population from the ends of the gradient over six developmental temperatures and examined reaction norm variation for five quantitative traits including thermal performance for fecundity, and reaction norms for thermotolerance, body size, viability and 23 transcript-level traits. Despite genetic variation for some quantitative traits, we found no differentiation between the populations for fecundity thermal optima and breadth, and the reaction norms for the other traits were largely parallel, supporting previous work suggesting that thermal evolution occurs by changes in trait means rather than by reaction norm shifts. We examined reaction norm variation in our expanded thermal regime for a gene set shown to previously exhibit GxE for expression plasticity in east Australian flies, as well as key heat-shock genes. Although there were differences in curvature between the populations suggesting a higher degree of thermal plasticity in expression patterns than for the quantitative traits, we found little evidence to support a role for genetic variation in maintaining expression plasticity.

The Interplay of Temperature and Genotype on Patterns of Alternative Splicing in Drosophila melanogaster

Alternative splicing is the highly regulated process of variation in the removal of introns from premessenger-RNA transcripts. The consequences of alternative splicing on the phenotype are well documented, but the impact of the environment on alternative splicing is not yet clear. We studied variation in alternative splicing among four different temperatures, 13, 18, 23, and 29°, in two Drosophila melanogaster genotypes. We show plasticity of alternative splicing with up to 10% of the expressed genes being differentially spliced between the most extreme temperatures for a given genotype. Comparing the two genotypes at different temperatures, we found <1% of the genes being differentially spliced at 18°. At extreme temperatures, however, we detected substantial differences in alternative splicing-with almost 10% of the genes having differential splicing between the genotypes: a magnitude similar to between species differences. Genes with differential alternative splicing between genotypes frequently exhibit dominant inheritance. Remarkably, the pattern of surplus of differences in alternative splicing at extreme temperatures resembled the pattern seen for gene expression intensity. Since different sets of genes were involved for the two phenotypes, we propose that purifying selection results in the reduction of differences at benign temperatures. Relaxed purifying selection at temperature extremes, on the other hand, may cause the divergence in gene expression and alternative splicing between the two strains in rarely encountered environments.

Phenotypic Plasticity through Transcriptional Regulation of the Evolutionary Hotspot Gene tan in Drosophila melanogaster

Phenotypic plasticity is the ability of a given genotype to produce different phenotypes in response to distinct environmental conditions. Phenotypic plasticity can be adaptive. Furthermore, it is thought to facilitate evolution. Although phenotypic plasticity is a widespread phenomenon, its molecular mechanisms are only beginning to be unravelled. Environmental conditions can affect gene expression through modification of chromatin structure, mainly via histone modifications, nucleosome remodelling or DNA methylation, suggesting that phenotypic plasticity might partly be due to chromatin plasticity. As a model of phenotypic plasticity, we study abdominal pigmentation of Drosophila melanogaster females, which is temperature sensitive. Abdominal pigmentation is indeed darker in females grown at 18°C than at 29°C. This phenomenon is thought to be adaptive as the dark pigmentation produced at lower temperature increases body temperature. We show here that temperature modulates the expression of tan (t), a pigmentation gene involved in melanin production. t is expressed 7 times more at 18°C than at 29°C in female abdominal epidermis. Genetic experiments show that modulation of t expression by temperature is essential for female abdominal pigmentation plasticity. Temperature modulates the activity of an enhancer of t without modifying compaction of its chromatin or level of the active histone mark H3K27ac. By contrast, the active mark H3K4me3 on the t promoter is strongly modulated by temperature. The H3K4 methyl-transferase involved in this process is likely Trithorax, as we show that it regulates t expression and the H3K4me3 level on the t promoter and also participates in female pigmentation and its plasticity. Interestingly, t was previously shown to be involved in inter-individual variation of female abdominal pigmentation in Drosophila melanogaster, and in abdominal pigmentation divergence between Drosophila species. Sensitivity of t expression to environmental conditions might therefore give more substrate for selection, explaining why this gene has frequently been involved in evolution of pigmentation.

Environmental conditions can strongly modulate the phenotype produced by a particular genotype. This process, called phenotypic plasticity, has major implications in medicine and agricultural sciences, and is thought to facilitate evolution. Phenotypic plasticity is observed in many animals and plants but its mechanisms are only partially understood. As a model of phenotypic plasticity, we study the effect of temperature on female abdominal pigmentation in the fruit fly Drosophila melanogaster. Here we show that temperature affects female abdominal pigmentation by modulating the expression of tan (t), a gene involved in melanin production, in female abdominal epidermis. This effect is mediated at least partly by a particular regulatory sequence of t, the t_MSE enhancer. However we detected no modulation of chromatin structure of t_MSE by temperature. By contrast, the level of the active chromatin mark H3K4me3 on the t promoter is strongly increased at lower temperature. We show that the H3K4 methyl-transferase Trithorax is involved in female abdominal pigmentation and its plasticity and regulates t expression and H3K4me3 level on the t promoter. Several studies have linked t to pigmentation evolution within and between Drosophila species. Our results suggest that sensitivity of t expression to temperature might facilitate its role in pigmentation evolution.

Spontaneous mutations and the origin and maintenance of quantitative genetic variation

Mutation and natural selection shape the genetic variation in natural populations. Here, we directly estimated the spontaneous mutation rate by sequencing new Drosophila mutation accumulation lines maintained with minimal natural selection. We inferred strong stabilizing natural selection on quantitative traits because genetic variation among wild-derived inbred lines was much lower than predicted from a neutral model and the mutational effects were much larger than allelic effects of standing polymorphisms. Stabilizing selection could act directly on the traits, or indirectly from pleiotropic effects on fitness. However, our data are not consistent with simple models of mutation-stabilizing selection balance; therefore, further empirical work is needed to assess the balance of evolutionary forces responsible for quantitative genetic variation.

DOI:http://dx.doi.org/10.7554/eLife.14625.001

A key challenge in evolutionary biology is to understand how genetic variation – differences in the DNA of individuals in a population – is generated and maintained to create the enormous diversity that exists in nature. Mutations to the DNA introduce new variation, but these are constantly removed from populations by two other evolutionary forces: natural selection and genetic drift.

Natural selection removes harmful genetic mutations that affect an organism’s fitness and reproduction, and genetic drift is the random increase in, or loss of, a genetic variant from a population over time. However, disentangling the effects of these evolutionary forces is challenging because the genetic variation we observe is often the final product of a long history of interaction between them.

Huang et al. have now investigated genetic variation by breeding fruit flies in the laboratory. Natural selection was minimized for these flies; genetic drift was therefore the main force that removed variation.

Huang et al. then sequenced the DNA of the flies to estimate the rate at which genetic mutations spontaneously occur. The sequences contained many more “high-impact” mutations (which directly affect how proteins in the fly’s cells work) than seen in sequences taken from a natural fly population.

Traits that are produced by the cumulative actions of many genes and the environment are known as quantitative traits. By examining how much variation genetic mutations introduced into the quantitative traits of each generation of the laboratory-grown flies, Huang et al. estimated how much variation should occur in a natural population whose quantitative traits evolved without natural selection. This estimate was much higher than the levels of genetic variation seen in nature, suggesting that natural selection acts to eliminate mutations that significantly affect quantitative traits.

Simple theoretical models cannot explain the relatively high spontaneous mutation rate and low genetic variation seen in the quantitative traits of natural populations. Therefore, further work is now required to understand more about the balance of evolutionary forces that maintain quantitative genetic variation.

DOI:http://dx.doi.org/10.7554/eLife.14625.002

Mitonuclear Epistasis for Development Time and Its Modification by Diet in Drosophila

Mitochondrial (mtDNA) and nuclear genes have to operate in a coordinated manner to maintain organismal function, and the regulation of this homeostasis presents a substantial source of potential epistatic (G × G) interactions. How these interactions shape the fitness landscape is poorly understood. Here we developed a novel mitonuclear epistasis model, using selected strains of the Drosophila Genetic Reference Panel (DGRP) and mitochondrial genomes from within Drosophila melanogaster and D. simulans to test the hypothesis that mtDNA × nDNA interactions influence fitness. In total we built 72 genotypes (12 nuclear backgrounds × 6 mtDNA haplotypes, with 3 from each species) to dissect the relationship between genotype and phenotype. Each genotype was assayed on four food environments. We found considerable variation in several phenotypes, including development time and egg-to-adult viability, and this variation was partitioned into genetic (G), environmental (E), and higher-order (G × G, G × E, and G × G × E) components. Food type had a significant impact on development time and also modified mitonuclear epistases, evidencing a broad spectrum of G × G × E across these genotypes. Nuclear background effects were substantial, followed by mtDNA effects and their G × G interaction. The species of mtDNA haplotype had negligible effects on phenotypic variation and there was no evidence that mtDNA variation has different effects on male and female fitness traits. Our results demonstrate that mitonuclear epistases are context dependent, suggesting the selective pressure acting on mitonuclear genotypes may vary with food environment in a genotype-specific manner.