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

Molecular tug-of-war: Plant immune recognition of herbivory

Plant defense responses against insect herbivores are induced through wound-induced signaling and the specific perception of herbivore-associated molecular patterns (HAMPs). In addition, herbivores can deliver effectors that suppress plant immunity. Here we review plant immune recognition of HAMPs and effectors, and argue that these initial molecular interactions upon a plant-herbivore encounter mediate and structure effective resistance. While the number of distinct HAMPs and effectors from both chewing and piercing-sucking herbivores has expanded rapidly with omics-enabled approaches, paired receptors and targets in the host are still not well characterized. Herbivore-derived effectors may also be recognized as HAMPs depending on the host plant species, potentially through the evolution of novel immune receptor functions. We compile examples of HAMPs and effectors where natural variation between species may inform evolutionary patterns and mechanisms of plant-herbivore interactions. Finally, we discuss the combined effects of wounding and HAMP recognition, and review potential signaling hubs, which may integrate both sensing functions. Understanding the precise mechanisms for plant sensing of herbivores will be critical for engineering resistance in agriculture.

The N-terminally truncated helper NLR NRG1C antagonizes immunity mediated by its full-length neighbors NRG1A and NRG1B

Both plants and animals utilize nucleotide-binding leucine-rich repeat immune receptors (NLRs) to perceive the presence of pathogen-derived molecules and induce immune responses. NLR genes are far more abundant and diverse in vascular plants than in animals. Truncated NLRs, which lack one or more of the canonical domains, are also commonly encoded in plant genomes. However, little is known about their functions, especially the N-terminally truncated ones. Here, we show that the Arabidopsis thaliana N-terminally truncated helper NLR (hNLR) gene N REQUIREMENT GENE1 (NRG1C) is highly induced upon pathogen infection and in autoimmune mutants. The immune response and cell death conferred by some Toll/interleukin-1 receptor-type NLRs (TNLs) were compromised in Arabidopsis NRG1C overexpression lines. Detailed genetic analysis revealed that NRG1C antagonizes the immunity mediated by its full-length neighbors NRG1A and NRG1B. Biochemical tests suggested that NRG1C might interfere with the EDS1-SAG101 complex, which functions in immunity signaling together with NRG1A/1B. Interestingly, Brassicaceae NRG1Cs are functionally exchangeable and that the Nicotiana benthamiana N-terminally truncated hNLR NRG2 also antagonizes NRG1 activity. Together, our study uncovers an unexpected negative role of N-terminally truncated hNLRs in immunity in different plant species.

Molecular innovations in plant TIR-based immunity signaling

A protein domain (Toll and Interleukin-1 receptor [TIR]-like) with homology to animal TIRs mediates immune signaling in prokaryotes and eukaryotes. Here, we present an overview of TIR evolution and the molecular versatility of TIR domains in different protein architectures for host protection against microbial attack. Plant TIR-based signaling emerges as being central to the potentiation and effectiveness of host defenses triggered by intracellular and cell-surface immune receptors. Equally relevant for plant fitness are mechanisms that limit potent TIR signaling in healthy tissues but maintain preparedness for infection. We propose that seed plants evolved a specialized protein module to selectively translate TIR enzymatic activities to defense outputs, overlaying a more general function of TIRs.

A tale of many families: calcium channels in plant immunity

Plants launch a concerted immune response to dampen potential infections upon sensing microbial pathogen and insect invasions. The transient and rapid elevation of the cytosolic calcium concentration [Ca2+]cyt is among the essential early cellular responses in plant immunity. The free Ca2+ concentration in the apoplast is far higher than that in the resting cytoplasm. Thus, the precise regulation of calcium channel activities upon infection is the key for an immediate and dynamic Ca2+ influx to trigger downstream signaling. Specific Ca2+ signatures in different branches of the plant immune system vary in timing, amplitude, duration, kinetics, and sources of Ca2+. Recent breakthroughs in the studies of diverse groups of classical calcium channels highlight the instrumental role of Ca2+ homeostasis in plant immunity and cell survival. Additionally, the identification of some immune receptors as noncanonical Ca2+-permeable channels opens a new view of how immune receptors initiate cell death and signaling. This review aims to provide an overview of different Ca2+-conducting channels in plant immunity and highlight their molecular and genetic mode-of-actions in facilitating immune signaling. We also discuss the regulatory mechanisms that control the stability and activity of these channels.

Recombination, selection and the evolution of tandem gene arrays

Multigene families—immunity genes or sensory receptors, for instance—are often subject to diversifying selection. Allelic diversity may be favored not only through balancing or frequency-dependent selection at individual loci but also by associating different alleles in multicopy gene families. Using a combination of analytical calculations and simulations, we explored a population genetic model of epistatic selection and unequal recombination, where a trade-off exists between the benefit of allelic diversity and the cost of copy abundance. Starting from the neutral case, where we showed that gene copy number is Gamma distributed at equilibrium, we derived also the mean and shape of the limiting distribution under selection. Considering a more general model, which includes variable population size and population substructure, we explored by simulations mean fitness and some summary statistics of the copy number distribution. We determined the relative effects of selection, recombination, and demographic parameters in maintaining allelic diversity and shaping the mean fitness of a population. One way to control the variance of copy number is by lowering the rate of unequal recombination. Indeed, when encoding recombination by a rate modifier locus, we observe exactly this prediction. Finally, we analyzed the empirical copy number distribution of 3 genes in human and estimated recombination and selection parameters of our model.

RFPDR: a random forest approach for plant disease resistance protein prediction

Background

Plant innate immunity relies on a broad repertoire of receptor proteins that can detect pathogens and trigger an effective defense response. Bioinformatic tools based on conserved domain and sequence similarity are within the most popular strategies for protein identification and characterization. However, the multi-domain nature, high sequence diversity and complex evolutionary history of disease resistance (DR) proteins make their prediction a real challenge. Here we present RFPDR, which pioneers the application of Random Forest (RF) for Plant DR protein prediction.

Methods

A recently published collection of experimentally validated DR proteins was used as a positive dataset, while 10x10 nested datasets, ranging from 400-4,000 non-DR proteins, were used as negative datasets. A total of 9,631 features were extracted from each protein sequence, and included in a full dimension (FD) RFPDR model. Sequence selection was performed, to generate a reduced-dimension (RD) RFPDR model. Model performances were evaluated using an 80/20 (training/testing) partition, with 10-cross fold validation, and compared to baseline, sequence-based and state-of-the-art strategies. To gain some insights into the underlying biology, the most discriminatory sequence-based features in the RF classifier were identified.

Results and Discussion

RD-RFPDR showed to be sensitive (86.4 ± 4.0%) and specific (96.9 ± 1.5%) for identifying DR proteins, while robust to data imbalance. Its high performance and robustness, added to the fact that RD-RFPDR provides valuable information related to DR proteins underlying properties, make RD-RFPDR an interesting approach for DR protein prediction, complementing the state-of-the-art strategies.

Discovery of Resistance Genes in Rye by Targeted Long-Read Sequencing and Association Genetics

