Next Article in Journal
Dissection of the Genetic Association between Anorexia Nervosa and Obsessive–Compulsive Disorder at the Network and Cellular Levels
Next Article in Special Issue
Improvement of a Yairipok Chujak Maize Landrace from North Eastern Himalayan Region for β-Carotene Content through Molecular Marker-Assisted Backcross Breeding
Previous Article in Journal
Characterization of Four Orphan Receptors (GPR3, GPR6, GPR12 and GPR12L) in Chickens and Ducks and Regulation of GPR12 Expression in Ovarian Granulosa Cells by Progesterone
Previous Article in Special Issue
Application of Genomics to Understand Salt Tolerance in Lentil
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Association Studies and Prediction of Tan Spot (Pyrenophora tritici-repentis) Infection in European Winter Wheat via Different Marker Platforms

1
Leibniz Institute of Plant Genetics and Crop Plant Research (IPK), Corrensstraße 3, D-06466 Stadt Seeland OT Gatersleben, Germany
2
European Wheat Breeding Center, BASF Agricultural Solutions GmbH, Am Schwabeplan 8, D-06466 Stadt Seeland OT Gatersleben, Germany
3
Julius-Kühn-Institute (JKI), D-38104 Braunschweig, Germany
4
TraitGenetics GmbH, Am Schwabeplan 1b, D-06466 Stadt Seeland OT Gatersleben, Germany
*
Author to whom correspondence should be addressed.
Genes 2021, 12(4), 490; https://doi.org/10.3390/genes12040490
Submission received: 18 February 2021 / Revised: 17 March 2021 / Accepted: 25 March 2021 / Published: 27 March 2021

Abstract

:
Tan spot, caused by the fungus Pyrenophora tritici-repentis (Ptr), is a severe foliar disease of wheat (Triticum aestivum L.). Improving genetic resistance is a durable strategy to reduce Ptr-related losses. Here, we dissected Ptr-infection’s genetic basis in 372 European wheat varieties via simple sequence repeats (SSRs) plus 35k and 90k single nucleotide polymorphism (SNP) marker platforms. In our phenotypic data analyses, Ptr infection showed a significant genotypic variance and a significant negative correlation with plant height. Genome-wide association studies revealed a highly quantitative nature of Ptr infection and identified two quantitative trait loci (QTL), viz., QTs.ipk-7A and QTs.ipk-7B, which imparted 21.23 and 5.84% of the genotypic variance, respectively. Besides, the Rht-D1 gene showed a strong allelic influence on the infection scores. Due to the complex genetic nature of the Ptr infection, the potential of genome-wide prediction (GP) was assessed via three different genetic models on individual and combined marker platforms. The GP results indicated that the marker density and marker platforms do not considerably impact prediction accuracy (~40–42%) and that higher-order epistatic interactions may not be highly pervasive. Our results provide a further understanding of Ptr-infection’s genetic nature, serve as a resource for marker-assisted breeding, and highlight the potential of genome-wide selection for improved Ptr resistance.

1. Introduction

Tan spot, also known as the yellow leaf spot, is a severe disease of wheat worldwide. Caused by the fungal pathogen Pyrenophora tritici-repentis, (Ptr; Died.) anamorph Drechslera tritici-repentis (Dtr; Died.) Shoem. (syn. Helminthosporium tritici-repentis), the Ptr infection is mainly diagnosed by tan-colored necrotic lesions with yellow margins that are often surrounded by chlorotic haloes on susceptible wheat leaves. Mature lesions have a dark area in the center. With time, the lesions become larger and often fuse, resulting in the decrease of leaves’ photosynthetic surface area. Consequently, dead leaf tissue areas translate to plant stress and eventually yield loss [1]. Ptr-associated yield losses prove more detrimental- especially at adult stages, e.g., between growth stages BBCH-45 and -65, i.e., mid of booting to mid of flowering [2]. The yield losses—that may reach up to 50%—are mainly attributed to the reduction in (1) leaf area index, (2) dry matter accumulation, and (3) the number of reproductive tillers [3,4]. Besides, reduced kernel size, kernel weight, and the number of kernels per ear were reported to be the main drivers of Ptr-associated yield losses [5]. The fungal spores overwinter in the previous wheat crop’s stubble residue and reproduce in the following spring and summer [6,7]. In high disease pressure, the tan spot can also infect ears and eventually kernels, leading to the seeds’ red- or pink-smudge disease [1]. No- or minimum-tillage practices were reported to result in high disease infestation where infected stubble or kernels from previous cropping seasons act as a disease inoculum [6,7,8]. The absence of cover crops and weedicide application coupled with susceptible wheat lines, a favorable environment (i.e., rainy summer), and no-tillage help the fungus flourish. The fungus produces at least three necrotrophic effectors (NEs), viz., Ptr-ToxA, Ptr-ToxB, and Ptr-ToxC (for reviews, see [9,10]). The NEs—previously called host-selective toxins—are recognized by host sensitivity (S) genes and lead to dominant susceptibility [9]. The lack of fungal NEs recognition by the host (wheat) results in an incompatible interaction and leads to resistant wheat lines. Based on the three NEs mentioned earlier, the Ptr isolates have been classified into eight races [11].
Farm or agronomic management practices, e.g., primary and secondary tillage, crop rotation, and cultivar mixtures, are suitable measures to prevent disease-associated damages [7]. However, the accompanying monetary demerits may prevent their continuous use, especially by smallholdings. On the other hand, the timely use of broad-spectrum foliar fungicides—especially in times of high pressure of multiple allied diseases, e.g., Septoria tritici blotch and Stagonospora nodorum blotch—can help prevent the disease spread and benefit economically by higher yields. Nevertheless, extensive fungicide applications may result in a high pathogen evolution rate and are not sustainable. Hence, improving the resistance by exploiting genetics is deemed as a durable strategy for sustainable gains.
Three dominant S-genes, viz., Tsn1, Tsc1, and Tsc2, have been identified on wheat chromosomes 5BL, 1AS, and 2BS, respectively [12,13,14]. Tsn1—the first and the only gene cloned thus far for the tan spot necrosis [15]—interacts with Ptr-ToxA, whereas Tsc1 and Tsc2—for tan spot chlorosis—were reported to interact with Ptr-ToxC and Ptr-ToxB, respectively. Besides, four tan spot resistance (tsr) qualitative genes, viz., tsr2–tsr5 (syn. tsn2–tsn5), were reported [16,17,18,19]. The presence of S-genes or absence of tsr-genes leads to cultivar susceptibility. The tan spot’s genetic architecture has been studied mainly via bi-parental mapping studies primarily to identify large-effect loci (reviewed in [10]), and, as a result, virtually tens of quantitative trait loci (QTL) have been identified, many of which correspond to the already identified S- or tsr-genes [20]. Genome-wide association studies (GWAS) that exploit the allelic diversity in diverse lines have also been performed to elucidate the tan spot’s genetic basis. Gurung et al. [21] were the first to show the potential of GWAS for tan spot to identify QTL in diverse spring wheat landraces. Since then, several studies report the QTL associated with both seedling and adult plant tan spot susceptibility and resistance in different panels comprising diverse landraces, breeding lines, and elite released varieties of both spring and winter wheat habitats [22,23,24,25,26,27,28]. Genome-wide prediction (GP) is a slightly different but related approach that exploits genome-wide markers’ effects rather than only the significant loci to predict the individual’s genetic merit for the trait under selection [29]. Recent GP studies on wheat diseases suggest its promising potential in breeding for improved quantitative resistance [30,31,32].
Here, we dissected the genetic basis of the Ptr infection in a diverse panel of recently registered 372 European wheat varieties previously studied only with the simple sequence repeats markers. We improved the molecular data by fingerprinting the varieties with high-density 35k and 90k single-nucleotide polymorphism marker arrays. We identified large-effect Ptr-associated QTL by combining all marker platforms suggesting the use of improved marker density. In addition, we studied the prospects of genome-wide selection (GS) by checking the efficiency of the individual marker platform to predict Ptr-infection’s genetic value. The GP accuracies showed that GS could be performed to improve quantitative genetic resistance and that marker platform, or marker density, does not substantially impact prediction accuracy.

