Clusters of resistance genes in plants evolve by divergent selection and a birth-and-death process

Phylogenetic, Structural, and Evolutionary Insights into Pepper NBS-LRR Resistance Genes

The comprehensive analysis of NBS-LRR resistance genes in the pepper (Capsicum annuum L.) genome reveals their structural diversity, evolutionary history, and functional importance in plant immunity. A total of 252 NBS-LRR genes were identified, distributed unevenly across all chromosomes, with 54% forming 47 gene clusters. These clusters, driven by tandem duplications and genomic rearrangements, underscore the dynamic evolution of resistance genes. Phylogenetic analysis demonstrated the dominance of the nTNL subfamily over the TNL subfamily, reflecting lineage-specific adaptations and evolutionary pressures. Structural analyses identified six conserved motifs (P-loop, RNBS-A, kinase-2, RNBS-B, RNBS-C, and GLPL) essential for ATP/GTP binding and resistance signaling. Subfamily-specific differences in motif composition and sequence similarity highlight their functional divergence and specialization. Comparative analyses across species further revealed a greater prevalence of nTNL genes in angiosperms, with significant losses of TNL genes in monocots. This study enhances our understanding of the evolution and diversification of plant-resistance genes and provides a foundation for developing disease-resistant crops through targeted breeding strategies.

Patterns of presence-absence variation of NLRs across populations of Solanum chilense are clade-dependent and mainly shaped by past demographic history

Understanding the evolution of pathogen resistance genes (nucleotide-binding site-leucine-rich repeats, NLRs) within a species requires a comprehensive examination of factors that affect gene loss and gain. We present a new reference genome of Solanum chilense, which leads to an increased number and more accurate annotation of NLRs. Using a target capture approach, we quantify the presence-absence variation (PAV) of NLR loci across 20 populations from different habitats. We build a rigorous pipeline to validate the identification of PAV of NLRs and then show that PAV is larger within populations than between populations, suggesting that maintenance of NLR diversity is linked to population dynamics. The amount of PAV appears not to be correlated with the NLR presence in gene clusters in the genome, but rather with the past demographic history of the species, with loss of NLRs in diverging (smaller) populations at the distribution edges. Finally, using a redundancy analysis, we find limited evidence of PAV being linked to environmental gradients. Our results suggest that random processes (genetic drift and demography) and weak positive selection for local adaptation shape the evolution of NLRs at the single nucleotide polymorphism and PAV levels in an outcrossing plant with high nucleotide diversity.

Characterization of acetolactate synthase genes and resistance mechanisms of multiple herbicide resistant Lolium multiflorum

Combining imidazolinone-tolerant wheat with imazamox presents an effective solution to combat weed resistance. However, Lolium multiflorum, a troublesome resistant weed infesting wheat fields, may have developed resistance to imazamox, and the potential resistance mechanisms are intriguing. In this study, we explored the susceptibility of L. multiflorum to imazamox and investigated the resistance mechanisms, including the contribution of the target enzyme acetolactate synthase (ALS) to resistance and the presence of non-target-site resistance (NTSR). Eight L. multiflorum populations suspected of being resistant to imazamox were collected, and six populations exhibited resistance, ranging from 2.45-fold to 16.32-fold. The LmALS1 gene from susceptible population D3 plants and multiple copies of the LmALS gene (LmALS1, LmALS2, LmALS2α, LmALS3, LmALS3α, LmALS3β) from resistant populations D5 and D8 plants were separately amplified. Two mutations (Pro/Gln197 to Thr, Trp574 to Leu) were found in LmALS1 in the resistant populations. Compared to D3, LmALS1 was overexpressed in D5 but not in D8. The presence of LmALS1 mutants (LmALS1-Thr197 and LmALS1- Leu574), along with LmALS2, LmALS3, and their subunits, contribute to the resistance phenotype by increasing bonding energies, weakening hydrogen bonds, or decreasing protein binding pocket volumes and surface area. Additionally, D5 and D8 populations exhibited multiple resistance (>40-fold) to three other ALS inhibitors: pyroxsulam, flucarbazone-sodium, and mesosulfuron-methyl. Pre-treatment with malathion and 4-chloro-7-nitrobenzoxadiazole (cytochrome P450 monooxygenase and glutathione S-transferase inhibitors respectively) reversed the resistance of the D8 population and partially reversed the resistance of the D5 population to imazamox. This study characterizes ALS genes and extends our knowledge into the ALS resistance mechanisms involved in L. multiflorum. It also deepens our understanding of the complex diversification resistance mechanisms, thereby facilitating advances in weed resistance management strategies in wheat fields.

Decoding resilience: ecology, regulation, and evolution of biosynthetic gene clusters

Secondary metabolism is crucial for plant survival and can generate chemistry with nutritional, therapeutic, and industrial value. Biosynthetic genes of selected secondary metabolites cluster within localised chromosomal regions. The arrangement of these biosynthetic gene clusters (BGCs) challenges the long-held model of random gene order in eukaryotes, raising questions about their regulation, ecological significance, and evolution. In this review, we address these questions by exploring the contribution of BGCs to ecologically relevant plant-biotic interactions, while also evaluating the molecular-(epi)genetic mechanisms controlling their coordinated stress- and tissue-specific expression. Based on evidence that BGCs have distinct chromatin signatures and are enriched with transposable elements (TEs), we integrate emerging hypotheses into an updated evolutionary model emphasising how stress-induced epigenetic processes have shaped BGC formation.

Genome-Wide Association Study for Resistance to Phytophthora sojae in Soybean [Glycine max (L.) Merr.]

Phytophthora sojae (Kauffman and Gerdemann) is an oomycete pathogen that threatens soybean (Glycine max L.) production worldwide. The development of soybean cultivars with resistance to this pathogen is of paramount importance for the sustainable management of the disease. The objective of this study was to identify genomic regions associated with resistance to P. sojae isolate 40468 through genome-wide association analyses of 983 soybean germplasms. To elucidate the genetic basis of resistance, three statistical models were employed: the compressed mixed linear model (CMLM), Bayesian-information and linkage disequilibrium iteratively nested keyway (BLINK), and fixed and random model circulating probability unification (FarmCPU). The three models consistently identified a genomic region (3.8-5.3 Mbp) on chromosome 3, which has been previously identified as an Rps cluster. A total of 18 single nucleotide polymorphisms demonstrated high statistical significance across all three models, which were distributed in eight linkage disequilibrium (LD) blocks within the aforementioned interval. Of the eight, LD3-2 exhibited the discernible segregation of phenotypic reactions by haplotype. Specifically, over 93% of accessions with haplotypes LD3-2-F or LD3-2-G displayed resistance, whereas over 91% with LD3-2-A, LD3-2-C, or LD3-2-D exhibited susceptibility. Furthermore, the BLINK and FarmCPU models identified new genomic variations significantly associated with the resistance on several other chromosomes, indicating that the resistance observed in this panel was due to the presence of different alleles of multiple Rps genes. These findings underscore the necessity for robust statistical models to accurately detect true marker-trait associations and provide valuable insights into soybean genetics and breeding.

Structural variation in the pangenome of wild and domesticated barley

