Natural genetic variation underlying the negative effect of elevated CO2 on ionome composition in Arabidopsis thaliana.

The elevation of atmospheric CO2 leads to a decline in plant mineral content, which might pose a significant threat to food security in coming decades. Although few genes have been identified for the negative effect of elevated CO2 on plant mineral composition, several studies suggest the existence of genetic factors. Here, we performed a large-scale study to explore genetic diversity of plant ionome responses to elevated CO2, using six hundred Arabidopsis thaliana accessions, representing geographical distributions ranging from worldwide to regional and local environments. We show that growth under elevated CO2 leads to a global decrease of ionome content, whatever the geographic distribution of the population. We observed a high range of genetic diversity, ranging from the most negative effect to resilience or even to a benefit in response to elevated CO2. Using genome-wide association mapping, we identified a large set of genes associated with this response, and we demonstrated that the function of one of these genes is involved in the negative effect of elevated CO2 on plant mineral composition. This resource will contribute to understand the mechanisms underlying the effect of elevated CO2 on plant mineral nutrition, and could help towards the development of crops adapted to a high-CO2 world.

A pan-genome of 69 Arabidopsis thaliana accessions reveals a conserved genome structure throughout the global species range.

Although originally primarily a system for functional biology, Arabidopsis thaliana has, owing to its broad geographical distribution and adaptation to diverse environments, developed into a powerful model in population genomics. Here we present chromosome-level genome assemblies of 69 accessions from a global species range. We found that genomic colinearity is very conserved, even among geographically and genetically distant accessions. Along chromosome arms, megabase-scale rearrangements are rare and typically present only in a single accession. This indicates that the karyotype is quasi-fixed and that rearrangements in chromosome arms are counter-selected. Centromeric regions display higher structural dynamics, and divergences in core centromeres account for most of the genome size variations. Pan-genome analyses uncovered 32,986 distinct gene families, 60% being present in all accessions and 40% appearing to be dispensable, including 18% private to a single accession, indicating unexplored genic diversity. These 69 new Arabidopsis thaliana genome assemblies will empower future genetic research.

Genome wide association study of Arabidopsis seed mucilage layers at a regional scale.

The myxospermous species Arabidopsis thaliana extrudes a polysaccharidic mucilage from the seed coat epidermis during imbibition. The whole seed mucilage can be divided into a seed-adherent layer and a fully soluble layer, both layers presenting natural genetic variations. The adherent mucilage is variable in size and composition, while the soluble mucilage is variable in composition and physical properties. Studies reporting both the genetic architecture and the putative selective agents acting on this natural genetic variation are scarce. In this study, we set up a Genome Wide Association study (GWAS) based on 424 natural accessions collected from 166 natural populations of A. thaliana located south-west of France and previously characterized for a very important number of abiotic and biotic factors. We identified an extensive genetic variation for both mucilage layers. The adherent mucilage was mainly related to precipitation and temperature whereas the non-adherent mucilage was unrelated to any environmental factors. By combining a hierarchical Bayesian model with a local score approach, we identified 55 and 28 candidate genes, corresponding to 26 and 10 QTLs for the adherent and non-adherent mucilages, respectively. Putative or characterized function and expression data available in the literature were used to filter the candidate genes. Only one gene among our set of candidate genes was already described as a seed mucilage actor, leaving a large set of new candidates putatively implicated inseed mucilage synthesis or release. The present study lay out foundation to understand the influence of regional ecological factors acting on seed mucilage in A. thaliana.

Genetic architecture of the response of Arabidopsis thaliana to a native plant-growth-promoting bacterial strain.

By improving plant nutrition and alleviating abiotic and biotic stresses, plant growth-promoting bacteria (PGPB) can help to develop eco-friendly and sustainable agricultural practices. Besides climatic conditions, soil conditions, and microbe-microbe interactions, the host genotype influences the effectiveness of PGPB. Yet, most GWAS conducted to characterize the genetic architecture of response to PGPB are based on non-native interactions between a host plant and PGPB strains isolated from the belowground compartment of other plants. In this study, a GWAS was set up under in vitro conditions to describe the genetic architecture of the response of Arabidopsis thaliana to the PGPB Pseudomonas siliginis, by inoculating seeds of 162 natural accessions from the southwest of France with one strain isolated from the leaf compartment in the same geographical region. Strong genetic variation of plant growth response to this native PGPB was observed at a regional scale, with the strain having a positive effect on the vegetative growth of small plants and a negative effect on the vegetative growth of large plants. The polygenic genetic architecture underlying this negative trade-off showed suggestive signatures of local adaptation. The main eco-evolutionary relevant candidate genes are involved in seed and root development.

