AbstractLocal adaptation frequently evolves in patches or environments that are connected via migration. In these cases, genomic regions that are linked to a locally adapted locus experience reduced effective migration rates. Via individual-based simulations of a two-patch system, we show that this reduced effective migration results in the accumulation of conditionally deleterious mutations, but not universally deleterious mutations, adjacent to adaptive loci. When there is redundancy in the genetic basis of local adaptation (i.e., genotypic redundancy), turnover of locally adapted polymorphisms allows conditionally deleterious mutation load to be purged. The amount of mutational load that accumulates adjacent to locally adapted loci is dependent on redundancy, recombination rate, migration rate, population size, strength of selection, and the phenotypic effect size of adaptive alleles. Our results highlight the need to be cautious when interpreting patterns of local adaptation at the level of phenotype or fitness, as the genetic basis of local adaptation can be transient, and evolution may confer a degree of maladaptation to nonlocal environments.

The ability of mutations to facilitate adaptation is central to evolution. To understand how mutations can lead to functional adaptation in a complex molecular machine, we created a defective version of the T4 clamp-loader complex, which is essential for DNA replication. This variant, which is ∼5,000-fold less active than the wild type, was made by replacing the catalytic domains with those from another phage. A directed-evolution experiment revealed that multiple substitutions to a single negatively charged residue in the chimeric clamp loader-Asp 86-restore fitness to within ∼20-fold of wild type. These mutations remove an adventitious electrostatic repulsive interaction between Asp 86 and the sliding clamp. Thus, the fitness decrease of the chimeric clamp loader is caused by a reduction in affinity between the clamp loader and the clamp. Deep mutagenesis shows that the reduced fitness of the chimeric clamp loader is also compensated for by lysine and arginine substitutions of several DNA-proximal residues in the clamp loader or the sliding clamp. Our results demonstrate that there is a latent capacity for increasing the affinity of the clamp loader for DNA and the sliding clamp, such that even single-point mutations can readily compensate for the loss of function due to suboptimal interactions elsewhere.

There is considerable evidence for local adaptation in nature, yet important questions remain regarding its genetic basis. How many loci are involved? What are their effect sizes? What is the relative importance of conditional neutrality versus genetic trade-offs? Here we address these questions in the self-pollinating, annual plant Arabidopsis thaliana. We used 400 recombinant inbred lines (RILs) derived from two locally adapted populations in Italy and Sweden, grew the RILs and parents at the parental locations, and mapped quantitative trait loci (QTL) for mean fitness (fruits/seedling planted). We previously published results from the first 3 years of the study, and here add five additional years, providing a unique opportunity to assess how temporal variation in selection might affect QTL detection and classification. We found 10 adaptive and one maladaptive QTL in Italy, and six adaptive and four maladaptive QTL in Sweden. The discovery of maladaptive QTL at both sites suggests that even locally adapted populations are not always at their genotypic optimum. Mean effect sizes for adaptive QTL, 0.97 and 0.55 fruits in Italy and Sweden, respectively, were large relative to the mean fitness of the RILs (approximately 8 fruits/seedling planted at both sites). Both genetic trade-offs (four cases) and conditional neutrality (seven cases) contribute to local adaptation in this system. The 8-year dataset provided greater power to detect QTL and to estimate their locations compared to our previous 3-year study, identifying one new genetic trade-off and resolving one genetic trade-off into two conditionally adaptive QTL.

Mutations are the cause of several diseases as well as the underlying force of evolution. A thorough understanding of their biophysical consequences is essential. We present a computational framework for evaluating different levels of mutual information (MI) and its dependence on mutation. We used molecular dynamics trajectories of the third PDZ domain and its different mutations. Nonlinear MI between all residue pairs are calculated by tensor Hermite polynomials up to the fifth order and compared with results from multivariate Gaussian distribution of joint probabilities. We show that MI is written as the sum of a Gaussian and a nonlinear component. Results for the PDZ domain show that the Gaussian term gives a sufficiently accurate representation of MI when compared with nonlinear terms up to the fifth order. Changes in MI between residue pairs show the characteristic patterns resulting from specific mutations. Emergence of new peaks in the MI versus residue index plots of mutated PDZ shows how mutation may change allosteric pathways. Triple correlations are characterized by evaluating MI between triplets of residues. We observed that certain triplets are strongly affected by mutation. Susceptibility of residues to perturbation is obtained by MI and discussed in terms of linear response theory.

Ecological differentiation can drive speciation but it is unclear how the genetic architecture of habitat-dependent fitness contributes to lineage divergence. We investigated the genetic architecture of cumulative flowering, a fitness component, in second-generation hybrids between Silene dioica and Silene latifolia transplanted into the natural habitat of each species. We used reduced-representation sequencing and Bayesian sparse linear mixed models (BSLMMs) to analyze the genetic control of cumulative flowering in each habitat. Our results point to a polygenic architecture of cumulative flowering. Allelic effects were mostly beneficial or deleterious in one habitat and neutral in the other. Positive-effect alleles often were derived from the native species, whereas negative-effect alleles, at other loci, tended to originate from the non-native species. We conclude that ecological differentiation is governed and maintained by many loci with small, habitat-dependent effects consistent with conditional neutrality. This pattern may result from differences in selection targets in the two habitats and from environmentally dependent deleterious load. Our results further suggest that selection for native alleles and against non-native alleles acts as a barrier to gene flow between species.

