Linking genotypes to phenotypes and fitness: how mechanistic biology can inform molecular ecology

The accessibility of new genomic resources, high‐throughput molecular technologies and analytical approaches such as genome scans have made finding genes contributing to fitness variation in natural populations an increasingly feasible task. Once candidate genes are identified, we argue that it is necessary to take a mechanistic approach and work up through the levels of biological organization to fully understand the impacts of genetic variation at these candidate genes. We demonstrate how this approach provides testable hypotheses about the causal links among levels of biological organization, and assists in designing relevant experiments to test the effects of genetic variation on phenotype, whole‐organism performance capabilities and fitness. We review some of the research programs that have incorporated mechanistic approaches when examining naturally occurring genetic and phenotypic variation and use these examples to highlight the value of developing a comprehensive understanding of the relationship between genotype and fitness. We give suggestions to guide future research aimed at uncovering and understanding the genetic basis of adaptation and argue that further integration of mechanistic approaches will help molecular ecologists better understand the evolution of natural populations.     

Epistasis for fitness-related quantitative traits in Arabidopsis thaliana grown in the field and in the greenhouse

The extent to which epistasis contributes to adaptation, population differentiation, and speciation is a long-standing and important problem in evolutionary genetics. Using recombinant inbred (RI) lines of Arabidopsis thaliana grown under natural field conditions, we have examined the genetic architecture of fitness-correlated traits with respect to epistasis; we identified both single-locus additive and two-locus epistatic QTL for natural variation in fruit number, germination, and seed length and width. For fruit number, we found seven significant epistatic interactions, but only two additive QTL. For seed germination, length, and width, there were from two to four additive QTL and from five to eight epistatic interactions. The epistatic interactions were both positive and negative. In each case, the magnitude of the epistatic effects was roughly double that of the effects of the additive QTL, varying from -41% to +29% for fruit number and from -5% to +4% for seed germination, length, and width. A number of the QTL that we describe participate in more than one epistatic interaction, and some loci identified as additive also may participate in an epistatic interaction; the genetic architecture for fitness traits may be a network of additive and epistatic effects. We compared the map positions of the additive and epistatic QTL for germination, seed width, and seed length from plants grown in both the field and the greenhouse. While the total number of significant additive and epistatic QTL was similar under the two growth conditions, the map locations were largely different. We found a small number of significant epistatic QTL x environment effects when we tested directly for them. Our results support the idea that epistatic interactions are an important part of natural genetic variation and reinforce the need for caution in comparing results from greenhouse-grown and field-grown plants.     

New insights into the genetic basis of natural chilling and cold shock tolerance in rice by genome‐wide association analysis

In order to understand cold adaptability and explore additional genetic resources for the cold tolerance improvement of rice, we investigated the genetic variation of 529 rice accessions under natural chilling and cold shock stress conditions at the seedling stage using genome‐wide association studies; a total of 132 loci were identified. Among them, 12 loci were common for both chilling and cold shock tolerance, suggesting that rice has a distinct and overlapping genetic response and adaptation to the two stresses. Haplotype analysis of a known gene OsMYB2, which is involved in cold tolerance, revealed indica–japonica differentiation and latitude tendency for the haplotypes of this gene. By checking the subpopulation and geographical distribution of accessions with tolerance or sensitivity under these two stress conditions, we found that the chilling tolerance group, which mainly consisted of japonica accessions, has a wider latitudinal distribution than the chilling sensitivity group. We conclude that the genetic basis of natural chilling stress tolerance in rice is distinct from that of cold shock stress frequently used for low‐temperature treatment in the laboratory and the cold adaptability of rice is associated with the subpopulation and latitudinal distribution.     

Genetic and epigenetic variations associated with adaptation to heterogeneous habitat conditions in a deciduous shrub

Environmentally induced phenotypic plasticity is thought to play an important role in the adaption of plant populations to heterogeneous habitat conditions, and yet the importance of epigenetic variation as a mechanism of adaptive plasticity in natural plant populations still merits further research. In this study, we investigated populations of Vitex negundo var. heterophylla (Chinese chastetree) from adjacent habitat types at seven sampling sites. Using several functional traits, we detected a significant differentiation between habitat types. With amplified fragment length polymorphisms (AFLP) and methylation‐sensitive AFLP (MSAP), we found relatively high levels of genetic and epigenetic diversity but very low genetic and epigenetic differences between habitats within sites. Bayesian clustering showed a remarkable habitat‐related differentiation and more genetic loci associated with the habitat type than epigenetic, suggesting that the adaptation to the habitat is genetically based. However, we did not find any significant correlation between genetic or epigenetic variation and habitat using simple and partial Mantel tests. Moreover, we found no correlation between genetic and ecologically relevant phenotypic variation and a significant correlation between epigenetic and phenotypic variation. Although we did not find any direct relationship between epigenetic variation and habitat environment, our findings suggest that epigenetic variation may complement genetic variation as a source of functional phenotypic diversity associated with adaptation to the heterogeneous habitat in natural plant populations.     