Corrigendum: Investigating genetic diversity within the most abundant and prevalent non-pathogenic leaf-associated bacteria interacting with Arabidopsis thaliana in natural habitats.

[This corrects the article DOI: 10.3389/fmicb.2022.984832.].

Epigenetic variations are more substantial than genetic variations in rapid adaptation of oyster to Pacific oyster mortality syndrome.

Disease emergence is accelerating with global changes. Understanding by which mechanisms host populations can rapidly adapt will be crucial for management practices. Pacific oyster mortality syndrome (POMS) imposes a substantial and recurrent selective pressure on oyster populations, and rapid adaptation may arise through genetics and epigenetics. In this study, we used (epi)genome-wide association mapping to show that oysters differentially exposed to POMS displayed genetic and epigenetic signatures of selection. Consistent with higher resistance to POMS, the genes targeted included many genes in several pathways related to immunity. By combining correlation, DNA methylation quantitative trait loci, and variance partitioning, we revealed that a third of phenotypic variation was explained by interactions between the genetic and epigenetic information, ~14% by the genome, and up to 25% by the epigenome alone. Similar to genetically based adaptation, epigenetic mechanisms notably governing immune responses can contribute substantially to the rapid adaptation of hosts to emerging infectious diseases.

Cycles of satellite and transposon evolution in Arabidopsis centromeres.

Centromeres are critical for cell division, loading CENH3 or CENPA histone variant nucleosomes, directing kinetochore formation and allowing chromosome segregation1,2. Despite their conserved function, centromere size and structure are diverse across species. To understand this centromere paradox3,4, it is necessary to know how centromeric diversity is generated and whether it reflects ancient trans-species variation or, instead, rapid post-speciation divergence. To address these questions, we assembled 346 centromeres from 66 Arabidopsis thaliana and 2 Arabidopsis lyrata accessions, which exhibited a remarkable degree of intra- and inter-species diversity. A. thaliana centromere repeat arrays are embedded in linkage blocks, despite ongoing internal satellite turnover, consistent with roles for unidirectional gene conversion or unequal crossover between sister chromatids in sequence diversification. Additionally, centrophilic ATHILA transposons have recently invaded the satellite arrays. To counter ATHILA invasion, chromosome-specific bursts of satellite homogenization generate higher-order repeats and purge transposons, in line with cycles of repeat evolution. Centromeric sequence changes are even more extreme in comparison between A. thaliana and A. lyrata. Together, our findings identify rapid cycles of transposon invasion and purging through satellite homogenization, which drive centromere evolution and ultimately contribute to speciation.

An atypical NLR gene confers bacterial wilt susceptibility in Arabidopsis.

Quantitative disease resistance (QDR) remains the most prevalent form of plant resistance in crop fields and wild habitats. Genome-wide association studies (GWAS) have proved to be successful in deciphering the quantitative genetic basis of complex traits such as QDR. To unravel the genetics of QDR to the devastating worldwide bacterial pathogen Ralstonia solanacearum, we performed a GWAS by challenging a highly polymorphic local mapping population of Arabidopsis thaliana with four R. solanacearum type III effector (T3E) mutants, identified as key pathogenicity determinants after a first screen on an A. thaliana core collection of 25 accessions. Although most quantitative trait loci (QTLs) were highly specific to the identity of the T3E mutant (ripAC, ripAG, ripAQ, and ripU), we finely mapped a common QTL located on a cluster of nucleotide-binding domain and leucine-rich repeat (NLR) genes that exhibited structural variation. We functionally validated one of these NLRs as a susceptibility factor in response to R. solanacearum, named it Bacterial Wilt Susceptibility 1 (BWS1), and cloned two alleles that conferred contrasting levels of QDR. Further characterization indicated that expression of BWS1 leads to suppression of immunity triggered by different R. solanacearum effectors. In addition, we showed a direct interaction between BWS1 and RipAC T3E, and BWS1 and SUPPRESSOR OF G2 ALLELE OF skp1 (SGT1b), the latter interaction being suppressed by RipAC. Together, our results highlight a putative role for BWS1 as a quantitative susceptibility factor directly targeted by the T3E RipAC, mediating negative regulation of the SGT1-dependent immune response.

The Genetic Architecture of Adaptation to Leaf and Root Bacterial Microbiota in Arabidopsis thaliana.

