Impacts of resistance gene genetics, function, and evolution on a durable future

Genome-wide characterization of the NBLRR gene family provides evolutionary and functional insights into blast resistance in pearl millet (Cenchrus americanus (L.) Morrone)

MAIN CONCLUSION

The investigation is the first report on genome-wide identification and characterization of NBLRR genes in pearl millet. We have shown the role of gene loss and purifying selection in the divergence of NBLRRs in Poaceae lineage and candidate CaNBLRR genes for resistance to Magnaporthe grisea infection. Plants have evolved multiple integral mechanisms to counteract the pathogens' infection, among which plant immunity through NBLRR (nucleotide-binding site, leucine-rich repeat) genes is at the forefront. The genome-wide mining in pearl millet (Cenchrus americanus (L.) Morrone) revealed 146 CaNBLRRs. The variation in the branch length of NBLRRs showed the dynamic nature of NBLRRs in response to evolving pathogen races. The orthology of NBLRRs showed a predominance of many-to-one orthologs, indicating the divergence of NBLRRs in the pearl millet lineage mainly through gene loss events followed by gene gain through single-copy duplications. Further, the purifying selection (Ka/Ks < 1) shaped the expansion of NBLRRs within the lineage of pear millet and other members of Poaceae. Presence of cis-acting elements, viz. TCA element, G-box, MYB, SARE, ABRE and conserved motifs annotated with P-loop, kinase 2, RNBS-A, RNBS-D, GLPL, MHD, Rx-CC and LRR suggests their putative role in disease resistance and stress regulation. The qRT-PCR analysis in pearl millet lines showing contrasting responses to Magnaporthe grisea infection identified CaNBLRR20, CaNBLRR33, CaNBLRR46 CaNBLRR51, CaNBLRR78 and CaNBLRR146 as putative candidates. Molecular docking showed the involvement of three and two amino acid residues of LRR domains forming hydrogen bonds (histidine, arginine and threonine) and salt bridges (arginine and lysine) with effectors. Whereas 14 and 20 amino acid residues of CaNBLRR78 and CaNBLRR20 showed hydrophobic interactions with 11 and 9 amino acid residues of effectors, Mg.00g064570.m01 and Mg.00g006570.m01, respectively. The present investigation gives a comprehensive overview of CaNBLRRs and paves the foundation for their utility in pearl millet resistance breeding through understanding of host-pathogen interactions.

Development of an NLR-ID Toolkit and Identification of Novel Disease-Resistance Genes in Soybean

The recognition of pathogen effectors through the nucleotide-binding leucine-rich repeat receptor (NLR) family is an important component of plant immunity. In addition to typical domains such as TIR, CC, NBS, and LRR, NLR proteins also contain some atypical integrated domains (IDs), the roles of which are rarely investigated. Here, we carefully screened the soybean (Glycine max) genome and identified the IDs that appeared in the soybean TNL-like proteins. Our results show that multiple IDs (36) are widely present in soybean TNL-like proteins. A total of 27 Gm-TNL-ID genes (soybean TNL-like gene encoding ID) were cloned and their antiviral activity towards the soybean mosaic virus (SMV)/tobacco mosaic virus (TMV) was verified. Two resistance (R) genes, SRA2 (SMV resistance gene contains AAA_22 domain) and SRZ4 (SMV resistance gene contains zf-RVT domain), were identified to possess broad-spectrum resistance characteristics towards six viruses including SMV, TMV, plum pox virus (PPV), cabbage leaf curl virus (CaLCuV), barley stripe mosaic virus (BSMV), and tobacco rattle virus (TRV). The effects of Gm-TNL-IDX (the domain of the Gm-TNL-ID gene after the TN domain) on the antiviral activity of a R protein SRC7TN (we previously reported the TN domain of the soybean broad-spectrum resistance gene SRC7) were validated, and most of Gm-TNL-IDX inhibits antiviral activity mediated by SRC7TN, possibly through intramolecular interactions. Yeast-two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) assays showed that seven Gm-TNL-IDX interacted with SMV-component proteins. Truncation analysis on a broad-spectrum antiviral protein SRZ4 indicated that SRZ4TIR is sufficient to mediate antiviral activity against SMV. Soybean cDNA library screening on SRZ4 identified 48 interacting proteins. In summary, our results indicate that the integration of IDs in soybean is widespread and frequent. The NLR-ID toolkit we provide is expected to be valuable for elucidating the functions of atypical NLR proteins in the plant immune system and lay the foundation for the development of engineering NLR for plant-disease control in the future.

Fine-mapping and evolutionary history of R-BPMV, a dominant resistance gene to Bean pod mottle virus in Phaseolus vulgaris L

KEY MESSAGE

R-BPMV is located within a recently expanded TNL cluster in the Phaseolus genus with suppressed recombination and known for resistance to multiple pathogens including potyviruses controlled by the I gene. Bean pod mottle virus (BPMV) is a comovirus that infects common bean and legumes in general. BPMV is distributed throughout the world and is a major threat on soybean, a closely related species of common bean. In common bean, BAT93 was reported to carry the R-BPMV resistance gene conferring resistance to BPMV and linked with the I resistance gene. To fine map R-BPMV, 182 recombinant inbred lines (RILs) derived from the cross BAT93 × JaloEEP558 were genotyped with polymerase chain reaction (PCR)-based markers developed using genome assemblies from G19833 and BAT93, as well as BAT93 BAC clone sequences. Analysis of RILs carrying key recombination events positioned R-BPMV to a target region containing at least 16 TIR-NB-LRR (TNL) sequences in BAT93. Because the I cluster presents a suppression of recombination and a large number of repeated sequences, none of the 16 TNLs could be excluded as R-BPMV candidate gene. The evolutionary history of the TNLs for the I cluster were reconstructed using microsynteny and phylogenetic analyses within the legume family. A single I TNL was present in Medicago truncatula and lost in soybean, mirroring the absence of complete BPMV resistance in soybean. Amplification of TNLs in the I cluster predates the divergence of the Phaseolus species, in agreement with the emergence of R-BPMV before the separation of the common bean wild centers of diversity. This analysis provides PCR-based markers useful in marker-assisted selection (MAS) and laid the foundation for cloning of R-BPMV resistance gene in order to transfer the resistance into soybean.

Exogenous dsRNA-Induced Silencing of the Phytophthora infestans Elicitin Genes inf1 and inf4 Suppresses Its Pathogenicity on Potato Plants

Phytophthora infestans, an Oomycete pathogen, has a devastating impact on potato agriculture, leading to the extensive use of chemical fungicides to prevent its outbreaks. Spraying double-stranded RNAs to suppress specific genes of the pathogen via the RNA interference (RNAi) pathway may provide an environmentally friendly alternative to chemicals. However, this novel approach will require various target genes and application strategies to be tested. Using the L4440 backbone, we have designed two plasmids to express dsRNA targeting inf1 and inf4 genes of P. infestans that are known to contribute to the disease development at different stages. The dsRNA produced by the bacteria was tested on potato explants and demonstrated a statistically significant reduction in lesions five days after inoculation compared to water treatment. The study results allow us to consider our approach to be promising for potato late blight control.

Induction and suppression of gene silencing in plants by nonviral microbes

Gene silencing is a conserved mechanism in eukaryotes that dynamically regulates gene expression. In plants, gene silencing is critical for development and for maintenance of genome integrity. Additionally, it is a critical component of antiviral defence in plants, nematodes, insects, and fungi. To overcome gene silencing, viruses encode effectors that suppress gene silencing. A growing body of evidence shows that gene silencing and suppression of silencing are also used by plants during their interaction with nonviral pathogens such as fungi, oomycetes, and bacteria. Plant-pathogen interactions involve trans-kingdom movement of small RNAs into the pathogens to alter the function of genes required for their development and virulence. In turn, plant-associated pathogenic and nonpathogenic microbes also produce small RNAs that move trans-kingdom into host plants to disrupt pathogen defence through silencing of plant genes. The mechanisms by which these small RNAs move from the microbe to the plant remain poorly understood. In this review, we examine the roles of trans-kingdom small RNAs and silencing suppressors produced by nonviral microbes in inducing and suppressing gene silencing in plants. The emerging model is that gene silencing and suppression of silencing play critical roles in the interactions between plants and their associated nonviral microbes.

EffectorO: Motif-Independent Prediction of Effectors in Oomycete Genomes Using Machine Learning and Lineage Specificity

