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

NRC Immune receptor networks show diversified hierarchical genetic architecture across plant lineages

Plants' complex immune systems include nucleotide-binding domain and leucine-rich repeat-containing (NLR) proteins, which help recognize invading pathogens. In solanaceous plants, the NRC (NLR required for cell death) family includes helper NLRs that form a complex genetic network with multiple sensor NLRs to provide resistance against pathogens. However, the evolution and function of NRC networks outside solanaceous plants are currently unclear. Here, we conducted phylogenomic and macroevolutionary analyses comparing NLRs identified from different asterid lineages and found that NRC networks expanded significantly in most lamiids but not in Ericales and campanulids. Using transient expression assays in Nicotiana benthamiana, we showed that NRC networks are simple in Ericales and campanulids, but have high complexity in lamiids. Phylogenetic analyses grouped the NRC helper NLRs into three NRC0 subclades that are conserved, and several family-specific NRC subclades of lamiids that show signatures of diversifying selection. Functional analyses revealed that members of the NRC0 subclades are partially interchangeable, whereas family-specific NRC members in lamiids lack interchangeability. Our findings highlight the distinctive evolutionary patterns of the NRC networks in asterids and provide potential insights into transferring disease resistance across plant lineages.

The NRC0 gene cluster of sensor and helper NLR immune receptors is functionally conserved across asterid plants

Nucleotide-binding domain and leucine-rich repeat-containing receptor (NLR) proteins can form complex receptor networks to confer innate immunity. An NLR-REQUIRED FOR CELL DEATH (NRC) is a phylogenetically related node that functions downstream of a massively expanded network of disease resistance proteins that protect against multiple plant pathogens. In this study, we used phylogenomic methods to reconstruct the macroevolution of the NRC family. One of the NRCs, termed NRC0, is the only family member shared across asterid plants, leading us to investigate its evolutionary history and genetic organization. In several asterid species, NRC0 is genetically clustered with other NLRs that are phylogenetically related to NRC-dependent disease resistance genes. This prompted us to hypothesize that the ancestral state of the NRC network is an NLR helper-sensor gene cluster that was present early during asterid evolution. We provide support for this hypothesis by demonstrating that NRC0 is essential for the hypersensitive cell death that is induced by its genetically linked sensor NLR partners in 4 divergent asterid species: tomato (Solanum lycopersicum), wild sweet potato (Ipomoea trifida), coffee (Coffea canephora), and carrot (Daucus carota). In addition, activation of a sensor NLR leads to higher-order complex formation of its genetically linked NRC0, similar to other NRCs. Our findings map out contrasting evolutionary dynamics in the macroevolution of the NRC network over the last 125 million years, from a functionally conserved NLR gene cluster to a massive genetically dispersed network.

Alternative splicing of a potato disease resistance gene maintains homeostasis between growth and immunity

Plants possess a robust and sophisticated innate immune system against pathogens and must balance growth with rapid pathogen detection and defense. The intracellular receptors with nucleotide-binding leucine-rich repeat (NLR) motifs recognize pathogen-derived effector proteins and thereby trigger the immune response. The expression of genes encoding NLR receptors is precisely controlled in multifaceted ways. The alternative splicing (AS) of introns in response to infection is recurrently observed but poorly understood. Here we report that the potato (Solanum tuberosum) NLR gene RB undergoes AS of its intron, resulting in 2 transcriptional isoforms, which coordinately regulate plant immunity and growth homeostasis. During normal growth, RB predominantly exists as an intron-retained isoform RB_IR, encoding a truncated protein containing only the N-terminus of the NLR. Upon late blight infection, the pathogen induces intron splicing of RB, increasing the abundance of RB_CDS, which encodes a full-length and active R protein. By deploying the RB splicing isoforms fused with a luciferase reporter system, we identified IPI-O1 (also known as Avrblb1), the RB cognate effector, as a facilitator of RB AS. IPI-O1 directly interacts with potato splicing factor StCWC15, resulting in altered localization of StCWC15 from the nucleoplasm to the nucleolus and nuclear speckles. Mutations in IPI-O1 that eliminate StCWC15 binding also disrupt StCWC15 re-localization and RB intron splicing. Thus, our study reveals that StCWC15 serves as a surveillance facilitator that senses the pathogen-secreted effector and regulates the trade-off between RB-mediated plant immunity and growth, expanding our understanding of molecular plant-microbe interactions.

Molecular dissection of the pseudokinase ZED1 expands effector recognition to the tomato immune receptor ZAR1

The highly conserved angiosperm immune receptor HOPZ-ACTIVATED RESISTANCE 1 (ZAR1) is a bacterial pathogen recognition hub that mediates resistance by guarding host kinases for modification by pathogen effectors. The pseudokinase HOPZ-ETI DEFICIENT 1 (ZED1) is the only known ZAR1-guarded protein that interacts directly with a pathogen effector, HopZ1a, from the bacterial pathogen Pseudomonas syringae, making it a promising system for rational design of effector recognition for plant immunity. Here, we conducted an in-depth molecular analysis of ZED1. We generated a library of 164 random ZED1 mutants and identified 50 mutants that could not recognize the effector HopZ1a when transiently expressed in Nicotiana benthamiana. Based on our random mutants, we generated a library of 27 point mutants and found evidence of minor functional divergence between Arabidopsis (Arabidopsis thaliana) and N. benthamiana ZAR1 orthologs. We leveraged our point mutant library to identify regions in ZED1 critical for ZAR1 and HopZ1a interactions and identified two likely ZED1-HopZ1a binding conformations. We explored ZED1 nucleotide and cation binding activity and showed that ZED1 is a catalytically dead pseudokinase, functioning solely as an allosteric regulator upon effector recognition. We used our library of ZED1 point mutants to identify the ZED1 activation loop regions as the most likely cause of interspecies ZAR1-ZED1 incompatibility. Finally, we identified a mutation that abolished ZAR1-ZED1 interspecies incompatibility while retaining the ability to mediate HopZ1a recognition, which enabled recognition of HopZ1a through tomato (Solanum lycopersicum) ZAR1. This provides an example of expanded effector recognition through a ZAR1 ortholog from a non-model species.

Redox signaling and oxidative stress in systemic acquired resistance

Plants fully depend on their immune systems to defend against pathogens. Upon pathogen attack, plants not only activate immune responses at the infection site but also trigger a defense mechanism known as systemic acquired resistance (SAR) in distal systemic tissues to prevent subsequent infections by a broad-spectrum of pathogens. SAR is induced by mobile signals produced at the infection site. Accumulating evidence suggests that reactive oxygen species (ROS) play a central role in SAR signaling. ROS burst at the infection site is one of the earliest cellular responses following pathogen infection and can spread to systemic tissues through membrane-associated NADPH oxidase-dependent relay production of ROS. It is well known that ROS ignite redox signaling and, when in excess, cause oxidative stress, damaging cellular components. In this review, we summarize current knowledge on redox regulation of several SAR signaling components. We discuss the ROS amplification loop in systemic tissues involving multiple SAR mobile signals. Moreover, we highlight the essential role of oxidative stress in generating SAR signals including azelaic acid and extracellular NAD(P) [eNAD(P)]. Finally, we propose that eNAD(P) is a damage-associated molecular pattern serving as a converging point of SAR mobile signals in systemic tissues.

Extensive immune receptor repertoire diversity in disease-resistant rice landraces