The majority of released rye cultivars are susceptible to leaf rust because of a low level of resistance in the predominant hybrid rye-breeding gene pools Petkus and Carsten. To discover new sources of leaf rust resistance, we phenotyped a diverse panel of inbred lines from the less prevalent Gülzow germplasm using six distinct isolates of Puccinia recondita f. sp. secalis and found that 55 out of 92 lines were resistant to all isolates. By performing a genome-wide association study using 261,406 informative SNP markers, we identified five resistance-associated QTLs on chromosome arms 1RS, 1RL, 2RL, 5RL and 7RS. To identify candidate Puccinia recondita (Pr) resistance genes in these QTLs, we sequenced the rye nucleotide-binding leucine-rich repeat (NLR) intracellular immune receptor complement using a Triticeae NLR bait-library and PacBio^® long-read single-molecule high-fidelity (HiFi) sequencing. Trait-genotype correlations across 10 resistant and 10 susceptible lines identified four candidate NLR-encoding Pr genes. One of these physically co-localized with molecular markers delimiting Pr3 on chromosome arm 1RS and the top-most resistance-associated QTL in the panel.

Genome sequences of three Aegilops species of the section Sitopsis reveal phylogenetic relationships and provide resources for wheat improvement

Aegilops is a close relative of wheat (Triticum spp.), and Aegilops species in the section Sitopsis represent a rich reservoir of genetic diversity for the improvement of wheat. To understand their diversity and advance their utilization, we produced whole-genome assemblies of Aegilops longissima and Aegilops speltoides. Whole-genome comparative analysis, along with the recently sequenced Aegilops sharonensis genome, showed that the Ae. longissima and Ae. sharonensis genomes are highly similar and are most closely related to the wheat D subgenome. By contrast, the Ae. speltoides genome is more closely related to the B subgenome. Haplotype block analysis supported the idea that Ae. speltoides genome is closest to the wheat B subgenome, and highlighted variable and similar genomic regions between the three Aegilops species and wheat. Genome-wide analysis of nucleotide-binding leucine-rich repeat (NLR) genes revealed species-specific and lineage-specific NLR genes and variants, demonstrating the potential of Aegilops genomes for wheat improvement.

New recognition specificity in a plant immune receptor by molecular engineering of its integrated domain

Plant nucleotide-binding and leucine-rich repeat domain proteins (NLRs) are immune sensors that recognize pathogen effectors. Here, we show that molecular engineering of the integrated decoy domain (ID) of an NLR can extend its recognition spectrum to a new effector. We relied for this on detailed knowledge on the recognition of the Magnaporthe oryzae effectors AVR-PikD, AVR-Pia, and AVR1-CO39 by, respectively, the rice NLRs Pikp-1 and RGA5. Both receptors detect their effectors through physical binding to their HMA (Heavy Metal-Associated) IDs. By introducing into RGA5_HMA the AVR-PikD binding residues of Pikp-1_HMA, we create a high-affinity binding surface for this effector. RGA5 variants carrying this engineered binding surface perceive the new ligand, AVR-PikD, and still recognize AVR-Pia and AVR1-CO39 in the model plant N. benthamiana. However, they do not confer extended disease resistance specificity against M. oryzae in transgenic rice plants. Altogether, our study provides a proof of concept for the design of new effector recognition specificities in NLRs through molecular engineering of IDs.

Methoxyacrylate Fungicide Candidate CL-15C Also Functions as a Plant Elicitor in Arabidopsis thaliana and Oryza sativa L

Unlike traditional fungicides targeting fungi, plant elicitors usually lack direct fungicidal activity but improve the plant immune system to resist fungi infection, which has gained increasing attention for better fungi resistance management and environment protection. (E)-methyl-2-(2-((((Z)-(amino-(3,4-dichloroisothiazol-5-yl)methylene)amino)oxy)methyl)phenyl)-2-(methoxyimino)acetate (CL-15C) was found to be a fungicide candidate with a broad spectrum. Here, we studied its immune-inducing ability and mechanism to strengthen the resistance of Arabidopsis thaliana against Pseudomonas syringae pv. tomato DC3000 (Pst DC3000) and Oryza sativa L. against Magnaporthe oryzae. CL-15C promoted a 2.20- and 1.47-fold increase in phenylalanine ammonia-lyase (PAL) activity in A. thaliana and O. sativa, respectively. It also facilitated a 1.89- and 1.32-fold increase in accumulation of salicylic acid (SA) in A. thaliana and O. sativa, respectively. Differential genes were clustered in the SA signaling pathway at 24 h after a CL-15C treatment in A. thaliana. Because PAL is a rate-limiting enzyme in the phenylalanine metabolic pathway, after a CL-15C treatment, a pal1(PAL 1) mutant was more susceptible to Pst DC3000 when compared with the wild type. Bacterial counts in leaves after a CL-15C treatment showed a 1.11-fold reduction in the pal1 mutant and a 1.54-fold reduction in the wild type. The effect of CL-15C on the PAL enzyme activity and SA content was attenuated in the pal1 mutant. Present experimental data implied that the immune-inducing activity of CL-15C was dependent on PAL gene-mediated synthesis of SA.

Chromosome-scale and haplotype-resolved genome assembly of a tetraploid potato cultivar

Potato is the most widely produced tuber crop worldwide. However, reconstructing the four haplotypes of its autotetraploid genome remained an unsolved challenge. Here, we report the 3.1 Gb haplotype-resolved (at 99.6% precision), chromosome-scale assembly of the potato cultivar ‘Otava’ based on high-quality long reads, single-cell sequencing of 717 pollen genomes and Hi-C data. Unexpectedly, ~50% of the genome was identical-by-descent due to recent inbreeding, which was contrasted by highly abundant structural rearrangements involving ~20% of the genome. Among 38,214 genes, only 54% were present in all four haplotypes with an average of 3.2 copies per gene. Taking the leaf transcriptome as an example, 11% of the genes were differently expressed in at least one haplotype, where 25% of them were likely regulated through allele-specific DNA methylation. Our work sheds light on the recent breeding history of potato, the functional organization of its tetraploid genome and has the potential to strengthen the future of genomics-assisted breeding.

Haplotype-resolved genome assembly of the tetraploid potato cultivar ‘Otava’ sheds light on functional organization of the tetraploid genome and provides the potential for genomics-assisted breeding.

Genome-wide identification of the NLR gene family in Haynaldia villosa by SMRT-RenSeq

Background

Nucleotide-binding and leucine-rich repeat (NLR) genes have attracted wide attention due to their crucial role in protecting plants from pathogens. SMRT-RenSeq, combining PacBio sequencing after resistance gene enrichment sequencing (RenSeq), is a powerful method for selectively capturing and sequencing full-length NLRs. Haynaldia villosa, a wild grass species with a proven potential for wheat improvement, confers resistance to multiple diseases. So, genome-wide identification of the NLR gene family in Haynaldia villosa by SMRT-RenSeq can facilitate disease resistance genes exploration.

