Heredity

Morales-Arce, Ana Y.; Sabin, Susanna J.; Stone, Anne C.; Jensen, Jeffrey D.

doi: 10.1038/s41437-020-00377-7pmid: 33060846

Recent progress in genomic sequencing from patient samples has allowed for the first detailed insight into the within-host genetic diversity of Mycobacterium tuberculosis (M.TB), revealing remarkably low levels of variation. While this has often been attributed to low mutation rates, other factors have been described, including resistance evolution (i.e., selective sweeps), widespread purifying and background selection, and, more recently, progeny skew. Here we review recent findings pertaining to the processes governing the evolutionary dynamics of M.TB, discuss their implications for improving our understanding of this important human pathogen, and make recommendations for future work. Significantly, this emerging evolutionary framework involving the joint estimation of demographic, selective, and reproductive processes is forming a new paradigm for the study of within-host pathogen evolution that will be widely applicable across organisms.

Wang, Xinxing; Li, Ao; Wang, Wei; Que, Huayong; Zhang, Guofan; Li, Li

doi: 10.1038/s41437-020-0351-7pmid: 32807851

Epigenetic mechanisms such as DNA methylation have the potential to affect organism acclimatization and adaptation to environmental changes by influencing their phenotypic plasticity; however, little is known about the role of methylation in the adaptive phenotypic divergence of marine invertebrates. Therefore, in this study, a typical intertidal species, the Pacific oyster (Crassostrea gigas), was selected to investigate the epigenetic mechanism of phenotypic plasticity in marine invertebrates. Intertidal and subtidal oysters subjected to one-generation common garden experiments and exhibited phenotypic divergence were used. The methylation landscape of both groups of oysters was investigated under temperate and high temperature. The two tidal oysters exhibited divergent methylation patterns, regardless of the temperature, which was mainly original environment-induced. Intertidal samples exhibited significant hypomethylation and more plasticity of methylation in response to heat shock, while subtidal samples showed hypermethylation and less plasticity. Combined with RNA-seq data, a positive relationship between methylation and expression in gene bodies was detected on a genome-wide scale. In addition, approximately 11% and 7% of differentially expressed genes showed significant methylation variation under high temperatures in intertidal and subtidal samples, respectively. Genes related to apoptosis and organism development may be regulated by methylation in response to high temperature in intertidal oysters, whereas oxidation-reduction and ion homeostasis-related genes were involved in subtidal oysters. The results also suggest that DNA methylation mediates phenotypic divergence in oysters adapting to different environments. This study provides new insight into the epigenetic mechanisms underlying phenotypic plasticity in adaptation to rapid climate change in marine organisms.

Fischer, Klaus; Kreyling, Jürgen; Beaulieu, Michaël; Beil, Ilka; Bog, Manuela; Bonte, Dries; Holm, Stefanie; Knoblauch, Sabine; Koch, Dustin; Muffler, Lena; Mouginot, Pierick; Paulinich, Maria; Scheepens, J. F.; Schiemann, Raijana; Schmeddes, Jonas;

Venney, Clare J.; Wellband, Kyle W.; Heath, Daniel D.

doi: 10.1038/s41437-020-0346-4pmid: 32699390

Genetic architecture and phenotypic plasticity are important considerations when studying trait variation within and among populations. Since environmental change can induce shifts in the genetic architecture and plasticity of traits, it is important to consider both genetic and environmental sources of phenotypic variation. While there is overwhelming evidence for environmental effects on phenotype, the underlying mechanisms are less clear. Variation in DNA methylation is a potential mechanism mediating environmental effects on phenotype due to its sensitivity to environmental stimuli, transgenerational inheritance, and influences on transcription. To characterize the effect of environment on methylation, we created two 6 × 6 (North Carolina II) Chinook salmon breeding crosses and reared the offspring in two environments: uniform hatchery tanks and seminatural stream channels. We sampled the fish twice during development, at the alevin (larval) and fry (juvenile) stages. We measured DNA methylation at 13 genes using a PCR-based bisulfite sequencing protocol. The genetic architecture of DNA methylation differed between rearing environments, with greater additive and nonadditive genetic variance in hatchery fish and greater maternal effects in seminatural channel fish, though gene-specific variation was evident. We observed plasticity in methylation across all assayed genes, as well as gene-specific effects at two genes in alevin and six genes in fry, indicating developmental stage-specific effects of rearing environment on methylation. Characterizing genetic and environmental influences on methylation is critical for future studies on DNA methylation as a potential mechanism for acclimation and adaptation.

