Resistance gene enrichment sequencing (RenSeq) enables reannotation of the NB‐LRR gene family from sequenced plant genomes and rapid mapping of resistance loci in segregating populations

RenSeq is a NB‐LRR (nucleotide binding‐site leucine‐rich repeat) gene‐targeted, Resistance gene enrichment and sequencing method that enables discovery and annotation of pathogen resistance gene family members in plant genome sequences. We successfully applied RenSeq to the sequenced potato Solanum tuberosum clone DM, and increased the number of identified NB‐LRRs from 438 to 755. The majority of these identified R gene loci reside in poorly or previously unannotated regions of the genome. Sequence and positional details on the 12 chromosomes have been established for 704 NB‐LRRs and can be accessed through a genome browser that we provide. We compared these NB‐LRR genes and the corresponding oligonucleotide baits with the highest sequence similarity and demonstrated that ~80% sequence identity is sufficient for enrichment. Analysis of the sequenced tomato S. lycopersicum ‘Heinz 1706’ extended the NB‐LRR complement to 394 loci. We further describe a methodology that applies RenSeq to rapidly identify molecular markers that co‐segregate with a pathogen resistance trait of interest. In two independent segregating populations involving the wild Solanum species S. berthaultii (Rpi‐ber2) and S. ruiz‐ceballosii (Rpi‐rzc1), we were able to apply RenSeq successfully to identify markers that co‐segregate with resistance towards the late blight pathogen Phytophthora infestans. These SNP identification workflows were designed as easy‐to‐adapt Galaxy pipelines.     

The Brassica napus receptor‐like protein RLM2 is encoded by a second allele of the LepR3/Rlm2 blackleg resistance locus

Leucine‐rich repeat receptor‐like proteins (LRR‐RLPs) are highly adaptable parts of the signalling apparatus for extracellular detection of plant pathogens. Resistance to blackleg disease of Brassica spp. caused by Leptosphaeria maculans is largely governed by host race‐specific R‐genes, including the LRR‐RLP gene LepR3. The blackleg resistance gene Rlm2 was previously mapped to the same genetic interval as LepR3. In this study, the LepR3 locus of the Rlm2 Brassica napus line ‘Glacier DH24287’ was cloned, and B. napus transformants were analysed for recovery of the Rlm2 phenotype. Multiple B. napus, B. rapa and B. juncea lines were assessed for sequence variation at the locus. Rlm2 was found to be an allelic variant of the LepR3 LRR‐RLP locus, conveying race‐specific resistance to L. maculans isolates harbouring AvrLm2. Several defence‐related LRR‐RLPs have previously been shown to associate with the RLK SOBIR1 to facilitate defence signalling. Bimolecular fluorescence complementation (BiFC) and co‐immunoprecipitation of RLM2‐SOBIR1 studies revealed that RLM2 interacts with SOBIR1 of Arabidopsis thaliana when co‐expressed in Nicotiana benthamiana. The interaction of RLM2 with AtSOBIR1 is suggestive of a conserved defence signalling pathway between B. napus and its close relative A. thaliana.     

Nucleotide‐binding resistance gene signatures in sugar beet,insights from a new reference genome

Nucleotide‐binding (NB‐ARC), leucine‐rich‐repeat genes (NLRs) account for 60.8% of resistance (R) genes molecularly characterized from plants. NLRs exist as large gene families prone to tandem duplication and transposition, with high sequence diversity among crops and their wild relatives. This diversity can be a source of new disease resistance, but difficulty in distinguishing specific sequences from homologous gene family members hinders characterization of resistance for improving crop varieties. Current genome sequencing and assembly technologies, especially those using long‐read sequencing, are improving resolution of repeat‐rich genomic regions and clarifying locations of duplicated genes, such as NLRs. Using the conserved NB‐ARC domain as a model, 231 tentative NB‐ARC loci were identified in a highly contiguous genome assembly of sugar beet, revealing diverged and truncated NB‐ARC signatures as well as full‐length sequences. The NB‐ARC‐associated proteins contained NLR resistance gene domains, including TIR, CC and LRR, as well as other integrated domains. Phylogenetic relationships of partial and complete domains were determined, and patterns of physical clustering in the genome were evaluated. Comparison of sugar beet NB‐ARC domains to validated R‐genes from monocots and eudicots suggested extensive Beta vulgaris‐specific subfamily expansions. The NLR landscape in the rhizomania resistance conferring Rz region of Chromosome 3 was characterized, identifying 26 NLR‐like sequences spanning 20 MB. This work presents the first detailed view of NLR family composition in a member of the Caryophyllales, builds a foundation for additional disease resistance work in B. vulgaris, and demonstrates an additional nucleic‐acid‐based method for NLR prediction in non‐model plant species.     

The tomato I gene for Fusarium wilt resistance encodes an atypical leucine‐rich repeat receptor‐like protein whose function is nevertheless dependent on SOBIR1 and SERK3/BAK1