Oomycete plant pathogens cause a wide variety of diseases, including late blight of potato, sudden oak death, and downy mildews of plants. These pathogens are major contributors to loss in numerous food crops. Oomycetes secrete effector proteins to manipulate their hosts to the advantage of the pathogen. Plants have evolved to recognize effectors, resulting in an evolutionary cycle of defense and counter-defense in plant-microbe interactions. This selective pressure results in highly diverse effector sequences that can be difficult to computationally identify using only sequence similarity. We developed a novel effector prediction tool, EffectorO, that uses two complementary approaches to predict effectors in oomycete pathogen genomes: i) a machine learning-based pipeline that predicts effector probability based on the biochemical properties of the N-terminal amino-acid sequence of a protein and ii) a pipeline based on lineage specificity to find proteins that are unique to one species or genus, a sign of evolutionary divergence due to adaptation to the host. We tested EffectorO on Bremia lactucae, which causes lettuce downy mildew, and Phytophthora infestans, which causes late blight of potato and tomato, and predicted many novel effector candidates while recovering the majority of known effector candidates. EffectorO will be useful for discovering novel families of oomycete effectors without relying on sequence similarity to known effectors. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

The grape powdery mildew resistance loci Ren2, Ren3, Ren4D, Ren4U, Run1, Run1.2b, Run2.1, and Run2.2 activate different transcriptional responses to Erysiphe necator

Multiple grape powdery mildew (PM) genetic resistance (R) loci have been found in wild grape species. Little is known about the defense responses associated with each R locus. In this study, we compare the defense mechanisms associated with PM resistance in interspecific crosses segregating for a single R locus from Muscadinia rotundifolia (Run1, Run1.2b, Run2.1, Run2.2), Vitis cinerea (Ren2), V. romanetii (Ren4D and Ren4U), and the interspecific hybrid Villard blanc (Ren3). By combining optical microscopy, visual scoring, and biomass estimation, we show that the eight R loci confer resistance by limiting infection at different stages. We assessed the defense mechanisms triggered in response to PM at 1 and 5 days post-inoculation (dpi) via RNA sequencing. To account for the genetic differences between species, we developed for each accession a diploid synthetic reference transcriptome by incorporating into the PN40024 reference homozygous and heterozygous sequence variants and de novo assembled transcripts. Most of the R loci exhibited a higher number of differentially expressed genes (DEGs) associated with PM resistance at 1 dpi compared to 5 dpi, suggesting that PM resistance is mostly associated with an early transcriptional reprogramming. Comparison of the PM resistance-associated DEGs showed a limited overlap between pairs of R loci, and nearly half of the DEGs were specific to a single R locus. The largest overlap of PM resistance-associated DEGs was found between Ren3 ⁺, Ren4D ⁺, and Ren4U ⁺ genotypes at 1 dpi, and between Ren4U ⁺ and Run1 ⁺ accessions at 5 dpi. The Ren3 ⁺, Ren4D ⁺, and Ren4U ⁺ were also found to have the highest number of R locus-specific DEGs in response to PM. Both shared and R locus-specific DEGs included genes from different defense-related categories, indicating that the presence of E. necator triggered distinct transcriptional responses in the eight R loci.

Genetic requirements for infection-specific responses in conferring disease resistance in Arabidopsis

Immunity in plants arises from defense regulatory circuits that can be conceptualized as modules. Both the types (and isolates) of pathogen and the repertoire of plant receptors may cause different modules to be activated and affect the magnitude of activation. Two major defense enzymes of Arabidopsis are ALD1 and ICS1/SID2. ALD1 is an aminotransferase needed for producing the metabolites pipecolic acid, hydroxy-pipecolic acid, and possibly other defense signals. ICS1/SID2 produces isochorismate, an intermediate in the synthesis of salicylic acid (SA) and SA-derivatives. Metabolites resulting from the activation of these enzymes are found in petiole exudates and may serve as priming signals for systemic disease resistance in Arabidopsis. Mutants lacking ALD1 are known to have reduced SA accumulation. To further investigate the role of ALD1 in relation to the SA-related module, immunity phenotypes of double mutants that disrupt ALD1 and ICS1/SID2 or SA perception by NPR1 were compared with each single mutant after infection by different Pseudomonas strains. Exudates collected from these mutants after infection were also evaluated for their ability to confer disease resistance when applied to wild-type plants. During infection with virulent or attenuated strains, the loss of ALD1 does not increase the susceptibility of npr1 or sid2 mutants, suggesting the main role of ALD1 in this context is in amplifying the SA-related module. In contrast, after an infection that leads to strong pathogen recognition via the cytoplasmic immune receptor RPS2, ALD1 acts additively with both NPR1 and ICS1/SID2 to suppress pathogen growth. The additive effects are observed in early basal defense responses as well as SA-related events. Thus, there are specific conditions that dictate whether the modules independently contribute to immunity to provide additive protection during infection. In the exudate experiments, intact NPR1 and ICS1/SID2, but not ALD1 in the donor plants were needed for conferring immunity. Mixing exudates showed that loss of SID2 yields exudates that suppress active exudates from wild-type or ald1 plants. This indicates that ICS1/SID2 may not only lead to positive defense signals, but also prevent a suppressive signal(s).

Haplotype-resolved powdery mildew resistance loci reveal the impact of heterozygous structural variation on NLR genes in Muscadinia rotundifolia

Muscadinia rotundifolia cv. Trayshed is a valuable source of resistance to grape powdery mildew. It carries 2 powdery mildew resistance-associated genetic loci, Run1.2 on chromosome 12 and Run2.2 on chromosome 18. The purpose of this study was to identify candidate resistance genes associated with each haplotype of the 2 loci. Both haplotypes of each resistance-associated locus were identified, phased, and reconstructed. Haplotype phasing allowed the identification of several structural variation events between haplotypes of both loci. Combined with a manual refinement of the gene models, we found that the heterozygous structural variants affected the gene content, with some resulting in duplicated or hemizygous nucleotide-binding leucine-rich repeat genes. Heterozygous structural variations were also found to impact the domain composition of some nucleotide-binding leucine-rich repeat proteins. By comparing the nucleotide-binding leucine-rich repeat proteins at Run1.2 and Run2.2 loci, we discovered that the 2 loci include different numbers and classes of nucleotide-binding leucine-rich repeat genes. To identify powdery mildew resistance-associated genes, we performed a gene expression profiling of the nucleotide-binding leucine-rich repeat genes at Run1.2b and Run2.2 loci with or without powdery mildew present. Several nucleotide-binding leucine-rich repeat genes were constitutively expressed, suggesting a role in powdery mildew resistance. These first complete, haplotype-resolved resistance-associated loci and the candidate nucleotide-binding leucine-rich repeat genes identified by this study are new resources that can aid the development of powdery mildew-resistant grape cultivars.

Mapping Locus RSC11K and predicting candidate gene resistant to Soybean mosaic virus strain SC11 through linkage analysis combined with genome resequencing of the parents in soybean

Soybean mosaic virus (SMV) strain SC11 was prevalent in middle China. Its resistance was controlled by a Mendelian single dominant gene RSC11K in soybean Kefeng-1. This study aimed at mapping RSC11K and identifying its candidate gene. RSC11K locus was mapped ~217 kb interval between two SNP-linkage-disequilibrium-blocks (Gm02_BLOCK_11273955_11464884 and Gm02_BLOCK_11486875_11491354) in W82.a1.v1 genome using recombinant inbred lines population derived from Kefeng-1 (Resistant) × NN1138-2 (Susceptible), but inserted with a ~245 kb segment in W82.a2.v1 genome. In the entire 462 kb RSC11K region, 429 SNPs, 142 InDels and 34 putative genes were identified with more SNPs/InDels distributed in non-functional regions. Thereinto, ten genes contained SNP/InDel variants with high and moderate functional impacts on proteins, among which Glyma.02G119700 encoded a typical innate immune receptor-like kinase involving in virus disease process and responded to SMV inoculation, therefore was recognized as RSC11K's candidate gene. The novel RSC11K locus and candidate genes may help developing SMV resistance germplasm.

Genome-wide association mapping reveals genomic regions frequently associated with lettuce field resistance to downy mildew

KEY MESSAGE

GWAS identified 63 QTLs for resistance to downy mildew. Though QTLs were distributed across all chromosomes, the genomic regions frequently associated with resistance were located on chromosomes 4 and 5. Lettuce downy mildew is one of the most economically important diseases of cultivated lettuce worldwide. We have applied the genome-wide association mapping (GWAS) approach to detect QTLs for field resistance to downy mildew in the panel of 496 accessions tested in 21 field experiments. The analysis identified 131 significant marker-trait associations that could be grouped into 63 QTLs. At least 51 QTLs were novel, while remaining 12 QTLs overlapped with previously described QTLs for lettuce field resistance to downy mildew. Unlike race-specific, dominant Dm genes that mostly cluster on three out of nine lettuce chromosomes, QTLs (qDMR loci) for polygenic resistance are randomly distributed across all nine chromosomes. The genomic regions frequently associated with lettuce field resistance to downy mildew are located on chromosomes 4 and 5 and could be used for detailed study of the mechanism of polygenic resistance. The most resistant accessions identified in the current study (cvs. Auburn, Grand Rapids, Romabella, PI 226514, and PI 249536) are being incorporated into our breeding program. Markers closely linked to the resistance QTLs could be potentially used for marker-assisted selection, or in combination with other markers in the genome, for a combined genomic and marker-assisted selection. Up to date this is the most comprehensive study of QTLs for field resistance to downy mildew and the first study that uses GWAS for mapping disease resistance loci in lettuce.