Results

In this study, SMRT-RenSeq was performed to identify the genome-wide NLR complement of H. villosa. In total, 1320 NLRs were annotated in 1169 contigs, including 772 complete NLRs. All the complete NLRs were phylogenetically analyzed and 11 main clades with special characteristics were derived. NLRs could be captured with high efficiency when aligned with cloned R genes, and cluster expansion in some specific gene loci was observed. The physical location of NLRs to individual chromosomes in H. villosa showed a perfect homoeologous relationship with group 1, 2, 3, 5 and 6 of other Triticeae species, however, NLRs physically located on 4VL were largely in silico predicted to be located on the homoeologous group 7. Fifteen types of integrated domains (IDs) were integrated in 52 NLRs, and Kelch and B3 NLR-IDs were found to have expanded in H. villosa, while DUF948, NAM-associated and PRT_C were detected as unique integrated domains implying the new emergence of NLR-IDs after H. villosa diverged from other species.

Conclusion

SMRT-RenSeq is a powerful tool to identify NLR genes from wild species using the baits of the evolutionary related species with reference sequences. The availability of the NLRs from H. villosa provide a valuable library for R gene mining and transfer of disease resistance into wheat.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08334-w.

NLR diversity and candidate fusiform rust resistance genes in loblolly pine

Resistance to fusiform rust disease in loblolly pine (Pinus taeda) is a classic gene-for-gene system. Early resistance gene mapping in the P. taeda family 10-5 identified RAPD markers for a major fusiform rust resistance gene, Fr1. More recently, single nucleotide polymorphism (SNP) markers associated with resistance were mapped to a full-length gene model in the loblolly pine genome encoding for a nucleotide-binding site leucine-rich repeat (NLR) protein. NLR genes are one of the most abundant gene families in plant genomes and are involved in effector-triggered immunity. Inter- and intraspecies studies of NLR gene diversity and expression have resulted in improved disease resistance. To characterize NLR gene diversity and discover potential resistance genes, we assembled de novo transcriptomes from 92 loblolly genotypes from across the natural range of the species. In these transcriptomes, we identified novel NLR transcripts that are not present in the loblolly pine reference genome and found significant geographic diversity of NLR genes providing evidence of gene family evolution. We designed capture probes for these NLRs to identify and map SNPs that stably cosegregate with resistance to the SC20-21 isolate of Cronartium quercuum f.sp. fusiforme (Cqf) in half-sib progeny of the 10-5 family. We identified 10 SNPs and 2 quantitative trait loci associated with resistance to SC20-21 Cqf. The geographic diversity of NLR genes provides evidence of NLR gene family evolution in loblolly pine. The SNPs associated with rust resistance provide a resource to enhance breeding and deployment of resistant pine seedlings.

Structural basis of NLR activation and innate immune signalling in plants

Animals and plants have NLRs (nucleotide-binding leucine-rich repeat receptors) that recognize the presence of pathogens and initiate innate immune responses. In plants, there are three types of NLRs distinguished by their N-terminal domain: the CC (coiled-coil) domain NLRs, the TIR (Toll/interleukin-1 receptor) domain NLRs and the RPW8 (resistance to powdery mildew 8)-like coiled-coil domain NLRs. CC-NLRs (CNLs) and TIR-NLRs (TNLs) generally act as sensors of effectors secreted by pathogens, while RPW8-NLRs (RNLs) signal downstream of many sensor NLRs and are called helper NLRs. Recent studies have revealed three dimensional structures of a CNL (ZAR1) including its inactive, intermediate and active oligomeric state, as well as TNLs (RPP1 and ROQ1) in their active oligomeric states. Furthermore, accumulating evidence suggests that members of the family of lipase-like EDS1 (enhanced disease susceptibility 1) proteins, which are uniquely found in seed plants, play a key role in providing a link between sensor NLRs and helper NLRs during innate immune responses. Here, we summarize the implications of the plant NLR structures that provide insights into distinct mechanisms of action by the different sensor NLRs and discuss plant NLR-mediated innate immune signalling pathways involving the EDS1 family proteins and RNLs.

Network organization of the plant immune system: from pathogen perception to robust defense induction

The plant immune system has been explored essentially through the study of qualitative resistance, a simple form of immunity, and from a reductionist point of view. The recent identification of genes conferring quantitative disease resistance revealed a large array of functions, suggesting more complex mechanisms. In addition, thanks to the advent of high-throughput analyses and system approaches, our view of the immune system has become more integrative, revealing that plant immunity should rather be seen as a distributed and highly connected molecular network including diverse functions to optimize expression of plant defenses to pathogens. Here, we review the recent progress made to understand the network complexity of regulatory pathways leading to plant immunity, from pathogen perception, through signaling pathways and finally to immune responses. We also analyze the topological organization of these networks and their emergent properties, crucial to predict novel immune functions and test them experimentally. Finally, we report how these networks might be regulated by environmental clues. Although system approaches remain extremely scarce in this area of research, a growing body of evidence indicates that the plant response to combined biotic and abiotic stresses cannot be inferred from responses to individual stresses. A view of possible research avenues in this nascent biology domain is finally proposed.

An Improved Assembly of the Albugo candida Ac2V Genome Reveals the Expansion of the "CCG" Class of Effectors

Albugo candida is an obligate oomycete pathogen that infects many plants in the Brassicaceae family. We resequenced the genome of isolate Ac2V using PacBio long reads and constructed an assembly augmented by Illumina reads. The Ac2VPB genome assembly is 10% larger and more contiguous compared with a previous version. Our annotation of the new assembly, aided by RNA-sequencing information, revealed a 175% expansion (40 to 110) in the CHxC effector class, which we redefined as "CCG" based on motif analysis. This class of effectors consist of arrays of phylogenetically related paralogs residing in gene sparse regions, and shows signatures of positive selection and presence/absence polymorphism. This work provides a resource that allows the dissection of the genomic components underlying A. candida adaptation and, particularly, the role of CCG effectors in virulence and avirulence on different hosts.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

OUP accepted manuscript

Potyviridae, the largest family of known RNA viruses (realm Riboviria), belongs to the picorna-like supergroup and has important agricultural and ecological impacts. Potyvirid genomes are translated into polyproteins, which are in turn hydrolyzed to release mature products. Recent sequencing efforts revealed an unprecedented number of potyvirids with a rich variability in gene content and genomic layouts. Here, we review the heterogeneity of non-core modules that expand the structural and functional diversity of the potyvirid proteomes. We provide a family-wide classification of P1 proteinases into the functional Types A and B, and discuss pretty interesting sweet potato potyviral ORF (PISPO), putative zinc fingers, and alkylation B (AlkB)—non-core modules found within P1 cistrons. The atypical inosine triphosphate pyrophosphatase (ITPase/HAM1), as well as the pseudo tobacco mosaic virus-like coat protein (TMV-like CP) are discussed alongside homologs of unrelated virus taxa. Family-wide abundance of the multitasking helper component proteinase (HC-pro) is revised. Functional connections between non-core modules are highlighted to support host niche adaptation and immune evasion as main drivers of the Potyviridae evolutionary radiation. Potential biotechnological and synthetic biology applications of potyvirid leader proteinases and non-core modules are finally explored.

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.

