Canalization of gene expression is a major signature of regulatory cold adaptation in temperate Drosophila melanogaster

Navigating the Cold: Integrative Transcriptome Sequencing Approach Reveals Ionoregulatory and Whole-Body Responses to Cold Acclimation in Drosophila ananassae

Understanding how species adapt to changing environments is a major goal in evolutionary biology and can elucidate the impact of climate change. Climate imposes inevitable effects on the geographical distribution of insects as their body temperature primarily depends on the environment. The vinegar fly Drosophila ananassae expanded from its tropical ancestral range to more temperate regions, which requires adaptation to colder climates. Transcriptome and genome-wide association studies focusing on the ancestral-range population identified the targets of selection related to ionoregulatory tissues. However, how cosmopolitan D. ananassae adapted to colder environments, where low temperatures last longer, is still unknown. Here, we present a study on the effect of long-term cold exposure on D. ananassae, examining the gene expression variation in the whole body and the ionoregulatory tissues, namely the hindgut and the Malpighian tubule. To elucidate molecular mechanisms of cold adaptation during species expansion, we included cold-tolerant and cold-sensitive strains from the ancestral species range and cold-tolerant strains from the derived species range. We show that cold acclimation improves cold tolerance and results in differential expression of more than half of the transcriptome in the ionoregulatory tissues and the whole body. Notably, we provide complementary insight into molecular processes at four levels: strains, populations, phenotypes, and tissues. By determining the biochemical pathways of phenotypic plasticity underlying cold tolerance, our results enhance our understanding of how environmental changes affect thermal adaptation in natural populations.

Removal of Exogenous Stimuli Reveals a Canalization of Circadian Physiology in a Vertically Migrating Copepod

Diel rhythms are observed across taxa and are important for maintaining synchrony between the environment and organismal physiology. A striking example of this is the diel vertical migration undertaken by zooplankton, some of which, such as the 5 mm-long copepod Pleuromamma xiphias (P. xiphias), migrate hundreds of meters daily between the surface ocean and deeper waters. Some of the molecular pathways that underlie the expressed phenotype at different stages of this migration are entrained by environmental variables (e.g., day length and food availability), while others are regulated by internal clocks. We identified a series of proteomic biomarkers that vary across ocean DVM and applied them to copepods incubated in 24 h of darkness to assess circadian control. The dark-incubated copepods shared some proteomic similarities to the ocean-caught copepods (i.e., increased abundance of carbohydrate metabolism proteins at night). Shipboard-incubated copepods demonstrated a clearer distinction between night and day proteomic profiles, and more proteins were differentially abundant than in the in situ copepods, even in the absence of the photoperiod and other environmental cues. This pattern suggests that there is a canalization of rhythmic diel physiology in P. xiphias that reflects likely circadian clock control over diverse molecular pathways.

Evidence that variation in root anatomy contributes to local adaptation in Mexican native maize

Mexican native maize (Zea mays ssp. mays) is adapted to a wide range of climatic and edaphic conditions. Here, we focus specifically on the potential role of root anatomical variation in this adaptation. Given the investment required to characterize root anatomy, we present a machine-learning approach using environmental descriptors to project trait variation from a relatively small training panel onto a larger panel of genotyped and georeferenced Mexican maize accessions. The resulting models defined potential biologically relevant clines across a complex environment that we used subsequently for genotype-environment association. We found evidence of systematic variation in maize root anatomy across Mexico, notably a prevalence of trait combinations favoring a reduction in axial hydraulic conductance in varieties sourced from cooler, drier highland areas. We discuss our results in the context of previously described water-banking strategies and present candidate genes that are associated with both root anatomical and environmental variation. Our strategy is a refinement of standard environmental genome-wide association analysis that is applicable whenever a training set of georeferenced phenotypic data is available.

Temporal regulation of temperature tolerances and gene expression in an arctic insect