Prediction of NB-LRR resistance genes based on full-length sequence homology

The activation of plant immunity is mediated by resistance (R)-gene receptors, also known as nucleotide-binding leucine-rich repeat (NB-LRR) genes, which in turn trigger the authentic defense response. R-gene identification is a crucial goal for both classic and modern plant breeding strategies for disease resistance. The conventional method identifies NB-LRR genes using a protein motif/domain-based search (PDS) within an automatically predicted gene set of the respective genome assembly. PDS proved to be imprecise since repeat masking prior to automatic genome annotation unwittingly prevented comprehensive NB-LRR gene detection. Furthermore, R-genes have diversified in a species-specific manner, so that NB-LRR gene identification cannot be universally standardized. Here, we present the full-length Homology-based R-gene Prediction (HRP) method for the comprehensive identification and annotation of a genome's R-gene repertoire. Our method has substantially addressed the complex genomic organization of tomato (Solanum lycopersicum) NB-LRR gene loci, proving to be more performant than the well-established RenSeq approach. HRP efficiency was also tested on three differently assembled and annotated Beta sp. genomes. Indeed, HRP identified up to 45% more full-length NB-LRR genes compared to previous approaches. HRP also turned out to be a more refined strategy for R-gene allele mining, testified by the identification of hitherto undiscovered Fom-2 homologs in five Cucurbita sp. genomes. In summary, our high-performance method for full-length NB-LRR gene discovery will propel the identification of novel R-genes towards development of improved cultivars.

RXam2, a NLR from cassava (Manihot esculenta) contributes partially to the quantitative resistance to Xanthomonas phaseoli pv. manihotis

The overexpression of RXam2, a cassava NLR (nucleotide-binding leucine-rich repeat) gene, by stable transformation and gene expression induction mediated by dTALEs, reduce cassava bacterial blight symptoms. Cassava (Manihot esculenta) is a tropical root crop affected by different pathogens including Xanthomonas phaseoli pv. manihotis (Xpm), the causal agent of cassava bacterial blight (CBB). Previous studies have reported resistance to CBB as a quantitative and polygenic character. This study sought to validate the functional role of a NLR (nucleotide-binding leucine-rich repeat) associated with a QTL to Xpm strain CIO151 called RXam2. Transgenic cassava plants overexpressing RXam2 were generated and analyzed. Plants overexpressing RXam2 showed a reduction in bacterial growth to Xpm strains CIO151, 232 and 226. In addition, designer TALEs (dTALEs) were developed to specifically bind to the RXam2 promoter region. The Xpm strain transformed with dTALEs allowed the induction of the RXam2 gene expression after inoculation in cassava plants and was associated with a diminution in CBB symptoms. These findings suggest that RXam2 contributes to the understanding of the molecular basis of quantitative disease resistance.

The barley leaf rust resistance gene Rph3 encodes a predicted membrane protein and is induced upon infection by avirulent pathotypes of Puccinia hordei

Leaf rust, caused by Puccinia hordei, is an economically significant disease of barley, but only a few major resistance genes to P. hordei (Rph) have been cloned. In this study, gene Rph3 was isolated by positional cloning and confirmed by mutational analysis and transgenic complementation. The Rph3 gene, which originated from wild barley and was first introgressed into cultivated Egyptian germplasm, encodes a unique predicted transmembrane resistance protein that differs from all known plant disease resistance proteins at the amino acid sequence level. Genetic profiles of diverse accessions indicated limited genetic diversity in Rph3 in domesticated germplasm, and higher diversity in wild barley from the Eastern Mediterranean region. The Rph3 gene was expressed only in interactions with Rph3-avirulent P. hordei isolates, a phenomenon also observed for transcription activator-like effector-dependent genes known as executors conferring resistance to Xanthomonas spp. Like known transmembrane executors such as Bs3 and Xa7, heterologous expression of Rph3 in N. benthamiana induced a cell death response. The isolation of Rph3 highlights convergent evolutionary processes in diverse plant-pathogen interaction systems, where similar defence mechanisms evolved independently in monocots and dicots.

Leaf rust is an economically significant disease of barley. Here the authors describe cloning of the barley Rph3 leaf rust resistance gene and reveal it encodes a predicted transmembrane protein that is expressed upon infection by Rph3-avirulent Puccinia hordei isolates.

Pathotype Complexity and Genetic Characterization of Phytophthora sojae Populations in Illinois, Indiana, Kentucky, and Ohio

Phytophthora sojae, the causal agent of Phytophthora root and stem rot of soybean, has been managed with single Rps genes since the 1960s but has subsequently adapted to many of these resistance genes, rendering them ineffective. The objective of this study was to examine the pathotype and genetic diversity of P. sojae from soil samples across Illinois, Indiana, Kentucky, and Ohio by assessing which Rps genes were still effective and identifying possible population clusters. There were 218 pathotypes identified from 473 P. sojae isolates with an average of 6.7 out of 15 differential soybean lines exhibiting a susceptible response for each isolate. Genetic characterization of 103 P. sojae isolates from across Illinois, Indiana, Kentucky, and Ohio with 19 simple sequence repeat markers identified 92 multilocus genotypes. There was a moderate level of population differentiation between these four states, with pairwise FST values ranging from 0.026 to 0.246. There were also moderate to high levels of differentiation between fields, with pairwise FST values ranging from 0.071 to 0.537. Additionally, cluster analysis detected the presence of P. sojae population structure across neighboring states. The level of pathotype and genetic diversity, in addition to the identification of population clusters, supports the hypothesis of occasional outcrossing events that allow an increase in diversity and the potential to select for a loss in avirulence to specific resistance genes within regions. The trend of suspected gene flow among neighboring fields is expected to be an ongoing issue with current agricultural practices.

A triple threat: the Parastagonospora nodorum SnTox267 effector exploits three distinct host genetic factors to cause disease in wheat

Parastagonospora nodorum is a fungal pathogen of wheat. As a necrotrophic specialist, it deploys effector proteins that target dominant host susceptibility genes to elicit programmed cell death (PCD). Here we identify and functionally validate the effector targeting the host susceptibility genes Snn2, Snn6 and Snn7. We utilized whole-genome sequencing, association mapping, gene-disrupted mutants, gain-of-function transformants, virulence assays, bioinformatics and quantitative PCR to characterize these interactions. A single proteinaceous effector, SnTox267, targeted Snn2, Snn6 and Snn7 to trigger PCD. Snn2 and Snn6 functioned cooperatively to trigger PCD in a light-dependent pathway, whereas Snn7-mediated PCD functioned in a light-independent pathway. Isolates harboring 20 SnTox267 protein isoforms quantitatively varied in virulence. The diversity and distribution of isoforms varied between populations, indicating adaptation to local selection pressures. SnTox267 deletion resulted in the upregulation of effector genes SnToxA, SnTox1 and SnTox3. We validated a novel effector operating in an inverse-gene-for-gene manner to target three genetically distinct host susceptibility genes and elicit PCD. The discovery of the complementary gene action of Snn2 and Snn6 indicates their potential function in a guard or decoy model. Additionally, differences in light dependency in the elicited pathways and upregulation of unlinked effectors sheds new light onto a complex fungal necrotroph-host interaction.

Genomic variation within the maize stiff-stalk heterotic germplasm pool

The stiff-stalk heterotic group in Maize (Zea mays L.) is an important source of inbreds used in U.S. commercial hybrid production. Founder inbreds B14, B37, B73, and, to a lesser extent, B84, are found in the pedigrees of a majority of commercial seed parent inbred lines. We created high-quality genome assemblies of B84 and four expired Plant Variety Protection (ex-PVP) lines LH145 representing B14, NKH8431 of mixed descent, PHB47 representing B37, and PHJ40, which is a Pioneer Hi-Bred International (PHI) early stiff-stalk type. Sequence was generated using long-read sequencing achieving highly contiguous assemblies of 2.13-2.18 Gbp with N50 scaffold lengths >200 Mbp. Inbred-specific gene annotations were generated using a core five-tissue gene expression atlas, whereas transposable element (TE) annotation was conducted using de novo and homology-directed methodologies. Compared with the reference inbred B73, synteny analyses revealed extensive collinearity across the five stiff-stalk genomes, although unique components of the maize pangenome were detected. Comparison of this set of stiff-stalk inbreds with the original Iowa Stiff Stalk Synthetic breeding population revealed that these inbreds represent only a proportion of variation in the original stiff-stalk pool and there are highly conserved haplotypes in released public and ex-Plant Variety Protection inbreds. Despite the reduction in variation from the original stiff-stalk population, substantial genetic and genomic variation was identified supporting the potential for continued breeding success in this pool. The assemblies described here represent stiff-stalk inbreds that have historical and commercial relevance and provide further insight into the emerging maize pangenome.

Large-scale gene gains and losses molded the NLR defense arsenal during the Cucurbita evolution

