Archaic lineages broaden our view on the history of Arabidopsis thaliana

Museum genomics reveals temporal genetic stasis and global genetic diversity in Arabidopsis thaliana

Global patterns of population genetic variation through time offer a window into evolutionary processes that maintain diversity. Over time, lineages may expand or contract their distribution, causing turnover in population genetic composition. At individual loci, migration, drift, environmental change (among other processes) may affect allele frequencies. Museum specimens of widely distributed species offer a unique window into the genetics of understudied populations and changes over time. Here, we sequenced genomes of 130 herbarium specimens and 91 new field collections of Arabidopsis thaliana and combined these with published genomes. We sought a broader view of genomic diversity across the species, and to test if population genomic composition is changing through time. We documented extensive and previously uncharacterized diversity in a range of populations in Africa, populations that are under threat from anthropogenic climate change. Through time, we did not find dramatic changes in genomic composition of populations. Instead, we found a pattern of genetic change every 100 years of the same magnitude seen when comparing Eurasian populations that are 185 km apart, potentially due to a combination of drift and changing selection. We found only mixed signals of polygenic adaptation at phenology and physiology QTL. We did find that genes conserved across eudicots show altered levels of directional allele frequency change, potentially due to variable purifying and background selection. Our study highlights how museum specimens can reveal new dimensions of population diversity and show how wild populations are evolving in recent history.

The bHLH transcription factor gene EGL3 accounts for the natural diversity in Arabidopsis fruit trichome pattern and morphology

The number and distribution of trichomes, i.e. the trichome pattern, in different plant organs show a conspicuous inter- and intraspecific diversity across Angiosperms that is presumably involved in adaptation to numerous environmental factors. The genetic and molecular mechanisms accounting for the evolution of trichome patterns have just begun to be elucidated. In this study, we aimed to identify and characterize MALAMBRUNO 1 (MAU1), a locus affecting trichome number in the fruits of Arabidopsis (Arabidopsis thaliana) natural populations. To this end, we developed introgression lines (ILs) from the hairy fruit accession Doñana (Don-0) in the genetic background of the Ler strain with glabrous fruits. Genetic analyses of ILs showed that MAU1 affects fruit trichome patterns through synergistic epistasis with the MYB genes TRICHOMELESS1 (TCL1), GLABRA1 (GL1), and TRIPTYCHON (TRY). In addition, fine mapping and characterization of transgenic lines demonstrated that MAU1 is the bHLH transcription factor gene EGL3, for which Don-0 carries a gain-of-function semidominant allele. Gene expression analyses did not detect differences between EGL3 alleles, thus supporting that a structural missense mutation is the causal nucleotide polymorphism of Don-0. Further phylogenetic analyses of EGL3 showed that most Arabidopsis populations with hairy fruits belong to 3 haplogroups, suggesting that additional EGL3 natural alleles account for fruit trichome development. Finally, the characterization of EGL3 pleiotropy indicates that Don-0 hyperfunction also increases stem trichome branching. We conclude that EGL3 interactions in the core gene regulatory network of trichome development explain the Arabidopsis natural diversity for fruit trichome pattern and morphology.

The trichome pattern diversity of Cardamine shares genetic mechanisms with Arabidopsis but differs in environmental drivers