Arabidopsis bioinformatics: tools and strategies

The sequencing of the Arabidopsis thaliana genome 21 years ago ushered in the genomics era for plant research. Since then, an incredible variety of bioinformatic tools permit easy access to large repositories of genomic, transcriptomic, proteomic, epigenomic and other '-omic' data. In this review, we cover some more recent tools (and highlight the 'classics') for exploring such data in order to help formulate quality, testable hypotheses, often without having to generate new experimental data. We cover tools for examining gene expression and co-expression patterns, undertaking promoter analyses and gene set enrichment analyses, and exploring protein-protein and protein-DNA interactions. We will touch on tools that integrate different data sets at the end of the article.

Recurrent expansions of B30.2-associated immune receptor families in fish

B30.2 domains, also known as PRY/SPRY, are key components of specific subsets of two large families of proteins involved in innate immunity: the tripartite motif proteins (TRIMs) and the Nod-like receptors (NLRs). TRIM proteins are important, often inducible factors of antiviral innate immunity, targeting multiple steps of viral cycles through a variety of mechanisms. NLRs prime and regulate systemic innate defenses, especially against bacteria, and control inflammation. Large TRIM and NLR subsets characterized by the presence of a B30.2 domain have been reported from a few fish species including zebrafish and seem to be strongly prone to gene duplication/expansion. Here, we performed a large-scale survey of these receptors across about 150 fish genomes, focusing on ray-finned fishes. We assessed the number and genomic distribution of domains and domain combinations associated with TRIMs, NLRs, and other genes containing B30.2 domains and looked for gene expansion patterns across fish groups. We then used a model to test the impact of taxonomy, genome size, and environmental variables on the copy numbers of these genes. Our findings reveal novel domain structures, clade-specific gains and losses. They also assist with the timing of the gene expansions, reveal patterns associated with the MHC, and lay the groundwork for further studies delving deeper into the forces that drive the copy number variation of immune genes on a species level.

Discovery of a novel powdery mildew (Blumeria graminis) resistance locus in rye (Secale cereale L.)

Powdery mildew is one of the most destructive diseases in the world, causing substantial grain yield losses and quality reduction in cereal crops. At present 23 powdery mildew resistance genes have been identified in rye, of which the majority are in wheat-rye translocation lines developed for wheat improvement. Here, we investigated the genetics underlying powdery mildew resistance in the Gülzow-type elite hybrid rye (Secale cereale L.) breeding germplasm. In total, 180 inbred breeding lines were genotyped using the state-of-the-art 600 K SNP array and phenotyped for infection type against three distinct field populations of B. graminis f. sp. secalis from Northern Germany (2013 and 2018) and Denmark (2020). We observed a moderate level of powdery mildew resistance in the non-restorer germplasm population, and by performing a genome-wide association study using 261,406 informative SNP markers, we identified a powdery mildew resistance locus, provisionally denoted PmNOS1, on the distal tip of chromosome arm 7RL. Using recent advances in rye genomic resources, we investigated whether nucleotide-binding leucine-rich repeat genes residing in the identified 17 Mbp block associated with PmNOS1 on recent reference genomes resembled known Pm genes.

Functional diversification gave rise to allelic specialization in a rice NLR immune receptor pair

Cooperation between receptors from the nucleotide-binding, leucine-rich repeats (NLR) superfamily is important for intracellular activation of immune responses. NLRs can function in pairs that, upon pathogen recognition, trigger hypersensitive cell death and stop pathogen invasion. Natural selection drives specialization of host immune receptors towards an optimal response, whilst keeping a tight regulation of immunity in the absence of pathogens. However, the molecular basis of co-adaptation and specialization between paired NLRs remains largely unknown. Here, we describe functional specialization in alleles of the rice NLR pair Pik that confers resistance to strains of the blast fungus Magnaporthe oryzae harbouring AVR-Pik effectors. We revealed that matching pairs of allelic Pik NLRs mount effective immune responses, whereas mismatched pairs lead to autoimmune phenotypes, a hallmark of hybrid necrosis in both natural and domesticated plant populations. We further showed that allelic specialization is largely underpinned by a single amino acid polymorphism that determines preferential association between matching pairs of Pik NLRs. These results provide a framework for how functionally linked immune receptors undergo co-adaptation to provide an effective and regulated immune response against pathogens. Understanding the molecular constraints that shape paired NLR evolution has implications beyond plant immunity given that hybrid necrosis can drive reproductive isolation.

A giant NLR gene confers broad-spectrum resistance to Phytophthora sojae in soybean

Phytophthora root and stem rot caused by P. sojae is a destructive soybean soil-borne disease found worldwide. Discovery of genes conferring broad-spectrum resistance to the pathogen is a need to prevent the outbreak of the disease. Here, we show that soybean Rps11 is a 27.7-kb nucleotide-binding site-leucine-rich repeat (NBS-LRR or NLR) gene conferring broad-spectrum resistance to the pathogen. Rps11 is located in a genomic region harboring a cluster of large NLR genes of a single origin in soybean, and is derived from rounds of unequal recombination. Such events result in promoter fusion and LRR expansion that may contribute to the broad resistance spectrum. The NLR gene cluster exhibits drastic structural diversification among phylogenetically representative varieties, including gene copy number variation ranging from five to 23 copies, and absence of allelic copies of Rps11 in any of the non-Rps11-donor varieties examined, exemplifying innovative evolution of NLR genes and NLR gene clusters.

The INCREASE project: Intelligent Collections of food-legume genetic resources for European agrofood systems

Food legumes are crucial for all agriculture-related societal challenges, including climate change mitigation, agrobiodiversity conservation, sustainable agriculture, food security and human health. The transition to plant-based diets, largely based on food legumes, could present major opportunities for adaptation and mitigation, generating significant co-benefits for human health. The characterization, maintenance and exploitation of food-legume genetic resources, to date largely unexploited, form the core development of both sustainable agriculture and a healthy food system. INCREASE will implement, on chickpea (Cicer arietinum), common bean (Phaseolus vulgaris), lentil (Lens culinaris) and lupin (Lupinus albus and L. mutabilis), a new approach to conserve, manage and characterize genetic resources. Intelligent Collections, consisting of nested core collections composed of single-seed descent-purified accessions (i.e., inbred lines), will be developed, exploiting germplasm available both from genebanks and on-farm and subjected to different levels of genotypic and phenotypic characterization. Phenotyping and gene discovery activities will meet, via a participatory approach, the needs of various actors, including breeders, scientists, farmers and agri-food and non-food industries, exploiting also the power of massive metabolomics and transcriptomics and of artificial intelligence and smart tools. Moreover, INCREASE will test, with a citizen science experiment, an innovative system of conservation and use of genetic resources based on a decentralized approach for data management and dynamic conservation. By promoting the use of food legumes, improving their quality, adaptation and yield and boosting the competitiveness of the agriculture and food sector, the INCREASE strategy will have a major impact on economy and society and represents a case study of integrative and participatory approaches towards conservation and exploitation of crop genetic resources.