Genome-wide annotation reveals that the gene birth-death process of the Cucurbita R family is associated with a species-specific diversification of TNL and CNL protein classes. The Cucurbitaceae family includes nearly 1000 plant species known universally as cucurbits. Cucurbita genus includes many economically important worldwide crops vulnerable to more than 200 pathogens. Therefore, the identification of pathogen-recognition genes is of utmost importance for this genus. The major class of plant-resistance (R) genes encodes nucleotide-binding site and leucine-rich repeat (NLR) proteins, and is divided into three sub-classes namely, TIR-NB-LRR (TNL), CC-NB-LRR (CNL) and RPW8-NB-LRR (RNL). Although the characterization of the NLR gene family has been carried out in important Cucurbita species, this information is still linked to the availability of sequenced genomes. In this study, we analyzed 40 de novo transcriptomes and 5 genome assemblies, which were explored to investigate the Cucurbita expressed-NLR (eNLR) and NLR repertoires using an ad hoc gene annotation approach. Over 1850 NLR-encoding genes were identified, finely characterized and compared to 96 well-characterized plant R-genes. The maximum likelihood analyses revealed an unusual diversification of CNL/TNL genes and a strong RNL conservation. Indeed, several gene gain and loss events have shaped the Cucurbita NLR family. Finally, to provide a first validation step Cucurbita, eNLRs were explored by real-time PCR analysis. The NLR repertories of the 12 Cucurbita species presented in this paper will be useful to discover novel R-genes.

A Single-Nucleotide Mutation in a GLUTAMATE RECEPTOR-LIKE Gene Confers Resistance to Fusarium Wilt in Gossypium hirsutum

Fusarium wilt (FW) disease of cotton, caused by the fungus Fusarium oxysporum f. sp. vasinfectum (Fov), causes severe losses in cotton production worldwide. Though significant advancements have been made in development of FW-resistant Upland cotton (Gossypium hirsutum) in resistance screening programs, the precise resistance genes and the corresponding molecular mechanisms for resistance to Fov remain unclear. Herein it is reported that Fov7, a gene unlike canonical plant disease-resistance (R) genes, putatively encoding a GLUTAMATE RECEPTOR-LIKE (GLR) protein, confers resistance to Fov race 7 in Upland cotton. A single nucleotide polymorphism (SNP) (C/A) in GhGLR4.8, resulting in an amino acid change (L/I), is associated with Fov resistance. A PCR-based DNA marker (GhGLR4.8SNP(A/C) ) is developed and shown to cosegregate with the Fov resistance. CRISPR/Cas9-mediated knockout of Fov7 results in cotton lines extremely susceptible to Fov race 7 with a loss of the ability to induce calcium influx in response to total secreted proteins (SEPs) of Fov. Furthermore, coinfiltration of SEPs with GhGLR4.8A results in a hypersensitive response. This first report of a GLR-encoding gene that functions as an R gene provides a new insight into plant-pathogen interactions and a new handle to develop cotton cultivars with resistance to Fov race 7.

Comparison and evolutionary analysis of Brassica nucleotide binding site leucine rich repeat (NLR) genes and importance for disease resistance breeding

The Brassica genus contains many agriculturally significant oilseed and vegetable crops, however the crop yield is threatened by a range of fungal and bacterial pathogens. Nucleotide Binding Site Leucine Rich Repeat (NLR) genes play important roles in plant innate immunity. The evolution of NLR genes is influenced by genomic processes and pathogen selection. At the whole genome level, whole genome duplications (WGDs) generate abundant gene copies, most of which are lost during genome fractionation. At sub-genomic levels, some retained copies undergo duplication forming clusters which facilitate rapid evolution through recombination. The number, distribution and genetic variations of the NLR genes vary among Brassica species and within populations suggesting differential selection pressure exerted by pathogen populations throughout the evolutionary history. A study of the evolution of disease resistance genes in agriculturally important plants such as Brassicas helps gain insights into their function and inform the identification of resistance genes for breeding of resistant lines.

The Tomato Interspecific NB-LRR Gene Arsenal and Its Impact on Breeding Strategies

Tomato (Solanum lycopersicum L.) is a model system for studying the molecular basis of resistance in plants. The investigation of evolutionary dynamics of tomato resistance (R)-loci provides unique opportunities for identifying factors that promote or constrain genome evolution. Nucleotide-binding domain and leucine-rich repeat (NB-LRR) receptors belong to one of the most plastic and diversified families. The vast amount of genomic data available for Solanaceae and wild tomato relatives provides unprecedented insights into the patterns and mechanisms of evolution of NB-LRR genes. Comparative analysis remarked a reshuffling of R-islands on chromosomes and a high degree of adaptive diversification in key R-loci induced by species-specific pathogen pressure. Unveiling NB-LRR natural variation in tomato and in other Solanaceae species offers the opportunity to effectively exploit genetic diversity in genomic-driven breeding programs with the aim of identifying and introducing new resistances in tomato cultivars. Within this motivating context, we reviewed the repertoire of NB-LRR genes available for tomato improvement with a special focus on signatures of adaptive processes. This issue is still relevant and not thoroughly investigated. We believe that the discovery of mechanisms involved in the generation of a gene with new resistance functions will bring great benefits to future breeding strategies.

Behind the curtain of the compartmentalization process: Exploring how xylem vessel diameter impacts vascular pathogen resistance

A key determinant of plant resistance to vascular infections lies in the ability of the host to successfully compartmentalize invaders at the xylem level. Growing evidence supports that the structural properties of the vascular system impact host vulnerability towards vascular pathogens. The aim of this study was to provide further insight into the impact of xylem vessel diameter on compartmentalization efficiency and thus vascular pathogen movement, using the interaction between Vitis and Phaeomoniella chlamydospora as a model system. We showed experimentally that an increased number of xylem vessels above 100 μm of diameter resulted in a higher mean infection level of host tissue. This benchmark was validated within and across Vitis genotypes. Although the ability of genotypes to restore vascular cambium integrity upon infection was highly variable, this trait did not correlate with their ability to impede pathogen movement at the xylem level. The distribution of infection severity of cuttings across the range of genotype's susceptibility suggests that a risk-based mechanism is involved. We used this experimental data to calibrate a mechanistic stochastic model of the pathogen spread and we provide evidence that the efficiency of the compartmentalization process within a given xylem vessel is a function of its diameter.

Genomic and Molecular Perspectives of Host-pathogen Interaction and Resistance Strategies against White Rust in Oilseed Mustard

Oilseed brassicas stand as the second most valuable source of vegetable oil and the third most traded one across the globe. However, the yield can be severely affected by infections caused by phytopathogens. White rust is a major oomycete disease of oilseed brassicas resulting in up to 60% yield loss globally. So far, success in the development of oomycete resistant Brassicas through conventional breeding has been limited. Hence, there is an imperative need to blend conventional and frontier biotechnological means to breed for improved crop protection and yield.

This review provides a deep insight into the white rust disease and explains the oomycete-plant molecular events with special reference to Albugo candida describing the role of effector molecules, A. candida secretome, and disease response mechanism along with nucleotide-binding leucine-rich repeat receptor (NLR) signaling. Based on these facts, we further discussed the recent progress and future scopes of genomic approaches to transfer white rust resistance in the susceptible varieties of oilseed brassicas, while elucidating the role of resistance and susceptibility genes. Novel genomic technologies have been widely used in crop sustainability by deploying resistance in the host. Enrichment of NLR repertoire, over-expression of R genes, silencing of avirulent and disease susceptibility genes through RNA interference and CRSPR-Cas are technologies which have been successfully applied against pathogen-resistance mechanism. The article provides new insight into Albugo and Brassica genomics which could be useful for producing high yielding and WR resistant oilseed cultivars across the globe.

Effector prediction and characterization in the oomycete pathogen Bremia lactucae reveal host-recognized WY domain proteins that lack the canonical RXLR motif

Pathogens that infect plants and animals use a diverse arsenal of effector proteins to suppress the host immune system and promote infection. Identification of effectors in pathogen genomes is foundational to understanding mechanisms of pathogenesis, for monitoring field pathogen populations, and for breeding disease resistance. We identified candidate effectors from the lettuce downy mildew pathogen Bremia lactucae by searching the predicted proteome for the WY domain, a structural fold found in effectors that has been implicated in immune suppression as well as effector recognition by host resistance proteins. We predicted 55 WY domain containing proteins in the genome of B. lactucae and found substantial variation in both sequence and domain architecture. These candidate effectors exhibit several characteristics of pathogen effectors, including an N-terminal signal peptide, lineage specificity, and expression during infection. Unexpectedly, only a minority of B. lactucae WY effectors contain the canonical N-terminal RXLR motif, which is a conserved feature in the majority of cytoplasmic effectors reported in Phytophthora spp. Functional analysis of 21 effectors containing WY domains revealed 11 that elicited cell death on wild accessions and domesticated lettuce lines containing resistance genes, indicative of recognition of these effectors by the host immune system. Only two of the 11 recognized effectors contained the canonical RXLR motif, suggesting that there has been an evolutionary divergence in sequence motifs between genera; this has major consequences for robust effector prediction in oomycete pathogens.