Pangenomes are collections of annotated genome sequences of multiple individuals of a species1. The structural variants uncovered by these datasets are a major asset to genetic analysis in crop plants2. Here we report a pangenome of barley comprising long-read sequence assemblies of 76 wild and domesticated genomes and short-read sequence data of 1,315 genotypes. An expanded catalogue of sequence variation in the crop includes structurally complex loci that are rich in gene copy number variation. To demonstrate the utility of the pangenome, we focus on four loci involved in disease resistance, plant architecture, nutrient release and trichome development. Novel allelic variation at a powdery mildew resistance locus and population-specific copy number gains in a regulator of vegetative branching were found. Expansion of a family of starch-cleaving enzymes in elite malting barleys was linked to shifts in enzymatic activity in micro-malting trials. Deletion of an enhancer motif is likely to change the developmental trajectory of the hairy appendages on barley grains. Our findings indicate that allelic diversity at structurally complex loci may have helped crop plants to adapt to new selective regimes in agricultural ecosystems.

QTL-seq and QTL mapping identify a new locus for Cercospora leaf spot (Cercospora canescens) resistance in mungbean (Vigna radiata) and a cluster of Receptor-like protein 12 (RLP12) genes as candidate genes for the resistance

KEY MESSAGE

QTL-seq, linkage mapping, and whole-genome resequencing revealed a new locus (qCLS5.1) controlling Cercospora canescens resistance in mungbean and Receptor-like protein 12 (RLP12) genes as candidate genes for the resistance. Cercospora leaf spot (CLS) disease, caused by Cercospora canescens, is a common disease of mungbean (Vigna radiata). In this study, the genetics of CLS resistance was investigated in a new source of resistance (accession V2817) and the resistance was finely mapped to identify candidate genes. F2 and F2:3 populations of the cross V1197 (susceptible) × V2718 and a BC1F1 population of the cross V1197 × (V1197 × V2817) were used in this study. Segregation analysis suggested that the resistance is controlled by a single dominant gene. QTL-seq using F2 individuals revealed that a single QTL (designated qCLS5.1) on chromosome 5 controlled the resistance. The qCLS5.1 was confirmed in the F2:3 and BC1F1 populations by QTL analysis. Fine mapping using 978 F2 individuals localized qCLS5.1 to a 48.94 Kb region containing three tandemly duplicated Receptor-like protein 12 (RLP12) genes. Whole-genome resequencing and alignment of V1197 and V2817 revealed polymorphisms causing amino acid changes and premature stop codons in the three RLP12 genes. Collectively, these results show that qCLS5.1 is a new locus for CLS resistance in mungbean, and a cluster of RLP12 genes are candidate genes for the resistance. The new locus qCLS5.1 will be useful for molecular breeding of durable CLS-resistant mungbean cultivars.

Majority of the Highly Variable NLRs in Maize Share Genomic Location and Contain Additional Target-Binding Domains

Nucleotide-binding, leucine-rich repeat (LRR) proteins (NLRs) are a major class of immune receptors in plants. NLRs include both conserved and rapidly evolving members; however, their evolutionary trajectory in crops remains understudied. Availability of crop pan-genomes enables analysis of the recent events in the evolution of this highly complex gene family within domesticated species. Here, we investigated the NLR complement of 26 nested association mapping (NAM) founder lines of maize. We found that maize has just four main subfamilies containing rapidly evolving highly variable NLR (hvNLR) receptors. Curiously, three of these phylogenetically distinct hvNLR lineages are located in adjacent clusters on chromosome 10. Members of the same hvNLR clade show variable expression and methylation across lines and tissues, which is consistent with their rapid evolution. By combining sequence diversity analysis and AlphaFold2 computational structure prediction, we predicted ligand-binding sites in the hvNLRs. We also observed novel insertion domains in the LRR regions of two hvNLR subfamilies that likely contribute to target recognition. To make this analysis accessible, we created NLRCladeFinder, a Google Colaboratory notebook, that accepts any newly identified NLR sequence, places it in the evolutionary context of the maize pan-NLRome, and provides an updated clade alignment, phylogenetic tree, and sequence diversity information for the gene of interest. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Genetic dissection of resistance to Phytophthora sojae using genome-wide association and linkage analysis in soybean [Glycine max (L.) Merr.]

KEY MESSAGE

Two novel and one known genomic regions associated with R-gene resistance to Phytophthora sojae were identified by genome-wide association analysis and linkage analysis in soybean. Phytophthora root and stem rot (PRR) caused by Phytophthora sojae is a severe disease that causes substantial economic losses in soybean [Glycine max (L.) Merr.]. The primary approach for successful disease management of PRR is using R-gene-mediated resistance. Based on the phenotypic evaluation of 376 cultivated soybean accessions for the R-gene type resistance to P. sojae (isolate 2457), a genome-wide association analysis identified two regions on chromosomes 3 and 8. The most significant genomic region (20.7-21.3 Mbp) on chromosome 8 was a novel resistance locus where no Rps gene was previously reported. Instead, multiple copies of the UDP-glycosyltransferase superfamily protein-coding gene, associated with disease resistance, were annotated in this new locus. Another genomic region on chromosome 3 was a well-known Rps cluster. Using the Daepung × Ilpumgeomjeong RIL population, a linkage analysis confirmed these two resistance loci and identified a resistance locus on chromosome 2. A unique feature of the resistance in Ilpumgeomjeong was discovered when phenotypic distribution was projected upon eight groups of RILs carrying different combinations of resistance alleles for the three loci. Interestingly, the seven groups carrying at least one resistance allele statistically differed from the other with none, regardless of the number of resistance alleles. This suggests that the respective three different resistance genes can confer resistance to P. sojae isolate 2457. Deployment of the three regions via marker-assisted selection will facilitate effectively improving resistance to particular P. sojae isolates in soybean breeding programs.

Phylogenomic insights into the diversity and evolution of RPW8-NLRs and their partners in plants

Plants use nucleotide-binding leucine-rich repeat receptors (NLRs) to sense pathogen effectors, initiating effector-triggered immunity (ETI). NLRs containing RESISTANCE TO POWDERY MILDEW 8 domain (RNLs) function as "helper" NLRs in flowering plants and support the immune responses mediated by "sensor" NLRs in cooperation with lipase-EP domain fused proteins (EP proteins). Despite their crucial roles in ETI, much remains unclear about the evolutionary trajectories of RNLs and their functional partners EP proteins. Here, we perform phylogenomic analyses of RNLs in 90 plants, covering the major diversity of plants, and identify the presence of RNLs in land plants and green algae, expanding the distribution of RNLs. We uncover a neglected major RNL group in gymnosperms, besides the canonical major group with NRG1s and ADR1s, and observe a drastic increase in RNL repertoire size in conifers. Phylogenetic analyses indicate that RNLs originated multiple times through domain shuffling, and the evolution of RNLs underwent a birth-and-death process. Moreover, we trace the origin of EP proteins back to the last common ancestor of vascular plants. We find that both RNLs and EP proteins evolve mainly under negative selection, revealing strong constraints on their function. Concerted losses and positive correlation in copy number are observed between RNL and EP sublineages, suggesting their cooperation in function. Together, our findings provide insights into the origin and evolution of plant helper NLRs, with implications for predicting novel innate immune signaling modules.