Ibañez, Verónica Noé; Masuelli, Ricardo Williams; Marfil, Carlos Federico

doi: 10.1038/s41437-020-00355-zpmid: 32801346

DNA methylation can be environmentally modulated and plays a role in phenotypic plasticity. To understand the role of environmentally induced epigenetic variation and its dynamics in natural populations and ecosystems, it is relevant to place studies in a real-world context. Our experimental model is the wild potato Solanum kurtzianum, a close relative of the cultivated potato S. tuberosum. It was evaluated in its natural habitat, an arid Andean region in Argentina characterised by spatial and temporal environmental fluctuations. The dynamics of phenotypic and epigenetic variability (with Methyl Sensitive Amplified Polymorphism markers, MSAP) were assayed in three genotypes across three growing seasons. These genotypes were cultivated permanently and also reciprocally transplanted between experimental gardens (EG) differing in ca. 1000 m of altitude. In two seasons, the genotypes presented differential methylation patterns associated to the EG. In the reciprocal transplants, a rapid epigenomic remodelling occurred according to the growing season. Phenotypic plasticity, both spatial (between EGs within season) and temporal (between seasons), was detected. The epigenetic and phenotypic variability was positively correlated. The lack of an evident mitotic epigenetic memory would be a common response to short-term environmental fluctuations. Thus, the environmentally induced phenotypic and epigenetic variation could contribute to populations persistence through time. These results have implications for understanding the great ecological diversity of wild potatoes.

Griffiths, Sarah M.; Butler, Mark J.; Behringer, Donald C.; Pérez, Thierry; Preziosi, Richard F.

doi: 10.1038/s41437-020-0344-6pmid: 32699391

Understanding population genetic structure can help us to infer dispersal patterns, predict population resilience and design effective management strategies. For sessile species with limited dispersal, this is especially pertinent because genetic diversity and connectivity are key aspects of their resilience to environmental stressors. Here, we describe the population structure of Ircinia campana, a common Caribbean sponge subject to mass mortalities and disease. Microsatellites were used to genotype 440 individuals from 19 sites throughout the Greater Caribbean. We found strong genetic structure across the region, and significant isolation by distance across the Lesser Antilles, highlighting the influence of limited larval dispersal. We also observed spatial genetic structure patterns congruent with oceanography. This includes evidence of connectivity between sponges in the Florida Keys and the southeast coast of the United States (>700 km away) where the oceanographic environment is dominated by the strong Florida Current. Conversely, the population in southern Belize was strongly differentiated from all other sites, consistent with the presence of dispersal-limiting oceanographic features, including the Gulf of Honduras gyre. At smaller spatial scales (<100 km), sites showed heterogeneous patterns of low-level but significant genetic differentiation (chaotic genetic patchiness), indicative of temporal variability in recruitment or local selective pressures. Genetic diversity was similar across sites, but there was evidence of a genetic bottleneck at one site in Florida where past mass mortalities have occurred. These findings underscore the relationship between regional oceanography and weak larval dispersal in explaining population genetic patterns, and could inform conservation management of the species.

Boyrie, Léa; Moreau, Corentin; Frugier, Florian; Jacquet, Christophe; Bonhomme, Maxime

doi: 10.1038/s41437-020-0349-1pmid: 32728044

The quest for signatures of selection using single nucleotide polymorphism (SNP) data has proven efficient to uncover genes involved in conserved and/or adaptive molecular functions, but none of the statistical methods were designed to identify interacting alleles as targets of selective processes. Here, we propose a statistical test aimed at detecting epistatic selection, based on a linkage disequilibrium (LD) measure accounting for population structure and heterogeneous relatedness between individuals. SNP-based (Trv\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$T_{r_v}$$\end{document}) and window-based (TcorPC1v\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$T_{corPC1_v}$$\end{document}) statistics fit a Student distribution, allowing to test the significance of correlation coefficients. As a proof of concept, we use SNP data from the Medicago truncatula symbiotic legume plant and uncover a previously unknown gene coadaptation between the MtSUNN (Super Numeric Nodule) receptor and the MtCLE02 (CLAVATA3-Like) signaling peptide. We also provide experimental evidence supporting a MtSUNN-dependent negative role of MtCLE02 in symbiotic root nodulation. Using human HGDP-CEPH SNP data, our new statistical test uncovers strong LD between SLC24A5 (skin pigmentation) and EDAR (hairs, teeth, sweat glands development) world-wide, which persists after correction for population structure and relatedness in Central South Asian populations. This result suggests that epistatic selection or coselection could have contributed to the phenotypic make-up in some human populations. Applying this approach to genome-wide SNP data will facilitate the identification of coadapted gene networks in model or non-model organisms.