Wheat Disease Resistance Genes and Their Diversification Through Integrated Domain Fusions

Plants are in a constant evolutionary arms race with their pathogens. At the molecular level, the plant nucleotide-binding leucine-rich repeat receptors (NLRs) family has coevolved with rapidly evolving pathogen effectors. While many NLRs utilize variable leucine-rich repeats (LRRs) to detect effectors, some have gained integrated domains (IDs) that may be involved in receptor activation or downstream signaling. The major objectives of this project were to identify NLR genes in wheat (Triticum aestivum L.) and assess IDs associated with immune signaling (e.g., kinase and transcription factor domains). We identified 2,151 NLR-like genes in wheat, of which 1,298 formed 547 gene clusters. Among the non-toll/interleukin-1 receptor NLR (non-TNL)-like genes, 1,552 encode LRRs, 802 are coiled-coil (CC) domain-encoding (CC-NBS-LRR or CNL) genes, and three encode resistance to powdery mildew 8 (RPW8) domains (RPW8-NBS-LRR or RNL). The expansion of the NLR gene family in wheat is attributable to its origin by recent polyploidy events. Gene clusters were likely formed by tandem duplications, and wheat NLR phylogenetic relationships were similar to those in barley and Aegilops. We also identified wheat NLR-ID fusion proteins as candidates for NLR functional diversification, often as kinase and transcription factor domains. Comparative analyses of the IDs revealed evolutionary conservation of more than 80% amino acid sequence similarity. Homology assessment indicates that these domains originated as functional non-NLR-encoding genes that were incorporated into NLR-encoding genes through duplication events. We also found that many of the NLR-ID genes encode alternative transcripts that include or exclude IDs, a phenomenon that seems to be conserved among species. To verify this, we have analyzed the alternative transcripts that include or exclude an ID of an NLR-ID from another monocotyledon species, rice (Oryza sativa). This indicates that plants employ alternative splicing to regulate IDs, possibly using them as baits, decoys, and functional signaling components. Genomic and expression data support the hypothesis that wheat uses alternative splicing to include and exclude IDs from NLR proteins.

Deciphering the Impact of a Bacterial Infection on Meiotic Recombination in Arabidopsis with Fluorescence Tagged Lines

Plants are under strong evolutionary pressure to maintain surveillance against pathogens. One major disease resistance mechanism is based on NB-LRR (NLR) proteins that specifically recognize pathogen effectors. The cluster organization of the NLR gene family could favor sequence exchange between NLR genes via recombination, favoring their evolutionary dynamics. Increasing data, based on progeny analysis, suggest the existence of a link between the perception of biotic stress and the production of genetic diversity in the offspring. This could be driven by an increased rate of meiotic recombination in infected plants, but this has never been strictly demonstrated. In order to test if pathogen infection can increase DNA recombination in pollen meiotic cells, we infected Arabidopsis Fluorescent Tagged Lines (FTL) with the virulent bacteria Pseudomonas syringae. We measured the meiotic recombination rate in two regions of chromosome 5, containing or not an NLR gene cluster. In all tested intervals, no significant difference in genetic recombination frequency between infected and control plants was observed. Although it has been reported that pathogen exposure can sometimes increase the frequency of recombinant progeny in plants, our findings suggest that meiotic recombination rate in Arabidopsis may be resilient to at least some pathogen attack. Alternative mechanisms are discussed.

Dissecting the Genetic Architecture of Aphanomyces Root Rot Resistance in Lentil by QTL Mapping and Genome-Wide Association Study

Lentil (Lens culinaris Medikus) is an important source of protein for people in developing countries. Aphanomyces root rot (ARR) has emerged as one of the most devastating diseases affecting lentil production. In this study, we applied two complementary quantitative trait loci (QTL) analysis approaches to unravel the genetic architecture underlying this complex trait. A recombinant inbred line (RIL) population and an association mapping population were genotyped using genotyping by sequencing (GBS) to discover novel single nucleotide polymorphisms (SNPs). QTL mapping identified 19 QTL associated with ARR resistance, while association mapping detected 38 QTL and highlighted accumulation of favorable haplotypes in most of the resistant accessions. Seven QTL clusters were discovered on six chromosomes, and 15 putative genes were identified within the QTL clusters. To validate QTL mapping and genome-wide association study (GWAS) results, expression analysis of five selected genes was conducted on partially resistant and susceptible accessions. Three of the genes were differentially expressed at early stages of infection, two of which may be associated with ARR resistance. Our findings provide valuable insight into the genetic control of ARR, and genetic and genomic resources developed here can be used to accelerate development of lentil cultivars with high levels of partial resistance to ARR.

Plants make galls to accommodate foreigners: some are friends, most are foes

At the colonization site of a foreign entity, plant cells alter their trajectory of growth and development. The resulting structure - a plant gall - accommodates various needs of the foreigner, which are phylogenetically diverse: viruses, bacteria, protozoa, oomycetes, true fungi, parasitic plants, and many types of animals, including rotifers, nematodes, insects, and mites. The plant species that make galls also are diverse. We assume gall production costs the plant. All is well if the foreigner provides a gift that makes up for the cost. Nitrogen-fixing nodule-inducing bacteria provide nutritional services. Gall wasps pollinate fig trees. Unfortunately for plants, most galls are made for foes, some of which are deeply studied pathogens and pests: Agrobacterium tumefaciens, Rhodococcus fascians, Xanthomonas citri, Pseudomonas savastanoi, Pantoea agglomerans, 'Candidatus' phytoplasma, rust fungi, Ustilago smuts, root knot and cyst nematodes, and gall midges. Galls are an understudied phenomenon in plant developmental biology. We propose gall inception for discovering unifying features of the galls that plants make for friends and foes, talk about molecules that plants and gall-inducers use to get what they want from each other, raise the question of whether plants colonized by arbuscular mycorrhizal fungi respond in a gall-like manner, and present a research agenda.

Comparative Transcriptome Profiling Reveals Compatible and Incompatible Patterns of Potato Toward Phytophthora infestans

Late blight, caused by Phytophthora infestans (P. infestans), is a devastating disease in potato worldwide. Our previous study revealed that the Solanum andigena genotype 03112-233 is resistant to P. infestans isolate 90128, but susceptible to the super race isolate, CN152. In this study, we confirmed by diagnostic resistance gene enrichment sequencing (dRenSeq) that the resistance of 03112-233 toward 90128 is most likely based on a distinct new R gene(s). To gain an insight into the mechanism that governs resistance or susceptibility in 03112-223, comparative transcriptomic profiling analysis based on RNAseq was initiated. Changes in transcription at two time points (24 h and 72 h) after inoculation with isolates 90128 or CN152 were analyzed. A total of 8,881 and 7,209 genes were differentially expressed in response to 90128 and CN152, respectively, and 1,083 differentially expressed genes (DEGs) were common to both time points and isolates. A substantial number of genes were differentially expressed in an isolate-specific manner with 3,837 genes showing induction or suppression following infection with 90128 and 2,165 genes induced or suppressed after colonization by CN152. Hierarchical clustering analysis suggested that isolates with different virulence profiles can induce different defense responses at different time points. Further analysis revealed that the compatible interaction caused higher induction of susceptibility genes such as SWEET compared with the incompatible interaction. The salicylic acid, jasmonic acid, and abscisic acid mediated signaling pathways were involved in the response against both isolates, while ethylene and brassinosteroids mediated defense pathways were suppressed. Our results provide a valuable resource for understanding the interactions between P. infestans and potato.

Inferring RPW8-NLRs's evolution patterns in seed plants: case study in Vitis vinifera

Genomic and transcriptomic studies in plants and, more in deep, in grapevine reveal that the disease-resistance RNL gene family is highly variable. RNLs (RPW8-NLRs) are a phylogenetically distinct class of nucleotide oligomerization domain (NOD)-like receptors (NLRs) identified in plants. Two RNLs, namely, the NRG1 (N Requirement Gene 1) and the ADR1 (Activated Disease Resistance 1), have been characterized; however, little is known about the RNL evolutionary history in higher plants. To trace the diversification of RNL gene subfamily, we scanned the NLR proteins of 73 plant genomes belonging to 29 taxa, revealing a noticeable diversification across species and within the same genus or botanic family together with a conspicuous expansion in important crop species. To explore the RNL variability in Vitis vinifera and gain information with respect to their structure, evolutionary diversification of five grape genomes ('Aglianico', 'Falanghina', 'Sultanina', 'Tannat', and 'Nebbiolo') has been compared to the reference genome ('Pinot Noir'). The number of RNLs ranged from 6 ('Sultanina') to 14 ('Nebbiolo'), in contrast to the 10 'Pinot Noir' RNLs. The phylogenetic study on grapevine RNLs revealed that all collapsed into NRG1-clade, rather than four. To investigate more in depth the means of intraspecific variability of grape RNL copies, a transcriptomic profiling in response to powdery mildew (PM) infection was carried out through qRT-PCRs and public databases interrogation. The RNL expression variability identified in transcriptome data sets supports the hypothesis of a functional expansion/contraction in grapevine varieties. Although no direct correlations between grapevine PM-resistance and RNL expression was identified, our work can provide good candidates for functional studies able to elucidate the putative "helper" role of RNLs in grape immune signalling.