Diverse epistatic effects in barley-powdery mildew interactions localize to host chromosome hotspots

Barley Mildew locus a (Mla) encodes a multi-allelic series of nucleotide-binding leucine-rich repeat (NLR) receptors that specify recognition to diverse cereal diseases. We exploited time-course transcriptome dynamics of barley and derived immune mutants infected with the powdery mildew fungus, Blumeria hordei (Bh), to infer gene effects governed by Mla6 and two other loci significant to disease development, Blufensin1 (Bln1), and Required for Mla6 resistance3 (rar3 = Sgt1 ΔKL308-309 ). Interactions of Mla6 and Bln1 resulted in diverse epistatic effects on the Bh-induced barley transcriptome, differential immunity to Pseudomonas syringae expressing the effector protease AvrPphB, and reaction to Bh. From a total of 468 barley NLRs, 115 were grouped under different gene effect models; genes classified under these models localized to host chromosome hotspots. The corresponding Bh infection transcriptome was classified into nine co-expressed modules, linking differential expression with pathogen structures, signifying that disease is regulated by an inter-organismal network that diversifies the response.

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.

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.

High-resolution bulked segregant analysis enables candidate gene identification for bacterial wilt resistance in Italian ryegrass (Lolium multiflorum Lam.)

Bacterial wilt, caused by Xanthomonas translucens pv. graminis (Xtg), is a serious disease of economically important forage grasses, including Italian ryegrass (Lolium multiflorum Lam.). A major QTL for resistance to Xtg was previously identified, but the precise location as well as the genetic factors underlying the resistance are yet to be determined. To this end, we applied a bulked segregant analysis (BSA) approach, using whole-genome deep sequencing of pools of the most resistant and most susceptible individuals of a large (n = 7484) biparental F2 population segregating for resistance to Xtg. Using chromosome-level genome assemblies as references, we were able to define a ~300 kb region highly associated with resistance on pseudo-chromosome 4. Further investigation of this region revealed multiple genes with a known role in disease resistance, including genes encoding for Pik2-like disease resistance proteins, cysteine-rich kinases, and RGA4- and RGA5-like disease resistance proteins. Investigation of allele frequencies in the pools and comparative genome analysis in the grandparents of the F2 population revealed that some of these genes contain variants with allele frequencies that correspond to the expected heterozygosity in the resistant grandparent. This study emphasizes the efficacy of combining BSA studies in very large populations with whole genome deep sequencing and high-quality genome assemblies to pinpoint regions associated with a binary trait of interest and accurately define a small set of candidate genes. Furthermore, markers identified in this region hold significant potential for marker-assisted breeding strategies to breed resistance to Xtg in Italian ryegrass cultivars more efficiently.

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.

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.

Haplotypes at the sorghum ARG4 and ARG5 NLR loci confer resistance to anthracnose

Sorghum anthracnose caused by the fungus Colletotrichum sublineola (Cs) is a damaging disease of the crop. Here, we describe the identification of ANTHRACNOSE RESISTANCE GENES (ARG4 and ARG5) encoding canonical nucleotide-binding leucine-rich repeat (NLR) receptors. ARG4 and ARG5 are dominant resistance genes identified in the sorghum lines SAP135 and P9830, respectively, that show broad-spectrum resistance to Cs. Independent genetic studies using populations generated by crossing SAP135 and P9830 with TAM428, fine mapping using molecular markers, comparative genomics and gene expression studies determined that ARG4 and ARG5 are resistance genes against Cs strains. Interestingly, ARG4 and ARG5 are both located within clusters of duplicate NLR genes at linked loci separated by ~1 Mb genomic region. SAP135 and P9830 each carry only one of the ARG genes while having the recessive allele at the second locus. Only two copies of the ARG5 candidate genes were present in the resistant P9830 line while five non-functional copies were identified in the susceptible line. The resistant parents and their recombinant inbred lines carrying either ARG4 or ARG5 are resistant to strains Csgl1 and Csgrg suggesting that these genes have overlapping specificities. The role of ARG4 and ARG5 in resistance was validated through sorghum lines carrying independent recessive alleles that show increased susceptibility. ARG4 and ARG5 are located within complex loci displaying interesting haplotype structures and copy number variation that may have resulted from duplication. Overall, the identification of anthracnose resistance genes with unique haplotype stucture provides a foundation for genetic studies and resistance breeding.

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.

Single amino acid change alters specificity of the multi-allelic wheat stem rust resistance locus SR9

Most rust resistance genes thus far isolated from wheat have a very limited number of functional alleles. Here, we report the isolation of most of the alleles at wheat stem rust resistance gene locus SR9. The seven previously reported resistance alleles (Sr9a, Sr9b, Sr9d, Sr9e, Sr9f, Sr9g, and Sr9h) are characterised using a synergistic strategy. Loss-of-function mutants and/or transgenic complementation are used to confirm Sr9b, two haplotypes of Sr9e (Sr9e_h1 and Sr9e_h2), Sr9g, and Sr9h. Each allele encodes a highly related nucleotide-binding site leucine-rich repeat (NB-LRR) type immune receptor, containing an unusual long LRR domain, that confers resistance to a unique spectrum of isolates of the wheat stem rust pathogen. The only SR9 protein effective against stem rust pathogen race TTKSK (Ug99), SR9H, differs from SR9B by a single amino acid. SR9B and SR9G resistance proteins are also distinguished by only a single amino acid. The SR9 allelic series found in the B subgenome are orthologs of wheat stem rust resistance gene Sr21 located in the A subgenome with around 85% identity in protein sequences. Together, our results show that functional diversification of allelic variants at the SR9 locus involves single and multiple amino acid changes that recognize isolates of wheat stem rust.

Chromosome-level genomes of three key Allium crops and their trait evolution

Allium crop breeding remains severely hindered due to the lack of high-quality reference genomes. Here we report high-quality chromosome-level genome assemblies for three key Allium crops (Welsh onion, garlic and onion), which are 11.17 Gb, 15.52 Gb and 15.78 Gb in size with the highest recorded contig N50 of 507.27 Mb, 109.82 Mb and 81.66 Mb, respectively. Beyond revealing the genome evolutionary process of Allium species, our pathogen infection experiments and comparative metabolomic and genomic analyses showed that genes encoding enzymes involved in the metabolic pathway of Allium-specific flavor compounds may have evolved from an ancient uncharacterized plant defense system widely existing in many plant lineages but extensively boosted in alliums. Using in situ hybridization and spatial RNA sequencing, we obtained an overview of cell-type categorization and gene expression changes associated with spongy mesophyll cell expansion during onion bulb formation, thus indicating the functional roles of bulb formation genes.

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.

Linked candidate genes of different functions for white mold resistance in common bean (Phaseolus vulgaris L) are identified by multiple QTL mapping approaches

