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.

Epigenomic Diversity in a Global Collection of Arabidopsis thaliana Accessions.

The epigenome orchestrates genome accessibility, functionality, and three-dimensional structure. Because epigenetic variation can impact transcription and thus phenotypes, it may contribute to adaptation. Here, we report 1,107 high-quality single-base resolution methylomes and 1,203 transcriptomes from the 1001 Genomes collection of Arabidopsis thaliana. Although the genetic basis of methylation variation is highly complex, geographic origin is a major predictor of genome-wide DNA methylation levels and of altered gene expression caused by epialleles. Comparison to cistrome and epicistrome datasets identifies associations between transcription factor binding sites, methylation, nucleotide variation, and co-expression modules. Physical maps for nine of the most diverse genomes reveal how transposons and other structural variants shape the epigenome, with dramatic effects on immunity genes. The 1001 Epigenomes Project provides a comprehensive resource for understanding how variation in DNA methylation contributes to molecular and non-molecular phenotypes in natural populations of the most studied model plant.

Small proteins modulate ion-channel-like ACD6 to regulate immunity in Arabidopsis thaliana.

The multi-pass transmembrane protein ACCELERATED CELL DEATH 6 (ACD6) is an immune regulator in Arabidopsis thaliana with an unclear biochemical mode of action. We have identified two loci, MODULATOR OF HYPERACTIVE ACD6 1 (MHA1) and its paralog MHA1-LIKE (MHA1L), that code for ∼7 kDa proteins, which differentially interact with specific ACD6 variants. MHA1L enhances the accumulation of an ACD6 complex, thereby increasing the activity of the ACD6 standard allele for regulating plant growth and defenses. The intracellular ankyrin repeats of ACD6 are structurally similar to those found in mammalian ion channels. Several lines of evidence link increased ACD6 activity to enhanced calcium influx, with MHA1L as a direct regulator of ACD6, indicating that peptide-regulated ion channels are not restricted to animals.

Species-wide genetic incompatibility analysis identifies immune genes as hot spots of deleterious epistasis.

Intraspecific genetic incompatibilities prevent the assembly of specific alleles into single genotypes and influence genome- and species-wide patterns of sequence variation. A common incompatibility in plants is hybrid necrosis, characterized by autoimmune responses due to epistatic interactions between natural genetic variants. By systematically testing thousands of F1 hybrids of Arabidopsis thaliana strains, we identified a small number of incompatibility hot spots in the genome, often in regions densely populated by nucleotide-binding domain and leucine-rich repeat (NLR) immune receptor genes. In several cases, these immune receptor loci interact with each other, suggestive of conflict within the immune system. A particularly dangerous locus is a highly variable cluster of NLR genes, DM2, which causes multiple independent incompatibilities with genes that encode a range of biochemical functions, including NLRs. Our findings suggest that deleterious interactions of immune receptors limit the combinations of favorable disease resistance alleles accessible to plant genomes.

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.

Mutation bias reflects natural selection in Arabidopsis thaliana.

Since the first half of the twentieth century, evolutionary theory has been dominated by the idea that mutations occur randomly with respect to their consequences1. Here we test this assumption with large surveys of de novo mutations in the plant Arabidopsis thaliana. In contrast to expectations, we find that mutations occur less often in functionally constrained regions of the genome-mutation frequency is reduced by half inside gene bodies and by two-thirds in essential genes. With independent genomic mutation datasets, including from the largest Arabidopsis mutation accumulation experiment conducted to date, we demonstrate that epigenomic and physical features explain over 90% of variance in the genome-wide pattern of mutation bias surrounding genes. Observed mutation frequencies around genes in turn accurately predict patterns of genetic polymorphisms in natural Arabidopsis accessions (r = 0.96). That mutation bias is the primary force behind patterns of sequence evolution around genes in natural accessions is supported by analyses of allele frequencies. Finally, we find that genes subject to stronger purifying selection have a lower mutation rate. We conclude that epigenome-associated mutation bias2 reduces the occurrence of deleterious mutations in Arabidopsis, challenging the prevailing paradigm that mutation is a directionless force in evolution.

