Next Article in Journal
Evaluation of Siderophores Generated by Pseudomonas Bacteria and Their Possible Application as Fe Biofertilizers
Next Article in Special Issue
Identification of Novel Loci Precisely Modulating Pre-Harvest Sprouting Resistance and Red Color Components of the Seed Coat in T. aestivum L.
Previous Article in Journal
Nitric Oxide (NO) Improves Wheat Growth under Dehydration Conditions by Regulating Phytohormone Levels and Induction of the Expression of the TADHN Dehydrin Gene
Previous Article in Special Issue
Analysis of the Structural Organization and Expression of the Vrn-D1 Gene Controlling Growth Habit (Spring vs. Winter) in Aegilops tauschii Coss.
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 12
Issue 23
10.3390/plants12234053
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Diversity and Genome-Wide Association Study for the Phenology Response of Winter Wheats of North America, Western Asia, and Europe

by
Adil El Baouchi
1,2,*,
Mohammed Ibriz
2,
Susanne Dreisigacker
3,
Marta S. Lopes
4,5 and
Miguel Sanchez Garcia
1,*
1
International Center for Agricultural Research in the Dry Areas (ICARDA), Rabat 10100, Morocco
2
Genetic and Biometric Laboratory, Faculty of Sciences, Ibn Tofail University, Kenitra BP. 242, Morocco
3
The International Maize and Wheat Improvement Center (CIMMYT), Texcoco 56237, Mexico
4
The International Maize and Wheat Improvement Center (CIMMYT), Ankara 3906511, Turkey
5
Sustainable Field Crops, Institute for Food and Agricultural Research and Technology (IRTA), 251981 Lleida, Spain
*
Authors to whom correspondence should be addressed.
Plants 2023, 12(23), 4053; https://doi.org/10.3390/plants12234053
Submission received: 7 October 2023 / Revised: 14 November 2023 / Accepted: 15 November 2023 / Published: 1 December 2023
(This article belongs to the Special Issue Cereal Genetics and Molecular Genetics)
Browse Figures
Review Reports Versions Notes

Abstract

:
Wheat is a staple food in many areas around the World. In the 20th century, breeders and scientists were able to boost wheat yield considerably. However, a yield plateau has become a concern and is threatening food security. Investments in cutting-edge technologies, including genomics and precision phenology measurements, can provide valuable tools to drive crop improvement. The objectives of this study were to (i) investigate the genetic diversity in a set of winter wheat lines, (ii) characterize their phenological response under different vernalization and photoperiod conditions, and (iii) identify effective markers associated with the phenological traits. A total of 249 adapted genotypes of different geographical origin were genotyped using the 35K Axiom® Wheat Breeder’s Array. A total of 11,476 SNPs were used for genetic analysis. The set showed an average polymorphism information content of 0.37 and a genetic diversity of 0.43. A population structure analysis revealed three distinct subpopulations mainly related to their geographical origin (Europe, North America, and Western Asia). The lines of CGIAR origin showed the largest diversity and the lowest genetic distance to all other subpopulations. The phenology of the set was studied under controlled conditions using four combinations of long (19 h light) and short photoperiod (13 h light) and long vernalization (49 days at 5 °C) and no vernalization. With this, phenological traits such as earliness per se (Eps), relative response to vernalization (RRV), and relative response to photoperiod (RRP) were calculated. The phenotypic variation of growing degree days was significant in all phenology combinations. RRV ranged from 0 to 0.56, while RRP was higher with an overall average of 0.25. The GWAS analysis detected 30 marker-trait associations linked to five phenological traits. The highest significant marker was detected on chromosome 2D with a value of −log10(p) = 11.69. Only four loci known to regulate flowering exceeded the Bonferroni correction threshold of −log10(p) > 5.1. These results outline a solid foundation to address global food security and offer tremendous opportunities for advancing crop improvement strategies.

1. Introduction

Bread wheat (Triticum aestivum L.) is a vital crop for food security with a global production of 760 M tons [1] annually. Wheat also occupies the largest land under cultivation with more than 200 million ha planted and is considered as one of humanity’s key sources of calories, contributing to up to 20% of its daily intake worldwide. To fulfill the demand for human food consumption by 2050, global wheat production need to increase by almost 70% [2,3]. Several studies show the importance of genetic diversity for breeding programs to be successful and highlight the importance of adaptive genes in the improvement of wheat [4,5,6,7,8,9]. Given the impact of climate change, breeding for earlier or later maturing lines as well as optimizing other phenology traits becomes crucial to ensure sustainable wheat adaptability and production. Continuous selection for genes or ideotypes can lead to the decrease and drift in allele frequency and therefore to a loss of genetic diversity. In breeding populations, genomic region consistently correlated with important traits reflect major drivers of genome divergence [10]. Ref. [11] noted that in modern wheat, a restricted number of loci reflect the phenotypic improvement, while the largest proportion of the genome remains unchanged. Vernalization, photoperiod sensitivity, and earliness per se are regulated by Vrn, Ppd, and Eps genes, respectively, and influence the capacity of wheat to adapt to a particular cropping region. All genes are key regulators of the shift from the vegetative stage to the reproductive development when the temperature is warm enough and rainfall is sufficient. Insights in the genetic variation for these genes (Vrn, Ppd, and Eps) in winter wheat provides breeding programs with a source of alleles to improve adaptation given the current requirement for adaptability to changing weather patterns.
Vernalization is a crucial process protecting the floral meristems from frost damage [12]. The process is controlled by four homoeologous loci: Vrn-1 [13], Vrn-2 [14], Vrn-3 [15], and Vrn-4 [16]. Winter wheat varieties carry recessive vernalization alleles [17]. After the vernalization requirement is met, the photoperiod response, which is primarily influenced by the Ppd-1 loci, governs and adjusts flowering timing for specific environments [18]. Earliness per se or “narrow-sense earliness” is a quantitative trait detected when vernalization and photoperiod requirements are satisfied [19]. The mayor genes at the Vrn and Ppd loci have been cloned and the effect of diverse alleles studied. Gene-specific markers established are commonly used in wheat breeding programs [20]. Conducting phenological studies under diverse controlled growth conditions presents an alternative for characterizing adaptability, particularly when simulating field conditions [21].
Genome-wide association studies (GWAS) have become an extremely popular and cost-effective approach for the identification of molecular marker related to complex traits of interest [22,23,24,25,26]. The further understanding of phenological traits under controlled conditions and the identification of additional molecular markers are beneficial for future wheat improvement. GWAS using winter wheat germplasm has already been employed to examine the genetic bases of several traits [27,28,29]. GWAS can also be used to validate and explore allele effects of known candidate genes [30]. The primary goal of this research was to (i) explore the genetic diversity and population structure of the winter wheat panel WWAGI, (ii) examine the phenological response of the panel under different controlled vernalization and photoperiod treatments, and (iii) validate and identify molecular markers associated with phenological traits.

2. Results

2.1. Genetic Diversity, Relationship, and Population Differentiation

A total of 11,476 polymorphic SNPs were used for genetic analyzes. The B genome held most polymorphic SNPs (50%), followed by the A (40%) and D genome (10%). The chromosome with most SNPs was chromosome 5B (815), while the lowest number of SNPs was noted for chromosomes 4D (32). The average value of genetic diversity (GD) was 0.43 and ranged from 0.38 (chromosome 2A) to 0.47 (chromosome 3B) (Table 1).
The PIC values ranged from 0.33 of chromosome 2A to 0.4 of 3B, with an average of 0.37. Lines from Western Asia exhibited the highest GD with a value of 0.42, followed by lines originating from the CGIAR and lines from North America with a value of 0.39. The lowest GD was observed within lines from Central Asia, which amounted to 0.27. Regarding PIC, the same pattern was identified, ranging from 0.22 for lines from Central Asia to 0.37 for lines derived from Western Asia with an average of 0.32. For the whole panel, the mean LD (R2) decreased as genome distance increased, indicating that the probability of LD was low across widely separated SNP pairs with a rate of 1 Mbp (Figure S1B).
The PCA revealed sparse grouping according to geographic and institutional origin. The first axis explained 9.94% of the total genotypic variance, the second axis 3.12%. The first axis differentiated between wheat from Central Asia and European lines, which are located on opposite sides of the axis. The second axis separated between the North American and the Western Asia wheat, positioning them on opposite sides of the axis (Figure 1A).
In addition, Western Asian lines were separated across the second axis, with Iranian lines predominantly on one side and those near European lines, principally from Turkey, on the other side. In agreement with the PCA, AMOVA results indicated that the variance within subpopulations explained most of the genetic variation (94%; p ˂ 0001), whereas variation among subpopulations was only 6% (p ˂ 0001) (Table 2). AMOVA results within groups showed that the highest proportion of genetic variance (26%) was observed within lines from the CGIAR (IWWIP, CIMMYT, and ICARDA) followed by the subpopulation of North America and Western Asia with 23% and 22%, respectively. However, the Western Asian subpopulation (mostly from Turkey and Iran) showed the largest genetic variation per entry as indicated by the mean squares (Table 2).
To determine the genetic diversity trend across geographic origin, the genetic distance among the six subpopulations was compared. The Nei’s genetic distance ranged from 0.012 to 0.1 with an overall mean distance of 0.05. The CGIAR subpopulation showed the lowest genetic distance with all other subpopulations with an average of 0.037 (Table 3).

2.2. Characterization Based on Gene Diversity for Ppd and Vrn

Regarding the major flowering time genes, 19% of the genotypes carried Ppd-B1a and 80% Ppd-B1b, while 61% carried Ppd-D1a and 31% the sensitive Ppd-D1b allele (Figure 2B). The Ppd-B1b/Ppd-D1a allelic combination was the most common in the WWAGI with 59% of the genotypes, followed by 22% for Ppd-B1b/Ppd-D1b. For the main vernalization genes, 89%, 76%, and 87% of the genotypes carried the winter type allele for Vrn-A1, Vrn-B1 and Vrn-D1, respectively, while only 8% carried the spring allele for Vrn-A1, 21% for Vrn-B1, and 9% for Vrn-D1 (Figure 2C). As expected, the triple winter Vrn alleles were the major component of the WWAGI with 66%, while all the other allelic combinations represented less than 16%. In addition, the AMOVA results suggest that the different allelic combinations of Vrn and Ppd had a limited effect on the population structure (Table 2).

2.3. Population Structure