Genome-Wide Association Analysis of Adaptation Using Environmentally Predicted Traits

Current methods for studying the genetic basis of adaptation evaluate genetic associations with ecologically relevant traits or single environmental variables, under the implicit assumption that natural selection imposes correlations between phenotypes, environments and genotypes. In practice, observed trait and environmental data are manifestations of unknown selective forces and are only indirectly associated with adaptive genetic variation. In theory, improved estimation of these forces could enable more powerful detection of loci under selection. Here we present an approach in which we approximate adaptive variation by modeling phenotypes as a function of the environment and using the predicted trait in multivariate and univariate genome-wide association analysis (GWAS). Based on computer simulations and published flowering time data from the model plant Arabidopsis thaliana, we find that environmentally predicted traits lead to higher recovery of functional loci in multivariate GWAS and are more strongly correlated to allele frequencies at adaptive loci than individual environmental variables. Our results provide an example of the use of environmental data to obtain independent and meaningful information on adaptive genetic variation.     

Quantifying effects of environmental and geographical factors on patterns of genetic differentiation

Elucidating the factors influencing genetic differentiation is an important task in biology, and the relative contribution from natural selection and genetic drift has long been debated. In this study, we used a regression-based approach to simultaneously estimate the quantitative contributions of environmental adaptation and isolation by distance on genetic variation in Boechera stricta, a wild relative of Arabidopsis. Patterns of discrete and continuous genetic differentiation coexist within this species. For the discrete differentiation between two major genetic groups, environment has larger contribution than geography, and we also identified a significant environment-by-geography interaction effect. Elsewhere in the species range, we found a latitudinal cline of genetic variation reflecting only isolation by distance. To further confirm the effect of environmental selection on genetic divergence, we identified the specific environmental variables predicting local genotypes in allopatric and sympatric regions. Water availability was identified as the possible cause of differential local adaptation in both geographical regions, confirming the role of environmental adaptation in driving and maintaining genetic differentiation between the two major genetic groups. In addition, the environment-by-geography interaction is further confirmed by the finding that water availability is represented by different environmental factors in the allopatric and sympatric regions. In conclusion, this study shows that geographical and environmental factors together created stronger and more discrete genetic differentiation than isolation by distance alone, which only produced a gradual, clinal pattern of genetic variation. These findings emphasize the importance of environmental selection in shaping patterns of species-wide genetic variation in the natural environment.     

Adaptive genetic diversity of trees for forest conservation in a future climate: a case study on Norway spruce in Austria

Genetic resources of forest trees are considered as a key factor for the persistence of forest ecosystems because the ability of tree species to survive under changing climate depends strongly on their intraspecific variation in climate response. Therefore, utilizing available genetic variation in climate response and planting alternative provenances suitable for future climatic conditions is considered as an important adaptation measure for forestry. On the other hand, the distribution of adaptive genetic diversity of many tree species is still unknown and the predicted shift of ecological zones and species’ distribution may threaten forest genetic resources that are important for adaptation. Here, we use Norway spruce in Austria as a case study to demonstrate the genetic variation in climate response and to analyse the existing network of genetic conservation units for its effectiveness to safeguard the hotspots of adaptive and neutral genetic diversity of this species. An analysis of the climate response of 480 provenances, clustered into 9 groups of climatically similar provenances, revealed high variation among provenance groups. The most productive and promising provenance clusters for future climates originate from three regions that today depict the warmest and driest areas of the natural spruce distribution in Austria. Gap analysis of the Austrian genetic conservation units in the EUFGIS Portal suggests adequate coverage of the genetic hotspots in southern parts of Austria, but not in eastern and northern Austria. Therefore conservation measures and sustainable utilization of the valuable genetic resources in these regions need to be expanded to cover their high adaptive genetic variation and local adaptation to a warmer climate. The study shows that current conservation efforts need to be evaluated for their effectiveness to protect genetic resources that are important for the survival of trees in a future climate.     