BjuWRR1, a CC-NB-LRR gene identified in Brassica juncea, confers resistance to white rust caused by Albugo candida

BjuWRR1, a CNL-type R gene, was identified from an east European gene pool line of Brassica juncea and validated for conferring resistance to white rust by genetic transformation. White rust caused by the oomycete pathogen Albugo candida is a significant disease of crucifer crops including Brassica juncea (mustard), a major oilseed crop of the Indian subcontinent. Earlier, a resistance-conferring locus named AcB1-A5.1 was mapped in an east European gene pool line of B. juncea-Donskaja-IV. This line was tested along with some other lines of B. juncea (AABB), B. rapa (AA) and B. nigra (BB) for resistance to six isolates of A. candida collected from different mustard growing regions of India. Donskaja-IV was found to be completely resistant to all the tested isolates. Sequencing of a BAC spanning the locus AcB1-A5.1 showed the presence of a single CC-NB-LRR protein encoding R gene. The genomic sequence of the putative R gene with its native promoter and terminator was used for the genetic transformation of a susceptible Indian gene pool line Varuna and was found to confer complete resistance to all the isolates. This is the first white rust resistance-conferring gene described from Brassica species and has been named BjuWRR1. Allelic variants of the gene in B. juncea germplasm and orthologues in the Brassicaceae genomes were studied to understand the evolutionary dynamics of the BjuWRR1 gene.

Genetic toolkits of the red alga Pyropia tenera against the three most common diseases in Pyropia farms

Disease outbreaks devastate Pyropia aquaculture farms every year. The three most common and serious diseases are Olpidiopsis-blight and red-rot disease caused by oomycete pathogens and green-spot disease caused by the PyroV1 virus. We hypothesized that a basic genetic profile of molecular defenses will be revealed by comparing and analyzing the genetic response of Pyropia tenera against the above three pathogens. RNAs isolated from infected thalli were hybridized onto an oligochip containing 15,115 primers designed from P. tenera expressed sequence tags (EST)s. Microarray profiles of the three diseases were compared and interpreted together with histochemical observation. Massive amounts of reactive oxygen species accumulated in P. tenera cells exposed to oomycete pathogens. Heat shock genes and serine proteases were the most highly up-regulated genes in all infection experiments. Genes involved in RNA metabolism, ribosomal proteins and antioxidant metabolism were also highly up-regulated. Genetic profiles of P. tenera in response to pathogens were most similar between the two biotrophic pathogens, Olpidiopsis pyropiae and PyroV1 virus. A group of plant resistance genes were specifically regulated against each pathogen. Our results suggested that disease response in P. tenera consists of a general constitutive defense and a genetic toolkit against specific pathogens.

Point Inoculation Method for Measuring Adult Plant Response of Wheat to Stripe Rust Infection

Depending on the pathogenicity of the stripe rust fungus Puccinia striiformis f. sp. tritici, the nature of resistance in the wheat host plant, and the environment, a broad range of disease phenotypes can be expressed. Therefore, the phenotyping of partial adult plant stripe rust resistance requires reliable and repeatable procedures, especially under controlled conditions. In this study, the development of a flag leaf point inoculation method, which resulted in a 100% initial infection rate, is reported. Flag leaf inoculations were achieved by placing 6-mm antibiotic test paper discs, dipped into a urediniospore and water suspension and covered with water-proof plastic tape, on the adaxial side of leaves. Results from independent trials allowed for the statistical comparison of stripe rust lesion expansion rate in wheat entries that differ in resistance. The technique is inexpensive, reliable, and applicable to routine screening for adult plant response type, quantitative comparison of stripe rust progress, environmental influences, and pathogenicity of different isolates.

Molecular mechanisms that limit the costs of NLR-mediated resistance in plants

Crop diseases cause significant yield losses, and the use of resistant cultivars can effectively mitigate these losses and control many plant diseases. Most plant resistance (R) genes encode immune receptors composed of nucleotide-binding and leucine-rich repeat (NLR) domains. These proteins mediate the specific recognition of pathogen avirulence effectors to induce defence responses. However, NLR-triggered immunity can be associated with a reduction in growth and yield, so-called 'fitness costs'. Recent data have shown that plants use an elaborate interplay of different mechanisms to control NLR gene transcript levels, as well as NLR protein abundance and activity, to avoid the associated cost of resistance in the absence of a pathogen. In this review, we discuss the different levels of NLR regulation (transcriptional, post-transcriptional and at the protein level). We address the apparent need for plants to maintain diverse modes of regulation. A recent model suggesting an equilibrium 'ON/OFF state' of NLR proteins, in the absence of a pathogen, provides the context for our discussion.

Pyramiding of nine transgenes in maize generates high-level resistance against necrotrophic maize pathogens

Key message Nine transgenes from different categories, viz. plant defense response genes and anti-apoptosis genes, played combined roles in maize to inhibit the necrotrophic pathogens Rhizoctonia solani and Bipolaris maydis. Maize sheath blight and southern corn leaf blight are major global threats to maize production. The management of these necrotrophic pathogens has encountered limited success due to the characteristics of their lifestyle. Here, we presented a transgenic pyramiding breeding strategy to achieve nine different resistance genes integrated in one transgenic maize line to combat different aspects of necrotrophic pathogens. These nine genes, selected from two different categories, plant defense response genes (Chi, Glu, Ace-AMP1, Tlp, Rs-AFP2, ZmPROPEP1 and Pti4), and anti-apoptosis genes (Iap and p35), were successfully transferred into maize and further implicated in resistance against the necrotrophic pathogens Rhizoctonia solani and Bipolaris maydis. Furthermore, the transgenic maize line 910, with high expression levels of the nine integrated genes, was selected from 49 lines. Under greenhouse and field trial conditions, line 910 showed significant resistance against maize sheath blight and southern corn leaf blight diseases. Higher-level resistance was obtained after the pyramiding of more resistance transgenes from different categories that function via different mechanisms. The present study provides a successful strategy for the management of necrotrophic pathogens.

Genome-Wide Identification of NBS-Encoding Resistance Genes in Sunflower (Helianthus annuus L.)

Nucleotide Binding Site—Leucine-Rich Repeat (NBS-LRR) genes encode disease resistance proteins involved in plants’ defense against their pathogens. Although sunflower is affected by many diseases, only a few molecular details have been uncovered regarding pathogenesis and resistance mechanisms. Recent availability of sunflower whole genome sequences in publicly accessible databases allowed us to accomplish a genome-wide identification of Toll-interleukin-1 receptor-like Nucleotide-binding site Leucine-rich repeat (TNL), Coiled Coil (CC)-NBS-LRR (CNL), Resistance to powdery mildew8 (RPW8)-NBS-LRR (RNL) and NBS-LRR (NL) protein encoding genes. Hidden Markov Model (HMM) profiling of 52,243 putative protein sequences from sunflower resulted in 352 NBS-encoding genes, among which 100 genes belong to CNL group including 64 genes with RX_CC like domain, 77 to TNL, 13 to RNL, and 162 belong to NL group. We also identified signal peptides and nuclear localization signals present in the identified genes and their homologs. We found that NBS genes were located on all chromosomes and formed 75 gene clusters, one-third of which were located on chromosome 13. Phylogenetic analyses between sunflower and Arabidopsis NBS genes revealed a clade-specific nesting pattern in CNLs, with RNLs nested in the CNL-A clade, and species-specific nesting pattern for TNLs. Surprisingly, we found a moderate bootstrap support (BS = 50%) for CNL-A clade being nested within TNL clade making both the CNL and TNL clades paraphyletic. Arabidopsis and sunflower showed 87 syntenic blocks with 1049 high synteny hits between chromosome 5 of Arabidopsis and chromosome 6 of sunflower. Expression data revealed functional divergence of the NBS genes with basal level tissue-specific expression. This study represents the first genome-wide identification of NBS genes in sunflower paving avenues for functional characterization and potential crop improvement.

Disease Resistance Mechanisms in Plants

Plants have developed a complex defense system against diverse pests and pathogens. Once pathogens overcome mechanical barriers to infection, plant receptors initiate signaling pathways driving the expression of defense response genes. Plant immune systems rely on their ability to recognize enemy molecules, carry out signal transduction, and respond defensively through pathways involving many genes and their products. Pathogens actively attempt to evade and interfere with response pathways, selecting for a decentralized, multicomponent immune system. Recent advances in molecular techniques have greatly expanded our understanding of plant immunity, largely driven by potential application to agricultural systems. Here, we review the major plant immune system components, state of the art knowledge, and future direction of research on plant⁻pathogen interactions. In our review, we will discuss how the decentralization of plant immune systems have provided both increased evolutionary opportunity for pathogen resistance, as well as additional mechanisms for pathogen inhibition of such defense responses. We conclude that the rapid advances in bioinformatics and molecular biology are driving an explosion of information that will advance agricultural production and illustrate how complex molecular interactions evolve.