Plants have powerful defense mechanisms and extensive immune receptor repertoires, yet crop monocultures are prone to epidemic diseases. Rice (Oryza sativa) is susceptible to many diseases, such as rice blast caused by Magnaporthe oryzae. Varietal resistance of rice to blast relies on intracellular nucleotide binding, leucine-rich repeat (NLR) receptors that recognize specific pathogen molecules and trigger immune responses. In the Yuanyang terraces in southwest China, rice landraces rarely show severe losses to disease whereas commercial inbred lines show pronounced field susceptibility. Here, we investigate within-landrace NLR sequence diversity of nine rice landraces and eleven modern varieties using complexity reduction techniques. We find that NLRs display high sequence diversity in landraces, consistent with balancing selection, and that balancing selection at NLRs is more pervasive in landraces than modern varieties. Notably, modern varieties lack many ancient NLR haplotypes that are retained in some landraces. Our study emphasizes the value of standing genetic variation that is maintained in farmer landraces as a resource to make modern crops and agroecosystems less prone to disease. The conservation of landraces is, therefore, crucial for ensuring food security in the face of dynamic biotic and abiotic threats.

Chromosome-Scale Assembly and Annotation of Eight Arabidopsis thaliana Ecotypes

The plant Arabidopsis thaliana is a model system used by researchers through much of plant research. Recent efforts have focused on discovering the genomic variation found in naturally occurring ecotypes isolated from around the world. These ecotypes have come from diverse climates and therefore have faced and adapted to a variety of abiotic and biotic stressors. The sequencing and comparative analysis of these genomes can offer insight into the adaptive strategies of plants. While there are a large number of ecotype genome sequences available, the majority were created using short-read technology. Mapping of short-reads containing structural variation to a reference genome bereft of that variation leads to incorrect mapping of those reads, resulting in a loss of genetic information and introduction of false heterozygosity. For this reason, long-read de novo sequencing of genomes is required to resolve structural variation events. In this article, we sequenced the genomes of eight natural variants of A. thaliana using nanopore sequencing. This resulted in highly contiguous assemblies with >95% of the genome contained within five contigs. The sequencing results from this study include five ecotypes from relict and African populations, an area of untapped genetic diversity. With this study, we increase the knowledge of diversity we have across A. thaliana ecotypes and contribute to ongoing production of an A. thaliana pan-genome.

Plant pangenomes for crop improvement, biodiversity and evolution

Plant genome sequences catalogue genes and the genetic elements that regulate their expression. Such inventories further research aims as diverse as mapping the molecular basis of trait diversity in domesticated plants or inquiries into the origin of evolutionary innovations in flowering plants millions of years ago. The transformative technological progress of DNA sequencing in the past two decades has enabled researchers to sequence ever more genomes with greater ease. Pangenomes - complete sequences of multiple individuals of a species or higher taxonomic unit - have now entered the geneticists' toolkit. The genomes of crop plants and their wild relatives are being studied with translational applications in breeding in mind. But pangenomes are applicable also in ecological and evolutionary studies, as they help classify and monitor biodiversity across the tree of life, deepen our understanding of how plant species diverged and show how plants adapt to changing environments or new selection pressures exerted by human beings.

Phylogenomics of Plant NLR Immune Receptors to Identify Functionally Conserved Sequence Motifs

In recent years, the increase in genome sequencing across diverse plant species has provided a significant advantage for phylogenomics studies, allowing the analysis of one of the most diverse gene families in plants: nucleotide-binding leucine-rich repeat receptors (NLRs). However, due to the sequence diversity of the NLR gene family, identifying key molecular features and functionally conserved sequence patterns is challenging through multiple sequence alignment. Here, we present a step-by-step protocol for a computational pipeline designed to identify evolutionarily conserved motifs in plant NLR proteins. In this protocol, we use a large-scale NLR dataset, including 1,862 NLR genes annotated from monocot and dicot species, to predict conserved sequence motifs, such as the MADA and EDVID motifs, within the coiled-coil (CC)-NLR subfamily. Our pipeline can be applied to identify molecular signatures that have remained conserved in the gene family over evolutionary time across plant species. Key features • Phylogenomics analysis of plant NLR immune receptor family. • Identification of functionally conserved sequence patterns among plant NLRs.

ATR2C ala2 from Arabidopsis-infecting downy mildew requires 4 TIR-NLR immune receptors for full recognition

Arabidopsis Col-0 RPP2A and RPP2B confer recognition of Arabidopsis downy mildew (Hyaloperonospora arabidopsidis [Hpa]) isolate Cala2, but the identity of the recognized ATR2Cala2 effector was unknown. To reveal ATR2Cala2, an F2 population was generated from a cross between Hpa-Cala2 and Hpa-Noks1. We identified ATR2Cala2 as a non-canonical RxLR-type effector that carries a signal peptide, a dEER motif, and WY domains but no RxLR motif. Recognition of ATR2Cala2 and its effector function were verified by biolistic bombardment, ectopic expression and Hpa infection. ATR2Cala2 is recognized in accession Col-0 but not in Ler-0 in which RPP2A and RPP2B are absent. In ATR2Emoy2 and ATR2Noks1 alleles, a frameshift results in an early stop codon. RPP2A and RPP2B are essential for the recognition of ATR2Cala2. Stable and transient expression of ATR2Cala2 under 35S promoter in Arabidopsis and Nicotiana benthamiana enhances disease susceptibility. Two additional Col-0 TIR-NLR (TNL) genes (RPP2C and RPP2D) adjacent to RPP2A and RPP2B are quantitatively required for full resistance to Hpa-Cala2. We compared RPP2 haplotypes in multiple Arabidopsis accessions and showed that all four genes are present in all ATR2Cala2-recognizing accessions.

A telomere-to-telomere gap-free reference genome assembly of avocado provides useful resources for identifying genes related to fatty acid biosynthesis and disease resistance

Avocado (Persea americana Mill.) is an economically valuable plant because of the high fatty acid content and unique flavor of its fruits. Its fatty acid content, especially the relatively high unsaturated fatty acid content, provides significant health benefits. We herein present a telomere-to-telomere gapless genome assembly (841.6 Mb) of West Indian avocado. The genome contains 40 629 predicted protein-coding genes. Repeat sequences account for 57.9% of the genome. Notably, all telomeres, centromeres, and a nucleolar organizing region are included in this genome. Fragments from these three regions were observed via fluorescence in situ hybridization. We identified 376 potential disease resistance-related nucleotide-binding leucine-rich repeat genes. These genes, which are typically clustered on chromosomes, may be derived from gene duplication events. Five NLR genes (Pa11g0262, Pa02g4855, Pa07g3139, Pa07g0383, and Pa02g3196) were highly expressed in leaves, stems, and fruits, indicating they may be involved in avocado disease responses in multiple tissues. We also identified 128 genes associated with fatty acid biosynthesis and analyzed their expression patterns in leaves, stems, and fruits. Pa02g0113, which encodes one of 11 stearoyl-acyl carrier protein desaturases mediating C18 unsaturated fatty acid synthesis, was more highly expressed in the leaves than in the stems and fruits. These findings provide valuable insights that enhance our understanding of fatty acid biosynthesis in avocado.

A chromosome-scale assembly of Brassica carinata (BBCC) accession HC20 containing resistance to multiple pathogens and an early generation assessment of introgressions into B. juncea (AABB)