Natural variation in trichome pattern (amount and distribution) is prominent among populations of many angiosperms. However, the degree of parallelism in the genetic mechanisms underlying this diversity and its environmental drivers in different species remain unclear. To address these questions, we analyzed the genomic and environmental bases of leaf trichome pattern diversity in Cardamine hirsuta, a relative of Arabidopsis (Arabidopsis thaliana). We characterized 123 wild accessions for their genomic diversity, leaf trichome patterns at different temperatures, and environmental adjustments. Nucleotide diversities and biogeographical distribution models identified two major genetic lineages with distinct demographic and adaptive histories. Additionally, C. hirsuta showed substantial variation in trichome pattern and plasticity to temperature. Trichome amount in C. hirsuta correlated positively with spring precipitation but negatively with temperature, which is opposite to climatic patterns in A. thaliana. Contrastingly, genetic analysis of C. hirsuta glabrous accessions indicated that, like for A. thaliana, glabrousness is caused by null mutations in ChGLABRA1 (ChGL1). Phenotypic genome-wide association studies (GWAS) further identified a ChGL1 haplogroup associated with low trichome density and ChGL1 expression. Therefore, a ChGL1 series of null and partial loss-of-function alleles accounts for the parallel evolution of leaf trichome pattern in C. hirsuta and A. thaliana. Finally, GWAS also detected other candidate genes (e.g. ChETC3, ChCLE17) that might affect trichome pattern. Accordingly, the evolution of this trait in C. hirsuta and A. thaliana shows partially conserved genetic mechanisms but is likely involved in adaptation to different environments.

The genomics and physiology of abiotic stressors associated with global elevational gradients in Arabidopsis thaliana

Phenotypic and genomic diversity in Arabidopsis thaliana may be associated with adaptation along its wide elevational range, but it is unclear whether elevational clines are consistent among different mountain ranges. We took a multi-regional view of selection associated with elevation. In a diverse panel of ecotypes, we measured plant traits under alpine stressors (low CO2 partial pressure, high light, and night freezing) and conducted genome-wide association studies. We found evidence of contrasting locally adaptive regional clines. Western Mediterranean ecotypes showed low water use efficiency (WUE)/early flowering at low elevations to high WUE/late flowering at high elevations. Central Asian ecotypes showed the opposite pattern. We mapped different candidate genes for each region, and some quantitative trait loci (QTL) showed elevational and climatic clines likely maintained by selection. Consistent with regional heterogeneity, trait and QTL clines were evident at regional scales (c. 2000 km) but disappeared globally. Antioxidants and pigmentation rarely showed elevational clines. High elevation east African ecotypes might have higher antioxidant activity under night freezing. Physiological and genomic elevational clines in different regions can be unique, underlining the complexity of local adaptation in widely distributed species, while hindering global trait-environment or genome-environment associations. To tackle the mechanisms of range-wide local adaptation, regional approaches are thus warranted.

Estimating dispersal rates and locating genetic ancestors with genome-wide genealogies

Spatial patterns in genetic diversity are shaped by individuals dispersing from their parents and larger-scale population movements. It has long been appreciated that these patterns of movement shape the underlying genealogies along the genome leading to geographic patterns of isolation-by-distance in contemporary population genetic data. However, extracting the enormous amount of information contained in genealogies along recombining sequences has, until recently, not been computationally feasible. Here, we capitalize on important recent advances in genome-wide gene-genealogy reconstruction and develop methods to use thousands of trees to estimate per-generation dispersal rates and to locate the genetic ancestors of a sample back through time. We take a likelihood approach in continuous space using a simple approximate model (branching Brownian motion) as our prior distribution of spatial genealogies. After testing our method with simulations we apply it to Arabidopsis thaliana. We estimate a dispersal rate of roughly 60 km2/generation, slightly higher across latitude than across longitude, potentially reflecting a northward post-glacial expansion. Locating ancestors allows us to visualize major geographic movements, alternative geographic histories, and admixture. Our method highlights the huge amount of information about past dispersal events and population movements contained in genome-wide genealogies.

Ecological trade-offs drive phenotypic and genetic differentiation of Arabidopsis thaliana in Europe