Pyramiding for Resistance Durability: Theory and Practice

Durable disease resistance is a key component of global food security, and combining resistance genes into "pyramids" is an important way to increase durability of resistance. The mechanisms by which pyramids impart durability are not well known. The traditional view of resistance pyramids considers the use of major resistance gene (R-gene) combinations deployed against pathogens that are primarily asexual. Interestingly, published examples of the successful use of pyramids in the traditional sense are rare. In contrast, most published descriptions of durable pyramids in practice are for cereal rusts, and tend to indicate an association between durability and cultivars combining major R-genes with incompletely expressed, adult plant resistance genes. Pyramids have been investigated experimentally for a diversity of pathogens, and many reduce disease levels below that of the single best gene. Resistance gene combinations have been identified through phenotypic reactions, molecular markers, and challenge against effector genes. As resistance genes do not express equally in all genetic backgrounds, however, a combination of genetic information and phenotypic analyses provide the ideal scenario for testing of putative pyramids. Not all resistance genes contribute equally to pyramids, and approaches have been suggested to identify the best genes and combinations of genes for inclusion. Combining multiple resistance genes into a single plant genotype quickly is a challenge that is being addressed through alternative breeding approaches, as well as through genomics tools such as resistance gene cassettes and gene editing. Experimental and modeling tests of pyramid durability are in their infancy, but have promise to help direct future studies of pyramids. Several areas for further work on resistance gene pyramids are suggested.

The TNL gene Rdr1 confers broad-spectrum resistance to Diplocarpon rosae

Black spot disease, which is caused by the ascomycete Diplocarpon rosae, is the most severe disease in field-grown roses in temperate regions and has been distributed worldwide, probably together with commercial cultivars. Here, we present data indicating that muRdr1A is the active Rdr1 gene, a single-dominant TIR-NBS-LRR (Toll/interleukin-1 receptor-nucleotide binding site-leucine rich repeat) (TNL)-type resistance gene against black spot disease, which acts against a broad range of pathogenic isolates independent of the genetic background of the host genotype. Molecular analyses revealed that, compared with the original donor genotype, the multiple integrations that are found in the primary transgenic clone segregate into different integration patterns in its sexual progeny and do not show any sign of overexpression. Rdr1 provides resistance to 13 different single-spore isolates belonging to six different races and broad field mixtures of conidia; thus far, Rdr1 is only overcome by two races. The expression of muRdr1A, the active Rdr1 gene, leads to interaction patterns that are identical in the transgenic clones and the non-transgenic original donor genotype. This finding indicates that the interacting avirulence (Avr) factor on the pathogen side must be widespread among the pathogen populations and may have a central function in the rose-black spot interaction. Therefore, the Rdr1 gene, pyramided with only a few other R genes by sexual crosses, might be useful for breeding roses that are resistant to black spot because the spread of new pathogenic races of the fungus appears to be slow.

Genomic and epigenomic immunity in common bean: the unusual features of NB-LRR gene family

In plants, a key class of genes comprising most of disease resistance (R) genes encodes Nucleotide-binding leucine-rich repeat (NL) proteins. Access to common bean (Phaseolus vulgaris) genome sequence provides unparalleled insight into the organization and evolution of this large gene family (∼400 NL) in this important crop. As observed in other plant species, most common bean NL are organized in cluster of genes. However, a particularity of common bean is that these clusters are often located in subtelomeric regions close to terminal knobs containing the satellite DNA khipu. Phylogenetically related NL are spread between different chromosome ends, suggesting frequent exchanges between non-homologous chromosomes. NL peculiar location, in proximity to heterochromatic regions, led us to study their DNA methylation status using a whole-genome cytosine methylation map. In common bean, NL genes displayed an unusual body methylation pattern since half of them are methylated in the three contexts, reminiscent of the DNA methylation pattern of repeated sequences. Moreover, 90 NL were also abundantly targeted by 24 nt siRNA, with 90% corresponding to methylated NL genes. This suggests the existence of a transcriptional gene silencing mechanism of NL through the RdDM (RNA-directed DNA methylation) pathway in common bean that has not been described in other plant species.

Wheat receptor-kinase-like protein Stb6 controls gene-for-gene resistance to fungal pathogen Zymoseptoria tritici

Deployment of fast-evolving disease-resistance genes is one of the most successful strategies used by plants to fend off pathogens^1,2. In gene-for-gene relationships, most cloned disease-resistance genes encode intracellular nucleotide-binding leucine-rich-repeat proteins (NLRs) recognizing pathogen-secreted isolate-specific avirulence (Avr) effectors delivered to the host cytoplasm^3,4. This process often triggers a localized hypersensitive response, which halts further disease development ⁵ . Here we report the map-based cloning of the wheat Stb6 gene and demonstrate that it encodes a conserved wall-associated receptor kinase (WAK)-like protein, which detects the presence of a matching apoplastic effector^6-8 and confers pathogen resistance without a hypersensitive response ⁹ . This report demonstrates gene-for-gene disease resistance controlled by this class of proteins in plants. Moreover, Stb6 is, to our knowledge, the first cloned gene specifying resistance to Zymoseptoria tritici, an important foliar fungal pathogen affecting wheat and causing economically damaging septoria tritici blotch (STB) disease^10-12.

Characterization of Phytophthora nicotianae Resistance Conferred by the Introgressed Nicotiana rustica Region, Wz, in Flue-Cured Tobacco

Black shank, caused by Phytophthora nicotianae, is one of the most important diseases affecting tobacco worldwide and is primarily managed through use of host resistance. An additional source of resistance to P. nicotianae, designated as Wz, has been introgressed into Nicotiana tabacum from N. rustica. The Wz gene region confers high levels of resistance to all races, but has not been characterized. Our study found Wz-mediated resistance is most highly expressed in the roots, with only a slight reduction in stem-lesion size in Wz genotypes compared with susceptible controls. No substantial relationships were observed between initial inoculum levels and disease development on Wz genotypes, which is generally consistent with qualitative or complete resistance. Isolates of P. nicotianae adapted for five host generations on plants with the Wz gene caused higher disease severity than isolates adapted on Wz plants for only one host generation. Wz-adapted isolates did not exhibit increased aggressiveness on genotypes with other sources of partial resistance, suggesting pathogen adaptation was specific to the Wz gene. To reduce potential for pathogen population shifts with virulence on Wz genotypes, Wz should be combined with other resistance sources and rotation of varying black shank resistance mechanisms is also recommended.

Genetic resistance against viruses in Phaseolus vulgaris L.: State of the art and future prospects

Viruses are obligate parasites that replicate intracellularly in many living organisms, including plants. Consequently, no chemicals are available that target only the virus without impacting host cells or vector organisms. The use of natural resistant varieties appears as the most reliable control strategy and remains the best and cheapest option in managing virus diseases, especially in the current ecological context of preserving biodiversity and environment in which the use of phytosanitary products becomes limited. Common bean is a grain legume cultivated mainly in Africa and Central-South America. Virus diseases of common bean have been extensively studied both by breeders to identify natural resistance genes in existing germplasms and by pathologists to understand the molecular bases of plant-virus interactions. Here we present a critical review in which we synthesize previous and recent information concerning 1) main viruses causing diseases in common bean, 2) genetic resistance to viruses in common bean, 3) the different resistance phenotypes observed and more particularly the effect of temperature, 4) the molecular bases of resistance genes to viruses in common bean, and 5) future prospects using transgenic-engineered resistant lines.

Identification and characterization of Sr13, a tetraploid wheat gene that confers resistance to the Ug99 stem rust race group

The Puccinia graminis f. sp. tritici (Pgt) Ug99 race group is virulent to most stem rust resistance genes currently deployed in wheat and poses a threat to global wheat production. The durum wheat (Triticum turgidum ssp. durum) gene Sr13 confers resistance to Ug99 and other virulent races, and is more effective at high temperatures. Using map-based cloning, we delimited a candidate region including two linked genes encoding coiled-coil nucleotide-binding leucine-rich repeat proteins designated CNL3 and CNL13. Three independent truncation mutations identified in each of these genes demonstrated that only CNL13 was required for Ug99 resistance. Transformation of an 8-kb genomic sequence including CNL13 into the susceptible wheat variety Fielder was sufficient to confer resistance to Ug99, confirming that CNL13 is Sr13CNL13 transcripts were slightly down-regulated 2-6 days after Pgt inoculation and were not affected by temperature. By contrast, six pathogenesis-related (PR) genes were up-regulated at high temperatures only when both Sr13 and Pgt were present, suggesting that they may contribute to the high temperature resistance mechanism. We identified three Sr13-resistant haplotypes, which were present in one-third of cultivated emmer and durum wheats but absent in most tested common wheats (Triticum aestivum). These results suggest that Sr13 can be used to improve Ug99 resistance in a large proportion of modern wheat cultivars. To accelerate its deployment, we developed a diagnostic marker for Sr13 The identification of Sr13 expands the number of Pgt-resistance genes that can be incorporated into multigene transgenic cassettes to control this devastating disease.