Brassica carinata (BBCC) commonly referred to as Ethiopian mustard is a natural allotetraploid containing the genomes of Brassica nigra (BB) and Brassica oleracea (CC). It is an oilseed crop endemic to the northeastern regions of Africa. Although it is under limited cultivation, B. carinata is valuable as it is resistant/highly tolerant to most of the pathogens affecting widely cultivated Brassica species of the U's triangle. We report a chromosome-scale genome assembly of B. carinata accession HC20 using long-read Oxford Nanopore sequencing and Bionano optical maps. The assembly has a scaffold N50 of ~39.8 Mb and covers ~1.11 Gb of the genome. We compared the long-read genome assemblies of the U's triangle species and found extensive gene collinearity between the diploids and allopolyploids with no evidence of major gene losses. Therefore, B. juncea (AABB), B. napus (AACC), and B. carinata can be regarded as strict allopolyploids. We cataloged the nucleotide-binding and leucine-rich repeat immune receptor (NLR) repertoire of B. carinata and, identified 465 NLRs, and compared these with the NLRs in the other Brassica species. We investigated the extent and nature of early-generation genomic interactions between the constituent genomes of B. carinata and B. juncea in interspecific crosses between the two species. Besides the expected recombination between the constituent B genomes, extensive homoeologous exchanges were observed between the A and C genomes. Interspecific crosses, therefore, can be used for transferring disease resistance from B. carinata to B. juncea and broadening the genetic base of the two allotetraploid species.

NbPTR1 confers resistance against Pseudomonas syringae pv. actinidiae in kiwifruit

Pseudomonas syringae pv. actinidiae biovar 3 (Psa3) causes a devastating canker disease in yellow-fleshed kiwifruit (Actinidia chinensis). The effector HopZ5, which is present in all isolates of Psa3 causing global outbreaks of pandemic kiwifruit canker disease, triggers immunity in Nicotiana benthamiana and is not recognised in susceptible A. chinensis cultivars. In a search for N. benthamiana nonhost resistance genes against HopZ5, we found that the nucleotide-binding leucine-rich repeat receptor NbPTR1 recognised HopZ5. RPM1-interacting protein 4 orthologues from N. benthamiana and A. chinensis formed a complex with NbPTR1 and HopZ5 activity was able to disrupt this interaction. No functional orthologues of NbPTR1 were found in A. chinensis. NbPTR1 transformed into Psa3-susceptible A. chinensis var. chinensis 'Hort16A' plants introduced HopZ5-specific resistance against Psa3. Altogether, this study suggested that expressing NbPTR1 in Psa3-susceptible kiwifruit is a viable approach to acquiring resistance to Psa3 and it provides valuable information for engineering resistance in otherwise susceptible kiwifruit genotypes.

Transcriptomics and metabolomics analyses provide insights into resistance genes of tree ferns

As ancient organisms, tree ferns play a crucial role as an evolutionary bridge between lower and higher plant species, providing various utilitarian benefits. However, they face challenges such as overexploitation, climate change, adverse environmental conditions, and insect pests, resulting in conservation concerns. In this study, we provide an overview of metabolic and transcriptomic resources of leaves in two typical tree ferns, A. spinulosa and A. metteniana, and explore the resistance genes for the first time. The landscape of metabolome showed that the compound skimmin may hold medicinal significance. A total of 111 differentially accumulated metabolites (DAMs) were detected, with pathway enrichment analysis highlighting 14 significantly enriched pathways, including 2-oxocarboxylic acid metabolism possibly associated with environmental adaptations. A total of 14,639 differentially expressed genes (DEGs) were found, among which 606 were resistance (R) genes. We identified BAM1 as a significantly differentially expressed R gene, which is one of the core genes within the R gene interaction network. Both the maximum-likelihood phylogenetic tree and the PPI network revealed a close relationship between BAM1, FLS2, and TMK. Moreover, BAM1 showed a significant positive correlation with neochlorogenic acid and kaempferol-7-O-glucoside. These metabolites, known for their antioxidant and anti-inflammatory properties, likely play a crucial role in the defense response of tree ferns. This research provides valuable insights into the metabolic and transcriptomic differences between A. spinulosa and A. metteniana, enhancing our understanding of resistance genes in tree ferns.

Copy number variation and population-specific immune genes in the model vertebrate zebrafish

Copy number variation in large gene families is well characterized for plant resistance genes, but similar studies are rare in animals. The zebrafish (Danio rerio) has hundreds of NLR immune genes, making this species ideal for studying this phenomenon. By sequencing 93 zebrafish from multiple wild and laboratory populations, we identified a total of 1513 NLRs, many more than the previously known 400. Approximately half of those are present in all wild populations, but only 4% were found in 80% or more of the individual fish. Wild fish have up to two times as many NLRs per individual and up to four times as many NLRs per population than laboratory strains. In contrast to the massive variability of gene copies, nucleotide diversity in zebrafish NLR genes is very low: around half of the copies are monomorphic and the remaining ones have very few polymorphisms, likely a signature of purifying selection.

Activating plant immunity: the hidden dance of intracellular Ca2+ stores

Calcium ion (Ca2+) serves as a versatile and conserved second messenger in orchestrating immune responses. In plants, plasma membrane-localized Ca2+-permeable channels can be activated to induce Ca2+ influx from extracellular space to cytosol upon pathogen infection. Notably, different immune elicitors can induce dynamic Ca2+ signatures in the cytosol. During pattern-triggered immunity, there is a rapid and transient increase in cytosolic Ca2+, whereas in effector-triggered immunity, the elevation of cytosolic Ca2+ is strong and sustained. Numerous Ca2+ sensors are localized in the cytosol or different intracellular organelles, which are responsible for detecting and converting Ca2+ signals. In fact, Ca2+ signaling coordinated by cytosol and subcellular compartments plays a crucial role in activating plant immune responses. However, the complete Ca2+ signaling network in plant cells is still largely ambiguous. This review offers a comprehensive insight into the collaborative role of intracellular Ca2+ stores in shaping the Ca2+ signaling network during plant immunity, and several intriguing questions for future research are highlighted.

Complete genome assembly provides a high-quality skeleton for pan-NLRome construction in melon

Melon (Cucumis melo L.), being under intensive domestication and selective breeding, displays an abundant phenotypic diversity. Wild germplasm with tolerance to stress represents an untapped genetic resource for discovery of disease-resistance genes. To comprehensively characterize resistance genes in melon, we generate a telomere-to-telomere (T2T) and gap-free genome of wild melon accession PI511890 (C. melo var. chito) with a total length of 375.0 Mb and a contig N50 of 31.24 Mb. The complete genome allows us to dissect genome architecture and identify resistance gene analogs. We construct a pan-NLRome using seven melon genomes, which include 208 variable and 18 core nucleotide-binding leucine-rich repeat receptors (NLRs). Multiple disease-related transcriptome analyses indicate that most up-regulated NLRs induced by pathogens are shell or cloud NLRs. The T2T gap-free assembly and the pan-NLRome not only serve as essential resources for genomic studies and molecular breeding of melon but also provide insights into the genome architecture and NLR diversity.

Population comparative genomics discovers gene gain and loss during grapevine domestication

Plant domestication are evolutionary experiments conducted by early farmers since thousands years ago, during which the crop wild progenitors are artificially selected for desired agronomic traits along with dramatic genomic variation in the course of moderate to severe bottlenecks. However, previous investigations are mainly focused on small-effect variants, while changes in gene contents are rarely investigated due to the lack of population-level assemblies for both the crop and its wild relatives. Here, we applied comparative genomic analyses to discover gene gain and loss during grapevine domestication using long-read assemblies of representative population samples for both domesticated grapevines (V. vinifera ssp. vinifera) and their wild progenitors (V. vinifera ssp. sylvestris). Only ∼7% of gene families were shared by 16 Vitis genomes while ∼8% of gene families were specific to each accession, suggesting dramatic variations of gene contents in grapevine genomes. Compared to wild progenitors, the domesticated accessions exhibited an increased presence of genes associated with asexual reproduction, while the wild progenitors showcased a higher abundance of genes related to pollination, revealing the transition from sexual reproduction to clonal propagation during domestication processes. Moreover, the domesticated accessions harbored fewer disease-resistance genes than wild progenitors. The SVs occurred frequently in aroma and disease-resistance related genes between domesticated grapevines and wild progenitors, indicating the rapid diversification of these genes during domestication. Our study provides insights and resources for biological studies and breeding programs in grapevine.