Genotype-by-environment interactions and adaptation to local temperature affect immunity and fecundity in Drosophila melanogaster

Natural populations of most organisms harbor substantial genetic variation for resistance to infection. The continued existence of such variation is unexpected under simple evolutionary models that either posit direct and continuous natural selection on the immune system or an evolved life history "balance" between immunity and other fitness traits in a constant environment. However, both local adaptation to heterogeneous environments and genotype-by-environment interactions can maintain genetic variation in a species. In this study, we test Drosophila melanogaster genotypes sampled from tropical Africa, temperate northeastern North America, and semi-tropical southeastern North America for resistance to bacterial infection and fecundity at three different environmental temperatures. Environmental temperature had absolute effects on all traits, but there were also marked genotype-by-environment interactions that may limit the global efficiency of natural selection on both traits. African flies performed more poorly than North American flies in both immunity and fecundity at the lowest temperature, but not at the higher temperatures, suggesting that the African population is maladapted to low temperature. In contrast, there was no evidence for clinal variation driven by thermal adaptation within North America for either trait. Resistance to infection and reproductive success were generally uncorrelated across genotypes, so this study finds no evidence for a fitness tradeoff between immunity and fecundity under the conditions tested. Both local adaptation to geographically heterogeneous environments and genotype-by-environment interactions may explain the persistence of genetic variation for resistance to infection in natural populations.     

Population of origin and environment interact to determine oomycete infections in spotted salamander populations

Spatial variation in disease risk in wild populations can depend both on environmental and genetic factors. Understanding the various contributions of each factor requires experimental manipulation of both the environment and genetic composition of populations under natural field conditions. We first examined natural patterns of oomycete composition and infection in the eggs of 13 populations of the spotted salamander Ambystoma maculatum. We then performed a fully factorial field transplant of the eggs of six populations to separate the contributions from population of origin and the environment on oomycete resistance in spotted salamanders. Among wild ponds, we found strong variation in oomycete infections in spotted salamander populations and differences in the composition of oomycete communities. In transplant experiments, salamander populations differed in their resistance to oomycete infections via a significant interaction between population of origin and environment. However, not all populations were locally adapted to local conditions. One population was significantly adapted to its home environment, and another one was significantly maladapted. These population effects could originate from differential adaptation of salamander populations to local oomycete communities or environmental conditions that mediate resistance, local adaptation and maladaptation of oomycetes to hosts, or from maternal transmission. Accounting for both environment and population of origin will often be necessary to understand disease dynamics in wild populations.     

Genetic variation and causes of genotype-environment interaction in the body size of blue tit (Parus caeruleus).

In several studies of natural populations of birds, the heritability of body size estimated by parent-offspring regression has been lower when offspring have developed in poor feeding regimens than when they developed in good feeding regimens. This has led to the suggestion that adaptation under poor regimens may be constrained by lack of genetic variation. We examined the influence of environmental conditions on expression of genetic variation in body size of nestling blue tits (Parus caeruleus) by raising full sibs in artificially reduced and enlarged broods, corresponding to good and poor feeding regimens, respectively. Individuals grown in the poor regimen attained smaller body size than their sibs grown in the good regimen. However, there was among-family variation in response to the treatments--i.e., genotype-environment interactions (GEIs). Partitioning the GEI variance into contributions attributable to (1) differences in the among-family genetic variance between the treatments and (2) imperfect correlation of genotypic values across treatments identified the latter as the main cause of the GEI. Parent-offspring regressions were not significantly different when offspring were reared in the good environment (h2 = 0.75) vs. when they were reared in the poor environment (h2 = 0.63). Thus, there was little evidence that genetic variance in body size was lower under the poor conditions than under the good conditions. These results do not support the view that the genetic potential for adaptation to poor feeding conditions is less than that for adaptation to good conditions, but they do suggest that different genotypes may be favored under the different conditions.     

Flowering time QTL in natural populations of Arabidopsis thaliana and implications for their adaptive value