White mold (WM) is a major disease in common bean (Phaseolus vulgaris L.), and its complex quantitative genetic control limits the development of WM resistant cultivars. WM2.2, one of the nine meta-QTL with a major effect on WM tolerance, explains up to 35% of the phenotypic variation and was previously mapped to a large genomic interval on Pv02. Our objective was to narrow the interval of this QTL using combined approach of classic QTL mapping and QTL-based bulk segregant analysis (BSA), and confirming those results with Khufu de novo QTL-seq. The phenotypic and genotypic data from two RIL populations, 'Raven'/I9365-31 (R31) and 'AN-37'/PS02-029C-20 (Z0726-9), were used to select resistant and susceptible lines to generate subpopulations for bulk DNA sequencing. The QTL physical interval was determined by considering overlapping interval of the identified QTL or peak region in both populations by three independent QTL mapping analyses. Our findings revealed that meta-QTL WM2.2 consists of three regions, WM2.2a (4.27-5.76 Mb; euchromatic), WM 2.2b (12.19 to 17.61 Mb; heterochromatic), and WM2.2c (23.01-25.74 Mb; heterochromatic) found in both populations. Gene models encoding for gibberellin 2-oxidase 8, pentatricopeptide repeat, and heat-shock proteins are the likely candidate genes associated with WM2.2a resistance. A TIR-NBS-LRR class of disease resistance protein (Phvul.002G09200) and LRR domain containing family proteins are potential candidate genes associated with WM2.2b resistance. Nine gene models encoding disease resistance protein [pathogenesis-related thaumatin superfamily protein and disease resistance-responsive (dirigent-like protein) family protein etc] found within the WM2.2c QTL interval are putative candidate genes. WM2.2a region is most likely associated with avoidance mechanisms while WM2.2b and WM2.2c regions trigger physiological resistance based on putative candidate genes.

The NLRomes of Zea mays NAM founder lines and Zea luxurians display presence-absence variation, integrated domain diversity, and mobility

Plant pathogens cause significant crop loss worldwide, and new resistance genes deployed to combat diseases can be overcome quickly. Understanding the existing resistance gene diversity within the germplasm of major crops, such as maize, is crucial for the development of new disease-resistant varieties. We analysed the nucleotide-binding leucine-rich repeat receptors (NLRs) of 26 recently sequenced diverse founder lines from the maize nested association mapping (NAM) population and compared them to the R gene complement present in a wild relative of maize, Zea luxurians. We found that NLRs in both species contain a large diversity of atypical integrated domains, including many domains that have not previously been found in the NLRs of other species. Additionally, the single Z. luxurians genome was found to have greater integrated atypical domain diversity than all 26 NAM founder lines combined, indicating that this species may represent a rich source of novel resistance genes. NLRs were also found to have very high sequence diversity and presence-absence variation among the NAM founder lines, with a large NLR cluster on Chr10 representing a diversity hotspot. Additionally, NLRs were shown to be mobile within maize genomes, with several putative interchromosomal translocations identified.

Investigation of Resistance Genes in Genus Vigna Reveals Highly Variable NLRome in Parallel Domesticated Member Species

Vigna is a unique genus that consist of multiple crop species that are domesticated in parallel fashion between 7-10 thousand years ago. Here we studied the evolution of nucleotide-binding site leucine-rich repeat receptor (NLR) genes across five crop species of genus Vigna. In total identified 286, 350, 234, 250, 108 and 161 NLR genes were from Phaseolous vulgaris, Vigna. unguiculata, Vigna mungo, Vigna radiata, Vigna angularis and Vigna umbellata respectively. Comprehensive phylogenetic and clusterization analysis reveals the presence of seven subgroups of Coiled coil like NLRs (CC-NLR) genes and four distinct lineages of Toll interleukin receptor like NLRs (TIR-NLR). Subgroup CC_G10-NLR shows large scale diversification among Vigna species suggesting genus specific distinct duplication pattern in Vigna species. Mainly birth of new NLR gene families and higher rate of terminal duplication is the major determinants for expansion of NLRome in genus Vigna. Recent expansion of NLRome in V. anguiculata and V. radiata was also observed which might suggest that domestication have supported their duplication of lineage specific NLR genes. In short, large scale difference in the architecture of NLRome were observed in diploid plant species. Our findings allowed us to hypothesized that independent parallel domestication is the major drivers of highly divergent evolution of NLRome in genus Vigna.

Puncta-localized TRAF domain protein TC1b contributes to the autoimmunity of snc1

Immune receptors play important roles in the perception of pathogens and initiation of immune responses in both plants and animals. Intracellular nucleotide-binding domain leucine-rich repeat (NLR)-type receptors constitute a major class of receptors in vascular plants. In the Arabidopsis thaliana mutant suppressor of npr1-1, constitutive 1 (snc1), a gain-of-function mutation in the NLR gene SNC1 leads to SNC1 overaccumulation and constitutive activation of defense responses. From a CRISPR/Cas9-based reverse genetics screen in the snc1 autoimmune background, we identified that mutations in TRAF CANDIDATE 1b (TC1b), a gene encoding a protein with four tumor necrosis factor receptor-associated factor (TRAF) domains, can suppress snc1 phenotypes. TC1b does not appear to be a general immune regulator as it is not required for defense mediated by other tested immune receptors. TC1b also does not physically associate with SNC1, affect SNC1 accumulation, or affect signaling of the downstream helper NLRs represented by ACTIVATED DISEASE RESISTANCE PROTEIN 1-L2 (ADR1-L2), suggesting that TC1b impacts snc1 autoimmunity in a unique way. TC1b can form oligomers and localizes to punctate structures of unknown function. The puncta localization of TC1b strictly requires its coiled-coil (CC) domain, whereas the functionality of TC1b requires the four TRAF domains in addition to the CC. Overall, we uncovered the TRAF domain protein TC1b as a novel positive contributor to plant immunity.

QTL study reveals candidate genes underlying host resistance in a Red Queen model system

Specific interactions of host and parasite genotypes can lead to balancing selection, maintaining genetic diversity within populations. In order to understand the drivers of such specific coevolution, it is necessary to identify the molecular underpinnings of these genotypic interactions. Here, we investigate the genetic basis of resistance in the crustacean host, Daphnia magna, to attachment and subsequent infection by the bacterial parasite, Pasteuria ramosa. We discover a single locus with Mendelian segregation (3:1 ratio) with resistance being dominant, which we call the F locus. We use QTL analysis and fine mapping to localize the F locus to a 28.8-kb region in the host genome, adjacent to a known resistance supergene. We compare the 28.8-kb region in the two QTL parents to identify differences between host genotypes that are resistant versus susceptible to attachment and infection by the parasite. We identify 13 genes in the region, from which we highlight eight biological candidates for the F locus, based on presence/absence polymorphisms and differential gene expression. The top candidates include a fucosyltransferase gene that is only present in one of the two QTL parents, as well as several Cladoceran-specific genes belonging to a large family that is represented in multiple locations of the host genome. Fucosyltransferases have been linked to resistance in previous studies of Daphnia-Pasteuria and other host-parasite systems, suggesting that P. ramosa spore attachment could be mediated by changes in glycan structures on D. magna cuticle proteins. The Cladoceran-specific candidate genes suggest a resistance strategy that relies on gene duplication. Our results add a new locus to a growing genetic model of resistance in the D. magna-P. ramosa system. The identified candidate genes will be used in future functional genetic studies, with the ultimate aim to test for cycles of allele frequencies in natural populations.