An NLR paralog Pit2 generated from tandem duplication of Pit1 fine-tunes Pit1 localization and function

NLR family proteins act as intracellular receptors. Gene duplication amplifies the number of NLR genes, and subsequent mutations occasionally provide modifications to the second gene that benefits immunity. However, evolutionary processes after gene duplication and functional relationships between duplicated NLRs remain largely unclear. Here, we report that the rice NLR protein Pit1 is associated with its paralogue Pit2. The two are required for the resistance to rice blast fungus but have different functions: Pit1 induces cell death, while Pit2 competitively suppresses Pit1-mediated cell death. During evolution, the suppression of Pit1 by Pit2 was probably generated through positive selection on two fate-determining residues in the NB-ARC domain of Pit2, which account for functional differences between Pit1 and Pit2. Consequently, Pit2 lost its plasma membrane localization but acquired a new function to interfere with Pit1 in the cytosol. These findings illuminate the evolutionary trajectory of tandemly duplicated NLR genes after gene duplication.

Are cereal grasses a single genetic system?

In 1993, a passionate and provocative call to arms urged cereal researchers to consider the taxon they study as a single genetic system and collaborate with each other. Since then, that group of scientists has seen their discipline blossom. In an attempt to understand what unity of genetic systems means and how the notion was borne out by later research, we survey the progress and prospects of cereal genomics: sequence assemblies, population-scale sequencing, resistance gene cloning and domestication genetics. Gene order may not be as extraordinarily well conserved in the grasses as once thought. Still, several recurring themes have emerged. The same ancestral molecular pathways defining plant architecture have been co-opted in the evolution of different cereal crops. Such genetic convergence as much as cross-fertilization of ideas between cereal geneticists has led to a rich harvest of genes that, it is hoped, will lead to improved varieties.

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

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

High allelic diversity in Arabidopsis NLRs is associated with distinct genomic features

Plants rely on Nucleotide-binding, Leucine-rich repeat Receptors (NLRs) for pathogen recognition. Highly variable NLRs (hvNLRs) show remarkable intraspecies diversity, while their low-variability paralogs (non-hvNLRs) are conserved between ecotypes. At a population level, hvNLRs provide new pathogen-recognition specificities, but the association between allelic diversity and genomic and epigenomic features has not been established. Our investigation of NLRs in Arabidopsis Col-0 has revealed that hvNLRs show higher expression, less gene body cytosine methylation, and closer proximity to transposable elements than non-hvNLRs. hvNLRs show elevated synonymous and nonsynonymous nucleotide diversity and are in chromatin states associated with an increased probability of mutation. Diversifying selection maintains variability at a subset of codons of hvNLRs, while purifying selection maintains conservation at non-hvNLRs. How these features are established and maintained, and whether they contribute to the observed diversity of hvNLRs is key to understanding the evolution of plant innate immune receptors.

Pattern recognition receptors as potential therapeutic targets for developing immunological engineered plants

There is an urgent need to combat pathogen infestations in crop plants to ensure food security worldwide. To counter this, plants have developed innate immunity mediated by Pattern Recognition Receptors (PRRs) that recognize pathogen-associated molecular patterns (PAMPs) and damage- associated molecular patterns (DAMPs). PRRs activate Pattern-Triggered Immunity (PTI), a defence mechanism involving intricate cell-surface and intracellular receptors. The diverse ligand-binding ectodomains of PRRs, including leucine-rich repeats (LRRs) and lectin domains, facilitate the recognition of MAMPs and DAMPs. Pathogen resistance is mediated by a variety of PTI responses, including membrane depolarization, ROS production, and the induction of defence genes. An integral part of intracellular immunity is the Nucleotide-binding Oligomerization Domain, Leucine-rich Repeat proteins (NLRs) which recognize and respond to effectors in a potent manner. Enhanced understanding of PRRs, their ligands, and downstream signalling pathways has contributed to the identification of potential targets for genetically modified plants. By transferring PRRs across plant species, it is possible to create broad-spectrum resistance, potentially offering innovative solutions for plant protection and global food security. The purpose of this chapter is to provide an update on PRRs involved in disease resistance, clarify the mechanisms by which PRRs recognize ligands to form active receptor complexes and present various applications of PRRs and PTI in disease resistance management for plants.

The plant immune system: From discovery to deployment

Plant diseases cause famines, drive human migration, and present challenges to agricultural sustainability as pathogen ranges shift under climate change. Plant breeders discovered Mendelian genetic loci conferring disease resistance to specific pathogen isolates over 100 years ago. Subsequent breeding for disease resistance underpins modern agriculture and, along with the emergence and focus on model plants for genetics and genomics research, has provided rich resources for molecular biological exploration over the last 50 years. These studies led to the identification of extracellular and intracellular receptors that convert recognition of extracellular microbe-encoded molecular patterns or intracellular pathogen-delivered virulence effectors into defense activation. These receptor systems, and downstream responses, define plant immune systems that have evolved since the migration of plants to land ∼500 million years ago. Our current understanding of plant immune systems provides the platform for development of rational resistance enhancement to control the many diseases that continue to plague crop production.

The NLR family of innate immune and cell death sensors

Nucleotide-binding oligomerization domain (NOD)-like receptors, also known as nucleotide-binding leucine-rich repeat receptors (NLRs), are a family of cytosolic pattern recognition receptors that detect a wide variety of pathogenic and sterile triggers. Activation of specific NLRs initiates pro- or anti-inflammatory signaling cascades and the formation of inflammasomes-multi-protein complexes that induce caspase-1 activation to drive inflammatory cytokine maturation and lytic cell death, pyroptosis. Certain NLRs and inflammasomes act as integral components of larger cell death complexes-PANoptosomes-driving another form of lytic cell death, PANoptosis. Here, we review the current understanding of the evolution, structure, and function of NLRs in health and disease. We discuss the concept of NLR networks and their roles in driving cell death and immunity. An improved mechanistic understanding of NLRs may provide therapeutic strategies applicable across infectious and inflammatory diseases and in cancer.

Auto-inhibition and activation of a short Argonaute-associated TIR-APAZ defense system

Short prokaryotic Ago accounts for most prokaryotic Argonaute proteins (pAgos) and is involved in defending bacteria against invading nucleic acids. Short pAgo associated with TIR-APAZ (SPARTA) has been shown to oligomerize and deplete NAD+ upon guide-mediated target DNA recognition. However, the molecular basis of SPARTA inhibition and activation remains unknown. In this study, we determined the cryogenic electron microscopy structures of Crenotalea thermophila SPARTA in its inhibited, transient and activated states. The SPARTA monomer is auto-inhibited by its acidic tail, which occupies the guide-target binding channel. Guide-mediated target binding expels this acidic tail and triggers substantial conformational changes to expose the Ago-Ago dimerization interface. As a result, SPARTA assembles into an active tetramer, where the four TIR domains are rearranged and packed to form NADase active sites. Together with biochemical evidence, our results provide a panoramic vision explaining SPARTA auto-inhibition and activation and expand understanding of pAgo-mediated bacterial defense systems.