The genetic basis of phenotypic traits is of great interest to evolutionary biologists, but their contribution to adaptation in nature is often unknown. To determine the genetic architecture of flowering time in ecologically relevant conditions, we used a recombinant inbred line population created from two locally adapted populations of Arabidopsis thaliana from Sweden and Italy. Using these RILs, we identified flowering time QTL in growth chambers that mimicked the natural temperature and photoperiod variation across the growing season in each native environment. We also compared the genomic locations of flowering time QTL to those of fitness (total fruit number) QTL from a previous three‐year field study. Ten total flowering time QTL were found, and in all cases, the Italy genotype caused early flowering regardless of the conditions. Two QTL were consistent across chamber environments, and these had the largest effects on flowering time. Five of the fitness QTL colocalized with flowering time QTL found in the Italy conditions, and in each case, the local genotype was favoured. In contrast, just two flowering time QTL found in the Sweden conditions colocalized with fitness QTL and in only one case was the local genotype favoured. This implies that flowering time may be more important for adaptation in Italy than Sweden. Two candidate genes (FLC and VIN3) underlying the major flowering time QTL found in the current study are implicated in local adaptation.     

Population genomic footprints of selection and associations with climate in natural populations of Arabidopsis halleri from the Alps

Natural genetic variation is essential for the adaptation of organisms to their local environment and to changing environmental conditions. Here, we examine genomewide patterns of nucleotide variation in natural populations of the outcrossing herb Arabidopsis halleri and associations with climatic variation among populations in the Alps. Using a pooled population sequencing (Pool‐Seq) approach, we discovered more than two million SNPs in five natural populations and identified highly differentiated genomic regions and SNPs using F_ST‐based analyses. We tested only the most strongly differentiated SNPs for associations with a nonredundant set of environmental factors using partial Mantel tests to identify topo‐climatic factors that may underlie the observed footprints of selection. Possible functions of genes showing signatures of selection were identified by Gene Ontology analysis. We found 175 genes to be highly associated with one or more of the five tested topo‐climatic factors. Of these, 23.4% had unknown functions. Genetic variation in four candidate genes was strongly associated with site water balance and solar radiation, and functional annotations were congruent with these environmental factors. Our results provide a genomewide perspective on the distribution of adaptive genetic variation in natural plant populations from a highly diverse and heterogeneous alpine environment.     

Caging and Uncaging Genetics

It is important for biology to understand if observations made in highly reductionist laboratory settings generalise to harsh and noisy natural environments in which genetic variation is sorted to produce adaptation. But what do we learn by studying, in the laboratory, a genetically diverse population that mirrors the wild? What is the best design for studying genetic variation? When should we consider it at all? The right experimental approach depends on what you want to know.     

Novel loci control variation in reproductive timing in Arabidopsis thaliana in natural environments

Molecular biologists are rapidly characterizing the genetic basis of flowering in model species such as Arabidopsis thaliana. However, it is not clear how the developmental pathways identified in controlled environments contribute to variation in reproductive timing in natural ecological settings. Here we report the first study of quantitative trait loci (QTL) for date of bolting (the transition from vegetative to reproductive growth) in A. thaliana in natural seasonal field environments and compare the results with those obtained under typical growth-chamber conditions. Two QTL specific to long days in the chamber were expressed only in spring-germinating cohorts in the field, and two loci specific to short days in the chamber were expressed only in fall-germinating cohorts, suggesting differential involvement of the photoperiod pathway in different seasonal environments. However, several other photoperiod-specific QTL with large effects in controlled conditions were undetectable in natural environments, indicating that expression of allelic variation at these loci was overridden by environmental factors specific to the field. Moreover, a substantial number of QTL with major effects on bolting date in one or more field environments were undetectable under controlled environment conditions. These novel loci suggest the involvement of additional genes in the transition to flowering under ecologically relevant conditions.     

Genomic signature of natural and anthropogenic stress in wild populations of the waterflea Daphnia magna: validation in space, time and experimental evolution

Natural populations are confronted with multiple selection pressures resulting in a mosaic of environmental stressors at the landscape level. Identifying the genetic underpinning of adaptation to these complex selection environments and assigning causes of natural selection within multidimensional selection regimes in the wild is challenging. The water flea Daphnia is a renowned ecological model system with its well-documented ecology, the possibility to analyse subfossil dormant egg banks and the short generation time allowing an experimental evolution approach. Capitalizing on the strengths of this model system, we here link candidate genome regions to three selection pressures, known to induce micro-evolutionary responses in Daphnia magna: fish predation, parasitism and land use. Using a genome scan approach in space, time and experimental evolution trials, we provide solid evidence of selection at the genome level under well-characterized environmental gradients in the wild and identify candidate genes linked to the three environmental stressors. Our study reveals differential selection at the genome level in Daphnia populations and provides evidence for repeatable patterns of local adaptation in a geographic mosaic of environmental stressors fuelled by standing genetic variation. Our results imply high evolutionary potential of local populations, which is relevant to understand the dynamics of trait changes in natural populations and their impact on community and ecosystem responses through eco-evolutionary feedbacks.     

