Genetic control and prospects of predictive breeding for European winter wheat's Zeleny sedimentation values and Hagberg-Perten falling number.

KEY MESSAGE: Sedimentation values and falling number in the last decades have helped maintain high baking quality despite rigorous selection for grain yield in wheat. Allelic combinations of major loci sustained the bread-making quality while improving grain yield. Glu-D1, Pinb-D1, and non-gluten proteins are associated with sedimentation values and falling number in European wheat. Zeleny sedimentation values (ZSV) and Hagberg-Perten falling number (HFN) are among the most important parameters that help determine the baking quality classes of wheat and, thus, influence the monetary benefits for growers. We used a published data set of 372 European wheat varieties evaluated in replicated field trials in multiple environments. ZSV and HFN traits hold a wide and significant genotypic variation and high broad-sense heritability. The genetic correlations revealed positive and significant associations of ZSV and HFN with each other, grain protein content (GPC) and grain hardness; however, they were all significantly negatively correlated with grain yield. Besides, GPC appeared to be the major predictor for ZSV and HFN. Our genome-wide association analyses based on high-quality SSR, SNP, and candidate gene markers revealed a strong quantitative genetic nature of ZSV and HFN by explaining their total genotypic variance as 41.49% and 38.06%, respectively. The association of known Glutenin (Glu-1) and Puroindoline (Pin-1) with ZSV provided positive analytic proof of our studies. We report novel candidate loci associated with globulins and albumins-the non-gluten monomeric proteins in wheat. In addition, predictive breeding analyses for ZSV and HFN suggest using genomic selection in the early stages of breeding programs with an average prediction accuracy of 81 and 59%, respectively.

Gabi wheat a panel of European elite lines as central stock for wheat genetic research.

In plant sciences, curation and availability of interoperable phenotypic and genomic data is still in its infancy and represents an obstacle to rapid scientific discoveries in this field. To that end, supplementing the efforts being made to generate open access wheat genome, pan wheat genome and other bioinformatic resources, we present the GABI-WHEAT panel of elite European cultivars comprising 358 winter and 14 summer wheat varieties released between 1975 to 2007. The panel has been genotyped with SNP arrays of increasing density to investigate several important agronomic, quality and disease resistance traits. The robustness of investigated traits and interoperability of genomic and phenotypic data was assessed in the current publication with the aim to transform this panel into a public data resource for future genetic research in wheat. Consecutively, the phenotypic data was formatted to comply with FAIR principles and linked to online databases to substantiate panel origin information and quality. Thus, we were able to make a valuable resource available for plant science in a sustainable way.

QTL mapping for resistance against cereal cyst nematode (Heterodera avenae Woll.) in wheat (Triticum aestivum L.).

The resistance to cereal cyst nematode (Heterodera avenae Woll.) in wheat (Triticum aestivum L.) was studied using 114 doubled haploid lines from a novel ITMI mapping population. These lines were screened for nematode infestation in a controlled environment for two years. QTL-mapping analyses were performed across two years (Y1 and Y2) as well as combining two years (CY) data. On the 114 lines that were screened, a total of 2,736 data points (genotype, batch or years, and replication combinations) were acquired. For QTL analysis, 12,093 markers (11,678 SNPs and 415 SSRs markers) were used, after filtering the genotypic data, for the QTL mapping. Composite interval mapping, using Haley-Knott regression (hk) method in R/QTL, was used for QTL analysis. In total, 19 QTLs were detected out of which 13 were novel and six were found to be colocalized or nearby to previously reported Cre genes, QTLs or MTAs for H. avenae or H. filipjevi. Nine QTLs were detected across all three groups (Y1, Y2 and CY) including a significant QTL "QCcn.ha-2D" on chromosome 2D that explains 23% of the variance. This QTL colocalized with a previously identified Cre3 locus. Novel QTL, QCcn.ha-2A, detected in the present study could be the possible unreported homeoloci to QCcn.ha-2D, QCcn.ha-2B.1 and QCcn.ha-2B.2. Six significant digenic epistatic interactions were also observed. In addition, 26 candidate genes were also identified including genes known for their involvement in PPNs (plant parasitic nematodes) resistance in different plant species. In-silico expression of putative candidate genes showed differential expression in roots during specific developmental stages. Results obtained in the present study are useful for wheat breeding to generate resistant genetic resources against H. avenae.