High-resolution mapping of Ryd4Hb, a major resistance gene to Barley yellow dwarf virus from Hordeum bulbosum

KEY MESSAGE

We mapped Ryd4Hb in a 66.5 kbp interval in barley and dissociated it from a sublethality factor. These results will enable a targeted selection of the resistance in barley breeding. Virus diseases are causing high yield losses in crops worldwide. The Barley yellow dwarf virus (BYDV) complex is responsible for one of the most widespread and economically important viral diseases of cereals. While no gene conferring complete resistance (immunity) has been uncovered in the primary gene pool of barley, sources of resistance were searched and identified in the wild relative Hordeum bulbosum, representing the secondary gene pool of barley. One such locus, Ryd4Hb, has been previously introgressed into barley, and was allocated to chromosome 3H, but is tightly linked to a sublethality factor that prevents the incorporation and utilization of Ryd4Hb in barley varieties. To solve this problem, we fine-mapped Ryd4Hb and separated it from this negative factor. We narrowed the Ryd4Hb locus to a corresponding 66.5 kbp physical interval in the barley 'Morex' reference genome. The region comprises a gene from the nucleotide-binding and leucine-rich repeat immune receptor family, typical of dominant virus resistance genes. The closest homolog to this Ryd4Hb candidate gene is the wheat Sr35 stem rust resistance gene. In addition to the fine mapping, we reduced the interval bearing the sublethality factor to 600 kbp in barley. Aphid feeding experiments demonstrated that Ryd4Hb provides a resistance to BYDV rather than to its vector. The presented results, including the high-throughput molecular markers, will permit a more targeted selection of the resistance in breeding, enabling the use of Ryd4Hb in barley varieties.

The genome and transcriptome of the snail Biomphalaria sudanica s.l.: immune gene diversification and highly polymorphic genomic regions in an important African vector of Schistosoma mansoni

BACKGROUND

Control and elimination of schistosomiasis is an arduous task, with current strategies proving inadequate to break transmission. Exploration of genetic approaches to interrupt Schistosoma mansoni transmission, the causative agent for human intestinal schistosomiasis in sub-Saharan Africa and South America, has led to genomic research of the snail vector hosts of the genus Biomphalaria. Few complete genomic resources exist, with African Biomphalaria species being particularly underrepresented despite this being where the majority of S. mansoni infections occur. Here we generate and annotate the first genome assembly of Biomphalaria sudanica sensu lato, a species responsible for S. mansoni transmission in lake and marsh habitats of the African Rift Valley. Supported by whole-genome diversity data among five inbred lines, we describe orthologs of immune-relevant gene regions in the South American vector B. glabrata and present a bioinformatic pipeline to identify candidate novel pathogen recognition receptors (PRRs).

RESULTS

De novo genome and transcriptome assembly of inbred B. sudanica originating from the shoreline of Lake Victoria (Kisumu, Kenya) resulted in a haploid genome size of ~ 944.2 Mb (6,728 fragments, N50 = 1.067 Mb), comprising 23,598 genes (BUSCO = 93.6% complete). The B. sudanica genome contains orthologues to all described immune genes/regions tied to protection against S. mansoni in B. glabrata, including the polymorphic transmembrane clusters (PTC1 and PTC2), RADres, and other loci. The B. sudanica PTC2 candidate immune genomic region contained many PRR-like genes across a much wider genomic region than has been shown in B. glabrata, as well as a large inversion between species. High levels of intra-species nucleotide diversity were seen in PTC2, as well as in regions linked to PTC1 and RADres orthologues. Immune related and putative PRR gene families were significantly over-represented in the sub-set of B. sudanica genes determined as hyperdiverse, including high extracellular diversity in transmembrane genes, which could be under pathogen-mediated balancing selection. However, no overall expansion in immunity related genes was seen in African compared to South American lineages.

CONCLUSIONS

The B. sudanica genome and analyses presented here will facilitate future research in vector immune defense mechanisms against pathogens. This genomic/transcriptomic resource provides necessary data for the future development of molecular snail vector control/surveillance tools, facilitating schistosome transmission interruption mechanisms in Africa.

Paired plant immune CHS3-CSA1 receptor alleles form distinct hetero-oligomeric complexes

Plant intracellular nucleotide-binding leucine-rich repeat receptors (NLRs) analyzed to date oligomerize and form resistosomes upon activation to initiate immune responses. Some NLRs are encoded in tightly linked co-regulated head-to-head genes whose products function together as pairs. We uncover the oligomerization requirements for different Arabidopsis paired CHS3-CSA1 alleles. These pairs form resting-state heterodimers that oligomerize into complexes distinct from NLRs analyzed previously. Oligomerization requires both conserved and allele-specific features of the respective CHS3 and CSA1 Toll-like interleukin-1 receptor (TIR) domains. The receptor kinases BAK1 and BIRs inhibit CHS3-CSA1 pair oligomerization to maintain the CHS3-CSA1 heterodimer in an inactive state. Our study reveals that paired NLRs hetero-oligomerize and likely form a distinctive "dimer of heterodimers" and that structural heterogeneity is expected even among alleles of closely related paired NLRs.

Two haplotype-resolved genome assemblies for AAB allotriploid bananas provide insights into banana subgenome asymmetric evolution and Fusarium wilt control

Bananas (Musa spp.) are one of the world's most important fruit crops and play a vital role in food security for many developing countries. Most banana cultivars are triploids derived from inter- and intraspecific hybridizations between the wild diploid ancestor species Musa acuminate (AA) and M. balbisiana (BB). We report two haplotype-resolved genome assemblies of the representative AAB-cultivated types, Plantain and Silk, and precisely characterize ancestral contributions by examining ancestry mosaics across the genome. Widespread asymmetric evolution is observed in their subgenomes, which can be linked to frequent homologous exchange events. We reveal the genetic makeup of triploid banana cultivars and verify that subgenome B is a rich source of disease resistance genes. Only 58.5% and 59.4% of Plantain and Silk genes, respectively, are present in all three haplotypes, with >50% of genes being differentially expressed alleles in different subgenomes. We observed that the number of upregulated genes in Plantain is significantly higher than that in Silk at one-week post-inoculation with Fusarium wilt tropical race 4 (Foc TR4), which confirms that Plantain can initiate defense responses faster than Silk. Additionally, we compared genomic and transcriptomic differences among the genes related to carotenoid synthesis and starch metabolism between Plantain and Silk. Our study provides resources for better understanding the genomic architecture of cultivated bananas and has important implications for Musa genetics and breeding.

Haplotype-resolved genome assembly provides insights into evolutionary history of the Actinidia arguta tetraploid

Actinidia arguta, known as hardy kiwifruit, is a widely cultivated species with distinct botanical characteristics such as small and smooth-fruited, rich in beneficial nutrients, rapid softening and tolerant to extremely low temperatures. It contains the most diverse ploidy types, including diploid, tetraploid, hexaploid, octoploid, and decaploid. Here we report a haplotype-resolved tetraploid genome (A. arguta cv. 'Longcheng No.2') containing four haplotypes, each with 40,859, 41,377, 39,833 and 39,222 protein-coding genes. We described the phased genome structure, synteny, and evolutionary analyses to identify and date possible WGD events. Ks calculations for both allelic and paralogous genes pairs throughout the assembled haplotypic individuals showed its tetraploidization is estimated to have formed ~ 1.03 Mya following Ad-α event occurred ~ 18.7 Mya. Detailed annotations of NBS-LRRs or CBFs highlight the importance of genetic variations coming about after polyploidization in underpinning ability of immune responses or environmental adaptability. WGCNA analysis of postharvest quality indicators in combination with transcriptome revealed several transcription factors were involved in regulating ripening kiwi berry texture. Taking together, the assembly of an A. arguta tetraploid genome provides valuable resources in deciphering complex genome structure and facilitating functional genomics studies and genetic improvement for kiwifruit and other crops.