Analysis of R Genes Related to Blackcurrant Reversion Virus Resistance in the Comparative Transcriptome of Ribes nigrum cv. Aldoniai

Blackcurrant reversion virus (BRV) is the most destructive mite-transmitted pathogen in blackcurrants. The understanding of the resistance to BRV is limited, hindering and delaying the selection process. To identify the resistance (R) gene for BRV resistance, a gene expression analysis based on de novo blackcurrant cv. Aldoniai comparative transcriptome analysis (mock- and BRV-inoculated samples at 2 and 4 days post-inoculation (dpi)) was performed. In this study, 111 annotated clusters associated with pathogenesis according to conservative R gene domains were identified. In virus-infected samples, only Cluster-12591.33361 showed significant expression at 4 dpi. The expression profiles of this cluster were significantly associated with the presence of BRV particles in plant tissues, making it a putative R gene in the dominant resistance strategy in the BRV-Ribes nigrum interaction. The newly identified gene R.nigrum_R belongs to the CC-NBS-LRR class and has 63.9% identity with RPM1 in Populus spp. This study provides new insights on dominant putative R genes related to resistance to BRV in R. nigrum, which could aid targeted research and genetic improvement in breeding programs of blackcurrants.

Predicting Cloned Disease Resistance Gene Homologs (CDRHs) in Radish, Underutilised Oilseeds, and Wild Brassicaceae Species

Brassicaceae crops, including Brassica, Camelina and Raphanus species, are among the most economically important crops globally; however, their production is affected by several diseases. To predict cloned disease resistance (R) gene homologs (CDRHs), we used the protein sequences of 49 cloned R genes against fungal and bacterial diseases in Brassicaceae species. In this study, using 20 Brassicaceae genomes (17 wild and 3 domesticated species), 3172 resistance gene analogs (RGAs) (2062 nucleotide binding-site leucine-rich repeats (NLRs), 497 receptor-like protein kinases (RLKs) and 613 receptor-like proteins (RLPs)) were identified. CDRH clusters were also observed in Arabis alpina, Camelina sativa and Cardamine hirsuta with assigned chromosomes, consisting of 62 homogeneous (38 NLR, 17 RLK and 7 RLP clusters) and 10 heterogeneous RGA clusters. This study highlights the prevalence of CDRHs in the wild relatives of the Brassicaceae family, which may lay the foundation for rapid identification of functional genes and genomics-assisted breeding to develop improved disease-resistant Brassicaceae crop cultivars.

Multilevel evolution shapes the function of NB-LRR encoding genes in plant innate immunity

A sophisticated innate immune system based on diverse pathogen receptor genes (PRGs) evolved in the history of plant life. To reconstruct the direction and magnitude of evolutionary trajectories of a given gene family, it is critical to detect the ancestral signatures. The rearrangement of functional domains made up the diversification found in PRG repertoires. Structural rearrangement of ancient domains mediated the NB-LRR evolutionary path from an initial set of modular proteins. Events such as domain acquisition, sequence modification and temporary or stable associations are prominent among rapidly evolving innate immune receptors. Over time PRGs are continuously shaped by different forces to find their optimal arrangement along the genome. The immune system is controlled by a robust regulatory system that works at different scales. It is important to understand how the PRG interaction network can be adjusted to meet specific needs. The high plasticity of the innate immune system is based on a sophisticated functional architecture and multi-level control. Due to the complexity of interacting with diverse pathogens, multiple defense lines have been organized into interconnected groups. Genomic architecture, gene expression regulation and functional arrangement of PRGs allow the deployment of an appropriate innate immunity response.

A novel locus conferring resistance to Puccinia hordei maps to the genomic region corresponding to Rph14 on barley chromosome 2HS

Barley leaf rust (BLR), caused by Puccinia hordei, is best controlled through genetic resistance. An efficient resistance breeding program prioritizes the need to identify, characterize, and map new sources of resistance as well as understanding the effectiveness, structure, and function of resistance genes. In this study, three mapping populations were developed by crossing Israelian barley lines "AGG-396," "AGG-397," and "AGG-403" (carrying unknown leaf rust resistance) with a susceptible variety "Gus" to characterize and map resistance. Genetic analysis of phenotypic data from rust testing F₃s with a P. hordei pathotype 5457 P+ revealed monogenic inheritance in all three populations. Targeted genotyping-by-sequencing of the three populations detected marker trait associations in the same genomic region on the short arm of chromosome 2H between 39 and 57 Mb (AGG-396/Gus), 44 and 64 Mb (AGG-397/Gus), and 31 and 58 Mb (AGG-403/Gus), suggesting that the resistance in all three lines is likely conferred by the same locus (tentatively designated RphAGG396). Two Kompetitive allele-specific PCR (KASP) markers, HvGBSv2-902 and HvGBSv2-932, defined a genetic distance of 3.8 cM proximal and 7.1 cM distal to RphAGG396, respectively. To increase the marker density at the RphAGG396 locus, 75 CAPS markers were designed between two flanking markers. Integration of marker data resulted in the identification of two critical recombinants and mapping RphAGG396 between markers- Mloc-28 (40.75 Mb) and Mloc-41 (41.92 Mb) narrowing the physical window to 1.17 Mb based on the Morex v2.0 reference genome assembly. To enhance map resolution, 600 F₂s were genotyped with markers- Mloc-28 and Mloc-41 and nine recombinants were identified, placing the gene at a genetic distance of 0.5 and 0.2 cM between the two markers, respectively. Two annotated NLR (nucleotide-binding domain leucine-rich repeat) genes (r2.2HG0093020 and r2.2HG0093030) were identified as the best candidates for RphAGG396. A closely linked marker was developed for RphAGG396 that can be used for marker-assisted selection.

A catalogue of resistance gene homologs and a chromosome-scale reference sequence support resistance gene mapping in winter wheat

A resistance gene atlas is an integral component of the breeder's arsenal in the fight against evolving pathogens. Thanks to high-throughput sequencing, catalogues of resistance genes can be assembled even in crop species with large and polyploid genomes. Here, we report on capture sequencing and assembly of resistance gene homologs in a diversity panel of 907 winter wheat genotypes comprising ex situ genebank accessions and current elite cultivars. In addition, we use accurate long-read sequencing and chromosome conformation capture sequencing to construct a chromosome-scale genome sequence assembly of cv. Attraktion, an elite variety representative of European winter wheat. We illustrate the value of our resource for breeders and geneticists by (i) comparing the resistance gene complements in plant genetic resources and elite varieties and (ii) conducting genome-wide associations scans (GWAS) for the fungal diseases yellow rust and leaf rust using reference-based and reference-free GWAS approaches. The gene content under GWAS peaks was scrutinized in the assembly of cv. Attraktion.

Genome-based high-resolution mapping of fusarium wilt resistance in sweet basil