We have identified the tomato I gene for resistance to the Fusarium wilt fungus Fusarium oxysporum f. sp. lycopersici (Fol) and show that it encodes a membrane‐anchored leucine‐rich repeat receptor‐like protein (LRR‐RLP). Unlike most other LRR‐RLP genes involved in plant defence, the I gene is not a member of a gene cluster and contains introns in its coding sequence. The I gene encodes a loopout domain larger than those in most other LRR‐RLPs, with a distinct composition rich in serine and threonine residues. The I protein also lacks a basic cytosolic domain. Instead, this domain is rich in aromatic residues that could form a second transmembrane domain. The I protein recognises the Fol Avr1 effector protein, but, unlike many other LRR‐RLPs, recognition specificity is determined in the C‐terminal half of the protein by polymorphic amino acid residues in the LRRs just preceding the loopout domain and in the loopout domain itself. Despite these differences, we show that I/Avr1‐dependent necrosis in Nicotiana benthamiana depends on the LRR receptor‐like kinases (RLKs) SERK3/BAK1 and SOBIR1. Sequence comparisons revealed that the I protein and other LRR‐RLPs involved in plant defence all carry residues in their last LRR and C‐terminal LRR capping domain that are conserved with SERK3/BAK1‐interacting residues in the same relative positions in the LRR‐RLKs BRI1 and PSKR1. Tyrosine mutations of two of these conserved residues, Q922 and T925, abolished I/Avr1‐dependent necrosis in N. benthamiana, consistent with similar mutations in BRI1 and PSKR1 preventing their interaction with SERK3/BAK1.     

Characterization of disease resistance genes in the Brassica napus pangenome reveals significant structural variation

Methods based on single nucleotide polymorphism (SNP), copy number variation (CNV) and presence/absence variation (PAV) discovery provide a valuable resource to study gene structure and evolution. However, as a result of these structural variations, a single reference genome is unable to cover the entire gene content of a species. Therefore, pangenomics analysis is needed to ensure that the genomic diversity within a species is fully represented. Brassica napus is one of the most important oilseed crops in the world and exhibits variability in its resistance genes across different cultivars. Here, we characterized resistance gene distribution across 50 B. napus lines. We identified a total of 1749 resistance gene analogs (RGAs), of which 996 are core and 753 are variable, 368 of which are not present in the reference genome (cv. Darmor‐bzh). In addition, a total of 15 318 SNPs were predicted within 1030 of the RGAs. The results showed that core R‐genes harbour more SNPs than variable genes. More nucleotide binding site‐leucine‐rich repeat (NBS‐LRR) genes were located in clusters than as singletons, with variable genes more likely to be found in clusters. We identified 106 RGA candidates linked to blackleg resistance quantitative trait locus (QTL). This study provides a better understanding of resistance genes to target for genomics‐based improvement and improved disease resistance.     

Overexpression of AtPAD4 in transgenic Brachypodium distachyon enhances resistance to Puccinia brachypodii

• Brachypodium distachyon (L.) has recently emerged as a model for temperate grasses for investigating the molecular basis of plant–pathogen interactions. Phytoalexin deficient 4 (PAD4) plays a regulatory role in mediating expression of genes involved in plant defence.
• In this research, we generated transgenic B. distachyon plants constitutively overexpressing AtPAD4. Two transgenic B. distachyon lines were verified using PCR and GUS phenotype.
• Constitutive expression of AtPAD4 in B. distachyon enhanced resistance to Puccinia brachypodii. P. brachypodii generated less urediniospores on transgenic than on wild‐type plants. AtPAD4 overexpression enhanced salicylic acid (SA) levels in B. distachyon‐infected tissues. qRT‐PCR showed that expression of pathogenesis‐related 1 (PR1) and other defence‐related genes were up‐regulated in transformed B. distachyon following infection with P. brachypodii.
• Our results indicate that AtPAD4 overexpression in B. distachyon plants led to SA accumulation and induced PR gene expression that reduced the rate of colonisation by P. brachypodii.

     

Rye Pm8 and wheat Pm3 are orthologous genes and show evolutionary conservation of resistance function against powdery mildew

The improvement of wheat through breeding has relied strongly on the use of genetic material from related wild and domesticated grass species. The 1RS chromosome arm from rye was introgressed into wheat and crossed into many wheat lines, as it improves yield and fungal disease resistance. Pm8 is a powdery mildew resistance gene on 1RS which, after widespread agricultural cultivation, is now widely overcome by adapted mildew races. Here we show by homology‐based cloning and subsequent physical and genetic mapping that Pm8 is the rye orthologue of the Pm3 allelic series of mildew resistance genes in wheat. The cloned gene was functionally validated as Pm8 by transient, single‐cell expression analysis and stable transformation. Sequence analysis revealed a complex mosaic of ancient haplotypes among Pm3‐ and Pm8‐like genes from different members of the Triticeae. These results show that the two genes have evolved independently after the divergence of the species 7.5 million years ago and kept their function in mildew resistance. During this long time span the co‐evolving pathogens have not overcome these genes, which is in strong contrast to the breakdown of Pm8 resistance since its introduction into commercial wheat 70 years ago. Sequence comparison revealed that evolutionary pressure acted on the same subdomains and sequence features of the two orthologous genes. This suggests that they recognize directly or indirectly the same pathogen effectors that have been conserved in the powdery mildews of wheat and rye.     

