Drought adaptation in Arabidopsis thaliana by extensive genetic loss-of-function
Abstract
Interdisciplinary syntheses are needed to scale up discovery of the environmental drivers and molecular basis of adaptation in nature. Here we integrated novel approaches using whole genome sequences, satellite remote sensing, and transgenic experiments to study natural loss-of-function alleles associated with drought histories in wild Arabidopsis thaliana. The genes we identified exhibit population genetic signatures of parallel molecular evolution, selection for loss-of-function, and shared associations with flowering time phenotypes in directions consistent with longstanding adaptive hypotheses seven times more often than expected by chance. We then confirmed predicted phenotypes experimentally in transgenic knockout lines. These findings reveal the importance of drought timing to explain the evolution of alternative drought tolerance strategies and further challenge popular assumptions about the adaptive value of genetic loss-of-function in nature. These results also motivate improved species-wide sequencing efforts to better identify loss-of-function variants and inspire new opportunities for engineering climate resilience in crops.
https://doi.org/10.7554/eLife.41038.001eLife digest
Water shortages caused by droughts lead to crop losses that affect billions of people around the world each year. By discovering how wild plants adapt to drought, it may be possible to identify traits and genes that help to improve the growth of crop plants when water is scarce. It has been suggested that plants have adapted to droughts by flowering at times of the year when droughts are less likely to occur. For example, if droughts are more likely to happen in spring, the plants may delay flowering until the summer.
Arabidopsis thaliana is a small plant that is found across Eurasia, Africa and North America, including in areas that are prone to drought at different times of the year. Individual plants of the same species may carry different versions of the same gene (known as alleles). Some of these alleles may not work properly and are referred to as loss-of-function alleles. Monroe et al. investigated whether A. thaliana plants carry any loss-of-function alleles that are associated with droughts happening in the spring or summer, and whether they are linked to when those plants will flower.
Monroe et al. analyzed satellite images collected over the last 30 years to measure when droughts have occurred. Next, they searched genome sequences of Arabidopsis thaliana for alleles that might help the plants to adapt to droughts in the spring or summer. Combining the two approaches revealed that loss-of-function alleles associated with spring droughts were strongly predicted to be associated with the plants flowering later in the year. Similarly, loss-of-function alleles associated with summer droughts were predicted to be associated with the plants flowering earlier in the year.
These findings support the idea that plants can adapt to drought by changing when they produce flowers, and suggest that loss-of-function alleles play a major role in this process. New techniques for editing genes mean it is easier than ever to generate new loss-of-function alleles in specific genes. Therefore, the results presented by Monroe et al. may help researchers to develop new varieties of crop plants that are better adapted to droughts.
https://doi.org/10.7554/eLife.41038.002Introduction
Discovering the environmental drivers and functional genetics of adaptation in nature is a key goal of evolutionary biology and valuable to advance applied genetics in agriculture. Understanding the genetics of drought adaptation in plants is particularly important as crop losses resulting from droughts affect billions of people each year, posing the greatest threat to global food stability. Because droughts also impose strong selection on natural plant populations, investigating drought adaptation in wild species is both useful for addressing fundamental questions of evolutionary biology, such as determining whether adaptation proceeds by few or many alleles, and informative for efforts to reverse engineer drought tolerance in crops (Mickelbart et al., 2015). Such an evolutionary research program is motivated by the need to understand adaptive drought tolerance strategies for different types of drought conditions, which can vary in severity and timing (Tardieu, 2012). Furthermore, previous limitations of single gene approaches have reinforced the necessity of developing methods to identify beneficial alleles at genomic scales and functional molecular resolutions (Dean and Thornton, 2007; Passioura, 2010).
Drought stress can occur throughout the year and drought timing is forecast to change over the next century (Trenberth et al., 2014). While dramatic evolutionary responses to drought events have been documented, (e.g. Franks et al., 2007), little is known about the relationship between drought timing and adaptation. However, the observation both in nature and agriculture that plants are particularly susceptible to drought while flowering (Nam et al., 2001; Dietrich and Smith, 2016) has contributed to the longstanding hypothesis that adaptive flowering time should reflect patterns in the seasonal timing of drought events (Passioura, 1996). Detailed studies of life history also reveal that locally adapted Arabidopsis thaliana (Arabidopsis hereafter) populations begin flowering in their home environments just prior to and after periods of increased historical drought frequency (Mojica et al., 2016).
Flowering time in Arabidopsis is correlated with other drought tolerance traits such as water use efficiency and can serve as a proxy for alternative drought tolerance strategies, with early flowering genotypes being associated with low water use efficiency (drought escape strategy) and late flowering genotypes with high water use efficiency (dehydration avoidance strategy) (McKay et al., 2003; Lovell et al., 2013; Kenney et al., 2014). Thus, the historical timing of drought experienced by locally adapted populations may explain the evolution of these strategies and the distribution of alleles responsible for natural flowering time variation. This hypothesis motivated our investigation to identify alleles associated with drought timing and test the prediction that they contribute to adaptive flowering time evolution.
Identifying functionally relevant genetic variation contributing to adaptation is needed to understand fundamental evolutionary processes. In contrast to early theoretical predictions and popular assumptions, loss-of-function (LoF) alleles, those that eliminate or ‘knockout’ a gene’s molecular function, are overrepresented among alleles reported as responsible for crop improvement and often produce adaptive phenotypes in wild species (Hoekstra et al., 2006; Rausher, 2008; Olsen and Wendel, 2013; Alonso-Blanco and Méndez-Vigo, 2014; Weigel and Nordborg, 2015b; Torkamaneh et al., 2018). Indeed, a number of individual genes exhibiting evidence of locally adaptive loss-of-function have been documented in Arabidopsis (Grant et al., 1998; Johanson et al., 2000; Kliebenstein, 2001; Kroymann et al., 2003; Mouchel et al., 2004; Aukerman, 1997; Hauser et al., 2001; Mauricio et al., 2003; Alonso-Blanco et al., 2005; Werner et al., 2005; Barboza et al., 2013; Xiang et al., 2014).
Discovering adaptive LoF alleles is particularly valuable for inspiring targeted molecular breeding because functionally similar mutations can be mined from the breeding pool or generated directly by non-transgenic native gene editing. Unfortunately, traditional genome-wide association scans based on the one-locus two-allele model perform poorly at detecting adaptive LoF alleles, which because of the large number of mutations that can create them, are likely to arise through parallel molecular evolution (Pennings and Hermisson, 2006; Barboza et al., 2013; Kerdaffrec et al., 2016). Species-wide whole genome sequences however, present the opportunity to advance beyond previous mapping and scanning methods that relied on linked polymorphisms by instead characterizing and contrasting functionally defined alleles.
Here, we combined long-term satellite-detected drought histories, whole genome sequence scans based on allele function, and transgenic knockout experiments in Arabidopsis to test historical predictions about how drought timing shapes the evolution of flowering time and outline a broadly scalable approach for discovering loss-of-function gene variants contributing to plant climate adaptation.
Results and discussion
To study global seasonal drought timing, satellite-detected measurements offer a valuable historical record. One such measurement, the Vegetative Health Index (VHI) has been used for decades to monitor drought, including in many places across the natural range of Arabidopsis (Kogan, 1997). Though primarily used as a tool to predict crop productivity, by quantifying drought induced vegetative stress this index also provides a resource for evolutionary ecologists to study seasonal patterns in drought-related episodes of natural selection. We analyzed 34 years of VHI data to characterize drought regimens at the home environments of Arabidopsis ecotypes (Figure 1, Supplementary file 1). We found that drought frequency during the spring (ß = 50.016, p < 2×10−16) and summer (ß = −28.035, p = 4.4×10−7) significantly predict flowering time among Arabidopsis ecotypes (Supplementary file 2A). We then generated a drought-timing index that quantifies the relative frequency of drought between spring and summer over the typical reproductive growing season and observed substantial differences in drought timing experienced by ecotypes (Figure 1—figure supplement 1). This environmental variation presented a useful cline to address classical hypotheses about the evolution of flowering time in relation to drought timing and identify LoF alleles potentially contributing to this evolution.
To identify candidate LoF alleles underlying drought adaptation and flowering time evolution, we analyzed whole genome sequences in Arabidopsis. We first surveyed the genomes of 1135 ecotypes (1001 Genomes Consortium, 2016) for LoF alleles in protein coding genes predicted to encode truncated amino acid sequences (Supplementary file 3A). To overcome the likely parallel evolutionary origins of LoF alleles that would have challenged previous methods, we classified alleles based functional allele state rather than individual polymorphisms for association testing. After filtering to reduce the likelihood of false positives (see materials and methods), we thus tested 2088 genes for LoF allele associations with drought timing (Figure 2A) and flowering time (Figure 2B). These analyses identified 247 genes in which LoF alleles are significantly associated with drought timing and/or flowering time after accounting for population structure and multiple testing (Supplementary file 3B). In contrast, when we performed these analyses on a permuted LoF genotype matrix, we found no genes that were significantly associated with drought timing or flowering time (Figure 1—figure supplement 1).
It should be noted that the 2088 genes tested for associations to flowering time and drought timing are not a complete representation of LoF alleles in Arabidopsis. In some cases, previously studied LoF alleles did not pass filtering steps (Supplementary file 3D,E). This was primarily because the frequency or quality of LoF allele calls in these genes fell below our filtering requirements (see materials and methods). In other cases, the Col-0 reference genome already has a documented LOF allele. Finally, we expect LoF alleles to be undetectable if they are the product of large insertions or deletions which cannot be properly identified with currently available resequencing data. Thus, while the methods used here are designed to minimize false positives (alleles classified as LoF, but which are actually functional), the likely occurrence of false negatives (undetected LoF alleles) in available data motivates the need for more sophisticated species wide genome sequencing efforts including a greater diversity of de-novo quality genomes for comprehensive detection of functionally relevant genetic variation across the species.
Associations to drought timing predicted associations of LoF alleles to flowering time directly. Together, summer drought and earlier flowering associated genes (Figure 2C), and spring drought and later flowering associated genes (Figure 2D) overlapped seven times more often than expected by chance (χ2=492, p < 2 × 10−16) and no shared associations were observed in the opposite direction. The strengths of the associations between LoF alleles and drought timing (P values) was also strongly correlated with the strengths of the associations to flowering time (r2 = 0.48. Figure 2—figure supplement 1E, Figure 2C,D). This result is comparable to overlapping peaks in a ‘Manhattan plot’ generated from a traditional genome wide association scan (e.g. Bosse et al., 2017). In contrast, these associations were weakly correlated when genotypes were permuted (r2 = 0.01 Figure 2—figure supplement 1F), indicating that the result is not simply explained as an artifact of allele frequencies or by the relationship between drought timing and flowering time (i.e. Supplementary file 1A). Thus, satellite-detected drought histories and a functional genome-wide scanning approach prove useful for predicting the direction and molecular targets of phenotypic evolution. Similar investigations with ecologically meaningful environmental variation could be valuable for discovering candidates underlying other important traits that are especially difficult to measure.
These results further support the classical hypothesis that the relationship between phenology and drought timing is the most important feature of plant drought tolerance (Passioura, 1996), indicating the evolution of ‘drought escape’ through earlier flowering in summer drought environments, and ‘dehydration avoidance’ by later flowering genotypes in spring drought environments. Because most Arabidopsis populations appear to exhibit a winter annual life habit, germinating in the fall and overwintering as a rosette (Ratcliffe, 1961; Thompson, 1994; Burghardt et al., 2015), late flowering genotypes in spring drought environments are expected to still encounter drought conditions. However, delayed flowering may ensure that droughts co-occur with vegetative growth rather than during the drought sensitive reproductive phase. This pattern is also consistent with hypotheses explaining the more water conservative water use and stomatal traits observed in late flowering genotypes (McKay et al., 2003; Lovell et al., 2013; Kenney et al., 2014; Kooyers, 2015) and those from spring drought environments (Dittberner et al., 2018). Future experimental work will be valuable to identify other plant physiological traits affected by the LoF alleles associated with drought timing.
These results provide new insight into the ecology and genetics of Arabidopsis life history evolution, but the complex ecological reality of these processes is undoubtedly beyond the scope of this study. We found that drought timing remains a significant predictor of allele associations to flowering time when controlling for allele associations with latitude and minimum temperature (slope estimate in multiple linear regression, p < 2×10−16, Supplementary file 2B). However, other unknown climatic variables or environmental interactions and non-linearities likely contribute to the flowering time adaptation as well. Flowering time is only one component of phenology and other adaptive life history transitions such a germination timing (Donohue, 2002) may also be influenced by drought timing and could change how drought timing affects the evolution of flowering time, a hypothesis that warrants further investigation. Furthermore, measuring flowering time in other environments, such alternate light regimes, may yield a different set of candidate genes using similar approaches.
Signatures of selection in the genes identified differ from the genome average and neutral expectations. As expected for genes harboring LoF alleles, these show parallel evolution of LoF and accelerated amino acid sequence evolution among Arabidopsis ecotypes (Figure 2—figure supplement 2A,B, Supplementary file 2C). We also found evidence of positive selection for LoF alleles in genes associated with drought timing and/or flowering time. While these genes have similar global frequencies of LoF alleles compared to genes not showing associations with drought timing and/or flowering time (Figure 2—figure supplement 2C), they tend to have significantly fewer unique LoF alleles (Figure 2—figure supplement 2D) and greater frequencies of each independent LoF allele (Figure 2E). This pattern is consistent with theoretical predictions and results from simulations of adaptation by parallel molecular evolution involving recurrent mutation combined with more rapid local fixation of alleles experiencing positive selection (Pennings and Hermisson, 2006). In cases where adaptation proceeds through the fixation of a single adaptive allele, traditional genome scanning approaches may be sufficient to detect causal loci. However, when genetic variation consists of multiple independent alleles, as is often the case for the genes examined here (Figure 2—figure supplement 2D), classifying alleles functionally before testing for associations is likely necessary.
The extent of LoF responsible for adaptive phenotypic evolution is much greater than once assumed (Smith, 1970; Albalat and Cañestro, 2016). LoF alleles identified were overwhelmingly associated with spring drought or later flowering rather than summer drought or earlier flowering (χ2 = 132, p < 2 × 10−16, Figure 2). Because the reference genome and gene models are from an early flowering Arabidopsis line, Col-0, this is consistent with the hypothesis that LoF alleles are particularly important in the evolution of phenotypic divergence (Rausher, 2008). This result also highlights the need to develop functional genomics resources informed by multiple de-novo quality reference genomes. We found that flowering time is strongly predicted by the accumulation of LoF alleles across the 214 candidate genes associated to spring drought and/or later flowering time (Figure 3A–E), estimating a 1 day increase for every three additional LoF alleles across these candidate genes (Figure 3F). This relationship is best represented as a simple linear regression; the addition of a non-linear quadratic predictor variable did not significantly improve the fit of the model (F = 0.7005, p = 0.4028). Importantly, we did not find a broader overabundance of LoF alleles in later flowering ecotypes or those from spring drought environments that would explain this relationship (e.g. Figure 2—figure supplement 3). Rather, these findings support a model of climate-associated evolution in complex traits that includes a substantial contribution from widespread genetic LoF and give promise to targeted LoF for directed phenotypic engineering.
Experimental knockout lines confirmed the later flowering times predicted from natural allele associations. To test phenotypic effects, we screened a panel of confirmed T-DNA insertion mutants representing a sample of candidate LoF alleles associated with spring drought and/or later flowering. As predicted by variation among Arabidopsis ecotypes (Figure 2D), the vast majority of knockout lines in these candidate genes (57 of 59, χ2 = 51, p = 8.045e-13) flowered later on average than the wild type genotype (Figure 3G, Supplementary file SF). LoF alleles identified through these analyses and experiments include those previously linked to flowering time (Cui et al., 2007) and drought responses (Aghdasi et al., 2012; Qin et al., 2017). Implementing a functional genome-wide association scan, we find that allele associations with ecologically meaningful environmental variation (drought timing) accurately predict associations with adaptive phenotypes directly (flowering time).
Together with validation in transgenic lines, these findings outline a scalable model for gaining deeper insights into the functional genomics of climate adaptation in nature. Combining large scale knockout experiments with functional genome wide association scans may be a valuable approach for future research to quantify the power to predict LoF allele effects. These results also further challenge historical assumptions about molecular adaptation that have implications for influencing evolutionary theory and public attitudes toward emerging molecular breeding approaches.
Groundbreaking yield increases during the green revolution of the 1960 s were largely attributable to semi-dwarf phenotypes caused by LoF alleles in both rice and barley (Spielmeyer et al., 2002; Jia et al., 2009). Later it was found that natural LoF alleles of the same gene in wild Arabidopsis produce similar phenotypes (Barboza et al., 2013), suggesting the potential to mine ecological species for information directly useful for crop improvement. Visions of a second green revolution powered and informed by such natural variation call for discoveries in evolutionary functional genomics at scales that have now become possible. The genes identified here could inspire future molecular breeding of climate resilient crops and this work more broadly highlights the value of integrating diverse disciplines to scale up the discovery of the climatic drivers of adaptation and functionally significant genetic variation at molecular resolutions.
Materials and methods
Satellite-Detected drought histories of Arabidopsis
Request a detailed protocolTo study patterns in historical drought, the remotely sensed Vegetative Health Index (VHI) was used, a satellite-detected drought measurement tool whose advantage is that it includes information about vegetative impacts of drought (Passioura, 1996; AghaKouchak et al., 2015). This index is based on multiple data sources from NOAA satellites, combining deviations from historic climatic (Temperature Condition Index derived from AVHRR-based observations in thermal bands) and vegetative conditions (Vegetative Condition Index derived from NDVI) to detect periods of ecological drought conditions and distinguish between other sources of vegetative stress such as cold (Kogan, 1997; Kogan et al., 2005; Rojas et al., 2011). VHI was collected weekly since 1981 at 16 km2 resolution on a scale from 0 to 100, where values below 40 reflect drought conditions (Kogan, 1997) (Figure 1A). The frequencies of observing drought conditions during photoperiodic spring (quarter surrounding spring equinox), summer (quarter surrounding summer solstice), fall (quarter surrounding fall equinox), and winter (quarter surrounding winter solstice) were calculated globally from 1981 to 2015 (Figure 1B) in R (R Core Development Team, 2017) using the raster package (Hijmans, 2016).
After removing ecotypes with missing location data or locations falling within pixels classified as water, seasonal drought frequencies and drought timing were calculated at the location of origin for 1,097 Arabidopsis ecotypes that were included as part of the 1001 Genomes Project (1001 Genomes Consortium, 2016) (Figure 1C, Supplementary file 1). Up to date global map files of seasonal drought frequency and the drought-timing index used here are available on Dryad and greymonroe.github.io/data alongside a brief tutorial showing how to extract data for points of interest in R. We tested whether seasonal drought frequencies significantly predicted with flowering time (flowering time described in subsequent section regarding LoF associations) by multiple linear regression (Supplementary file 2A)
To characterize the seasonal timing of droughts during an important period of Arabidopsis’ life history, a univariate drought-timing index was generated that quantifies whether the historical frequency of drought increases or decreases over the course of the typical Arabidopsis reproductive growing season (Ratcliffe, 1961; Thompson, 1994; Burghardt et al., 2015). Specifically, this index is equal to the natural log transformed ratio between spring and summer drought frequency. More negative values reflect environments where drought frequency increases from spring to summer and are referred to here as ‘summer drought environments,’ (e.g. Figure 1B left). Conversely, more positive values reflect environments where drought frequency decreases from spring to summer and are referred to here as ‘spring drought environments,’ (e.g. Figure 1B right).
Loss-of-Function (LoF) Alleles in Arabidopsis genomes
Request a detailed protocolTo identify functionally definitive gene variants (Hoekstra and Coyne, 2007; Weigel and Nordborg, 2015a; Byers et al., 2017), LoF alleles (Albalat and Cañestro, 2016) were identified from whole genome sequence data of 1,135 Arabidopsis accessions (Olson, 1999; Cutter and Jovelin, 2015; 1001 Genomes Consortium, 2016) using R scripts. First, genes were filtered to those containing at least 5% frequency of predicted frameshift or premature stop mutations and less than 5% missing allele calls from results generated by the 1,001 Genomes Consortium (1001 Genomes Consortium, 2016) using ‘SnpEff’ (Cingolani et al., 2012). To reduce instances where exon skipping might ameliorate LoF mutations (Gan et al., 2011), genes were filtered to those with a single predicted gene model (Lamesch et al., 2012). Additionally, to preclude false LoF calls for cases where compensatory mutations restore gene function or in which an insignificant portion of the final protein product is affected by putative LoF mutations (MacArthur et al., 2012), coding regions were translated into predicted amino acid sequences from which lengths from start to stop codon were calculated in R. LoF alleles were defined as those producing protein products with at least 10% lost because of late start codons and/or prematurely truncated translation. Allelic heterogeneity expected to mask these genes from traditional GWAS (Remington, 2015; Monroe et al., 2016; Flood and Hancock, 2017) was corrected for by classifying all alleles as either functional (0) or non-functional (1). A final frequency filter was re-applied (5% global LoF allele frequency), resulting in 2088 genes for downstream association analyses (Supplementary file 3B). Finally, to compare the results of this pipeline to genes known to harbor natural LoF alleles (Mouchel et al., 2004; Shindo et al., 2008; Gujas et al., 2012; Kliebenstein, 2001; Kroymann et al., 2003; Grant et al., 1998; Tian et al., 2003; Mauricio et al., 2003; Werner et al., 2005; Aukerman, 1997; Flowers et al., 2009; Xiang et al., 2014; Xiang et al., 2016; Amiguet-Vercher et al., 2015; Johanson et al., 2000; Le Corre et al., 2002; McKay et al., 2003; Stinchcombe et al., 2004; Shindo et al., 2005; Flowers et al., 2009; Méndez-Vigo et al., 2011; Lovell et al., 2013; Hauser et al., 2001; Bloomer et al., 2012; Alonso-Blanco et al., 2005; Zhen and Ungerer, 2008; Kang et al., 2013; Monroe et al., 2016; Zhu et al., 2015; Barboza et al., 2013), we manually performed this functional allele calling approach on a set of 16 genes (Supplementary file D,E)
LoF associations to drought timing and flowering time
Request a detailed protocolTo identify candidate LoF alleles responsible for climate adaptation and phenotypic evolution, the relationships between functional allele state and drought timing and between functional allele state and flowering time were evaluated for each of the 2088 genes that passed preceding filtering steps. Specifically, the association between functional allele state among Arabidopsis ecotypes and historical drought timing at their locations of origin was tested by logistic regression in a generalized linear model in R (R Core Development Team, 2017). This association study differs from traditional GWAS in several respects. First, because the alleles studied here are functionally defined, they are expected to be more likely to have a phenotypic impact than random SNPs. Second, the scope of our analyses were restricted to a subset of the genome - 2088 genes with high confidence LoF allele calls that passed previous filtering steps, rather than tens of thousands to millions of SNPs. Finally, in contrast to traditional GWAS, which is designed to identify associated chromosomal regions rather than functionally definitive genetic variations, our approach is motivated by the ability to identify alleles at molecular resolutions whose functional relevance can be tested empirically. Thus, the balance of opportunity costs related to trade-offs between false positive and false negative associations that generally challenge GWAS are shifted to reduce false negatives rather than minimizing false positives. For these reasons, we implemented analyses based on (Price et al., 2006) to balance false positives and false negatives. Population structure was accounted for by performing a principal component analysis on the kinship matrix among all ecotypes and including in each model the first three resulting principal components, which explain >75% of variance in relatedness between ecotypes (Price et al., 2006). The P-values (Pdrought timing) of the slope estimates (βdrought timing) for drought timing in these models were adjusted to account for multiple tests by a Bonferroni correction to identify those significantly associated (Supplementary file 3C).
Summer drought genes were identified as those in which LoF alleles are found in ecotypes that experience a significantly (βdrought timing <0 and Pdrought timing <0.05) more negative drought-timing index (summer drought environments where drought frequency increases over the course of the reproductive growing season, Figure 1B left and Figure 2A top). Conversely, spring drought genes were identified as those in which LoF alleles are found in ecotypes that experience a significantly (βdrought timing >0 and Pdrought timing <0.05) more positive drought-timing index (spring drought environments where drought frequency decreases over the course of the reproductive growing season, Figure 1B right and Figure 2A bottom).
The above analytical approach was repeated to test whether functional allele state is associated with the reported common garden flowering times of Arabidopsis ecotypes (Alonso-Blanco and Méndez-Vigo, 2014) (Supplementary file 1). See Alonso-Blanco et al. (Alonso-Blanco and Méndez-Vigo, 2014) for details, but in brief, flowering time was measured in growth chambers at 10°C (considerably less missing data than experiment at 16°C) under 16 hour days. Earlier flowering genes were identified as those in which LoF alleles are found in ecotypes that flower significantly (βflowering time <0 and Pflowering time <0.05) earlier than ecotypes with a functional allele (Figure 2B top). Later flowering genes were identified as those in which LoF alleles are found in ecotypes that flower significantly (βflowering time >0 and Pflowering time <0.05) later than ecotypes with a functional allele (Figure 2B bottom). The preceding analyses revealed considerable overlap between genes associated with both drought timing and flowering time. To assess whether this result was an artifact of the binary LoF allele calls, we randomly permuted the genotype matrix and repeated the analyses described above, testing for significant associations between allele states and drought timing and/or flowering time. Quantile-quantile plots of P values were visualized using qqPlot in the GWASTools package in R (Gogarten et al., 2012) (Figure 2—figure supplement 1A–D)
Overlap between drought timing and flowering time associated genes
Request a detailed protocolTo address the longstanding hypothesis that flowering time reflects adaptation to drought timing (Fox, 1990; Passioura, 1996; Kooyers, 2015), and to test the corresponding prediction that alleles associated with drought timing are also associated with flowering time, the groups of genes identified with significant associations to drought timing or flowering time were compared (Figure 2C and D). Deviation from the null hypothesis of independent associations to drought timing and flowering time was evaluated by a chi-squared test (Expected number of co-associated genes = 12, Observed = 83, χ2 = 492, p = 2×10−16).
The magnitude of P-values have historically served as the basis of selecting candidate loci for further examination toward their contribution to environmental adaptation or phenotypic evolution in quantitative trait locus mapping and genome wide association scans [e.g. (Bosse et al., 2017). To test whether associations to environment (drought timing) can be used to identify loci associated with phenotypes (flowering time) directly, the correlation between log transformed P-values describing allele associations with drought timing (Pdrought timing) and with flowering time (Pflowering time) was calculated (Figure 2—figure supplement 1E, r2 = 0.48,) and visualized separately for genes associated to summer drought/earlier flowering (Figure 2C) and to spring drought/later flowering (Figure 2D). To control for the possibility that allele frequencies or the relationship between drought timing and flowering time explained these observations, we also tested whether allele associations were correlated when generated from association analyses using a matrix of randomly permuted genotypes with the same allele frequencies (Figure 2—figure supplement 1F, r2 = 0.01).
Finally, to control for the possibility that correlated LoF allele associations were explained by confounding environmental variables we tested whether the LoF allele associations to drought timing remained predictive while accounting for LoF allele associations with latitude and minimum temperature of the coldest month (Hijmans et al., 2005) using a multiple linear regression in R (Supplementary file 3B). To do so, we repeated the association analyses described in the previous section but instead tested for LoF allele associations with latitude and minimum temperatures. We then included these P values (Supplementary file 2B) in a multiple linear regression where the strength of the association to flowering time was predicted by the associations to drought timing, latitude, and minimum temperature simultaneously.
Signatures of selection
Request a detailed protocolTo assess whether histories of selection for genes identified differ from the genome wide expectation, measures of amino acid sequence evolution were evaluated for 122 genes in which loss-of-function is associated with drought timing or flowering time and for which there are orthologs identified between A. lyrata and A. thaliana (Goodstein et al., 2012). For each gene, sequences were aligned using MAFFT (Katoh and Standley, 2013), codons with gaps removed, and the number of non-synonymous and synonymous polymorphisms among A. thaliana accessions (PN and PS) as well as synonymous and non-synonymous divergence (DN and DS) from A. lyrata were measured using mkTest.rb (https://github.com/kern-lab/). The ratios PN/PS and DN/DS were then calculated to measure the proportion of variants predicted to affect amino acid sequences that are segregating among ecotypes and diverged from A. lyrata, respectively. These calculations were also performed for genes not associated to drought timing or flowering time (n = 912) and the remaining genes across the A. thaliana genome (n = 20373) with orthologs between A. lyrata and A. thaliana. To test whether genes identified show evidence of accelerated protein sequence evolution, comparisons were made to genes associated with drought timing or flowering time for both PN/PS (Figure 2—figure supplement 2A) and DN/DS ((Figure 2—figure supplement 2A,B) by two-sided students t-tests (α = 0.05) in R (R Core Development Team, 2017).
Because theory predicts adaptation by loss-of-function to proceed through multiple independent alleles, but to exhibit a fewer number of different alleles than in neutral loci at similar LoF allele frequencies (Pennings and Hermisson, 2006; Ralph and Coop, 2010; Ralph and Coop, 2015), the number of unique LoF alleles was estimated by protein length in the genes that passed preceding filtering steps. To address the hypothesis that genes in which LoF alleles are associated to drought history or flowering time are likely to reflect positive selection compared to genes in which LoF are random with respect to drought history or flowering time, the total number of unique LoF alleles between these groups was compared using a two-sided students t-test (log10 transformed, p = 5.8×10−7, (Figure 2—figure supplement 2D). To control for the possibility that this result in an artifact of reduced frequency of LoF alleles in genes identified, the global frequency of LoF was also compared between these groups (log10 transformed, two-sided students t-test, p = 0.11, (Figure 2—figure supplement 2C). Finally, to further test the prediction that LoF alleles in genes identified have increased in frequency because of more positive selection, the frequency per specific LoF allele was compared between groups (log10 transformed, two-sided students t-test, p = 3.4×10−7, Figure 2E).
Candidate genes contributing to later flowering time by widespread LoF
Request a detailed protocolThe significance of the tendency for LoF associations to spring drought/later flowering time (Figure 2D) was tested by chi-squared tests (spring drought vs. summer drought, p < 2×10−16; later vs. earlier flowering, p < 2×10−16, spring drought/later flowering vs. summer drought/earlier flowering, p < 2×10−16). The chromosomal locations of candidate genes (those associated to spring drought/later flowering time) were mapped onto the Arabidopsis genome (Lamesch et al., 2012) (Figure 3A). To address the hypothesis that widespread LoF contributes to later flowering time phenotypes, the total number of LoF in candidate genes for each ecotype was calculated and the correlation between this value and flowering time evaluated (Figure 3F, r2 = 0.39, p < 2×10−16). We also tested whether a model which included a non-linear predictor (squared value of the total number of LoF in candidate genes) was a better fit than the simple linear model by an analysis of variance (F = 0.7005, p = 0.4028).
Experimental testing of predicted phenotypes in gene knockout lines
Request a detailed protocolThe preceding analyses provided compelling evidence of LoF in candidate genes as important in the evolution of later flowering time phenotypes. To test the prediction that non-functionalization of these genes causes increased flowering time, phenotypes were measured in transgenic lines in a subsample of candidate genes showing a significant association between loss-of-function and spring drought environments and/or later flowering time. Motivated by the general need to develop a high throughput approach of studying naturally adaptive LoF, knockout lines from the Arabidopsis Biological Resource Center were chosen from a collection created by the SALK Institute in which a T-DNA insertion in an exon of candidate genes has already been identified and confirmed to be homozygous (O'Malley and Ecker, 2010; Rutter et al., 2017). These T-DNA knockout lines were generated by the SALK institute (Supplementary file 3F) and exist in a common genetic background (Columbia) (Alonso et al., 2003). Seeds were planted in 2’ pots containing wet potting soil and stratified for 5 days at 4°C. Seedlings were thinned to a single plant per pot one week after stratification. Plants were grown (59 T-DNA knockout lines, 10 reps of each line and 30 reps Columbia) in a stratified (by shelf), randomized design in growth chambers (Conviron ATC60, Controlled Environments, Winnipeg, MB) under 16 hr of light at 20°C. Flowering time was measured as days after planting to the emergence of the first open flower, based on the definition of flowering time used by the 1,001 Genomes Consortium (1001 Genomes Consortium, 2016). We calculated the least squares mean (lsmean from ‘lsmeans’ package in R) flowering time for each line from a mixed model where shelf and tray were included as random effects (Supplementary file 3F). We tested the prediction that knockout lines would flower later (have higher lsmean flowering time estimates) than the wild type Columbia genotype by a chi-squared test (p = 8.1×10−13).
Data availability
All data generated or analyzed during this study are included in the manuscript and supporting files.
-
1001 Genomes Data CenterGMI-MPI Arabidopsis thaliana genomes.
References
-
Remote sensing of drought: Progress, challenges and opportunitiesReviews of Geophysics 53:452–480.https://doi.org/10.1002/2014RG000456
-
Cloning and expression analysis of Arabidopsis TRR14 gene under salt and drought stressJournal of Cell and Molecular Research 4:1–10.https://doi.org/10.22067/jcmr.v4i1.12269
-
Genetic architecture of naturally occurring quantitative traits in plants: an updated synthesisCurrent Opinion in Plant Biology 18:37–43.https://doi.org/10.1016/j.pbi.2014.01.002
-
Mechanistic approaches to the study of evolution: the functional synthesisNature Reviews Genetics 8:675–688.https://doi.org/10.1038/nrg2160
-
The genomic basis of adaptation in plantsCurrent Opinion in Plant Biology 36:88–94.https://doi.org/10.1016/j.pbi.2017.02.003
-
Population genomics of the Arabidopsis thaliana flowering time gene networkMolecular Biology and Evolution 26:2475–2486.https://doi.org/10.1093/molbev/msp161
-
Drought and the evolution of flowering time in desert annualsAmerican Journal of Botany 77:1508–1518.https://doi.org/10.1002/j.1537-2197.1990.tb12563.x
-
Phytozome: a comparative platform for green plant genomicsNucleic Acids Research 40:D1178–D1186.https://doi.org/10.1093/nar/gkr944
-
Very high resolution interpolated climate surfaces for global land areasInternational Journal of Climatology 25:1965–1978.https://doi.org/10.1002/joc.1276
-
raster: Geographic Data Analysis and Modelingraster: Geographic Data Analysis and Modeling, https://rdrr.io/cran/raster/.
-
GA-20 oxidase as a candidate for the semidwarf gene sdw1/denso in barleyFunctional & Integrative Genomics 9:255–262.https://doi.org/10.1007/s10142-009-0120-4
-
MAFFT multiple sequence alignment software version 7: improvements in performance and usabilityMolecular Biology and Evolution 30:772–780.https://doi.org/10.1093/molbev/mst010
-
Global drought watch from spaceBulletin of the American Meteorological Society 78:621–636.https://doi.org/10.1175/1520-0477(1997)078<0621:GDWFS>2.0.CO;2
-
Modelling corn production in China using AVHRR‐based vegetation health indicesInternational Journal of Remote Sensing 26:2325–2336.https://doi.org/10.1080/01431160500034235
-
The Arabidopsis Information Resource (TAIR): improved gene annotation and new toolsNucleic Acids Research 40:D1202–D1210.https://doi.org/10.1093/nar/gkr1090
-
Pleiotropy of FRIGIDA enhances the potential for multivariate adaptationProceedings of the Royal Society B: Biological Sciences 280:20131043.https://doi.org/10.1098/rspb.2013.1043
-
Natural selection for polymorphism in the disease resistance gene Rps2 of Arabidopsis thalianaGenetics 163:735–746.
-
Genetic mechanisms of abiotic stress tolerance that translate to crop yield stabilityNature Reviews Genetics 16:237–251.https://doi.org/10.1038/nrg3901
-
Effect of timing of drought stress on growth and grain yield of extra-short-duration pigeonpea linesThe Journal of Agricultural Science 136:179–189.https://doi.org/10.1017/S0021859601008607
-
A bountiful harvest: genomic insights into crop domestication phenotypesAnnual Review of Plant Biology 64:47–70.https://doi.org/10.1146/annurev-arplant-050312-120048
-
When less is more: gene loss as an engine of evolutionary changeThe American Journal of Human Genetics 64:18–23.https://doi.org/10.1086/302219
-
Scaling up: the essence of effective agricultural researchFunctional Plant Biology 37:585–591.https://doi.org/10.1071/FP10106
-
SoftwareR: A language and environment for statistical computingR Foundation for Statistical Computing, Vienna, Austria.
-
Adaptation to habitat in a group of annual plantsThe Journal of Ecology 49:187–203.https://doi.org/10.2307/2257433
-
Evolutionary transitions in floral ColorInternational Journal of Plant Sciences 169:7–21.https://doi.org/10.1086/523358
-
Assessing drought probability for agricultural areas in Africa with coarse resolution remote sensing imageryRemote Sensing of Environment 115:343–352.https://doi.org/10.1016/j.rse.2010.09.006
-
Parallel Loss-of-Function at the RPM1 bacterial resistance locus in arabidopsis thalianaFrontiers in Plant Science 3:287.https://doi.org/10.3389/fpls.2012.00287
-
Fitness effects of mutation: testing genetic redundancy in Arabidopsis thalianaJournal of Evolutionary Biology 30:1124–1135.https://doi.org/10.1111/jeb.13081
-
Any trait or trait-related allele can confer drought tolerance: just design the right drought scenarioJournal of Experimental Botany 63:25–31.https://doi.org/10.1093/jxb/err269
-
Global warming and changes in droughtNature Climate Change 4:17–22.https://doi.org/10.1038/nclimate2067
-
Population genomics for understanding adaptation in wild plant speciesAnnual Review of Genetics 49:315–338.https://doi.org/10.1146/annurev-genet-120213-092110
-
Population genomics for understanding adaptation in wild plant speciesAnnual Review of Genetics 49:315–338.https://doi.org/10.1146/annurev-genet-120213-092110
-
Relaxed selection on the CBF/DREB1 regulatory genes and reduced freezing tolerance in the southern range of Arabidopsis thalianaMolecular Biology and Evolution 25:2547–2555.https://doi.org/10.1093/molbev/msn196
Article and author information
Author details
Funding
National Science Foundation (1701918)
- John Grey Monroe
U.S. Department of Agriculture (2014- 38420-21801)
- John Grey Monroe
National Science Foundation (IOS-1402393)
- John T Lovell
National Science Foundation (1022196)
- John K McKay
National Science Foundation (1556262)
- John K McKay
Cargill (Research support)
- John K McKay
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
E Buckler, D Des Marais, A Henry, J Lasky, T Mitchell-Olds, J Ross-Ibarra, and D Sloan provided valuable feedback and insightful discussion that improved this work. This study was financially supported by NSF Awards DEB 1022196 and 1556262 to JKM, NSF Award 1701918 and USDA-NIFA Award 2014-38420-21801 to JGM, as well as generous funding from Cargill, Inc. The work conducted by the US Department of Energy Joint Genome Institute, a DOE Office of Science User Facility, is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231. Data used are included in the main text, supplementary materials, and public repositories.
Copyright
© 2018, Monroe et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 8,571
- views
-
- 968
- downloads
-
- 64
- 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
- Genetics and Genomics
Evolutionary arms races can arise at the contact surfaces between host and viral proteins, producing dynamic spaces in which genetic variants are continually pursued. However, the sampling of genetic variation must be balanced with the need to maintain protein function. A striking case is given by protein kinase R (PKR), a member of the mammalian innate immune system. PKR detects viral replication within the host cell and halts protein synthesis to prevent viral replication by phosphorylating eIF2α, a component of the translation initiation machinery. PKR is targeted by many viral antagonists, including poxvirus pseudosubstrate antagonists that mimic the natural substrate, eIF2α, and inhibit PKR activity. Remarkably, PKR has several rapidly evolving residues at this interface, suggesting it is engaging in an evolutionary arms race, despite the surface’s critical role in phosphorylating eIF2α. To systematically explore the evolutionary opportunities available at this dynamic interface, we generated and characterized a library of 426 SNP-accessible nonsynonymous variants of human PKR for their ability to escape inhibition by the model pseudosubstrate inhibitor K3, encoded by the vaccinia virus gene K3L. We identified key sites in the PKR kinase domain that harbor K3-resistant variants, as well as critical sites where variation leads to loss of function. We find K3-resistant variants are readily available throughout the interface and are enriched at sites under positive selection. Moreover, variants beneficial against K3 were also beneficial against an enhanced variant of K3, indicating resilience to viral adaptation. Overall, we find that the eIF2α-binding surface of PKR is highly malleable, potentiating its evolutionary ability to combat viral inhibition.
-
- Ecology
- Evolutionary Biology
Seasonal polyphenism enables organisms to adapt to environmental challenges by increasing phenotypic diversity. Cacopsylla chinensis exhibits remarkable seasonal polyphenism, specifically in the form of summer-form and winter-form, which have distinct morphological phenotypes. Previous research has shown that low temperature and the temperature receptor CcTRPM regulate the transition from summer-form to winter-form in C. chinensis by impacting cuticle content and thickness. However, the underling neuroendocrine regulatory mechanism remains largely unknown. Bursicon, also known as the tanning hormone, is responsible for the hardening and darkening of the insect cuticle. In this study, we report for the first time on the novel function of Bursicon and its receptor in the transition from summer-form to winter-form in C. chinensis. Firstly, we identified CcBurs-α and CcBurs-β as two typical subunits of Bursicon in C. chinensis, which were regulated by low temperature (10 °C) and CcTRPM. Subsequently, CcBurs-α and CcBurs-β formed a heterodimer that mediated the transition from summer-form to winter-form by influencing the cuticle chitin contents and cuticle thickness. Furthermore, we demonstrated that CcBurs-R acts as the Bursicon receptor and plays a critical role in the up-stream signaling of the chitin biosynthesis pathway, regulating the transition from summer-form to winter-form. Finally, we discovered that miR-6012 directly targets CcBurs-R, contributing to the regulation of Bursicon signaling in the seasonal polyphenism of C. chinensis. In summary, these findings reveal the novel function of the neuroendocrine regulatory mechanism underlying seasonal polyphenism and provide critical insights into the insect Bursicon and its receptor.