Fusarium wilt of basil is a disease of sweet basil (Ocimum basilicum L.) plants caused by the fungus Fusarium oxysporum f. sp. basilici (FOB). Although resistant cultivars were released > 20 years ago, the underlying mechanism and the genes controlling the resistance remain unknown. We used genetic mapping to elucidate FOB resistance in an F₂ population derived from a cross between resistant and susceptible cultivars. We performed genotyping by sequencing of 173 offspring and aligning the data to the sweet basil reference genome. In total, 23,411 polymorphic sites were detected, and a single quantitative trait locus (QTL) for FOB resistance was found. The confidence interval was < 600 kbp, harboring only 60 genes, including a cluster of putative disease-resistance genes. Based on homology to a fusarium resistance protein from wild tomato, we also investigated a candidate resistance gene that encodes a transmembrane leucine-rich repeat - receptor-like kinase - ubiquitin-like protease (LRR-RLK-ULP). Sequence analysis of that gene in the susceptible parent vs. the resistant parent revealed multiple indels, including an insertion of 20 amino acids next to the transmembrane domain, which might alter its functionality. Our findings suggest that this LRR-RLK-ULP might be responsible for FOB resistance in sweet basil and demonstrate the usefulness of the recently sequenced basil genome for QTL mapping and gene mining.

Chromosome-level genome assembly of Mentha longifolia L. reveals gene organization underlying disease resistance and essential oil traits

Mentha longifolia (L.) Huds., a wild, diploid mint species, has been developed as a model for mint genetic and genomic research to aid breeding efforts that target Verticillium wilt disease resistance and essential oil monoterpene composition. Here, we present a near-complete, chromosome-scale mint genome assembly for M. longifolia USDA accession CMEN 585. This new assembly is an update of a previously published genome draft, with dramatic improvements. A total of 42,107 protein-coding genes were annotated and placed on 12 chromosomal scaffolds. One hundred fifty-three genes contained conserved sequence domains consistent with nucleotide binding site-leucine-rich-repeat plant disease resistance genes. Homologs of genes implicated in Verticillium wilt resistance in other plant species were also identified. Multiple paralogs of genes putatively involved in p-menthane monoterpenoid biosynthesis were identified and several cases of gene clustering documented. Heterologous expression of candidate genes, purification of recombinant target proteins, and subsequent enzyme assays allowed us to identify the genes underlying the pathway that leads to the most abundant monoterpenoid volatiles. The bioinformatic and functional analyses presented here are laying the groundwork for using marker-assisted selection in improving disease resistance and essential oil traits in mints.

Genome-Wide Identification, Characterization, and Comparative Analysis of NLR Resistance Genes in Coffea spp

The largest family of disease resistance genes in plants are nucleotide-binding site leucine-rich repeat genes (NLRs). The products of these genes are responsible for recognizing avirulence proteins (Avr) of phytopathogens and triggering specific defense responses. Identifying NLRs in plant genomes with standard gene annotation software is challenging due to their multidomain nature, sequence diversity, and clustered genomic distribution. We present the results of a genome-wide scan and comparative analysis of NLR loci in three coffee species (Coffea canephora, Coffea eugenioides and their interspecific hybrid Coffea arabica). A total of 1311 non-redundant NLR loci were identified in C. arabica, 927 in C. canephora, and 1079 in C. eugenioides, of which 809, 562, and 695 are complete loci, respectively. The NLR-Annotator tool used in this study showed extremely high sensitivities and specificities (over 99%) and increased the detection of putative NLRs in the reference coffee genomes. The NLRs loci in coffee are distributed among all chromosomes and are organized mostly in clusters. The C. arabica genome presented a smaller number of NLR loci when compared to the sum of the parental genomes (C. canephora, and C. eugenioides). There are orthologous NLRs (orthogroups) shared between coffee, tomato, potato, and reference NLRs and those that are shared only among coffee species, which provides clues about the functionality and evolutionary history of these orthogroups. Phylogenetic analysis demonstrated orthologous NLRs shared between C. arabica and the parental genomes and those that were possibly lost. The NLR family members in coffee are subdivided into two main groups: TIR-NLR (TNL) and non-TNL. The non-TNLs seem to represent a repertoire of resistance genes that are important in coffee. These results will support functional studies and contribute to a more precise use of these genes for breeding disease-resistant coffee cultivars.

Evolutionary footprint of plant immunity

There are pieces of evidence from genomic footprints and fossil records indicating that plants have co-evolved with microbes after terrestrialization for more than 407 million years. Therefore, to truly comprehend plant evolution, we need to understand the co-evolutionary process and history between plants and microbes. Recent developments in genomes and transcriptomes of a vast number of plant species as well as microbes have greatly expanded our knowledge of the evolution of the plant immune system. In this review, we summarize recent advances in the co-evolution between plants and microbes with emphasis on the plant side and point out future research needed for understanding plant-microbial co-evolution. Knowledge of the evolution and variation of the plant immune system will better equip us on designing crops with boosted performance in agricultural fields.

Mining of Cloned Disease Resistance Gene Homologs (CDRHs) in Brassica Species and Arabidopsis thaliana

Simple Summary

Developing cultivars with resistance genes (R genes) is an effective strategy to support high yield and quality in Brassica crops. The availability of clone R gene and genomic sequences in Brassica species and Arabidopsis thaliana provide the opportunity to compare genomic regions and survey R genes across genomic databases. In this paper, we aim to identify genes related to cloned genes through sequence identity, providing a repertoire of species-wide related R genes in Brassica crops. The comprehensive list of candidate R genes can be used as a reference for functional analysis.

Abstract

Various diseases severely affect Brassica crops, leading to significant global yield losses and a reduction in crop quality. In this study, we used the complete protein sequences of 49 cloned resistance genes (R genes) that confer resistance to fungal and bacterial diseases known to impact species in the Brassicaceae family. Homology searches were carried out across Brassica napus, B. rapa, B. oleracea, B. nigra, B. juncea, B. carinata and Arabidopsis thaliana genomes. In total, 660 cloned disease R gene homologs (CDRHs) were identified across the seven species, including 431 resistance gene analogs (RGAs) (248 nucleotide binding site-leucine rich repeats (NLRs), 150 receptor-like protein kinases (RLKs) and 33 receptor-like proteins (RLPs)) and 229 non-RGAs. Based on the position and distribution of specific homologs in each of the species, we observed a total of 87 CDRH clusters composed of 36 NLR, 16 RLK and 3 RLP homogeneous clusters and 32 heterogeneous clusters. The CDRHs detected consistently across the seven species are candidates that can be investigated for broad-spectrum resistance, potentially providing resistance to multiple pathogens. The R genes identified in this study provide a novel resource for the future functional analysis and gene cloning of Brassicaceae R genes towards crop improvement.

Evolution of resistance (R) gene specificity

Plant resistance (R) genes are members of large gene families with significant within and between species variation. It has been hypothesised that a variety of processes have shaped R gene evolution and the evolution of R gene specificity. In this review, we illustrate the main mechanisms that generate R gene diversity and provide examples of how they can change R gene specificity. Next, we explain which evolutionary mechanisms are at play and how they determine the fate of new R gene alleles and R genes. Finally, we place this in a larger context by comparing the diversity and evolution of R gene specificity within and between species scales.

Mapping and characterization of the Rx3 gene for resistance to Xanthomonas euvesicatoria pv. euvesicatoria race T1 in tomato

KEY MESSAGE