The NLR immune receptor ADR1 and lipase-like proteins EDS1 and PAD4 mediate stomatal immunity in Nicotiana benthamiana and Arabidopsis

In the presence of pathogenic bacteria, plants close their stomata to prevent pathogen entry. Intracellular nucleotide-binding leucine-rich repeat (NLR) immune receptors recognize pathogenic effectors and activate effector-triggered immune responses. However, the regulatory and molecular mechanisms of stomatal immunity involving NLR immune receptors are unknown. Here, we show that the Nicotiana benthamiana RPW8-NLR central immune receptor ACTIVATED DISEASE RESISTANCE 1 (NbADR1), together with the key immune proteins ENHANCED DISEASE SUSCEPTIBILITY 1 (NbEDS1) and PHYTOALEXIN DEFICIENT 4 (NbPAD4), plays an essential role in bacterial pathogen- and flg22-induced stomatal immunity by regulating the expression of salicylic acid (SA) and abscisic acid (ABA) biosynthesis or response-related genes. NbADR1 recruits NbEDS1 and NbPAD4 in stomata to form a stomatal immune response complex. The transcription factor NbWRKY40e, in association with NbEDS1 and NbPAD4, modulates the expression of SA and ABA biosynthesis or response-related genes to influence stomatal immunity. NbADR1, NbEDS1, and NbPAD4 are required for the pathogen infection-enhanced binding of NbWRKY40e to the ISOCHORISMATE SYNTHASE 1 promoter. Moreover, the ADR1-EDS1-PAD4 module regulates stomatal immunity in Arabidopsis (Arabidopsis thaliana). Collectively, our findings show the pivotal role of the core intracellular immune receptor module ADR1-EDS1-PAD4 in stomatal immunity, which enables plants to limit pathogen entry.

Genetic basis of Arabidopsis thaliana responses to infection by naïve and adapted isolates of turnip mosaic virus

Plant viruses account for enormous agricultural losses worldwide, and the most effective way to combat them is to identify genetic material conferring plant resistance to these pathogens. Aiming to identify genetic associations with responses to infection, we screened a large panel of Arabidopsis thaliana natural inbred lines for four disease-related traits caused by infection by A. thaliana-naïve and -adapted isolates of the natural pathogen turnip mosaic virus (TuMV). We detected a strong, replicable association in a 1.5 Mb region on chromosome 2 with a 10-fold increase in relative risk of systemic necrosis. The region contains several plausible causal genes as well as abundant structural variation, including an insertion of a Copia transposon into a Toll/interleukin receptor (TIR-NBS-LRR) coding for a gene involved in defense, that could be either a driver or a consequence of the disease-resistance locus. When inoculated with TuMV, loss-of-function mutant plants of this gene exhibited different symptoms than wild-type plants. The direction and severity of symptom differences depended on the adaptation history of the virus. This increase in symptom severity was specific for infections with the adapted isolate. Necrosis-associated alleles are found worldwide, and their distribution is consistent with a trade-off between resistance during viral outbreaks and a cost of resistance otherwise, leading to negative frequency-dependent selection.

Differential expansion and retention patterns of LRR-RLK genes across plant evolution

To maximize overall fitness, plants must accurately respond to a host of growth, developmental, and environmental signals throughout their life. Many of these internal and external signals are perceived by the leucine-rich repeat receptor-like kinases, which play roles in regulating growth, development, and immunity. This largest family of receptor kinases in plants can be divided into subfamilies based on the conservation of the kinase domain, which demonstrates that shared evolutionary history often indicates shared molecular function. Here we investigate the evolutionary history of this family across the evolution of 112 plant species. We identify lineage-specific expansions of the malectin-domain containing subfamily LRR subfamily I primarily in the Brassicales and bryophytes. Most other plant lineages instead show a large expansion in LRR subfamily XII, which in Arabidopsis is known to contain key receptors in pathogen perception. This striking asymmetric expansion may reveal a dichotomy in the evolutionary history and adaptation strategies employed by plants. A greater understanding of the evolutionary pressures and adaptation strategies acting on members of this receptor family offers a way to improve functional predictions for orphan receptors and simplify the identification of novel stress-related receptors.

The structured coalescent in the context of gene copy number variation

The Structured Coalescent was introduced to describe the coalescent process in spatially subdivided populations with migration. Here, we re-interpret migration routes of individuals in the original model as "migration routes" of single genes in tandemly arranged gene arrays. A gene copy may change its position within the array via unequal recombination. Hence, in a coalescent framework, two copies sampled from two chromosomes may coalesce only if they are at exactly homologous positions. Otherwise, one or multiple recombination events have to occur before they can coalesce, thereby increasing mean coalescence time and expected genetic diversity among the copies in a gene array. We explicitly calculate the transition probabilities on these routes backward in time. We simulate the structured coalescent with migration and coalescence rates informed by the unequal recombination process of gene copies. With this novel interpretation of population structure models we determine coalescence times and expected genetic diversity in samples of orthologous and paralogous copies from a gene family. As a case study, we discuss the site frequency spectrum of a small gene family in the two scenarios of high and of no gene copy number variation among individuals. These examples underline the significance of our model, since standard test-statistics may lead to misinterpretations when analyzing sequence data of multi-copy genes due to their different expected genetic diversity.

The plant immune receptor SNC1 monitors helper NLRs targeted by a bacterial effector

Plants deploy intracellular receptors to counteract pathogen effectors that suppress cell-surface-receptor-mediated immunity. To what extent pathogens manipulate intracellular receptor-mediated immunity, and how plants tackle such manipulation, remains unknown. Arabidopsis thaliana encodes three similar ADR1 class helper nucleotide-binding domain leucine-rich repeat receptors (ADR1, ADR1-L1, and ADR1-L2), which are crucial in plant immunity initiated by intracellular receptors. Here, we report that Pseudomonas syringae effector AvrPtoB suppresses ADR1-L1- and ADR1-L2-mediated cell death. ADR1, however, evades such suppression by diversifying into two ubiquitination sites targeted by AvrPtoB. The intracellular sensor SNC1 interacts with and guards the CCR domains of ADR1-L1/L2. Removal of ADR1-L1/L2 or delivery of AvrPtoB activates SNC1, which then signals through ADR1 to trigger immunity. Our work elucidates the long-sought-after function of SNC1 in defense, and also how plants can use dual strategies, sequence diversification, and a multi-layered guard-guardee system, to counteract pathogen's attack on core immunity functions.

The genome and transcriptome of the snail Biomphalaria sudanica s.l.: Immune gene diversification and highly polymorphic genomic regions in an important African vector of Schistosoma mansoni

Background

Control and elimination of schistosomiasis is an arduous task, with current strategies proving inadequate to break transmission. Exploration of genetic approaches to interrupt Schistosoma mansoni transmission, the causative agent for human intestinal schistosomiasis in sub-Saharan Africa and South America, has led to genomic research of the snail vector hosts of the genus Biomphalaria. Few complete genomic resources exist, with African Biomphalaria species being particularly underrepresented despite this being where the majority of S. mansoni infections occur. Here we generate and annotate the first genome assembly of Biomphalaria sudanica sensu lato, a species responsible for S. mansoni transmission in lake and marsh habitats of the African Rift Valley. Supported by whole-genome diversity data among five inbred lines, we describe orthologs of immune-relevant gene regions in the South American vector B. glabrata and present a bioinformatic pipeline to identify candidate novel pathogen recognition receptors (PRRs).