Building pangenome graphs.

Pangenome graphs can represent all variation between multiple reference genomes, but current approaches to build them exclude complex sequences or are based upon a single reference. In response, we developed the PanGenome Graph Builder, a pipeline for constructing pangenome graphs without bias or exclusion. The PanGenome Graph Builder uses all-to-all alignments to build a variation graph in which we can identify variation, measure conservation, detect recombination events and infer phylogenetic relationships.

Fast-forward genetics identifies plant CPL phosphatases as regulators of miRNA processing factor HYL1.

MicroRNAs (miRNAs) are processed from primary transcripts that contain partially self-complementary foldbacks. As in animals, the core microprocessor in plants is a Dicer protein, DICER-LIKE1 (DCL1). Processing accuracy and strand selection is greatly enhanced through the RNA binding protein HYPONASTIC LEAVES 1 (HYL1) and the zinc finger protein SERRATE (SE). We have combined a luciferase-based genetic screen with whole-genome sequencing for rapid identification of new regulators of miRNA biogenesis and action. Among the first six mutants analyzed were three alleles of C-TERMINAL DOMAIN PHOSPHATASE-LIKE 1 (CPL1)/FIERY2 (FRY2). In the miRNA processing complex, SE functions as a scaffold to mediate CPL1 interaction with HYL1, which needs to be dephosphorylated for optimal activity. In the absence of CPL1, HYL1 dephosphorylation and hence accurate processing and strand selection from miRNA duplexes are compromised. Our findings thus define a new regulatory step in plant miRNA biogenesis.

The EDS1-PAD4-ADR1 node mediates Arabidopsis pattern-triggered immunity.

Plants deploy cell-surface and intracellular leucine rich-repeat domain (LRR) immune receptors to detect pathogens1. LRR receptor kinases and LRR receptor proteins at the plasma membrane recognize microorganism-derived molecules to elicit pattern-triggered immunity (PTI), whereas nucleotide-binding LRR proteins detect microbial effectors inside cells to confer effector-triggered immunity (ETI). Although PTI and ETI are initiated in different host cell compartments, they rely on the transcriptional activation of similar sets of genes2, suggesting pathway convergence upstream of nuclear events. Here we report that PTI triggered by the Arabidopsis LRR receptor protein RLP23 requires signalling-competent dimers of the lipase-like proteins EDS1 and PAD4, and of ADR1 family helper nucleotide-binding LRRs, which are all components of ETI. The cell-surface LRR receptor kinase SOBIR1 links RLP23 with EDS1, PAD4 and ADR1 proteins, suggesting the formation of supramolecular complexes containing PTI receptors and transducers at the inner side of the plasma membrane. We detected similar evolutionary patterns in LRR receptor protein and nucleotide-binding LRR genes across Arabidopsis accessions; overall higher levels of variation in LRR receptor proteins than in LRR receptor kinases are consistent with distinct roles of these two receptor families in plant immunity. We propose that the EDS1-PAD4-ADR1 node is a convergence point for defence signalling cascades, activated by both surface-resident and intracellular LRR receptors, in conferring pathogen immunity.

Transposon dynamics in the emerging oilseed crop Thlaspi arvense.