2. Materials and Methods

2.1. Collection and Analyses of the Phenotypic Data

A panel (GABI) of European wheat lines comprising 372 varieties (358 winter type; 14 spring type) was evaluated for tan spot (Ptr) infection/resistance. The Ptr-infection’s phenotypic data were gathered from three replications in two environments, with each environment considered a location-by-year combination. The inoculation was performed by using a mixture of various German tan spot field isolates. Ten flag and ten first leaves were evaluated from every genotype in each replication for the percentage of Ptr infected area. The average percent Ptr infected area from all leaves was taken to represent each variety’s overall Ptr score in each replication. A detailed protocol for inoculation at various growth stages and disease scoring methodology is provided in Kollers et al. [23]. The field trials were conducted in α-lattice design. More details about the field trials, agronomic practices, climatic conditions, and calculation of the across-replications arithmetic entry means of each genotype in individual environments have been described previously [23]. Since disease data are generally skewed, we performed the square-root transformation on the individual environment’s data to improve the statistical normality. The normality of the phenotypic data was assessed via the Shapiro–Wilk test at P = 0.001 .
To compute the across-environment individual variance components of the genotype, environment, and the residuals, the following linear mixed-effect model was used by assuming all effects except the intercept as random:
y i j = μ + g i + e j + ε i j
where y i j is the phenotypic value (arithmetic mean) of the i t h genotype in the j t h environment, μ is the common intercept term, g i is the effect of the i t h genotype, e j is the effect of the j t h environment, and ε i j is the corresponding residual term as ε ~ N ( 0 , I σ ε 2 ) with I and σ ε 2 being the identity matrix and residual variance. The broad-sense heritability ( H 2 ) was calculated as:
H 2 = σ g 2 σ g 2 + ( σ ε 2 n E )
where σ g 2 and σ ε 2 denote the genotype and residuals’ variance components, respectively, and n E represents the number of environments. The best linear unbiased estimations (BLUEs) across environments were calculated by assuming the intercept and genotype effects fixed in Equation (1). Since plant height (PH; cm) and heading date (HD; the number of days counted after 1st January) are purposed as morphological escape traits for various diseases [32,33], we retrieved data from previously published multiple-environment studies on the same panel [34,35]. We calculated the genetic correlations among all the traits based on their BLUEs computed across environments.

2.2. Collection and Analyses of the Genotypic Data

All 372 wheat varieties were genotyped with marker platforms, viz., microsatellites (simple sequence repeats; SSRs), and single nucleotide polymorphism (SNP) arrays. In total, the SSR genotyping resulted in 732 markers with 782 scorable genetic loci representing 3178 (2581 mapped and 597 unmapped) alleles, as described previously [23]. For SNP genotyping of the panel, two state-of-the-art marker platforms, viz., 35k Affymatrix breeders’ array and 90k Illumina iSELECT array were employed which generated 35,143 and 81,587 markers ( p ) , respectively. Besides, we genotyped the whole panel with functional markers for the candidate genes, such as photoperiodism (Ppd-D1) and reduced height (Rht-B1 and Rht-D1). Detailed information about the primer design for the candidate genes is given in Kollers et al. [23]. The SSR markers’ genetic positions were taken from the International Triticeae Mapping Initiative (ITMI) DH mapping population described in Sorrells et al. [36]. On the other hand, SNP markers from both 35k and 90k arrays were anchored onto the physical map of wheat (RefSeq v1.1), and the physical position of the markers and their corresponding information, e.g., location, gene-ID, and gene-length (start and end positions) were retrieved from Sun et al. [37]. In total, of the 35k and 90k SNP arrays, 26,236 (74.65%) and 60,638 (74.32%) makers were physically mapped onto the chromosomes. The SNP markers from both arrays plus the SSRs and candidate-gene markers’ scores were combined, which resulted in an n × p matrix of 372 × 119,966 and subjected to the quality check. The quality criteria were implemented to remove the markers with a minimum of 0.05 minor allele frequency and >5% missing or heterozygous calls; the remaining missing or heterozygous calls were imputed with the mean value of both alleles.

2.3. Genome-Wide Association Studies

