Pm21 from Haynaldia Villosa Encodes a CC-NBS-LRR Protein Conferring Powdery Mildew Resistance in Wheat
Liping Xing,Ping Hu,Jiaqian Liu,Kamil Witek,Shuang Zhou,Jiefei Xu,Weihao Zhou,Li Gao,Zhenpu Huang,Ruiqi Zhang,Xiue Wang,Peidu Chen,Haiyan Wang,Jonathan D. G. Jones,Miroslava Karafiatova,Jan Vrana,Jan Bartos,Jaroslav Dolezel,Yuanchun Tian,Yufeng Wu,Aizhong Cao
DOI: https://doi.org/10.1016/j.molp.2018.02.013
IF: 27.5
2018-01-01
Molecular Plant
Abstract:Wheat powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), is a destructive disease of wheat throughout the world. One of the most important environmental-friendly and economical methods to reduce wheat loss caused by Bgt is to develop highly resistant varieties (Kuraparthy et al., 2007Kuraparthy V. Chhuneja P. Dhaliwal H.S. Kaur S. Bowden R.L. Gill B.S. Characterization and mapping of cryptic alien introgression from Aegilops geniculata with new leaf rust and stripe rust resistance genes Lr57 and Yr40 in wheat.Theor. Appl. Genet. 2007; 114: 1379-1389Crossref PubMed Scopus (154) Google Scholar). Pm21 from the wild species Haynaldia villosa (also known as Dasypyrum villosum) confers high resistance to Bgt in wheat throughout all growth stages. It has now become one of the most highly effective genetic loci introgressed into wheat from wild species, and the commercial varieties harboring Pm21 have been widely used in wheat production with more than 4 million hectares in China. Owing to the recombination inhibition of alien chromosomes, the development of cytogenetic stocks has allowed the gradual narrowing of the physical location of the Pm21-containing chromosomal region and the identification of candidate genes. Development of the whole-arm translocation line T6VS·6AL led to the identification of the Pm21 locus on the 6VS arm (Chen et al., 1995Chen P.D. Qi L.L. Zhou B. Zhang S.Z. Liu D.J. Development and molecular cytogenetic analysis of wheat-Haynaldia villosa 6VS/6AL translocation lines specifying resistance to powdery mildew.Theor. Appl. Genet. 1995; 91: 1125-1128Crossref PubMed Scopus (223) Google Scholar), after which development of the resistant deletion stock del.6VS-1 (fraction length [FL] 0.58) and susceptible deletion stock del.6VS-2 (FL 0.45) subsequently resulted in the location of its chromosome region from FL 0.45 to FL 0.58 (Cao et al., 2011Cao A.Z. Xing L.P. Wang X.Y. Yang X.M. Wang W. Sun Y.L. Qian C. Ni J.N. Chen Y.P. Liu D.J. et al.Serine/threonine kinase gene Stpk-V, a key member of powdery mildew resistance gene Pm21, confers powdery mildew resistance in wheat.Proc. Natl. Acad. Sci. USA. 2011; 108: 7727-7732Crossref PubMed Scopus (249) Google Scholar). The development of two translocation lines, NAU418 and NAU419, which involves small segments of 6VS, then resulted in the definition of Pm21 to chromosomal region between molecular markers CINAU273 and CINAU276 (Chen et al., 2013Chen P.D. You C.F. Hu Y. Chen S.W. Zhou B. Cao A.Z. Wang X.E. Radiation-induced translocations with reduced Haynaldia villosa chromatin at the Pm21 locus for powdery mildew resistance in wheat.Mol. Breed. 2013; 31: 477-484Crossref Scopus (62) Google Scholar). We previously identified the Stpk-V gene as a key member of Pm21 gene locus; this gene was located in the 6VS FL 0.45-FL 0.58 region and also in the introgressed segment in NAU418 and NAU419. Overexpression of Stpk-V resulted in enhanced broad-spectrum resistance to powdery mildew, while Stpk-V silencing significantly but not completely compromised the resistance of T6VS·6AL (Cao et al., 2011Cao A.Z. Xing L.P. Wang X.Y. Yang X.M. Wang W. Sun Y.L. Qian C. Ni J.N. Chen Y.P. Liu D.J. et al.Serine/threonine kinase gene Stpk-V, a key member of powdery mildew resistance gene Pm21, confers powdery mildew resistance in wheat.Proc. Natl. Acad. Sci. USA. 2011; 108: 7727-7732Crossref PubMed Scopus (249) Google Scholar). The alien segment identified by genomic in situ hybridization (GISH) in the previous cytogenetic stocks was too large for mapping Pm21 gene; moreover, a lack of both reference data and loss-of-function mutants throughout the process of cloning this alien gene has limited its genetic elucidation. Therefore, whether a cluster of genes at the Pm21 gene locus contribute to the powdery mildew resistance remains unknown. In this study, comprehensive approaches including cytogenetic stock development, mutagenesis, resistance gene enrichment and sequencing (RenSeq)-PacBio, long-range assembly, and functional analysis were subsequently used to successfully clone the Pm21 gene. To precisely map the physical location of Pm21, we performed a long-term study involving cytogenetics approach, molecular marker analysis, and powdery mildew resistance evaluations to identify cryptic alien introgression lines among a large irradiated population of T6VS·6AL. NAU427 was identified as a useful stock; GISH using H. villosa genomic DNA as a probe could not identify the introgressed 6VS segment (Supplemental Figure 1A–1C). In addition, among all the 50 molecular markers located on the 6VS segments introgressed in NAU418, one previously developed molecular marker, 6EST258 (Supplemental Figure 1D), and 25 newly developed molecular markers were detected in NAU427 (Supplemental Figures 2 and 3). Two backcross populations were used to further validate the strong resistance effect of the introgressed segment in NAU427 (Supplemental Figure 4). This technique allowed us to narrow the Pm21 locus to a small physical region and provided effective material for precisely locating candidate genes. Most disease resistance (R) genes cloned in cereals are characterized as nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes (NLRs). Several genes on 6VS have been reported to be related to resistance against powdery mildew; however, none of these genes were NLRs, and silencing of these genes could only partially compromise the resistance of T6VS·6AL. Thus, we attempted to identify the expressed NLRs that were located in the introgressed segment of NAU427 and determined the relationship of these genes with Pm21. The RenSeq approach provided a rapid method for cloning full-length NLR genes from the whole genome (Sánchez-Martín et al., 2016Sánchez-Martín J. Steuernagel B. Ghosh S. Herren G. Hurni S. Adamski N. Vrána J. Kubaláková M. Krattinger S.G. Wicker T. et al.Rapid gene isolation in barley and wheat by mutant chromosome sequencing.Genome Biol. 2016; 17: 221Crossref PubMed Scopus (179) Google Scholar, Steuernagel et al., 2016Steuernagel B. Periyannan S.K. Hernández-Pinzón I. Witek K. Rouse M.N. Yu G. Hatta A. Ayliffe M. Bariana H. Jones J.D. et al.Rapid cloning of disease-resistance genes in plants using mutagenesis and sequence capture.Nat. Biotechnol. 2016; 34: 652-655Crossref PubMed Scopus (238) Google Scholar, Witek et al., 2016Witek K. Jupe F. Witek A.I. Baker D. Clark M.D. Jones J.D. Accelerated cloning of a potato late blight-resistance gene using RenSeq and SMRT sequencing.Nat. Biotechnol. 2016; 34: 656-660Crossref PubMed Scopus (151) Google Scholar). The RenSeq-PacBio results of the H. villosa genome data produced 107,153 reads of insert (ROIs) in 4.5 kb on average, then de novo-assembled NLR-containing ROIs resulted in the identification of 1509 contigs that harbored 482 unique full-length NLRs. To screen the expressed NLRs in H. villosa we annotated the de novo-assembled RNA-sequencing (RNA-seq) transcripts, and identified 457 expressed NLRs. In total, 45 individual NLR-carrying contigs were identified as the putative NLRs located on the 6V chromosome by collinearity analysis. Five of the expressed NLRs were mapped to the 6VS segment that was introgressed into NAU418; however, only two NLRs were located on the 6VS segment that was introgressed into NAU427 (Figure 1A). Therefore, two full-length NLR1-V and NLR2-V, both of which were located on the same RenSeq-PacBio Hv-contig 39 and identified as expressed NLRs in RNA-seq transcripts, were determined to be physically located at the Pm21 locus (Supplemental Figure 5A). The protein sequence similarity between NLR1-V and NLR2-V is 63.8% (Supplemental Figure 5B and 5C). To identify the putative candidate of Pm21, we performed a screen for powdery mildew-susceptible mutants of T6VS·6AL derived from ethyl methanesulfonate (EMS)-treated seeds, and identified eight individual mutants, namely SM-1 to SM-8 (Figure 1B). The sequences of NLR1-V and NLR2-V were compared between the T6VS·6AL and the mutant lines. Non-synonymous sequence changes were found in the NLR1-V coding sequence in all the mutants but not in the wild-type T6VS·6AL allele (Figure 1D and Supplemental Figure 6). Nevertheless, NLR2-V sequences were identical in the mutants and the wild-type plants. Therefore, NLR1-V was most likely a candidate gene of Pm21. The NLR1-V gene contains three exons and two introns. The genomic sequence of NLR1-V spans 3698 bp from the translation initiation (ATG) to the stop (TAA) codons, and its open reading frame (ORF) sequence is 2730 bp, corresponding to a predicted 909-amino-acid protein (Figure 1C and 1D). Domain analysis showed that NLR1-V contains the conserved coiled-coil (CC), NBS, and LRR domains belonging to a typical CC-NLR type R protein (Figure 1D). Phylogenetic analysis showed that NLR1-V was most closely related to the barley Rpg5 (Arora et al., 2013Arora D. Gross T. Brueggeman R. Allele characterization of genes required for rpg4-mediated wheat stem rust resistance identifies Rpg5 as the R gene.Phytopathology. 2013; 103: 1153-1161Crossref PubMed Scopus (21) Google Scholar), while distantly related to the broad-spectrum R proteins Pigm-R and Pigm-S from rice (Deng et al., 2017Deng Y.W. Zhai K.R. Xie Z. Yang D.Y. Zhu X.D. Liu J.Z. Wang X. Qin P. Yang Y. Zhang G. et al.Epigenetic regulation of antagonistic receptors confers rice blast resistance with yield balance.Science. 2017; 355: 962-965Crossref PubMed Scopus (313) Google Scholar) (Supplemental Figure 7). Single-cell transient overexpression system and Bgt infection assay were used to test the function of NLR1-V in wheat leaf epidermal cells. The Haustorium Index (HI) in SM-1 that co-transformed with the NLR1-V gene and the activity of β-glucuronidase (GUS) was significantly reduced compared with the HI in SM-1 transiently expressing the GUS gene, but it did not significantly differ from the HI in the T6VS·6AL lines overexpressing GUS (Figure 1E). These results indicate that NLR1-V could restore the resistance of SM-1 to the same level as could Pm21 and that, by preventing Bgt-induced haustoria formation, NLR1-V positively contributes to powdery mildew resistance. Barley stripe mosaic virus-induced gene silencing (BSMV-VIGS) targeting both the NBS domain and the LRR domain was used to further evaluate the function of NLR1-V gene in T6VS·6AL (Figure 1D and Supplemental Table 2). Both groups of the NLR1-V-silenced individuals were completely susceptible to Bgt; many powdery mildew colonies were clearly visible, and normal Bgt development was observed on NLR1-V-silenced leaves (Figure 1F and Supplemental Figure 8), indicating that silencing of NLR1-V could completely compromise the resistance of T6VS·6AL. Expression of NLR1-V was rapidly induced in the seedlings of resistant T6VS·6AL during Bgt infection (Supplemental Figure 9). The native promoter of NLR1-V was cloned and fused to the GUS gene, which was also strongly responsive to Bgt (Figure 1G). Stable transgenic plants of SM-1 carrying the NLR1-V gene driven by its native promoter were subsequently generated. Six primary independent transgenic plants displayed complete resistance phenotypes, as observed for the Pm21-harboring line T6VS·6AL (Figure 1H). To further reveal the physical structure of the Pm21-containing segment in NAU427, we sorted (Supplemental Figure 10), sequenced, assembled, and scaffolded the translocation chromosome of 6VS·6AL via Chicago long-range linkage by Dovetail Genomics (Thind et al., 2017Thind A.K. Wicker T. Šimková H. Fossati D. Moullet O. Brabant C. Vrána J. Doležel J. Krattinger S.G. Rapid cloning of genes in hexaploid wheat using cultivar-specific long-range chromosome assembly.Nat. Biotechnol. 2017; 35: 793-796Crossref PubMed Scopus (123) Google Scholar) (Supplemental Table 3). A 2.58-Mb region encompassing 19 ORFs was extracted by using the aforementioned 25 new co-segregating markers of the Pm21 gene (Figure 1I and Supplemental Figure 11). Micro-collinearity analysis revealed that the gene order in the Pm21 locus is relatively conserved in the analyzed grass species, whereas the diversity of NLR1-V and its orthologs among different species is higher than that of its flanking genes and their corresponding orthologs. A structural change with insertion leading to gene disruption in NLR1-V orthologs was observed in the subgenomes of common wheat (Supplemental Figure 12). Wild species constitute useful germplasms for identifying resistance genes to enhance the genetic diversity of wheat, whereas the lack of genomic references for wild species greatly limits the genetic elucidation of explored beneficial genes. In this study, comprehensive strategies were used to clone the Pm21 gene from H. villosa. Development of cryptic introgression lines represents an alternative strategy for targeting wild genes, and recently developed innovative sequencing technologies could be used to clone genes from wild species that lack genomic references. Previously, several genes have been reported as components involved in the Pm21-mediated resistance. NLR1-V is a typical CC-NLR resistance gene and provides high resistance to Bgt. Further studies are required to elucidate the underlying resistance mechanism and to reveal whether NLR1-V sufficiently contributes to broad-spectrum resistance, as does T6VS·6AL. This work was supported by the Important National Science and Technology Specific Projects of Transgenic Research (Grant No. 2018ZX0800963B), the Natural Science Foundation of China (Grant No. 31671685, 31771779), the Fundamental Research Funds for the National Central Universities (Grant No. KYZ201601), the Czech Science Foundation (Grant No. P501/12/G090), and the Ministry of Education, Youth and Sports of the Czech Republic (Grant No. LO1204 from the National Program of Sustainability I).