Results

De novo genome and transcriptome assembly of inbred B. sudanica originating from the shoreline of Lake Victoria (Kisumu, Kenya) resulted in a haploid genome size of ~944.2 Mb (6732 fragments, N50=1.067 Mb), comprising 23,598 genes (BUSCO=93.6% complete). The B. sudanica genome contains orthologues to all described immune genes/regions tied to protection against S. mansoni in B. glabrata. The B. sudanica PTC2 candidate immune genomic region contained many PRR-like genes across a much wider genomic region than has been shown in B. glabrata, as well as a large inversion between species. High levels of intra-species nucleotide diversity were seen in PTC2, as well as in regions linked to PTC1 and RADres orthologues. Immune related and putative PRR gene families were significantly over-represented in the sub-set of B. sudanica genes determined as hyperdiverse, including high extracellular diversity in transmembrane genes, which could be under pathogen-mediated balancing selection. However, no overall expansion in immunity related genes were seen in African compared to South American lineages.

Conclusions

The B. sudanica genome and analyses presented here will facilitate future research in vector immune defense mechanisms against pathogens. This genomic/transcriptomic resource provides necessary data for the future development of molecular snail vector control/surveillance tools, facilitating schistosome transmission interruption mechanisms in Africa.

SMRT-AgRenSeq-d in potato (Solanum tuberosum) as a method to identify candidates for the nematode resistance Gpa5

Potato is the third most important food crop in the world. Diverse pathogens threaten sustainable crop production but can be controlled, in many cases, through the deployment of disease resistance genes belonging to the family of nucleotide-binding, leucine-rich-repeat (NLR) genes. To identify effective disease resistance genes in established varieties, we have successfully established SMRT-AgRenSeq in tetraploid potatoes and have further enhanced the methodology by including dRenSeq in an approach that we term SMR-AgRenSeq-d. The inclusion of dRenSeq enables the filtering of candidates after the association analysis by establishing a presence/absence matrix across resistant and susceptible varieties that is translated into an F1 score. Using a SMRT-RenSeq-based sequence representation of the NLRome from the cultivar Innovator, SMRT-AgRenSeq-d analyses reliably identified the late blight resistance benchmark genes Rpi-R1, Rpi-R2-like, Rpi-R3a, and Rpi-R3b in a panel of 117 varieties with variable phenotype penetrations. All benchmark genes were identified with an F1 score of 1, which indicates absolute linkage in the panel. This method also identified nine strong candidates for Gpa5 that controls the potato cyst nematode (PCN) species Globodera pallida (pathotypes Pa2/3). Assuming that NLRs are involved in controlling many types of resistances, SMRT-AgRenSeq-d can readily be applied to diverse crops and pathogen systems.

Global expression patterns of R-genes in tomato and potato

Introduction

As key-players of plant immunity, the proteins encoded by resistance genes (R-genes) recognize pathogens and initiate pathogen-specific defense responses. The expression of some R-genes carry fitness costs and therefore inducible immune responses are likely advantageous. To what degree inducible resistance driven by R-genes is triggered by pathogen infection is currently an open question.

Methods

In this study we analyzed the expression of 940 R-genes of tomato and potato across 315 transcriptome libraries to investigate how interspecific interactions with microbes influence R-gene expression in plants.

Results

We found that most R-genes are expressed at a low level. A small subset of R-genes had moderate to high levels of expression and were expressed across many independent libraries, irrespective of infection status. These R-genes include members of the class of genes called NRCs (NLR required for cell death). Approximately 10% of all R-genes were differentially expressed during infection and this included both up- and down-regulation. One factor associated with the large differences in R-gene expression was host tissue, reflecting a considerable degree of tissue-specific transcriptional regulation of this class of genes.

Discussion

These results call into question the widespread view that R-gene expression is induced upon pathogen attack. Instead, a small core set of R-genes is constitutively expressed, imparting upon the plant a ready-to-detect and defend status.

Plant immune receptors interact with hemibiotrophic pathogens to activate plant immunity

Phytopathogens pose a devastating threat to the productivity and yield of crops by causing destructive plant diseases in natural and agricultural environments. Hemibiotrophic pathogens have a variable-length biotrophic phase before turning to necrosis and are among the most invasive plant pathogens. Plant resistance to hemibiotrophic pathogens relies mainly on the activation of innate immune responses. These responses are typically initiated after the plant plasma membrane and various plant immune receptors detect immunogenic signals associated with pathogen infection. Hemibiotrophic pathogens evade pathogen-triggered immunity by masking themselves in an arms race while also enhancing or manipulating other receptors to promote virulence. However, our understanding of plant immune defenses against hemibiotrophic pathogens is highly limited due to the intricate infection mechanisms. In this review, we summarize the strategies that different hemibiotrophic pathogens interact with host immune receptors to activate plant immunity. We also discuss the significant role of the plasma membrane in plant immune responses, as well as the current obstacles and potential future research directions in this field. This will enable a more comprehensive understanding of the pathogenicity of hemibiotrophic pathogens and how distinct plant immune receptors oppose them, delivering valuable data for the prevention and management of plant diseases.

Crop adaptation to climate change: An evolutionary perspective

The disciplines of evolutionary biology and plant and animal breeding have been intertwined throughout their development, with responses to artificial selection yielding insights into the action of natural selection and evolutionary biology providing statistical and conceptual guidance for modern breeding. Here we offer an evolutionary perspective on a grand challenge of the 21st century: feeding humanity in the face of climate change. We first highlight promising strategies currently under way to adapt crops to current and future climate change. These include methods to match crop varieties with current and predicted environments and to optimize breeding goals, management practices, and crop microbiomes to enhance yield and sustainable production. We also describe the promise of crop wild relatives and recent technological innovations such as speed breeding, genomic selection, and genome editing for improving environmental resilience of existing crop varieties or for developing new crops. Next, we discuss how methods and theory from evolutionary biology can enhance these existing strategies and suggest novel approaches. We focus initially on methods for reconstructing the evolutionary history of crops and their pests and symbionts, because such historical information provides an overall framework for crop-improvement efforts. We then describe how evolutionary approaches can be used to detect and mitigate the accumulation of deleterious mutations in crop genomes, identify alleles and mutations that underlie adaptation (and maladaptation) to agricultural environments, mitigate evolutionary trade-offs, and improve critical proteins. Continuing feedback between the evolution and crop biology communities will ensure optimal design of strategies for adapting crops to climate change.

Towards integrative plant pathology

The field of plant pathology has revealed many of the mechanisms underlying the arms race, providing crucial knowledge and genetic resources for improving plant health. Although the host-microbe interaction seemingly favors rapidly evolving pathogens, it has also generated a vast evolutionary history of largely unexplored plant immunodiversity. We review studies that characterize the scope and distribution of genetic and ecological diversity in model and non-model systems with specific reference to pathogen effector diversity, plant immunodiversity in both cultivated species and their wild relatives, and diversity in the plant-associated microbiota. We show how the study of evolutionary and ecological processes can reveal patterns of genetic convergence, conservation, and diversification, and that this diversity is increasingly tractable in both experimental and translational systems. Perhaps most importantly, these patterns of diversity provide largely untapped resources that can be deployed for the rational engineering of durable resistance for sustainable agriculture.

Jurassic NLR: Conserved and dynamic evolutionary features of the atypically ancient immune receptor ZAR1