Terrestrial arthropods in the Arctic are exposed to highly variable temperatures that frequently reach cold and warm extremes. Yet, ecophysiological studies on arctic insects typically focus on the ability of species to tolerate low temperatures, whereas studies investigating physiological adaptations of species to periodically warm and variable temperatures are few. In this study, we investigated temporal changes in thermal tolerances and the transcriptome in the Greenlandic seed bug Nysius groenlandicus, collected in the field across different times and temperatures in Southern Greenland. We found that plastic changes in heat and cold tolerances occurred rapidly (within hours) and at a daily scale in the field, and that these changes are correlated with diurnal temperature variation. Using RNA sequencing, we provide molecular underpinnings of the rapid adjustments in thermal tolerance across ambient field temperatures and in the laboratory. We show that transcriptional responses are sensitive to daily temperature changes, and days characterized by high temperature variation induced markedly different expression patterns than thermally stable days. Further, genes associated with laboratory-induced heat responses, including expression of heat shock proteins and vitellogenins, were shared across laboratory and field experiments, but induced at time points associated with lower temperatures in the field. Cold stress responses were not manifested at the transcriptomic level.

Tropical super flies: Integrating Cas9 into Drosophila ananassae and its phenotypic effects

Ectotherms such as insects are animals whose body temperature largely depends on ambient temperature and temperature variations provide a selection pressure affecting the geographical distribution of these species. However, over the course of evolution, some insect species managed to colonize environments characterized by various temperature ranges. Therefore, insects provide an excellent study system to investigate the basis of adaptation to temperature changes and extremes. We are generally using the vinegar fly Drosophila ananassae as a model system to investigate the genetic basis of cold tolerance. This species has expanded from its tropical ancestral range to more temperate regions resulting in a cosmopolitan, domestic distribution. Previously, we identified candidate genes significantly associated with cold tolerance in this species. We now established molecular genetic tools to assess the function of these genes. Using CRISPR/Cas9 methodology for genome editing and the PiggyBac system, the Cas9 enzyme was successfully integrated into the genome of three fly strains with different levels of cold tolerance. We further report on preliminary findings that the Cas9 integration itself did not have a consistent effect on tolerance to cold. In conclusion, we offer with our study the molecular tools that allow studying stress-related candidate genes in D. ananassae in the future. In addition, we point out and provide guidance on the challenges that come with genome editing in a non-model species.

Genetic assimilation of ancestral plasticity during parallel adaptation to zinc contamination in Silene uniflora

Phenotypic plasticity in ancestral populations is hypothesized to facilitate adaptation, but evidence is piecemeal and often contradictory. Further, whether ancestral plasticity increases the probability of parallel adaptive changes has not been explored. The most general finding is that ancestral responses to a new environment are reversed following adaptation (known as reversion). We investigated the contribution of ancestral plasticity to adaptive evolution of gene expression in two independently evolved lineages of zinc-tolerant Silene uniflora. We found that the general pattern of reversion is driven by the absence of a widespread stress response in zinc-adapted plants compared with zinc-sensitive plants. We show that ancestral plasticity that moves expression closer to the optimum value in the new environment influences the evolution of gene expression among genes that are likely to be involved in adaptation and increases the chance that genes are recruited repeatedly during adaptation. However, despite convergence in gene expression levels between independently adapted lineages, ancestral plasticity does not influence how similar expression values of adaptive genes become. Surprisingly, we also observed that ancestral plasticity that increases fitness often becomes genetically determined and fixed, that is, genetically assimilated. These results emphasize the important role of ancestral plasticity in parallel adaptation.

Genomic shifts, phenotypic clines and fitness costs associated with cold-tolerance in the Asian tiger mosquito

Climatic variation is a key driver of genetic differentiation and phenotypic traits evolution, and local adaptation to temperature is expected in widespread species. We investigated phenotypic and genomic changes in the native range of the Asian tiger mosquito, Aedes albopictus. We first refine the phylogeographic structure based on genome-wide regions (1,901 ddRAD SNPs) from 41 populations. We then explore the patterns of cold adaptation using phenotypic traits measured in common garden (wing size and cold tolerance) and genotype–temperature associations at targeted candidate regions (51,706 exon capture SNPs) from 9 populations. We confirm the existence of three evolutionary lineages including clades A (Malaysia, Thailand, Cambodia, and Laos), B (China and Okinawa), and C (South Korea and Japan). We identified temperature-associated differentiation in fifteen out of 221 candidate regions but none in ddRAD regions, supporting the role of directional selection in detected genes. These include genes involved in lipid metabolism and a circadian clock gene. Most outlier SNPs are differently fixed between clades A and C, while clade B has an intermediate pattern. Females are larger at higher latitude yet produce no more eggs, which might favor the storage of energetic reserves in colder climate. Non-diapausing eggs from temperate populations survive better to cold exposure than those from tropical populations, suggesting they are protected from freezing damages but this cold tolerance has a fitness cost in terms of egg viability. Altogether, our results provide strong evidence for the thermal adaptation of A. albopictus across its wide temperature range.

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.

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.