Exploring genomic variation associated with drought stress in Picea mariana populations

Predicted increases in drought and heat stress will likely induce shifts in species bioclimatic envelopes. Genetic variants adapted to water limitation may prove pivotal for species response under scenarios of increasing drought. In this study, we aimed to explore this hypothesis by investigating genetic variation in 16 populations of black spruce (Picea mariana) in relation to climate variables in Alaska. A total of 520 single nucleotide polymorphisms (SNPs) were genotyped for 158 trees sampled from areas of contrasting climate regimes. We used multivariate and univariate genotype‐by‐environment approaches along with available gene annotations to investigate the relationship between climate and genetic variation among sampled populations. Nine SNPs were identified as having a significant association with climate, of which five were related to drought stress response. Outlier SNPs with respect to the overall environment were significantly overrepresented for several biological functions relevant for coping with variable hydric regimes, including osmotic stress response. This genomic imprint is consistent with local adaptation of black spruce to drought stress. These results suggest that natural selection acting on standing variation prompts local adaptation in forest stands facing water limitation. Improved understanding of possible adaptive responses could inform our projections about future forest dynamics and help prioritize populations that harbor valuable genetic diversity for conservation.     

DNA methylation as a possible mechanism affecting ability of natural populations to adapt to changing climate

Environmentally induced epigenetic variation has been recently recognized as a possible mechanism allowing plants to rapidly adapt to novel conditions. Despite increasing evidence on the topic, little is known on how epigenetic variation affects responses of natural populations to changing climate. We studied the effects of experimental demethylation (DNA methylation is an important mediator of heritable control of gene expression) on performance of a clonal grass, Festuca rubra, coming from localities with contrasting temperature and moisture regimes. We compared performance of demethylated and control plants from different populations under two contrasting climatic scenarios and explored whether the response to demethylation depended on genetic relatedness of the plants. Demethylation significantly affected plant performance. Its effects interacted with population of origin and partly with conditions of cultivation. The effects of demethylation also varied between distinct genotypes with more closely related genotypes showing more similar response to demethylation. For belowground biomass, demethylated plants showed signs of adaptation to drought that were not apparent in plants that were naturally methylated. The results suggest that DNA methylation may modify the response of this species to moisture. DNA methylation may thus affect the ability of clonal plants to adapt to novel climatic conditions. Whether this variation in DNA methylation may also occur under natural conditions, however, remains to be explored. Despite the significant interactions between population of origin and demethylation, our data do not provide clear evidence that DNA methylation enabled adaptation to different environments. In fact, we obtained stronger evidence of local adaptation in demethylated than in naturally‐methylated plants. As changes in DNA methylation may be quite dynamic, it is thus possible that epigenetic variation can mask plant adaptations to conditions of their origin due to pre‐cultivation of the plants under standardized conditions. This possibility should be considered in future experiments exploring plant adaptations.     

Fitness costs associated with low genetic variation are reduced in a harsher environment in amphibian island populations

A basic premise of conservation geneticists is that low levels of genetic variation are associated with fitness costs in terms of reduced survival and fecundity. These fitness costs may frequently vary with environmental factors and should increase under more stressful conditions. However, there is no consensus on how fitness costs associated with low genetic variation change under natural conditions in relation to the stressfulness of the environment. On the Swedish west coast, natterjack toad Bufo calamita populations show a strong population genetic structure and large variation in the amount of within-population genetic variation. We experimentally examined the survival of natterjack larvae from six populations with different genetic variation in three thermal environments corresponding to (a) the mean temperature of natural ponds (stable, laboratory), (b) a high temperature environment occurring in desiccating ponds (stable, laboratory) and (c) an outdoor treatment mimicking the natural, variable thermal conditions (fluctuating, semi-natural). We found that larvae in the outdoor treatment had poorer survival than larvae in the stable environments suggesting that the outdoor treatment was more stressful. Overall, populations with higher genetic variation had higher larval survival. However, a significant interaction between treatments and genetic variation indicated that fitness costs associated with low genetic variation were less severe in the outdoor treatment. Thus, we found no support for the hypothesis that fitness costs associated with low genetic variation increase under more stressful conditions. Our results suggest that natural thermal stress may mask fitness losses associated with low genetic variation in these populations.