Moises Exposito-Alonso

1,135 Genomes Reveal the Global Pattern of Polymorphism in Arabidopsis thaliana

Summary Arabidopsis thaliana serves as a model organism for the study of fundamental physiological, cellular, and molecular processes. It has also greatly advanced our understanding of intraspecific genome variation. We present a detailed map of variation in 1,135 high-quality re-sequenced natural inbred lines representing the native Eurasian and North African range and recently colonized North America. We identify relict populations that continue to inhabit ancestral habitats, primarily in the Iberian Peninsula. They have mixed with a lineage that has spread to northern latitudes from an unknown glacial refugium and is now found in a much broader spectrum of habitats. Insights into the history of the species and the fine-scale distribution of genetic diversity provide the basis for full exploitation of A. thaliana natural variation through integration of genomes and epigenomes with molecular and non-molecular phenotypes.

Effector-Triggered Immune Response in Arabidopsis thaliana Is a Quantitative Trait

We identified loci responsible for natural variation in Arabidopsis thaliana (Arabidopsis) responses to a bacterial pathogen virulence factor, HopAM1. HopAM1 is a type III effector protein secreted by the virulent Pseudomonas syringae strain Pto DC3000. Delivery of HopAM1 from disarmed Pseudomonas strains leads to local cell death, meristem chlorosis, or both, with varying intensities in different Arabidopsis accessions. These phenotypes are not associated with differences in bacterial growth restriction. We treated the two phenotypes as quantitative traits to identify host loci controlling responses to HopAM1. Genome-wide association (GWA) of 64 Arabidopsis accessions identified independent variants highly correlated with response to each phenotype. Quantitative trait locus (QTL) mapping in a recombinant inbred population between Bur-0 and Col-0 accessions revealed genetic linkage to regions distinct from the top GWA hits. Two major QTL associated with HopAM1-induced cell death were also associated with HopAM1-induced chlorosis. HopAM1-induced changes in Arabidopsis gene expression showed that rapid HopAM1-dependent cell death in Bur-0 is correlated with effector-triggered immune responses. Studies of the effect of mutations in known plant immune system genes showed, surprisingly, that both cell death and chlorosis phenotypes are enhanced by loss of EDS1, a regulatory hub in the plant immune-signaling network. Our results reveal complex genetic architecture for response to this particular type III virulence effector, in contrast to the typical monogenic control of cell death and disease resistance triggered by most type III effectors.

The rate and potential relevance of new mutations in a colonizing plant lineage

By following the evolution of populations that are initially genetically homogeneous, much can be learned about core biological principles. For example, it allows for detailed studies of the rate of emergence of de novo mutations and their change in frequency due to drift and selection. Unfortunately, in multicellular organisms with generation times of months or years, it is difficult to set up and carry out such experiments over many generations. An alternative is provided by “natural evolution experiments” that started from colonizations or invasions of new habitats by selfing lineages. With limited or missing gene flow from other lineages, new mutations and their effects can be easily detected. North America has been colonized in historic times by the plant Arabidopsis thaliana, and although multiple intercrossing lineages are found today, many of the individuals belong to a single lineage, HPG1. To determine in this lineage the rate of substitutions—the subset of mutations that survived natural selection and drift–, we have sequenced genomes from plants collected between 1863 and 2006. We identified 73 modern and 27 herbarium specimens that belonged to HPG1. Using the estimated substitution rate, we infer that the last common HPG1 ancestor lived in the early 17th century, when it was most likely introduced by chance from Europe. Mutations in coding regions are depleted in frequency compared to those in other portions of the genome, consistent with purifying selection. Nevertheless, a handful of mutations is found at high frequency in present-day populations. We link these to detectable phenotypic variance in traits of known ecological importance, life history and growth, which could reflect their adaptive value. Our work showcases how, by applying genomics methods to a combination of modern and historic samples from colonizing lineages, we can directly study new mutations and their potential evolutionary relevance.

The rate and effect of de novo mutations in a colonizing lineage of Arabidopsis thaliana

Because colonizations and invasions are often associated with genetic bottlenecks, they offer an opportunity to directly observe de novo mutations and their subsequent fate. North America has recently been colonized by Arabidopsis thaliana, and many of the individuals found today belong to a single lineage, HPG1. To determine substitution rates under natural conditions in this lineage, we have sequenced 100 HPG1 genomes from plants collected between 1863 and 2006. We infer that the last common HPG1 ancestor lived in the early 17th century, most likely the time when HPG1 began to colonize N. America. Demographic reconstructions infer substantial population size fluctuations during the past four centuries. Even though changing demographics can undermine the effect of natural selection, we observed that mutations at coding sites were at lower frequency than mutations at other sites, consistent with the effect of purifying selection. Exceptionally, some mutations rose to high frequency and some had measurable effects in root development, consistent with positive selection acting over mutations with an adaptive value. Our work showcases how by applying genomics methods to a combination of modern and historic samples we can learn about plant colonisations and invasions and observe “evolution in action”.

Adaptive diversification of growth allometry in the plant Arabidopsis thaliana