Genetic dissection of a TIR‐NB‐LRR locus from the wild North American grapevine species Muscadinia rotundifolia identifies paralogous genes conferring resistance to major fungal and oomycete pathogens in cultivated grapevine

The most economically important diseases of grapevine cultivation worldwide are caused by the fungal pathogen powdery mildew (Erysiphe necator syn. Uncinula necator) and the oomycete pathogen downy mildew (Plasmopara viticola). Currently, grapegrowers rely heavily on the use of agrochemicals to minimize the potentially devastating impact of these pathogens on grape yield and quality. The wild North American grapevine species Muscadinia rotundifolia was recognized as early as 1889 to be resistant to both powdery and downy mildew. We have now mapped resistance to these two mildew pathogens in M. rotundifolia to a single locus on chromosome 12 that contains a family of seven TIR‐NB‐LRR genes. We further demonstrate that two highly homologous (86% amino acid identity) members of this gene family confer strong resistance to these unrelated pathogens following genetic transformation into susceptible Vitis vinifera winegrape cultivars. These two genes, designated r esistance to P lasmopara v iticola (MrRPV1) are the first resistance genes to be cloned from a grapevine species. Both MrRUN1 and MrRPV1 were found to confer resistance to multiple powdery and downy mildew isolates from France, North America and Australia; however, a single powdery mildew isolate collected from the south‐eastern region of North America, to which M. rotundifolia is native, was capable of breaking MrRUN1‐mediated resistance. Comparisons of gene organization and coding sequences between M. rotundifolia and the cultivated grapevine V. vinifera at the MrRUN1/MrRPV1 locus revealed a high level of synteny, suggesting that the TIR‐NB‐LRR genes at this locus share a common ancestor.     

Assembly and annotation of a draft genome sequence for Glycine latifolia,a perennial wild relative of soybean

Glycine latifolia (Benth.) Newell & Hymowitz (2n = 40), one of the 27 wild perennial relatives of soybean, possesses genetic diversity and agronomically favorable traits that are lacking in soybean. Here, we report the 939‐Mb draft genome assembly of G. latifolia (PI 559298) using exclusively linked‐reads sequenced from a single Chromium library. We organized scaffolds into 20 chromosome‐scale pseudomolecules utilizing two genetic maps and the Glycine max (L.) Merr. genome sequence. High copy numbers of putative 91‐bp centromere‐specific tandem repeats were observed in consecutive blocks within predicted pericentromeric regions on several pseudomolecules. No 92‐bp putative centromeric repeats, which are abundant in G. max, were detected in G. latifolia or Glycine tomentella. Annotation of the assembled genome and subsequent filtering yielded a high confidence gene set of 54 475 protein‐coding loci. In comparative analysis with five legume species, genes related to defense responses were significantly overrepresented in Glycine‐specific orthologous gene families. A total of 304 putative nucleotide‐binding site (NBS)‐leucine‐rich‐repeat (LRR) genes were identified in this genome assembly. Different from other legume species, we observed a scarcity of TIR‐NBS‐LRR genes in G. latifolia. The G. latifolia genome was also predicted to contain genes encoding 367 LRR‐receptor‐like kinases, a family of proteins involved in basal defense responses and responses to abiotic stress. The genome sequence and annotation of G. latifolia provides a valuable source of alternative alleles and novel genes to facilitate soybean improvement. This study also highlights the efficacy and cost‐effectiveness of the application of Chromium linked‐reads in diploid plant genome de novo assembly.     

Genetic and Molecular Mapping of Stripe Rust Resistance Genes in Wheat Cultivar Zhongliang 12

Stripe rust, caused by Puccinia striiformis f.sp. tritici (Pst), is one of the most widespread and destructive diseases of wheat worldwide. Resistance breeding is constantly pursued for decades to tackle the variations of prevalent Pst races. Zhongliang 12 has strong resistance to abiotic stresses, wide adaptability, higher resistance to stripe rust and excellent biological characteristics. To identify the resistance gene(s) against stripe rust, Zhongliang 12 was crossed with stripe rust susceptible genotype Mingxian 169, and F₁, F₂, F_2 : 3 and BC₁ progenies were tested with Chinese Pst race CYR30 and CYR31 in seedling stage in greenhouse. Zhongliang 12 possessed different dominant genes for resistance to each race. Linkage maps were constructed with four simple sequence repeats (SSRs) markers, Xwmc695, Xcfd20, Xbarc121 and Xbarc49, for the gene on wheat chromosome 7AL conferring resistance to CYR30 (temporarily designated as Yrzhong12‐1) with genetic distance ranging from 3.1 to 10.8 cM and four SSR markers, Xpsp3003, Xcfd2129, Xwmc673 and Xwmc51, for the gene on wheat chromosome 1AL conferring resistance to CYR31 (temporarily designated as Yrzhong12‐2) with genetic distance ranging from 3.9 cM to 9.3 cM. The molecular markers closely linked to each gene should be useful in marker‐assisted selection in breeding programmes for against stripe rust.