Two gap-free reference genomes and a global view of the centromere architecture in rice

Rice (Oryza sativa), a major staple throughout the world and a model system for plant genomics and breeding, was the first crop genome sequenced almost two decades ago. However, reference genomes for all higher organisms to date contain gaps and missing sequences. Here, we report the assembly and analysis of gap-free reference genome sequences for two elite O. sativa xian/indica rice varieties, Zhenshan 97 and Minghui 63, which are being used as a model system for studying heterosis and yield. Gap-free reference genomes provide the opportunity for a global view of the structure and function of centromeres. We show that all rice centromeric regions share conserved centromere-specific satellite motifs with different copy numbers and structures. In addition, the similarity of CentO repeats in the same chromosome is higher than across chromosomes, supporting a model of local expansion and homogenization. Both genomes have over 395 non-TE genes located in centromere regions, of which ∼41% are actively transcribed. Two large structural variants at the end of chromosome 11 affect the copy number of resistance genes between the two genomes. The availability of the two gap-free genomes lays a solid foundation for further understanding genome structure and function in plants and breeding climate-resilient varieties.

A new comprehensive annotation of leucine-rich repeat-containing receptors in rice

Oryza sativa (rice) plays an essential food security role for more than half of the world's population. Obtaining crops with high levels of disease resistance is a major challenge for breeders, especially today, given the urgent need for agriculture to be more sustainable. Plant resistance genes are mainly encoded by three large leucine-rich repeat (LRR)-containing receptor (LRR-CR) families: the LRR-receptor-like kinase (LRR-RLK), LRR-receptor-like protein (LRR-RLP) and nucleotide-binding LRR receptor (NLR). Using lrrprofiler, a pipeline that we developed to annotate and classify these proteins, we compared three publicly available annotations of the rice Nipponbare reference genome. The extended discrepancies that we observed for LRR-CR gene models led us to perform an in-depth manual curation of their annotations while paying special attention to nonsense mutations. We then transferred this manually curated annotation to Kitaake, a cultivar that is closely related to Nipponbare, using an optimized strategy. Here, we discuss the breakthrough achieved by manual curation when comparing genomes and, in addition to 'functional' and 'structural' annotations, we propose that the community adopts this approach, which we call 'comprehensive' annotation. The resulting data are crucial for further studies on the natural variability and evolution of LRR-CR genes in order to promote their use in breeding future resilient varieties.

Evolutionary trade-offs at the Arabidopsis WRR4A resistance locus underpin alternate Albugo candida race recognition specificities

The oomycete Albugo candida causes white rust of Brassicaceae, including vegetable and oilseed crops, and wild relatives such as Arabidopsis thaliana. Novel White Rust Resistance (WRR) genes from Arabidopsis enable new insights into plant/parasite co-evolution. WRR4A from Arabidopsis accession Columbia (Col-0) provides resistance to many but not all white rust races, and encodes a nucleotide-binding, leucine-rich repeat immune receptor. Col-0 WRR4A resistance is broken by AcEx1, an isolate of A. candida. We identified an allele of WRR4A in Arabidopsis accession Øystese-0 (Oy-0) and other accessions that confers full resistance to AcEx1. WRR4AOy-0 carries a C-terminal extension required for recognition of AcEx1, but reduces recognition of several effectors recognized by the WRR4ACol-0 allele. WRR4AOy-0 confers full resistance to AcEx1 when expressed in the oilseed crop Camelina sativa.

Regulation of Cell Death and Signaling by Pore-Forming Resistosomes

Nucleotide-binding leucine-rich repeat receptors (NLRs) are the largest class of immune receptors in plants. They play a key role in the plant surveillance system by monitoring pathogen effectors that are delivered into the plant cell. Recent structural biology and biochemical analyses have uncovered how NLRs are activated to form oligomeric resistosomes upon the recognition of pathogen effectors. In the resistosome, the signaling domain of the NLR is brought to the center of a ringed structure to initiate immune signaling and regulated cell death (RCD). The N terminus of the coiled-coil (CC) domain of the NLR protein HOPZ-ACTIVATED RESISTANCE 1 likely forms a pore in the plasma membrane to trigger RCD in a way analogous to animal pore-forming proteins that trigger necroptosis or pyroptosis. NLRs that carry TOLL-INTERLEUKIN1-RECEPTOR as a signaling domain may also employ pore-forming resistosomes for RCD execution. In addition, increasing evidence supports intimate connections between NLRs and surface receptors in immune signaling. These new findings are rapidly advancing our understanding of the plant immune system.

One Hundred Years of Hybrid Necrosis: Hybrid Autoimmunity as a Window into the Mechanisms and Evolution of Plant-Pathogen Interactions

Hybrid necrosis in plants refers to a genetic autoimmunity syndrome in the progeny of interspecific or intraspecific crosses. Although the phenomenon was first documented in 1920, it has been unequivocally linked to autoimmunity only recently, with the discovery of the underlying genetic and biochemical mechanisms. The most common causal loci encode immune receptors, which are known to differ within and between species. One mechanism can be explained by the guard hypothesis, in which a guard protein, often a nucleotide-binding site-leucine-rich repeat protein, is activated by interaction with a plant protein that mimics standard guardees modified by pathogen effector proteins. Another surprising mechanism is the formation of inappropriately active immune receptor complexes. In this review, we summarize our current knowledge of hybrid necrosis and discuss how its study is not only informing the understanding of immune gene evolution but also revealing new aspects of plant immune signaling.

Genomic-Assisted Marker Development Suitable for CsCvy-1 Selection in Cucumber Breeding

Cucumber is a widely grown vegetable crop plant and a host to many different plant pathogens. Cucumber vein yellowing virus (CVYV) causes economic losses on cucumber crops in Mediterranean countries and in some part of India such as West Bengal and in African countries such as Sudan. CVYV is an RNA potyvirus transmitted mechanically and by whitefly (Bemisia tabaci) in a semipersistent manner. Control of this virus is heavily dependent on the management of the insect vector and breeding virus-resistant lines. DNA markers have been used widely in conventional plant breeding programs via marker-assisted selection (MAS). However, very few resistance sources against CVYV in cucumber exist, and also the lack of tightly linked molecular markers to these sources restricts the rapid generation of resistant lines. In this work, we used genomics coupled with the bulked segregant analysis method and generated the MAS-friendly Kompetitive allele specific PCR (KASP) markers suitable for CsCvy-1 selection in cucumber breeding using a segregating F2 mapping population and commercial plant lines. Variant analysis was performed to generate single-nucleotide polymorphism (SNP)-based markers for mapping the population and genotyping the commercial lines. We fine-mapped the region by generating new markers down to 101 kb with eight genes. We provided SNP data for this interval, which could be useful for breeding programs and cloning the candidate genes.