Genome evolution is partly driven by the mobility of transposable elements (TEs) which often leads to deleterious effects, but their activity can also facilitate genetic novelty and catalyze local adaptation. We explored how the intraspecific diversity of TE polymorphisms might contribute to the broad geographic success and adaptive capacity of the emerging oil crop Thlaspi arvense (field pennycress). We classified the TE inventory based on a high-quality genome assembly, estimated the age of retrotransposon TE families and comprehensively assessed their mobilization potential. A survey of 280 accessions from 12 regions across the Northern hemisphere allowed us to quantify over 90,000 TE insertion polymorphisms (TIPs). Their distribution mirrored the genetic differentiation as measured by single nucleotide polymorphisms (SNPs). The number and types of mobile TE families vary substantially across populations, but there are also shared patterns common to all accessions. Ty3/Athila elements are the main drivers of TE diversity in T. arvense populations, while a single Ty1/Alesia lineage might be particularly important for transcriptome divergence. The number of retrotransposon TIPs is associated with variation at genes related to epigenetic regulation, including an apparent knockout mutation in BROMODOMAIN AND ATPase DOMAIN-CONTAINING PROTEIN 1 (BRAT1), while DNA transposons are associated with variation at the HSP19 heat shock protein gene. We propose that the high rate of mobilization activity can be harnessed for targeted gene expression diversification, which may ultimately present a toolbox for the potential use of transposition in breeding and domestication of T. arvense.

Activation of an atypical plant NLR with an N-terminal deletion initiates cell death at the vacuole.

Plants evolve nucleotide-binding leucine-rich repeat receptors (NLRs) to induce immunity. Activated coiled-coil (CC) domain containing NLRs (CNLs) oligomerize and form apparent cation channels promoting calcium influx and cell death, with the alpha-1 helix of the individual CC domains penetrating the plasma membranes. Some CNLs are characterized by putative N-myristoylation and S-acylation sites in their CC domain, potentially mediating permanent membrane association. Whether activated Potentially Membrane Localized NLRs (PMLs) mediate cell death and calcium influx in a similar way is unknown. We uncovered the cell-death function at the vacuole of an atypical but conserved Arabidopsis PML, PML5, which has a significant deletion in its CCG10/GA domain. Active PML5 oligomers localize in Golgi membranes and the tonoplast, alter vacuolar morphology, and induce cell death, with the short N-terminus being sufficient. Mutant analysis supports a potential role of PMLs in plant immunity. PML5-like deletions are found in several Brassicales paralogs, pointing to the evolutionary importance of this innovation. PML5, with its minimal CC domain, represents the first identified CNL utilizing vacuolar-stored calcium for cell death induction.

Argonaute10 as a miRNA locker.

The stability and translation efficiency of many messenger RNAs is regulated by microRNAs (miRNAs), which exert their effects through associated Argonaute proteins. In this issue, Zhu, Zhang, and colleagues reveal that plants also exploit miRNA binding by Argonautes as a sequestering mechanism that prevents miRNAs from fulfilling their normal roles.

Standing genetic variation fuels rapid evolution of herbicide resistance in blackgrass.

Repeated herbicide applications in agricultural fields exert strong selection on weeds such as blackgrass (Alopecurus myosuroides), which is a major threat for temperate climate cereal crops. This inadvertent selection pressure provides an opportunity for investigating the underlying genetic mechanisms and evolutionary processes of rapid adaptation, which can occur both through mutations in the direct targets of herbicides and through changes in other, often metabolic, pathways, known as non-target-site resistance. How much target-site resistance (TSR) relies on de novo mutations vs. standing variation is important for developing strategies to manage herbicide resistance. We first generated a chromosome-level reference genome for A. myosuroides for population genomic studies of herbicide resistance and genome-wide diversity across Europe in this species. Next, through empirical data in the form of highly accurate long-read amplicons of alleles encoding acetyl-CoA carboxylase (ACCase) and acetolactate synthase (ALS) variants, we showed that most populations with resistance due to TSR mutations-23 out of 27 and six out of nine populations for ACCase and ALS, respectively-contained at least two TSR haplotypes, indicating that soft sweeps are the norm. Finally, through forward-in-time simulations, we inferred that TSR is likely to mainly result from standing genetic variation, with only a minor role for de novo mutations.

Transcriptional control of gene expression by microRNAs.