The population structure of the WWAGI was examined using the sNMF function, and the results showed that the cross-entropy criteria did not reveal a minimum value or distinct plateau and continuously dropped as k increased. However, significant results were observed at K2 and K3. At each K level, genotypes were allocated to different subpopulations, and genotypes with an ancestry coefficient less than 70% were classified as admixed and mapped to principal component analysis (PCA) clusters. At the K = 2 level, 20% were assigned to subpopulation 1 and 26% to subpopulation 2, while the remaining 54% were admixture. Most of subpopulation 1 comprised lines from Europe and Western Asia (64% and 22%). Subpopulation 2 contained lines from three groups, Western Asia, CGIAR, and North America, representing 91%. At K = 3, 64% of the lines were admixed, while the 36% remaining were divided between three subpopulations with 17%, 10%, and 8% for subpopulations 1, 2, and 3, respectively. Subpopulation 1 was composed mainly of the European lines, while subpopulation 2 was largely made of North American lines. For subpopopulation3, 57% were representing Western Asia and 33% represented lines from CGIAR (Figure 1C).
In the geographical region, Vrn-1 and Ppd-1 genes are lined up against the phylogenetic trees (in three different layers). The phylogenetic can be interpreted to split the dataset into four clades (Figure 2A). The subgroup I was mainly composed of genotypes coming from North America with few lines from Europe and CGIAR, and the second subgroup II principally combines Lines from Europe with few from Central, Western Asia, and CGIAR, while subgroup III comprises genotypes coming from North America, CGIAR, and Europe. The fourth subgroup is primarily composed of lines from Western Asia and CGIAR. The first and the second subgroup (I and II) englobed essentially triple winter types for the three copies of the Vrn gene (Figure 2A). For the Ppd gene, the first subgroup is mainly composed of the double sensitive alleles at Ppd-B1b and Ppd-D1b, and the sensitive Ppd-B1b and insensitive Ppd-D1a. At the same time, the second subgroup is predominantly made of the sensitive/insensitive Ppd-B1b/Ppd-D1a. For the subgroup III, the Vrn gene is represented by the triple winter type and double recessive winter type, while for the Ppd gene, the subgroup III englobed the sensitive/insensitive Ppd-B1b-/Ppd-D1a mainly, followed by the insensitive/sensitive Ppd-B1a/Ppd-D1b. The last subgroup IV was represented mainly by the double winter type and double recessive winter type for the Vrn gene, and by the sensitive/insensitive Ppd-B1b/Ppd-D1a for the Ppd gene.

2.4. Response to Photoperiod and Vernalization

The results of the phenology experiment (Table 4) showed that the genotype effect was significant in all treatments for the GDD to heading date. GDD ranged from 1013.9 to 1629.42 with an average of 1273.13 for full pre-vernalization and a long photoperiod (V+P+), whereas for no vernalization and a long photoperiod (V−P+), it ranged from 896.85 to 2241.87 with an average of 1417.13. The highest average was found for V+P− with 1705.05, which ranged from 1270 to 2150. Across all treatments, the GDD showed high heritabilities (Table 5).
The highest heritability was detected for V+P+ (h = 0.92) followed by V−P+ with 0.91 and V+P− by 0.81. On the other hand, there was a strong significant difference for GDD across treatments and the interaction of genotype × treatment (Table 6).
For the response to vernalization and photoperiod treatments of the WWAGI, the experiment revealed that among all lines, only 35% (85 lines) were able to flower under no vernalization, whereas 70% of the lines (179) under short photoperiod were able to flower. The results showed that the presence of one spring allele is sufficient to trigger flowering for 78% of the lines carrying one (Figure 3A). However, just 36% for Vrn-A1, 27% for Vrn-B1, and 36% for Vrn-D1 with winter alleles were able to flower. At the same time, the percentage of flowering was almost similar for the lines carrying one or two spring alleles at the Vrn loci with 83% and 85, respectively, while it was only 15% for the triple winter type (Figure 3B). Most of the winter wheat lines that flowered without vernalization still showed a quantitative response to vernalization, and this response was associated with the Vrn allele they carried. Thus, for the average extra days to flowering under no vernalization, the lines carrying the three winter alleles showed the highest value with 33 extra days as compared to the vernalized treatment. The lines carrying winter alleles at Vrn-B1 showed a delay in flowering under no vernalization of 35 days, while the genotypes with spring alleles revealed a lower number of extra days to flowering with only 9 extra days for Vrn-A1 and 22 days for Vrn-B1 as compared to growing under 49 days of vernalization. At the same time, Vrn-D1 did not show a clear effect as compared to its spring version (Figure 3C). This result indicated that the Vrn-A1 gene contributed the most to flowering time under non-vernalization. On the other hand, Vrn-D1 showed the weakest effect on flowering, and therefore, the presence of its spring or winter allele did not delay the flowering time significantly. Furthermore, considering any winter allele, there was a noticeable increase in the average extra days to flowering when comparing genotypes with a single winter allele to those with double or triple winter alleles. Genotypes with a single winter allele exhibit a delay of 14 days, while it was 22 and 49 days for double and triple winter allele genotypes, respectively (Figure 3D). The number of winter alleles increased the extra days to flowering by 71% from single to triple. Regarding the photoperiod treatment (Figure 3E), the percentage of lines that achieved flowering under short photoperiod was 100% for lines with both insensitive Ppd alleles, while the percentage was reduced by only 9% for those carrying the sensitive Ppd-B1b. However, for those carrying the sensitive Ppd-D1b allele, the percentage of flowering was decreased to 56%, while it was only 11% for the lines carrying both sensitive Ppd alleles (Ppd-B1b/Ppd-D1b). Most lines that flowered under a short photoperiod also showed a quantitative response to a long photoperiod, which was related to the Ppd allele they carried (Figure 3E). In addition, the highest value of extra days under short photoperiod was detected for the Ppd-B1b/Ppd-D1b combination with 21 days alongside with Ppd-B1a/Ppd-D1b with 22 days, while Ppd-B1b/Ppd-D1a and Ppd-B1a/Ppd-D1a combinations showed a lower number of days with 15 and 17 days, respectively (Figure 3F). This suggested that Ppd-D1 is the strongest regulator of flowering time under a short photoperiod, while no significant effect of Ppd-B1 on flowering was detected.
The phenology experiment revealed that out of the 249 lines used, only 31 lines needed both vernalization and a long photoperiod to flower, while 102 lines required only vernalization, 20 lines needed a long photoperiod, and 65 genotypes required neither (Figure 4).
In terms of region-specific distribution (Figure 5), the CGIAR lines were mostly composed of two major groups: lines that need vernalization but not a long photoperiod (49%), and those that did not require vernalization but required a long photoperiod (44%). This was similar for the North America lines with 45% and 40%, respectively. In the European group, most lines tested required vernalization but not a long photoperiod (76% of the lines), while 13% required a long photoperiod without vernalization. For Central Asia, three groups with identical percentages (33% each) were present. However, almost 88% of the lines coming from the Western Asia region either did not require vernalization and a long photoperiod, or they needed vernalization but not a long photoperiod (Figure 5).
The quantitative response to the treatments allowed the calculation of the relative response to vernalization (RRV) and relative response to photoperiod (RRP) indices, quantifying the sensitivity of the WWAGI set to these. Across the WWAGI, and in absolute terms, the RRV ranged from 0 to 0.56 with an average of 0.2, while RRP was higher with an overall average of 0.25, but less variable. In comparison to Eps, both RRP and RRV showed a significant moderate correlation with R2 of 0.33 for RRP and of 0.12 for RRV (Figure 6A,B). However, no significant correlation was detected between RRP and RRV themselves (Figure 6C).
Furthermore, the results revealed that Eps was significantly higher for the lines that did not require vernalization to flower than for those who did. No difference could be found in Eps between the lines requiring or not a long photoperiod to flower. The genotypes that required vernalization to flower were mainly originated from Europe, North America, and CGIAR, whereas those from Central Asia and Western Asia (Iranian group) did not, in general, require vernalization (Figure 7E). Furthermore, 77% were also photoperiod-insensitive, while the 23% remaining were photoperiod-sensitive. In addition, most North America lines and 6% from the CGIAR group exhibited photoperiod sensitivity, whereas the remaining were found to be insensitive to the photoperiod. This suggests that there may be a potential relationship between geographic origin and photoperiod sensitivity in these lines (Figure 7F). Regarding the vernalization and photoperiod experiment (Figure 7G), PCA revealed two groups clustered with their respective vernalization and photoperiod requirements. Furthermore, no clear pattern was detected in the relationship between RRP, RRV, Eps, and population structure. However, a significant differentiation was apparent in the case of strict vernalization and photosensitivity. Specifically, most of the European and North American lines exhibited a strict vernalization requirement, while lines originating from Asia and CGIAR did not require vernalization. Regarding strict photosensitivity, only a subset of lines from North America and CGIAR required a long photoperiod, indicating that the population structure may be influenced by strict vernalization and photoperiod requirements.

2.5. Genome-Wide Association Study

GWAS for the traits studied identified 17 significant markers (6 in genomes A and B and 5 in genome D; −log10(p) > 3) located on chromosomes 1A, 2D, 5A, 5B, 5D, 7B, and 7D. The percentage of phenotypic variation explained by significant markers ranged from 9% to 49%. The highest percentage was observed for RRV with an average of 43%, while the lowest value was detected for RRP. The percentage of phenotypic variation explained by functional genes ranged from 12% to 44%, while from 9% to 45%, it was explained by the newly identified markers. Four significant markers were detected for strict photosensitive, distributed across three chromosomes and in charge of 25–41% of the phenotypic variation (Table 7).
The Ppd-D1 in chromosome 2D was the most significant association (LOD = 11.8) and responsible for the highest proportion (41%) of the phenotypic variation for strict photosensitivity. The lowest number of MTAs were identified for Eps with only 1 marker Vrn-B1 in chromosome 5B (LOD = 3.7). For strict vernalization, five significant markers were identified, distributed in chromosomes 1A, 5A, and 5B. The mean of the phenotypic variation explained by these markers was 29%. Allele Vrn-B1 explained the highest phenotypic variation with 36% and was present in 54 lines; among them, 78% were characterized by strict vernalization, and all of them carry the winter type allele of Vrn-B1. Out of the seven MTAs, the remaining five were identified for RRV and two for RRP. Associated SNPs identified for each trait are listed with their respective LOD in Table 7. Only one marker was associated with two different traits: Vrn-B1, with strict vernalization and Eps. The marker AX-94796479 was present in only 17 lines, of which 11 showed no relative response to vernalization, and the remaining showed a response below 0.13 (Figure 8).

3. Discussion

3.1. Genetic Diversity and Population Structure