Dynamic Diversity of NLR Genes in Triticum and Mining of Promising NLR Alleles for Disease Resistance

Bread wheat is an essential crop with the second-highest global production after maize. Currently, wheat diseases are a serious threat to wheat production. Therefore, efficient breeding for disease resistance is extremely urgent in modern wheat. Here, we identified 2012 NLR genes from hexaploid wheat, and Ks values of paired syntenic NLRs showed a significant peak at 3.1–6.3 MYA, which exactly coincided with the first hybridization event between A and B genome lineages at ~5.5 MYA. We provided a landscape of dynamic diversity of NLRs from Triticum and Aegilops and found that NLR genes have higher diversity in wild progenitors and relatives. Further, most NLRs had opposite diversity patterns between genic and 2 Kb-promoter regions, which might respectively link sub/neofunctionalization and loss of duplicated NLR genes. Additionally, we identified an alien introgression of chromosome 4A in tetraploid emmer wheat, which was similar to that in hexaploid wheat. Transcriptome data from four experiments of wheat disease resistance helped to profile the expression pattern of NLR genes and identified promising NLRs involved in broad-spectrum disease resistance. Our study provided insights into the diversity evolution of NLR genes and identified beneficial NLRs to deploy into modern wheat in future wheat disease-resistance breeding.

Polymorphisms in plants to restrict losses to pathogens: From gene family expansions to complex network evolution

Genetic polymorphisms are the basis of the natural diversity seen in all life on earth, also in plant-pathogen interactions. Initially, studies on plant-pathogen interaction focused on reporting phenotypic variation in resistance properties and on the identification of underlying major genes. Nowadays, the field of plant-pathogen interactions is moving from focusing on families of single dominant genes involved in gene-for-gene interactions to an understanding of the plant immune system in the context of a much more complex signaling network and quantitative resistance. Simultaneously, studies on pathosystems from the wild and genome analyses advanced, revealing tremendous variation in natural plant populations. It is now imperative to place studies on genetic diversity and evolution of plant-pathogen interactions in the appropriate molecular biological, as well as evolutionary, context.

The ETS-ETI cycle: evolutionary processes and metapopulation dynamics driving the diversification of pathogen effectors and host immune factors

The natural diversity of pathogen effectors and host immune components represents a snapshot of the underlying evolutionary processes driving the host-pathogen arms race. In plants, this arms race is manifested by an ongoing cycle of disease and resistance driven by pathogenic effectors that promote disease (effector-triggered susceptibility; ETS) and plant resistance proteins that recognize effector activity to trigger immunity (effector-triggered immunity; ETI). Here we discuss how this ongoing ETS-ETI cycle has shaped the natural diversity of both plant resistance proteins and pathogen effectors. We focus on the evolutionary forces that drive the diversification of the molecules that determine the outcome of plant-pathogen interactions and introduce the concept of metapopulation dynamics (i.e., the introduction of genetic variation from conspecific organisms in different populations) as an alternative mechanism that can introduce and maintain diversity in both host and pathogen populations.

A playbook for developing disease-resistant crops through immune receptor identification and transfer

Plants are resistant to most pathogens because of an immune system that perceives invading microbes and activates defense. A large repertoire of innate immune receptors mediates specific direct or indirect recognition of pathogen-derived molecules. Disease is often a consequence of insufficient immune surveillance, and the transfer of immune receptor genes from resistant plants to susceptible crop varieties is an effective strategy for combating disease outbreaks. We discuss approaches for identifying intracellular and cell surface immune receptors, with particular focus on recently developed and emerging methodologies. We also review considerations for the transfer of immune receptor genes into crop species, including additional host factors that may be required for immune receptor function. Together, these concepts lay out a broadly applicable playbook for developing crop varieties with durable disease resistance.

Two NLR immune receptors acquired high-affinity binding to a fungal effector through convergent evolution of their integrated domain

A subset of plant NLR immune receptors carry unconventional integrated domains in addition to their canonical domain architecture. One example is rice Pik-1 that comprises an integrated heavy metal-associated (HMA) domain. Here, we reconstructed the evolutionary history of Pik-1 and its NLR partner, Pik-2, and tested hypotheses about adaptive evolution of the HMA domain. Phylogenetic analyses revealed that the HMA domain integrated into Pik-1 before Oryzinae speciation over 15 million years ago and has been under diversifying selection. Ancestral sequence reconstruction coupled with functional studies showed that two Pik-1 allelic variants independently evolved from a weakly binding ancestral state to high-affinity binding of the blast fungus effector AVR-PikD. We conclude that for most of its evolutionary history the Pik-1 HMA domain did not sense AVR-PikD, and that different Pik-1 receptors have recently evolved through distinct biochemical paths to produce similar phenotypic outcomes. These findings highlight the dynamic nature of the evolutionary mechanisms underpinning NLR adaptation to plant pathogens.

Two NLR immune receptors acquired high-affinity binding to a fungal effector through convergent evolution of their integrated domain

A subset of plant NLR immune receptors carry unconventional integrated domains in addition to their canonical domain architecture. One example is rice Pik-1 that comprises an integrated heavy metal-associated (HMA) domain. Here, we reconstructed the evolutionary history of Pik-1 and its NLR partner, Pik-2, and tested hypotheses about adaptive evolution of the HMA domain. Phylogenetic analyses revealed that the HMA domain integrated into Pik-1 before Oryzinae speciation over 15 million years ago and has been under diversifying selection. Ancestral sequence reconstruction coupled with functional studies showed that two Pik-1 allelic variants independently evolved from a weakly binding ancestral state to high-affinity binding of the blast fungus effector AVR-PikD. We conclude that for most of its evolutionary history the Pik-1 HMA domain did not sense AVR-PikD, and that different Pik-1 receptors have recently evolved through distinct biochemical paths to produce similar phenotypic outcomes. These findings highlight the dynamic nature of the evolutionary mechanisms underpinning NLR adaptation to plant pathogens.

Multi-omics approach highlights differences between RLP classes in Arabidopsis thaliana

BACKGROUND

The Leucine rich-repeat (LRR) receptor-like protein (RLP) family is a complex gene family with 57 members in Arabidopsis thaliana. Some members of the RLP family are known to be involved in basal developmental processes, whereas others are involved in defence responses. However, functional data is currently only available for a small subset of RLPs, leaving the remaining ones classified as RLPs of unknown function.

RESULTS

Using publicly available datasets, we annotated RLPs of unknown function as either likely defence-related or likely fulfilling a more basal function in plants. Then, using these categories, we can identify important characteristics that differ between the RLP subclasses. We found that the two classes differ in abundance on both transcriptome and proteome level, physical clustering in the genome and putative interaction partners. However, the classes do not differ in the genetic di versity of their individual members in accessible pan-genome data.

CONCLUSIONS