Plant diversity is shaped by trade-offs between traits related to competitive ability, propagule dispersal, and stress resistance. However, we still lack a clear understanding of how these trade-offs influence species distribution and population dynamics. In Arabidopsis thaliana, recent genetic analyses revealed a group of cosmopolitan genotypes that successfully recolonized Europe from its center after the last glaciation, excluding older (relict) lineages from the distribution except for their north and south margins. Here, we tested the hypothesis that cosmopolitans expanded due to higher colonization ability, while relicts persisted at the margins due to higher tolerance to competition and/or stress. We compared the phenotypic and genetic differentiation between 71 European genotypes originating from the center, and the south and north margins. We showed that a trade-off between plant fecundity and seed mass shapes the differentiation of A. thaliana in Europe, suggesting that the success of the cosmopolitan groups could be explained by their high dispersal ability. However, at both north and south margins, we found evidence of selection for alleles conferring low dispersal but highly competitive and stress-resistance abilities. This study sheds light on the role of ecological trade-offs as evolutionary drivers of the distribution and dynamics of plant populations.

Climate biogeography of Arabidopsis thaliana: linking distribution models and individual variation

AIM

Patterns of individual variation are key to testing hypotheses about the mechanisms underlying biogeographic patterns. If species distributions are determined by environmental constraints, then populations near range margins may have reduced performance and be adapted to harsher environments. Model organisms are potentially important systems for biogeographical studies, given the available range-wide natural history collections, and the importance of providing biogeographical context to their genetic and phenotypic diversity.

LOCATION

Global.

TAXON

Arabidopsis thaliana ("Arabidopsis").

METHODS

We fit occurrence records to climate data, and then projected the distribution of Arabidopsis under last glacial maximum, current, and future climates. We confronted model predictions with individual performance measured on 2,194 herbarium specimens, and we asked whether predicted suitability was associated with life-history and genomic variation measured on ~900 natural accessions.

RESULTS

The most important climate variables constraining the Arabidopsis distribution were winter cold in northern and high elevation regions and summer heat in southern regions. Herbarium specimens from regions with lower habitat suitability in both northern and southern regions were smaller, supporting the hypothesis that the distribution of Arabidopsis is constrained by climate-associated factors. Climate anomalies partly explained interannual variation in herbarium specimen size, but these did not closely correspond to local limiting factors identified in the distribution model. Late-flowering genotypes were absent from the lowest suitability regions, suggesting slower life histories are only viable closer to the center of the realized niche. We identified glacial refugia farther north than previously recognized, as well as refugia concordant with previous population genetic findings. Lower latitude populations, known to be genetically distinct, are most threatened by future climate change. The recently colonized range of Arabidopsis was well-predicted by our native-range model applied to certain regions but not others, suggesting it has colonized novel climates.

MAIN CONCLUSIONS

Integration of distribution models with performance data from vast natural history collections is a route forward for testing biogeographical hypotheses about species distributions and their relationship with evolutionary fitness across large scales.

Above and belowground phenotypic response to exogenous auxin across Arabidopsis thaliana mutants and natural accessions varies from seedling to reproductive maturity

Background

Plant hormones influence phenology, development, and function of above and belowground plant structures. In seedlings, auxin influences the initiation and development of lateral roots and root systems. How auxin-related genes influence root initiation at early life stages has been investigated from numerous perspectives. There is a gap in our understanding of how these genes influence root size through the life cycle and in mature plants. Across development, the influence of a particular gene on plant phenotypes is partly regulated by the addition of a poly-A tail to mRNA transcripts via alternative polyadenylation (APA). Auxin related genes have documented variation in APA, with auxin itself contributing to APA site switches. Studies of the influence of exogenous auxin on natural plant accessions and mutants of auxin pathway gene families exhibiting variation in APA are required for a more complete understanding of genotype by development by hormone interactions in whole plant and fitness traits.

Methods

We studied Arabidopsis thaliana homozygous mutant lines with inserts in auxin-related genes previously identified to exhibit variation in number of APA sites. Our growth chamber experiment included wildtype Col-0 controls, mutant lines, and natural accession phytometers. We applied exogenous auxin through the life cycle. We quantified belowground and aboveground phenotypes in 14 day old, 21 day old seedlings and plants at reproductive maturity. We contrasted root, rosette and flowering phenotypes across wildtype, auxin mutant, and natural accession lines, APA groups, hormone treatments, and life stages using general linear models.

Results

The root systems and rosettes of mutant lines in auxin related genes varied in response to auxin applications across life stages and varied between genotypes within life stages. In seedlings, exposure to auxin decreased size, but increased lateral root density, whereas at reproductive maturity, plants displayed greater aboveground mass and total root length. These differences may in part be due to a shift which delayed the reproductive stage when plants were treated with auxin. Root traits of auxin related mutants depended on the number of APA sites of mutant genes and the plant's developmental stage. Mutants with inserts in genes with many APA sites exhibited lower early seedling belowground biomass than those with few APA sites but only when exposed to exogenous auxin. As we observed different responses to exogenous auxin across the life cycle, we advocate for further studies of belowground traits and hormones at reproductive maturity. Studying phenotypic variation of genotypes across life stages and hormone environments will uncover additional shared patterns across traits, assisting efforts to potentially reach breeding targets and enhance our understanding of variation of genotypes in natural systems.

Pan-European study of genotypes and phenotypes in the Arabidopsis relative Cardamine hirsuta reveals how adaptation, demography, and development shape diversity patterns

We study natural DNA polymorphisms and associated phenotypes in the Arabidopsis relative Cardamine hirsuta. We observed strong genetic differentiation among several ancestry groups and broader distribution of Iberian relict strains in European C. hirsuta compared to Arabidopsis. We found synchronization between vegetative and reproductive development and a pervasive role for heterochronic pathways in shaping C. hirsuta natural variation. A single, fast-cycling ChFRIGIDA allele evolved adaptively allowing range expansion from glacial refugia, unlike Arabidopsis where multiple FRIGIDA haplotypes were involved. The Azores islands, where Arabidopsis is scarce, are a hotspot for C. hirsuta diversity. We identified a quantitative trait locus (QTL) in the heterochronic SPL9 transcription factor as a determinant of an Azorean morphotype. This QTL shows evidence for positive selection, and its distribution mirrors a climate gradient that broadly shaped the Azorean flora. Overall, we establish a framework to explore how the interplay of adaptation, demography, and development shaped diversity patterns of 2 related plant species.

Ancestral self-compatibility facilitates the establishment of allopolyploids in Brassicaceae

Self-incompatibility systems based on self-recognition evolved in hermaphroditic plants to maintain genetic variation of offspring and mitigate inbreeding depression. Despite these benefits in diploid plants, for polyploids who often face a scarcity of mating partners, self-incompatibility can thwart reproduction. In contrast, self-compatibility provides an immediate advantage: a route to reproductive viability. Thus, diploid selfing lineages may facilitate the formation of new allopolyploid species. Here, we describe the mechanism of establishment of at least four allopolyploid species in Brassicaceae (Arabidopsis suecica, Arabidopsis kamchatica, Capsella bursa-pastoris, and Brassica napus), in a manner dependent on the prior loss of the self-incompatibility mechanism in one of the ancestors. In each case, the degraded S-locus from one parental lineage was dominant over the functional S-locus of the outcrossing parental lineage. Such dominant loss-of-function mutations promote an immediate transition to selfing in allopolyploids and may facilitate their establishment.

The genome sequence of thale cress, Arabidopsis thaliana (Heynh., 1842)

We present a genome assembly of an Arabidopsis thaliana specimen (thale cress; Streptophyta; Magnoliopsida; Brassicales; Brassicaceae). The genome sequence spans 138 megabases. Most of the assembly (98.76%) is scaffolded into five chromosomal pseudomolecules. The mitochondrial and plastid genomes were also assembled and are 368.8 and 154.5 kilobases in length respectively.

Parent-of-Origin Effects on Seed Size Modify Heterosis Responses in Arabidopsis thaliana

Parent-of-origin effects arise when a phenotype depends on whether it is inherited maternally or paternally. Parent-of-origin effects can exert a strong influence on F1 seed size in flowering plants, an important agronomic and life-history trait that can contribute to biomass heterosis. Here we investigate the natural variation in the relative contributions of the maternal and paternal genomes to F1 seed size across 71 reciprocal pairs of F1 hybrid diploids and the parental effect on F1 seed size heterosis. We demonstrate that the paternally derived genome influences F1 seed size more significantly than previously appreciated. We further demonstrate (by disruption of parental genome dosage balance in F1 triploid seeds) that hybridity acts as an enhancer of genome dosage effects on F1 seed size, beyond that observed from hybridity or genome dosage effects on their own. Our findings indicate that interactions between genetic hybridity and parental genome dosage can enhance heterosis effects in plants, opening new avenues for boosting heterosis breeding in crop plants.

Multiple Sources of Introduction of North American Arabidopsis thaliana from across Eurasia

Large-scale movement of organisms across their habitable range, or migration, is an important evolutionary process that can shape genetic diversity and influence the adaptive spread of alleles. Although human migrations have been studied in great detail with modern and ancient genomes, recent anthropogenic influence on reducing the biogeographical constraints on the migration of nonnative species has presented opportunities in several study systems to ask the questions about how repeated introductions shape genetic diversity in the introduced range. We present an extensive overview of population structure of North American Arabidopsis thaliana by studying a set of 500 whole-genome sequenced and over 2,800 RAD-seq genotyped individuals in the context of global diversity represented by Afro-Eurasian genomes. We use methods based on haplotype and rare-allele sharing as well as phylogenetic modeling to identify likely sources of introductions of extant N. American A. thaliana from the native range in Africa and Eurasia. We find evidence of admixture among the introduced lineages having increased haplotype diversity and reduced mutational load. We also detect signals of selection in immune-system-related genes that may impart qualitative disease resistance to pathogens of bacterial and oomycete origin. We conclude that multiple introductions to a nonnative range can rapidly enhance the adaptive potential of a colonizing species by increasing haplotypic diversity through admixture. Our results lay the foundation for further investigations into the functional significance of admixture.

Large tandem duplications affect gene expression, 3D organization, and plant-pathogen response

Rapid plant genome evolution is crucial to adapt to environmental changes. Chromosomal rearrangements and gene copy number variation (CNV) are two important tools for genome evolution and sources for the creation of new genes. However, their emergence takes many generations. In this study, we show that in Arabidopsis thaliana, a significant loss of ribosomal RNA (rRNA) genes with a past history of a mutation for the chromatin assembly factor 1 (CAF1) complex causes rapid changes in the genome structure. Using long-read sequencing and microscopic approaches, we have identified up to 15 independent large tandem duplications in direct orientation (TDDOs) ranging from 60 kb to 1.44 Mb. Our data suggest that these TDDOs appeared within a few generations, leading to the duplication of hundreds of genes. By subsequently focusing on a line only containing 20% of rRNA gene copies (20rDNA line), we investigated the impact of TDDOs on 3D genome organization, gene expression, and cytosine methylation. We found that duplicated genes often accumulate more transcripts. Among them, several are involved in plant–pathogen response, which could explain why the 20rDNA line is hyper-resistant to both bacterial and nematode infections. Finally, we show that the TDDOs create gene fusions and/or truncations and discuss their potential implications for the evolution of plant genomes.

MSH2 shapes the meiotic crossover landscape in relation to interhomolog polymorphism in Arabidopsis

During meiosis, DNA double-strand breaks undergo interhomolog repair to yield crossovers between homologous chromosomes. To investigate how interhomolog sequence polymorphism affects crossovers, we sequenced multiple recombinant populations of the model plant Arabidopsis thaliana. Crossovers were elevated in the diverse pericentromeric regions, showing a local preference for polymorphic regions. We provide evidence that crossover association with elevated diversity is mediated via the Class I crossover formation pathway, although very high levels of diversity suppress crossovers. Interhomolog polymorphism causes mismatches in recombining molecules, which can be detected by MutS homolog (MSH) mismatch repair protein heterodimers. Therefore, we mapped crossovers in a msh2 mutant, defective in mismatch recognition, using multiple hybrid backgrounds. Although total crossover numbers were unchanged in msh2 mutants, recombination was remodelled from the diverse pericentromeres towards the less-polymorphic sub-telomeric regions. Juxtaposition of megabase heterozygous and homozygous regions causes crossover remodelling towards the heterozygous regions in wild type Arabidopsis, but not in msh2 mutants. Immunostaining showed that MSH2 protein accumulates on meiotic chromosomes during prophase I, consistent with MSH2 regulating meiotic recombination. Our results reveal a pro-crossover role for MSH2 in regions of higher sequence diversity in A. thaliana.