MicroRNAs (miRNAs) control gene expression in animals and plants. Like another class of small RNAs, siRNAs, they affect gene expression posttranscriptionally. While siRNAs in addition act in transcriptional gene silencing, a role of miRNAs in transcriptional regulation has been less clear. We show here that in moss Physcomitrella patens mutants without a DICER-LIKE1b gene, maturation of miRNAs is normal but cleavage of target RNAs is abolished and levels of these transcripts are drastically reduced. These mutants accumulate miRNA:target-RNA duplexes and show hypermethylation of the genes encoding target RNAs, leading to gene silencing. This pathway occurs also in the wild-type upon hormone treatment. We propose that initiation of epigenetic silencing by DNA methylation depends on the ratio of the miRNA and its target RNA.

Publisher Correction: Natural selection on the Arabidopsis thaliana genome in present and future climates.

An Amendment to this paper has been published and can be accessed via a link at the top of the paper.

Natural selection on the Arabidopsis thaliana genome in present and future climates.

Through the lens of evolution, climate change is an agent of natural selection that forces populations to change and adapt, or face extinction. However, current assessments of the risk of biodiversity associated with climate change1 do not typically take into account how natural selection influences populations differently depending on their genetic makeup2. Here we make use of the extensive genome information that is available for Arabidopsis thaliana and measure how manipulation of the amount of rainfall affected the fitness of 517 natural Arabidopsis lines that were grown in Spain and Germany. This allowed us to directly infer selection along the genome3. Natural selection was particularly strong in the hot-dry location in Spain, where 63% of lines were killed and where natural selection substantially changed the frequency of approximately 5% of all genome-wide variants. A significant portion of this climate-driven natural selection of variants was predictable from signatures of local adaptation (R2 = 29-52%), as genetic variants that were found in geographical areas with climates more similar to the experimental sites were positively selected. Field-validated predictions across the species range indicated that Mediterranean and western Siberian populations-at the edges of the environmental limits of this species-currently experience the strongest climate-driven selection. With more frequent droughts and rising temperatures in Europe4, we forecast an increase in directional natural selection moving northwards from the southern end of Europe, putting many native A. thaliana populations at evolutionary risk.

Arabidopsis FLL2 promotes liquid-liquid phase separation of polyadenylation complexes.

An important component of cellular biochemistry is the concentration of proteins and nucleic acids in non-membranous compartments1,2. These biomolecular condensates are formed from processes that include liquid-liquid phase separation. The multivalent interactions necessary for liquid-liquid phase separation have been extensively studied in vitro1,3. However, the regulation of this process in vivo is poorly understood. Here we identify an in vivo regulator of liquid-liquid phase separation through a genetic screen targeting factors required for Arabidopsis RNA-binding protein FCA function. FCA contains prion-like domains that phase-separate in vitro, and exhibits behaviour in vivo that is consistent with phase separation. The mutant screen identified a functional requirement for FLL2, a coiled-coil protein, in the formation of FCA nuclear bodies. FCA reduces transcriptional read-through by promoting proximal polyadenylation at many sites in the Arabidopsis genome3,4. FLL2 was required to promote this proximal polyadenylation, but not the binding of FCA to target RNA. Ectopic expression of FLL2 increased the size and number of FCA nuclear bodies. Crosslinking with formaldehyde captured in vivo interactions between FLL2, FCA and the polymerase and nuclease modules of the RNA 3'-end processing machinery. These 3' RNA-processing components colocalized with FCA in the nuclear bodies in vivo, which indicates that FCA nuclear bodies compartmentalize 3'-end processing factors to enhance polyadenylation at specific sites. Our findings show that coiled-coil proteins can promote liquid-liquid phase separation, which expands our understanding of the principles that govern the in vivo dynamics of liquid-like bodies.

Contrasting patterns of microbial dominance in the Arabidopsis thaliana phyllosphere.

