Genome-wide association mapping within a local Arabidopsis thaliana population more fully reveals the genetic architecture for defensive metabolite diversity

A Genome-Wide Association Screen for Genes Affecting Leaf Trichome Development and Epidermal Metal Accumulation in Arabidopsis

To identify novel genes engaged in plant epidermal development, we characterized the phenotypic variability of rosette leaf epidermis of 310 sequenced Arabidopsis thaliana accessions, focusing on trichome shape and distribution, compositional characteristics of the trichome cell wall, and histologically detectable metal ion distribution. Some of these traits correlated with cLimate parameters of our accession's locations of origin, suggesting environmental selection. A novel metal deposition pattern in stomatal guard cells was observed in some accessions. Subsequent GWAS analysis identified 1546 loci with protein sequence-altering SNPs associated with one or more traits, including 5 genes with previously reported relevant mutant phenotypes and 80 additional genes with known or predicted roles in relevant developmental and cellular processes. Some candidates, including GFS9/TT9, exhibited environmentally correlated allele distribution. Several large gene famiLies, namely DUF674, DUF784, DUF1262, DUF1985, DUF3741, cytochrome P450, receptor-Like kinases, Cys/His-rich C1 domain proteins and formins were overrepresented among the candidates for various traits, suggesting epidermal development-related functions. A possible participation of formins in guard cell metal deposition was supported by observations in available loss of function mutants. Screening of candidate gene lists against the STRING interactome database uncovered several predominantly nuclear protein interaction networks with possible novel roles in epidermal development.

Study design and the sampling of deleterious rare variants in biobank-scale datasets

One key component of study design in population genetics is the "geographic breadth" of a sample (i.e., how broad a region across which individuals are sampled). How the geographic breadth of a sample impacts observations of rare, deleterious variants is unclear, even though such variants are of particular interest for biomedical and evolutionary applications. Here, in order to gain insight into the effects of sample design on ascertained genetic variants, we formulate a stochastic model of dispersal, genetic drift, selection, mutation, and geographically concentrated sampling. We use this model to understand the effects of the geographic breadth of sampling effort on the discovery of negatively selected variants. We find that samples which are more geographically broad will discover a greater number variants as compared geographically narrow samples (an effect we label "discovery"); though the variants will be detected at lower average frequency than in narrow samples (e.g. as singletons, an effect we label "dilution"). Importantly, these effects are amplified for larger sample sizes and moderated by the magnitude of fitness effects. We validate these results using both population genetic simulations and empirical analyses in the UK Biobank. Our results are particularly important in two contexts: the association of large-effect rare variants with particular phenotypes and the inference of negative selection from allele frequency data. Overall, our findings emphasize the importance of considering geographic breadth when designing and carrying out genetic studies, especially at biobank scale.

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.

Inference of Host-Pathogen Interaction Matrices from Genome-Wide Polymorphism Data

Host-pathogen coevolution is defined as the reciprocal evolutionary changes in both species due to genotype × genotype (G×G) interactions at the genetic level determining the outcome and severity of infection. While co-analyses of hosts and pathogen genomes (co-genome-wide association studies) allow us to pinpoint the interacting genes, these do not reveal which host genotype(s) is/are resistant to which pathogen genotype(s). The knowledge of this so-called infection matrix is important for agriculture and medicine. Building on established theories of host-pathogen interactions, we here derive four novel indices capturing the characteristics of the infection matrix. These indices can be computed from full genome polymorphism data of randomly sampled uninfected hosts, as well as infected hosts and their pathogen strains. We use these indices in an approximate Bayesian computation method to pinpoint loci with relevant G×G interactions and to infer their underlying interaction matrix. In a combined single nucleotide polymorphism dataset of 451 European humans and their infecting hepatitis C virus (HCV) strains and 503 uninfected individuals, we reveal a new human candidate gene for resistance to HCV and new virus mutations matching human genes. For two groups of significant human-HCV (G×G) associations, we infer a gene-for-gene infection matrix, which is commonly assumed to be typical of plant-pathogen interactions. Our model-based inference framework bridges theoretical models of G×G interactions with host and pathogen genomic data. It, therefore, paves the way for understanding the evolution of key G×G interactions underpinning HCV adaptation to the European human population after a recent expansion.

Polygenic architecture of flowering time and its relationship with local environments in the grass Brachypodium distachyon

Synchronizing the timing of reproduction with the environment is crucial in the wild. Among the multiple mechanisms, annual plants evolved to sense their environment, the requirement of cold-mediated vernalization is a major process that prevents individuals from flowering during winter. In many annual plants including crops, both a long and short vernalization requirement can be observed within species, resulting in so-called early-(spring) and late-(winter) flowering genotypes. Here, using the grass model Brachypodium distachyon, we explored the link between flowering-time-related traits (vernalization requirement and flowering time), environmental variation, and diversity at flowering-time genes by combining measurements under greenhouse and outdoor conditions. These experiments confirmed that B. distachyon natural accessions display large differences regarding vernalization requirements and ultimately flowering time. We underline significant, albeit quantitative effects of current environmental conditions on flowering-time-related traits. While disentangling the confounding effects of population structure on flowering-time-related traits remains challenging, population genomics analyses indicate that well-characterized flowering-time genes may contribute significantly to flowering-time variation and display signs of polygenic selection. Flowering-time genes, however, do not colocalize with genome-wide association peaks obtained with outdoor measurements, suggesting that additional genetic factors contribute to flowering-time variation in the wild. Altogether, our study fosters our understanding of the polygenic architecture of flowering time in a natural grass system and opens new avenues of research to investigate the gene-by-environment interaction at play for this trait.