Costa-Neto, Germano; Fritsche-Neto, Roberto; Crossa, José

doi: 10.1038/s41437-020-00353-1pmid: 32855544

Modern whole-genome prediction (WGP) frameworks that focus on multi-environment trials (MET) integrate large-scale genomics, phenomics, and envirotyping data. However, the more complex the statistical model, the longer the computational processing times, which do not always result in accuracy gains. We investigated the use of new kernel methods and modeling structures involving genomics and nongenomic sources of variation in two MET maize data sets. Five WGP models were considered, advancing in complexity from a main-effect additive model (A) to more complex structures, including dominance deviations (D), genotype × environment interaction (AE and DE), and the reaction-norm model using environmental covariables (W) and their interaction with A and D (AW + DW). A combination of those models built with three different kernel methods, Gaussian kernel (GK), Deep kernel (DK), and the benchmark genomic best linear-unbiased predictor (GBLUP/GB), was tested under three prediction scenarios: newly developed hybrids (CV1), sparse MET conditions (CV2), and new environments (CV0). GK and DK outperformed GB in prediction accuracy and reduction of computation time (~up to 20%) under all model–kernel scenarios. GK was more efficient in capturing the variation due to A + AE and D + DE effects and translated it into accuracy gains (~up to 85% compared with GB). DK provided more consistent predictions, even for more complex structures such as W + AW + DW. Our results suggest that DK and GK are more efficient in translating model complexity into accuracy, and more suitable for including dominance and reaction-norm effects in a biologically accurate and faster way.

Chebib, Jobran; Jackson, Benjamin C.; López-Cortegano, Eugenio; Tautz, Diethard; Keightley, Peter D.

doi: 10.1038/s41437-020-00361-1pmid: 32868871

For over a century, inbred mice have been used in many areas of genetics research to gain insight into the genetic variation underlying traits of interest. The generalizability of any genetic research study in inbred mice is dependent upon all individual mice being genetically identical, which in turn is dependent on the breeding designs of companies that supply inbred mice to researchers. Here, we compare whole-genome sequences from individuals of four commonly used inbred strains that were procured from either the colony nucleus or from a production colony (which can be as many as ten generations removed from the nucleus) of a large commercial breeder, in order to investigate the extent and nature of genetic variation within and between individuals. We found that individuals within strains are not isogenic, and there are differences in the levels of genetic variation that are explained by differences in the genetic distance from the colony nucleus. In addition, we employ a novel approach to mutation rate estimation based on the observed genetic variation and the expected site frequency spectrum at equilibrium, given a fully inbred breeding design. We find that it provides a reasonable per nucleotide mutation rate estimate when mice come from the colony nucleus (~7.9 × 10−9 in C3H/HeN), but substantially inflated estimates when mice come from production colonies.

Zimmerman, Shawna J.; Aldridge, Cameron L.; Langin, Kathryn M.; Wann, Gregory T.; Scott Cornman, R.; Oyler-McCance, Sara J.

doi: 10.1038/s41437-020-0352-6pmid: 32807852

The warming climate will expose alpine species adapted to a highly seasonal, harsh environment to novel environmental conditions. A species can shift their distribution, acclimate, or adapt in response to a new climate. Alpine species have little suitable habitat to shift their distribution, and the limits of acclimation will likely be tested by climate change in the long-term. Adaptive genetic variation may provide the raw material for species to adapt to this changing environment. Here, we use a genomic approach to describe adaptive divergence in an alpine-obligate species, the white-tailed ptarmigan (Lagopus leucura), a species distributed from Alaska to New Mexico, across an environmentally variable geographic range. Previous work has identified genetic structure and morphological, behavioral, and physiological differences across the species’ range; however, those studies were unable to determine the degree to which adaptive divergence is correlated with local variation in environmental conditions. We used a genome-wide dataset generated from 95 white-tailed ptarmigan distributed throughout the species’ range and genotype–environment association analyses to identify the genetic signature and environmental drivers of local adaptation. We detected associations between multiple environmental gradients and candidate adaptive loci, suggesting ptarmigan populations may be locally adapted to the plant community composition, elevation, local climate, and to the seasonality of the environment. Overall, our results suggest there may be groups within the species’ range with genetic variation that could be essential for adapting to a changing climate and helpful in guiding conservation action.