Our work has several implications for work related to functional studies on RLPs as well as for the understanding of RLP gene family evolution. Using our annotations, we can make suggestions on which RLPs can be identified as potential immune receptors using genetics tools and thereby complement disease studies. The lack of differences in nucleotide diversity between the two RLP subclasses further suggests that non-synonymous diversity of gene sequences alone cannot distinguish defence from developmental genes. By contrast, differences in transcript and protein abundance or clustering at genomic loci might also allow for functional annotations and characterisation in other plant species.

Genotype-Specific Expression and NLR Repertoire Contribute to Phenotypic Resistance Diversity in Plantago lanceolata

High levels of phenotypic variation in resistance appears to be nearly ubiquitous across natural host populations. Molecular processes contributing to this variation in nature are still poorly known, although theory predicts resistance to evolve at specific loci driven by pathogen-imposed selection. Nucleotide-binding leucine-rich repeat (NLR) genes play an important role in pathogen recognition, downstream defense responses and defense signaling. Identifying the natural variation in NLRs has the potential to increase our understanding of how NLR diversity is generated and maintained, and how to manage disease resistance. Here, we sequenced the transcriptomes of five different Plantago lanceolata genotypes when inoculated by the same strain of obligate fungal pathogen Podosphaera plantaginis. A de novo transcriptome assembly of RNA-sequencing data yielded 24,332 gene models with N50 value of 1,329 base pairs and gene space completeness of 66.5%. The gene expression data showed highly varying responses where each plant genotype demonstrated a unique expression profile in response to the pathogen, regardless of the resistance phenotype. Analysis on the conserved NB-ARC domain demonstrated a diverse NLR repertoire in P. lanceolata consistent with the high phenotypic resistance diversity in this species. We find evidence of selection generating diversity at some of the NLR loci. Jointly, our results demonstrate that phenotypic resistance diversity results from a crosstalk between different defense mechanisms. In conclusion, characterizing the architecture of resistance in natural host populations may shed unprecedented light on the potential of evolution to generate variation.

Creation and judicious application of a wheat resistance gene atlas

Disease-resistance (R) gene cloning in wheat (Triticum aestivum) has been accelerated by the recent surge of genomic resources, facilitated by advances in sequencing technologies and bioinformatics. However, with the challenges of population growth and climate change, it is vital not only to clone and functionally characterize a few handfuls of R genes, but also to do so at a scale that would facilitate the breeding and deployment of crops that can recognize the wide range of pathogen effectors that threaten agroecosystems. Pathogen populations are continually changing, and breeders must have tools and resources available to rapidly respond to those changes if we are to safeguard our daily bread. To meet this challenge, we propose the creation of a wheat R-gene atlas by an international community of researchers and breeders. The atlas would consist of an online directory from which sources of resistance could be identified and deployed to achieve more durable resistance to the major wheat pathogens, such as wheat rusts, blotch diseases, powdery mildew, and wheat blast. We present a costed proposal detailing how the interacting molecular components governing disease resistance could be captured from both the host and the pathogen through biparental mapping, mutational genomics, and whole-genome association genetics. We explore options for the configuration and genotyping of diversity panels of hexaploid and tetraploid wheat, as well as their wild relatives and major pathogens, and discuss how the atlas could inform a dynamic, durable approach to R-gene deployment. Set against the current magnitude of wheat yield losses worldwide, recently estimated at 21%, this endeavor presents one route for bringing R genes from the lab to the field at a considerable speed and quantity.

Building pan-genome infrastructures for crop plants and their use in association genetics

Pan-genomic studies aim at representing the entire sequence diversity within a species to provide useful resources for evolutionary studies, functional genomics and breeding of cultivated plants. Cost reductions in high-throughput sequencing and advances in sequence assembly algorithms have made it possible to create multiple reference genomes along with a catalogue of all forms of genetic variations in plant species with large and complex or polyploid genomes. In this review, we summarize the current approaches to building pan-genomes as an in silico representation of plant sequence diversity and outline relevant methods for their effective utilization in linking structural with phenotypic variation. We propose as future research avenues (i) transcriptomic and epigenomic studies across multiple reference genomes and (ii) the development of user-friendly and feature-rich pan-genome browsers.

Novel genomic approaches to study antagonistic coevolution between hosts and parasites

Host-parasite coevolution is ubiquitous, shaping genetic and phenotypic diversity and the evolutionary trajectory of interacting species. With the advances of high throughput sequencing technologies applicable to model and non-model organisms alike, it is now feasible to study in greater detail (a) the genetic underpinnings of coevolution, (b) the speed and type of dynamics at coevolving loci, and (c) the genomic consequences of coevolution. This review focuses on three recently developed approaches that leverage information from host and parasite full genome data simultaneously to pinpoint coevolving loci and draw inference on the coevolutionary history. First, co-genome-wide association study (co-GWAS) methods allow pinpointing the loci underlying host-parasite interactions. These methods focus on detecting associations between genetic variants and the outcome of experimental infection tests or on correlations between genomes of naturally infected hosts and their infecting parasites. Second, extensions to population genomics methods can detect genes under coevolution and infer the coevolutionary history, such as fitness costs. Third, correlations between host and parasite population size in time are indicative of coevolution, and polymorphism levels across independent spatially distributed populations of hosts and parasites can reveal coevolutionary loci and infer coevolutionary history. We describe the principles of these three approaches and discuss their advantages and limitations based on coevolutionary theory. We present recommendations for their application to various host (prokaryotes, fungi, plants, and animals) and parasite (viruses, bacteria, fungi, and macroparasites) species. We conclude by pointing out methodological and theoretical gaps to be filled to extract maximum information from full genome data and thereby to shed light on the molecular underpinnings of coevolution.

Machine-learning predicts genomic determinants of meiosis-driven structural variation in a eukaryotic pathogen

Species harbor extensive structural variation underpinning recent adaptive evolution. However, the causality between genomic features and the induction of new rearrangements is poorly established. Here, we analyze a global set of telomere-to-telomere genome assemblies of a fungal pathogen of wheat to establish a nucleotide-level map of structural variation. We show that the recent emergence of pesticide resistance has been disproportionally driven by rearrangements. We use machine learning to train a model on structural variation events based on 30 chromosomal sequence features. We show that base composition and gene density are the major determinants of structural variation. Retrotransposons explain most inversion, indel and duplication events. We apply our model to Arabidopsis thaliana and show that our approach extends to more complex genomes. Finally, we analyze complete genomes of haploid offspring in a four-generation pedigree. Meiotic crossover locations are enriched for new rearrangements consistent with crossovers being mutational hotspots. The model trained on species-wide structural variation accurately predicts the position of >74% of newly generated variants along the pedigree. The predictive power highlights causality between specific sequence features and the induction of chromosomal rearrangements. Our work demonstrates that training sequence-derived models can accurately identify regions of intrinsic DNA instability in eukaryotic genomes.

Comparative analysis of de novo genomes reveals dynamic intra-species divergence of NLRs in pepper

BACKGROUND

Peppers (Capsicum annuum L.) containing distinct capsaicinoids are the most widely cultivated spices in the world. However, extreme genomic diversity among species represents an obstacle to breeding pepper.

RESULTS