Natural variation in the transcriptional response of Drosophila melanogaster to oxidative stress

Broadly distributed species must cope with diverse and changing environmental conditions, including various forms of stress. Cosmopolitan populations of Drosophila melanogaster are more tolerant to oxidative stress than those from the species’ ancestral range in sub-Saharan Africa, and the degree of tolerance is associated with an insertion/deletion polymorphism in the 3′ untranslated region of the Metallothionein A (MtnA) gene that varies clinally in frequency. We examined oxidative stress tolerance and the transcriptional response to oxidative stress in cosmopolitan and sub-Saharan African populations of D. melanogaster, including paired samples with allelic differences at the MtnA locus. We found that the effect of the MtnA polymorphism on oxidative stress tolerance was dependent on the genomic background, with the deletion allele increasing tolerance only in a northern, temperate population. Genes that were differentially expressed under oxidative stress included MtnA and other metallothioneins, as well as those involved in glutathione metabolism and other genes known to be part of the oxidative stress response or the general stress response. A gene coexpression analysis revealed further genes and pathways that respond to oxidative stress including those involved in additional metabolic processes, autophagy, and apoptosis. There was a significant overlap among the genes induced by oxidative and cold stress, which suggests a shared response pathway to these two stresses. Interestingly, the MtnA deletion was associated with consistent changes in the expression of many genes across all genomic backgrounds, regardless of the expression level of the MtnA gene itself. We hypothesize that this is an indirect effect driven by the loss of microRNA binding sites within the MtnA 3′ untranslated region.

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.

Small heat shock protein Hsp67Bc plays a significant role in Drosophila melanogaster cold stress tolerance

Hsp67Bc in Drosophila melanogaster is a member of the small heat shock protein family, the main function of which is to prevent the aggregation of misfolded or damaged proteins. Hsp67Bc interacts with Starvin and Hsp23, which are known to be a part of the cold stress response in the fly during the recovery phase. In this study, we investigated the role of the Hsp67Bc gene in the cold stress response. We showed that in adult Drosophila, Hsp67Bc expression increases after cold stress and decreases after 1.5 h of recovery, indicating the involvement of Hsp67Bc in short-term stress recovery. We also implemented a deletion in the D. melanogaster Hsp67Bc gene using imprecise excision of a P-element, and analysed the cold tolerance of Hsp67Bc-null mutants at different developmental stages. We found that Hsp67Bc-null homozygous flies are viable and fertile but display varying cold stress tolerance throughout the stages of ontogenesis: the survival after cold stress is slightly impaired in late third instar larvae, unaffected in pupae, and notably affected in adult females. Moreover, the recovery from chill coma is delayed in Hsp67Bc-null adults of both sexes. In addition, the deletion in the Hsp67Bc gene caused more prominent up-regulation of Hsp70 following cold stress, suggesting the involvement of Hsp70 in compensation of the lack of the Hsp67Bc protein. Taken together, our results suggest that Hsp67Bc is involved in the recovery of flies from a comatose state and contributes to the protection of the fruit fly from cold stress.

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.

Three Quantitative Trait Loci Explain More than 60% of Variation for Chill Coma Recovery Time in a Natural Population of Drosophila ananassae