Rx3 encodes a typical CC-NBS-LRR resistance protein and confers the resistance to Xanthomonas euvesicatoria pv. euvesicatoria race T1 causing bacterial spot in tomato. Bacterial spot caused by at least four species of Xanthomonas is an epidemic disease severely affecting tomato production worldwide. The use of resistant cultivars is an economical and effective approach to control the disease. An unimproved tomato breeding line Hawaii 7988 has been considered as the most reliable source for resistance to X. euvesicatoria pv. euvesicatoria race T1, and the Rx3 locus located at a 4.53-Mb region on chromosome 5 (SL4.0) is the major locus for resistance to race T1 in this line. In the current study, the Rx3 locus was firstly located to a 1.05-Mb region based on comparisons of marker polymorphisms between the susceptible line Ohio 88119 and resistant lines Hawaii 7998, Ohio 9834 and FG02-7530. Using recombinant inbred lines (F5:6, F6:7, and F7:8) derived from a cross between Ohio 88119 and Ohio 9834, the Rx3 locus was finally mapped to a 64.3-kb interval between markers MG-Rx3-4 and MG-Rx3-A6. The Solyc05g053980 gene, designated as Rx3, encoding a coiled-coil nucleotide-binding leucine-rich repeat protein was considered as the candidate for the Rx3 locus. Expression of the gene could be induced by the infection of race T1 strain. Knockout of the Solyc05g053980 gene through CRISPR/Cas9 editing system in the resistant line FG02-7530 decreased resistance to race T1 strain. These results provide a close step for understanding the resistance mechanism to race T1 in Hawaii 7998 and guide tomato breeders accordingly to improve bacterial spot disease resistance in tomato.

The origin and evolution of a plant resistosome

The nucleotide-binding, leucine-rich receptor (NLR) protein HOPZ-ACTIVATED RESISTANCE 1 (ZAR1), an immune receptor, interacts with HOPZ-ETI-DEFICIENT 1 (ZED1)-related kinases (ZRKs) and AVRPPHB SUSCEPTIBLE 1-like proteins to form a pentameric resistosome, triggering immune responses. Here, we show that ZAR1 emerged through gene duplication and that ZRKs were derived from the cell surface immune receptors wall-associated protein kinases (WAKs) through the loss of the extracellular domain before the split of eudicots and monocots during the Jurassic period. Many angiosperm ZAR1 orthologs, but not ZAR1 paralogs, are capable of oligomerization in the presence of AtZRKs and triggering cell death, suggesting that the functional ZAR1 resistosome might have originated during the early evolution of angiosperms. Surprisingly, inter-specific pairing of ZAR1 and AtZRKs sometimes results in the formation of a resistosome in the absence of pathogen stimulation, suggesting within-species compatibility between ZAR1 and ZRKs as a result of co-evolution. Numerous concerted losses of ZAR1 and ZRKs occurred in angiosperms, further supporting the ancient co-evolution between ZAR1 and ZRKs. Our findings provide insights into the origin of new plant immune surveillance networks.

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.

Integrating Phenotypic and Gene Expression Linkage Mapping to Dissect Rust Resistance in Chickling Pea

Rusts are among the most important foliar biotrophic fungal diseases in legumes. Lathyrus cicera crop can be severely damaged by Uromyces pisi, to which partial resistance has been identified. Nevertheless, the underlying genetic basis and molecular mechanisms of this resistance are poorly understood in L. cicera. To prioritise the causative variants controlling partial resistance to rust in L. cicera, a recombinant inbred line (RIL) population, segregating for response to this pathogen, was used to combine the detection of related phenotypic- and expression-quantitative trait loci (pQTLs and eQTLs, respectively). RILs' U. pisi disease severity (DS) was recorded in three independent screenings at seedling (growth chamber) and in one season of exploratory screening at adult plant stage (semi-controlled field conditions). A continuous DS range was observed in both conditions and used for pQTL mapping. Different pQTLs were identified under the growth chamber and semi-controlled field conditions, indicating a distinct genetic basis depending on the plant developmental stage and/or the environment. Additionally, the expression of nine genes related to U. pisi resistance in L. cicera was quantified for each RIL individual and used for eQTL mapping. One cis-eQTL and one trans-eQTL were identified controlling the expression variation of one gene related to rust resistance - a member of glycosyl hydrolase family 17. Integrating phenotyping, gene expression and linkage mapping allowed prioritising four candidate genes relevant for disease-resistance precision breeding involved in adaptation to biotic stress, cellular, and organelle homeostasis, and proteins directly involved in plant defence.

Meta-analysis reveals consensus genomic regions associated with multiple disease resistance in wheat (Triticum aestivum L.)

In wheat, meta-QTLs (MQTLs) and candidate genes (CGs) were identified for multiple disease resistance (MDR). For this purpose, information was collected from 58 studies for mapping QTLs for resistance to one or more of the five diseases. As many as 493 QTLs were available from these studies, which were distributed in five diseases as follows: septoria tritici blotch (STB) 126 QTLs; septoria nodorum blotch (SNB), 103 QTLs; fusarium head blight (FHB), 184 QTLs; karnal bunt (KB), 66 QTLs; and loose smut (LS), 14 QTLs. Of these 493 QTLs, only 291 QTLs could be projected onto a consensus genetic map, giving 63 MQTLs. The CI of the MQTLs ranged from 0.04 to 15.31 cM with an average of 3.09 cM per MQTL. This is a ~ 4.39 fold reduction from the CI of QTLs, which ranged from 0 to 197.6 cM, with a mean of 13.57 cM. Of 63 MQTLs, 60 were anchored to the reference physical map of wheat (the physical interval of these MQTLs ranged from 0.30 to 726.01 Mb with an average of 74.09 Mb). Thirty-eight (38) of these MQTLs were verified using marker-trait associations (MTAs) derived from genome-wide association studies. As many as 874 CGs were also identified which were further investigated for differential expression using data from five transcriptome studies, resulting in 194 differentially expressed candidate genes (DECGs). Among the DECGs, 85 genes had functions previously reported to be associated with disease resistance. These results should prove useful for fine mapping and cloning of MDR genes and marker-assisted breeding.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01282-z.

Identification and characterization of genome-wide resistance gene analogs (RGAs) of durian (Durio zibethinus L.)

BACKGROUND

Durian (Durio zibethinus L.) is a tropical fruit crop which is popular in Southeast Asia but recently gaining popularity in other parts of the world. In this study, we analyzed the resistance gene analogs (RGAs) of durian through mining of the currently available reference genome of its 'Musang King' cultivar (PRJNA400310).

RESULTS

A total of 2586 RGAs were identified in the durian genome consisting of 47 nucleotide binding site proteins (NBS), 158 NBS-leucine rich repeat proteins (NL), 400 coiled-coil NBS-LRR (CNL), 72 toll/interleukin-1 receptor NBS-LRR (TNL), 54 coiled-coil NBS (CN), 10 toll/interleukin-1 receptor NBS (TN), 19 toll/interleukin-1 receptor with unknown domain (TX), 246 receptor-like proteins (RLP), 1,377 receptor-like kinases (RLK), 185 TM-CC, and 18 other NBS-containing proteins with other domains. These RGAs were functionally annotated and characterized via gene ontology (GO) analysis. Among the RGAs with the highest copies in durian genome include the putative disease resistance RPP13-like protein 1, disease resistance protein At4g27190, disease resistance protein RPS6, Probable disease resistance protein At4g27220, and putative disease resistance protein RGA3, while 35 RGAs were found to be novel. Phylogenetic analyses revealed that the genome-wide RGAs were broadly clustered into four major clades based on their domain classification.