Here, we report de novo genome assemblies of Capsicum annuum 'Early Calwonder (non-pungent, ECW)' and 'Small Fruit (pungent, SF)' along with their annotations. In total, we assembled 2.9 Gb of ECW and SF genome sequences, representing over 91% of the estimated genome sizes. Structural and functional annotation of the two pepper genomes generated about 35,000 protein-coding genes each, of which 93% were assigned putative functions. Comparison between newly and publicly available pepper gene annotations revealed both shared and specific gene content. In addition, a comprehensive analysis of nucleotide-binding and leucine-rich repeat (NLR) genes through whole-genome alignment identified five significant regions of NLR copy number variation (CNV). Detailed comparisons of those regions revealed that these CNVs were generated by intra-specific genomic variations that accelerated diversification of NLRs among peppers.

CONCLUSIONS

Our analyses unveil an evolutionary mechanism responsible for generating CNVs of NLRs among pepper accessions, and provide novel genomic resources for functional genomics and molecular breeding of disease resistance in Capsicum species.

Plant NLR diversity: the known unknowns of pan-NLRomes

Plants and pathogens constantly adapt to each other. As a consequence, many members of the plant immune system, and especially the intracellular nucleotide-binding site leucine-rich repeat receptors, also known as NOD-like receptors (NLRs), are highly diversified, both among family members in the same genome, and between individuals in the same species. While this diversity has long been appreciated, its true extent has remained unknown. With pan-genome and pan-NLRome studies becoming more and more comprehensive, our knowledge of NLR sequence diversity is growing rapidly, and pan-NLRomes provide powerful platforms for assigning function to NLRs. These efforts are an important step toward the goal of comprehensively predicting from sequence alone whether an NLR provides disease resistance, and if so, to which pathogens.

Analysis of intraspecies diversity reveals a subset of highly variable plant immune receptors and predicts their binding sites

The evolution of recognition specificities by the immune system depends on the generation of receptor diversity and on connecting the binding of new antigens with the initiation of downstream signaling. In plant immunity, the innate Nucleotide-Binding Leucine-Rich Repeat (NLR) receptor family enables antigen binding and immune signaling. In this study, we surveyed the NLR complements of 62 ecotypes of Arabidopsis thaliana and 54 lines of Brachypodium distachyon and identified a limited number of NLR subfamilies that show high allelic diversity. We show that the predicted specificity-determining residues cluster on the surfaces of Leucine-Rich Repeat domains, but the locations of the clusters vary among NLR subfamilies. By comparing NLR phylogeny, allelic diversity, and known functions of the Arabidopsis NLRs, we formulate a hypothesis for the emergence of direct and indirect pathogen-sensing receptors and of the autoimmune NLRs. These findings reveal the recurring patterns of evolution of innate immunity and can inform NLR engineering efforts.

Genetic and environmental modulation of transposition shapes the evolutionary potential of Arabidopsis thaliana

BACKGROUND

How species can adapt to abrupt environmental changes, particularly in the absence of standing genetic variation, is poorly understood and a pressing question in the face of ongoing climate change. Here we leverage publicly available multi-omic and bio-climatic data for more than 1000 wild Arabidopsis thaliana accessions to determine the rate of transposable element (TE) mobilization and its potential to create adaptive variation in natural settings.

RESULTS

We demonstrate that TE insertions arise at almost the same rate as base substitutions. Mobilization activity of individual TE families varies greatly between accessions, in association with genetic and environmental factors as well as through complex gene-environment interactions. Although the distribution of TE insertions across the genome is ultimately shaped by purifying selection, reflecting their typically strong deleterious effects when located near or within genes, numerous recent TE-containing alleles show signatures of positive selection. Moreover, high rates of transposition appear positively selected at the edge of the species' ecological niche. Based on these findings, we predict through mathematical modeling higher transposition activity in Mediterranean regions within the next decades in response to global warming, which in turn should accelerate the creation of large-effect alleles.

CONCLUSIONS

Our study reveals that TE mobilization is a major generator of genetic variation in A. thaliana that is finely modulated by genetic and environmental factors. These findings and modeling indicate that TEs may be essential genomic players in the demise or rescue of native populations in times of climate crises.

Disentangling cause and consequence: genetic dissection of the DANGEROUS MIX2 risk locus, and activation of the DM2h NLR in autoimmunity

Nucleotide-binding domain-leucine-rich repeat-type immune receptors (NLRs) protect plants against pathogenic microbes through intracellular detection of effector proteins. However, this comes at a cost, as NLRs can also induce detrimental autoimmunity in genetic interactions with foreign alleles. This may occur when independently evolved genomes are combined in inter- or intraspecific crosses, or when foreign alleles are introduced by mutagenesis or transgenesis. Most autoimmunity-inducing NLRs are encoded within highly variable NLR gene clusters with no known immune functions, which were termed autoimmune risk loci. Whether risk NLRs differ from sensor NLRs operating in natural pathogen resistance and how risk NLRs are activated in autoimmunity is unknown. Here, we analyzed the DANGEROUS MIX2 risk locus, a major autoimmunity hotspot in Arabidopsis thaliana. By gene editing and heterologous expression, we show that a single gene, DM2h, is necessary and sufficient for autoimmune induction in three independent cases of autoimmunity in accession Landsberg erecta. We focus on autoimmunity provoked by an EDS1-yellow fluorescent protein (YFP)NLS fusion protein to characterize DM2h functionally and determine features of EDS1-YFPNLS activating the immune receptor. Our data suggest that risk NLRs function in a manner reminiscent of sensor NLRs, while autoimmunity-inducing properties of EDS1-YFPNLS in this context are unrelated to the protein's functions as an immune regulator. We propose that autoimmunity, at least in some cases, may be caused by spurious, stochastic interactions of foreign alleles with coincidentally matching risk NLRs.

Insights into amphicarpy from the compact genome of the legume Amphicarpaea edgeworthii

Amphicarpy (seed heteromorphy) is a unique and fascinating reproductive strategy wherein a single plant produces both aerial and subterranean fruits. This strategy is believed to be an adaptation to life under stressful or uncertain environments. Here, we sequenced and de novo assembled a chromosome-level genome assembly of the legume Amphicarpaea edgeworthii Benth. The 299-Mb A. edgeworthii genome encodes 27 899 protein-coding genes and is the most compact sequenced legume genome reported until date. Its reduced genome size may be attributed to the reduced long-terminal repeat retrotransposon content, which stems from the unequal homologous recombination. Gene families related to immunity and stress resistance have been contracted in A. edgeworthii, which is consistent with the notion that the amphicarpic reproductive strategy may be a complementary mechanism for its weak environmental-adaptation ability. We demonstrated the 'ABCE' model for the differentiation of chasmogamous and cleistogamous flowers. In addition, the characteristics of aerial and subterranean seeds in hard-seededness were explored. Thus, we suggest that the A. edgeworthii genome, which is the first of an amphicarpic plant, offers significant insights into its unusual reproductive strategy that is a key resource towards comprehending the evolution of angiosperms.