Mutation accumulation (MA) has profound ecological and evolutionary consequences. One example is that accumulation of conditionally neutral mutations leads to fitness trade-offs among heterogenous habitats which cause population divergence. Here we suggest that temperature, which controls the rates of all biochemical and biophysical processes, should play a crucial role for determining mutational effects. Particularly, warmer temperatures may mitigate the effects of some, not all, deleterious mutations and cause stronger environmental dependence in MA effects.
We experimentally tested the above hypothesis by measuring the growth performance of ten Escherichia coli genotypes on six carbon resources across ten temperatures, where the ten genotypes were derived from a single ancestral strain and accumulated spontaneous mutations. We analyzed resource dependence of MA consequences for growth yields. The MA genotypes typically showed reduced growth yields relative to the ancestral type; and the magnitude of reduction was smaller at intermediate temperatures. Stronger resource dependence in MA consequences for growth performance was observed at higher temperatures. Specifically, the MA genotypes were more likely to show impaired growth performance on all the six carbon resources when grown at lower temperatures; but suffered growth performance loss only on some, not all the six, carbon substrates at higher temperatures.
Higher temperatures increase the chance that MA causes conditionally neutral fitness effects while MA is more likely to cause fitness loss regardless of available resources at lower temperatures. This finding has implications for understanding how geographic patterns in population divergence may emerge, and how conservation practices, particularly protection of diverse microhabitats, may mitigate the impacts of global warming.

Ionomics measures elemental concentrations in biological organisms and provides a snapshot of physiology under different conditions. In this study, we evaluate genetic variation of the ionome in outbred, perennial switchgrass in three environments across the species\' native range, and explore patterns of genotype-by-environment interactions. We grew 725 clonally replicated genotypes of a large full sib family from a four-way linkage mapping population, created from deeply diverged upland and lowland switchgrass ecotypes, at three common gardens. Concentrations of 18 mineral elements were determined in whole post-anthesis tillers using ion coupled plasma mass spectrometry (ICP-MS). These measurements were used to identify quantitative trait loci (QTL) with and without QTL-by-environment interactions (QTLxE) using a multi-environment QTL mapping approach. We found that element concentrations varied significantly both within and between switchgrass ecotypes, and GxE was present at both the trait and QTL level. Concentrations of 14 of the 18 elements were under some genetic control, and 77 QTL were detected for these elements. 74% of QTL colocalized multiple elements, half of QTL exhibited significant QTLxE, and roughly equal numbers of QTL had significant differences in magnitude and sign of their effects across environments. The switchgrass ionome is under moderate genetic control and by loci with highly variable effects across environments.

Toxicant resistance is a complex trait, affected both by genetics and the environment. Like most complex traits, it can exhibit sexual dimorphism, yet sex is often overlooked as a factor in studies of toxicant resistance. Paraquat, one such toxicant, is a commonly used herbicide and is known to produce mitochondrial oxidative stress, decrease dopaminergic neurons and dopamine (DA) levels, and decrease motor ability. While the main effects of paraquat are well-characterized, less is known about the naturally occurring variation in paraquat susceptibility. The purpose of this study was to map the genes contributing to low-dose paraquat susceptibility in Drosophila melanogaster, and to determine if susceptibility differs between the sexes. One hundred of the Drosophila Genetic Reference Panel (DGRP) lines were scored for susceptibility via climbing ability and used in a genome-wide association study (GWAS). Variation in seventeen genes in females and thirty-five genes in males associated with paraquat susceptibility. Only two candidate genes overlapped between the sexes despite a significant positive correlation between male and female susceptibilities. Many associated polymorphisms had significant interactions with sex, with most having conditionally neutral effects. Conditional neutrality between the sexes probably stems from sex-biased expression which may result from partial resolution of sexual conflict. Candidate genes were verified with RNAi knockdowns, gene expression analyses, and DA quantification. Several of these genes are novel associations with paraquat susceptibility. This research highlights the importance of assessing both sexes when studying toxicant susceptibility.

The central abundance hypothesis predicts that local adaptation is a function of the distance to the centre of a species\' geographic range. To test this hypothesis, we gathered genomic diversity data from 49 populations, 646 individuals and 33,464 SNPs of two wild relatives of maize, the teosintes Zea mays ssp. parviglumis and Zea. mays. ssp. mexicana. We examined the association between the distance to their climatic and geographic centroids and the enrichment of SNPs bearing signals of adaptation. We identified candidate adaptive SNPs in each population by combining neutrality tests and cline analyses. By applying linear regression models, we found that the number of candidate SNPs is positively associated with niche suitability, while genetic diversity is reduced at the limits of the geographic distribution. Our results suggest that overall, populations located at the limit of the species\' niches are adapting locally. We argue that local adaptation to this limit could initiate ecological speciation processes and facilitate adaptation to global change.

A longstanding idea in evolutionary physiology is that an enzyme cannot jointly optimize performance at both high and low temperatures due to a trade-off between stability and activity. Although a stability-activity trade-off has been observed for well-characterized examples, such a trade-off is not imposed by any physical chemical constraint. To better understand the pervasiveness of this trade-off, I investigated the stability-activity relationship for comparative biochemical studies of purified orthologous enzymes identified by a literature search. The nature of this relationship varied greatly among studies. Notably, studies of enzymes with low mean synonymous nucleotide sequence divergence were less likely to exhibit the predicted negative correlation between stability and activity. Similarly, a survey of directed evolution investigations of the stability-activity relationship indicated that these traits are often uncoupled among nearly identical yet phenotypically divergent enzymes. This suggests that the presumptive trade-off often reported for investigations of enzymes with high mean sequence divergence may in some cases instead be a consequence of the degeneration over time of enzyme function in unselected environments, rather than a direct effect of thermal adaptation. The results caution against the general assertion of a stability-activity trade-off during enzyme adaptation.