Understanding the role of the host genome in modulating microbiota variation is a need to shed light on the holobiont theory and overcome the current limits on the description of host-microbiota interactions at the genomic and molecular levels. However, the host genetic architecture structuring microbiota is only partly described in plants. In addition, most association genetic studies on microbiota are often carried out outside the native habitats where the host evolves and the identification of signatures of local adaptation on the candidate genes has been overlooked. To fill these gaps and dissect the genetic architecture driving adaptive plant-microbiota interactions, we adopted a genome-environment association (GEA) analysis on 141 whole-genome sequenced natural populations of Arabidopsis thaliana characterized in situ for their leaf and root bacterial communities in fall and spring, and a large range of nonmicrobial ecological factors (i.e., climate, soil, and plant communities). A much higher fraction of among-population microbiota variance was explained by the host genetics than by nonmicrobial ecological factors. Importantly, the relative importance of host genetics and nonmicrobial ecological factors in explaining the presence of particular operational taxonomic units (OTUs) differs between bacterial families and genera. In addition, the polygenic architecture of adaptation to bacterial communities was highly flexible between plant compartments and seasons. Relatedly, signatures of local adaptation were stronger on quantitative trait loci (QTLs) of the root microbiota in spring. Finally, plant immunity appears as a major source of adaptive genetic variation structuring bacterial assemblages in A. thaliana.

Genome-wide association studies in plant pathosystems: success or failure?

Harnessing natural genetic variation is an established alternative to artificial genetic variation for investigating the molecular dialog between partners in plant pathosystems. Herein, we review the successes of genome-wide association studies (GWAS) in both plants and pathogens. While GWAS in plants confirmed that the genetic architecture of disease resistance is polygenic, dynamic during the infection kinetics, and dependent on the environment, GWAS shortened the time of identification of quantitative trait loci (QTLs) and revealed both complex epistatic networks and a genetic architecture dependent upon the geographical scale. A similar picture emerges from the few GWAS in pathogens. In addition, the ever-increasing number of functionally validated QTLs has revealed new molecular plant defense mechanisms and pathogenicity determinants. Finally, we propose recommendations to better decode the disease triangle.

The ecologically relevant genetics of plant-plant interactions.

Interactions among plants have been long recognized as a major force driving plant community dynamics and crop yield. Surprisingly, our knowledge of the ecological genetics associated with variation of plant-plant interactions remains limited. In this opinion article by scientists from complementary disciplines, the international PLANTCOM network identified four timely questions to foster a better understanding of the mechanisms mediating plant assemblages. We propose that by identifying the key relationships among phenotypic traits involved in plant-plant interactions and the underlying adaptive genetic and molecular pathways, while considering environmental fluctuations at diverse spatial and time scales, we can improve predictions of genotype-by-genotype-by-environment interactions and modeling of productive and stable plant assemblages in wild habitats and crop fields.

Investigating genetic diversity within the most abundant and prevalent non-pathogenic leaf-associated bacteria interacting with Arabidopsis thaliana in natural habitats.

Microbiota modulates plant health and appears as a promising lever to develop innovative, sustainable and eco-friendly agro-ecosystems. Key patterns of microbiota assemblages in plants have been revealed by an extensive number of studies based on taxonomic profiling by metabarcoding. However, understanding the functionality of microbiota is still in its infancy and relies on reductionist approaches primarily based on the establishment of representative microbial collections. In Arabidopsis thaliana, most of these microbial collections include one strain per OTU isolated from a limited number of habitats, thereby neglecting the ecological potential of genetic diversity within microbial species. With this study, we aimed at estimating the extent of genetic variation between strains within the most abundant and prevalent leaf-associated non-pathogenic bacterial species in A. thaliana located south-west of France. By combining a culture-based collection approach consisting of the isolation of more than 7,000 bacterial colonies with an informative-driven approach, we isolated 35 pure strains from eight non-pathogenic bacterial species. We detected significant intra-specific genetic variation at the genomic level and for growth rate in synthetic media. In addition, significant host genetic variation was detected in response to most bacterial strains in in vitro conditions, albeit dependent on the developmental stage at which plants were inoculated, with the presence of both negative and positive responses on plant growth. Our study provides new genetic and genomic resources for a better understanding of the plant-microbe ecological interactions at the microbiota level. We also highlight the need of considering genetic variation in both non-pathogenic bacterial species and A. thaliana to decipher the genetic and molecular mechanisms involved in the ecologically relevant dialog between hosts and leaf microbiota.

A Genome-Wide Association study in Arabidopsis thaliana to decipher the adaptive genetics of quantitative disease resistance in a native heterogeneous environment.