Wheat breeding is a continuous process and exploitation of genetic diversity is crucial. For attaining specific breeding goals not only is the discovery of genes important, but also the contribution of these genes needs to be analyzed and dissected, and their association with specific traits of interest needs to be assessed and quantified. This process of gene discovery will ultimately contribute to a deeper knowledge of the wheat genome and its function. Within this framework, we investigated the genetic diversity of a worldwide winter wheat association genetics panel, originated from seven different regions and 32 countries. This population showed high GD (0.47) and PIC (0.33), suggesting a higher levels of polymorphism as compared to previous studies [31,32,33]. Moreover, the B genome showed the highest GD and PIC followed by the A genome and finally the D genome, similarly to previous studies [33,34].
The results from the two clustering methods (PCA, and NJ tree) demonstrated the diversity of the 249 WWAGI lines grouped into three subpopulations (Europe, North America, and Western Asia) and the presence of additional subgroups characterized by common related origins (CGIAR). This can be attributed to the CGIAR’s utilization of the entire spectrum of genetic diversity available among modern varieties and the increased adoption of CGIAR advanced lines through worldwide distribution of international nurseries. Similarly, Ref. [35] studied a panel of 1000 wheat accessions from 91 countries and found that the population was clustered based on geographical origin in three main groups, European, Asian, and American. Geographic origin has previously been proven to be a major driver of population clustering. For instance, Ref. [36] showed in a panel of 3200 wheat accessions that the population was structured based on geographic origin and growth habit. Also, Ref. [37] revealed that using 635 wheat accessions that the geographic origin and growth habit were the primary basis for clustering. AMOVA analysis based on geographical origin revealed that more than 6% of the overall genetic variations were identified among the origin groups. These observations are comparable to those reported by Ref. [38] with a collection of 290 winter wheat varieties and by Ref. [39] with a panel of 753 U.S. wheat varieties including 517 winter wheat varieties, and are lower than those reported by Ref. [40], where 11% of variations were captured among countries. The modest diversity among the modern winter wheat panel (6% based on regions) displayed in Table 2 was most likely because of the large and regular sharing of parents among the wheat community. Nevertheless, these findings indicate that there is genetic diversity available for use within wheat breeding programs. The clustering patterns shown by PCA and NJ tree offer strong evidence that the panel contains a pool of unused genetic potential. By confirming the presence of distinct subpopulations, such as those associated with CGIAR lines, consequently, it is possible that these lines have valuable traits that can be used to enhance the wheat breeding program effort.
A large number of studies have used Vrn and Ppd alleles to explain genetic diversity and population structure [12,41,42]. In our study, we explored population structure taking into consideration the gene-specific molecular markers Vrn and Ppd. For the vernalization genes, Vrn-A1a was almost absent in the WWAGI with only 21 lines, while vrn-A1 was present in 222 lines. Lines from CGIAR, Europe, and Western Asia lines have been characterized with photoperiod sensitivity at Ppd-B1 and insensitive Ppd-D1 as the most common allelic combination for Ppd, while the most common Central Asia was the double sensitive and the insensitive/sensitive Ppd-B1a/ppd-D1b for North America. However, the Central Asian genotypes was notably the smallest among the studied groups, which may raise a statistical limitation or potential bias due to the limited sample size. These outcomes are comparable to those discussed by [43] with 257 wheat cultivars from CIMMYT where 91% carried Ppd-D1a. However, Ref. [44] described that out of 245 European cultivars, 224 carried the sensitive allele Ppd-D1b and only 8.5% were carried the photoperiod-insensitive allele Ppd-D1a, found mainly in southern Europe, while Ref. [45] showed that 58% out of 521 wheat cultivars from Europe carried the insensitive allele Ppd-D1a.
For Vrn, and as expected, the most common allelic combination was the winter type for all regions except Central Asian lines that were double winter type for Vrn-A1 and Vrn-B1 and mostly carried the spring type for Vrn-D1.

3.2. Phenology Variability and Marker Trait Association

Flowering time in wheat has previously been linked to variation in the Vrn and Ppd genes such as insertions, deletions, and mutations [46]. In this study, the panel was characterized in a growth chamber with controlled light and temperature, simulating different and extreme vernalization and photoperiod conditions, to determine the phenology response. In general terms, it has been established that spring wheat growth habit does not require any cold treatment to transition to flowering [47]. However, for winter wheat growth habit, as vernalization is a quantitative response, increasing the duration of vernalization until requirements are met result in accelerated flowering [48]. Under no vernalization, we observed that the percentage of flowering was more than 78% for lines carrying spring alleles at the Vrn loci, while it was less than 36% for those carrying winter Vrn alleles. All the lines that flowered under no vernalization still revealed a response to low-temperature exposure, and this response was linked to the allelic combination at Vrn loci. Specifically, lines carrying triple winter Vrn alleles had the lowest percentage of flowering (less than 15%) as compared to double and single Vrn type (more than 85%), reflecting the cumulative influence of Vrn genes on flowering. This highlights that Vrn genes impact significantly flowering as described by large body of research [48,49,50]. As expected, this allelic variation significantly impacted the extra days to flower, for Vrn-A1 and Vrn-B1 we found that the genotypes carrying the spring alleles showed the lowest average extra days to flowering, as compared to their response under long vernalization, from less than 22 days to more than 33 days for winter lines. Furthermore, we found that the difference in extra days to flowering between winter and spring alleles is higher at Vrn-A1 with 24 days as compared to 13 days for Vrn-B1. This is in agreement with several research suggesting that Vrn-A1 is the strongest driver of vernalization requirement [14,51,52], while no difference was detected between lines carrying the Vrn-D1 spring allele and those carrying the winter allele. This finding is in agreement with results published by [53]. However, lines carrying the triple winter combination can flower under no vernalization, derived mainly from the CGIAR and Western Asia regions. This is probably due to another source of variation in vernalization response via copy number variation (CNV) as reported by [48,50] and/or additional Vrn genes, e.g., Vrn-3. CNV for Vrn-A1 was reported as the highest among the Vrn loci [50,54], which can be from one to four copies and therefore may exhibit multiple effects on flowering, such that the higher the number of copies, the longer the vernalization duration required [55]. This suggests that the use of genetic markers associated with a Vrn allele in breeding programs may be enough to select for gross phenological adaptation, but for fine-tuning it, markers for CNV, either using associated haplotypes or through RT-PCR, are necessary.
Wheat photoperiodism is primarily regulated by the Ppd-D1 locus on chromosome 2D, which strongly controls flowering time by separating lines into photoperiod-sensitive and -insensitive. As a result, it influences the proportion of lines that flower under a short photoperiod and therefore the extra days required for flowering [56,57]. Analogously, lines carrying the insensitive allele of Ppd-D1a showed the highest percentage of flowering under short photoperiod conditions with more than 91%, as compared to less than 56% for lines carrying the sensitive allele of Ppd-D1. Moreover, the average extra days under short photoperiod was lower for the Ppd-D1-insensitive allele (less than 17 days), while it was more than 3 weeks for the sensitive Ppd-D1 allele. In addition, no difference was detected in the extra days to flowering between lines carrying the sensitive and the insensitive Ppd-B1 allele but still impacting the absolute sensitivity to photoperiod. We found that lines carrying the sensitive allele at Ppd-D1 exhibit a delay in flowering by 5.5 days. Similarly, Ref. [58] reported in a panel of 1100 winter wheat cultivars that the strongest association with phenology including heading date was observed for Ppd-D1, thus delaying flowering by 4.8 days [58]. It was demonstrated in a large panel of winter wheat that Ppd-D1 was the primary switch to flowering and that the copy number at Ppd-B1 was the second most important source of variation, although [59] reported that Ppd-B1 could be as strong as Ppd-D1, and both of them significantly accelerate spike development under a short photoperiod. This is in agreement with previous observations [50,60,61]. Our results on the photoperiod response demonstrated that the sensitive allelic combination of both Ppd-B1 and Ppd-D1 plays a significant role in blocking flowering under a short photoperiod.
The geographic pattern distribution of the vernalization and photoperiod response experiment in our WWAGI revealed that most of the European lines (76%) required long vernalization to flower but could still do it under short days. Most probably, this is because most of them carry the insensitive Ppd-D1 allele. Refs. [56,62] reported that southern and eastern European wheats were usually photoperiod-insensitive, while wheat from the northern latitudes of Europe were photoperiod-sensitive. Analogously, the majority of our European lines (90%) originated from eastern and southern Europe. Regarding North American lines, the experiment demonstrated that the vast majority (85%) need long vernalization to flower, and 40% of these lines require long photoperiods as well. These results concur with those reported by [63] in Canadian winter wheat material with a diversity panel of 203 winter wheat genotypes and [64] using 299 winter wheat genotypes from USA. Almost all the CGIAR lines (93%) were photoperiod-insensitive with 49% of them requiring long vernalization, suggesting that all these lines were adapted to lower latitudes. This is expected given that the vast majority of CGIAR wheat have insensitive photoperiod genes in their ancestors [49]. The Central Asia lines showed the highest proportion of photoperiod-sensitive lines (66%), and this is probably due to the high latitude at which they are grown. Nevertheless, it is important to mention that these genotypes represented the smallest group among the studied groups, which may introduce potential statistical limitations. The distribution across the Western Asian lines was similar to CGIAR lines with slight differences in those that do not require vernalization and long photoperiod. The proportion of the Western Asian lines that did not require vernalization and long photoperiod was 8% higher than the CGIAR lines, while they were less with 14% for genotypes that needed vernalization but not a photoperiod. Strict vernalization and photoperiod requirements exhibited a major impact on population structure since the European and North American Lines were clustered together as strict vernalization, while Asian and CGIAR lines were not. Similarly for strict photosensitivity, one cluster from North America and CGIAR were sensitive to a long photoperiod, confirming that selection based on vernalization and photoperiod requirements play an important role in population structure.
To identify unique and novel MTAs for specific traits and under specific conditions, genome-wide association studies have been used extensively in wheat [23,25,65,66,67,68]. Previous studies have revealed that genes involved in flowering are mostly found on chromosomes 1A, 2B, 3A, 3B, 5A, 6A, 6B, 7A, 7B, and 7D [69,70,71]. In the current study, a total of 30 significant markers associated with five phenological traits were identified. The high significant functional marker Ppd-D1 was found for the strict photosensitive which is consistent with previous genetic research that identified the Ppd-D1 gene is the most efficient photoperiod response gene, followed by Ppd-B1 and Ppd-A1 [23,60]. For the strict vernalization, we identified two significant MTAs of Vrn-B1 located on chromosome 5B; these identified genes serve as the primary regulators of plant growth in the vernalization pathway [13,72], and are overlapping with several studies previously reported [12,16,41,45,73]. The Vrn-B1 gene involved considerably in the regulation of flowering, and it has been found to be impacted by copy number variation [48,50]. During the investigation, we found that lines with the Vrn-B1a allele flowered 12 days earlier than those with Vrn-B1. This finding is in good agreement with the results published by [23], which revealed the same earliness of flowering for cultivars with allele Vrn-B1a but with just 3 days. We have observed that the LOD score of Ppd was higher than Vrn, and our result demonstrated that Ppd-D1 gene had the highest effect on flowering followed by Vrn-B1. The authors of [23] revealed that Ppd-D1 had the most important effect on heading and flowering then Vrn-B1 in Chinese wheat. This is due to the fact that Vrn genes are positioned in a more dynamic region of the wheat genome as compared to Ppd genes; as a result, the Vrn genes have more allelic variation and genetic diversity than Ppd genes. Furthermore, it is demonstrated that Eps genes are located near the centromere of 5B chromosome [74]. Until now, Eps-Am1 in Triticum monococcum on chromosome 1A is the best identified Eps locus in wheat [19,75]. Our results indicated that gene Vrn-B1 on chromosome 5B which is related to vernalization, has a significant impact on Eps and therefore on the flowering time in wheat.

4. Materials and Methods

4.1. Plant Material

For this study, a panel of 249 advanced lines and varieties from the Winter Wheat Association Genetics Initiative (WWAGI) developed by the International Maize and Wheat Improvement Center (CIMMYT) was used (Table S1). The lines originated from five different regions—Central and Western Asia, Europe, and North America—or originated from the CGIAR International Winter Wheat Improvement Program (IWWYP) a joint venture of the Turkish Government, the International Center for Maize and Wheat improvement (CIMMYT), and the International Center for Agricultural Research in the Dry Areas (ICARDA). The genetic material used are representative of the International Winter Wheat Improvement Program germplasm dissemination work (International Nurseries) and its international collaboration. Constituting a set of relevant breeding material representative of a diverse set of winter wheat breeding programs. Only a smaller number of lines derived from other regions. The largest number of lines (65) originated from the CGIAR (CGR), followed by United States (61) and Iran (27), Turkey (22), Russia (10), Ukraine (9), Bulgaria (6), and Moldova (5) (Figure 9A,B).