Plant nucleotide-binding leucine-rich repeat (NLR) immune receptors generally exhibit hallmarks of rapid evolution, even at the intraspecific level. We used iterative sequence similarity searches coupled with phylogenetic analyses to reconstruct the evolutionary history of HOPZ-ACTIVATED RESISTANCE1 (ZAR1), an atypically conserved NLR that traces its origin to early flowering plant lineages ∼220 to 150 million yrs ago (Jurassic period). We discovered 120 ZAR1 orthologs in 88 species, including the monocot Colocasia esculenta, the magnoliid Cinnamomum micranthum, and most eudicots, notably the Ranunculales species Aquilegia coerulea, which is outside the core eudicots. Ortholog sequence analyses revealed highly conserved features of ZAR1, including regions for pathogen effector recognition and cell death activation. We functionally reconstructed the cell death activity of ZAR1 and its partner receptor-like cytoplasmic kinase (RLCK) from distantly related plant species, experimentally validating the hypothesis that ZAR1 evolved to partner with RLCKs early in its evolution. In addition, ZAR1 acquired novel molecular features. In cassava (Manihot esculenta) and cotton (Gossypium spp.), ZAR1 carries a C-terminal thioredoxin-like domain, and in several taxa, ZAR1 duplicated into 2 paralog families, which underwent distinct evolutionary paths. ZAR1 stands out among angiosperm NLR genes for having experienced relatively limited duplication and expansion throughout its deep evolutionary history. Nonetheless, ZAR1 also gave rise to noncanonical NLRs with integrated domains and degenerated molecular features.

An atypical NLR gene confers bacterial wilt susceptibility in Arabidopsis

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

Co-evolution within the plant holobiont drives host performance

Plants interact with a diversity of microorganisms that influence their growth and resilience, and they can therefore be considered as ecological entities, namely "plant holobionts," rather than as singular organisms. In a plant holobiont, the assembly of above- and belowground microbiota is ruled by host, microbial, and environmental factors. Upon microorganism perception, plants activate immune signaling resulting in the secretion of factors that modulate microbiota composition. Additionally, metabolic interdependencies and antagonism between microbes are driving forces for community assemblies. We argue that complex plant-microbe and intermicrobial interactions have been selected for during evolution and may promote the survival and fitness of plants and their associated microorganisms as holobionts. As part of this process, plants evolved metabolite-mediated strategies to selectively recruit beneficial microorganisms in their microbiota. Some of these microbiota members show host-adaptation, from which mutualism may rapidly arise. In the holobiont, microbiota members also co-evolved antagonistic activities that restrict proliferation of microbes with high pathogenic potential and can therefore prevent disease development. Co-evolution within holobionts thus ultimately drives plant performance.

Solanum americanum genome-assisted discovery of immune receptors that detect potato late blight pathogen effectors

Potato (Solanum tuberosum) and tomato (Solanum lycopersicon) crops suffer severe losses to late blight caused by the oomycete pathogen Phytophthora infestans. Solanum americanum, a relative of potato and tomato, is globally distributed and most accessions are highly blight resistant. We generated high-quality reference genomes of four S. americanum accessions, resequenced 52 accessions, and defined a pan-NLRome of S. americanum immune receptor genes. We further screened for variation in recognition of 315P. infestans RXLR effectors in 52 S. americanum accessions. Using these genomic and phenotypic data, we cloned three NLR-encoding genes, Rpi-amr4, R02860 and R04373, that recognize cognate P. infestans RXLR effectors PITG_22825 (AVRamr4), PITG_02860 and PITG_04373. These genomic resources and methodologies will support efforts to engineer potatoes with durable late blight resistance and can be applied to diseases of other crops.

Ptr1 and ZAR1 immune receptors confer overlapping and distinct bacterial pathogen effector specificities

Some nucleotide-binding and leucine-rich repeat receptors (NLRs) indirectly detect pathogen effectors by monitoring their host targets. In Arabidopsis thaliana, RIN4 is targeted by multiple sequence-unrelated effectors and activates immune responses mediated by RPM1 and RPS2. These effectors trigger cell death in Nicotiana benthamiana, but the corresponding NLRs have yet not been identified. To identify N. benthamiana NLRs (NbNLRs) that recognize Arabidopsis RIN4-targeting effectors, we conducted a rapid reverse genetic screen using an NbNLR VIGS library. We identified that the N. benthamiana homolog of Ptr1 (Pseudomonas tomato race 1) recognizes the Pseudomonas effectors AvrRpt2, AvrRpm1, and AvrB. We demonstrated that recognition of the Xanthomonas effector AvrBsT and the Pseudomonas effector HopZ5 is conferred independently by the N. benthamiana homolog of Ptr1 and ZAR1. Interestingly, the recognition of HopZ5 and AvrBsT is contributed unequally by Ptr1 and ZAR1 in N. benthamiana and Capsicum annuum. In addition, we showed that the RLCK XII family protein JIM2 is required for the NbZAR1-dependent recognition of AvrBsT and HopZ5. The recognition of sequence-unrelated effectors by NbPtr1 and NbZAR1 provides an additional example of convergently evolved effector recognition. Identification of key components involved in Ptr1 and ZAR1-mediated immunity could reveal unique mechanisms of expanded effector recognition.

Vitis labrusca genome assembly reveals diversification between wild and cultivated grapevine genomes

Wild grapevines are important genetic resources in breeding programs to confer adaptive fitness traits and unique fruit characteristics, but the genetics underlying these traits, and their evolutionary origins, are largely unknown. To determine the factors that contributed to grapevine genome diversification, we performed comprehensive intragenomic and intergenomic analyses with three cultivated European (including the PN40024 reference genome) and two wild North American grapevine genomes, including our newly released Vitis labrusca genome. We found the heterozygosity of the cultivated grapevine genomes was twice as high as the wild grapevine genomes studied. Approximately 30% of V. labrusca and 48% of V. vinifera Chardonnay genes were heterozygous or hemizygous and a considerable number of collinear genes between Chardonnay and V. labrusca had different gene zygosity. Our study revealed evidence that supports gene gain-loss events in parental genomes resulted in the inheritance of hemizygous genes in the Chardonnay genome. Thousands of segmental duplications supplied source material for genome-specific genes, further driving diversification of the genomes studied. We found an enrichment of recently duplicated, adaptive genes in similar functional pathways, but differential retention of environment-specific adaptive genes within each genome. For example, large expansions of NLR genes were discovered in the two wild grapevine genomes studied. Our findings support variation in transposable elements contributed to unique traits in grapevines. Our work revealed gene zygosity, segmental duplications, gene gain-and-loss variations, and transposable element polymorphisms can be key driving forces for grapevine genome diversification.

Extrachromosomal circular DNA and structural variants highlight genome instability in Arabidopsis epigenetic mutants

Abundant extrachromosomal circular DNA (eccDNA) is associated with transposable element (TE) activity. However, how the eccDNA compartment is controlled by epigenetic regulations and what is its impact on the genome is understudied. Here, using long reads, we sequence both the eccDNA compartment and the genome of Arabidopsis thaliana mutant plants affected in DNA methylation and post-transcriptional gene silencing. We detect a high load of TE-derived eccDNA with truncated and chimeric forms. On the genomic side, on top of truncated and full length TE neo-insertions, we detect complex structural variations (SVs) notably at a disease resistance cluster being a natural hotspot of SV. Finally, we serendipitously identify large tandem duplications in hypomethylated plants, suggesting that SVs could have been overlooked in epigenetic mutants. We propose that a high eccDNA load may alter DNA repair pathways leading to genome instability and the accumulation of SVs, at least in plants.