A SNP-based genetic dissection of versatile traits in bread wheat (Triticum aestivum L.).

The continuous increase in global population prompts increased wheat production. Future wheat (Triticum aestivum L.) breeding will heavily rely on dissecting molecular and genetic bases of wheat yield and related traits which is possible through the discovery of quantitative trait loci (QTLs) in constructed populations, such as recombinant inbred lines (RILs). Here, we present an evaluation of 92 RILs in a bi-parental RIL mapping population (the International Triticeae Mapping Initiative Mapping Population [ITMI/MP]) using newly generated phenotypic data in 3-year experiments (2015), older phenotypic data (1997-2009), and newly created single nucleotide polymorphism (SNP) marker data based on 92 of the original RILs to search for novel and stable QTLs. Our analyses of more than 15 unique traits observed in multiple experiments included analyses of 46 traits in three environments in the USA, 69 traits in eight environments in Germany, 149 traits in 10 environments in Russia, and 28 traits in four environments in India (292 traits in 25 environments) with 7584 SNPs (292 × 7584 = 2 214 528 data points). A total of 874 QTLs were detected with limit of detection (LOD) scores of 2.01-3.0 and 432 QTLs were detected with LOD > 3.0. Moreover, 769 QTLs could be assigned to 183 clusters based on the common markers and relative proximity of related QTLs, indicating gene-rich regions throughout the A, B, and D genomes of common wheat. This upgraded genotype-phenotype information of ITMI/MP can assist breeders and geneticists who can make crosses with suitable RILs to improve or investigate traits of interest.

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