Natural variation at FLM splicing has pleiotropic effects modulating ecological strategies in Arabidopsis thaliana

Investigating the evolution of complex phenotypes and the underlying molecular bases of their variation is critical to understand how organisms adapt to their environment. Applying classical quantitative genetics on a segregating population derived from a Can-0xCol-0 cross, we identify the MADS-box transcription factor FLOWERING LOCUS M (FLM) as a player of the phenotypic variation in plant growth and color. We show that allelic variation at FLM modulates plant growth strategy along the leaf economics spectrum, a trade-off between resource acquisition and resource conservation, observable across thousands of plant species. Functional differences at FLM rely on a single intronic substitution, disturbing transcript splicing and leading to the accumulation of non-functional FLM transcripts. Associations between this substitution and phenotypic and climatic data across Arabidopsis natural populations, show how noncoding genetic variation at a single gene might be adaptive through pleiotropic effects.

FLOWERING LOCUS M (FLM) is known as a repressor of Arabidopsis flowering. Here, the authors show that a single intronic substitution of FLM modulates leaf color and plant growth strategy along the leaf economics spectrum, as well as plays a role in plant adaptation.

A Species-Wide Inventory of NLR Genes and Alleles in Arabidopsis thaliana

Infectious disease is both a major force of selection in nature and a prime cause of yield loss in agriculture. In plants, disease resistance is often conferred by nucleotide-binding leucine-rich repeat (NLR) proteins, intracellular immune receptors that recognize pathogen proteins and their effects on the host. Consistent with extensive balancing and positive selection, NLRs are encoded by one of the most variable gene families in plants, but the true extent of intraspecific NLR diversity has been unclear. Here, we define a nearly complete species-wide pan-NLRome in Arabidopsis thaliana based on sequence enrichment and long-read sequencing. The pan-NLRome largely saturates with approximately 40 well-chosen wild strains, with half of the pan-NLRome being present in most accessions. We chart NLR architectural diversity, identify new architectures, and quantify selective forces that act on specific NLRs and NLR domains. Our study provides a blueprint for defining pan-NLRomes.

Evolution of development of pollen performance

With the origin of pollination in ancient seed plants, the male gametophyte ("pollen") began to evolve a new and unique life history stage, the progamic phase, a post-pollination period in which pollen sexual maturation occurs in interaction with sporophyte-derived tissues. Pollen performance traits mediate the timing of the fertilization process, often in competition with other pollen, via the speed of pollen germination, sperm development, and pollen tube growth. Studies of pollen development rarely address the issue of performance or its evolution, which involves linking variation in developmental rates to relative fitness within populations or to adaptations on a macroevolutionary scale. Modifications to the pollen tube pathway and changes in the intensity of pollen competition affect the direction and strength of selection on pollen performance. Hence, pollen developmental evolution is always contextual-it involves both the population biology of pollen reaching stigmas and the co-evolution of sporophytic traits, such as the pollen tube pathway and mating system. For most species, performance evolution generally reflects a wandering history of periods of directional selection and relaxed selection, channeled by developmental limitations, a pattern that favors the accumulation of diversity and redundancy in developmental mechanisms and the genetic machinery. Developmental biologists are focused on finding universal mechanisms that underlie pollen function, and these are largely mechanisms that have evolved through their effects on performance. Here, we suggest ways in which studies of pollen performance or function could progress by cross-fertilization between the "evo" and "devo" fields.