Ectothermic species such as insects are particularly vulnerable to climatic fluctuations. Nevertheless, many insects that evolved and diversified in the tropics have successfully colonized temperate regions all over the globe. To shed light on the genetic basis of cold tolerance in such species, we conducted a quantitative trait locus (QTL) mapping experiment for chill coma recovery time (CCRT) in Drosophila ananassae, a cosmopolitan species that has expanded its range from tropical to temperate regions. We created a mapping population of recombinant inbred advanced intercross lines (RIAILs) from two founder strains with diverging CCRT phenotypes. The RIAILs were phenotyped for their CCRT and, together with the founder strains, genotyped for polymorphic markers with double-digest restriction site-associated DNA (ddRAD) sequencing. Using a hierarchical mapping approach that combined standard interval mapping and a multiple-QTL model, we mapped three QTL which altogether explained 64% of the phenotypic variance. For two of the identified QTL, we found evidence of epistasis. To narrow down the list of cold tolerance candidate genes, we cross-referenced the QTL intervals with genes that we previously identified as differentially expressed in response to cold in D. ananassae, and with thermotolerance candidate genes of D. melanogaster. Among the 58 differentially expressed genes that were contained within the QTL, GF15058 showed a significant interaction of the CCRT phenotype and gene expression. Further, we identified the orthologs of four D. melanogaster thermotolerance candidate genes, MtnA, klarsicht, CG5246 (D.ana/GF17132) and CG10383 (D.ana/GF14829) as candidates for cold tolerance in D. ananassae.

Metabolically similar cohorts of bacteria exhibit strong cooccurrence patterns with diet items and eukaryotic microbes in lizard guts

Gut microbiomes perform essential services for their hosts, including helping them to digest food and manage pathogens and parasites. Performing these services requires a diverse and constantly changing set of metabolic functions from the bacteria in the microbiome. The metabolic repertoire of the microbiome is ultimately dependent on the outcomes of the ecological interactions of its member microbes, as these interactions in part determine the taxonomic composition of the microbiome. The ecological processes that underpin the microbiome's ability to handle a variety of metabolic challenges might involve rapid turnover of the gut microbiome in response to new metabolic challenges, or it might entail maintaining sufficient diversity in the microbiome that any new metabolic demands can be met from an existing set of bacteria. To differentiate between these scenarios, we examine the gut bacteria and resident eukaryotes of two generalist-insectivore lizards, while simultaneously identifying the arthropod prey each lizard was digesting at the time of sampling. We find that the cohorts of bacteria that occur significantly more or less often than expected with arthropod diet items or eukaryotes include bacterial species that are highly similar to each other metabolically. This pattern in the bacterial microbiome could represent an early step in the taxonomic shifts in bacterial microbiome that occur when host lineages change their diet niche over evolutionary timescales.

Genetic Basis of Phenotypic Differences Between Chinese Yunling Black Goats and Nubian Goats Revealed by Allele-Specific Expression in Their F1 Hybrids

Chinese Yunling black goats and African Nubian goats are divergent breeds showing significant differences in body size, milk production, and environmental adaptation. However, the genetic mechanisms underlying these phenotypic differences remain to be elucidated. In this report, we provide a detailed portrait of allele-specific expression (ASE) from 54 RNA-Seq analyses across six tissues from nine F1 hybrid offspring generated by crossing the two breeds combined with 13 genomes of the two breeds. We identified a total of 524 genes with ASE, which are involved in bone development, muscle cell differentiation, and the regulation of lipid metabolic processes. We further found that 38 genes with ASE were also under directional selection by comparing 13 genomes of the two breeds; these 38 genes play important roles in metabolism, immune responses, and the adaptation to hot and humid environments. In conclusion, our study shows that the exploration of genes with ASE in F1 hybrids provides an efficient way to understand the genetic basis underlying the phenotypic differences of two diverse goat breeds.

Transcriptome Analysis Reveals Candidate Genes for Cold Tolerance in Drosophila ananassae

Coping with daily and seasonal temperature fluctuations is a key adaptive process for species to colonize temperate regions all over the globe. Over the past 18,000 years, the tropical species Drosophila ananassae expanded its home range from tropical regions in Southeast Asia to more temperate regions. Phenotypic assays of chill coma recovery time (CCRT) together with previously published population genetic data suggest that only a small number of genes underlie improved cold hardiness in the cold-adapted populations. We used high-throughput RNA sequencing to analyze differential gene expression before and after exposure to a cold shock in cold-tolerant lines (those with fast chill coma recovery, CCR) and cold-sensitive lines (slow CCR) from a population originating from Bangkok, Thailand (the ancestral species range). We identified two candidate genes with a significant interaction between cold tolerance and cold shock treatment: GF14647 and GF15058. Further, our data suggest that selection for increased cold tolerance did not operate through the increased activity of heat shock proteins, but more likely through the stabilization of the actin cytoskeleton and a delayed onset of apoptosis.