CONCLUSION

To our knowledge, this is the most comprehensive analysis of durian RGAs which provides a valuable resource for genetic, agronomic, and other biological research of this important tropical fruit crop.

Fitness benefits play a vital role in the retention of the Pi-ta susceptible alleles

In plants, large numbers of R genes, which segregate as loci with alternative alleles conferring different levels of disease resistance to pathogens, have been maintained over a long period of evolution. The reason why hosts harbor susceptible alleles in view of their null contribution to resistance is unclear. In rice, a single copy gene, Pi-ta, segregates for 2 expressed clades of alleles, 1 resistant and the other susceptible. We simulated loss-of-function of the Pi-ta susceptible allele using the CRISPR/Cas9 system to detect subsequent fitness changes and obtained insights into fitness effects related to the retention of the Pi-ta susceptible allele. Our creation of an artificial knockout of the Pi-ta susceptible allele suffered fitness-related trait declines of up to 49% in terms of filled grain yield upon the loss of Pi-ta function. The Pi-ta susceptible alleles might serve as an off-switch to downstream immune signaling, thus contributing to the fine-tuning of plant defense responses. The results demonstrated that the susceptible Pi-ta alleles should have evolved pleiotropic functions, facilitating their retention in populations. As Pi-ta is a single copy gene with no paralogs in the genome, its function cannot be compensated by an alternative gene; whereas most other R genes form gene clusters by tandem duplications, and the function could be compensated by paralogs with high sequence similarity. This attempt to evaluate the fitness effects of the R gene in crops indicates that not all disease resistance genes incur fitness costs, which also provides a plausible explanation for how host genomes can tolerate the possible genetic load associated with a vast repertoire of R genes.

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.

A cluster of atypical resistance genes in soybean confers broad-spectrum antiviral activity

Soybean mosaic virus (SMV) is a severe soybean (Glycine max) pathogen. Here we characterize a soybean SMV resistance cluster (SRC) that comprises five resistance (R) genes. SRC1 encodes a Toll/interleukin-1 receptor and nucleotide-binding site (TIR-NBS [TN]) protein, SRC4 and SRC6 encode TIR proteins with a short EF-hand domain, while SRC7 and SRC8 encode TNX proteins with a noncanonical basic secretory protein (BSP) domain at their C-termini. We mainly studied SRC7, which contains a noncanonical BSP domain and gave full resistance to SMV. SRC7 possessed broad-spectrum antiviral activity toward several plant viruses including SMV, plum pox virus, potato virus Y, and tobacco mosaic virus. The TIR domain alone was both necessary and sufficient for SRC7 immune signaling, while the NBS domain enhanced its activity. Nuclear oligomerization via the interactions of both TIR and NBS domains was essential for SRC7 function. SRC7 expression was transcriptionally inducible by SMV infection and salicylic acid (SA) treatment, and SA was required for SRC7 triggered virus resistance. SRC7 expression was posttranscriptionally regulated by miR1510a and miR2109, and the SRC7-miR1510a/miR2109 regulatory network appeared to contribute to SMV-soybean interactions in both resistant and susceptible soybean cultivars. In summary, we report a soybean R gene cluster centered by SRC7 that is regulated at both transcriptional and posttranscriptional levels, possesses a yet uncharacterized BSP domain, and has broad-spectrum antiviral activities. The SRC cluster is special as it harbors several functional R genes encoding atypical TIR-NBS-LRR (TNL) type R proteins, highlighting its importance in SMV-soybean interaction and plant immunity.

Epigenetics: a catalyst of plant immunity against pathogens

The plant immune system protects against pests and diseases. The recognition of stress-related molecular patterns triggers localised immune responses, which are often followed by longer-lasting systemic priming and/or up-regulation of defences. In some cases, this induced resistance (IR) can be transmitted to following generations. Such transgenerational IR is gradually reversed in the absence of stress at a rate that is proportional to the severity of disease experienced in previous generations. This review outlines the mechanisms by which epigenetic responses to pathogen infection shape the plant immune system across expanding time scales. We review the cis- and trans-acting mechanisms by which stress-inducible epigenetic changes at transposable elements (TEs) regulate genome-wide defence gene expression and draw particular attention to one regulatory model that is supported by recent evidence about the function of AGO1 and H2A.Z in transcriptional control of defence genes. Additionally, we explore how stress-induced mobilisation of epigenetically controlled TEs acts as a catalyst of Darwinian evolution by generating (epi)genetic diversity at environmentally responsive genes. This raises questions about the long-term evolutionary consequences of stress-induced diversification of the plant immune system in relation to the long-held dichotomy between Darwinian and Lamarckian evolution.

The barley immune receptor Mla recognizes multiple pathogens and contributes to host range dynamics

Crop losses caused by plant pathogens are a primary threat to stable food production. Stripe rust (Puccinia striiformis) is a fungal pathogen of cereal crops that causes significant, persistent yield loss. Stripe rust exhibits host species specificity, with lineages that have adapted to infect wheat and barley. While wheat stripe rust and barley stripe rust are commonly restricted to their corresponding hosts, the genes underlying this host specificity remain unknown. Here, we show that three resistance genes, Rps6, Rps7, and Rps8, contribute to immunity in barley to wheat stripe rust. Rps7 cosegregates with barley powdery mildew resistance at the Mla locus. Using transgenic complementation of different Mla alleles, we confirm allele-specific recognition of wheat stripe rust by Mla. Our results show that major resistance genes contribute to the host species specificity of wheat stripe rust on barley and that a shared genetic architecture underlies resistance to the adapted pathogen barley powdery mildew and non-adapted pathogen wheat stripe rust.

Genomic regions associated with resistance to Fusarium wilt in castor identified through linkage and association mapping approaches

Fusarium wilt, caused by Fusarium oxysporum f. sp. ricini, is the most destructive disease in castor. Host plant resistance is the best strategy for the management of wilt. Identification of molecular markers linked to wilt resistance will enhance the efficiency and effectiveness of breeding for wilt resistance. In the present study, genomic regions linked to wilt resistance were mapped using a bi-parental population of 185 F₆-RILs and a genetically diverse panel of 300 germplasm accessions. Quantitative trait loci (QTL) analysis performed using a linkage map consisting of 1090 SNP markers identified a major QTL on chromosome 7 with an LOD score of 18.7, which explained 44% of the phenotypic variance. The association mapping performed using genotypic data from 3465 SNP loci revealed 69 significant associations (p < 1 × 10^-4) for wilt resistance. The phenotypic variance explained by the individual SNPs ranged from 0.063 to 0.210. The QTL detected in the bi-parental mapping population was not identified in the association analysis. Thus, the results of this study indicate the possibility of vast gene diversity for Fusarium wilt resistance in castor.