Significance Are there biological constants unifying phenotypic diversity across scales? Metabolic scaling theory (MST) predicts mathematical regularity and constancy in the allometric scaling of growth rate with body size across species. Here we show that adaptation to climate in Arabidopsis thaliana is associated with local strains that substantially deviate from the values predicted by MST. This deviation can be linked to increased stress tolerance at the expense of seed production, and it occurs through selection on genes that are involved in the abiotic stress response and are geographically correlated with climatic conditions. This highlights the evolutionary role of allometric diversification and helps establish the physiological bases of plant adaptation to contrasting environments. Seed plants vary tremendously in size and morphology; however, variation and covariation in plant traits may be governed, at least in part, by universal biophysical laws and biological constants. Metabolic scaling theory (MST) posits that whole-organismal metabolism and growth rate are under stabilizing selection that minimizes the scaling of hydrodynamic resistance and maximizes the scaling of resource uptake. This constrains variation in physiological traits and in the rate of biomass accumulation, so that they can be expressed as mathematical functions of plant size with near-constant allometric scaling exponents across species. However, the observed variation in scaling exponents calls into question the evolutionary drivers and the universality of allometric equations. We have measured growth scaling and fitness traits of 451 Arabidopsis thaliana accessions with sequenced genomes. Variation among accessions around the scaling exponent predicted by MST was correlated with relative growth rate, seed production, and stress resistance. Genomic analyses indicate that growth allometry is affected by many genes associated with local climate and abiotic stress response. The gene with the strongest effect, PUB4, has molecular signatures of balancing selection, suggesting that intraspecific variation in growth scaling is maintained by opposing selection on the trade-off between seed production and abiotic stress resistance. Our findings suggest that variation in allometry contributes to local adaptation to contrasting environments. Our results help reconcile past debates on the origin of allometric scaling in biology and begin to link adaptive variation in allometric scaling to specific genes.

Spatio-temporal variation in fitness responses to contrasting environments in Arabidopsis thaliana : FITNESS RESPONSES TO NOVEL ENVIRONMENTS

The evolutionary response of organisms to global climate change is expected to be strongly conditioned by preexisting standing genetic variation. In addition, natural selection imposed by global climate change on fitness‐related traits can be heterogeneous over time. We estimated selection of life‐history traits of an entire genetic lineage of the plant Arabidopsis thaliana occurring in north‐western Iberian Peninsula that were transplanted over multiple years into two environmentally contrasting field sites in southern Spain, as southern environments are expected to move progressively northwards with climate change in the Iberian Peninsula. The results indicated that natural selection on flowering time prevailed over that on recruitment. Selection favored early flowering in six of eight experiments and late flowering in the other two. Such heterogeneity of selection for flowering time might be a powerful mechanism for maintaining genetic diversity in the long run. We also found that north‐western A. thaliana accessions from warmer environments exhibited higher fitness and higher phenotypic plasticity for flowering time in southern experimental facilities. Overall, our transplant experiments suggested that north‐western Iberian A. thaliana has the means to cope with increasingly warmer environments in the region as predicted by trends in global climate change models.

A map of climate change-driven natural selection in Arabidopsis thaliana

Through the lens of evolution, climate change is a directional selection acting on populations, forcing them to change and adapt, or face extinction. We studied climate-driven selection at the genetic level by measuring fitness of 517 Arabidopsis thaliana lines grown in rainfall-manipulation experiments in Spain and Germany. Natural selection in the hot-dry Mediterranean location was the strongest, killing 63% of the lines and significantly changing the frequency of about 5 % of all genome-wide variants. Rainfall treatment dictated the direction and magnitude of selection: variants that increased in frequency in a low-rainfall environment, decreased in frequency under high rainfall, and vice versa. A significant proportion of this selection was predictable from the annotation of variants and from the climate in the areas where different variants are found (R2 = 29–52 %). Extrapolating our predictions to other areas of the species range suggested that selection driven by local climate would be strongest in the Mediterranean and Western Siberia — the edge of the species’ environmental limits — and weakest in Central Europe. With rapidly increasing droughts and rising temperatures in Europe, we forecast a wave of directional selection moving North, consequently decimating and causing a genetic turnover of some native A. thaliana populations.

A rainfall-manipulation experiment with 517 Arabidopsis thaliana accessions

The gold standard for studying natural selection and adaptation in the wild is to quantify lifetime fitness of individuals from natural populations that have been grown together in a common garden, or that have been reciprocally transplanted. By combining fitness values with species traits and genome sequences, one can infer selection coefficients at the genetic level. Here we present a rainfall-manipulation experiment with 517 whole-genome sequenced natural accessions of the plant Arabidopsis thaliana spanning the global distribution of the species. The experiments were conducted in two field stations in contrasting climates, in the Mediterranean and in Central Europe, where we built rainout shelters and simulated high and low rainfall. Using custom image analysis we quantified fitness- and phenology-related traits for 23,154 pots, which contained about 14,500 plants growing independently, and over 310,000 plants growing in small populations (max. 30 plants). This large field experiment dataset, which associates fitness and ecologically-relevant traits with genomes, will provide an important resource to test eco-evolutionary genetic theories and to understand the potential evolutionary impacts of future climates on an important plant model species.