Mef2 induction of the immediate early gene Hr38/Nr4a is terminated by Sirt1 to promote ethanol tolerance

Drug naïve animals given a single dose of ethanol show changed responses to subsequent doses, including the development of ethanol tolerance and ethanol preference. These simple forms of behavioral plasticity are due in part to changes in gene expression and neuronal properties. Surprisingly little is known about how ethanol initiates changes in gene expression or what the changes do. Here we demonstrate a role in ethanol plasticity for Hr38, the sole Drosophila homolog of the mammalian Nr4a1/2/3 class of immediate early response transcription factors. Acute ethanol exposure induces transient expression of Hr38 and other immediate early neuronal activity genes. Ethanol activates the Mef2 transcriptional activator to induce Hr38, and the Sirt1 histone/protein deacetylase is required to terminate Hr38 induction. Loss of Hr38 decreases ethanol tolerance and causes precocious but short‐lasting ethanol preference. Similarly, reduced Mef2 activity in all neurons or specifically in the mushroom body α/β neurons decreases ethanol tolerance; Sirt1 promotes ethanol tolerance in these same neurons. Genetically decreasing Hr38 expression levels in Sirt1 null mutants restores ethanol tolerance, demonstrating that both induction and termination of Hr38 expression are important for behavioral plasticity to proceed. These data demonstrate that Hr38 functions as an immediate early transcription factor that promotes ethanol behavioral plasticity.

Are Nonsense Alleles of Drosophila melanogaster Genes under Any Selection?

A gene which carries a bona fide loss-of-function mutation effectively becomes a functionless pseudogene, free from selective constraint. However, there is a number of molecular mechanisms that may lead to at least a partial preservation of the function of genes carrying even drastic alleles. We performed a direct measurement of the strength of negative selection acting on nonsense alleles of protein-coding genes in the Zambian population of Drosophila melanogaster. Within those exons that carry nonsense mutations, negative selection, assayed by the ratio of missense over synonymous nucleotide diversity levels, appears to be absent, consistent with total loss of function. In other exons of nonsense alleles, negative selection was deeply relaxed but likely not completely absent, and the per site number of missense alleles declined significantly with the distance from the premature stop codon. This pattern may be due to alternative splicing which preserves function of some isoforms of nonsense alleles of genes.

A transcriptomics assessment of oxygen-temperature interactions reveals novel candidate genes underlying variation in thermal tolerance and survival

While single stress responses are fairly well researched, multiple, interactive stress responses are not-despite the obvious importance thereof. Here, using D. melanogaster, we investigated the effects of simultaneous exposures to low O₂ (hypoxia) and varying thermal conditions on mortality rates, estimates of thermal tolerance and the transcriptome. We used combinations of 21 (normoxia), 10 or 5kPa O₂ with control (23°C), cold (4°C) or hot (31°C) temperature exposures before assaying chill coma recovery time (CCRT) and heat knock down time (HKDT) as measures of cold and heat tolerance respectively. We found that mortality was significantly affected by temperature, oxygen partial pressure (PO₂) and the interaction between the two. Cold treatments resulted in low mortality (<5%), regardless of PO₂ treatment; while hot treatments resulted in higher mortality (∼20%), especially at 5kPa O₂ which was lethal for most flies (∼80%). Both CCRT and HKDT were significantly affected by temperature, but not PO₂, of the treatments, and the interaction of temperature and PO₂ was non-significant. Hot treatments led to significantly longer CCRT, and shorter HKDT in comparison to cold treatments. Global gene expression profiling provided the first transcriptome level response to the combined stress of PO₂ and temperature, showing that stressful treatments resulted in higher mortality and induced transcripts that were associated with protein kinases, catabolic processes (proteases, hydrolases, peptidases) and membrane function. Several genes and pathways that may be responsible for the protective effects of combined PO₂ and cold treatments were identified. We found that urate oxidase was upregulated in all three cold treatments, regardless of the PO₂. Small heat shock proteins Hsp22 and Hsp23 were upregulated after both 10 and 21kPa O₂-hot treatments. Collectively, the data from PO₂-hot treatments suggests that hypoxia does exacerbate heat stress, through an as yet unidentified mechanism. Hsp70B and an unannotated transcript (CG6733) were significantly differentially expressed after 5kPa O₂-cold and 10kPa O₂-hot treatments relative to their controls. Downregulation of these transcripts was correlated with reduced thermal tolerance (longer CCRT and shorter HKDT), suggesting that these genes may be important candidates for future research.

