Intergenerational adaptations to stress are evolutionarily conserved, stress-specific, and have deleterious trade-offs
Abstract
Despite reports of parental exposure to stress promoting physiological adaptations in progeny in diverse organisms, there remains considerable debate over the significance and evolutionary conservation of such multigenerational effects. Here, we investigate four independent models of intergenerational adaptations to stress in C. elegans - bacterial infection, eukaryotic infection, osmotic stress and nutrient stress - across multiple species. We found that all four intergenerational physiological adaptations are conserved in at least one other species, that they are stress-specific, and that they have deleterious trade-offs in mismatched environments. By profiling the effects of parental bacterial infection and osmotic stress exposure on progeny gene expression across species we established a core set of 587 genes that exhibited a greater than 2-fold intergenerational change in expression in response to stress in C. elegans and at least one other species, as well as a set of 37 highly conserved genes that exhibited a greater than 2-fold intergenerational change in expression in all four species tested. Furthermore, we provide evidence suggesting that presumed adaptive and deleterious intergenerational effects are molecularly related at the gene expression level. Lastly, we found that none of the effects we detected of these stresses on C. elegans F1 progeny gene expression persisted transgenerationally three generations after stress exposure. We conclude that intergenerational responses to stress play a substantial and evolutionarily conserved role in regulating animal physiology and that the vast majority of the effects of parental stress on progeny gene expression are reversible and not maintained transgenerationally.
Data availability
RNA-seq data that support the findings of this study have been deposited at NCBI GEO and are available under the accession code GSE173987.
-
The time-resolved transcriptome of C. elegansNCBI Sequence Read Archive - Supplemental Table 1.
Article and author information
Author details
Funding
Centre Trophoblast Research (Next Generation fellowship)
- Nicholas O Burton
National Institutes of Health
- Kinsey Fisher
- L Ryan Baugh
National Institutes of Health (GM117408)
- L Ryan Baugh
Natural Sciences and Engineering Research Council of Canada (Grant #522691522691)
- Alexandra Willis
- Aaron W Reinke
Alfred P Sloan Research Fellowship (FG2019-12040)
- Aaron W Reinke
Cancer Research UK (C13474/A18583)
- Eric A Miska
Cancer Research UK (C6946/A14492)
- Eric A Miska
Wellcome Trust (104640/Z/14/Z)
- Eric A Miska
Wellcome Trust (092096/Z/10/Z)
- Eric A Miska
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Copyright
© 2021, Burton et al.
This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,560
- views
-
- 363
- downloads
-
- 32
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Evolutionary Biology
Spatial patterns in genetic diversity are shaped by individuals dispersing from their parents and larger-scale population movements. It has long been appreciated that these patterns of movement shape the underlying genealogies along the genome leading to geographic patterns of isolation-by-distance in contemporary population genetic data. However, extracting the enormous amount of information contained in genealogies along recombining sequences has, until recently, not been computationally feasible. Here, we capitalize on important recent advances in genome-wide gene-genealogy reconstruction and develop methods to use thousands of trees to estimate per-generation dispersal rates and to locate the genetic ancestors of a sample back through time. We take a likelihood approach in continuous space using a simple approximate model (branching Brownian motion) as our prior distribution of spatial genealogies. After testing our method with simulations we apply it to Arabidopsis thaliana. We estimate a dispersal rate of roughly 60 km2/generation, slightly higher across latitude than across longitude, potentially reflecting a northward post-glacial expansion. Locating ancestors allows us to visualize major geographic movements, alternative geographic histories, and admixture. Our method highlights the huge amount of information about past dispersal events and population movements contained in genome-wide genealogies.
-
- Evolutionary Biology
Understanding the genomic basis of natural variation in plant pest resistance is an important goal in plant science, but it usually requires large and labor-intensive phenotyping experiments. Here, we explored the possibility that non-target reads from plant DNA sequencing can serve as phenotyping proxies for addressing such questions. We used data from a whole-genome and -epigenome sequencing study of 207 natural lines of field pennycress (Thlaspi arvense) that were grown in a common environment and spontaneously colonized by aphids, mildew, and other microbes. We found that the numbers of non-target reads assigned to the pest species differed between populations, had significant SNP-based heritability, and were associated with climate of origin and baseline glucosinolate contents. Specifically, pennycress lines from cold and thermally fluctuating habitats, presumably less favorable to aphids, showed higher aphid DNA load, i.e., decreased aphid resistance. Genome-wide association analyses identified genetic variants at known defense genes but also novel genomic regions associated with variation in aphid and mildew DNA load. Moreover, we found several differentially methylated regions associated with pathogen loads, in particular differential methylation at transposons and hypomethylation in the promoter of a gene involved in stomatal closure, likely induced by pathogens. Our study provides first insights into the defense mechanisms of Thlaspi arvense, a rising crop and model species, and demonstrates that non-target whole-genome sequencing reads, usually discarded, can be leveraged to estimate intensities of plant biotic interactions. With rapidly increasing numbers of large sequencing datasets worldwide, this approach should have broad application in fundamental and applied research.