Quantitative Resistance to Plant Pathogens in Pyramiding Strategies for Durable Crop Protection

Quantitative resistance has gained interest in plant breeding for pathogen control in low-input cropping systems. Although quantitative resistance frequently has only a partial effect and is difficult to select, it is considered more durable than major resistance (R) genes. With the exponential development of molecular markers over the past 20 years, resistance QTL have been more accurately detected and better integrated into breeding strategies for resistant varieties with increased potential for durability. This review summarizes current knowledge on the genetic inheritance, molecular basis, and durability of quantitative resistance. Based on this knowledge, we discuss how strategies that combine major R genes and QTL in crops can maintain the effectiveness of plant resistance to pathogens. Combining resistance QTL with complementary modes of action appears to be an interesting strategy for breeding effective and potentially durable resistance. Combining quantitative resistance with major R genes has proven to be a valuable approach for extending the effectiveness of major genes. In the plant genomics era, improved tools and methods are becoming available to better integrate quantitative resistance into breeding strategies. Nevertheless, optimal combinations of resistance loci will still have to be identified to preserve resistance effectiveness over time for durable crop protection.

A wheat NBS-LRR gene TaRGA19 participates in Lr19-mediated resistance to Puccinia triticina

Wheat leaf rust, caused by Puccinia triticina (Pt), is one of the most severe fungal diseases on wheat globally. Rational utilization of wheat leaf rust resistance (Lr) genes is still the best choice for control this disease. Wheat seedlings carrying Lr19 showed a high resistance phenotype to all Pt races in China. So far, all the cloned seedling Lr genes including Lr1, Lr10 and Lr21, encode protein with NBS-LRR domain. In this study, a wheat gene with NBS-LRR domain from previously established Lr19-resistance-related cDNA library was cloned and designated as TaRGA19. Full length of this gene was amplified by rapid amplification of cDNA ends (RACE). By blast against IWGSC wheat genome database, we have noticed that TaRGA19 was located on chromosome 2DS, which was different from Lr19 located on chromosome 7DL. Compared with susceptible Thatcher line, expression level of TaRGA19 was upregulated in wheat isogenic lines carrying Lr19 (TcLr19) after inoculation of Pt race THTS. By particle bombardment, TaRGA19-GFP fused protein was localized on plasma membrane of epidermal cells. Using virus-induced gene silencing (VIGS), TaRGA19-knockdown plants of TcLr19 showed reduced resistance and few sporulation phenotype upon Pt challenge. Further histological observation indicated that Pt hyphal growth at the infection sites was less suppressed in the TaRGA19-knockdown plants. In conclusion, we speculate this TaRGA19 gene was involved in the Lr19-mediated resistance to wheat leaf rust along with other components.

Comparative Genomics of Non-TNL Disease Resistance Genes from Six Plant Species

Disease resistance genes (R genes), as part of the plant defense system, have coevolved with corresponding pathogen molecules. The main objectives of this project were to identify non-Toll interleukin receptor, nucleotide-binding site, leucine-rich repeat (nTNL) genes and elucidate their evolutionary divergence across six plant genomes. Using reference sequences from Arabidopsis, we investigated nTNL orthologs in the genomes of common bean, Medicago, soybean, poplar, and rice. We used Hidden Markov Models for sequence identification, performed model-based phylogenetic analyses, visualized chromosomal positioning, inferred gene clustering, and assessed gene expression profiles. We analyzed 908 nTNL R genes in the genomes of the six plant species, and classified them into 12 subgroups based on the presence of coiled-coil (CC), nucleotide binding site (NBS), leucine rich repeat (LRR), resistance to Powdery mildew 8 (RPW8), and BED type zinc finger domains. Traditionally classified CC-NBS-LRR (CNL) genes were nested into four clades (CNL A-D) often with abundant, well-supported homogeneous subclades of Type-II R genes. CNL-D members were absent in rice, indicating a unique R gene retention pattern in the rice genome. Genomes from Arabidopsis, common bean, poplar and soybean had one chromosome without any CNL R genes. Medicago and Arabidopsis had the highest and lowest number of gene clusters, respectively. Gene expression analyses suggested unique patterns of expression for each of the CNL clades. Differential gene expression patterns of the nTNL genes were often found to correlate with number of introns and GC content, suggesting structural and functional divergence.

Coevolution between Mutualists and Parasites in Symbiotic Communities May Lead to the Evolution of Lower Virulence

Most eukaryotes harbor a diverse community of parasitic, mutualistic, and commensal microbial symbionts. Although the diversity of these microbial symbiotic communities has recently drawn considerable attention, theory regarding the evolution of interactions among symbionts and with the host is still in its nascent stages. Here we evaluate the role of interactions among coinfecting symbionts in the evolution of symbiont virulence toward the host. To do so, we place the virulence-transmission trade-off into a community context and model the evolution of symbiont trophic modes along the continuum from parasitism (virulence) to mutualism (negative virulence). We establish a framework for studying multiple infections of a host by the same symbiont species and coinfection by multiple species, using a concept of shared costs, wherein the negative consequences of virulence (or harm) toward the host are shared among symbionts. Our results show that mutualism can be maintained under infection by multiple symbionts when shared costs are sufficiently low, while greater virulence and parasitism toward the host are more likely when shared costs are high. Last, for coinfection by more than one species, we show that if the presence of a mutualist ameliorates some of the costs of pathogen virulence, then the symbiotic community may more often evolve to a more commensal state and maintain mutualisms.

Xylem Vessel Diameter Affects the Compartmentalization of the Vascular Pathogen Phaeomoniella chlamydospora in Grapevine

Fungal wilt diseases are a threat to global food safety. Previous studies in perennial crops showed that xylem vessel diameter affects disease susceptibility. We tested the hypothesis that xylem vessel diameter impacts occlusion processes and pathogen compartmentalization in Vitis vinifera L. We studied the interaction between four grape commercial cultivars with the vascular wilt pathogen Phaeomoniella chlamydospora. We used qPCR and wood necrotic lesion length to measure fungal colonization coupled with histological studies to assess differences in xylem morphology, pathogen compartmentalization, and fungal colonization strategy. We provided evidence that grape cultivar with wide xylem vessel diameter showed increased susceptibility to P. chlamydospora. The host response to pathogen included vessel occlusion with tyloses and gels, deposition of non-structural phenolic compounds and suberin in vessel walls and depletion of starch in parenchyma cells. Pathogen compartmentalization was less efficient in wide xylem vessels than in narrow diameter vessels. Large vessels displayed higher number of tyloses and gel pockets, which provided substrate for P. chlamydospora growth and routes to escape occluded vessels. We discuss in which capacity xylem vessel diameter is a key determinant of the compartmentalization process and in turn grape cultivar resistance to disease caused by P. chlamydospora.

Inheritance, fine-mapping, and candidate gene analyses of resistance to soybean mosaic virus strain SC5 in soybean

Soybean mosaic virus (SMV) is one of the most devastating pathogens for soybeans in China. Among the country-wide 22 strains, SC5 dominates in Huang-Huai and Changjiang valleys. For controlling its damage, the resistance gene was searched through Mendelian inheritance study, gene fine-mapping, and candidate gene analysis combined with qRT-PCR (quantitative real-time polymerase chain reaction) analysis. The parents F1, F2, and RILs (recombinant inbred lines) of the cross Kefeng-1 (Resistance, R) × NN1138-2 (Susceptible, S) were used to examine the inheritance of SC5-resistance. The F1 was resistant and the F2 and RILs segregated in a 3R:1S and 1R:1S ratio, respectively, indicating a single dominant gene conferring the Kefeng-1 resistance. Subsequently, the genomic region conferring the resistance was found in "Bin 352-Bin353 with 500 kb" on Chromosome 2 using the phenotyping data of the 427 RILs and a high-density genetic map with 4703 bin markers. In the 500 kb genomic region, 38 putative genes are contained. The association analysis between the SNPs in a putative gene and the resistance phenotype for the 427 RILs prioritized 11 candidate genes using Chi-square criterion. The expression levels of these genes were tested by qRT-PCR. On infection with SC5, 7 out of the 11 genes had differential expression in Kefeng-1 and NN1138-2. Furthermore, integrating SNP-phenotype association analysis with qRT-PCR expression profiling analysis, Glyma02g13495 was found the most possible candidate gene for SC5-resistance. This finding can facilitate the breeding for SC5-resistance through marker-assisted selection and provide a platform to gain a better understanding of SMV-resistance gene system in soybean.