Original Research ARTICLE
Genome-Wide Association Study for Spot Blotch Resistance in Hard Winter Wheat
- Department of Agronomy, Horticulture, and Plant Science, South Dakota State University, Brookings, SD, United States
Spot blotch (SB) caused by Cochliobolus sativus (anamorph: Bipolaris sorokiniana) is an economically important disease of wheat worldwide. Under a severe epidemic condition, the disease can cause yield losses up to 70%. Previous approaches like bi-parental mapping for identifying SB resistant genes/QTLs exploited only a limited portion of the available genetic diversity with a lower capacity to detect polygenic traits, and had a lower marker density. In this study, we performed genome-wide association study (GWAS) for SB resistance in hard winter wheat association mapping panel (HWWAMP) of 294 genotypes. The HWWAMP was evaluated for response to B. sorokiniana (isolate SD40), and a range of reactions was observed with 10 resistant, 38 moderately resistant, 120 moderately resistant- moderately susceptible, 111 moderately susceptible, and 15 susceptible genotypes. GWAS using 15,590 high-quality SNPs and 294 genotypes we identified six QTLs (p = <0.001) on chromosomes 2D, 3A, 4A, 4B, 5A, and 7B that collectively explained 30% of the total variation for SB resistance. Highly associated SNPs were identified for all six QTLs, QSb.sdsu-2D.1 (SNP: Kukri_c31121_1460, R2 = 4%), QSb.sdsu-3A.1 (SNP: Excalibur_c46082_440, R2 = 4%), QSb.sdsu-4A.1 (SNP: IWA8475, R2 = 5.5%), QSb.sdsu-4B.1 (SNP: Excalibur_rep_c79414_306, R2 = 4%), QSb.sdsu-5A.1 (SNP: Kukri_rep_c104877_2166, R2 = 6%), and QSb.sdsu-7B.1 (SNP: TA005844-0160, R2 = 6%). Our study not only validates three (2D, 5A, and 7B) genomic regions identified in previous studies but also provides highly associated SNP markers for marker assisted selection. In addition, we identified three novel QTLs (QSb.sdsu-3A.1, QSb.sdsu-4A.1, and QSb.sdsu-4B.1) for SB resistance in wheat. Gene annotation analysis of the candidate regions identified nine NBS-LRR and 38 other plant defense-related protein families across multiple QTLs, and these could be used for fine mapping and further characterization of SB resistance in wheat. Comparative analysis with barley indicated the SB resistance locus on wheat chromosomes 2D, 3A, 5A, and 7B identified in our study are syntenic to the previously identified SB resistance locus on chromosomes 2H, 3H, 5H, and 7H in barley. The 10 highly resistant genotypes and SNP markers identified in our study could be very useful resources for breeding of SB resistance in wheat.
Wheat (Triticum aestivum L) is an important cereal crop grown worldwide and remains a vital source for human food (FAO, 2017). However, its production is continuously challenged by a number of environmental and biological factors (Sharma et al., 2007b; Gupta et al., 2018). Spot blotch (SB) caused by Bipolaris sorokiniana (Sacc.) Shoem syn. Drechslera sorokiniana (Sacc.) Subrm and Jain (syn. Helminthosporium sativum, teleomorph Cochliobolus sativus) is one of the destructive fungal diseases that affects wheat and several other small grains worldwide (Dubin and Rajaram, 1996; Duveiller and Dubin, 2002; Joshi and Chand, 2002; Sharma et al., 2007b; Singh and Singh, 2007; Gurung et al., 2009; Chowdhury et al., 2013). It has a wide range of hosts within wild and cultivated Poaceae species (Kumar et al., 2002; Pandey et al., 2005; O’Boyle et al., 2014). In susceptible lines, SB symptoms are characterized by small, dark brown lesions that extend 1–2 mm long without chlorotic margins during initial infection (Chand et al., 2003). Later, the leaves are killed when the light brown to dark brown colored oval to elongated blotches extend and merge very quickly. In addition to leaves, the fungus causes common root rot (Wildermuth et al., 1997), seedling blight and seed rot or black point on the embryo (Kumar et al., 2002; Hudec and Muchova, 2008). Average yield loss of 15–20% due to SB has been reported from several countries, but under suitable climatic conditions the losses in yield can reach up to 70% in susceptible genotypes, in addition to the reduction in seed quality (Mehta et al., 1992; Lemerle et al., 1996; Fernandez et al., 1998, 2014; Wang et al., 2002; Fernandez and Jefferson, 2004; Sharma et al., 2004, 2007b; Siddique et al., 2006; Sharma and Duveiller, 2007; Acharya et al., 2011).
Breeding for resistance is one of the most economical and sustainable component of integrated disease management (Crute and Pink, 1996; Bartoš et al., 2002; Duveiller et al., 2007; Gupta et al., 2010; Khan et al., 2010; Chowdhury et al., 2013; Ban et al., 2016; Vasistha et al., 2016; Kumar et al., 2017). However, our ability to develop SB resistant genotypes depends on an understanding of the mechanism of resistance present in the host and identification of the resistant genes responsible for SB resistance. Given the challenges in large-scale germplasm screening (Leng et al., 2016; Osman et al., 2016; Kumar et al., 2017), development of molecular markers linked to disease resistance genes can facilitate marker-assisted selection (MAS) and increase the efficiency of breeding for disease resistance in wheat (Collard et al., 2005; Gupta et al., 2010; Miedaner and Korzun, 2012; Müller et al., 2018).
With the availability of large numbers of molecular markers (Miedaner and Korzun, 2012; Korte and Farlow, 2013) more efficient mapping techniques like genome-wide association studies (GWAS) have become popular for analyzing an unlimited number of traits in genetically identical materials across a wide range of environments (Gupta et al., 2005; Ersoz et al., 2007; Miedaner and Korzun, 2012; Korte and Farlow, 2013; Ogura and Busch, 2015; Kushwaha et al., 2017). GWAS has been used to characterize disease resistance in many crop species: blast resistance genes in rice (Raboin et al., 2016), maize (Xiao et al., 2017), SB resistance in wild barley (Roy et al., 2010), resistance to multiple leaf spot diseases of spring wheat (Gurung et al., 2014), resistance to bacterial leaf streak and SB in spring wheat (Adhikari et al., 2012), Fusarium head blight resistance in wheat (Arruda et al., 2016), tan spot resistance in European winter wheat (Kollers et al., 2014), mapping for resistance to leaf and stripe rust in winter-habit hexaploid wheat landraces (Sun et al., 2015).
Complex quantitative inheritance (Dubin and Rajaram, 1996; Joshi et al., 2004b; Kumar et al., 2015b) of SB resistance in wheat has slowed the progress in breeding for SB resistance. Many studies, using methods of both bi-parental mapping and association mapping (AM) have reported several SB resistance QTLs on chromosome 1A, 1B, 2A, 2B, 2D, 3B, 5A, 5B, 6B, 6D, 7A, 7B, and 7D (Neupane et al., 2007; Sharma et al., 2007a; Gonzalez-Hernandez et al., 2009; Kumar et al., 2009, 2010, 2015a, 2017; Adhikari et al., 2012; Lillemo et al., 2013; Zhuang et al., 2013; Gurung et al., 2014; Zhu et al., 2014; Lu et al., 2016; Zhang et al., 2015; Gupta et al., 2018). However, only three major QTLs designated as Sb1 on 7D (Lillemo et al., 2013), Sb2 on 5B (Kumar et al., 2015a), and Sb3 on 3B (Lu et al., 2016) are well described. Most of these studies have been focused on hard spring wheat, and relatively few studies characterized SB resistance in hard winter wheat germplasm.
Our ability to deploy and develop SB resistant winter wheat cultivars depends on identification of resistant QTLs responsible for the traits. Therefore, the objectives of this study were (i) to identify winter wheat genotypes carrying resistance genes against B. sorokiniana and (ii) to locate putative QTLs and identify SNP markers useful for MAS. This work will contribute toward the development of genome-wide breeding strategies and MAS for SB resistance in wheat.
Materials and Methods
A set of 294 genotypes (Supplementary Table 1) from the hard winter wheat association mapping panel (HWWAMP) of 299 genotypes developed under the USDA TCAP project (Guttieri et al., 2015, 2017) were selected for this study due to availability of sufficient seed. The genotypes were comprised of landraces, advanced breeding lines, and varieties released since the 1940s representing United States Great Plains areas. The HWWAMP panel was designed to include genotypes that represent the existing germplasm of the wheat growing regions of the United States such as Colorado, Kansas, Michigan, Montana, Nebraska, North Dakota, Oklahoma, South Dakota, and Texas (Grogan et al., 2016; Guttieri et al., 2017).
Evaluating HWWAMP for Spot Blotch
The screening of 294 genotypes and the susceptible check, Glenlea and resistant check, Salamouni against B. sorokiniana was conducted at the South Dakota State University Young Brothers Seed Technology greenhouse complex in Brookings, SD, United States. The experiment was conducted in a randomized complete block design (RCBD) with three replications and repeated three times. Three seeds of each genotype were planted in a single root trainer container (Ray Leach “Cone-tainer”TM Single Cell System) and arranged within each tray (Stuewe & Sons, Tangent, OR, United States). Plants were grown in greenhouse at 24°C/16°C (day/night) with 14-h photoperiod and 32% relative humidity. Soluble fertilizer was added every 2 weeks after planting and watering was scheduled every 2 days.
One of the B. sorokiniana isolates, SD40, was used as inoculum in all experiments. SD40 is one of the most virulent isolate predominantly found in United States Great Plains areas and routinely used to screen breeding materials (provided by Dr. Shaukat Ali, SDSU). The isolate was originally derived from a single spore and method of its isolation and cultivation is well described by Kumar et al. (2007). This fungus is easily identifiable based on the color morphology and number of nuclei in lab conditions (Shoemaker, 1959). The stored conidia of the fungus were streaked on V8-PDA (Lamari and Bernier, 1989) plates using a sterile rod. The cultures were grown at room temperature under continuous darkness and harvested 5 days later when the pink colony in the plate began to darken. The spores in the plates were scraped with a flamed microscope slide, diluted with distilled water and transferred to beaker covered by cheesecloth to filter out mycelium. The spore suspension was adjusted to 3500/ml prior to use.
A fully expanded third-leaf of all the seedlings were sprayed with spore suspension using a hand held sprayer (Preval, Chicago Aerosol, Coal City, IL, United States). The spore suspension for inoculations was prepared by adding a 100 μl/L Tween-20 (polyoxyethylene-20-sorbitan monolaurate (Sigma-Aldrich, St. Louis, MO, United States) as a dispersing agent. Approximately 0.5 ml of inoculum suspension was applied to each plant using the sprayer. After inoculation, seedlings were incubated in darkness for 12 h at 20°C in a mist chamber near 100% RH by setting the humidifiers to release the mist for 2 min every 30 min to maintain a humid environment for disease development. Subsequently, plants were transferred to the greenhouse bench where the temperature was regulated at 24°C/16°C (day/night) until the plants were rated for symptom development 10 days post-inoculation.
The infection responses (IRs) of each genotype for SB were assessed based on the five-class (1–5) rating scale used by Lamari and Bernier (1989). IRs were based on the type (presence of necrosis and/or chlorosis) and relative size of lesions observed on the third leaves of the seedlings (Supplementary Figure 1). The third fully expanded inoculated leaf during inoculation was rated as follows: (1) small dark brown to black spot without any surrounding chlorosis or necrosis (resistant), (2) small dark brown to black spot with very little necrosis or chlorosis (moderately resistant), (3) small dark brown to black spot completely surrounded by distinct chlorosis or tan necrosis ring, lesions generally not coalescing (moderately resistant to moderately susceptible), (4) small dark brown to black spot completely surrounded by distinct chlorotic or tan necrotic zone and some of the lesions coalescing (moderately susceptible), and (5) the dark brown or black centers may not be distinguishable, most lesions consist of coalescing chlorotic or tan necrotic zones (susceptible). The infection type corresponding to each score in our experiment is shown in Supplementary Figure 1. In our study, we rated the plants based on necrosis symptoms as SB does not produce chlorotic symptoms in wheat.
Distribution of the phenotypic data for SB was visualized using the histogram. The Shapiro–Wilk test was conducted in R to check the normality of the untransformed and transformed data based on square root method. The homogeneity of the variance across experiments was checked using Bartlett’s test. Data from three repeated experiments, each with three replications, were combined and the overall mean of all experiments was used for analysis of GWAS if the experiment was homogenous. The phenotype data was analyzed using linear mixed model (LMM) approaches with a randomized group-based jackknife technique using R version 3.3.3 (R Core Team, 2016). Broad sense heritability (H2) was estimated by dividing genetic variance over the combined sum of error variance and genetic variance.
The genotype data were obtained from the wheat T3 Toolbox, a public repository.1 The genotyping was conducted using Illumina iSelect 90K under the USDA-TCAP (Guttieri et al., 2017). To avoid spurious marker-trait associations, SNP markers with minor allele frequency (MAF) <0.05 and missing data >10% were excluded from further analyses. The genetic and physical positions of SNP markers from the wheat 90K array were obtained from the consensus map with 46,977 SNPs developed using a combination of eight mapping populations (Wang S. et al., 2014) and the International Wheat Genome Sequencing Consortium website.2 After filtering the high-quality polymorphic SNPs, markers were imputed using TASSEL Version 5.0 software (Bradbury et al., 2007).
Population Structure and Kinship
In order to avoid the distortion of population structure and linkage disequilibrium (LD), SNP markers were first thinned into 0.0005 cM apart to retain only markers with high pairwise correlation using TASSEL v5 software (Bradbury et al., 2007). After keeping informative SNPs in the analysis and eliminating redundant information, we analyzed the genetic stratification, i.e., population structure (Q) within the HWWAMP with STRUCTURE v2.3.4 (Pritchard et al., 2000) using a model-based clustering method. STRUCTURE runs were performed for each specified K values (number of subpopulations, from 2 to 6) using the default setting of the admixture model for the ancestry of individuals and correlated allele frequencies. Burn-in period and a number of Markov Chain Monte Carlo (MCMC) iterations under Linux environment were set to 20,000 and 50,000, respectively (Evanno et al., 2005). The best fit number of clusters was calculated according to Evanno et al. (2005) using STRUCTURE HARVESTER (Duncan et al., 2017). The likely number of population structure was chosen from principal coordinates (PCO) plot, i.e., K vs. ΔK where the rate of change in the log probability between successive K values was the highest.
Linkage Disequilibrium Estimation
Linkage disequilibrium is defined as the non-random association of alleles at different loci in a given population and is represented by the square of the correlation coefficient (r2) between markers. Markers which were in perfect LD (r2 = 1) with another markers were removed before the LD analysis. The r2 between intra- and inter-chromosomal SNP markers were estimated using TASSEL v5 (Bradbury et al., 2007). LD estimates expressed as r2 and based on a sliding window of 100 markers throughout the genome, were calculated and plotted against genetic distance. From the unlinked loci, two markers were considered significant when LD P < 0.001. We plotted the intra-chromosomal r2 values against the genetic distance using excel to see how rapidly the LD decay occurs. The distance at which the smooth curve intercepts the critical r2 was drawn using logarithmic trended smooth lines as described by Hao et al., 2012. A critical value of r2 (basal LD) was estimated using 95% percentile of non-syntenic (inter-chromosomal) r2 distribution below which relationship between two pairs of loci are assumed not to be caused by physical linkage (Laidò et al., 2014). The distance at which the LD decays to 0.7 cM was considered as the critical distance up to which a QTL region can extend (Zhao et al., 2005).
Genome-Wide Association Analysis
Genomic regions associated with SB resistance were identified using TASSEL v.5.0 (Bradbury et al., 2007) and the new enhanced version of the genome association and prediction integrated tool (GAPIT) (Tang et al., 2016) in R version 3.3.3 (R Core Team, 2016). The enhanced version of GAPIT implements computationally powerful statistical approaches such as general linear model (GLM), mixed linear model (MLM) (Zhang Z. et al., 2010), compressed mixed linear model (CMLM) (Li et al., 2014), enhanced compressed mixed linear model (ECMLM) (Li et al., 2014), factored spectrally transformed LMMs (FaST-LMM Select) (Listgarten et al., 2013) and SUPER (Settlement of MLM Under Progressively Exclusive Relationship) (Wang Q. et al., 2014). In GLMs, marker data, disease data and the PCA matrix were integrated as covariates to correct for the effects of population substructure. Unlike GLM, MLM accounts for both population structure and individual kinship as a covariate to reduce type-I error.
We selected the MLM method for our data by comparing the statistical power they have and type I error they produce (results not presented). The MLM for GWAS were represented by y = Xβ + Qv + Ku +e, where y is the vector of the phenotypic values, β is fixed effects due to marker, v is fixed effects due to population structures, e is the vector of residuals, and u is a vector of random effects due to the portion of breeding values not accounted by the marker. X, Q, and K represent matrices from the marker, population structure estimated from the structure or principal component analysis and kinship, respectively. The variance of u is derived as, Var (u) = 2 KVg, where K represents the relative kinship matrix inferred from genotypes of the HWWAMP based on the proportion of shared alleles and Vg is the additive portion of the genetic variance.
In Silico Annotation of SNPs and Syntenic Regions
The sequences of the significant markers associated with SB were extracted from the Infinium iSELECT 90K (Cavanagh et al., 2013) and were BLASTN searched against the Chinese Spring wheat RefSeq v1.0 (IWGSC 20183). The search was limited to the top hit with an E-value cut off of 1E-50 with an identity higher than 75%. Being allohexaploid species (2n = 6x = 42) wheat chromosomes are found in homeologous status (A, B, D) which shared similarities among the homeologous chromosomes. Therefore, we removed SNPs that were mapped to multiple chromosomes.
We further identified the target region for each of the QTLs on pseudomolecule that co-localized with the significant markers contained in each LD block. Next, the sequence segments were BLASTN searched against the wheat coding DNA sequence (CDS) and followed by BLASTX search against the wheat protein database (Duncan et al., 2017). Out of several lists of an annotated protein family, those related to previously described disease resistance protein families were further identified by a BLASTP search against Pfam database (Pfam 31.04). We then compared the candidate regions with barley for comparative analysis of SB resistance genes across related species. We used flanking markers from QTLs identified in our study to BLASTN search against the known barley SB resistance QTLs (Roy et al., 2010; Zhou and Steffenson, 2013) and placed them on the most recent barley genome assembly (Mascher et al., 2017) and then produced synteny representation using circos (Krzywinski et al., 2009).
Variations in Seedling Infection Response
Seedlings showed a range of infection types (Supplementary Figure 1) within 294 hard winter wheat genotypes when inoculated with B. sorokiniana (Figure 1 and Supplementary Table 1). As expected, the resistant check, Salamouni and susceptible check, Glenlea exhibited a mean disease score of 2 and 5, respectively (Figure 1). Out of 294 genotypes, a total of 48 were resistant whereas other 240 genotypes appeared to be in the moderately resistant – moderately susceptible, moderately susceptible and susceptible categories in all three experiments. Out of the 48 resistant genotypes, 10 showed highly resistant (score 1) and 38 showed modrately resistant reaction (score 2) to the SB. The 10 highly resistant genotypes could be potential sources for SB resistance breeding (Table 1). Of the other 248 accessions, 120 showed either moderately resistant – moderately susceptible (score 3) response, whereas, 111 genotypes showed a moderately susceptible (score 4) and 15 showed a highly susceptible (score 5) response to SB of wheat (Figure 1).
FIGURE 1. Frequency distribution of the mean spot blotch infection response of 294 HWWAMP genotypes. The x-axis exhibits 1–5 scores of mean infection response of each genotype. The y-axis represents the number of genotypes (also numbered on the bar) that exhibited the infection response. Salamouni and Glenlea were the resistant and susceptible checks of the experiment, respectively.
TABLE 1. List of 10 genotypes from the HWWAMP that showed highly resistant reactions (R) against B. sorokiniana.
Cultivars and breeding lines that exhibited a highly resistant response to the SB pathogens were: Duster, Colt, Custer, Intrada, MT0495, NE99495, OK04525, OK05122, OK05723W, and Venango (Table 1). The variation among data sets of three different experiments and among the genotypes was analyzed using LMM approaches with a randomized group-based jackknife resampling technique (Supplementary Table 2). At α = 0.05, there were highly significant differences among the predicted genotypic effects (P = 1.33e-14) when compared to the population mean (μ) (Supplementary Table 2). On the other hand, the variations among the disease score of repeated experiments were not significant. There was a high repeatability (96%) among the experiments, and the phenotypic variation explained by the genotypes (broad sense heritability H2) was estimated to be 80% (Supplementary Table 2).
Genotyping, Marker Distribution and Allele Frequency
We obtained genotypic data for 294 HWWAMP constituting 21,555 SNPs from the T3 wheat database.5 We removed 5,487 markers by filtering markers with a (MAF less than 5%. Further, 458 markers with unknown chromosome positions were also removed. Finally, 15,590 high-quality SNP markers across 294 accessions of HWWAMP were used for structure and GWAS analysis (Supplementary Table 3). The map position of 15,590 SNP loci were obtained from Wang S. et al. (2014). The 15,590 SNPs covered map length of 3,661 cM on all 21 wheat chromosomes giving an average interval distance of 0.54 cM (Supplementary Table 3). The number of SNP markers assigned to the A, B, and D chromosomes were 6,211, 7,630, and 1,749, respectively. The individual length chromosomes on genetic map ranged from 119 to 244 cM. The average number of markers per cM on genome A (4.97) and B (6.6) were relatively higher when compared to D genome (1.5) (Supplementary Table 3). The chromosome 4D (52) and 7D (133) harbored the lowest number of informative markers. Each locus was characterized by the presence of major and minor allele with a frequency between 0.50 to 0.95, and 0.05 to 0.50, respectively.
Linkage disequilibrium is the non-random co-segregation of alleles at two or more loci on the same chromosome or between loci on different chromosomes. Out of 15,990 markers used for AM, only 1,842 markers were used for LD analysis taking out non-informative markers. Among a total of 91,307 locus pairs detected, 13,076 locus pairs (14.3%) were in LD at P < 0.001 of which 7,744 locus pairs (8.5%) were found at r2 > 0.1 and P < 0.001 (Supplementary Table 4). However, the distance at which LD starts decaying depends on meiotic events and/or genetic drifts occured in the population. We estimated the LD decay distances in the whole genome and within each genome of winter wheat association panel using logarithmic trended smooth lines developed from scatter plots of syntenic r2 vs. the genetic distance (cM) between pairs of two-locus (Supplementary Figure 2). We estimated the distance point where LD value (r2) decreases below 0.1 or half strength of D’ (D’ = 0.5) based on the curve of the nonlinear logarithmic trend line. The LD (r2 > 0.1) decay distance of about 4.5 cM was estimated for the whole genome. Whereas, LD decay distances for individual A, B, and D genomes were approximately 3.4, 3.6, and 14.2 cM, respectively (Supplementary Figure 2).
Population Stratification and Kinship
In order to avoid distortion due to population structure, SNP markers were first thinned to retain only 1,842 markers that were at least 0.5 cM apart in TASSEL 5. A STRUCTURE analysis indicated the presence of four subpopulations (K1–K4) in the HWWAMP where the four clusters contained 45, 37, 99, and 113 genotypes, respectively (Supplementary Figures 3, 4). Average distances (expected heterozygosity) between individuals within each cluster (K1–K4) were 0.06, 0.17, 0.27, and 0.23, respectively. Similarly, the net nucleotide distance among structures, i.e., the average probability that a pair of alleles were different among K1 vs. K2, K3, K4 were 0.28, 0.19, 0.16; K2 vs. K3, K4 were 0.24, 0.21; and K3 vs. K4 were 0.10. We further divided the individuals based on their inferred ancestry and made a principal component analysis. Individuals within and across different structures shared genetic similarities, and were considered admixed when the cumulative shared ancestry across the structures was above 40% or they retained greater than 60% ancestry within their respective structure. From a PCA analysis, the proportion and cumulative variances of the first four PCA were 0.43, 0.29, 0.28, 0.00 and 0.43, 0.72, 1.00, 1.00, respectively (results not shown). Similarly, a matrix of kinships among individual genotypes was calculated using all 15,990 SNPs. The heat map calculated using the classical equation from VanRaden (2008) showed a high kinship relationship among individuals (Supplementary Figure 5).
Model Comparison for Marker-Trait Association
Six statistical models such as GLM, MLM, CMLM, ECMLM, FaST-LMM, and SUPER were compared to select the one which reduces the type-I error and increases the power of SNP discovery. Our analysis indicated that MLM, CMLM, and ECMLM similarly reduced the type-I error and increased power when compared to others. The analysis showed no differences in the number of significant SNPs discovered using all the models except GLM and FAST. However, due to its power and robustness, we selected MLM (Zhang Z. et al., 2010) with optimum compression as the model of choice for subsequent analysis. The quantile-quantile (Q–Q) plots for the test statistics (Supplementary Figure 6) indicated the absence of inflation of statistics or overall systematic bias caused by the population stratification when MLM was used than GLM. The number of significant markers associated with the SB response with the GLM were much greater than with the MLM (Supplementary Table 5). Indeed, at P < 0.001 there were 13 and 26 significant makers trait associations that were identified by MLM and GLM, respectively. The 13 markers identified by MLM were also significant by GLM (Supplementary Table 5).
Markers Associated With the Spot Blotch QTLs
GWAS analysis identified several genomic regions linked to SB resistance including some in genomic regions identified in earlier studies (Joshi et al., 2004a; Lillemo et al., 2013; Kumar et al., 2015a; Lu et al., 2016). We identified six genomic regions showing highly significant (P < 0.001) markers associated with SB resistance on chromosomes 2D, 3A, 4A, 4B, 5A, and 7B (Figure 2). A total 30% of the variation was explained by the most significant SNP from each genomic region. QTLs for SB resistance have been reported in the regions similar to three of the six QTLs identified in our study on chromosomes 2D (QSb.sdsu-2D.1), 5A (QSb.sdsu-5A.1), and 7B (QSb.sdsu-7B.1). We identified three novel regions that contribute to SB resistance. The QTLs, QSb.sdsu-3A.1, QSb.sdsu-4A.1, and Q.Sb.sdsu-4B.1 (P < 0.001) explained 4, 6, and 4% of the variation, respectively (Table 2). We identified several co-localized SNPs markers that were associated to each QTL region (Supplementary Tables 6, 7). We further identified the QTL spanning region based on the significant SNP markers in LD. The SNP marker Kukri_rep_c104877_2166 which was at 59.1 cM (480,285,174 bp) on chromosome 5A had the highest association (P = 4.02E-05 and R2 = 5.9%). Likewise, group of 7B markers such as Excalibur_c5700_670, Kukri_c21628_1215, Tdurum_contig9966_646, Kukri_c22495_552, Excalibur_c5700_527, Excalibur_c58742_144, TA005844-0160, Excalibur_c5700_705, BS00075332_51, and Tdurum_contig90495_232 were found at about 86.1 cM were high associated with SB resistance and explained similar phenotypic variation (R2 = 6.3%) (Supplementary Table 6).
FIGURE 2. Genome wide association scan for SB resistance in hard winter wheat association mapping panel (HWWAMP). Manhattan plot developed using mixed linear model (MLM) (A) and general linear model (GLM) (B) in TASSEL v.5. The -log10 (P) values from a genome-wide scan are plotted against positions on each of the 21 wheat chromosomes. Horizontal lines indicate genome-wide significance thresholds.
TABLE 2. Summary of SNP markers linked to significant SB resistance QTLs detected from genome-wide association analysis of 294 winter wheat genotypes.
QTL Effects on SB Resistance
QTLs with negative additive effect had favorable alleles and showed a reduction in SB whereas QTLs with positive additive effect had alleles showing unfavorable effect on SB reaction (Table 2 and Supplementary Tables 8, 9). The combination of individual QTL effects produces total resistance effects that the genotype exhibit. The cultivar, Colt, for instance, was one of the most resistant genotypes and harbored desirable alleles on chromosome 2D, 3A, 4A, and 7B. The combination of the four alleles reduced the disease by a score 1.73 when compared to the mean score of the HRWWAMP and hence exhibited the highest resistance (Table 2 and Supplementary Tables 8, 9). Similarly, OK05723 and Duster had desirable allele on chromosome 4A and 7B, and 3A, 4A, and 7B, respectively, that gave higher protection against SB. QTLs QSb.sdsu-2D.1, QSb.sdsu-3A.1, QSb.sdsu-4A.1, and QSb.sdsu-7B.1 could give the highest protection against the disease if a genotype harbored a combination of the following desirable alleles: CC, CC, GG, and CC, respectively (Supplementary Table 9). For QTLs QSb.sdsu-4B.1 and QSb.sdsu-5A.1 the best approach would be to have genotypes with alleles AA and TT, respectively (Supplementary Table 9).
In Silico Functional Annotation of the Candidate Region
Plant reactions to diseases are very complex, and involve the activation of sets of genes, encoding for different proteins. To facilitate the identification of additional markers to identify candidate proteins that may be involved in pathogenesis-related (PR) response we annotated the coding sequences in the candidate regions (Figure 3 and Supplementary Table 10). The chromosome regions spanning through the set of SNP markers that showed P < 0.05 up-and downstream of the highest significant marker was identified for further analysis (Figure 3). After BLASTN (1e-50) search of markers sequence to reference wheat genome (RefSeq v1.0), a 2.07, 56.7, 6.58, 1.31, 10.03, and 3.00 Mb segment of DNA were identified on chromosomes 2D, 3A, 4A, 4B, 5A, and 7B, respectively and were used for functional annotation. Out of the total of 129 proteins annotated in the candidate regions, we found 9, 22, 5, 2, 9, and 3 disease-related proteins in candidate regions on chromosomes 2D, 3A, 4A, 4B, 5A, and 7B, respectively (Supplementary Table 10). We further studied QTLs QSb.sdsu-7B.1 and QSb.sdsu-2D.1 in detail as they explained the maximum variation of SB resistance in our study and were flanked in relatively smaller regions. Figure 3 shows the SNP associations with SB response and is plotted as Manhattan plot (Figures 3A-I,B-I) along the consensus genetic map (Wang S. et al., 2014). Within the selected region of the up-and-downstream of the most significant markers, we compared the local r2 LD pattern and two haplotype blocks were identified on each of chromosome 7B (Figure 3A-II) and 2D (Figure 3B-II). The haplotype block DNA segment on 7B was 2.6 Mb (1.1 cM) (Figure 3A-III) whereas 2D haplotype block was 2.1 Mb (0.7 cM) (Figure 3B-III) and the functions of plant disease defense-related genes in the haplotype blocks (Figures 3A-IV,B-IV) were identified.
FIGURE 3. Gene annotation of QTLs identified on chromosome 7B (A) and 2D (B) for spot blotch resistance in hard winter wheat. The far left image is a Manhattan plot indicating the level of marker association with the trait (I). Next is a visualization of linkage disequilibrium (black is a D’ value of 100%, another color is a D’ value of less than 100%), significant haplotype blocks are outlined in black (II). Names of the markers in the region of interest along with their cM position (III), and physical position (IV). Most significant marker is highlighted in red. The far right image is a physical map of candidate genes on 7B and 2D chromosome segments spanning from 608.7 to 611.7 Mb and 606.9 to 608.8 Mb, respectively (IV and V). Genes in bold code Pathogenesis-related (PR) proteins. The cM position is based on Wang S. et al. (2014) and the physical position is based on IWGSC 2018.
Our gene annotation analysis identified several NB-ARC and NBS-LRR proteins containing N-terminal nucleotide-binding site (NBS) and C-terminal leucine-rich repeat (LRRs) in QTLs spanning regions of QSb.sdsu-2D.1, QSb.sdsu-3A.1, QSb.sdsu-4A.1, and QSb.sdsu-7B.1. Similarly, Chitinase class I protein was observed in the QTL region QSb.sdsu-5A.1. The protein description, functions, ID and other related information are presented in Supplementary Table 10. Proteins related to cell death and response to oxidative stress such as peroxidase (POX) superfamily and zinc-binding dehydrogenase protein family were observed on QSb.sdsu-7B.1 and QSb.sdsu-5A.1, respectively. In addition to NBS-LRR proteins QSb.sdsu-2D.1 also harbors zinc finger domain and RING/U-box superfamily proteins that may play important role in resistance to fungal pathogens. Further, protein families reported as transporters (EamA-like transporter family), endopeptidase inhibitors (serine protease inhibitor – SERPIN), kinase activity (protein kinase superfamily and tyrosine kinase), and detoxification (ATP-binding cassette (ABC) transporter C and cytochrome P450) were also found in QTL candidate regions. Besides the main families, proteins with unknown function or indirectly involved in disease resistance like senescence and dehydration-associated protein, ABC transporters were identified. Similarly, protein families such as glycosyl hydrolase, mitochondrial glycoprotein, auxin signaling F-box 2 (AFB2), pentatricopeptide repeat (PPR) superfamily, proteins were found in the candidate regions of few QTLs.
Comparison of SB Resistance QTLs Between Wheat and Barley
Shared synteny can be one of the most reliable criteria for establishing the orthology of genomic regions in different species. We performed a comparative analysis of wheat and barley for the candidate regions of six QTLs identified in our study (Figure 4). Synteny analysis indicated that four QTLs (QSb.sdsu-2D.1, QSb.sdsu-3A.1, QSb.sdsu-5A.1, and QSb.sdsu-7B.1) on chromosomes 2D, 3A, 5A, and 7B carrying SB resistance QTLs in wheat corresponds to 2H, 3H, 5H, and 7H chromosomes in barley (Figure 4). SB resistance genes/QTLs have been reported in these syntenic regions in barley (Roy et al., 2010; Zhou and Steffenson, 2013) further validating the QTLs identified in our study.
FIGURE 4. Comparative analysis showing synteny in spot blotch resistance QTLs mapped on wheat (TGACv1 2014) and barley (Mascher et al., 2017) chromosomes. The map was made using CIRCOS (Circular Genome Data Visualization). Each color indicates a different chromosome. Arcs in bold shows the corresponding marker related to the respective QTL.
Phenotypic Variability for Resistance Against Spot Blotch of Wheat
Developing SB resistant wheat cultivars is likely the most economical and durable strategy for minimizing loss due to SB when compared to the overall usage of fungicides (Duveiller and Gilchrist, 1994; Duveiller and Sharma, 2009; McMullen and Adhikari, 2009). Several sources of SB resistance have been reported in the Indian national collection (Kumar et al., 2017), CIMMYT germplasm and derivatives of CIMMYT primary synthetic bread wheat (Mikhailova et al., 2004; Zhu et al., 2014), multi-resistant cultivars from Nepal (Mahto et al., 2011) and some Brazilian varieties (Mehta, 1998). However, most of these studies have been conducted in spring wheat. In winter wheat, modern European winter wheat cultivars and breeding lines (Liatukas and Ruzgas, 2012) have been evaluated, but the United States germplasm has not been extensively exploited. In the present study, the 294 HWW genotypes showed a variable reaction against SB indicating the existence of genetic variability among them similar to the observations reported in spring wheat (Khan and Chowdhury, 2011; Kumar et al., 2017; Ojha et al., 2017). We observed a repeatability of 96% and a heritability of 80% demonstrating the reliability of the SB resistance evaluation under greenhouse conditions. However, studies have shown both positive and negative correlation for response to SB wheat between in vivo and in vitro conditions (Rosyara et al., 2009; Liatukas and Ruzgas, 2012). We identified 10 genotypes (Table 1) highly resistant to SB that can be utilized for developing SB resistance in winter wheat cultivars. Colt [Agate sib: (NE69441)//(Tx65A1503-1) 391-56-D8/Kaw] was one of the top cultivars that exhibited a reproducible and highly resistant reaction to the SB pathogen (Table 1). It was among the first semi-dwarf wheat ever released in Nebraska in 1983. Another cultivar Duster, one of the most popular variety in the Oklahoma since 2006, was found to be highly resistant and could be a promising source of SB resistance. It has a pedigree background of W0405D/NE78488//W7469C/TX81V6187. Similarly, Custer which had been a top performer in Oklahoma in 1987, and some of its parents (F-29-76/TAM-105//Chisholm) were among the HWWAMP accessions were also evaluated in this experiment. Both F-29-76 and TAM-105 showed susceptible reaction whereas Chisholm exhibited resistant reaction suggesting that Chisholm served as sources of resistance in Custer. Of the 168 winter wheat cultivars released in the United States since the green revolution of 1970s and evaluated in our study, only 16% of the genotypes showed a resistant and moderately resistant response (Supplementary Table 1). A similar percentage (16%) of the evaluated breeding lines were also resistant and can be employed for SB resistance breeding.
SNP Marker Distribution and LD Decay
The distribution and density of informative markers reflects the overall genetic richness and diversity of the wheat genome. The D genome had the least average marker density ranging from 0.3–4.2 markers/cM as compared to A (3.9–6.4) and B (3.8–7.2) genomes. Further, the total markers used in for AM from D genome (11.2%) were very small when compared to A (39.8%) and B (48.9%) genomes suggesting a lower genetic diversity and lower level of effective recombination in the D genome as observed in prevoius studies Barnes et al. (2012) and Nielsen et al. (2014). LD is an integral part of AM and the estimation of LD decay distances can help in determining the power of AM. We estimated a critical value of r2 (basal LD = 0.1) from 95th percentile distribution of the inter-chromosomal LD (Francki et al., 2006; Laidò et al., 2014) below which we assumed the absence of physical relationship within pairs of loci and hence they are physically present on different chromosomes. Our results also demonstrated that the D genome had high LD decay distance (14 cM) when compared to A and B genomes (Supplementary Figure 2). Whole genome LD decay distance we observed in HWWAMP was similar to whole genome LD decay observed of in Chinese winter wheat using SSR marker (Chao et al., 2010; Barnes et al., 2012). However, a lower rate of average LD decay (higher distance) was observed in European hexaploid wheat (23 cM) (Nielsen et al., 2014) and US Elite hard red winter wheat (10 cM) at r2 > 0.1(Zhang D. et al., 2010). When comparing the wheat landraces to modern cultivars, the LD decay distance (5 cM) was lower in landraces as compared to modern winter wheat cultivars (5–10 cM) (Hao et al., 2011) suggesting a possible reduction in diversity in the modern wheat cultivars in China. The LD decay analysis in HWWAMP showed variation among genomes and within the genome itself, indicating variability in recombination hot spots, differences in selection pressure imposed on alleles of wheat genomes, and an evidence of the recombination events in past breeding history.
QTLs for SB Resistance
We identified six genomic locations on chromosomes 2D, 3A, 4A, 4B, 5A, and 7B associated with SB resistance in HWWAMP and the linked SNPs (TA005844-0160, Kukri_rep_c104877_2166, IWA8475, Kukri_c31121_1460, Excalibur_rep_c79414_306, and Excalibur_c46082_440) explained about 30% of the total variation. We compared the MTA with MLM and GLM (Figure 2), and observed that GLM analysis exhibited low power and a higher risk of false-positive detection (result not presented). In the GLM, which only accounts for population structure as a covariate, the additive variance and error variance could not be separated as GLM uses maximum compression (compression = n) with all taxa as a single group. Unlike GLM, however, MLM (Zhang Z. et al., 2010) takes account of population structure and individual kinship in association analysis to reduce type I error instigated due to relatedness and population structure. We observed similar results with MLM, cMLM (Tang et al., 2016) and ECMLM (Li et al., 2014). All subpopulations showed a good degree of admixture which could be a result of frequently exchanged germplasm among the winter wheat breeding programs (Supplementary Figure 3).
We identified three novel QTLs for SB resistance (QSb.sdsu-3A.1, QSb.sdsu-4A.1, and Q.Sb.sdsu-4B.1) in addition to three regions on chromosomes 2D, 5A, and 7B reported earlier to be significantly associated with SB resistance. Previously, Gurung et al. (2014) reported significant QTLs associated with SB on chromosomes 1B, 5A, 5B, 6B, and 7B, whereas Adhikari et al. (2012) detected QTLs on chromosomes 1A, 3B, 7B, and 7D, and Lillemo et al. (2013) detected QTLs on chromosomes 5B, 7A, and 7D. Similarly, SB resistance QTLs were reported on wheat chromosomes 1B, 3B, and 5A (Zhu et al., 2014); 2AL, 2BS, 5BL, and 6DL (Kumar et al., 2009); 2BS, 2DS, 3BS, 7BS and 7DS (Kumar et al., 2010); 5B, 6A, and 6D (Sharma et al., 2007b) and 7B and 7D on (Ban et al., 2016). With dense marker coverage, we not only validate the QTLs QSb.sdsu-2D.1, QSb.sdsu-5A.1, and QSb.sdsu-7B.1 on chromosomes 2D, 5A, and 7B but also provided highly significant associated SNP markers that could be used for MAS for SB resistance. However, with a limited number of common markers between previous SSR or DArT based studies, comparisons could be distorted and allelism tests will be required to determine which of the QTLs identified in our study have been previously identified.
In our analysis, none of the three major QTLs, Sb1 (Joshi et al., 2004a; Lillemo et al., 2013), Sb2 (Kumar et al., 2015a), and Sb3 (Lu et al., 2016) contributing to SB resistance that have been characterized using simple sequence repeat (SSR) markers were significant at P < 0.001 but all these QTLs showed a peak and were significant at P < 0.005 suggesting the presence of these QTLs in the HWWAMP. Each of the three QTLs was explaining ∼3% variation in HWWAMP but was lower than was reported in previous studies (Joshi et al., 2004a,b, Lillemo et al., 2013; Kumar et al., 2015a,b; Lu et al., 2016). The high density of marker coverage in our study and the availability of wheat genome sequences permitted comparison with barley (Figure 4). Identification of several R gene clusters and mapping of SB QTLs in syntenic regions in barley suggest four wheat QTLs could be orthologous to barley and this information could support the beginning for the search of candidate genes in wheat and understanding the mechanism of SB resistance.
Functional Annotations of Candidate Regions
Plant defense systems can be categorized into ever existing constitutive defense systems that are triggered by pathogen-associated molecular patterns (PAMPs) or the temporarily induced system that targets to defend an attacked area of the plant. The genes encoding the specificity determinants of effector-triggered immunity are known as resistance (R) genes (DeYoung and Innes, 2006). The production of PR proteins in response to pathogens are the primary mechanisms in induced plant’s self-defense system. Numerous PR proteins have been characterized in recent years, and largely classified into at least 17 protein families and several pathogenesis-related proteins that do not constitute a superfamily of proteins (Dangl and Jones, 2001; Sels et al., 2008). The genomic regions spanning the SB QTLs identified in our study harbored many genes. However, not all genes are equally important in the regulation of quantitative traits like SB resistance. The candidate genes that were commonly found across multiple QTLs reported are more likely the ones that determine the trait (Swamy et al., 2011). In our study, we found NB-ARC and NBS-LRR in many of the annotated QTL regions (Figure 3 and Supplementary Table 10). The NBS-LRR are the most common R-genes by which highly conserved NBS domains can bind and hydrolyze ATP or GTP, whereas the LRR motif is typically involved in protein-protein interactions and is responsible for recognition specificity (Wan et al., 2012; Zhang et al., 2016).
Peroxidase superfamily protein that plays a role in self-defense (Hiraga et al., 2001) by catalyzing oxidoreduction between H2O2 and various reductants was one of the classical enzymes that was observed in the QSb.sdsu-7B.1 region on chromosome 7B. Chitinase class I was observed in the QSb.sdsu-5A.1 region. Plant chitinases are the proteins that hydrolyzes the N-acetylglucosamine polymer chitin and also takes part in pathogenesis related activities (Punja and Zhang, 1993). Another disease related protein family identified in the candidate regions was protein kinase that not only play major role in phosphorylation in plants (Mulekar and Huq, 2013) but are also involved in signal transduction and activation of plant defense mechanisms (Goff and Ramonell, 2007). Many plant defense related genes such as OsMPK6 in rice (Ueno et al., 2015) and Pto in tomato (Martin et al., 1993) have been reported to translate protein kinase proteins. In our study protein kinase protiens were observed in QTL region QSb.sdsu-3A.1, QSb.sdsu-4A.1, and QSb.sdsu-5A.1. Further, the zinc finger proteins identified in some of the candidate regions in our study have been reported to regulate resistance mechanism through its active involvement in sequence-specific binding to DNA/RNA and contribution in protein-protein recognitions (Gupta et al., 2012). All the above described disease-related genes identified in the candidate regions could help in the development of a new markers and further characterization of the SB resistance QTLs in wheat.
Implications of GWAS for SB Resistance Breeding in Wheat
The ultimate goal of characterizing SB resistance genes is to find closely linked markers for assisting in the selection and further understanding the underlying network of genes and their interactions to achieve resistance responses. This comprehensive understanding will help in developing durable disease resistant cultivars. We identified groups of SNP markers associated with six QTLs that had different levels of effects (Table 2 and Supplementary Table 6). A few genotypes used in our study encompass multiple resistant alleles and showed highest level of resistance responses (Table 1 and Supplementary Table 1). Past inheritance studies on resistance to SB suggest polygenic types of resistance that appears to be based on many minor genes with small individual effects (Dubin and Rajaram, 1996; Joshi et al., 2004b; Gurung et al., 2014; Kumar et al., 2015b). Backcrosses of two parents harboring Qsb.bhu-2A on chromosome 2A and Qsb.bhu-5B on chromosome 5B in suitable parents achieved higher resistance in susceptible cultivar HUW 234 in India (Vasistha et al., 2016). This suggested that assembling multiple desirable alleles for SB QTLs from multiple parents into a single genotype through the use of marker-assisted gene pyramiding could help in improving disease resistance in wheat varieties. We identified that winter wheat cv. Custer (OK) carries desirable alleles for SB QTLs on chromosomes 2D, 3A, 4A, and 7B and also demonstrated a very high level of SB resistance. The significant SNP markers linked to these SB QTLs (Table 2 and Supplementary Table 7) can be used to develop a Kompetitive Allele Specific PCR (KASP) assay and used for MAS, as it is well evidenced in several crops.
We identified 10 winter wheat genotypes highly resistant to SB and six genomic regions associated with SB resistance along with tightly linked SNPs. Genotypes with multiple SB resistance QTLs could be used for future breeding and the linked SNP markers could facilitate quick mobilization of the SB resistance QTLs.
GA, SA, and SS designed the experiments. GA and SA performed the SB evaluations. GA, JS, and SS performed the data analysis. GA and SS wrote the manuscript. JS, SA, JGH, and BT contributed to the interpretation of results and also critically reviewed the manuscript.
This project was collectively funded by the USDA hatch project SD00H538-15, USDA NACA 58-3020-5-025, and South Dakota Wheat Commission 3X7262. The genotyping of the hard winter wheat association mapping panel was performed under Agriculture and Food Research Initiative Competitive Grants 2011-68002-30029 (Triticeae-CAP) and 2017-67007-25939 (Wheat-CAP) from the USDA National Institute of Food and Agriculture. The funders had no role in the study design, data collection and analysis, decision to publish or preparation of the manuscript.
Conflict of Interest Statement
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.
We would like to thank the South Dakota Agriculture Experimental Station (Brookings, SD, United States) for providing the resources to conduct the experiments. We would also like to thank IWGSC for providing pre-publication access to RefSeq v1.0 for functional analysis of SB QTL regions.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.00926/full#supplementary-material
- ^ http://triticeaetoolbox.org
- ^ http://www.wheatgenome.org
- ^ https://wheat-urgi.versailles.inra.fr/Seq-Repository/Assemblies
- ^ https://xfam.wordpress.com/2017
- ^ https://triticeaetoolbox.org/wheat
Adhikari, T. B., Gurung, S., Hansen, J. M., Jackson, E. W., and Bonman, J. M. (2012). Association mapping of quantitative trait loci in spring wheat landraces conferring resistance to bacterial leaf streak and spot blotch. Plant Genome 5, 1–16. doi: 10.3835/plantgenome2011.12.0032
Arruda, M. P., Brown, P., Brown-Guedira, G., Krill, A. M., Thurber, C., Merrill, K. R., et al. (2016). Genome-wide association mapping of Fusarium head blight resistance in wheat using genotyping-by-sequencing. Plant Genome 9, 1–14. doi: 10.3835/plantgenome2015.04.0028
Ban, T., Singh, V., Singh, G., Chaudhury, A., Ojha, A., Tyagi, B. S., et al. (2016). Phenotyping at hot spots and tagging of QTLs conferring spot blotch resistance in bread wheat. Mol. Biol. Rep. 43, 1293–1303. doi: 10.1007/s11033-016-4066-z
Barnes, J., Chen, X., Min, D., Yasir, T. A., and Hu, Y. G. (2012). Genetic diversity, population structure and linkage disequilibrium in elite Chinese winter wheat investigated with SSR markers. Int. J. Mol. Sci. 7:e44510.
Bartoš, P., Šip, V., Chrpová, J., Vacke, J., Stuchliková, E., Blažková, V., et al. (2002). Achievements and prospects of wheat breeding for disease resistance. Czech J. Genet. Plant Breed. 38, 16–28. doi: 10.17221/6107-CJGPB
Bradbury, P. J., Zhang, Z., Kroon, D. E., Casstevens, T. M., Ramdoss, Y., and Buckler, E. S. (2007). TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics 23, 2633–2635. doi: 10.1093/bioinformatics/btm308
Cavanagh, C. R., Chao, S., Wang, S., Huang, B. E., Stephen, S., Kiani, S., et al. (2013). Genome-wide comparative diversity uncovers multiple targets of selection for improvement in hexaploid wheat landraces and cultivars. Proc. Natl. Acad. Sci. U.S.A. 110, 8057–8062. doi: 10.1073/pnas.1217133110
Chand, R., Pandey, S. P., Singh, H. V., Kumar, S., and Joshi, A. K. (2003). Variability and its probable cause in natural populations of spot blotch pathogen Bipolaris sorokiniana of wheat (T. aestivum L.) in India. J. Plant Dis. Prot. 110, 27–35.
Chao, S., Dubcovsky, J., Dvorak, J., Luo, M. C., Baenziger, S. P., Matnyazov, R., et al. (2010). Population- and genome-specific patterns of linkage disequilibrium and SNP variation in spring and winter wheat (Triticum aestivum L.). BMC Genomics 11:727. doi: 10.1186/1471-2164-11-727
Collard, B. C. Y., Jahufer, M. Z. Z., Brouwer, J. B., and Pang, E. C. K. (2005). An introduction to markers, quantitative trait loci (QTL) mapping and marker-assisted selection for crop improvement: the basic concepts. Euphytica 142, 169–196. doi: 10.1007/s10681-005-1681-5
Duveiller, E., and Gilchrist, L. (1994). “Production constraints due to Bipolaris sorokiniana in wheat: current situation and future prospects,” in Proceedings of the Wheat in Warm Area, Rice–Wheat Farming Systems, Dinajpur.
Duveiller, E., Singh, R. P., and Nicol, J. M. (2007). The challenges of maintaining wheat productivity: pests, diseases, and potential epidemics. Euphytica 157, 417–430. doi: 10.1007/s10681-007-9380-z
Duveiller, E. M., and Sharma, R. C. (2009). Genetic improvement and crop management strategies to minimize yield losses in warm non-traditional wheat growing areas due to spot blotch pathogen Cochliobolus sativus. J. Phytopathol. 157, 521–534. doi: 10.1111/j.1439-0434.2008.01534.x
Ersoz, E. S., Yu, J., and Buckler, E. S. (2007). “Applications of linkage disequilibrium and association mapping in crop plants,” in Genomics-assisted Crop Improvements, eds R. J. Varshney, and R. Tuberosa, (Berlin: Springer), 97–119.
Evanno, G., Regnaut, S., and Goudet, J. (2005). Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol. Ecol. 14, 2611–2620. doi: 10.1111/j.1365-294X.2005.02553.x
Fernandez, M. R., Fox, S. L., Hucl, P., and Singh, A. K. (2014). Leaf spotting reaction of spring common, durum and spelt wheat, and Kamut under organic management in western Canada. Can. J. Plant Sci. 94, 929–935. doi: 10.4141/cjps2013-310
Fernandez, M. R., Zentner, R. P., McConkey, B. G., and Campbell, C. A. (1998). Effects of crop rotations and fertilizer management on leaf spotting diseases of spring wheat in southwestern Saskatchewan. Can. J. Plant Sci. 78, 489–496. doi: 10.4141/P97-082
Francki, M., Breseghello, F., and Sorrells, M. E. (2006). Association mapping of kernel size and milling quality in wheat (Triticum aestivum L.) cultivars. Genetics 172, 1165–1177. doi: 10.1534/genetics.105.044586
Gonzalez-Hernandez, J. L., Singh, P. K., Mergoum, M., Adhikari, T. B., Kianian, S. F., Simsek, S., et al. (2009). A quantitative trait locus on chromosome 5B controls resistance of Triticum turgidum (L.) var. diccocoides to Stagonospora nodorum blotch. Euphytica 166, 199–206. doi: 10.1007/s10681-008-9825-z
Grogan, S. M., Anderson, J., Baenziger, P. S., Frels, K., Guttieri, M. J., Haley, S. D., et al. (2016). Phenotypic plasticity of winter wheat heading date and grain yield across the US Great Plains. Crop Sci. 56, 2223–2236. doi: 10.2135/cropsci2015.06.0357
Gupta, P. K., Chand, R., Vasistha, N. K., Pandey, S. P., Kumar, U., Mishra, V. K., et al. (2018). Spot blotch disease of wheat: the current status of research on genetics and breeding. Plant Pathol. 67, 508–531. doi: 10.1111/ppa.12781
Gupta, P. K., Rustgi, S., and Kulwal, P. L. (2005). Linkage disequilibrium and association studies in higher plants: present status and future prospects. Plant Mol. Biol. 57, 461–485. doi: 10.1007/s11103-005-0257-z
Gupta, S. K., Rai, A. K., Kanwar, S. S., and Sharma, T. R. (2012). Comparative analysis of zinc finger proteins involved in plant disease resistance. PLoS One 7:e42578. doi: 10.1371/journal.pone.0042578
Gurung, S., Mamidi, S., Bonman, J. M., Xiong, M., Brown-Guedira, G., and Adhikari, T. B. (2014). Genome-wide association study reveals novel quantitative trait loci associated with resistance to multiple leaf spot diseases of spring wheat. PLoS One 9:e108179. doi: 10.1371/journal.pone.0108179
Guttieri, M. J., Baenziger, P. S., Frels, K., Carver, B., Arnall, B., and Waters, B. M. (2015). Variation for grain mineral concentration in a diversity panel of current and historical Great Plains hard winter wheat germplasm. Crop Sci. 55, 1035–1052. doi: 10.2135/cropsci2014.07.0506
Guttieri, M. J., Frels, K., Regassa, T., Waters, B. M., and Baenziger, P. S. (2017). Variation for nitrogen use efficiency traits in current and historical great plains hard winter wheat. Euphytica 213:87. doi: 10.1007/s10681-017-1869-5
Hao, C., Wang, L., Ge, H., Dong, Y., and Zhang, X. (2011). Genetic diversity and linkage disequilibrium in Chinese bread wheat (Triticum aestivum L.) revealed by SSR markers. PLoS One 6:e17279. doi: 10.1371/journal.pone.0017279
Hao, D., Cheng, H., Yin, Z., Cui, S., Zhang, D., Wang, H., et al. (2012). Identification of single nucleotide polymorphisms and haplotypes associated with yield and yield components in soybean (Glycine max) landraces across multiple environments. Theor. Appl. Genet. 124, 447–458. doi: 10.1007/s00122-011-1719-0
Joshi, A. K., and Chand, R. (2002). Variation and inheritance of leaf angle, and its association with spot blotch (Bipolaris sorokiniana) severity in wheat (Triticum aestivum). Euphytica 124, 283–291. doi: 10.1023/a:1015773404694
Joshi, A. K., Kumar, S., Chand, R., and Ortiz-Ferrara, G. (2004b). Inheritance of resistance to spot blotch caused by Bipolaris sorokiniana in spring wheat. Plant Breed. 123, 213–219. doi: 10.1111/j.1439-0523.2004.00954.x
Khan, H., and Chowdhury, S. (2011). Identification of resistance source in wheat germplasm against spot blotch disease caused by Bipolaris sorokiniana. Arch. Phytopathol. Plant Prot. 44, 840–844. doi: 10.1080/03235400903345273
Kollers, S., Rodemann, B., Ling, J., Korzun, V., Ebmeyer, E., Argillier, O., et al. (2014). Genome-wide association mapping of tan spot resistance (Pyrenophora tritici-repentis) in European winter wheat. Mol. Breed. 34, 363–371. doi: 10.1007/s11032-014-0039-x
Krzywinski, M., Schein, J., Birol, I., Connors, J., Gascoyne, R., Horsman, D., et al. (2009). Circos: an information aesthetic for comparative genomics. Genome Res. 19, 1639–1645. doi: 10.1101/gr.092759.109
Kumar, D., Chand, R., Prasad, L. C., and Joshi, A. K. (2007). A new technique for monoconidial culture of the most aggressive isolate in a given population of Bipolaris sorokiniana, cause of foliar spot blotch in wheat and barley. World J. Microbiol. Biotechnol. 23, 1647–1651. doi: 10.1007/s11274-007-9410-y
Kumar, J., Schäfer, P., Hückelhoven, R., Langen, G., Baltruschat, H., Stein, E., et al. (2002). Bipolaris sorokiniana, a cereal pathogen of global concern: cytological and molecular approaches towards better control. Mol. Plant Pathol. 3, 185–195. doi: 10.1046/j.1364-3703.2002.00120.x
Kumar, S., Archak, S., Tyagi, R. K., Kumar, J., Vikas, V. K., Jacob, S. R., et al. (2017). Evaluation of 19,460 wheat accessions conserved in the Indian national Genebank to identify new sources of resistance to rust and spot blotch diseases. PLoS One 12:e0175610. doi: 10.1371/journal.pone.0175610
Kumar, S., Roder, M. S., Tripathi, S. B., Chand, R., Joshi, A. K., and Kumar, U. (2015a). Mendelization and fine mapping of a bread wheat spot blotch disease resistance QTL. Mol. Breed. 35:218. doi: 10.1007/s11032-015-0411-5
Kumar, S., Tripathi, S. B., and Kumar, U. (2015b). Dissection of wheat spot blotch disease resistance QTLs in to single Mendelian genes. Indian J. Genet. Plant Breed. 75, 434–439. doi: 10.5958/0975-6906.2015.00070.x
Kumar, U., Joshi, A. K., Kumar, S., Chand, R., and Roder, M. (2009). Mapping of resistance to spot blotch disease caused by Bipolaris sorokiniana in spring wheat. Theor. Appl. Genet. 118, 783–792. doi: 10.1007/s00122-008-0938-5
Kumar, U., Joshi, A. K., Kumar, S., Chand, R., Roder, M. S., and Röder, M. S. (2010). Quantitative trait loci for resistance to spot blotch caused by Bipolaris sorokiniana in wheat (T. aestivum L.) lines “Ning 8201” and “Chirya 3”. Mol. Breed. 26, 477–491. doi: 10.1007/s11032-009-9388-2
Laidò, G., Marone, D., Russo, M. A., Colecchia, S. A., Mastrangelo, A. M., De Vita, P., et al. (2014). Linkage disequilibrium and genome-wide association mapping in tetraploid wheat (Triticum turgidum L.). PLoS One 9:e95211. doi: 10.1371/journal.pone.0095211
Lamari, L., and Bernier, C. C. (1989). Evaluation of wheat lines and cultivars to tan spot [Pyrenophora tritici-repentis] based on lesion type. Can. J. Plant Pathol. 11, 49–56. doi: 10.1080/07060668909501146
Lemerle, D., Yuan, T. H., Murray, G. M., and Morris, S. (1996). Survey of weeds and diseases in cereal crops in the southern wheat belt of New South Wales. Aust. J. Exp. Agric. 36, 545–554. doi: 10.1071/ea9960545
Leng, Y. Q., Wang, R., Ali, S., Zhao, M. X., and Zhong, S. B. (2016). Sources and genetics of spot blotch resistance to a new pathotype of Cochliobolus sativus in the USDA National small grains collection. Plant Dis. 100, 1988–1993. doi: 10.1094/pdis-02-16-0152-re
Li, M., Liu, X., Bradbury, P., Yu, J., Zhang, Y.-M., Todhunter, R. J., et al. (2014). Enrichment of statistical power for genome-wide association studies. BMC Biol. 12:73. doi: 10.1186/s12915-014-0073-5
Lillemo, M., Joshi, A. K., Prasad, R., Chand, R., and Singh, R. P. (2013). QTL for spot blotch resistance in bread wheat line Saar co-locate to the biotrophic disease resistance loci Lr34 and Lr46. Theor. Appl. Genet. 126, 711–719. doi: 10.1007/s00122-012-2012-6
Lu, P., Liang, Y., Li, D., Wang, Z., Li, W., Wang, G., et al. (2016). Fine genetic mapping of spot blotch resistance gene Sb3 in wheat (Triticum aestivum). Theor. Appl. Genet. 129, 577–589. doi: 10.1007/s00122-015-2649-z
Mahto, B. N., Gurung, S., and Adhikari, T. B. (2011). Assessing genetic resistance to spot blotch, Stagonospora nodorum blotch and tan spot in wheat from Nepal. Eur. J. Plant Pathol. 131, 249–260. doi: 10.1007/s10658-011-9803-5
Martin, G. B., Brommonschenkel, S. H., Chunwongse, J., Frary, A., Ganal, M. W., Spivey, R., et al. (1993). Map-based cloning of a protein kinase gene conferring disease resistance in tomato. Science 262, 1432–1436. doi: 10.1126/science.7902614
Mascher, M., Gundlach, H., Himmelbach, A., Beier, S., Twardziok, S. O., Wicker, T., et al. (2017). A chromosome conformation capture ordered sequence of the barley genome. Nature 544, 427–433. doi: 10.1038/nature22043
McMullen, M., and Adhikari, T. (2009). Fungal Leaf Spot Diseases of Wheat: Tan Spot, Stagonospora nodorum blotch and Septoria tritici blotch. North Dakota State University Extension Service, 1249. Available at: https://www.ag.ndsu.edu/ndipm/publications/wheat/documents/pp1249.pdf
Mehta, Y. R. (1998). “Constraints on the integrated management of spot blotch of wheat,” in Helminthosporium Blight Wheat Spot Blotch Tan Spot, eds E. Duveiller, H. J. Dubin, J. Reeves, and A. McNab (Mexico, DF: CIMMYT), 18–27.
Mehta, Y. R., Riede, C. R., Campos, L. A. C., and Kohli, M. M. (1992). Integrated management of major wheat diseases in brazil - an example for the southern cone region of latin-America. Crop Prot. 11, 517–524. doi: 10.1016/0261-2194(92)90168-5
Mikhailova, L. A., Lianfa, S., Gogoleva, S. G., and Gultyaeva, E. I. (2004). Sources and donors of resistance to wheat spot blotch caused by Bipolaris sorokiniana Shoem.(Cochliobolus sativus Drechs. ex Dastur). Arch. Phytopathol. Plant Prot. 37, 161–167. doi: 10.1080/03235400310001631909
Müller, T., Schierscher-Viret, B., Fossati, D., Brabant, C., Schori, A., Keller, B., et al. (2018). Unlocking the diversity of genebanks: whole-genome marker analysis of Swiss bread wheat and spelt. Theor. Appl. Genet. 131, 407–416. doi: 10.1007/s00122-017-3010-5
Neupane, R. B., Sharma, R. C., Duveiller, E., Ortiz-Ferrara, G., Ojha, B. R., Rosyara, U. R., et al. (2007). Major gene controls of field resistance to spot blotch in wheat genotypes “Milan/Shanghai #7” and “Chirya.3”. Plant Dis. 91, 692–697. doi: 10.1094/pdis-91-6-0692
Nielsen, N. H., Backes, G., Stougaard, J., Andersen, S. U., and Jahoor, A. (2014). Genetic diversity and population structure analysis of European hexaploid bread wheat (Triticum aestivum L.) varieties. PLoS One 9:e94000. doi: 10.1371/journal.pone.0094000
O’Boyle, P. D., Brooks, W. S., Barnett, M. D., Berger, G. L., Steffenson, B. J., Stromberg, E. L., et al. (2014). Mapping net blotch resistance in “Nomini” and CIho 2291 Barley. Crop Sci. 54, 2596–2602. doi: 10.2135/cropsci2013.08.0514
Ogura, T., and Busch, W. (2015). From phenotypes to causal sequences: using genome wide association studies to dissect the sequence basis for variation of plant development. Curr. Opin. Plant Biol. 23, 98–108. doi: 10.1016/j.pbi.2014.11.008
Ojha, A., Singh, G., Tyagi, B., Singh, V. A., Rajita, and Kumar, P. (2017). Screening of resistance source against spot blotch disease caused By Bipolaris Sorokiniana in Triticum aestivum L. Int. J. Adv. Res. 5, 23–28. doi: 10.21474/IJAR01/2716
Osman, M., He, X. Y., Benedettelli, S., Ali, S., and Singh, P. K. (2016). Identification of new sources of resistance to fungal leaf and head blight diseases of wheat. Eur. J. Plant Pathol. 145, 305–320. doi: 10.1007/s10658-015-0843-0
Pandey, S. P., Kumar, S., Kumar, U., Chand, R., and Joshi, A. K. (2005). Sources of inoculum and reappearance of spot blotch of wheat in rice–wheat cropping. Eur. J. Plant Pathol. 111, 47–55. doi: 10.1007/s10658-004-2404-9
Raboin, L.-M., Ballini, E., Tharreau, D., Ramanantsoanirina, A., Frouin, J., Courtois, B., et al. (2016). Association mapping of resistance to rice blast in upland field conditions. Rice 9:59. doi: 10.1186/s12284-016-0131-4
Rosyara, U. R., Khadka, K., Subedi, S., Sharma, R. C., and Duveiller, E. (2009). Field resistance to spot blotch is not associated with undesirable physio-morphological traits in three spring wheat populations. J. Plant Pathol. 91, 113–122.
Roy, J. K., Smith, K. P., Muehlbauer, G. J., Chao, S., Close, T. J., and Steffenson, B. J. (2010). Association mapping of spot blotch resistance in wild barley. Mol. Breed. 26, 243–256. doi: 10.1007/s11032-010-9402-8
Sels, J., Mathys, J., De Coninck, B. M. A., Cammue, B. P. A., and De Bolle, M. F. C. (2008). Plant pathogenesis-related (PR) proteins: a focus on PR peptides. Plant Physiol. Biochem. 46, 941–950. doi: 10.1016/j.plaphy.2008.06.011
Sharma, R. C., Duveiller, E., and Jacquemin, J. M. (2007a). Microsatellite markers associated with spot blotch resistance in spring wheat. J. Phytopathol. 155, 316–319. doi: 10.1111/j.1439-0434.2007.01238.x
Sharma, R. C., Duveiller, E., and Ortiz-Ferrara, G. (2007b). Progress and challenge towards reducing wheat spot blotch threat in the Eastern Gangetic Plains of South Asia: Is climate change already taking its toll? Field Crop. Res. 103, 109–118. doi: 10.1016/j.fcr.2007.05.004
Sharma, R. C., Sah, S. N., and Duveiller, E. (2004). Combining ability analysis of resistance to Helminthosporium leaf blight in spring wheat. Euphytica 136, 341–348. doi: 10.1023/B:EUPH.0000032752.75504.bd
Sun, D., Kertho, A., Mamidi, S., Bonman, J. M., McClean, P. E., and Acevedo, M. (2015). Genome-wide association mapping for resistance to leaf and stripe rust in winter-habit hexaploid wheat landraces. PLoS One 10:e0129580. doi: 10.1371/journal.pone.0129580
Swamy, B. P. M., Vikram, P., Dixit, S., Ahmed, H. U., and Kumar, A. (2011). Meta-analysis of grain yield QTL identified during agricultural drought in grasses showed consensus. BMC Genomics 12:319. doi: 10.1186/1471-2164-12-319
Tang, Y., Liu, X., Wang, J., Li, M., Wang, Q., Tian, F., et al. (2016). GAPIT Version 2: an enhanced integrated tool for genomic association and prediction. Plant Genome 9, 1–9. doi: 10.3835/plantgenome2015.11.0120
Ueno, Y., Yoshida, R., Kishi-Kaboshi, M., Matsushita, A., Jiang, C.-J., Goto, S., et al. (2015). Abiotic stresses antagonize the rice defence pathway through the tyrosine-dephosphorylation of OsMPK6. PLoS Pathog. 11:e1005231. doi: 10.1371/journal.ppat.1005231
Vasistha, N. K., Balasubramaniam, A., Mishra, V. K., Chand, R., Srinivasa, J., Yadav, P. S., et al. (2016). Enhancing spot blotch resistance in wheat by marker-aided backcross breeding. Euphytica 207, 119–133. doi: 10.1007/s10681-015-1548-3
Wan, H., Yuan, W., Ye, Q., Wang, R., Ruan, M., Li, Z., et al. (2012). Analysis of TIR-and non-TIR-NBS-LRR disease resistance gene analogous in pepper: characterization, genetic variation, functional divergence and expression patterns. BMC Genomics 13:502. doi: 10.1186/1471-2164-13-502
Wang, H., Fernandez, M. R., Clarke, F. R., DePauw, R. M., and Clarke, J. M. (2002). Effect of leaf spotting diseases on grain yield and seed traits of wheat in southern Saskatchewan. Can. J. Plant Sci. 82, 507–512. doi: 10.4141/P01-128
Wang, S., Wong, D., Forrest, K., Allen, A., Chao, S., Huang, B. E., et al. (2014). Characterization of polyploid wheat genomic diversity using a high-density 90 000 single nucleotide polymorphism array. Plant Biotechnol. J. 12, 787–796. doi: 10.1111/pbi.12183
Wildermuth, G. B., Thomas, G. A., Radford, B. J., McNamara, R. B., and Kelly, A. (1997). Crown rot and common root rot in wheat grown under different tillage and stubble treatments in southern Queensland, Australia. Soil Tillage Res. 44, 211–224. doi: 10.1016/s0167-1987(97)00054-8
Zhang, D., Bai, G., Zhu, C., Yu, J., and Carver, B. F. (2010). Genetic diversity, population structure, and linkage disequilibrium in US elite winter wheat. Plant Genome 3, 117–127. doi: 10.3835/plantgenome2010.03.0004
Zhang, K., Nair, S. K., Babu, R., Magorokosho, C., Mahuku, G., Semagn, K., et al. (2015). Fine mapping of Msv1, a major QTL for resistance to Maize Streak Virus leads to development of production markers for breeding pipelines. J. Genet. Genomics 128, 1839–1854.
Zhang, Y., Xia, R., Kuang, H., and Meyers, B. C. (2016). The diversification of plant NBS-LRR defense genes directs the evolution of microRNAs that target them. Mol. Biol. Evol. 33, 2692–2705. doi: 10.1093/molbev/msw154
Zhang, Z., Ersoz, E., Lai, C.-Q., Todhunter, R. J., Tiwari, H. K., Gore, M. A., et al. (2010). Mixed linear model approach adapted for genome-wide association studies. Nat. Genet. 42, 355–360. doi: 10.1038/ng.546
Zhao, H., Nettleton, D., Soller, M., and Dekkers, J. C. M. (2005). Evaluation of linkage disequilibrium measures between multi-allelic markers as predictors of linkage disequilibrium between markers and QTL. Genet. Res. 86, 77–87. doi: 10.1017/S001667230500769X
Zhou, H., and Steffenson, B. (2013). Genome-wide association mapping reveals genetic architecture of durable spot blotch resistance in US barley breeding germplasm. Mol. Breed. 32, 139–154. doi: 10.1007/s11032-013-9858-4
Zhu, Z. W., Bonnett, D., Ellis, M., Singh, P., Heslot, N., Dreisigacker, S., et al. (2014). Mapping resistance to spot blotch in a CIMMYT synthetic-derived bread wheat. Mol. Breed. 34, 1215–1228. doi: 10.1007/s11032-014-0111-6
Zhuang, Y., Gala, A., and Yen, Y. (2013). Identification of functional genic components of major Fusarium head blight resistance quantitative trait loci in wheat cultivar Sumai 3. Mol. Plant Microbe Interact. 26, 442–450. doi: 10.1094/mpmi-10-12-0235-r
Keywords: Triticum aestivum, hard winter wheat, spot blotch, GWAS, QTLs, SNPs, marker-assisted selection
Citation: Ayana GT, Ali S, Sidhu JS, Gonzalez Hernandez JL, Turnipseed B and Sehgal SK (2018) Genome-Wide Association Study for Spot Blotch Resistance in Hard Winter Wheat. Front. Plant Sci. 9:926. doi: 10.3389/fpls.2018.00926
Received: 07 May 2018; Accepted: 11 June 2018;
Published: 06 July 2018.
Edited by:Anna Maria Mastrangelo, Centro di Ricerca per l’Orticoltura (CRA), Italy
Reviewed by:Ivan A. Matus, Instituto de Investigaciones Agropecuarias (INIA), Chile
Marion S. Röder, Leibniz-Institut für Pflanzengenetik und Kulturpflanzenforschung (IPK), Germany
Copyright © 2018 Ayana, Ali, Sidhu, Gonzalez Hernandez, Turnipseed and Sehgal. 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: Sunish K. Sehgal, email@example.com