Pathogens are often the main selective agents acting in plant communities, thereby influencing the distribution of polymorphism at loci affecting resistance within and among natural plant populations. In addition, the outcome of plant-pathogen interactions can be drastically affected by abiotic and biotic factors at different spatial and temporal grains. The characterization of the adaptive genetic architecture of disease resistance in native heterogeneous environments is however still missing. In this study, we conducted an in situ Genome-Wide Association study in the spatially heterogeneous native habitat of a highly genetically polymorphic local mapping population of Arabidopsis thaliana, to unravel the adaptive genetic architecture of quantitative disease resistance. Disease resistance largely differed among three native soils and was affected by the presence of the grass Poa annua. The observation of strong crossing reactions norms among the 195 A. thaliana genotypes for disease resistance among micro-habitats, combined with a negative fecundity-disease resistance relationship in each micro-habitat, suggest that alternative local genotypes of A. thaliana are favored under contrasting environmental conditions at the scale of few meters. A complex genetic architecture was detected for disease resistance and fecundity. However, only few QTLs were common between these two traits. Heterogeneous selection in this local population should therefore promote the maintenance of polymorphism at only few candidate resistance genes.

APOK3, a pollen killer antidote in Arabidopsis thaliana.

The principles of heredity state that the two alleles carried by a heterozygote are equally transmitted to the progeny. However, genomic regions that escape this rule have been reported in many organisms. It is notably the case of genetic loci referred to as gamete killers, where one allele enhances its transmission by causing the death of the gametes that do not carry it. Gamete killers are of great interest, particularly to understand mechanisms of evolution and speciation. Although being common in plants, only a few, all in rice, have so far been deciphered to the causal genes. Here, we studied a pollen killer found in hybrids between two accessions of Arabidopsis thaliana. Exploring natural variation, we observed this pollen killer in many crosses within the species. Genetic analyses revealed that three genetically linked elements are necessary for pollen killer activity. Using mutants, we showed that this pollen killer works according to a poison-antidote model, where the poison kills pollen grains not producing the antidote. We identified the gene encoding the antidote, a chimeric protein addressed to mitochondria. De novo genomic sequencing in 12 natural variants with different behaviors regarding the pollen killer revealed a hyper variable locus, with important structural variations particularly in killer genotypes, where the antidote gene recently underwent duplications. Our results strongly suggest that the gene has newly evolved within A. thaliana. Finally, we identified in the protein sequence polymorphisms related to its antidote activity.

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

A paradoxical finding from genome-wide association studies (GWAS) in plants is that variation in metabolite profiles typically maps to a small number of loci, despite the complexity of underlying biosynthetic pathways. This discrepancy may partially arise from limitations presented by geographically diverse mapping panels. Properties of metabolic pathways that impede GWAS by diluting the additive effect of a causal variant, such as allelic and genetic heterogeneity and epistasis, would be expected to increase in severity with the geographical range of the mapping panel. We hypothesized that a population from a single locality would reveal an expanded set of associated loci. We tested this in a French Arabidopsis thaliana population (less than 1 km transect) by profiling and conducting GWAS for glucosinolates, a suite of defensive metabolites that have been studied in depth through functional and genetic mapping approaches. For two distinct classes of glucosinolates, we discovered more associations at biosynthetic loci than the previous GWAS with continental-scale mapping panels. Candidate genes underlying novel associations were supported by concordance between their observed effects in the TOU-A population and previous functional genetic and biochemical characterization. Local populations complement geographically diverse mapping panels to reveal a more complete genetic architecture for metabolic traits. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.

Molecular and Genomic Characterization of the Pseudomonas syringae Phylogroup 4: An Emerging Pathogen of Arabidopsis thaliana and Nicotiana benthamiana.

Environmental fluctuations such as increased temperature, water availability, and air CO2 concentration triggered by climate change influence plant disease dynamics by affecting hosts, pathogens, and their interactions. Here, we describe a newly discovered Pseudomonas syringae strain found in a natural population of Arabidopsis thaliana collected from the southwest of France. This strain, called Psy RAYR-BL, is highly virulent on natural Arabidopsis accessions, Arabidopsis model accession Columbia 0, and tobacco plants. Despite the severe disease phenotype caused by the Psy RAYR-BL strain, we identified a reduced repertoire of putative Type III virulence effectors by genomic sequencing compared to P. syringae pv tomato (Pst) DC3000. Furthermore, hopBJ1Psy is found exclusively on the Psy RAYR-BL genome but not in the Pst DC3000 genome. The plant expression of HopBJ1Psy induces ROS accumulation and cell death. In addition, HopBJ1Psy participates as a virulence factor in this plant-pathogen interaction, likely explaining the severity of the disease symptoms. This research describes the characterization of a newly discovered plant pathogen strain and possible virulence mechanisms underlying the infection process shaped by natural and changing environmental conditions.

