Identification of candidate genes for YrZM175
After anchoring the genomic interval of YrZM175 according to the Chinese Spring genome sequence we found that the target region contained a cluster of putative resistance genes. We cloned six putative disease resistance genes in this interval; four had no sequence difference between Zhongmai 175 and Avocet S or Mingxian 169; NLR gene TraesCS2A02G560600 had a SNP in exon 2 between Zhongmai 175 and Mingxian 169, and no PCR product for kinase gene TraesCS2A02G562800 was obtained from Zhongmai 175 and Avocet S or Mingxian 169. We performed transgenic assays of TraesCS2A02G560600 using two transgenic vectors; one was the CDS domain driven by Ubi promoter and another was the 6,251-bp ORF (including coding and intron region) with a 2,770-bp presumed native promoter and 1,500-bp terminator (3’-UTR). All the positive T1 plants were susceptible to race CYR32 in seedling stage tests, indicating that TraesCS2A02G560600 was not a gene candidate for YrZM175. We amplified five and two genes upstream and downstream, respectively, of TraesCS2A02G562800 according to the Chinese Spring genome and again failed to obtain PCR products, indicating great sequence differences of this genomic region between Zhongmai 175 and the Chinese Spring reference genome. With the release of reference genomes such as T. turgidum ssp. durum cv. Svevo (Maccaferri et al. 2019), T. aestivum cv. Zang1817 (Guo et al. 2020), 10 + Genome (Walkowiak et al. 2020), and T. aestivum cv. Fielder (Sato et al. 2021) we were able to extend our investigation. Since TraesCS2A02G562800 was the only undetermined putative resistance gene in the 636.4-kb target interval on chromosome arm 2AL of Zhongmai 175 according to the Chinese Spring genome, we performed genome walking to obtain the unknown sequence corresponding to TraesCS2A02G562800 in Zhongmai 175. After sequencing and blasting, we found that products corresponding to TraesCS2A02G562800 on 2AL in Zhongmai 175 were almost the same as a genome sequence (TRITD2Bv1G265970) on chromosome arm 2BL of durum wheat Svevo. Tests of relative TRITD2Bv1G265970 expression in infected leaf samples of Zhongmai 175 and Mingxian 169 at 0–7 dpi showed that relative expression of TRITD2Bv1G265970 in Zhongmai 175 was 3.8-fold higher than in Mingxian 169 at 1 dpi. However, the linkage distance between TRITD2Bv1G265970 and the resistance gene was 0.21 cM, indicating that TRITD2Bv1G265970 was not the gene candidate for YrZM175.
We amplified and sequenced ten putative disease resistance genes in the collinear chromosome arm 2AL interval in five Triticum genomes in Zhongmai 175 and obtained no candidate gene sequences. Finally, we assessed TRITD2Bv1G264480 on chromosome arm 2BL in T. turgidum as a candidate gene based on a co-segregating marker, sequence alignment of resistant and susceptible parents and gene expression patterns. Our research indicated that Zhongmai 175 has complicated genetic origins that caused difficulties in mapping the target genes according to single reference genome. In addition to the Chinese Spring genome sequence, it is also very necessary to refer other available genome sequences of wheat and related species for fine mapping and cloning of disease resistance genes in bread wheat.
Candidate gene TRITD2Bv1G264480
The putative resistance gene TRITD2Bv1G264480 encodes an NBS-LRR protein that belongs to two overlapping homologous superfamilies: P-loop containing nucleoside triphosphate hydrolase, and apoptotic protease-activating factors (helical domain) (http://www.ebi.ac.uk/interpro/entry/InterPro/IPR002182/). The predicted protein of TRITD2Bv1G264480 is closest to TRIDC2BG090350 on chromosome 2B from T. dicoccoides (Fig. S9). The paralogous genes on chromosome 2A in genomes of seven bread wheat varieties clustered together phylogenetically (Fig. S9), with 89% and 66% sequence similarities in the CDS domain to TRITD2Bv1G264480 on chromosome 2B and the paralogous gene TRITD2Av1G294930 on chromosome 2A in T. turgidum, respectively, indicating much closer phylogenetic relationship with TRITD2Bv1G264480 from chromosome 2B in T. turgidum, whereas the orthologous genes on chromosome 2B in T. aestivum cv. CDC Landmark, cv. CDC Stanley, cv. Fielder and cv. Julius showed 67.4% similarity with TRITD2Bv1G264480. In addition, paralogous genes on chromosomes 2A in T. Urartu and 2D in Ae. tauschii and orthologous genes on chromosome 2B in T. spelta also had quite distant phylogenetic relationships with TRITD2Bv1G264480 (Fig. S9). Therefore, we inferred TRITD2Bv1G264480 was transferred from chromosome 2B to 2A during evolution from an old crossover event of bread wheat. The CDS domain of candidate gene in Zhongmai 175 showed 97.8% and 96.8% similarities with TRITD2Bv1G264480 (T. turgidum) and TRIDC2BG090350 (T. dicoccoides), respectively, whereas the sequence similarities with homologues genes on chromosome 2A in other bread wheat varieties (Zang1817, Jagger, Fielder, Mace, CDC Landmark and Norin 61) were relatively low (all below 88.3%), indicating multiple allelic variations or haplotypes of the candidate gene present in bread wheat during evolution.
The relationship between YrZM175 and other Yr genes on chromosome 2A
Zhongmai 175 has been widely grown in the North Winter Wheat Zone of China since 2008. It was derived from the cross between Jing 411 and BPM27 (Fig. S10). BPM27 is a disease-resistant line developed by Prof. Zuomin Yang at China Agricultural University in the 1980s; one of its parents was VPM1, a French line derived from crosses of Ae. ventricosa, T. turgidum L. var. carthlicum (T. persicum) and common wheat cv. Marne (Fig. S10). The 2NS translocation fragment with Yr17 in Zhongmai 175 came from Ae. ventricosa.
Lu et al. (2016) mapped YrZM175 on chromosome arm 2AS using 344 F2 plants and 147 F2:3 lines derived from cross Lunxuan 987/Zhongmai 175 tested with CYR29. The present study used many more F2 plants from a cross of Avocet S and Zhongmai 175 but employed race CYR32. The discrepancy of gene locations can be attributed to different mapping populations, molecular markers and PST races used in the two studies. Unfortunately, the PST race CYR29 used previously (Lu et al. 2016) was no longer available to test the present population, thus we are not very sure whether the resistance genes are the same or not. The race CYR32 is virulent to VPM1 (Yr17) (Table S9). Screening of the 1,793 susceptible F2 plants in the Avocet S/Zhongmai 175 population with the marker URIC/LN2 (Helguera et al. 2003) as proxy for Yr17 indicated linkage of 55.2 cM. YrZM175 is located at the distal end (763,452,916–764,089,317 bp) of chromosome arm 2AL, whereas Yr17 is located at the distal (25.0–38.0 cM linkage interval) of chromosome arm 2AS (Cruz et al. 2016) (Fig. S11). Lu et al. (2016) performed multi-pathotype test on VPM1 (Yr17) and Zhongmai 175 (Yr17 and YrZM175), the two lines showed different reactions to 7 PST races at the seedling stage (Table S9). A test with the marker URIC/LN2 in 92 susceptible F2 plants indicated that the linkage distance between Yr17 and YrZM175 was 39.0 cM (Lu et al. 2016), being similar to the present study.
To date, three Yr genes have been mapped on chromosome arm 2AL, i.e., Yr1, Yr32 and YrJ22 (Lupton and Macer 1962; Eriksen et al. 2004; Chen et al. 2016). Yr1 is linked to Xgwm382 (Fig. 5B), with a genetic distance of 5.6 cM (Zheng et al. 2017; Bansal et al. 2009). Xgwm382 (at 772.96 Mb) was tested on all 1,793 susceptible F2 plants from Avocet S/Zhongmai 175; it was linked to YrZM175 (763,452,916–764,089,317 bp) with a genetic distance of 5.1 cM (Fig. 5C). The EST-SSR marker BU099658 (763.69 Mb) may accurately show the presence of Yr1 (Hasancebi et al. 2014), which locates at almost the same physical position as YrZM175. Because Zhongmai 175 and Chinese 166 have different pedigrees and different responses to PST races (Lu et al. 2016; Chen et al. 2016) (Table S9), we believe YrZM175 and Yr1 are different. Based on the linkage maps and physical positions of linked markers, the physical position of YrZM175 is close to Yr1, possibly indicating allelism. To clarify their relationship, tests of allelism and gene cloning should be carried out. Xgwm382 is also closely linked to YrJ22 with a genetic distance of 2.5 cM (Fig. 5D) (Chen et al. 2016). The physical position of YrJ22 is at 768.0–769.0 Mb on 2AL (Personal communication with Dr. Can Chen), a 5 Mb physical distance from YrZM175. As Jimai 22 and Zhongmai 175 have different pedigrees and different arrays of reactions to PST races (Table S9) we assume that YrZM175 and YrJ22 are different genes or alleles.
Eriksen et al. (2004) found that the marker Xcdo678 linked with Yr32 in a genetic distance of 35 cM, while Xcdo678 is co-segregating with Pm4 (Ma et al. 1994). Pm4 and Yr1 are closely linked to each other with a genetic distance of 2.0 cM (McIntosh and Arts 1996). As a result, the genetic distance between Yr32 and Yr1 is about 33 cM (Fig. 5A) (Yang et al. 2019), indicating that Yr32 is far away from YrZM175 and Yr1. Based on linkage maps and molecular marker analyses we conclude that YrZM175 is different from Yr1, Yr17, Yr32 and YrJ22.
Complexity Of The Distal Genome Region Of Chromosome 2al
Structural genomic variation such as presence-absence variations (PAVs) and copy number variations (CNVs) have been recognized to have the potential to generate phenotypic variation in maize (Springer et al. 2009) and barley (Muñoz-Amatriaín et al. 2013). Megabase-scale PAVs of Tripsacum origin were confirmed to be under selection during maize domestication and adaptation (Huang et al. 2021). Rimbert et al. (2018) also found high frequencies of PAVs in the distal regions of wheat chromosomes. Akhunov et al. (2003) found that duplicated loci were most frequently located in the distal regions of chromosomes, and their distribution was positively correlated with recombination rate. The higher recombination rate indicates that these regions are fast-evolving in adapting to biotic and abiotic stresses. Besides, there are higher polymorphisms in the distal regions of chromosomes with the decreased levels of synteny between homoeologous chromosomes as distance from centromeres increase (Akhunov et al. 2003). Our experimental results on chromosome arm 2AL in Zhongmai 175 were in agreement with these related reports mentioned above; this is reflected by the result that the candidate gene TRITD2Bv1G264480 for YrZM175 likely originated from chromosome arm 2BL of Triticum turgidum. The distal region of chromosome has great structural genomic variations, therefore is easier to evolve new alleles that trigger phenotypic variation.
In conclusion, YrZM175 is a stripe rust resistance gene that confers moderate to high resistance to stripe rust in the field. We fine mapped YrZM175 and identified a candidate gene based on a high-resolution linkage map, collinearity analysis and gene expression analysis. Transgenic assays are ongoing to validate the functions of gene candidate for YrZM175. The identification and fine mapping of YrZM175 provide options for deployment in combination with other effective resistance genes.