A diverse panel of 755 bread wheat accessions harbors untapped genetic diversity in landraces and reveals novel genetic regions conferring powdery mildew resistance

KEY MESSAGE

A bread wheat panel reveals rich genetic diversity in Turkish, Pakistani and Iranian landraces and novel resistance loci to diverse powdery mildew isolates via subsetting approaches in association studies. Wheat breeding for disease resistance relies on the availability and use of diverse genetic resources. More than 800,000 wheat accessions are globally conserved in gene banks, but they are mostly uncharacterized for the presence of resistance genes and their potential for agriculture. Based on the selective reduction of previously assembled collections for allele mining for disease resistance, we assembled a trait-customized panel of 755 geographically diverse bread wheat accessions with a focus on landraces, called the LandracePLUS panel. Population structure analysis of this panel based on the TaBW35K SNP array revealed an increased genetic diversity compared to 632 landraces genotyped in an earlier study and 17 high-quality sequenced wheat accessions. The additional genetic diversity found here mostly originated from Turkish, Iranian and Pakistani landraces. We characterized the LandracePLUS panel for resistance to ten diverse isolates of the fungal pathogen powdery mildew. Performing genome-wide association studies and dividing the panel further by a targeted subsetting approach for accessions of distinct geographical origin, we detected several known and already cloned genes, including the Pm2a gene. In addition, we identified 22 putatively novel powdery mildew resistance loci that represent useful sources for resistance breeding and for research on the mildew-wheat pathosystem. Our study shows the value of assembling trait-customized collections and utilizing a diverse range of pathogen races to detect novel loci. It further highlights the importance of integrating landraces of different geographical origins into future diversity studies.

SHI family transcription factors regulate an interspecific barrier

Pre-zygotic interspecies incompatibility in angiosperms is an important mechanism to prevent unfavourable hybrids between species. Here we report our identification of STIGMATIC PRIVACY 2 (SPRI2), a transcription factor that has a zinc-finger domain and regulates interspecies barriers in Arabidopsis thaliana, via genome-wide association study. Knockout analysis of SPRI2/SRS7 and its paralogue SPRI2-like/SRS5 demonstrated their necessity in rejecting male pollen from other species within female pistils. Additionally, they govern mRNA transcription of xylan O-acetyltransferases (TBL45 and TBL40) related to cell wall modification, alongside SPRI1, a pivotal transmembrane protein for interspecific pollen rejection. SPRI2/SRS7 is localized as condensed structures in the nucleus formed via liquid-liquid phase separation (LLPS), and a prion-like sequence in its amino-terminal region was found to be responsible for the formation of the condensates. The LLPS-regulated SPRI2/SRS7 discovered in this study may contribute to the establishment of interspecific reproductive barriers through the transcriptional regulation of cell wall modification genes and SPRI1.

Indirect genetic effects are shaped by demographic history and ecology in Arabidopsis thaliana

The phenotype of an individual can be affected by the genes of its conspecifics through indirect genetic effects (IGEs). IGEs have been studied across different organisms including wild and domesticated animals and plants, but little is known about their genetic architecture. Here, in a large-scale intraspecific interaction experiment, we show that the contribution of IGEs to the biomass variation of Arabidopsis thaliana is comparable to values classically reported in animals. Moreover, we identify 11 loci explaining 85.1% of the variability in IGEs. We find that positive IGE alleles (that is, those with positive effects on neighbour biomass) occur both in relict accessions from southern Eurasia and in post-glacial colonizers from northern Scandinavia, and that they are likely to have two divergent origins: for nine loci, they evolved in the post-glacial colonizers independently from the relicts, while the two others were introgressed in the post-glacial colonizer from the relicts. Finally, we find that variation in IGEs probably reflects divergent adaptations to the contrasting environments of the edges and the centre of the native range of the species. These findings reveal a surprisingly tractable genetic basis of IGEs in A. thaliana that is shaped by the ecology and the demographic history of the species.

The design of mapping populations: Impacts of geographic scale on genetic architecture and mapping efficacy for defense and immunity

Genome-wide association studies (GWAS) have yielded tremendous insight into the genetic architecture of trait variation. However, the collections of loci they uncover are far from exhaustive. As many of the complicating factors that confound or limit the efficacy of GWAS are exaggerated over broad geographic scales, a shift toward more analyses using mapping panels sampled from narrow geographic localities ("local" populations) could provide novel, complementary insights. Here, we present an overview of the major complicating factors, review mounting evidence from genomic analyses that these factors are pervasive, and synthesize theoretical and empirical evidence for the power of GWAS in local populations.

Genetic basis of speciation and adaptation: from loci to causative mutations

Does evolution proceed in small steps or large leaps? How repeatable is evolution? How constrained is the evolutionary process? Answering these long-standing questions in evolutionary biology is indispensable for both understanding how extant biodiversity has evolved and predicting how organisms and ecosystems will respond to changing environments in the future. Understanding the genetic basis of phenotypic diversification and speciation in natural populations is key to properly answering these questions. The leap forward in genome sequencing technologies has made it increasingly easier to not only investigate the genetic architecture but also identify the variant sites underlying adaptation and speciation in natural populations. Furthermore, recent advances in genome editing technologies are making it possible to investigate the functions of each candidate gene in organisms from natural populations. In this article, we discuss how these recent technological advances enable the analysis of causative genes and mutations and how such analysis can help answer long-standing evolutionary biology questions.

This article is part of the theme issue ‘Genetic basis of adaptation and speciation: from loci to causative mutations’.