Parallel reduction in flowering time from de novo mutations enable evolutionary rescue in colonizing lineages.

Understanding how populations adapt to abrupt environmental change is necessary to predict responses to future challenges, but identifying specific adaptive variants, quantifying their responses to selection and reconstructing their detailed histories is challenging in natural populations. Here, we use Arabidopsis from the Cape Verde Islands as a model to investigate the mechanisms of adaptation after a sudden shift to a more arid climate. We find genome-wide evidence of adaptation after a multivariate change in selection pressures. In particular, time to flowering is reduced in parallel across islands, substantially increasing fitness. This change is mediated by convergent de novo loss of function of two core flowering time genes: FRI on one island and FLC on the other. Evolutionary reconstructions reveal a case where expansion of the new populations coincided with the emergence and proliferation of these variants, consistent with models of rapid adaptation and evolutionary rescue.

Draft Genome Sequence of Pseudomonas syringae RAYR-BL, a Strain Isolated from Natural Accessions of Arabidopsis thaliana Plants.

Here, we report the genome sequence of the P. syringae strain RAYR-BL, isolated from natural accessions of Arabidopsis plants. The draft genome sequence consists of 5.85 Mbp assembled in 110 contigs. The study of P. syringae RAYR-BL is a valuable tool to investigate molecular features of plant-pathogen interaction under environmental conditions.

Study of natural diversity in response to a key pathogenicity regulator of Ralstonia solanacearum reveals new susceptibility genes in Arabidopsis thaliana.

Ralstonia solanacearum gram-negative phytopathogenic bacterium exerts its virulence through a type III secretion system (T3SS) that translocates type III effectors (T3Es) directly into the host cells. T3E secretion is finely controlled at the posttranslational level by helper proteins, T3SS control proteins, and type III chaperones. The HpaP protein, one of the type III secretion substrate specificity switch (T3S4) proteins, was previously highlighted as a virulence factor on Arabidopsis thaliana Col-0 accession. In this study, we set up a genome-wide association analysis to explore the natural diversity of response to the hpaP mutant of two A. thaliana mapping populations: a worldwide collection and a local population. Quantitative genetic variation revealed different genetic architectures in both mapping populations, with a global delayed response to the hpaP mutant compared to the GMI1000 wild-type strain. We have identified several quantitative trait loci (QTLs) associated with the hpaP mutant inoculation. The genes underlying these QTLs are involved in different and specific biological processes, some of which were demonstrated important for R. solanacearum virulence. We focused our study on four candidate genes, RKL1, IRE3, RACK1B, and PEX3, identified using the worldwide collection, and validated three of them as susceptibility factors. Our findings demonstrate that the study of the natural diversity of plant response to a R. solanacearum mutant in a key regulator of virulence is an original and powerful strategy to identify genes directly or indirectly targeted by the pathogen.

The Genomic Architecture of Competitive Response of Arabidopsis thaliana Is Highly Flexible Among Plurispecific Neighborhoods.

Plants are daily challenged by multiple abiotic and biotic stresses. A major biotic constraint corresponds to competition with other plant species. Although plants simultaneously interact with multiple neighboring species throughout their life cycle, there is still very limited information about the genetics of the competitive response in the context of plurispecific interactions. Using a local mapping population of Arabidopsis thaliana, we set up a genome wide association study (GWAS) to estimate the extent of genetic variation of competitive response in 12 plant species assemblages, based on three competitor species (Poa annua, Stellaria media, and Veronica arvensis). Based on five phenotypic traits, we detected strong crossing reaction norms not only between the three bispecific neighborhoods but also among the plurispecific neighborhoods. The genetic architecture of competitive response was highly dependent on the identity and the relative abundance of the neighboring species. In addition, most of the enriched biological processes underlying competitive responses largely differ among neighborhoods. While the RNA related processes might confer a broad range response toolkit for multiple traits in diverse neighborhoods, some processes, such as signaling and transport, might play a specific role in particular assemblages. Altogether, our results suggest that plants can integrate and respond to different species assemblages depending on the identity and number of each neighboring species, through a large range of candidate genes associated with diverse and unexpected processes leading to developmental and stress responses.