Sphingomonas is one of the most abundant bacterial genera in the phyllosphere of wild Arabidopsis thaliana, but relative to Pseudomonas, the ecology of Sphingomonas and its interaction with plants is poorly described. We analyzed the genomic features of over 400 Sphingomonas isolates collected from local A. thaliana populations, which revealed much higher intergenomic diversity than for the considerably more uniform Pseudomonas isolates found in the same host populations. Variation in Sphingomonas plasmid complements and additional genomic features suggest high adaptability of this genus, and the widespread presence of protein secretion systems hints at frequent biotic interactions. While some of the isolates showed plant-protective phenotypes in lab tests, this was a rare trait. To begin to understand the extent of strain sharing across alternate hosts, we employed amplicon sequencing and a bulk-culturing metagenomics approach on both A. thaliana and neighboring plants. Our data reveal that both Sphingomonas and Pseudomonas thrive on other diverse plant hosts, but that Sphingomonas is a poor competitor in dying or dead leaves.

Predictable and stable epimutations induced during clonal plant propagation with embryonic transcription factor.

Clonal propagation is frequently used in commercial plant breeding and biotechnology programs because it minimizes genetic variation, yet it is not uncommon to observe clonal plants with stable phenotypic changes, a phenomenon known as somaclonal variation. Several studies have linked epigenetic modifications induced during regeneration with this newly acquired phenotypic variation. However, the factors that determine the extent of somaclonal variation and the molecular changes underpinning this process remain poorly understood. To address this gap in our knowledge, we compared clonally propagated Arabidopsis thaliana plants derived from somatic embryogenesis using two different embryonic transcription factors- RWP-RK DOMAIN-CONTAINING 4 (RKD4) or LEAFY COTYLEDON2 (LEC2) and from two epigenetically distinct founder tissues. We found that both the epi(genetic) status of the explant and the regeneration protocol employed play critical roles in shaping the molecular and phenotypic landscape of clonal plants. Phenotypic variation in regenerated plants can be largely explained by the inheritance of tissue-specific DNA methylation imprints, which are associated with specific transcriptional and metabolic changes in sexual progeny of clonal plants. For instance, regenerants were particularly affected by the inheritance of root-specific epigenetic imprints, which were associated with an increased accumulation of salicylic acid in leaves and accelerated plant senescence. Collectively, our data reveal specific pathways underpinning the phenotypic and molecular variation that arise and accumulate in clonal plant populations.

miR472 deficiency enhances Arabidopsis thaliana defense without reducing seed production.

After having co-existed in plant genomes for at least 200 million years, the products of microRNA (miRNA) and Nucleotide-Binding Leucine Rich Repeat protein (NLR) genes formed a regulatory relationship in the common ancestor of modern gymnosperms and angiosperms. From then on, DNA polymorphisms occurring at miRNA target sequences within NLR transcripts must have been compensated by mutations in the corresponding mature miRNA sequence. The potential evolutionary advantage of such regulation remains largely unknown and might be related to two non-exclusive scenarios: miRNA-dependent regulation of NLR levels might prevent defense mis-activation with negative effects on plant growth and reproduction; or reduction of active miRNA levels in response to pathogen derived molecules (PAMPS and silencing suppressors) might rapidly release otherwise silent NLR transcripts for rapid translation and thereby enhance defense. Here, we used Arabidopsis thaliana plants deficient for miR472 function to study the impact of releasing its NLR targets on plant growth and reproduction and on defense against the fungal pathogen Plectospharaella cucumerina. We show that miR472 regulation has a dual role, participating both in the tight regulation of plant defense and growth. MIM472 lines, with reduced active miR472, are more resistant to pathogens and, correlatively, have reduced relative growth compared to wild-type plants although the end of their reproductive phase is delayed, exhibiting higher adult biomass and similar seed yield as the wild-type. Our study highlights how negative consequences of defense activation might be compensated by changes in phenology and that miR472 reduction is an integral part of plant defense responses.