Tan spot, caused by the fungus Pyrenophoratritici-repentis (Ptr), is a severe foliar disease of wheat (Triticumaestivum 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 single sequence repeats (SSR) in addition to 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 resistance 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.

Linkage mapping identifies a non-synonymous mutation in FLOWERING LOCUS T (FT-B1) increasing spikelet number per spike.

Total spikelet number per spike (TSN) is a major component of spike architecture in wheat (Triticum aestivum L.). A major and consistent quantitative trait locus (QTL) was discovered for TSN in a doubled haploid spring wheat population grown in the field over 4 years. The QTL on chromosome 7B explained up to 20.5% of phenotypic variance. In its physical interval (7B: 6.37-21.67 Mb), the gene FLOWERING LOCUS T (FT-B1) emerged as candidate for the observed effect. In one of the parental lines, FT-B1 carried a non-synonymous substitution on position 19 of the coding sequence. This mutation modifying an aspartic acid (D) into a histidine (H) occurred in a highly conserved position. The mutation was observed with a frequency of ca. 68% in a set of 135 hexaploid wheat varieties and landraces, while it was not found in other plant species. FT-B1 only showed a minor effect on heading and flowering time (FT) which were dominated by a major QTL on chromosome 5A caused by segregation of the vernalization gene VRN-A1. Individuals carrying the FT-B1 allele with amino acid histidine had, on average, a higher number of spikelets (15.1) than individuals with the aspartic acid allele (14.3) independent of their VRN-A1 allele. We show that the effect of TSN is not mainly related to flowering time; however, the duration of pre-anthesis phases may play a major role.

DEFECTIVE ENDOSPERM-D1 (Dee-D1) is crucial for endosperm development in hexaploid wheat.

Hexaploid wheat (Triticum aestivum L.) is a natural allopolyploid and provides a usable model system to better understand the genetic mechanisms that underlie allopolyploid speciation through the hybrid genome doubling. Here we aimed to identify the contribution of chromosome 1D in the development and evolution of hexaploid wheat. We identified and mapped a novel DEFECTIVE ENDOSPERM-D1 (Dee-D1) locus on 1DL that is involved in the genetic control of endosperm development. The absence of Dee-D1 leads to non-viable grains in distant crosses and alters grain shape, which negatively affects grain number and thousand-grain weight. Dee-D1 can be classified as speciation locus with a positive effect on the function of genes which are involved in endosperm development in hybrid genomes. The presence of Dee-D1 is necessary for the normal development of endosperm, and thus play an important role in the evolution and improvement of grain yield in hexaploid wheat.

Comparison Between Core Set Selection Methods Using Different Illumina Marker Platforms: A Case Study of Assessment of Diversity in Wheat.

Collections of plant genetic resources stored in genebanks are an important source of genetic diversity for improvement in plant breeding programs and for conservation of natural variation. The establishment of reduced representative collections from a large set of genotypes is a valuable tool that provides cost-effective access to the diversity present in the whole set. Software like Core Hunter 3 is available to generate high quality core sets. In addition, general clustering approaches, e.g., k-medoids, are available to subdivide a large data set into small groups with maximum genetic diversity between groups. Illumina genotyping platforms are a very efficient tool for the assessment of genetic diversity of plant genetic resources. The accumulation of genotyping data over time using commercial genotyping platforms raises the question of how such huge amount of information can be efficiently used for creating core collections. In the present study, after developing a 15K wheat Infinium array with 12,908 SNPs and genotyping a set of 479 hexaploid winter wheat lines (Triticum aestivum), a larger data set was created by merging 411 lines previously genotyped with the 90K iSelect array. Overlaying the markers from the 15K and 90K arrays enabled the identification of a common set of 12,806 markers, suggesting that the 15K array is a valuable and cost-effective resource for plant breeding programs. Finally, we selected genetically diverse core sets out of these 890 wheat genotypes derived from five collections based on the common markers from the 15K and 90K SNP arrays. Two different approaches, k-medoids and Core Hunter 3 were compared,and k-medoids was identified as an efficient method for selecting small core sets out of a large collection of genotypes while retaining the genetic diversity of the original population.

Prospects of GWAS and predictive breeding for European winter wheat's grain protein content, grain starch content, and grain hardness.

Grain quality traits determine the classification of registered wheat (Triticum aestivum L.) varieties. Although environmental factors and crop management practices exert a considerable influence on wheat quality traits, a significant proportion of the variance is attributed to the genetic factors. To identify the underlying genetic factors of wheat quality parameters viz., grain protein content (GPC), grain starch content (GSC), and grain hardness (GH), we evaluated 372 diverse European wheat varieties in replicated field trials in up to eight environments. We observed that all of the investigated traits hold a wide and significant genetic variation, and a significant negative correlation exists between GPC and GSC plus grain yield. Our association analyses based on 26,694 high-quality single nucleotide polymorphic markers revealed a strong quantitative genetic nature of GPC and GSC with associations on groups 2, 3, and 6 chromosomes. The identification of known Puroindoline-b gene for GH provided a positive analytic proof for our studies. We report that a locus QGpc.ipk-6A controls both GPC and GSC with opposite allelic effects. Based on wheat's reference and pan-genome sequences, the physical characterization of two loci viz., QGpc.ipk-2B and QGpc.ipk-6A facilitated the identification of the candidate genes for GPC. Furthermore, by exploiting additive and epistatic interactions of loci, we evaluated the prospects of predictive breeding for the investigated traits that suggested its efficient use in the breeding programs.

Identification and validation of QTL for spike fertile floret and fruiting efficiencies in hexaploid wheat (Triticum aestivum L.).

KEY MESSAGE: This study identified and validated two QTL associated with spike fertile floret and fruiting efficiencies. They represent two new loci to use in MAS to improve wheat yield potential. The spike fruiting efficiency (FE-grains per unit spike dry weight at anthesis, GN/SDW) is a promising trait to improve wheat yield potential. It depends on fertile floret efficiency (fertile florets per unit SDW-FFE, FF/SDW) and grain set (grains per fertile floret-GST). Given its difficult measurement, it is often estimated as the grains per unit of nongrain spike dry weight at maturity (FEm). In this study, quantitative trait loci (QTL) were mapped using a double haploid population (Baguette 19/BIOINTA 2002, with high and low FE, respectively) genotyped with the iSelect 90 K SNP array and evaluated in five environments. We identified 37 QTL, but two were major with an R2 > 10% and stable for being at least present in three environments: the QFEm.perg-3A (on Chr. 3A, 51.6 cM, 685.12 Mb) for FEm and the QFFE.perg-5A (on Chr. 5A, 42.1 cM, 461.49 Mb) for FFE, FE and FEm. Both QTL were validated using two independent F2 populations and KASP markers. For the most promising QTL, QFFE.perg-5A, the presence of the allele for high FFE resulted in + 4% FF, + 9% GN, + 13% GST, + 16% yield gSDW-1 and + 5% yield spike-1. QFEm.perg-3A and QFFE.perg-5A represent two new loci to use in MAS to improve wheat yield potential.

Negative dominance and dominance-by-dominance epistatic effects reduce grain-yield heterosis in wide crosses in wheat.

The genetics underlying heterosis, the difference in performance of crosses compared with midparents, is hypothesized to vary with relatedness between parents. We established a unique germplasm comprising three hybrid wheat sets differing in the degree of divergence between parents and devised a genetic distance measure giving weight to heterotic loci. Heterosis increased steadily with heterotic genetic distance for all 1903 hybrids. Midparent heterosis, however, was significantly lower in the hybrids including crosses between elite and exotic lines than in crosses among elite lines. The analysis of the genetic architecture of heterosis revealed this to be caused by a higher portion of negative dominance and dominance-by-dominance epistatic effects. Collectively, these results expand our understanding of heterosis in crops, an important pillar toward global food security.

The potential of hybrid breeding to enhance leaf rust and stripe rust resistance in wheat.

KEY MESSAGE: Hybrid wheat breeding is a promising strategy to improve the level of leaf rust and stripe rust resistance in wheat. Leaf rust and stripe rust belong to the most important fungal diseases in wheat production. Due to a dynamic development of new virulent races, epidemics appear in high frequency and causes significant losses in grain yield and quality. Therefore, research is needed to develop strategies to breed wheat varieties carrying highly efficient resistances. Stacking of dominant resistance genes through hybrid breeding is such an approach. Within this study, we investigated the genetic architecture of leaf rust and stripe rust resistance of 1750 wheat hybrids and their 230 parental lines using a genome-wide association study. We observed on average a lower rust susceptibility for hybrids in comparison to their parental inbred lines and some hybrids outperformed their better parent with up to 56%. Marker-trait associations were identified on chromosome 3D and 4A for leaf rust and on chromosome 2A, 2B, and 6A for stripe rust resistance by using a genome-wide association study with a Bonferroni-corrected threshold of P < 0.10. Detected loci on chromosomes 4A and 2A were located within previously reported genomic regions affecting leaf rust and stripe rust resistance, respectively. The degree of dominance was for most associations favorable in the direction of improved resistance. Thus, resistance can be increased in hybrid wheat breeding by fixing complementary leaf rust and stripe rust resistance genes with desired dominance effects in opposite parental pools.

Exome association analysis sheds light onto leaf rust (Puccinia triticina) resistance genes currently used in wheat breeding (Triticum aestivum L.).

Resistance breeding is crucial for a sustainable control of leaf rust (Puccinia triticina) in wheat (Triticum aestivum L.) while directly targeting functional variants is the Holy Grail for efficient marker-assisted selection and map-based cloning. We assessed the limits and prospects of exome association analysis for severity of leaf rust in a large hybrid wheat population of 1574 single-crosses plus their 133 parents. After imputation and quality control, exome sequencing revealed 202 875 single-nucleotide polymorphisms (SNPs) covering 19.7% of the high-confidence annotated gene space. We performed intensive data mining and found significant associations for 2171 SNPs corresponding to 50 different loci. Some of these associations mapped in the proximity of the already known resistance genes Lr21, Lr34-B, Lr1 and Lr10, while other associated genomic regions, such as those on chromosomes 1A and 3D, harboured several annotated genes putatively involved in resistance. Validation with an independent population helped to narrow down the list of putative resistance genes that should be targeted by fine-mapping. We expect that the proposed strategy of intensive data mining coupled with validation will significantly influence research in plant genetics and breeding.

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.

We dissected the genetic basis of total spikelet number (TSN) along with other traits, viz. spike length (SL) and flowering time (FT) in a panel of 518 elite European winter wheat varieties. Genome-wide association studies (GWAS) based on 39,908 SNP markers revealed highly significant quantitative trait loci (QTL) for TSN on chromosomes 2D, 7A, and 7B, for SL on 5A, and FT on 2D, with 2D-QTL being the functional marker for the gene Ppd-D1. The physical region of the 7A-QTL for TSN revealed the presence of a wheat ortholog (TaAPO-A1) to APO1-a rice gene that positively controls the spikelet number on the panicles. Interspecific analyses of the TaAPO-A1 orthologs showed that it is a highly conserved gene important for floral development and present in a wide range of terrestrial plants. Intraspecific studies of the TaAPO-A1 across wheat genotypes revealed a polymorphism in the conserved F-box domain, defining two haplotypes. A KASP marker developed on the polymorphic site showed a highly significant association of TaAPO-A1 with TSN, explaining 23.2% of the total genotypic variance. Also, the TaAPO-A1 alleles showed weak but significant differences for SL and grain yield. Our results demonstrate the importance of wheat sequence resources to identify candidate genes for important traits based on genetic analyses.

Genome-wide Association Mapping and Prediction of Adult Stage Septoria tritici Blotch Infection in European Winter Wheat via High-Density Marker Arrays.

blotch (STB) caused by the fungus is a devastating foliar disease of wheat ( L.) that can lead to substantial yield losses. Quantitative genetic resistance has been proposed as a durable strategy for STB control. In this study, we dissected the genetic basis of STB infection in 371 European wheat varieties based on 35k and 90k single nucleotide polymorphism marker arrays. The phenotypic data analyses suggested that large genetic variance exists for STB infection with a broad-sense heritability of 0.78. Genome-wide association studies (GWAS) propose the highly quantitative nature of STB infection with potential associations on chromosomes 1A, 1B, 2D, 4A, 5A, 6A, 6D, 7A, and 7B. Increased marker density in GWAS by combining markers from both arrays helped to detect additional markers explaining increased genotypic variance. Linkage disequilibrium analyses revealed genes with a possible role in disease resistance. The potential of genomic prediction (GP) assessed via two models accounting for additive effects and additive plus epistatic interactions among the loci suggested the possibility of genomic selection for improved STB resistance. Genomic prediction results also indicated that the higher-order epistatic interactions are not abundant and that both marker platforms are equally suitable for GP of STB infection. Our results provide further understanding of the quantitative genetic nature of STB infection, serve as a resource for marker-assisted breeding, and highlight the potential of genomic selection for improved STB resistance.

Whole-Genome Association Mapping and Genomic Prediction for Iron Concentration in Wheat Grains.

Malnutrition of iron (Fe) affects two billion people worldwide. Therefore, enhancing grain Fe concentration (GFeC) in wheat (Triticum aestivum L.) is an important goal for breeding. Here we study the genetic factors underlying GFeC trait by genome-wide association studies (GWAS) and the prediction abilities using genomic prediction (GP) in a panel of 369 European elite wheat varieties which was genotyped with 15,523 mapped single-nucleotide polymorphism markers (SNP) and a subpanel of 183 genotypes with 44,233 SNP markers. The resulting means of GFeC from three field experiments ranged from 24.42 to 52.42 µg·g-1 with a broad-sense heritability (H²) equaling 0.59 over the years. GWAS revealed 41 and 137 significant SNPs in the whole and subpanel, respectively, including significant marker-trait associations (MTAs) for best linear unbiased estimates (BLUEs) of GFeC on chromosomes 2A, 3B and 5A. Putative candidate genes such as NAC transcription factors and transmembrane proteins were present on chromosome 2A (763,689,738⁻765,710,113 bp). The GP for a GFeC trait ranged from low to moderate values. The current study reported GWAS of GFeC for the first time in hexaploid wheat varieties. These findings confirm the utility of GWAS and GP to explore the genetic architecture of GFeC for breeding programs aiming at the improvement of wheat grain quality.

Identifying Candidate Genes for Enhancing Grain Zn Concentration in Wheat.

Wheat (Triticum aestivum L.) is one of the major staple food crops worldwide. Despite efforts in improving wheat quality, micronutrient levels are still below the optimal range for human nutrition. In particular, zinc (Zn) deficiency is a widespread problem in human nutrition in countries relying mainly on a cereal diet; hence improving Zn accumulation in grains is an imperative need. This study was designed to understand the genetic architecture of Zn grain concentrations in wheat grains. We performed a genome-wide association study (GWAS) for grain Zn concentrations in 369 European wheat genotypes, using field data from 3 years. The complete wheat panel was genotyped by high-density arrays of single nucleotide polymorphic (SNP) markers (90k iSELECT Infinium and 35k Affymetrix arrays) resulting in 15,523 polymorphic markers. Additionally, a subpanel of 183 genotypes was analyzed with a novel 135k Affymetrix marker array resulting in 28,710 polymorphic SNPs for high-resolution mapping of the potential genomic regions. The mean grain Zn concentration of the genotypes ranged from 25.05-52.67 µg g-1 dry weight across years with a moderate heritability value. Notably, 40 marker-trait associations (MTAs) were detected in the complete panel of varieties on chromosomes 2A, 3A, 3B, 4A, 4D, 5A, 5B, 5D, 6D, 7A, 7B, and 7D. The number of MTAs in the subpanel was increased to 161 MTAs whereas the most significant and consistent associations were located on chromosomes 3B (723,504,241-723,611,488 bp) and 5A (462,763,758-466,582,184 bp) having major effects. These genomic regions include newly identified putative candidate genes, which are related to Zn uptake and transport or represent bZIP and mitogen-activated protein kinase genes. These findings provide the basis for understanding the genetic background of Zn concentration in wheat grains that in turn may help breeders to select high Zn-containing genotypes to improve human health and grain quality.

Manipulation and prediction of spike morphology traits for the improvement of grain yield in wheat.

In wheat (Triticum spp.), modifying inflorescence (spike) morphology can increase grain number and size and thus improve yield. Here, we demonstrated the potential for manipulating and predicting spike morphology, based on 44 traits. In 12 wheat cultivars, we observed that detillering (removal of branches), which alters photosynthate distribution, changed spike morphology. Our genome-wide association study detected close associations between carbon partitioning (e.g. tiller number, main shoot dry weight) and spike morphology (e.g. spike length, spikelet density) traits in 210 cultivars. Most carbon-partitioning traits (e.g. tiller dry weight, harvest index) demonstrated high prediction abilities (>0.5). For spike morphology, some traits (e.g. total and fertile spikelet number, spike length) displayed high prediction abilities (0.3-0.5), but others (e.g. spikelet fertility, spikelet density) exhibited low prediction abilities (<0.2). Grain size traits were closely correlated in field and greenhouse experiments. Stepwise regression analysis suggests that significantly associated traits in the greenhouse explain 35.35% of the variation in grain yield and 67.63% of the variation in thousand-kernel weight in the field. Therefore, the traits identified in this study affect spike morphology; these traits can be used to predict and improve plant architecture and thus increase yield.

Genetic dissection of pre-anthesis sub-phase durations during the reproductive spike development of wheat.

Flowering time is an important factor affecting grain yield in wheat. In this study, we divided reproductive spike development into eight sub-phases. These sub-phases have the potential to be delicately manipulated to increase grain yield. We measured 36 traits with regard to sub-phase durations, determined three grain yield-related traits in eight field environments and mapped 15 696 single nucleotide polymorphism (SNP, based on 90k Infinium chip and 35k Affymetrix chip) markers in 210 wheat genotypes. Phenotypic and genetic associations between grain yield traits and sub-phase durations showed significant consistency (Mantel test; r = 0.5377, P < 0.001). The shared quantitative trait loci (QTLs) revealed by the genome-wide association study suggested a close association between grain yield and sub-phase duration, which may be attributed to effects on spikelet initiation/spikelet number (double ridge to terminal spikelet stage, DR-TS) and assimilate accumulation (green anther to anthesis stage, GA-AN). Moreover, we observed that the photoperiod-sensitivity allele at the Ppd-D1 locus on chromosome 2D markedly extended all sub-phase durations, which may contribute to its positive effects on grain yield traits. The dwarfing allele at the Rht-D1 (chromosome 4D) locus altered the sub-phase duration and displayed positive effects on grain yield traits. Data for 30 selected genotypes (from among the original 210 genotypes) in the field displayed a close association with that from the greenhouse. Most importantly, this study demonstrated specific connections to grain yield in narrower time windows (i.e. the eight sub-phases), rather than the entire stem elongation phase as a whole.

Genome-wide association analyses of plant growth traits during the stem elongation phase in wheat.

One of the primary objectives of wheat breeding was to increase grain yield. Floral abortion during the stem elongation phase (SEP) leads to a loss of more than 50% of the grain number potential. In this study, we quantified 75 plant growth-associated traits at seven stages during the SEP and mapped 15 696 single nucleotide polymorphism (SNP) markers in 210 accessions of wheat (Triticum aestivum). Our genomewide association study identified trait-associated SNPs that are shared among various stages of the SEP, as well as SNPs that are shared between plant growth traits and grain yield in the field. The genomic selection analysis shows variation among the prediction abilities of various traits and stages. Furthermore, we found that the allelic variants of Ppd-D1 (chromosome 2D) and Rht-D1 (chromosome 4D) loci affect some plant growth traits (e.g. leaf area and spike length). These results have identified a narrow time window within the SEP in which plant growth traits can be manipulated to alter grain yield. This suggests that there may be multiple ways to regulate plant growth during the SEP, to ultimately influence grain number in wheat.