4.2. Genotyping

The WWAGI was genotyped with the 35K Axiom® Wheat Breeder’s Array (Affymetrix UK Ltd., High Wycombe, UK) as per manufacturer’s guidelines [76]. The array comprised a total of 35,143 SNP markers. Genomic DNA was isolated from plants at the growth stage. Modified CTAB method was applied for the DNA isolation [77]. Vrn-1 alleles (Vrn-A1, Vrn-B1, and Vrn-B3) were characterized to determine the growth habit of each genotype. The gene-specific markers described by Refs. [17,78] were used to identify variation in the promoter and intron-1 region of the Vrn-1 genes. In addition, the presence of a SNP, located in Exon 4 of Vrn-A1 was tested [52]. Ppd-1 alleles which confer the sensitivity to day lengths were genotyped using gene-specific markers related to the junction between intact Ppd-B1 copies in the ‘Sonora64’ allele and the 2089 bp deletion in the Ppd-D1 gene [79].
Polymerase chain reaction (PCR) amplifications for STS markers were performed in a total volume of 10 μL containing a final concentration of 1× Buffer with Green Dye (Promega Corp. Madison, Wisconsin, USA), 200 μM deoxynucleotide triphosphates (dNTPs), 1.2 mM magnesium chloride (MgCl2), 0.25 μM of each primer, 1U of DNA polymerase (GoTaq®Flexi, Promega Corp., Cat. # M8295), and 50 ng of DNA template. PCR conditions were performed using the following temperature profile: 94 °C for 2 min followed by 30 cycles of 94 °C for 1 min, 54–60 °C for 2 min (dependent on the primer), and 72 °C for 2 min. PCR amplified products were separated by electrophoresis on 1.2% agarose gels using 1× TAE buffer, visualized under UV light. The SNP polymorphisms were scored using KASP reagents. Reactions contained 2.5 mL of water, 2.5 mL of 2× KASPar reaction mix, 0.07 mL of assay mix, and 50 ng of dried DNA with a PCR profile of 94 °C for 15 min activation time, followed by 20 cycles of 94 °C for 10 s, 57 °C for 5 s, and 72 °C for 10 s, followed by 18 cycles of 94 °C for 10 s, 57 °C for 20 s, and 72 °C for 40 s. Fluorescence was read as an end point reading at 25 °C.
After removing markers with minor allele frequencies (MAF) of <5% and markers with more than 10% missing data, a refined set of 11,476 Axiom SNPs markers remained. The distribution of the filtered SNPs across the seven wheat chromosomes (Figure S1A) was drawn using the Cmplot package for R 4.1.1 statistical software (R Core Team, 2019, R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org/) [80].

4.3. Population Structure

The Remington approach was used to conduct a linkage disequilibrium (LD) study [81]. The population structure was assessed using Neighbor joining (NJ) tree and principal component analysis (PCA). NJ analysis was conducted using the “BiocManager” and “ggtree” R packages [82], while LD and PCA were performed using TASSEL version 5.0 [83]. To further uncover population structure of WWAGI, a sparse “Non-negative Matrix Factorization” (sNMF) algorithm in LEA package was performed in R [84]. The method consists of estimation of admixture coefficient for each individual separately without taking into consideration a previous population membership supposition. The cross-entropy criterion from the sNMF function was used to check various potential population sets ranging from K = 2 to K = 10. Admixed lines were set to have less than 70% ancestral membership to any subpopulation.

4.4. Diversity Analysis

Major allele frequency (MAF), genetic diversity (GD), and polymorphism information content (PIC) for all chromosomes and individual regions and origin of winter wheat lines were calculated using PowerMarker version 3.25 software [85]. The analysis of molecular variance (AMOVA) and Nei’s genetic distance was performed using the same software to estimate genetic diversity between and within the predefined subpopulations. The optimal distances were assessed by a permutation test of 1000 [85].

4.5. Growth Chamber Experiment and Phenology Measurements

Three independent experiments were carried out in a climate-controlled growth chamber at 25 °C ± 1 °C temperature with a combination of two photoperiod (long photoperiod ‘P+’ and short photoperiod ‘P−’) and two vernalization durations (‘V+’ for full pre-vernalization and ‘V−’ for no vernalization). The details of the temperature and light treatments are provided in Table 4.
The experiments were set up growing one plant per genotype in a 1 L pot filled with peat moss and arranged in an augmented design with two repeated check cultivars NEKOTA and SERI, repeated 15 times per experiment. For each experiment both, days from germination to heading and growing degree days (GDD) for the same period were computed adopting the method of [72,86]. Extra days to flowering were assessed by deducting days required to flower under non-vernalization or under short photoperiod from days needed to flowering under vernalization or long photoperiod. Three phenological traits were calculated: earliness per se (Eps) was obtained using the rate of development (R) under long photoperiod and full vernalization, and relative response to vernalization (RRV) and relative response to photoperiod (RRP) were calculated according to [21,87] as follows:
R = 1 F                   R R V = 1 − R N V R F V                   R R P = 1 − R S P R L P
F = growing degree days from germination to heading date. RNV = rate of development under long photoperiod and no vernalization. RFV = rate of development under long photoperiod and full vernalization. RSP = rate of development under short photoperiod and full vernalization. RLP = rate of development under long photoperiod and full vernalization.

4.6. Statistical Analysis

For each experiment, a linear mixed model with spatial adjustment was built and used to determine adjusted means. The R package “statgenSTA” [88] was used to fit the model, which comprised a fixed effect for genotype and block, and an additional spatial component derived using Spatial Analysis of Field Trials with Splines (SpATS), which utilizes two-dimensional smoothing with P-splines as reported by [89]. Heritability (h2) was computed as a function of the effective dimension (ED) to the genetic component (EDg) h2 = EDg/(ng−l), where ng is the number of genotypes [90]. BLUEs (best linear unbiased estimators) were estimated with genotypes considered as fixed effect, heritability was treated as random [91].

4.7. Genome-Wide Association Study (GWAS)

The filtered SNPs, gene-specific marker results, and the BLUEs for each quantitative phenological trait were used for the GWAS analysis using GAPIT V3 (Genomic Association and Prediction Integrated Tool) in R [92]. Three models were performed: GLM [93], MLM [94], and CMLM [95] (https://github.com/jiabowang/GAPIT3 (accessed on 8 September 2023)). In addition to PCA, a relative kinship analysis was carried out to identify the genetic relationship between genotypes using the VanRadden method implemented in GAPIT [96,97]. Generalized linear mixed-model association test (GMMAT) was used for testing the vernalization requirement and photoperiod sensitivity (strict vernalization and photoperiod) [98]. Marker trait associations (MTA) were declared if their −log10(p) value was higher than 3.0. The significance threshold for MTA was elaborated based on Bonferroni correction of p ≤ 0.05/n, where n is the number of SNP markers used. Due to the conservative nature of the Bonferroni adjustment and to avoid false negatives, markers with an LOD > 3 are also highlighted and reported. A Q–Q plot with a predicted vs. observed log10(p) value was used to determine the quality of the fitted GWAS model. Manhattan plots was drawn using the CMplot R package (https://cran.r-project.org/web/packages/CMplot/index.html) to visualize the MTA.

5. Conclusions

Wheat is a crucial crop for global food security. In the present study, the WWAGI panel has been studied. The characterization and genetic dissection of the vernalization and photoperiod diversity of this panel added more value to our understanding of linkages, diversity, and population structure, while also making it easier to identify the best sources of genetic variability. As well as their relationship with phenology traits, quantitative and qualitative vernalization and photoperiod requirement in controlled condition can provide a valuable alternative to fields trials as it appears to avoid numerous constraints that may affect response accuracy. Additionally, this information will support the breeding program to improve the adaptation of the new winter wheat varieties. Genome-wide association for phenology in the WWAGI showed that a number of regulatory vernalization and photoperiod genes, such as Ppd-D1, Vrn-A1, and Vrn-B1, have a significant impact on the genetic architecture of flowering time. However, the distribution of significant MTAs led to the conclusion that there are many additional important genetic loci involved.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants12234053/s1, Figure S1: (A) Distribution in the winter wheat genome based on the physical map of SNPs markers used in the winter wheat association genetics initiative (WWAGI) of 249 advanced lines and varieties from all over the world. Color legend represents marker density. (B) Linkage disequilibrium (LD) decay expressed genetic map distance and dashed line indicates the derived threshold for LD due to linkage. Table S1: Genotypes information used in the winter wheat association genetics initiative (WWAGI) of 249 advanced lines and varieties. Table S2: AMOVA of the winter wheat association genetics initiative (WWAGI) based on Vrn and Ppd Allelic combination.

Author Contributions

A.E.B., M.S.L. and M.S.G. conceived the study; M.S.L. assembled the panel; A.E.B. and M.S.G. did the phenotyping; S.D. contribute the genotyping; A.E.B. and M.S.G. did computational analysis and visualization; A.E.B. performed the research work and prepared the first draft of the manuscript; M.I., M.S.L., S.D. and M.S.G. contributed to the final draft of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by a European Union project entitled “Addressing the Challenges of Climate Change for Sustainable Food Security in Turkey, Iran and Morocco, through the creation and dissemination of an international database to promote the use of wheat genetic resources and increase genetic gains.” CFP 2014/2015-W3B-PR-18-Turkey.

Data Availability Statement

The data presented in this study are available in the article. The raw MS files are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. FAO. Crop Prospects and Food Situation; FAO: Rome, Italy, 2019; ISBN 9789251358580. [Google Scholar]
  2. Mohammadi, R. Breeding for Increased Drought Tolerance in Wheat: A Review. Crop Pasture Sci. 2018, 69, 223–241. [Google Scholar] [CrossRef]
  3. Semenov, M.A.; Stratonovitch, P.; Alghabari, F.; Gooding, M.J. Adapting Wheat in Europe for Climate Change. J. Cereal Sci. 2014, 59, 245–256. [Google Scholar] [CrossRef]
  4. Andrés, F.; Coupland, G. The Genetic Basis of Flowering Responses to Seasonal Cues. Nat. Rev. Genet. 2012, 13, 627–639. [Google Scholar] [CrossRef]
  5. Gegas, V.C.; Nazari, A.; Griffiths, S.; Simmonds, J.; Fish, L.; Orford, S.; Sayers, L.; Doonan, J.H.; Snape, J.W. A Genetic Framework for Grain Size and Shape Variation in Wheat. Plant Cell 2010, 22, 1046–1056. [Google Scholar] [CrossRef]
  6. Khadka, K.; Torkamaneh, D.; Kaviani, M.; Belzile, F.; Raizada, M.N.; Navabi, A. Population Structure of Nepali Spring Wheat (Triticum aestivum L.) Germplasm. BMC Plant Biol. 2020, 20, 530. [Google Scholar] [CrossRef]
  7. Mazzucotelli, E.; Sciara, G.; Mastrangelo, A.M.; Desiderio, F.; Xu, S.S.; Faris, J.; Hayden, M.J.; Tricker, P.J.; Ozkan, H.; Echenique, V.; et al. The Global Durum Wheat Panel (GDP): An International Platform to Identify and Exchange Beneficial Alleles. Front. Plant Sci. 2020, 11, 569905. [Google Scholar] [CrossRef]
  8. Pascual, L.; Ruiz, M.; López-Fernández, M.; Pérez-Penã, H.; Benavente, E.; Vázquez, J.F.; Sansaloni, C.; Giraldo, P. Genomic Analysis of Spanish Wheat Landraces Reveals Their Variability and Potential for Breeding. BMC Genom. 2020, 21, 122. [Google Scholar] [CrossRef]
  9. Rufo, R.; Alvaro, F.; Royo, C.; Soriano, J.M. From Landraces to Improved Cultivars: Assessment of Genetic Diversity and Population Structure of Mediterranean Wheat Using SNP Markers. PLoS ONE 2019, 14, e0219867. [Google Scholar] [CrossRef]
  10. Cavanagh, C.R.; Chao, S.; Wang, S.; Huang, B.E.; Stephen, S.; Kiani, S.; Forrest, K.; Saintenac, C.; Brown-Guedira, G.L.; Akhunova, A.; et al. Genome-Wide Comparative Diversity Uncovers Multiple Targets of Selection for Improvement in Hexaploid Wheat Landraces and Cultivars. Proc. Natl. Acad. Sci. USA 2013, 110, 8057–8062. [Google Scholar] [CrossRef]
  11. Ali, M.; Danting, S.; Wang, J.; Sadiq, H.; Rasheed, A.; He, Z.; Li, H. Genetic Diversity and Selection Signatures in Synthetic-Derived Wheats and Modern Spring Wheat. Front. Plant Sci. 2022, 13, 877496. [Google Scholar] [CrossRef]
  12. Muterko, A.; Kalendar, R.; Salina, E. Allelic Variation at the VERNALIZATION-A1, VRN-B1, VRN-B3, and PHOTOPERIOD-A1 Genes in Cultivars of Triticum Durum Desf. Planta 2016, 244, 1253–1263. [Google Scholar] [CrossRef]
  13. Yan, L.; Loukoianov, A.; Tranquilli, G.; Helguera, M.; Fahima, T.; Dubcovsky, J. Positional Cloning of the Wheat Vernalization Gene VRN1. Proc. Natl. Acad. Sci. USA 2003, 100, 6263–6268. [Google Scholar] [CrossRef]
  14. Yan, L.; Loukoianov, A.; Blechl, A.; Tranquilli, G.; Ramakrishna, W.; SanMiguel, P.; Bennetzen, J.L.; Echenique, V.; Dubcovsky, J. The Wheat VRN2 Gene Is a Flowering Repressor Down-Regulated by Vernalization. Science 2004, 303, 1640–1644. [Google Scholar] [CrossRef]
  15. Yan, L.; Fu, D.; Li, C.; Blechl, A.; Tranquilli, G.; Bonafede, M.; Sanchez, A.; Valarik, M.; Yasuda, S.; Dubcovsky, J. The Wheat and Barley Vernalization Gene VRN3 Is an Orthologue of FT. Proc. Natl. Acad. Sci. USA 2006, 103, 19581–19586. [Google Scholar] [CrossRef]
  16. Kippes, N.; Debernardi, J.M.; Vasquez-Gross, H.A.; Akpinar, B.A.; Budak, H.; Kato, K.; Chao, S.; Akhunov, E.; Dubcovsky, J. Identification of the VERNALIZATION 4 Gene Reveals the Origin of Spring Growth Habit in Ancient Wheats from South Asia. Proc. Natl. Acad. Sci. USA 2015, 112, E5401–E5410. [Google Scholar] [CrossRef]
  17. Yan, L.; Helguera, M.; Kato, K.; Fukuyama, S.; Sherman, J.; Dubcovsky, J. Allelic Variation at the VRN-1 Promoter Region in Polyploid Wheat. Theor. Appl. Genet. 2004, 109, 1677–1686. [Google Scholar] [CrossRef]
  18. Fjellheim, S.; Boden, S.; Trevaskis, B. The Role of Seasonal Flowering Responses in Adaptation of Grasses to Temperate Climates. Front. Plant Sci. 2014, 5, 431. [Google Scholar] [CrossRef]
  19. Zikhali, M.; Wingen, L.U.; Griffiths, S. Delimitation of the Earliness per Se D1 (Eps-D1) Flowering Gene to a Subtelomeric Chromosomal Deletion in Bread Wheat (Triticum aestivum). J. Exp. Bot. 2016, 67, 287–299. [Google Scholar] [CrossRef]
  20. Nazim Ud Dowla, M.A.N.; Edwards, I.; O’Hara, G.; Islam, S.; Ma, W. Developing Wheat for Improved Yield and Adaptation Under a Changing Climate: Optimization of a Few Key Genes. Engineering 2018, 4, 514–522. [Google Scholar] [CrossRef]
  21. Herndl, M.; White, J.W.; Hunt, L.A.; Graeff, S.; Claupein, W. Field-Based Evaluation of Vernalization Requirement, Photoperiod Response and Earliness per Se in Bread Wheat (Triticum aestivum L.). Field Crop. Res. 2008, 105, 193–201. [Google Scholar] [CrossRef]
  22. Borrill, P.; Adamski, N.; Uauy, C. Genomics as the Key to Unlocking the Polyploid Potential of Wheat. New Phytol. 2015, 208, 1008–1022. [Google Scholar] [CrossRef]
  23. Zhang, X.; Chen, J.; Yan, Y.; Yan, X.; Shi, C.; Zhao, L.; Chen, F. Genome-Wide Association Study of Heading and Flowering Dates and Construction of Its Prediction Equation in Chinese Common Wheat. Theor. Appl. Genet. 2018, 131, 2271–2285. [Google Scholar] [CrossRef]
  24. Alemu, A.; El Baouchi, A.; El Hanafi, S.; Kehel, Z.; Eddakhir, K.; Tadesse, W. Genetic Analysis of Grain Protein Content and Dough Quality Traits in Elite Spring Bread Wheat (Triticum aestivum) Lines through Association Study. J. Cereal Sci. 2021, 100, 103214. [Google Scholar] [CrossRef]
  25. Shewabez, E.; Bekele, E.; Alemu, A.; Mugnai, L.; Tadesse, W. Genetic Characterization and Genome-Wide Association Mapping for Stem Rust Resistance in Spring Bread Wheat. BMC Genomic Data 2022, 23, 11. [Google Scholar] [CrossRef]
  26. Devate, N.B.; Kumar Parmeshwarappa, S.V.; Kumar Manjunath, K.; Chauhan, D.; Singh, S.; Bahadur Singh, J.; Kumar, M.; Patil, R.; Khan, H.; Jain, N.; et al. Genome-Wide Association Mapping for Component Traits of Drought and Heat Tolerance in Wheat. Front. Plant Sci. 2022, 13, 943033. [Google Scholar] [CrossRef]
  27. Breseghello, F.; Sorrells, M.E. Association Mapping of Kernel Size and Milling Quality in Wheat (Triticum aestivum L.) Cultivars. Genetics 2006, 172, 1165–1177. [Google Scholar] [CrossRef]
  28. Poland, J.A.; Bradbury, P.J.; Buckler, E.S.; Nelson, R.J. Genome-Wide Nested Association Mapping of Quantitative Resistance to Northern Leaf Blight in Maize. Proc. Natl. Acad. Sci. USA 2011, 108, 6893–6898. [Google Scholar] [CrossRef]
  29. Würschum, T.; Langer, S.M.; Longin, C.F.H.; Korzun, V.; Akhunov, E.; Ebmeyer, E.; Schachschneider, R.; Schacht, J.; Kazman, E.; Reif, J.C. Population Structure, Genetic Diversity and Linkage Disequilibrium in Elite Winter Wheat Assessed with SNP and SSR Markers. Theor. Appl. Genet. 2013, 126, 1477–1486. [Google Scholar] [CrossRef] [PubMed]
  30. Sehgal, D.; Mondal, S.; Guzman, C.; Garcia Barrios, G.; Franco, C.; Singh, R.; Dreisigacker, S. Validation of Candidate Gene-Based Markers and Identification of Novel Loci for Thousand-Grain Weight in Spring Bread Wheat. Front. Plant Sci. 2019, 10, 1189. [Google Scholar] [CrossRef]
  31. Liu, J.; He, Z.; Rasheed, A.; Wen, W.; Yan, J.; Zhang, P.; Wan, Y.; Zhang, Y.; Xie, C.; Xia, X. Genome-Wide Association Mapping of Black Point Reaction in Common Wheat (Triticum aestivum L.). BMC Plant Biol. 2017, 17, 220. [Google Scholar] [CrossRef]
  32. Eltaher, S.; Sallam, A.; Belamkar, V.; Emara, H.A.; Nower, A.A.; Salem, K.F.M.; Poland, J.; Baenziger, P.S. Genetic Diversity and Population Structure of F3:6 Nebraska Winter Wheat Genotypes Using Genotyping-by-Sequencing. Front. Genet. 2018, 9, 76. [Google Scholar] [CrossRef]
  33. Sheoran, S.; Jaiswal, S.; Kumar, D.; Raghav, N.; Sharma, R.; Pawar, S.; Paul, S.; Iquebal, M.A.; Jaiswar, A.; Sharma, P.; et al. Uncovering Genomic Regions Associated with 36 Agro-Morphological Traits in Indian Spring Wheat Using GWAS. Front. Plant Sci. 2019, 10, 527. [Google Scholar] [CrossRef] [PubMed]
  34. Ain, Q.U.; Rasheed, A.; Anwar, A.; Mahmood, T.; Imtiaz, M.; Mahmood, T.; Xia, X.; He, Z.; Quraishi, U.M. Genome-Wide Association for Grain Yield under Rainfed Conditions in Historical Wheat Cultivars from Pakistan. Front. Plant Sci. 2015, 6, 743. [Google Scholar] [CrossRef]
  35. Muleta, K.T.; Bulli, P.; Zhang, Z.; Chen, X.; Pumphrey, M. Unlocking Diversity in Germplasm Collections via Genomic Selection: A Case Study Based on Quantitative Adult Plant Resistance to Stripe Rust in Spring Wheat. Plant Genome 2017, 10, plantgenome2016.12.0124. [Google Scholar] [CrossRef]
  36. Bonman, J.M.; Babiker, E.M.; Cuesta-marcos, A.; Esvelt-klos, K.; Brown-guedira, G.; Chao, S.; See, D.; Chen, J.; Akhunov, E.; Zhang, J.; et al. Genetic Diversity among Wheat Accessions from the USDA National Small Grains Collection. Corp. Sci. 2015, 55, 1243–1253. [Google Scholar] [CrossRef]
  37. Balfourier, F.; Bouchet, S.; Robert, S.; De Oliveira, R.; Rimbert, H.; Kitt, J.; Choulet, F.; Wheat, I.; Sequencing, G. Worldwide Phylogeography and History of Wheat Genetic Diversity. Sci. Adv. 2019, 5, eaav0536. [Google Scholar] [CrossRef]
  38. Lehnert, H.; Berner, T.; Lang, D.; Beier, S.; Stein, N.; Himmelbach, A.; Kilian, B.; Keilwagen, J. Insights into Breeding History, Hotspot Regions of Selection, and Untapped Allelic Diversity for Bread Wheat Breeding. Plant J. 2022, 112, 897–918. [Google Scholar] [CrossRef] [PubMed]
  39. Sthapit, S.R.; Ruff, T.M.; Hooker, M.A.; See, D.R. Population Structure and Genetic Diversity of U.S. Wheat Varieties. Plant Genome 2022, 15, e20196. [Google Scholar] [CrossRef]
  40. Novoselović, D.; Bentley, A.R.; Šimek, R.; Dvojković, K.; Sorrells, M.E.; Gosman, N.; Horsnell, R.; Drezner, G.; Šatović, Z. Characterizing Croatian Wheat Germplasm Diversity and Structure in a European Context by DArT Markers. Front. Plant Sci. 2016, 7, 184. [Google Scholar] [CrossRef] [PubMed]
  41. Kiss, T.; Balla, K.; Cseh, A.; Berki, Z.; Horváth, Á.; Vida, G.; Veisz, O.; Karsai, I. Assessment of the Genetic Diversity, Population Structure and Allele Distribution of Major Plant Development Genes in Bread Wheat Cultivars Using DArT and Gene-Specific Markers. Cereal Res. Commun. 2021, 49, 549–557. [Google Scholar] [CrossRef]
  42. Iqbal, M.; Shahzad, A.; Ahmed, I. Allelic Variation at the Vrn-A1, Vrn-B1, Vrn-D1, Vrn-B3 and Ppd-D1a Loci of Pakistani Spring Wheat Cultivars. Electron. J. Biotechnol. 2011, 14, 1–8. [Google Scholar] [CrossRef]
  43. Yang, F.P.; Zhang, X.K.; Xia, X.C.; Laurie, D.A.; Yang, W.X.; He, Z.H. Distribution of the Photoperiod Insensitive Ppd-D1a Allele in Chinese Wheat Cultivars. Euphytica 2009, 165, 445–452. [Google Scholar] [CrossRef]
  44. Andrey, B.S.; Börner, A.; Salina, E.A. Effect of VRN-1 and PPD-D1 Genes on Heading Time in European Bread Wheat Cultivars. Plant Breed. 2015, 134, 49–55. [Google Scholar] [CrossRef]
  45. Kiss, T.; Balla, K.; Veisz, O.; Láng, L.; Bedő, Z.; Griffiths, S.; Isaac, P.; Karsai, I. Allele Frequencies in the VRN-A1, VRN-B1 and VRN-D1 Vernalization Response and PPD-B1 and PPD-D1 Photoperiod Sensitivity Genes, and Their Effects on Heading in a Diverse Set of Wheat Cultivars (Triticum aestivum L.). Mol. Breed. 2014, 34, 297–310. [Google Scholar] [CrossRef]
  46. Sanchez-Garcia, M.; Bentley, A. Global Journeys of Adaptive Wheat Genes; Woodhead Publishing: Cambridge, UK, 2019; ISBN 9780081021637. [Google Scholar]
  47. Royo, C.; Dreisigacker, S.; Soriano, J.M.; Lopes, M.S.; Ammar, K.; Villegas, D. Allelic Variation at the Vernalization Response (Vrn-1) and Photoperiod Sensitivity (Ppd-1) Genes and Their Association with the Development of Durum Wheat Landraces and Modern Cultivars. Front. Plant Sci. 2020, 11, 838. [Google Scholar] [CrossRef] [PubMed]
  48. Dixon, L.E.; Karsai, I.; Kiss, T.; Adamski, N.M.; Liu, Z.; Ding, Y.; Allard, V.; Boden, S.A.; Griffiths, S. VERNALIZATION1 Controls Developmental Responses of Winter Wheat under High Ambient Temperatures. Development 2019, 146, dev172684. [Google Scholar] [CrossRef] [PubMed]
  49. Ortiz Ferrara, G.; Mosaad, M.G.; Mahalakshmi, V.; Rajaram, S. Photoperiod and Vernalisation Response of Mediterranean Wheats, and Implications for Adaptation. Euphytica 1998, 100, 377–384. [Google Scholar] [CrossRef]
  50. Díaz, A.; Zikhali, M.; Turner, A.S.; Isaac, P.; Laurie, D.A. Copy Number Variation Affecting the Photoperiod-B1 and Vernalization-A1 Genes Is Associated with Altered Flowering Time in Wheat (Triticum aestivum). PLoS ONE 2012, 7, e33234. [Google Scholar] [CrossRef]
  51. Chen, A.; Dubcovsky, J. Wheat TILLING Mutants Show That the Vernalization Gene VRN1 Down-Regulates the Flowering Repressor VRN2 in Leaves but Is Not Essential for Flowering. PLoS Genet. 2012, 8, e1003134. [Google Scholar] [CrossRef]
  52. Eagles, H.; Cane, K.; Trevaskisc, B. Veery Wheats Carry an Allele of Vrn-A1 That Has Implications for Freezing Tolerance in Winter Wheats. Plant Breed. 2011, 130, 413–418. [Google Scholar] [CrossRef]
  53. Zhang, X.K.; Xiao, Y.G.; Zhang, Y.; Xia, X.C.; Dubcovsky, J.; He, Z.H. Allelic Variation at the Vernalization Genes Vrn-A1, Vrn-B1, Vrn-D1, and Vrn-B3 in Chinese Wheat Cultivars and Their Association with Growth Habit. Crop Sci. 2008, 48, 458–470. [Google Scholar] [CrossRef]
  54. Würschum, T.; Boeven, P.H.G.; Langer, S.M.; Longin, C.F.H.; Leiser, W.L. Multiply to Conquer: Copy Number Variations at Ppd-B1 and Vrn-A1 Facilitate Global Adaptation in Wheat. BMC Genet. 2015, 16, 96. [Google Scholar] [CrossRef]
  55. Li, G.; Yu, M.; Fang, T.; Cao, S.; Carver, B.F.; Yan, L. Vernalization Requirement Duration in Winter Wheat Is Controlled by TaVRN-A1 at the Protein Level. Plant J. 2013, 76, 742–753. [Google Scholar] [CrossRef]
  56. Langer, S.M.; Longin, C.F.H.; Würschum, T. Flowering Time Control in European Winter Wheat. Front. Plant Sci. 2014, 5, 537. [Google Scholar] [CrossRef] [PubMed]
  57. Kamran, A.; Iqbal, M.; Spaner, D. Flowering Time in Wheat (Triticum aestivum L.): A Key Factor for Global Adaptability. Euphytica 2014, 197, 1–26. [Google Scholar] [CrossRef]
  58. Würschum, T.; Langer, S.M.; Longin, C.F.H.; Tucker, M.R.; Leiser, W.L. A Three-Component System Incorporating Ppd-D1, Copy Number Variation at Ppd-B1, and Numerous Small-Effect Quantitative Trait Loci Facilitates Adaptation of Heading Time in Winter Wheat Cultivars of Worldwide Origin. Plant Cell Environ. 2018, 41, 1407–1416. [Google Scholar] [CrossRef] [PubMed]
  59. Tanio, M.; Kato, K. Development of Near-Isogenic Lines for Photoperiod-Insensitive Genes, Ppd-B1 and Ppd-D1, Carried by the Japanese Wheat Cultivars and Their Effect on Apical Development. Breed. Sci. 2007, 57, 65–72. [Google Scholar] [CrossRef]
  60. Worland, A.J. The Influence of Flowering Time Genes on Environmental Adaptability in European Wheats. Euphytica 1996, 89, 49–57. [Google Scholar] [CrossRef]
  61. Bentley, A.R.; Horsnell, R.; Werner, C.P.; Turner, A.S.; Rose, G.A.; Bedard, C.; Howell, P.; Wilhelm, E.P. Short, Natural, and Extended Photoperiod Response in BC 2 F 4 Lines of Bread Wheat with Different Photoperiod-1 (Ppd-1) Alleles. J. Exp. Bot. 2013, 64, 1783–1793. [Google Scholar] [CrossRef] [PubMed]
  62. Worland, A.J.; Bo, A. The Influence of Photoperiod Genes on the Adaptability of European Winter Wheats. Euphytica 1998, 100, 385–394. [Google Scholar] [CrossRef]
  63. Whittal, A.; Kaviani, M.; Graf, R.; Humphreys, G.; Navabi, A. Allelic Variation of Vernalization and Photoperiod Response Genes in a Diverse Set of North American High Latitude Winter Wheat Genotypes. PLoS ONE 2018, 13, e0209543. [Google Scholar] [CrossRef] [PubMed]
  64. Grogan, S.M.; Brown-guedira, G.; Haley, S.D.; Mcmaster, G.S.; Reid, D.; Smith, J.; Byrne, P.F. Allelic Variation in Developmental Genes and Effects on Winter Wheat Heading Date in the U.S. Great Plains. PLoS ONE 2016, 11, e0152852. [Google Scholar] [CrossRef] [PubMed]
  65. Tadesse, W.; Ogbonnaya, F.C.; Jighly, A.; Sanchez-Garcia, M.; Sohail, Q.; Rajaram, S.; Baum, M. Genome-Wide Association Mapping of Yield and Grain Quality Traits in Winter Wheat Genotypes. PLoS ONE 2015, 10, e0141339. [Google Scholar] [CrossRef] [PubMed]
  66. Monostori, I.; Szira, F.; Tondelli, A.; ÁrendÁs, T.; Gierczik, K.; Cattivelli, L.; Galiba, G.; VÁgújfalvi, A. Genome-Wide Association Study and Genetic Diversity Analysis on Nitrogen Use Efficiency in a Central European Winter Wheat (Triticum aestivum L.) Collection. PLoS ONE 2017, 12, e0189265. [Google Scholar] [CrossRef]
  67. Li, Z.; Lhundrup, N.; Guo, G.; Dol, K.; Chen, P.; Gao, L.; Chemi, W.; Zhang, J.; Wang, J.; Nyema, T.; et al. Characterization of Genetic Diversity and Genome-Wide Association Mapping of Three Agronomic Traits in Qingke Barley (Hordeum vulgare L.) in the Qinghai-Tibet Plateau. Front. Genet. 2020, 11, 638. [Google Scholar] [CrossRef]
  68. El Hanafi, S.; Backhaus, A.; Bendaou, N.; Sanchez-Garcia, M.; Al-Abdallat, A.; Tadesse, W. Genome-Wide Association Study for Adult Plant Resistance to Yellow Rust in Spring Bread Wheat (Triticum aestivum L.). Euphytica 2021, 217, 87. [Google Scholar] [CrossRef]
  69. Kobayashi, F.; Tanaka, T.; Kanamori, H.; Wu, J.; Katayose, Y.; Handa, H. Characterization of a Mini Core Collection of Japanese Wheat Varieties Using Singlenucleotide Polymorphisms Generated by Genotyping-by-Sequencing. Breed. Sci. 2016, 66, 213–225. [Google Scholar] [CrossRef] [PubMed]
  70. Ogbonnaya, F.C.; Rasheed, A.; Okechukwu, E.C.; Jighly, A.; Makdis, F.; Wuletaw, T.; Hagras, A.; Uguru, M.I.; Agbo, C.U. Genome-Wide Association Study for Agronomic and Physiological Traits in Spring Wheat Evaluated in a Range of Heat Prone Environments. Theor. Appl. Genet. 2017, 130, 1819–1835. [Google Scholar] [CrossRef] [PubMed]
  71. Sheoran, S.; Jaiswal, S.; Raghav, N.; Sharma, R.; Sabhyata; Gaur, A.; Jaisri, J.; Tandon, G.; Singh, S.; Sharma, P.; et al. Genome-Wide Association Study and Post-Genome-Wide Association Study Analysis for Spike Fertility and Yield Related Traits in Bread Wheat. Front. Plant Sci. 2022, 12, 820761. [Google Scholar] [CrossRef]
  72. Aslam, M.A.; Ahmed, M.; Stöckle, C.O.; Higgins, S.S.; ul Hassan, F.; Hayat, R. Can Growing Degree Days and Photoperiod Predict Spring Wheat Phenology? Front. Environ. Sci. 2017, 5, 57. [Google Scholar] [CrossRef]
  73. Zhang, Y.; Liu, W.C.; Li, J.; Wei, H.T.; Hu, X.R.; Li, Y.J.; Lu, B.R.; Yang, W.Y. Distribution and Selective Effects of Vrn-A1, Vrn-B1, and Vrn-D1 Genes in Derivative Varieties from Four Cornerstone Breeding Parents of Wheat in China. Agric. Sci. China 2010, 9, 1389–1399. [Google Scholar] [CrossRef]
  74. Kamran, A.; Iqbal, M.; Navabi, A.; Randhawa, H.; Pozniak, C.; Spaner, D. Earliness per Se QTLs and Their Interaction with the Photoperiod Insensitive Allele Ppd-D1a in the Cutler × AC Barrie Spring Wheat Population. Theor. Appl. Genet. 2013, 126, 1965–1976. [Google Scholar] [CrossRef] [PubMed]
  75. Faricelli, M.E.; Valárik, M.; Dubcovsky, J. Control of Flowering Time and Spike Development in Cereals: The Earliness per Se Eps-1 Region in Wheat, Rice, and Brachypodium. Funct. Integr. Genom. 2010, 10, 293–306. [Google Scholar] [CrossRef] [PubMed]
  76. Allen, A.M.; Winfield, M.O.; Burridge, A.J.; Downie, R.C.; Benbow, H.R.; Barker, G.L.A.; Wilkinson, P.A.; Coghill, J.; Waterfall, C.; Davassi, A.; et al. Characterization of a Wheat Breeders’ Array Suitable for High-Throughput SNP Genotyping of Global Accessions of Hexaploid Bread Wheat (Triticum aestivum). Plant Biotechnol. J. 2017, 15, 390–401. [Google Scholar] [CrossRef]
  77. Dreisigacke, S.; Sukumaran, S.; Guzmán, C.; He, X.; Lan, C.; Bonnett, D.; Crossa, J. Molecular Marker-Based Selection Tools in Spring Bread Wheat Improvement: CIMMYT Experience and Prospects; Springer: Berlin/Heidelberg, Germany, 2016; Volume 11, ISBN 978-3-319-27088-3. [Google Scholar]
  78. Fu, D.; Szűcs, P.; Yan, L.; Helguera, M.; Skinner, J.S.; Von Zitzewitz, J.; Hayes, P.M.; Dubcovsky, J. Large Deletions within the First Intron in VRN-1 Are Associated with Spring Growth Habit in Barley and Wheat. Mol. Genet. Genomics 2005, 273, 54–65. [Google Scholar] [CrossRef]
  79. Beales, J.; Turner, A.; Griffiths, S.; Snape, J.W.; Laurie, D.A. A Pseudo-Response Regulator Is Misexpressed in the Photoperiod Insensitive Ppd-D1a Mutant of Wheat (Triticum aestivum L.). Theor. Appl. Genet. 2007, 115, 721–733. [Google Scholar] [CrossRef]
  80. Lilin-yin Package ‘CMplot’. 2022. pp. 1–8. Available online: https://github.com/YinLiLin/CMplot (accessed on 1 October 2023).
  81. Remington, D.L.; Thornsberry, J.M.; Matsuoka, Y.; Wilson, L.M.; Whitt, S.R.; Doebley, J.; Kresovich, S.; Goodman, M.M.; Buckler IV, E.S. Structure of Linkage Disequilibrium and Phenotypic Associations in the Maize Genome. Proc. Natl. Acad. Sci. USA 2001, 98, 11479–11484. [Google Scholar] [CrossRef]
  82. Yu, G.; Smith, D.K.; Zhu, H.; Guan, Y.; Lam, T.T.Y. Ggtree: An R Package for Visualization and Annotation of Phylogenetic Trees with Their Covariates and Other Associated Data. Methods Ecol. Evol. 2017, 8, 28–36. [Google Scholar] [CrossRef]
  83. Bradbury, P.J.; Zhang, Z.; Kroon, D.E.; Casstevens, T.M.; Ramdoss, Y.; Buckler, E.S. TASSEL: Software for Association Mapping of Complex Traits in Diverse Samples. Bioinformatics 2007, 23, 2633–2635. [Google Scholar] [CrossRef]
  84. Frichot, E.; François, O. LEA: An R Package for Landscape and Ecological Association Studies. Methods Ecol. Evol. 2015, 6, 925–929. [Google Scholar] [CrossRef]
  85. Liu, K.; Muse, S.V. PowerMaker: An Integrated Analysis Environment for Genetic Maker Analysis. Bioinformatics 2005, 21, 2128–2129. [Google Scholar] [CrossRef]
  86. Wang, E.; Engel, T. Simulation of Phenological Development of Wheat Crops. Agric. Syst. 1998, 58, 1–24. [Google Scholar] [CrossRef]
  87. White, J.W.; Laing, D.R. Photoperiod Response of Flowering in Diverse Genotypes of Common Bean (Phaseolus vulgaris). F. Crops Res. 1989, 22, 113–128. [Google Scholar] [CrossRef]
  88. Van Rossum, B.-J. StatgenGxE: Genotype by Environment (GxE) Analysis, Version 1.0.5. 2022. Available online: https://cran.r-project.org/web/packages/statgenGxE/statgenGxE.pdf (accessed on 1 October 2023).
  89. Rodríguez-Álvarez, M.X.; Boer, M.P.; van Eeuwijk, F.A.; Eilers, P.H.C. Correcting for Spatial Heterogeneity in Plant Breeding Experiments with P-Splines. Spat. Stat. 2018, 23, 52–71. [Google Scholar] [CrossRef]
  90. Rodríguez-Álvarez, M.X.; Boer, M.P.; van Eeuwijk, F.A.; Eilers, P.H.C. Spatial Models for Field Trials. arXiv 2016, arXiv:1607.08255. [Google Scholar]
  91. Tossim, H.A.; Nguepjop, J.R.; Diatta, C.; Sambou, A.; Seye, M.; Sane, D.; Rami, J.F.; Fonceka, D. Assessment of 16 Peanut (Arachis hypogaea L.) CSSLs Derived from an Interspecific Cross for Yield and Yield Component Traits: QTL Validation. Agronomy 2020, 10, 583. [Google Scholar] [CrossRef]
  92. Wang, J.; Zhang, Z. GAPIT Version 3: Boosting Power and Accuracy for Genomic Association and Prediction. Genom. Proteom. Bioinforma. 2021, 19, 629–640. [Google Scholar] [CrossRef]
  93. Price, A.L.; Patterson, N.J.; Plenge, R.M.; Weinblatt, M.E.; Shadick, N.A.; Reich, D. Principal Components Analysis Corrects for Stratification in Genome-Wide Association Studies. Nat. Genet. 2006, 38, 904–909. [Google Scholar] [CrossRef]
  94. Zhang, Z.; Ersoz, E.; Lai, C.Q.; Todhunter, R.J.; Tiwari, H.K.; Gore, M.A.; Bradbury, P.J.; Yu, J.; Arnett, D.K.; Ordovas, J.M.; et al. Mixed Linear Model Approach Adapted for Genome-Wide Association Studies. Nat. Genet. 2010, 42, 355–360. [Google Scholar] [CrossRef]
  95. Li, M.; Liu, X.; Bradbury, P.; Yu, J.; Zhang, Y.M.; Todhunter, R.J.; Buckler, E.S.; Zhang, Z. Enrichment of Statistical Power for Genome-Wide Association Studies. BMC Biol. 2014, 12, 73. [Google Scholar] [CrossRef]
  96. VanRaden, P.M. Efficient Methods to Compute Genomic Predictions. J. Dairy Sci. 2008, 91, 4414–4423. [Google Scholar] [CrossRef] [PubMed]
  97. Lipka, A.E.; Tian, F.; Wang, Q.; Peiffer, J.; Li, M.; Bradbury, P.J.; Gore, M.A.; Buckler, E.S.; Zhang, Z. GAPIT: Genome Association and Prediction Integrated Tool. Bioinformatics 2012, 28, 2397–2399. [Google Scholar] [CrossRef] [PubMed]
  98. Chen, H.; Wang, C.; Conomos, M.P.; Stilp, A.M.; Li, Z.; Sofer, T.; Szpiro, A.A.; Chen, W.; Brehm, J.M.; Celedón, J.C.; et al. Control for Population Structure and Relatedness for Binary Traits in Genetic Association Studies via Logistic Mixed Models. Am. J. Hum. Genet. 2016, 98, 653–666. [Google Scholar] [CrossRef] [PubMed]
Plants 12 04053 g001
Figure 1. The genetic relationships and population structure of the Winter Wheat Association Genetics Initiative (WWAGI) with 249 advanced lines and varieties using principal component analysis (PCA) of SNP markers concluded via PCA and sparse Non-negative Matrix Factorization (sNMF) analyses. Spatial distribution of the genotypes is based on geographic origin (A) and assignment to sNMF groups at K = 2 (B), K = 3 (C). The phylogenetic tree of the WWAGI represents a clustering of individuals, through respective layers (from inside to outside), geographical region, and assignment to sNMF groups at K = 2 and K = 3. Genotypes with the highest ancestry coefficient below 70% are colored in gray (D).
Figure 1. The genetic relationships and population structure of the Winter Wheat Association Genetics Initiative (WWAGI) with 249 advanced lines and varieties using principal component analysis (PCA) of SNP markers concluded via PCA and sparse Non-negative Matrix Factorization (sNMF) analyses. Spatial distribution of the genotypes is based on geographic origin (A) and assignment to sNMF groups at K = 2 (B), K = 3 (C). The phylogenetic tree of the WWAGI represents a clustering of individuals, through respective layers (from inside to outside), geographical region, and assignment to sNMF groups at K = 2 and K = 3. Genotypes with the highest ancestry coefficient below 70% are colored in gray (D).
Plants 12 04053 g001
Plants 12 04053 g002
Figure 2. (A) Phylogeny-based clustering tree represents a clustering of individuals, through respective layers (from inside to outside), geographical region, and Vrn and Ppd allelic combination. (B,C) Characterization of the genotypes used based on gene diversity for Ppd and Vrn.
Figure 2. (A) Phylogeny-based clustering tree represents a clustering of individuals, through respective layers (from inside to outside), geographical region, and Vrn and Ppd allelic combination. (B,C) Characterization of the genotypes used based on gene diversity for Ppd and Vrn.
Plants 12 04053 g002
Plants 12 04053 g003
Figure 3. Unraveling regulatory flowering genes, vernalization response, and photoperiod sensitivity of the Winter Wheat Association Genetics Initiative (WWAGI). (A) Impact of Vrn Genes on percentage of flowering in spring and winter Alleles. (B) Impact of different winter types alleles on percentage of flowering. (C) Impact of Vrn Genes on average extra days to flowering in Spring and Winter Alleles under no vernalization. (D) Impact of different winter types alleles on average extra days to flowering under no vernalization. (E) Impact of Ppd-B1/Ppd-D1 combination alleles on percentage of flowering under short photoperiod. (F) Impact of Ppd-B1/Ppd-D1 combination alleles on average extra days to flowering under short photoperiod.
Figure 3. Unraveling regulatory flowering genes, vernalization response, and photoperiod sensitivity of the Winter Wheat Association Genetics Initiative (WWAGI). (A) Impact of Vrn Genes on percentage of flowering in spring and winter Alleles. (B) Impact of different winter types alleles on percentage of flowering. (C) Impact of Vrn Genes on average extra days to flowering in Spring and Winter Alleles under no vernalization. (D) Impact of different winter types alleles on average extra days to flowering under no vernalization. (E) Impact of Ppd-B1/Ppd-D1 combination alleles on percentage of flowering under short photoperiod. (F) Impact of Ppd-B1/Ppd-D1 combination alleles on average extra days to flowering under short photoperiod.
Plants 12 04053 g003
Plants 12 04053 g004
Figure 4. Phenology characterization of the winter wheat panel based on vernalization–photoperiod response requirement.
Figure 4. Phenology characterization of the winter wheat panel based on vernalization–photoperiod response requirement.
Plants 12 04053 g004
Plants 12 04053 g005
Figure 5. Phenology characterization of the winter wheat panel based on region and vernalization–photoperiod response requirement.
Figure 5. Phenology characterization of the winter wheat panel based on region and vernalization–photoperiod response requirement.
Plants 12 04053 g005
Plants 12 04053 g006
Figure 6. Relationship between relative response to vernalization (RRV), relative response to photoperiod (RRP), and earliness per se (Eps). (A) Regression between Eps and RRV. (B) Regression between Eps and RRP. (C) Regression between RRV and RRP. (D) Boxplot analysis of the Eps under vernalization treatment. (E) Boxplot analysis of the Eps under photoperiod treatment.
Figure 6. Relationship between relative response to vernalization (RRV), relative response to photoperiod (RRP), and earliness per se (Eps). (A) Regression between Eps and RRV. (B) Regression between Eps and RRP. (C) Regression between RRV and RRP. (D) Boxplot analysis of the Eps under vernalization treatment. (E) Boxplot analysis of the Eps under photoperiod treatment.
Plants 12 04053 g006
Plants 12 04053 g007
Figure 7. PCA based on geographic region (A), earliness per se (B), relative response to photoperiod (C), relative response to vernalization (D), complete vernalization (E), complete photosensitive (F), and vernalization and photoperiod requirement (G).
Figure 7. PCA based on geographic region (A), earliness per se (B), relative response to photoperiod (C), relative response to vernalization (D), complete vernalization (E), complete photosensitive (F), and vernalization and photoperiod requirement (G).
Plants 12 04053 g007
Plants 12 04053 g008
Figure 8. Genome-wide association scan for RRV, RRP, Eps, complete photosensitive, and complete winter in WWAGI. The plots show Manhattan plots with the names of the significant SNPs. The chromosomes are shown on the X-axis, and the genome-wide scan −log10 (p-values) values are shown on the Y-axis using common threshold (−log10p > 3) in black dotted line and Bonferroni correction (−log10p > 5.2) in red dotted line.
Figure 8. Genome-wide association scan for RRV, RRP, Eps, complete photosensitive, and complete winter in WWAGI. The plots show Manhattan plots with the names of the significant SNPs. The chromosomes are shown on the X-axis, and the genome-wide scan −log10 (p-values) values are shown on the Y-axis using common threshold (−log10p > 3) in black dotted line and Bonferroni correction (−log10p > 5.2) in red dotted line.
Plants 12 04053 g008
Plants 12 04053 g009
Figure 9. Characterization of the Winter Wheat Association Genetics Initiative (WWAGI) panel with 249 advanced lines and varieties used based on region and geographic origin (A,B).
Figure 9. Characterization of the Winter Wheat Association Genetics Initiative (WWAGI) panel with 249 advanced lines and varieties used based on region and geographic origin (A,B).
Plants 12 04053 g009
Table 1. Major allele frequency, sample size, genetic diversity, and the polymorphism information content (PIC) of Chromosome, Region of the Winter Wheat Association Genetics Initiative (WWAGI) with 249 advanced lines and varieties.
Table 1. Major allele frequency, sample size, genetic diversity, and the polymorphism information content (PIC) of Chromosome, Region of the Winter Wheat Association Genetics Initiative (WWAGI) with 249 advanced lines and varieties.
Major Allele
Frequency		Number
of SNPs		Genetic
Diversity		PIC
Chromosome
1A0.705340.430.37
1B0.707100.420.36
1D0.702210.410.35
2A0.736610.380.33
2B0.687870.460.39
2D0.682340.430.36
3A0.704510.430.37
3B0.677770.470.40
3D0.651060.470.40
4A0.723830.40.35
4B0.703080.430.37
4D0.67320.440.36
5A0.665960.460.39
5B0.728150.410.36
5D0.691370.420.36
6A0.695030.440.38
6B0.686740.440.38
6D0.72760.40.35
7A0.705250.440.38
7B0.684700.450.38
7D0.72830.40.34
Major Allele
Frequency		Sample
Size		Genetic
Diversity		PIC
Region
Central Asia0.8190.270.22
CGIAR0.72650.390.35
Europe0.74480.380.34
North America0.72620.390.35
Other0.75110.330.28
Western Asia0.70540.420.37
Table 2. AMOVA of the Winter Wheat Association Genetics Initiative (WWAGI) with 249 advanced lines and varieties based on geographical region.
Table 2. AMOVA of the Winter Wheat Association Genetics Initiative (WWAGI) with 249 advanced lines and varieties based on geographical region.
Region
SourcedfMean SqSum SqPercentage
Among Region528,330141,6526%
Within Central Asia8684354,7482%
Within CGIAR648921570,93126%
Within Europe478685408,21118%
Within North America618413513,16423%
Within Other10383038,2982%
Within Western Asia539118483,26722%
Table 3. Population differentiation of the Winter Wheat Association Genetics Initiative (WWAGI) with 249 advanced lines and varieties calculated as the Nei genetic distance according to geographical region.
Table 3. Population differentiation of the Winter Wheat Association Genetics Initiative (WWAGI) with 249 advanced lines and varieties calculated as the Nei genetic distance according to geographical region.
RegionCentral AsiaCGIAREuropeNorth AmericaOtherWestern Asia
Central Asia0.000.050.070.080.100.06
CGIAR0.050.000.030.030.060.01
Europe0.070.030.000.040.060.03
North America0.080.030.040.000.050.04
Other0.100.060.060.050.000.06
Western Asia0.060.010.030.040.060.00
Table 4. Phenology experiment conditions during the three vernalization and photoperiod treatments.
Table 4. Phenology experiment conditions during the three vernalization and photoperiod treatments.
ExperimentVernalization (Days at 5 °C)Photoperiod (Hour in Days/Night)
V+P+4920/4
V−P+020/4
V+P−4914/10
Table 5. Descriptive statistic of growing degree days (GDD) for the WWAGI used based on vernalization and photoperiod requirement experiment.
Table 5. Descriptive statistic of growing degree days (GDD) for the WWAGI used based on vernalization and photoperiod requirement experiment.
TreatmentMeanMinMaxVarianceGenetic VarianceHeritability
V+P+1273.131013.91629.4218,120.5916,994.020.92
V−P+1417.13896.852241.87138,772.3124,085.40.91
V+P−1705.051270215035,540.1525,163.10.81
Table 6. ANOVA of growing degree days (GDD) for the WWAGI genotypes used based on vernalization and photoperiod requirement experiment (treatment).
Table 6. ANOVA of growing degree days (GDD) for the WWAGI genotypes used based on vernalization and photoperiod requirement experiment (treatment).
DfSum sqMean sqF ValuePr (>F)
Treatment220,430,31310,215,1572681.48<2.2 × 10−16 ***
Genotype23412,550,47353,63514.0794.459 × 10−16 ***
GenotypeXTreatment25511,103,56243,54311.431.820 × 10−14 ***
Residuals40152,3813810
*** significant at the probability of 0.01%.
Table 7. Significant marker associations for each trait under study.
Table 7. Significant marker associations for each trait under study.
Trait 1MarkerChromosomePosition−log10(p)R2
Strict photosensitivePpd-D1 *2D3620749711.820.41
Strict photosensitiveAX-944069572D350691434.50.31
Strict photosensitiveAX-946569887D1976657533.460.25
Strict photosensitiveAX-943834047B1374853543.120.25
Strict vernalizationVrn-B1 *5B5770152007.750.36
Strict vernalizationAX-951663975B5835352233.760.27
Strict vernalizationAX-951946331A5006795313.370.28
Strict vernalizationAX-949130415A5981903023.310.27
Strict vernalizationAX-949187695B5897625863.30.26
EpsVrn-B15B5770152003.740.12
RRPAX-945234055D3980198823.680.1
RRPAX-946697685A3775572153.040.09
RRVAX-947964795A5902259395.090.49
RRVVrn-A15A5903995784.050.44
RRVAX-944042195A5907029923.480.42
RRVAX-949264195B5787194443.310.41
RRVAX-946796922D5729337853.080.4
1 Relative response to vernalization (RRV), relative response to photoperiod (RRP), and earliness per se (Eps). * Significant marker that passes Bonferroni correction.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

El Baouchi, A.; Ibriz, M.; Dreisigacker, S.; Lopes, M.S.; Garcia, M.S. Genetic Diversity and Genome-Wide Association Study for the Phenology Response of Winter Wheats of North America, Western Asia, and Europe. Plants 2023, 12, 4053. https://doi.org/10.3390/plants12234053

AMA Style

El Baouchi A, Ibriz M, Dreisigacker S, Lopes MS, Garcia MS. Genetic Diversity and Genome-Wide Association Study for the Phenology Response of Winter Wheats of North America, Western Asia, and Europe. Plants. 2023; 12(23):4053. https://doi.org/10.3390/plants12234053

Chicago/Turabian Style

El Baouchi, Adil, Mohammed Ibriz, Susanne Dreisigacker, Marta S. Lopes, and Miguel Sanchez Garcia. 2023. "Genetic Diversity and Genome-Wide Association Study for the Phenology Response of Winter Wheats of North America, Western Asia, and Europe" Plants 12, no. 23: 4053. https://doi.org/10.3390/plants12234053

APA Style

El Baouchi, A., Ibriz, M., Dreisigacker, S., Lopes, M. S., & Garcia, M. S. (2023). Genetic Diversity and Genome-Wide Association Study for the Phenology Response of Winter Wheats of North America, Western Asia, and Europe. Plants, 12(23), 4053. https://doi.org/10.3390/plants12234053

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