Genome-wide association studies (GWAS) were performed on data taken from the individual environment and markers (both SSRs and SNPs) passing the quality criteria and the functional-gene markers. Let n be the varieties and p the predictor marker genotypes. Following Yu et al. [38], a standard linear mixed-effect model was used to perform GWAS as:
y = 1 μ + E τ + X β + P v + Z u + ε
where y is the column vector of adjusted means of each genotype calculated in the individual environment, μ is the common intercept, τ ,   β , v , u , and ε are the vectors of the individual environment, markers, population structure (principal components), polygenic background, and the error effects, respectively; E , X , P , and Z are the corresponding design matrices. In the model, μ , τ ,   β , and v were assumed to be fixed while u and ε as random with u ~ N ( 0 , G σ u 2 ) and ε ~ N ( 0 , I σ ε 2 ) . The n × n variance-covariance additive relationship matrix ( G ) was calculated from an n × p matrix W = ( w i k ) of marker genotypes (being 0, 1, or 2) as:
G = k = 1 p ( w i k 2 p k ) ( w j k 2 p k ) 2 k = 1 p p k ( 1 p k )
where w i k and w j k are the profiles of the k t h marker for the i t h and j t h variety, respectively; p k is the estimated frequency of one allele in k t h marker, described as a second solution in VanRaden [39]. Since population stratification and familial relatedness can severely impact the power to detect the real marker-trait associations (MTA) in GWAS, different methods were used to correct for population stratification and relatedness viz., (1) multiple linear regression (naïve), (2) correction of population structure by the first three principal components (PC[1–3]), (3) correction of familial relatedness via genomic relationship matrix G , and (4) correction of both population structure and familial relatedness by PC[1–3] and G . It is expected that using both PCs and G in the model better corrects for the false positives. The models described above were compared by plotting expected versus observed log 10 ( P ) values in the form of a quantile-quantile (qq) plot. The most conservative model was determined by checking how well the observed log 10 ( P ) values aligned with the expected.
To declare the MTA, a liberal false discovery rate (FDR) to account for multiple testing was applied at α = 0.20 [40]. As described in Utz et al. [41], the percentage of the adjusted genotypic variance ( p G ) explained by all QTL was determined as:
p G = ( R a d j 2 H 2 ) × 100
where, R a d j 2 was calculated as R a d j 2 = R 2 ( z N z 1 ) ( 1 R 2 ) by fitting the MTA ( z ) in the order of their descending P-values in a multiple linear regression model; R 2 ,   N , and H 2 denote the regression coefficient, number of observations, and the broad-sense heritability calculated in Equation (2), respectively. The p G explained by the individual MTA was accordingly calculated from their sum of squares. The identified QTL were named based on recommended rules for gene or QTL symbolization in wheat (available online: https://wheat.pw.usda.gov/ggpages/wgc/98/Intro.htm, accessed on 18 January 2021).

2.4. Genome-Wide Predictions

Genome-wide prediction (GP) studies were performed by using three different models with different assumptions, viz., genomic best linear unbiased predictions (GBLUP), Bayesian alphabet B (BayesB), and reproducing kernel Hilbert space regression (RKHSR). GBLUP is a standard robust parametric procedure to predict the total genetic value of the trait under consideration by exploiting additive effects of the markers assuming equal variances [39,42]. It is a linear model of the form:
y = 1 μ + g + ε
where y is the column vector of BLUEs calculated across environments in Equation (1), μ is a common intercept, and g = Z u ; the Z u and ε are explained in Equation (3).
Since the distribution of marker variances across loci is not always equal, the BayesB model, which is of the form:
y = 1 μ + X β + ε
utilizes a scaled inverse Chi-squared ( χ 2 ) distribution on the marker variances. This circumvents the problem of equal variance by assuming a prior distribution ( π ; the prior proportion of non-zero effects) that yields a scaled t-distribution for marker effects by using both shrinkage and variable selection methods. Here, y is explained in Equation (6), and X , β and ε are explained in Equation (3). Following Pérez and de los Campos [43], the prior distribution can be modeled as:
p ( β j , σ β 2 , π ) = { k [ π N ( β j k | 0 , σ β 2 ) + ( 1 π ) 1 ( β j k = 0 ) ] χ 2 ( σ β j k 2 | d f β , S β ) } × B ( ( π | p 0 , π 0 ) × G ( S β | r , s )
where N and B denote normal and beta densities; β and σ β 2 represent the vector of regression coefficients and respective variance. To set the hyper-parameters, we implemented the built-in procedures of the BGLR statistical package [43].
The RKHSR is a semiparametric method that accounts for the additive as well as epistatic interactions among loci [44]. It is of the same form as GBLUP (Equation (6)) with the assumption that g = K α , and thus—by using Gaussian kernel—can be represented as:
y = 1 μ + K α + ε
where y , μ , and ε are the same as described in Equation (6), and α is the vector of random effects with α ~ N ( 0 ,   K σ α 2 ) . Here, K is n × n symmetric positive-definite matrix and is defined as K i j = e ( h × d i j 2 p ) where, K i j represents the measured relationship between the i t h and j t h variety based on their marker profiles, d i j 2 is the Euclidean distance between the i t h and j t h variety and h is the bandwidth parameter. To determine the optimum h , three different values as h = 0.5 × ( 1 / 5 , 1 , 5 ) were tested in a five-fold cross-validation scenario, and the value representing the highest accuracy was chosen.
We evaluated the genome-wide prediction accuracy ( r G P ) of all models using a five-fold cross-validation scenario, as described in Muqaddasi et al. 2020. Briefly, the varieties were randomly divided into five subsets; four were used as the training set to estimate the remaining test set’s genetic values. The accuracy of prediction was defined as the Pearson’s product-moment correlation between the observed ( y ) and predicted ( y ^ ) genetic values standardized by the square-root of the broad-sense heritability as r G P = c o r ( y ,   y ^ ) H . Since the cross-validation runs were repeated for 100 cycles, mean and standard deviation values were calculated to show the individual prediction model’s performance. Unless stated otherwise, all calculations were performed in R software [45] mainly by using lme4 [46] and rrBLUP [47] packages.

3. Results

3.1. Phenotypic Data Analyses Reveal Significant Genetic Variation and a Strong Negative Correlation of Tan Spot Infection with Plant Height

The tan spot (Ptr) infection assessment on 372 wheat varieties registered primarily for European markets was performed in replicated field trials. The phenotypic data from the individual environment was square-root transformed (Table S1). We observed a moderate ( r = 0.20 ) but significant Pearson’s product-moment correlation between both environments’ adjusted means (Figure S1a). The ANOVA showed that both genotype and environment variance was significantly ( P < 0.001 ) larger than zero (Table 1). The best linear unbiased estimations (BLUEs) calculated across environments approximated a statistically normal distribution (Shapiro–Wilk P = 0.003 ) and ranged from 1.58 to 3.97 with a mean of 2.51 and median of 2.48; the 1st and the 3rd quantiles amounted to 2.23 and 2.77, respectively (Figure 1a and Figure S1b). The broad-sense heritability amounted to 0.33, suggesting a sizeable environmental variance; this is expected due to uneven disease pressures in different environments.
We retrieved data for plant height and heading date from previously published studies to observe their influence on the tan spot infection. We observed a highly significant negative Pearson’s product-moment correlation of tan spot infection with PH while a moderate negative correlation with HD (Figure S1c). This indicates that taller and later heading plants—on average—escape Ptr infection and that shorter plants are more susceptible to the disease infestation.

3.2. GWAS Reveals Medium- to Large-Effect Loci Controlling the Tan Spot

We performed GWAS based on environment-specific phenotypic scores and the genotypic matrix comprising the full set of quality markers ( p = 28 , 114 ) that were combined from SSRs, two SNP arrays, and candidate-gene markers. It was shown earlier that, on this panel, increasing the marker density results in improved detection of the marker-trait associations (MTA) [32]. In this study, the GWAS model correcting both population structure and genomic relationships (Figure S2) could sufficiently control spurious MTA detection (Figure 1b,c). Our GWAS resulted in the detection of two quantitative trait loci (QTL) and, in total, identified 28 MTA, of which 19 were distributed on chromosome 7A (QTs.ipk-7A) and 1 on chromosome 7B (QTs.ipk-7B). The remaining eight MTA were unmapped and, therefore, no chromosomal and physical position was assigned to them (Table 2 and Table S2). Since QTs.ipk-7A harbored several MTA, only one marker with the highest log 10 ( P ) value and genotypic variance ( p G a d j ) —hereafter termed as a representative marker—was taken to represent the QTL. The representative markers of the QTL, viz., QTs.ipk-7A, and QTs.ipk-7B imparted p G a d j = 21.23 and 5.84%, respectively. The total p G a d j imparted by all MTA amounted to 25.79%.
As observed in the phenotypic data analyses, a highly significant negative correlation of Ptr infection was observed with the plant height, suggesting that taller plants escape the disease infestation. Nevertheless, our GWAS—albeit setting a liberal MTA detection threshold (FDR) of 0.20—did not identify the Rht genes. The FDR value for Rht-D1 was, however, 0.25 and, therefore, being close to the threshold and frequent (Rht-D1a = 0.41; Rht-D1b = 0.59) in the European wheat germplasm, we investigated its genetic/allelic influence on the tan spot. The findings concurred with the phenotypic analyses where the impact of Rht-D1a (wild-type; tall allele) was significantly greater than Rht-D1b (dwarfing allele; short allele) in terms of reducing the Ptr infection (Figure 1b and Figure 2a). Similarly, allele-wise phenotypic distribution showed a significant difference between the varieties harboring the reference (major) and variant (minor) allele of the representative marker for the 7A-QTL QTs.ipk-7A (gene-ID = TraesCS7A02G264300; Table 2, Figure 2b). The same was true for another small- or medium-effect locus QTs.ipk-7B on chromosome 7B (gene-ID = TraesCS7B02G444900; Table 2; Figure 2c). Since the physical interval of the QTs.ipk-7A is large (~36-Mb), it is difficult to identify a single causative gene. Nevertheless, the large effect of 7A-QTL explaining >20% genotypic variance merits its future use for gene cloning and downstream molecular and functional analyses.
The extent of linkage disequilibrium (LD), the non-random association between different loci, plays a vital role in GWAS. The panel under investigation has been previously examined for the LD via different marker platforms [48,49]. In addition, the population structure and related parameters have been published earlier [49]. Here, to observe the alleles’ distribution in the investigated germplasm, we performed PC analyses based on singular value decomposition, as described previously [49]. The first ten PCs accounted for 29.2% of the total variation (Figure 3a). A two-dimensional scatterplot of the first two PCs for Rht-D1 alleles showed a clear distribution of both alleles on the opposite sides of the central axis (Figure 3b). However, the large-effect QTL on chromosome 7A (i.e., QTs.ipk-7A)—the minor allele of which was present in only 6.7% of the varieties—showed no clear pattern (Table 2; Figure 3c).

3.3. Genome-Wide Prediction Studies Show That Marker Density, Marker Platform, and Genetic Models Do Not Substantially Influence the Prediction Accuracies

To observe the influence of individual marker platforms on the genome-wide prediction accuracies of Ptr infection, we tested three different models making different assumptions in this study, thus creating four scenarios as (1) SSR alleles, (2) 35k SNP array, (3) 90k SNP array, and (4) the full set of markers altogether. In every scenario, we incorporated the functional candidate-gene markers as well. The mean prediction accuracies resulting from the five-fold cross-validation scenario of Ptr infection generally produced similar results (~40%) across all three tested model scenarios, i.e., the GBLUP model that accounted for the main additive effects of markers assuming equal variances, BayesB by assuming unequal marker variances, and RKHSR that accounted for both additive and the epistatic interactions among the loci (Figure 4a–d). Overall, the 90k platform outperformed every other scenario with higher (1–2%) prediction accuracies. The RKHSR resulted in relatively better prediction accuracy than the GBLUP and BayesB, suggesting—albeit not highly prevalent—the presence of epistatic interactions for the Ptr infection.

4. Discussion

4.1. A Parallel Exploitation of Genetic Variation and Morphological Escape Traits Can Help Improve the Tan Spot Resistance in Wheat

A significant genetic variation for the traits under selection provides a substantial impetus in improving breeding programs’ genetic gains. However, especially for disease traits, besides the genotypic variation, the presence of a large and significant genotype-by-environment interaction is virtually a norm mainly because (1) the disease pressures are uneven across environments, and (2) the environmental effects are very unpredictable. We evaluated 372 registered wheat varieties in replicated field trials and observed significant genotypic variation for tan spot (Ptr) infection. However, due to large and significant genotype-by-environment interaction, we observed a moderate broad-sense heritability that amounted to 0.33. Recently, based on multiple environment trials, Juliana et al. [30] reported similar moderate broad-sense heritability estimates for the tan spot adult plant resistance in wheat.
Coupled with significant genetic variation, certain easy-to-score morphological traits have been purposed to escape disease infestations not only for the tan spot but also for other diseases, e.g., Fusarium head blight and Septoria tritici blotch [32,33,48]. We observed a highly significant negative correlation of plant height and a moderate negative but significant correlation of heading date with tan spot infection. Based on previous studies and this study, it seems that the major genes for plant height (Rht) or photoperiodism (Ppd) may show a pleiotropic effect on the disease traits.

4.2. The Influence of Rht-D1 and QTs.ipk-7A on the Tan Spot for Marker-Assisted Selection

We identified two significant QTL associated with tan spot on chromosomes 7A and 7B at the 263.18 and 709.08-Mb positions. Although previous studies on both bi-parental and diverse populations have reported tan spot-associated loci on chromosomes 7A and 7B, none of them resulted in identifying QTL imparting >20% of the genotypic variance [22,23,24,25,26,27,28]. Recently, Liu et al., [20] in a meta-QTL study, identified one QTL on chromosomes 7A (116.1–133.2-Mb) and two on 7B (21.0–34.0-Mb and 614.2–622.8-Mb). The physical distances of the QTL identified in our study from the meta-QTL are large and, given an extensive linkage disequilibrium in wheat, they may be considered novel. Also, the comparison of markers and their corresponding positions is not possible, mainly due to different marker systems and maps (physical and/or genetic). Marker-assisted selection (MAS) is profitable per unit time and cost only when the trait-tagged markers impart considerable genotypic variation. Therefore, due to sizeable genotypic variance, i.e., 21.23%, the 7A-QTL is of interest for MAS. The 7B-QTL explained 5.84% of the genotypic variance and can be considered a second target for MAS.
Besides, as shown in Figure 1b and Figure 2a, the functional marker for the candidate gene Rht-D1—although it did not pass the significance threshold—showed a relatively large effect on the tan spot infection score. This indicates that the MAS based on the Rht-D1a, i.e., the wild-type or tall allele for plant height, may help select for tan spot resistance. Since the tan spot infection is more lethal at later plant growth stages (e.g., BBCH-45–65) [2,4], the relatively taller plant selection should help escape the disease infestation. Consistent with this observation, the genotypes harboring Rht-D1a were more resistant than those bearing Rht-D1b (Figure 2a). Semi-dwarf or short-statured plants are, on the other hand, desired in breeding programs to achieve higher stem/stand strength. This warrants the use of genes other than Rht-D1a to tailor plant height. The frequency of Rht-B1b—as reported previously—is shallow in European varieties [34], which was perhaps why it was not identified as significantly associated with the tan spot in our study. Nonetheless, the selection of Rht-B1b to reduce height may not seem advisable given the similar effects of both genes on several other traits. For this purpose, other Rht-genes (e.g., Rht8 or Rht24) may be used to fine-tune the plant height for improved lodging resistance in breeding programs [50,51].

4.3. Genome-Wide Prediction Accuracy Reveals the Prospects of Genome-Wide Selection for Tan Spot Resistance

Improving qualitative disease resistance by selecting for or against major genes or QTL is a resource and time-efficient measure. However, most disease genes are responsive only against one or a few pathogen races and lack a broad-spectrum application [52]. Moreover, the practical difficulty and costs become co-extensive while pyramiding several QTL in an elite background and, therefore, inadvertently affect the breeding operations. Also, relying on only one or a few large-effect genes can result in the acceleration of pathogen evolution. For long-term sustainable genetic gains, improving quantitative resistance is deemed a durable strategy. Therefore, instead of concentrating on only large-effect loci, using the total genetic value predicted by both small- and large-effect loci helps select lines with relatively broad-spectrum resistance.
In our study, we observed a highly quantitative genetic nature of the tan spot where the significantly associated markers—albeit considerably improving marker density—in total, explained only 25.79% of the genotypic variance. Markers that do not cross significance thresholds in GWAS, e.g., Rht-D1, are usually not used for MAS. Therefore, instead of concentrating only on large-effect loci, genome-wide prediction (GP) of the total genetic value of tan spot based on small- and large-effect markers is a holistic tactic to improve the broad-spectrum resistance. We evaluated the GP accuracy for tan spot resistance by modeling the loci’s additive effects assuming equal variances, unequal variances, and epistatic interaction. In line with a previous study, the mean GP accuracies calculated across 100 cycles and five-fold cross-validation scenarios amounted to ~40%, with virtually no statistically significant difference between the models [30]. Although epistatic interactions were previously reported to be pervasive in self-pollinating species like wheat [53], we observed only a slight increase in the prediction accuracy by modeling epistatic interactions.
Based on the hypothesis that marker platform and thus marker density influences the outcome of GP accuracy, we evaluated all three models on all marker platforms, viz., SSR, 35k, and 90k SNP arrays individually and marker loci combined from all platforms. However, increasing the marker density did not result in any significant increase in GP accuracy. This finding is in line with previous reports where GP accuracy was not influenced above a certain number of markers [31,32], underlining that all marker platforms are almost equally efficient to predict tan spot infection. In practical breeding, nevertheless, the usefulness of GP may be hampered by shifts in the virulence spectrum of the pathogen in different environments or breeding target zones.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/genes12040490/s1, Figure S1: Distribution and correlation of tan spot (Ptr) resistance with plant height and heading date, Figure S2: Genomic relationship matrix, Table S1: List of genotypes, BLUEs, and associated markers, Table S2: List of tan spot resistance associated QTL.

Author Contributions

Conceptualization, Q.H.M. and M.S.R.; methodology, Q.H.M.; formal analysis, Q.H.M. and R.K.; visualization, Q.H.M.; genotypic data curation, M.W.G.; phenotypic data curation, B.R.; writing—original draft preparation, Q.H.M.; editing, Q.H.M., R.K., V.M., and J.C.R.; funding acquisition, M.S.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded within the framework of the project GABI-Wheat by the Federal Ministry of Education and Research (BMBL), grant number 0315067.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article and supplementary files.

Acknowledgments

We thank A. Fliefer, P. Jaquim, K. Wendehake, and J. Plieske for genotyping the varieties. We are grateful to three anonymous reviews whose comments helped to improve the manuscript.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results. Q.H.M. (presently), V.M., and M.W.G. are members of various companies. However, this does limit the availability or sharing of data and materials.

References

  1. McMullen, M.; Adhikari, T. Fungal Leaf Spot Diseases of Wheat: Tan Spot, Stagonospora nodorum Blotch and Septoria tritici Blotch; North Dakota State University: Fargo, ND, USA, 2009. [Google Scholar]
  2. De Wolf, E.; Effertz, R.; Ali, S.; Francl, L. Vistas of tan spot research. Can. J. Plant Pathol. 1998, 20, 349–370. [Google Scholar] [CrossRef]
  3. Rees, R.; Platz, G.; Mayer, R. Yield losses in wheat from yellow spot: Comparison of estimates derived from single tillers and plots. Aust. J. Agric. Res. 1982, 33, 899–908. [Google Scholar] [CrossRef]
  4. Rees, R.; Platz, G. Effects of yellow spot on wheat: Comparison of epidemics at different stages of crop development. Aust. J. Agric. Res. 1983, 34, 39–46. [Google Scholar] [CrossRef]
  5. Shabeer, A.; Bockus, W. Tan spot effects on yield and yield components relative to growth stage in winter wheat. Plant Dis. 1988, 72, 599–602. [Google Scholar] [CrossRef]
  6. Rees, R.; Platz, G. The occurrence and control of yellow spot of wheat in north-eastern Australia. Aust. J. Exp. Agric. 1979, 19, 369–372. [Google Scholar] [CrossRef]
  7. Bockus, W.W.; Claassen, M.M. Effects of crop rotation and residue management practices on severity of tan spot of winter wheat. Plant Dis. 1992, 76, 633–636. [Google Scholar] [CrossRef]
  8. Sutton, J.; Vyn, T. Crop sequences and tillage practices in relation to diseases of winter wheat in Ontario. Can. J. Plant Pathol. 1990, 12, 358–368. [Google Scholar] [CrossRef]
  9. Ciuffetti, L.M.; Manning, V.A.; Pandelova, I.; Betts, M.F.; Martinez, J.P. Host-selective toxins, Ptr ToxA and Ptr ToxB, as necrotrophic effectors in the Pyrenophora tritici-repentis—Wheat interaction. New Phytol. 2010, 187, 911–919. [Google Scholar] [CrossRef]
  10. Faris, J.D.; Liu, Z.; Xu, S.S. Genetics of tan spot resistance in wheat. Theor. Appl. Genet. 2013, 126, 2197–2217. [Google Scholar] [CrossRef] [PubMed]
  11. Strelkov, S.; Lamari, L. Host–parasite interactions in tan spot [Pyrenophora tritici-repentis] of wheat. Can. J. Plant Pathol. 2003, 25, 339–349. [Google Scholar] [CrossRef]
  12. Faris, J.; Anderson, J.A.; Francl, L.; Jordahl, J. Chromosomal location of a gene conditioning insensitivity in wheat to a necrosis-inducing culture filtrate from Pyrenophora tritici-repentis. Phytopathology 1996, 86, 459–463. [Google Scholar] [CrossRef]
  13. Effertz, R.; Anderson, J.; Francl, L. Restriction fragment length polymorphism mapping of resistance to two races of Pyrenophora tritici-repentis in adult and seedling wheat. Phytopathology 2001, 91, 572–578. [Google Scholar] [CrossRef] [PubMed]
  14. Abeysekara, N.S.; Friesen, T.L.; Liu, Z.; McClean, P.E.; Faris, J.D. Marker development and saturation mapping of the tan spot Ptr ToxB sensitivity locus Tsc2 in hexaploid wheat. Plant Genome 2010, 3. [Google Scholar] [CrossRef]
  15. Faris, J.D.; Zhang, Z.; Lu, H.; Lu, S.; Reddy, L.; Cloutier, S.; Fellers, J.P.; Meinhardt, S.W.; Rasmussen, J.B.; Xu, S.S. A unique wheat disease resistance-like gene governs effector-triggered susceptibility to necrotrophic pathogens. Proc. Natl. Acad. Sci. USA 2010, 107, 13544–13549. [Google Scholar] [CrossRef]
  16. Singh, P.; Gonzalez-Hernandez, J.; Mergoum, M.; Ali, S.; Adhikari, T.; Kianian, S.; Elias, E.; Hughes, G. Identification and molecular mapping of a gene conferring resistance to Pyrenophora tritici-repentis race 3 in tetraploid wheat. Phytopathology 2006, 96, 885–889. [Google Scholar] [CrossRef]
  17. Tadesse, W.; Hsam, S.L.; Wenzel, G.; Zeller, F.J. Identification and monosomic analysis of tan spot resistance genes in synthetic wheat lines (Triticum turgidum L. × Aegilops tauschii Coss.). Crop Sci. 2006, 46, 1212–1217. [Google Scholar] [CrossRef]
  18. Tadesse, W.; Hsam, S.; Zeller, F. Evaluation of common wheat cultivars for tan spot resistance and chromosomal location of a resistance gene in the cultivar ‘Salamouni’. Plant Breed. 2006, 125, 318–322. [Google Scholar] [CrossRef]
  19. Singh, P.; Mergoum, M.; Gonzalez-Hernandez, J.; Ali, S.; Adhikari, T.; Kianian, S.; Elias, E.; Hughes, G. Genetics and molecular mapping of resistance to necrosis inducing race 5 of Pyrenophora tritici-repentis in tetraploid wheat. Mol. Breed. 2008, 21, 293–304. [Google Scholar] [CrossRef]
  20. Liu, Y.; Salsman, E.; Wang, R.; Galagedara, N.; Zhang, Q.; Fiedler, J.D.; Liu, Z.; Xu, S.; Faris, J.D.; Li, X. Meta-QTL analysis of tan spot resistance in wheat. Theor. Appl. Genet. 2020, 133, 2363–2375. [Google Scholar] [CrossRef]
  21. Gurung, S.; Mamidi, S.; Bonman, J.; Jackson, E.; Del Rio, L.; Acevedo, M.; Mergoum, M.; Adhikari, T. Identification of novel genomic regions associated with resistance to Pyrenophora tritici-repentis races 1 and 5 in spring wheat landraces using association analysis. Theor. Appl. Genet. 2011, 123, 1029–1041. [Google Scholar] [CrossRef] [PubMed]
  22. Patel, J.S.; Mamidi, S.; Bonman, J.M.; Adhikari, T.B. Identification of QTL in spring wheat associated with resistance to a novel isolate of Pyrenophora tritici-repentis. Crop Sci. 2013, 53, 842–852. [Google Scholar] [CrossRef]
  23. Kollers, S.; Rodemann, B.; Ling, J.; Korzun, V.; Ebmeyer, E.; Argillier, O.; Hinze, M.; Plieske, J.; Kulosa, D.; Ganal, M.W. Genome-wide association mapping of tan spot resistance (Pyrenophora tritici-repentis) in European winter wheat. Mol. Breed. 2014, 34, 363–371. [Google Scholar] [CrossRef]
  24. Liu, Z.; El-Basyoni, I.; Kariyawasam, G.; Zhang, G.; Fritz, A.; Hansen, J.; Marais, F.; Friskop, A.; Chao, S.; Akhunov, E.; et al. Evaluation and association mapping of resistance to tan spot and Stagonospora nodorum blotch in adapted winter wheat germplasm. Plant Dis. 2015, 99, 1333–1341. [Google Scholar] [CrossRef]
  25. Juliana, P.; Singh, R.P.; Singh, P.K.; Poland, J.A.; Bergstrom, G.C.; Huerta-Espino, J.; Bhavani, S.; Crossa, J.; Sorrells, M.E. Genome-wide association mapping for resistance to leaf rust, stripe rust and tan spot in wheat reveals potential candidate genes. Theor. Appl. Genet. 2018, 131, 1405–1422. [Google Scholar] [CrossRef] [PubMed]
  26. Dinglasan, E.G.; Singh, D.; Shankar, M.; Afanasenko, O.; Platz, G.; Godwin, I.D.; Voss-Fels, K.P.; Hickey, L.T. Discovering new alleles for yellow spot resistance in the Vavilov wheat collection. Theor. Appl. Genet. 2019, 132, 149–162. [Google Scholar] [CrossRef] [PubMed]
  27. Galagedara, N.; Liu, Y.; Fiedler, J.; Shi, G.; Chiao, S.; Xu, S.S.; Faris, J.D.; Li, X.; Liu, Z. Genome-wide association mapping of tan spot resistance in a worldwide collection of durum wheat. Theor. Appl. Genet. 2020, 133, 2227–2237. [Google Scholar] [CrossRef] [PubMed]
  28. Kokhmetova, A.; Sehgal, D.; Ali, S.; Atishova, M.; Kumarbayeva, M.; Leonova, I.; Dreisigacker, S. Genome-wide association study of tan spot resistance in a hexaploid wheat collection from Kazakhstan. Front. Genet. 2020, 11, 581214. [Google Scholar] [CrossRef]
  29. Meuwissen, T.H.; Hayes, B.J.; Goddard, M.E. Prediction of total genetic value using genome-wide dense marker maps. Genetics 2001, 157, 1819–1829. [Google Scholar]
  30. Juliana, P.; Singh, R.P.; Singh, P.K.; Crossa, J.; Rutkoski, J.E.; Poland, J.A.; Bergstrom, G.C.; Sorrells, M.E. Comparison of models and whole-genome profiling approaches for genomic-enabled prediction of Septoria tritici blotch, Stagonospora nodorum blotch, and tan spot resistance in wheat. Plant Genome 2017, 10. [Google Scholar] [CrossRef]
  31. Jiang, Y.; Zhao, Y.; Rodemann, B.; Plieske, J.; Kollers, S.; Korzun, V.; Ebmeyer, E.; Argillier, O.; Hinze, M.; Ling, J.; et al. Potential and limits to unravel the genetic architecture and predict the variation of Fusarium head blight resistance in European winter wheat (Triticum aestivum L.). Heredity 2015, 114, 318–326. [Google Scholar] [CrossRef]
  32. Muqaddasi, Q.H.; Zhao, Y.; Rodemann, B.; Plieske, J.; Ganal, M.W.; Röder, M.S. Genome-wide association mapping and prediction of adult stage Septoria tritici blotch infection in European winter wheat via high-density marker arrays. Plant Genome 2019, 12, 180029. [Google Scholar] [CrossRef]
  33. Srinivasachary; Gosman, N.; Steed, A.; Hollins, T.; Bayles, R.; Jennings, P.; Nicholson, P. Semi-dwarfing Rht-B1 and Rht-D1 loci of wheat differ significantly in their influence on resistance to Fusarium head blight. Theor. Appl. Genet. 2009, 118, 695. [Google Scholar] [CrossRef]
  34. Zanke, C.D.; Ling, J.; Plieske, J.; Kollers, S.; Ebmeyer, E.; Korzun, V.; Argillier, O.; Stiewe, G.; Hinze, M.; Neumann, K.; et al. Whole genome association mapping of plant height in winter wheat (Triticum aestivum L.). PLoS ONE 2014, 9, e113287. [Google Scholar] [CrossRef]
  35. Zanke, C.; Ling, J.; Plieske, J.; Kollers, S.; Ebmeyer, E.; Korzun, V.; Argillier, O.; Stiewe, G.; Hinze, M.; Beier, S.; et al. Genetic architecture of main effect QTL for heading date in European winter wheat. Front. Plant Sci. 2014, 5, 217. [Google Scholar] [CrossRef] [PubMed]
  36. Sorrells, M.E.; Gustafson, J.P.; Somers, D.; Chao, S.; Benscher, D.; Guedira-Brown, G.; Huttner, E.; Kilian, A.; McGuire, P.E.; Ross, K.; et al. Reconstruction of the Synthetic W7984 × Opata M85 wheat reference population. Genome 2011, 54, 875–882. [Google Scholar] [CrossRef]
  37. Sun, C.; Dong, Z.; Zhao, L.; Ren, Y.; Zhang, N.; Chen, F. The Wheat 660K SNP array demonstrates great potential for marker-assisted selection in polyploid wheat. Plant Biotechnol. J. 2020, 18, 1354–1360. [Google Scholar] [CrossRef] [PubMed]
  38. Yu, J.; Pressoir, G.; Briggs, W.H.; Bi, I.V.; Yamasaki, M.; Doebley, J.F.; McMullen, M.D.; Gaut, B.S.; Nielsen, D.M.; Holland, J.B.; et al. A unified mixed-model method for association mapping that accounts for multiple levels of relatedness. Nat. Genet. 2006, 38, 203. [Google Scholar] [CrossRef] [PubMed]
  39. VanRaden, P.M. Efficient methods to compute genomic predictions. J. Dairy Sci. 2008, 91, 4414–4423. [Google Scholar] [CrossRef]
  40. Benjamini, Y.; Hochberg, Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B Methodol. 1995, 57, 289–300. [Google Scholar] [CrossRef]
  41. Utz, H.F.; Melchinger, A.E.; Schön, C.C. Bias and sampling error of the estimated proportion of genotypic variance explained by quantitative trait loci determined from experimental data in maize using cross validation and validation with independent samples. Genetics 2000, 154, 1839–1849. [Google Scholar]
  42. Habier, D.; Fernando, R.L.; Dekkers, J.C. The impact of genetic relationship information on genome-assisted breeding values. Genetics 2007, 177, 2389–2397. [Google Scholar] [CrossRef]
  43. Pérez, P.; de Los Campos, G. Genome-wide regression and prediction with the BGLR statistical package. Genetics 2014, 198, 483–495. [Google Scholar] [CrossRef]
  44. Gianola, D.; Fernando, R.L.; Stella, A. Genomic-assisted prediction of genetic value with semiparametric procedures. Genetics 2006, 173, 1761–1776. [Google Scholar] [CrossRef] [PubMed]
  45. Team, R.C. R: A Language and Environment for Statistical Computing; R Core Team: Diepoldsau, Switzerland, 2013. [Google Scholar]
  46. Bates, D.; Mächler, M.; Bolker, B.; Walker, S. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 2015, 67, 1–48. [Google Scholar] [CrossRef]
  47. Endelman, J.B. Ridge regression and other kernels for genomic selection with R package rrBLUP. Plant Genome 2011, 4, 200–255. [Google Scholar] [CrossRef]
  48. Kollers, S.; Rodemann, B.; Ling, J.; Korzun, V.; Ebmeyer, E.; Argillier, O.; Hinze, M.; Plieske, J.; Kulosa, D.; Ganal, M.W.; et al. Whole genome association mapping of Fusarium head blight resistance in European winter wheat (Triticum aestivum L.). PLoS ONE 2013, 8, e57500. [Google Scholar] [CrossRef]
  49. Muqaddasi, Q.H.; Brassac, J.; Koppolu, R.; Plieske, J.; Ganal, M.W.; Röder, M.S. TaAPO-A1, an ortholog of rice ABERRANT PANICLE ORGANIZATION 1, is associated with total spikelet number per spike in elite European hexaploid winter wheat (Triticum aestivum L.) varieties. Sci. Rep. 2019, 9, 1–12. [Google Scholar] [CrossRef]
  50. Korzun, V.; Röder, M.; Ganal, M.; Worland, A.; Law, C. Genetic analysis of the dwarfing gene (Rht8) in wheat. Part I. Molecular mapping of Rht8 on the short arm of chromosome 2D of bread wheat (Triticum aestivum L.). Theor. Appl. Genet. 1998, 96, 1104–1109. [Google Scholar] [CrossRef]
  51. Würschum, T.; Langer, S.M.; Longin, C.F.H.; Tucker, M.R.; Leiser, W.L. A modern Green Revolution gene for reduced height in wheat. Plant J. 2017, 92, 892–903. [Google Scholar] [CrossRef] [PubMed]
  52. Figueroa, M.; Hammond-Kosack, K.E.; Solomon, P.S. A review of wheat diseases—A field perspective. Mol. Plant Pathol. 2018, 19, 1523–1536. [Google Scholar] [CrossRef] [PubMed]
  53. Heslot, N.; Yang, H.P.; Sorrells, M.E.; Jannink, J.L. Genomic selection in plant breeding: A comparison of models. Crop Sci. 2012, 52, 146–160. [Google Scholar] [CrossRef]
Figure 1. Phenotypic distribution and summary of the genome-wide association studies (GWAS) of tan spot infection in European winter wheat varieties. (a) Histogram of the square-root transformed tan spot Pyrenophora tritici-repentis (Ptr) infection scores on wheat varieties. (b) Manhattan plot showing the distribution of marker significance [ log 10 ( P values ) ] along the wheat chromosomes. The correction for population stratification was performed by using the first three principal components ( P C [ 1 3 ] ) and an additive relationship matrix (G) in the linear mixed-effect model. The red dashed line marks the false discovery rate (FDR) threshold ( α = 0.20 ) to detect marker-trait associations. (c) Quantile-quantile (qq) plot showing the distribution of observed versus expected (red dashed line) log 10 ( P value ) based on the naïve model (the general linear model without correction for population structure), P C [ 1 3 ] model (the population structure corrected with first three PCs), the G model (population structure corrected with a genomic relationship matrix), and the P C [ 1 3 ] + G model (population structure corrected with both PCs and the G matrix). The pink highlighted markers in the Manhattan plot designate the representative markers. The color code of different models is given in the figure legend. n = the number of varieties; p = significance value of the Shapiro–Wilk normality test; H 2 = broad-sense heritability; P = the number of quality marker loci; unm = the unmapped markers.
Figure 1. Phenotypic distribution and summary of the genome-wide association studies (GWAS) of tan spot infection in European winter wheat varieties. (a) Histogram of the square-root transformed tan spot Pyrenophora tritici-repentis (Ptr) infection scores on wheat varieties. (b) Manhattan plot showing the distribution of marker significance [ log 10 ( P values ) ] along the wheat chromosomes. The correction for population stratification was performed by using the first three principal components ( P C [ 1 3 ] ) and an additive relationship matrix (G) in the linear mixed-effect model. The red dashed line marks the false discovery rate (FDR) threshold ( α = 0.20 ) to detect marker-trait associations. (c) Quantile-quantile (qq) plot showing the distribution of observed versus expected (red dashed line) log 10 ( P value ) based on the naïve model (the general linear model without correction for population structure), P C [ 1 3 ] model (the population structure corrected with first three PCs), the G model (population structure corrected with a genomic relationship matrix), and the P C [ 1 3 ] + G model (population structure corrected with both PCs and the G matrix). The pink highlighted markers in the Manhattan plot designate the representative markers. The color code of different models is given in the figure legend. n = the number of varieties; p = significance value of the Shapiro–Wilk normality test; H 2 = broad-sense heritability; P = the number of quality marker loci; unm = the unmapped markers.
Genes 12 00490 g001
Figure 2. Allelic influence of the tan spot-associated quantitative trait loci (QTL) in European wheat. (a) Distribution of the best linear unbiased estimations (BLUEs; left panel) to compare the allele-wise distribution of the Rht-D1 alleles. (b) Allelic distribution of the representative marker for tan spot-associated locus QTs.ipk-7A. (c) Allelic distribution of the representative marker for tan spot-associated locus QTs.ipk-7B. The first, second, and third rows in the figure header show the QTL name, the gene-ID corresponding to the most significant marker in the QTL, and the QTL’s chromosome and physical position of the representative marker. The x-axis shows the representative marker names and their alleles. n = number of varieties harbored by the corresponding panel; P = significance value of Welch two-sample t-test; R = reference (major) allele; and V = variant (minor) allele.
Figure 2. Allelic influence of the tan spot-associated quantitative trait loci (QTL) in European wheat. (a) Distribution of the best linear unbiased estimations (BLUEs; left panel) to compare the allele-wise distribution of the Rht-D1 alleles. (b) Allelic distribution of the representative marker for tan spot-associated locus QTs.ipk-7A. (c) Allelic distribution of the representative marker for tan spot-associated locus QTs.ipk-7B. The first, second, and third rows in the figure header show the QTL name, the gene-ID corresponding to the most significant marker in the QTL, and the QTL’s chromosome and physical position of the representative marker. The x-axis shows the representative marker names and their alleles. n = number of varieties harbored by the corresponding panel; P = significance value of Welch two-sample t-test; R = reference (major) allele; and V = variant (minor) allele.
Genes 12 00490 g002
Figure 3. Principal component (PC) analysis on the wheat marker loci combined from the SSR alleles, functional candidate-gene markers, 35k, and 90k single nucleotide polymorphism arrays. (a) Scree plot showing first ten PCs and their corresponding proportion of variance. (b) Two-dimensional scatterplot showing the absence of pronounced clustering among the varieties, except those based on the Rht-D1 alleles. (c) Scatterplot showing the lack of clustering among the varieties based on the QTs.ipk-7A representative marker alleles. (d) Scatterplot showing the absence of clustering among the varieties based on the QTs.ipk-7B representative marker alleles. Color codes are given in the respective sub-figure’s legend/s. n and p denote the number of varieties and the marker genotypes used in the analysis, respectively.
Figure 3. Principal component (PC) analysis on the wheat marker loci combined from the SSR alleles, functional candidate-gene markers, 35k, and 90k single nucleotide polymorphism arrays. (a) Scree plot showing first ten PCs and their corresponding proportion of variance. (b) Two-dimensional scatterplot showing the absence of pronounced clustering among the varieties, except those based on the Rht-D1 alleles. (c) Scatterplot showing the lack of clustering among the varieties based on the QTs.ipk-7A representative marker alleles. (d) Scatterplot showing the absence of clustering among the varieties based on the QTs.ipk-7B representative marker alleles. Color codes are given in the respective sub-figure’s legend/s. n and p denote the number of varieties and the marker genotypes used in the analysis, respectively.
Genes 12 00490 g003
Figure 4. Accuracy of genome-wide prediction (GP) of tan spot (Ptr) infection in wheat. The figure header represents four different GP scenarios, viz., (a) GP based on SSR alleles and candidate genes, (b) GP based on quality 35k SNPs and candidate genes, (c) GP based on quality 90k SNPs and candidate genes, and (d) GP based on markers combined from every platform. The GP accuracy assessment is based on three models, viz., genomic best linear unbiased prediction (GBLUP), Bayesian alphabet B (BayesB), and reproducing kernel Hilbert space regression (RKHSR). The GP accuracies were evaluated through 100 random five-fold cross-validation cycles. Symbols μ and σ denote the mean accuracy and standard deviation of the corresponding model, respectively.
Figure 4. Accuracy of genome-wide prediction (GP) of tan spot (Ptr) infection in wheat. The figure header represents four different GP scenarios, viz., (a) GP based on SSR alleles and candidate genes, (b) GP based on quality 35k SNPs and candidate genes, (c) GP based on quality 90k SNPs and candidate genes, and (d) GP based on markers combined from every platform. The GP accuracy assessment is based on three models, viz., genomic best linear unbiased prediction (GBLUP), Bayesian alphabet B (BayesB), and reproducing kernel Hilbert space regression (RKHSR). The GP accuracies were evaluated through 100 random five-fold cross-validation cycles. Symbols μ and σ denote the mean accuracy and standard deviation of the corresponding model, respectively.
Genes 12 00490 g004
Table 1. ANOVA for tan spot in European wheat varieties.
Table 1. ANOVA for tan spot in European wheat varieties.
DfSum SqMean SqF-ValuePr(>F)Sig. σ 2 SD
Genotype371113.030.31.496.48 × 10−5***0.0500.224
Env.1260.42260.41273.89<2.00 × 10−16***0.7000.836
Residuals37175.840.2 0.2040.452
Df = degree of freedom; Sq = squares; Sig. = significance code; σ 2 = variance; SD = standard deviation; Env. = environment; *** = significant at the 0.001 probability (P) level.
Table 2. The quantitative trait loci (QTL) and markers associated with the wheat tan spot infection from the combined set of candidate genes, simple sequence repeats (SSRs), 35k, and 90k marker platforms.
Table 2. The quantitative trait loci (QTL) and markers associated with the wheat tan spot infection from the combined set of candidate genes, simple sequence repeats (SSRs), 35k, and 90k marker platforms.
QTLMarkerChr.Pos. (bp)|log10(P)|MAF p G a d j
QTs.ipk-7AEx_c37521_6707A2462583333.780.0618.57
QTs.ipk-7Awsnp_Ra_c12708_202814397A2462584373.780.0618.57
QTs.ipk-7ABS00067759_517A2476627563.780.0618.57
QTs.ipk-7AAX_944208107A2476627583.820.0618.75
QTs.ipk-7AEx_c6348_12057A2492798813.780.0618.57
QTs.ipk-7AExcalibur_c15904_13317A2500371453.780.0618.57
QTs.ipk-7Awsnp_Ex_c4819_86006187A2521761343.780.0618.57
QTs.ipk-7ARAC875_c109483_5237A2543869613.780.0618.57
QTs.ipk-7Awsnp_Ex_c26560_358032107A2543869613.780.0618.57
QTs.ipk-7Awsnp_be471272A_Ta_2_17A2609056263.780.0618.57
QTs.ipk-7AExcalibur_rep_c67182_5847A2631772874.250.0721.23
QTs.ipk-7ATdurum_contig10174_1557A2631772874.250.0721.23
QTs.ipk-7ARAC875_rep_c81362_3447A2723395443.780.0618.57
QTs.ipk-7AAX_947067827A2729977983.820.0618.75
QTs.ipk-7Awsnp_Ex_rep_c102317_875126607A2729977983.780.0618.57
QTs.ipk-7ABS00036422_517A2755456433.780.0618.57
QTs.ipk-7Awsnp_Ra_c13009_206907357A2759197163.780.0618.57
QTs.ipk-7Awsnp_Ex_c5330_94221067A2759204693.780.0618.57
QTs.ipk-7ABS00066695_517A2823685773.780.0618.57
QTs.ipk-7BAX_946525617B7090824224.070.405.84
Chr. = chromosome name; Pos. (bp) = physical position of the corresponding marker in base-pairs; |log10(P)| = negative log transformed significance (P) value of the corresponding marker; MAF = minor allele frequency; p G a d j = percentage of adjusted genotypic variance imparted by the corresponding marker. The bold markers are the representative QTL descrbied in the text.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Muqaddasi, Q.H.; Kamal, R.; Mirdita, V.; Rodemann, B.; Ganal, M.W.; Reif, J.C.; Röder, M.S. Genome-Wide Association Studies and Prediction of Tan Spot (Pyrenophora tritici-repentis) Infection in European Winter Wheat via Different Marker Platforms. Genes 2021, 12, 490. https://doi.org/10.3390/genes12040490

AMA Style

Muqaddasi QH, Kamal R, Mirdita V, Rodemann B, Ganal MW, Reif JC, Röder MS. Genome-Wide Association Studies and Prediction of Tan Spot (Pyrenophora tritici-repentis) Infection in European Winter Wheat via Different Marker Platforms. Genes. 2021; 12(4):490. https://doi.org/10.3390/genes12040490

Chicago/Turabian Style

Muqaddasi, Quddoos H., Roop Kamal, Vilson Mirdita, Bernd Rodemann, Martin W. Ganal, Jochen C. Reif, and Marion S. Röder. 2021. "Genome-Wide Association Studies and Prediction of Tan Spot (Pyrenophora tritici-repentis) Infection in European Winter Wheat via Different Marker Platforms" Genes 12, no. 4: 490. https://doi.org/10.3390/genes12040490

APA Style

Muqaddasi, Q. H., Kamal, R., Mirdita, V., Rodemann, B., Ganal, M. W., Reif, J. C., & Röder, M. S. (2021). Genome-Wide Association Studies and Prediction of Tan Spot (Pyrenophora tritici-repentis) Infection in European Winter Wheat via Different Marker Platforms. Genes, 12(4), 490. https://doi.org/10.3390/genes12040490

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop