Skip to main content

ORIGINAL RESEARCH article

Front. Plant Sci., 12 December 2022
Sec. Plant Breeding

Identification of new resistance loci against wheat sharp eyespot through genome-wide association study

Xujiang Wu,,,&#x;Xujiang Wu1,2,3,4†Junchan Wang&#x;Junchan Wang1†Di WuDi Wu1Wei JiangWei Jiang1Zhifu GaoZhifu Gao1Dongsheng LiDongsheng Li1Rongling WuRongling Wu1Derong Gao,,,Derong Gao1,2,3,4Yong Zhang,,,*Yong Zhang1,2,3,4*
  • 1Institute of Agricultural Science of the Lixiahe District in Jiangsu Province, Yangzhou, China
  • 2Key Laboratory of Wheat Biology and Genetic Improvement on Low and Middle Yangtze River Valley Wheat Region, Ministry of Agriculture, Yangzhou, China
  • 3Jiangsu Co-Innovation Center for Modern Production Technology of Grain Crops, Yangzhou University, Yangzhou, China
  • 4Jiangsu Key Laboratory of Crop Genomics and Molecular Breeding, Yangzhou University, Yangzhou, China

Introduction: Wheat sharp eyespot caused by Rhizoctonia cerealis is a serious pathogenic disease affecting plants. The effective strategy for controlling this disease is breeding resistant cultivar. However, to date, no wheat varieties are fully resistant to sharp eyespot, and only a few quantitative trait loci (QTLs) have been shown to be associated with sharp eyespot resistance.

Methods: To understand the genetic basis of this disease, a genome-wide association study (GWAS) of sharp eyespot resistance in 262 varieties from all China wheat regions was conducted.

Results: After cultivation for three years, only 6.5% of the varieties were resistant to sharp eyespot. Notably, the varieties from the middle and lower Yangtze River displayed higher sharp eyespot resistance than those from Huanghuai wheat zone. Only two varieties had the same resistance level to the control Shanhongmai. The results of GWAS showed that 5 single nucleotide polymorphism (SNP) loci were markedly related to sharp eyespot resistance in the three years repeatedly, and two QTLs, qSE-6A and qSE-7B, on chromosome 6A and 7B were identified. Based on the ‘CG’ haplotypes of significant SNPs, we found that the two QTLs exhibited additive effects on attenuating sharp eyespot resistance.

Discussion: These results provide novel insights into the genetic basis of sharp eyespot resistance in China wheat varieties. The SNPs related to sharp eyespot resistance can be applied for marker-assisted selection in plant breeding.

Introduction

Bread wheat (Triticum aestivum L.) is an essential cereal crop cultivated worldwide, and is the main source of calories for the growth of global population (Gupta et al., 2020). Sharp eyespot leads to wheat yield losses, thereby posing serious threats to global food security. The extensive uses of nematicides and fumigants are not recommended because of environmental concerns. Breeding resistant varieties is considered the most effective and economical strategy for controlling plant diseases compared with pesticides and nematicides. To date, no wheat germplasm is fully resistant to sharp eyespot, and only a few varieties demonstrate partial resistance, including CI12633, Luck and AQ24788-83 (Chen et al., 2013; Wu et al., 2017).

The development of sharp eyespot disease was affected by the environmental conditions (Glynne, 1950). Several environmental factors are related to sharp eyespot, such as temperature, sowing time, soil-related factors, nutrient availability and agronomic practices (Hamada et al., 2011). Most scholars agreed that the high disease incidence in early sown fields can be attributed to the prolonged period of Rhizoctonia cerealis infection before winter. The favorable temperature of sharp eyespot is 16-28°C. If the temperature gets too high or too low, it will slow down the incidence of the disease. Some studies found that the disease severity was affected by soil fertilizers (Colbach et al., 1997). At first, many elements (the environments, the experiments repeat or the standard of disease sore) affected the phenotype of the disease, some researchers thought that the sharp eyespot was controlled by one or two major genes (Cai et al., 2006).

An accurate method for detecting quantitative traits loci (QTLs) in sharp eyespot resistance is a vital precondition for wheat sharp eyespot resistance studies. To date, the major methods for assessing sharp eyespot resistance include seeded inoculation, colonized wheat kernels on the soil surface (Cai et al., 2006; Wu et al., 2017), screening a disease nursery (Li et al., 2013), two-step colonized wheat kernels (Liu et al., 2020) and toothpick inoculation (Cai et al., 2006; Liu et al., 2020). In this study, we chose the spreading of colonized wheat kernels on the soil surface inoculation method, and moisturized by fog. Although this method greatly increases the repeatability, and is labor-intensive and time-consuming, it can increase the degree of accuracy for the phenotyping (Wu et al., 2017).

With the development types and quantity of molecule markers, most scientists discovered that sharp eyespot in wheat was a standard quantitative trait, which can be controlled by QTLs or multiple genes (Chen et al., 2013; Wu et al., 2017; Guo et al., 2017; Liu et al., 2021; Su et al., 2021). Ren Lijuan (Lijuan et al., 2004) first identified sharp eyespot resistance QTLs using SSR markers. Following that, approximately 30 QTLs of sharp eyespot resistance have been discovered on the wheat 21 chromosomes, except 2A, using recombinant inbred lines (RILs) mapping populations. In addition, Chen Jiang et al. (Chen et al., 2013) performed an association mapping with 549 simple-sequence repeat DNA markers using 241 RILs, and identified seven QTLs significantly contributed to sharp eyespot resistance. Wu Xujiang et al. (Wu et al., 2017) detected five QTLs significantly associated with sharp eyespot resistance through association analysis with the Illumina iSelect 90K SNP wheat chip and 101 SSR markers using 224 RILs. Liu Caiyun et al. (Liu et al., 2021) identified two stable QTL and three novel stable QTL through association analysis with the polymorphic Affymetrix Axiom™ 55K SNPs using 215 F8 lines from Niavt14/Xuzhou25. Over a period of twenty years, the discovered resistance genes were about thirty. However, only a few of these QTLs can be applied for sharp eyespot resistance. Moreover, no QTL for sharp eyespot resistance is fine-mapped and has been cloned so far by the map-based cloning method.

At present, all the mapping population of sharp eyespot research is RILs. However, this approach is time consuming, as the QTL mapping may require a large biparental population. Genome-wide association studies (GWAS) have become a crucial method to detect putative genes for complex traits in plants. Compared to the method in sharp eyespot study, GWAS employs high-density single nucleotide polymorphism (SNP) markers and more diverse natural populations, which can identify QTLs much closer to the putative genes and determine favorable haplotypes of sharp eyespot resistance among tested varieties.

In this research, 262 wheat cultivars were genotyped using the 90K SNP array, involving a large number of wheat cultivars from all China wheat regions. According to the multi-environmental trail data, GWAS was employed to identify favorable alleles of sharp eyespot. The results provide new insights into marker-assisted selection in plant breeding.

Materials and methods

Wheat materials

An association panel comprising 262 wheat cultivars was used for SNP genotyping and 3-year (2018, 2019 and 2020) sharp eyespot resistance phenotyping (Supplemental Table 1). Of them, 43 cultivars were collected from the middle and lower reaches of Yangtze River, 21 from Huanghuai wheat zone, 5 cultivars were derived from foreign regions, and 7 cultivars were obtained from Shanxi provinces in the Northern China Winter Wheat Region and 15 cultivars from Southern China. Six wheat varieties with known sharp eyespot resistance levels, CI12633 (moderately resistant, the disease index was between 44.8-50.27), shanhongmai (moderately resistant), Yangmai158 (highly susceptible), and Yangmai9 (high susceptible) were served as references in the analysis. The experiments were carried out at Wantou Experimental Station, Institute of Agricultural Sciences of the Lixiahe, Yangzhou, Jiangsu Province, China (latitude 32.24°N, altitude 8 m, annual precipitation approximately 1000 mm, from early November to the next May) for three years.

Sharp eyespot fungal inoculation

The R. cerealis isolate R0301, a predominant strain in Jiangsu province, China, was obtained from Prof. Shibin Cai and Prof. Huaigu Chen (Jiangsu Academy of Agricultural Sciences, China). Before use, R0301 was inoculated on potato-dextrose agar (PDA). Inoculum was generated on the sterilized wheat kernels as described previously (Chen et al., 2013).

Randomized complete block designed with 3 replicates/environment was used in this study. The cultivars (50 seeds/row) were sown in one row of 80*40 cm for each replication. Some seedlings were removed to obtain a density of 40 cultivars/row. The inoculum method was the same with Wu Xujiang et al. (Wu et al., 2017). In early March, when the weather became warm and wet, the plants were inoculated by placing R0301-colonized wheat kernels on the soil surface close to the plants, but not touching them. The kernels were then covered with earth and water was sprinkled three times a day for the first month in non-raining days. Overall, the weather was very suitable for sharp eyespot outbreak two of three seasons with consistent rainfall occurring throughout the experimental period until the final disease recording.

Phenotyping

In spring next year, when the susceptible mocks Yangmai 9 began to succumb to R. cerealis infection, all sample plants except the edge plants were dug up and individually assessed for disease (approximately 200-250 wheat stems of each line). The infection types were categorized qualitatively from 0 to 5: 0, no lesion; 1, the lesion appeared on the sheaths rather than stems; 2, the width of the lesion was <50% of the infected stem perimeter; 3, the width of the lesion is >50 and <75% of the infected stem perimeter; 4, the width of the disease lesion is >75% of the infected stem perimeter; 5, white spike or dead plant. Disease index =((0×X0+1×X1+2×X2+3×X3+4×X4+5×X5)/[(X0+X1+X2+X3 +X4+X5)×5]) ×100, where X0-X5 indicated plants with infection types 0-5, respectively (Wu et al., 2017).

DNA extraction and genotyping

A modified CTAB method described in the CIMMYT laboratory protocol was used to extract genomic DNA from fresh leaf samples (Clarke, 2009). The RNA yield was assessed using a NanoDrop 8000 spectrophotometer. The Illumina iSelect 90K SNP technology (Chinese academy of agricultural sciences, Beijing) was used for genotyping. SNPs with minor allele frequency (MAF) <0.10 and/or with >10% of missing data were eliminated. The SNPs were mapped on wheat chromosomes in accordance with Chinese Spring reference genome sequences (International Wheat Genome Sequencing Consortium (IWGSC), 2018).

Population structure and linkage disequilibrium analyses

PS was assessed by Bayesian cluster analysis using the Structure 2.3.4 software used 26517 SNPs. There were 34039 SNPs distributed on all 21 wheat chromosomes with r2<0.2 (Najafi et al., 2019). LD was conducted with TASSEL 5.0 software using the expected versus observed allele frequencies of the SNPs (Liu et al., 2017).

GWAS analysis

The identification of SNPs related to sharp eyespot resistance was performed as described previously (Zhou et al., 2016). The GWAS panel consisted of 90k SNP markers (Supplemental Table 2), which was analyzed based on the mixed linear model (MLM) and TASSEL 5.0 software (Bradbury et al., 2007). In the MLM, the kinship matrix is used jointly with PS (Henderson, 1975). To reduce type I error, the region with SNPs (P<1*10-4) in a 200-kb genomic window was selected for further analyses. PerL was used to construct Manhattan maps (Christiansen et al., 2012). A standard linear mixed model was fitted by EMMAX (Kang et al., 2010), and then added with 3 principal component covariates. The Manhattan and quantile-quantile plots were generated using the R package (https://cran.rproject.org/web/pachages/qqman/).

Statistical analysis

Statistical analyses, including ANOVA, heritability, correlation, and descriptive statistics, were conducted using the R statistical package (R Core Team, 2019). The broad-sense heritability for the traits was estimated by the formula H2=VG/(VG+VE) where VG and VE represent the estimates of genetic and environmental variance, respectively (Arora et al., 2017).

Results

Phenotypic variation, correlation, and heritability of sharp eyespot

Phenotypic variations of sharp eyespot among the 262 wheat accessions were evaluated at the experimental field in Yangzhou during the winter seasons of 2018-2020 in MLM. Phenotypic evaluation showed the presence of sharp eyespot resistant in three years. The mean of sharp eyespot disease index ranged from 25.86% to 95.23% (Supplemental Table 3), with a normal distribution pattern. The standard deviation and coefficient of variation (CV) were 13.5 and 4.61%, respectively. The frequency of sharp eyespot disease index was normally distributed in each year (Figure 1), implying the nature of the disease.

FIGURE 1
www.frontiersin.org

Figure 1 Frequency distributions of the sharp eyespot disease index in 2018, 2019 and 2020.

Significant positive correlations were observed among the three years for sharp eyespot disease index (r=0.75-0.95), suggesting that the traits occur consistently (Table 1). The severity of the disease index was most high in 2018 growing season (Supplemental Figure 1). Wheat cultivars from different wheat regions in China had varying degree of resistance to sharp eyespot. The cultivars from Jiangsu province were consistently resistant to sharp eyespot in three seasons, while those from Henan province had the greatest vulnerability. The ANOVA results revealed significant differences among the genotypes for sharp eyespots resistance (Table 2). High broad-sense heritability was observed for the trait, with 0.86, indicating that the trait was largely influenced by genetic factors.

TABLE 1
www.frontiersin.org

Table 1 Correlation coefficients among the three trial environments for the sharp eyespot disease index.

TABLE 2
www.frontiersin.org

Table 2 The ANOVA of wheat sharp eyespot disease scores of the cultivars.

The cultivars from different wheat regions had varying degree of resistance to sharp eyespot (Table 3). Those from Huanghuai region and foreign were consistently resistant to sharp eyespot in the three seasons, while from the wheat region in the middle and lower reaches of Yangtze River exhibited the greatest vulnerability.

TABLE 3
www.frontiersin.org

Table 3 Sharp eyespot disease scores of different wheat regions.

Population structure and LD analyses

PC analysis is compulsory in many populations as the presence of a large number of genotypes in the research can result in a higher quantity of false associations between the phenotypes and unlinked markers. A total of 34039 SNPs (Supplemental Table 4), of which 19975 were mapped to a relative position on a single chromosome. Among them, 6751, 8248 and 4587 were found in A, B and D genome chromosome, respectively, but 389 SNPs could not be presented in genome chromosome. From the kinship analysis, a break in the slope was observed at K = 2, followed by a flattening of the curve (Supplemental Figure 2A); hence, it was revealed that the most likely number of sub-populations was two (K = 2; Supplemental Figure 2B). In addition, the results were verified by principle component analysis (PCA), on account of the standardized covariance of genetic distances of SNP markers (Supplemental Figure 2C).

For LD analysis, only SNPs with a known position and a MAF >5% were selected. The LD half-width (the point at which LD decayed to 48% of the peak) was 10.5 cM. The r2 values of A, B and D sub-genomes gradually declined with increasing genetic distance (Supplemental Figure 3).

Identification of new QTLs for sharp eyespot resistance

According to the 34039 SNP dataset and the disease index of the selected varieties, 5, 19 and 14 SNPs were markedly related to sharp eyespot resistance in 2018, 2019 and 2020, respectively. The highest number of SNPs in 2019 may be attributed to a large phenotypic variation. The contribution of each SNP on phenotypic variance ranged from 2.7% to 8.5%. These SNPs were distributed mainly on chromosomes 6A, 7A and 7B. Among them, two SNP markers were repeatedly observed in all the three experiments, and three SNPs were detected in two years (Supplemental Table 5). These five SNPs were primarily located on chromosome 6A. Due to the high degree of LD in the wheat genome, the SNP clusters identified on chromosome 6A (333.95-351.08 Mb) might represent the chromosome regions with sharp eyespot loci. Although many QTLs were detected, only two regions were found to be large effective QTL for sharp eyespot resistance, which were designated as qSE-6A and qSE-7B, respectively (Figure 2).

FIGURE 2
www.frontiersin.org

Figure 2 GWAS analysis of SNP loci related to sharp eyespot resistance in the three years. (A–D), Manhattan plots of SNPs related to sharp eyespot resistance on 21 wheat chromosomes in three years and mean of three years. The x-axis and y-axis represent the genomic coordinates and LOD score for each SNP, respectively. (E–H), Quantile-quantile (Q-Q) plots for the GWAS data in 2018 (E), 2019 (F), 2020 (G) and mean of three years (H).

The additive effect of SNP haplotypes related to QTLs qSE-6A and qSE-7B related SNP haplotypes on sharp eyespot resistance

To assess the effect of QTL on sharp eyespot resistance, the most significantly associated SNPs were used to distinguish varieties with/without the resistant alleles. For qSE-6A, BS00078658_51 was the most significant SNP locus, in which ‘C’ and ‘T’ were the favorable and unfavorable alleles, respectively. The LOD score of each dot represents a transformed P value, -lg(P). A total of 222 varieties containing ‘C’ and 40 varieties containing ‘T’ were designated C-type and T-type, respectively. From the three years results, we found that the mean disease score of C-type varieties was remarkably decreased compared to the T-type varieties, implying the significant effect of qSE-6A on sharp eyespot resistance (Figure 3A). For qSE-7B, wsnp_Ex_c3501_6407527 was the most significant SNP locus, in which the favorable and unfavorable haplotypes were ‘G’ and ‘A’ halpotypes, respectively. In total, 114 and 148 varieties harbored the resistant ‘G’ and susceptible ‘A’ haplotypes, respectively. The mean disease score of ‘G’ haplotype variety was markedly reduced compared to ‘A’ haplotype variety (Figure 3B). In addition, the mean disease score of ‘CG’ haplotype variety (n = 97) harboring favorable alleles in both qSE-6A and qSE-7B loci was even significantly lower than that of ‘TA’ haplotype variety (n = 22) harboring unfavorable alleles. Based on the difference in mean diseases score, the resistance effects of qSE-6A, qSE-7B and qSE-6A&qSE-7B were determined to be 59.47, 57.22 and 56.82, respectively, and without the two genes was 63.17, indicating the additive effects of the two QTLs (Figure 3C). These findings demonstrate that both qSE-6A and qSE-7B have significant pyramiding effects on wheat sharp eyespot resistance.

FIGURE 3
www.frontiersin.org

Figure 3 Resistance effects and ‘CG’ haplotype distribution of qSE-6A and qSE-7B. (A), Effects of the favorable allele ‘C’ of SNP BS00078658_51 (represent qSE-6A) versus the unfavorable allele ‘T’. (B), Effects of the favorable allele ‘G’ of SNP wsnp_Ex_c3501_6407527 (represent qSE-7B) versus the unfavorable allele ‘A’. (C), Effects of the favorable haplotype ‘CG’ of SNP BS00078658_51 and wsnp_Ex_c3501_6407527 compared to the unfavorable haplotype ‘TA’. (D) the ratio of ‘CG’ haplotype varieties in each wheat region.

Furthermore, the varieties were mainly located in the middle and lower reaches of the Yangtze River, Huanghuai and southwest wheat growing regions. Only a few were introduced from abroad, which harbored the favorable haplotypes “CG”. The reaction to wheat sharp eyespot in all wheat region was close, except for those from abroad. In addition, all known resistance varieties, such as Shanhongmai, CI12633 and Limai16, belonged to the “CG” haplotypes. Hence, the associated SNPs in the qSE-6A and qSE-7B regions can be used for marker-assisted breeding program to increase the resistance to wheat sharp eyespot (Figure 3D).

Discussion

The genetic architecture of wheat sharp eyespot resistance remains largely unknown due to the lack of appropriate mapping populations and reliable disease evaluation methods. Here, the wheat sharp eyespot resistance of 262 varieties was evaluated in three years, and the results showed that the overall resistance of the varieties was considerably less, and as low as 3% of these varieties were classified as the resistant. The results of GWAS analysis showed that 15 SNP loci were markedly related to wheat sharp eyespot resistance, and two SNP markers on chromosomes 6A and 7B were repeatedly observed in the three years. Two QTLs were designated and named as qSE-6A and qSE-7B, and the ‘CG’ haplotype varieties could decrease sharp eyespot disease index by about 20% and 10%, respectively. The QTLs of sharp eyespot resistance on chromosome 6A and 7B were mapped to the wheat reference Chinese spring genome based on their flanking or associated markers. The location region of qSE-6A (physical location was 351086145bp) was not overlapped with previous QTLs (physical location was 77690500bp) on chromosome 6A, implying that the qSE-6A is a novel sharp eyespot resistance locus.

No studies have been reported on the use of GWAS with high-density SNP markers for identifying loci related to sharp eyespot resistance. To further verify qSE-6A and qSE-7B, some LD analysis-related populations should be established via a backcrossing or hybridization strategy.

One of the advantages of GWAS is that it can be used to analyze the target varieties with favorable/unfavorable alleles, which is unattainable by the LD analysis of two parental alleles. Hence, GWAS data can be successfully adopted in the breeding practices. Based on the results of SNP markers related to sharp eyespot resistance, the suitable parents can be selected, and we also can attain rapid transfer of the target loci by performing genomic breeding and marker-assisted selection. However, it may impractical to enhance the overall resistance of wheat to sharp eyespot through marker-assisted polygenic polymerization, as this disease is controlled by multiple genes. Thus, it is of imperative importance to alter the frequencies of wheat resistant to sharp eyespot during conventional breeding. In this research, we observed that some abroad imported varieties and landraces were significantly more resistant than other region varieties, which may be attributed to the fact that the majority of abroad varieties and landraces harbored “CG” haplotypes in qSE-6A and qSE-7B regions.

In this study, the two QTLs, qSE-6A and qSE-7B, were delimited into ≤0.8 Mb regions with only 15 annotated high-confidence SNPs. Thus, the putative genes in some of the QTL regions can be identified based on the wheat reference genome annotation (International Wheat Genome Sequencing Consortium (IWGSC), 2018). The disease resistance mechanism in plants is relatively complex, which involves the activation of disease resistance genes, pathogenesis-related genes, transcription factors, mitogen-activated protein kinase pathways, and recognition of pathogen-associated patterns or effectors. Thus far, 9 resistance genes have been cloned for wheat resistance to sharp eyespot, and the types of resistance genes are diversified, such as COMT (Wang et al., 2018), pathogen-induced ERF1 (Wei et al., 2016), TaAGC1 (Zhu et al., 2015), RcCUT1 (Lu et al., 2018), TaCAD12 (Rong et al., 2016), PP2A (Zhu et al., 2018), TaCPK7-D (Wei et al., 2016), TaRCR1 (Zhu et al., 2017), TaRIM1 (Shan et al., 2016), and DmAMP1W (Su et al., 2020). These genes all contributed to sharp eyespot resistance and can serve as potential QTLs for plant disease resistance.

In the QTLs associated with sharp eyespot resistance, 193 genes were annotated, and 5 of them were resistant genes from 2 QTLs (Supplemental Table 6). All of the five genes were detected in three years, implying that these genes can be used for identifying causal variants in these QTLs. Based on the transcriptomic data from our previous research, a putative gene that activated by sharp eyespot pathogen was identified in this study. qSE-7B on chromosome 7B was significant under the three environments. This QTL region contains a carboxypeptidase gene, TraesCS7B01G054900.1, which participates in lignin synthesis. In the future, CRISPR-Cas9-based gene editing will be conducted to validate these candidate genes.

In this study, GWAS analysis was performed to assess sharp eyespot resistance based on 34039 SNPs in 262 wheat accessions. Two QTLs were correlated with sharp eyespot in three years. The two QTLs were stably detected in the three years, and the QTL in chromosome 6A was newly discovered. Functional analysis revealed that 8 core genes were associated with sharp eyespot resistance. The robustness of multilocus GWAS approach in studying wheat sharp eyespot resistance was evaluated. This study allows us to better understand the mechanism of wheat sharp eyespot resistance. The identified QTLs and putative SNPs can be employed to characterize genes related to wheat sharp eyespot resistance and further improve this disease through marker-assisted breeding.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.

Author contributions

XW and JW: Conceptualization, Supervision, Resources, writing review and editing and Funding acquisition. YZ: Supervision, Resources and Funding acquisition. DW and WJ performed the experiments; ZG and DL provided additional resources and suggestions; RW and DG provided guidance throughout the entire study. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by Natural Science Foundation of Jiangsu Province (BK20211117) and Jiangsu Agriculture Science and Technology Innovation fund (CX(20)3009), and the China Agricultural Research System (CARS-03-8).

Acknowledgments

The authors would like to express their gratitude to EditSprings (https://www.editsprings.cn) for the expert linguistic services provided. The authors acknowledge the technical assistance of Hao Chenyang coming from CAAS.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.1056935/full#supplementary-material

Supplementary Table 2 | The information of 90K SNP

Supplementary Table 3 | Phenotypic variation for sharp eyespt in three years (2018-2020)

Supplementary Table 4 | The raw analysis data

Supplementary Table 6 | Annotated genes within the identified QTL region

References

Arora, S., Singh, N., Kaur, S., Bains, N. S., Uauy, C., Poland, J., et al. (2017). Genome-wide association study of grain architecture in wild wheat aegilops tauschii. Front. Plant Sci. 8, 886. doi: 10.3389/fpls.2017.00886

PubMed Abstract | CrossRef Full Text | Google Scholar

Bradbury, P. J., Zhang, Z., Kroon, D. E., Casstevens, T. M., Ramdoss, Y., Buckler, E. S. (2007). TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics 23, 2633–2635. doi: 10.1093/bioinformatics/btm308

PubMed Abstract | CrossRef Full Text | Google Scholar

Cai, S-b., Ren, L-j., Yan, W., Wu, J.-Z., Chen, H.-G., Wu, X.-Y., et al. (2006). Germplasm development and mapping of resistance to sharp eyespot (Rhizoctoniacerealis) in wheat. Scientia Agricult. Sin. 39, 928–934.

Google Scholar

Chen, J., Li, G. H., Du, Z. Y., Quan, W., Zhang, H. Y., Che, M. Z., et al. (2013). Mapping of QTL conferring resistance to sharp eyespot (Rhizoctonia cerealis) in bread wheat at the adult plant growth stage. Theor. Appl. Genet. 126, 2865–2878. doi: 10.1007/s00122-013-2178-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Christiansen, T., Foy, B. D., Wall, L., Orwant, J. (2012). Programming Perl: Unmatched power for text processing and scripting (Sebastopol, CA, USA: O’Reilly).

Google Scholar

Clarke, J. D. (2009). Cetyltrimethyl ammonium bromide (CTAB) DNA miniprep for plant DNA isolation. Cold Spring Harb. Protoc. 2009, 5177. doi: 10.1101/pdb.prot5177

CrossRef Full Text | Google Scholar

Colbach, N., Lucas, P., Cavelier, N., Cavelier, A. (1997). Influence of cropping system on sharp eyespot in winter wheat. Crop Prot. 16, 415–422. doi: 10.1016/S0261-2194(97)00018-5

CrossRef Full Text | Google Scholar

Glynne, M. D. (1950). Sharp eyespot as a severe disease of oats. Nature 166, 232. doi: 10.1038/166232a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, Y., Du, Z., Chen, J., Zhang, Z. (2017). QTL mapping of wheat plant architectural characteristics and their genetic relationship with seven QTLs conferring resistance to sheath blight. PLos One 12, e0174939. doi: 10.1371/journal.pone.0174939

PubMed Abstract | CrossRef Full Text | Google Scholar

Gupta, P. K., Balyan, H. S., Sharma, S., Kumar, R. (2020). Genetics of yield, abiotic stress tolerance and biofortification in wheat (Triticum aestivum l.). Theor. Appl. Genet. 133, 1569–1602. doi: 10.1007/s00122-020-03583-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Hamada, M. S., Yin, Y., Chen, H., Ma, Z. (2011). The escalating threat of rhizoctonia cerealis, the causal agent of sharp eyespot in wheat. Pest Manag. Sci. 67, 1411–1419. doi: 10.1002/ps.2236

PubMed Abstract | CrossRef Full Text | Google Scholar

Henderson, C. R. (1975). Best linear unbiased estimation and prediction under a selection model. Biometrics 31, 423–447. doi: 10.2307/2529430

PubMed Abstract | CrossRef Full Text | Google Scholar

International Wheat Genome Sequencing Consortium (IWGSC) (2018). Shifting the limits in wheat research and breeding using a fully annotated reference genome. Science 361 (6403), eaar7191. doi: 10.1126/science.aar7191

PubMed Abstract | CrossRef Full Text | Google Scholar

Kang, H. M., Sul, J. H., Service, S. K., Zaitlen, N. A., Kong, S-y., Freimer, N. B., et al. (2010). Variance component model to account for sample structure in genome-wide association studies. Nat. Genet. 42, 348–354. doi: 10.1038/ng.548

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, J., Anderson, N. P., Mundt, C. C. (2020). Methods for screening wheat genotypes for resistance to sharp eyespot in the field and greenhouse. Plant Dis. 104(12), 3192–3196. doi: 10.1094/PDIS-12-19-2652-RE

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C., Guo, W., Zhang, Q., Fu, B., Yang, Z., Sukumaran, S., et al. (2021). Genetic dissection of adult plant resistance to sharp eyespot using an updated genetic map of Niavt14 x Xuzhou25 winter wheat recombinant inbred line population. Plant Dis. 105 (4), 997–1005. doi: 10.1094/PDIS-09-20-1924-RE

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, J., He, Z., Rasheed, A., Wen, W., Yan, J., Zhang, P., et al. (2017). Genome-wide association mapping of black point reaction in common wheat (Triticum aestivum l.). BMC Plant Biol. 17 (1), 220. doi: 10.1186/s12870-017-1167-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., Wang, X., Chen, H., Li, W., Zhang, H. (2013). Reaction of wheat-progenies and wheat germplasm to sharp eyespot. Acta Agronomica Sin. 39, 999. doi: 10.3724/SP.J.1006.2013.00999

CrossRef Full Text | Google Scholar

Lu, L., Rong, W., Massart, S., Zhang, Z. (2018). Genome-wide identification and expression analysis of cutinase gene family in rhizoctonia cerealis and functional study of an active cutinase RcCUT1 in the fungal-wheat interaction. Front. Microbiol. 9, 1813. doi: 10.3389/fmicb.2018.01813

PubMed Abstract | CrossRef Full Text | Google Scholar

Najafi, A., Janghorbani, S., Motahari, S. A., Fatemizadeh, E. (2019). “Statistical association mapping of population-structured genetic data,” in IEEE/ACM Trans Comput Biol Bioinform, 16 (2), 638–649.

Google Scholar

R Core Team (2019). R: A language and environment for statistical computing (Vienna, Austria: R Foundation for Statistical Computing).

Google Scholar

Ren, L. J., Cai, S. B., Tang, T., Zhong, W. J., Zhou, M. P., Yan, W., et al. (2004). SSR markers linked resistance QTLs to sharp eyespot (Rhizoctonia cerealis) in wheat. J. YangZhou Univ. (Agricult. Life Sci. Edition) 25, 16–19.

Google Scholar

Rong, W., Luo, M., Shan, T., Wei, X., Du, L., Xu, H., et al. (2016). A wheat cinnamyl alcohol dehydrogenase TaCAD12 contributes to host resistance to the sharp eyespot disease. Front. Plant Sci. 7, 1723. doi: 10.3389/fpls.2016.01723

PubMed Abstract | CrossRef Full Text | Google Scholar

Shan, T., Rong, W., Xu, H., Du, L., Liu, X., Zhang, Z. (2016). The wheat R2R3-MYB transcription factor TaRIM1 participates in resistance response against the pathogen rhizoctonia cerealis infection through regulating defense genes. Sci. Rep. 6, 28777. doi: 10.1038/srep28777

PubMed Abstract | CrossRef Full Text | Google Scholar

Su, Q., Wang, K., Zhang, Z. (2020). Ecotopic expression of the antimicrobial peptide DmAMP1W improves resistance of transgenic wheat to two diseases: Sharp eyespot and common root rot. Int. J. Mol. Sci. 21 (12), 647. doi: 10.3390/ijms21020647

PubMed Abstract | CrossRef Full Text | Google Scholar

Su, J., Zhao, J., Zhao, S., Li, M., Pang, S., Kang, Z., et al. (2021). Genetics of resistance to common root rot (Spot blotch), fusarium crown rot, and sharp eyespot in wheat. Front. Genet. 12, 699342. doi: 10.3389/fgene.2021.699342

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, M., Zhu, X., Wang, K., Lu, C., Luo, M., Shan, T., et al. (2018). A wheat caffeic acid 3-o-methyltransferase TaCOMT-3D positively contributes to both resistance to sharp eyespot disease and stem mechanical strength. Sci. Rep. 8 (1), 6543. doi: 10.1038/s41598-018-24884-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, X., Shen, F., Hong, Y., Rong, W., Du, L., Liu, X., et al. (2016). The wheat calcium-dependent protein kinase TaCPK7-d positively regulates host resistance to sharp eyespot disease. Mol. Plant Pathol. 17 (8), 1252–1264. doi: 10.1111/mpp.12360

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, X., Cheng, K., Zhao, R., Zang, S., Bie, T., Jiang, Z., et al. (2017). Quantitative trait loci responsible for sharp eyespot resistance in common wheat CI12633. Sci. Rep. 7, 11799. doi: 10.1038/s41598-017-12197-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, L., Liu, S., Wu, W., Chen, D., Zhan, X., Zhu, A., et al. (2016). Dissection of genetic architecture of rice plant height and heading date by multiple-strategy-based association studies. Sci. Rep. 6, 29718. doi: 10.1038/srep29718

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, X., Lu, C., Du, L., Ye, X., Liu, X., Coules, A., et al. (2017). The wheat NB-LRR gene TaRCR1 is required for host defence response to the necrotrophic fungal pathogen rhizoctonia cerealis. Plant Biotechnol. J. 15 (6), 674–687. doi: 10.1111/pbi.12665

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, X., Wang, Y., Su, Z., Lv, L., Zhang, Z. (2018). Silencing of the wheat protein phosphatase 2A catalytic subunit TaPP2Ac enhances host resistance to the necrotrophic pathogen rhizoctonia cerealis. Front. Plant Sci. 9, 1437. doi: 10.3389/fpls.2018.01437

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, X., Yang, K., Wei, X., Zhang, Q., Rong, W., Du, L., et al. (2015). The wheat AGC kinase TaAGC1 is a positive contributor to host resistance to the necrotrophic pathogen rhizoctonia cerealis. J. Exp. Bot. 66 (21), 6591–6603. doi: 10.1093/jxb/erv367

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: wheat, sharp eyespot, genome-wide association study (GWAS), quantitative trait loci (QTL), single nucleotide polymorphism (SNP)

Citation: Wu X, Wang J, Wu D, Jiang W, Gao Z, Li D, Wu R, Gao D and Zhang Y (2022) Identification of new resistance loci against wheat sharp eyespot through genome-wide association study. Front. Plant Sci. 13:1056935. doi: 10.3389/fpls.2022.1056935

Received: 29 September 2022; Accepted: 28 November 2022;
Published: 12 December 2022.

Edited by:

Xuecai Zhang, International Maize and Wheat Improvement Center, Mexico

Reviewed by:

Yuanfeng Hao, Institute of Crop Sciences (CAAS), China
Caixia Lan, Huazhong Agricultural University, China

Copyright © 2022 Wu, Wang, Wu, Jiang, Gao, Li, Wu, Gao and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yong Zhang, zy@wheat.org.cn

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.