The roles of plasticity and evolutionary change in shaping gene expression variation in natural populations of extremophile fish

The notorious plasticity of gene expression responses and the complexity of environmental gradients complicate the identification of adaptive differences in gene regulation among populations. We combined transcriptome analyses in nature with common-garden and exposure experiments to establish cause-effect relationships between the presence of a physiochemical stressor and expression differences, as well as to test how evolutionary change and plasticity interact to shape gene expression variation in natural systems. We studied two evolutionarily independent population pairs of an extremophile fish (Poecilia mexicana) living in toxic, hydrogen sulphide (H₂ S)-rich springs and adjacent nontoxic habitats and assessed genomewide expression patterns of wild-caught and common-garden-raised individuals exposed to different concentrations of H₂ S. We found that 7.7% of genes that were differentially expressed between sulphidic and nonsulphidic ecotypes remained differentially expressed in the laboratory, indicating that sources of selection other than H₂ S-or plastic responses to other environmental factors-contribute substantially to gene expression patterns observed in the wild. Concordantly differentially expressed genes in the wild and the laboratory were primarily associated with H₂ S detoxification, sulphur processing and metabolic physiology. While shared, ancestral plasticity played a minor role in shaping gene expression variation observed in nature, we documented evidence for evolved population differences in the constitutive expression as well as the H₂ S inducibility of candidate genes. Mechanisms underlying gene expression variation also varied substantially across the two ecotype pairs. These results provide a springboard for studying evolutionary modifications of gene regulatory mechanisms that underlie expression variation in locally adapted populations.

Allele specific expression analysis identifies regulatory variation associated with stress-related genes in the Mexican highland maize landrace Palomero Toluqueño

BACKGROUND

Gene regulatory variation has been proposed to play an important role in the adaptation of plants to environmental stress. In the central highlands of Mexico, farmer selection has generated a unique group of maize landraces adapted to the challenges of the highland niche. In this study, gene expression in Mexican highland maize and a reference maize breeding line were compared to identify evidence of regulatory variation in stress-related genes. It was hypothesised that local adaptation in Mexican highland maize would be associated with a transcriptional signature observable even under benign conditions.

METHODS

Allele specific expression analysis was performed using the seedling-leaf transcriptome of an F1 individual generated from the cross between the highland adapted Mexican landrace Palomero Toluqueño and the reference line B73, grown under benign conditions. Results were compared with a published dataset describing the transcriptional response of B73 seedlings to cold, heat, salt and UV treatments.

RESULTS

A total of 2,386 genes were identified to show allele specific expression. Of these, 277 showed an expression difference between Palomero Toluqueño and B73 alleles under benign conditions that anticipated the response of B73 cold, heat, salt and/or UV treatments, and, as such, were considered to display a prior stress response. Prior stress response candidates included genes associated with plant hormone signaling and a number of transcription factors. Construction of a gene co-expression network revealed further signaling and stress-related genes to be among the potential targets of the transcription factors candidates.

DISCUSSION

Prior activation of responses may represent the best strategy when stresses are severe but predictable. Expression differences observed here between Palomero Toluqueño and B73 alleles indicate the presence of cis-acting regulatory variation linked to stress-related genes in Palomero Toluqueño. Considered alongside gene annotation and population data, allele specific expression analysis of plants grown under benign conditions provides an attractive strategy to identify functional variation potentially linked to local adaptation.