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.

Dissecting Bread Wheat Heterosis through the Integration of Agronomic and Physiological Traits.

To meet the challenge of feeding almost 10 billion people by 2050, wheat yield has to double by 2050. However, over the past 20 years, yield increase has slowed down and even stagnated in the main producing countries. Following the example of maize, hybrids have been suggested as a solution to overcome yield stagnation in wheat. However, wheat heterosis is still limited and poorly understood. Gaining a better understanding of hybrid vigor holds the key to breed for better varieties. To this aim, we have developed and phenotyped for physiological and agronomic traits an incomplete factorial design consisting of 91 hybrids and their nineteen female and sixteen male parents. Monitoring the plant development with normalized difference vegetation index revealed that 89% of the hybrids including the five higher yielding hybrids had a longer grain filling phase with a delayed senescence that results in larger grain size. This average increase of 7.7% in thousand kernel weight translated to a positive mid-parent heterosis for grain yield for 86% of hybrids. In addition, hybrids displayed a positive grain protein deviation leading to a +4.7% heterosis in protein yield. These results shed light on the physiological bases underlying yield heterosis in wheat, paving new ways to breed for better wheat hybrids.

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.

The roles of pleiotropy and close linkage as revealed by association mapping of yield and correlated traits of wheat (Triticum aestivum L.).

Grain yield (GY) of bread wheat (Triticum aestivum L.) is quantitatively inherited. Correlated GY-syndrome traits such as plant height (PH), heading date (HD), thousand grain weight (TGW), test weight (TW), grains per ear (GPE), and ear weight (EW) influence GY. Most quantitative genetics studies assessed the multiple-trait (MT) complex of GY-syndrome using single-trait approaches, and little is known about its underlying pleiotropic architecture. We investigated the pleiotropic architecture of wheat GY-syndrome through MT association mapping (MT-GWAS) using 372 varieties phenotyped in up to eight environments and genotyped with 18 832 single nucleotide polymorphisms plus 24 polymorphic functional markers. MT-GWAS revealed a total of 345 significant markers spread genome wide, representing 8, 40, 11, 40, 34, and 35 effective GY-PH, GY-HD, GY-TGW, GY-TW, GY-GPE, and GY-EW associations, respectively. Among them, pleiotropic roles of Rht-B1 and TaGW2-6B loci were corroborated. Only one marker presented simultaneous associations for three traits (i.e. GY-TGW-TW). Close linkage was difficult to differentiate from pleiotropy; thus, the pleiotropic architecture of GY-syndrome was dissected more as a cause of pleiotropy rather than close linkage. Simulations showed that minor allele frequencies, along with sizes and distances between quantitative trait loci for two traits, influenced the ability to distinguish close linkage from pleiotropy.

Genome-wide association mapping of resistance to eyespot disease (Pseudocercosporella herpotrichoides) in European winter wheat (Triticum aestivum L.) and fine-mapping of Pch1.

KEY MESSAGE: Genotypes with recombination events in the Triticum ventricosum introgression on chromosome 7D allowed to fine-map resistance gene Pch1, the main source of eyespot resistance in European winter wheat cultivars. Eyespot (also called Strawbreaker) is a common and serious fungal disease of winter wheat caused by the necrotrophic fungi Oculimacula yallundae and Oculimacula acuformis (former name Pseudocercosporella herpotrichoides). A genome-wide association study (GWAS) for eyespot was performed with 732 microsatellite markers (SSR) and 7761 mapped SNP markers derived from the 90 K iSELECT wheat array using a panel of 168 European winter wheat varieties as well as three spring wheat varieties and phenotypic evaluation of eyespot in field tests in three environments. Best linear unbiased estimations (BLUEs) were calculated across all trials and ranged from 1.20 (most resistant) to 5.73 (most susceptible) with an average value of 4.24 and a heritability of H 2 = 0.91. A total of 108 SSR and 235 SNP marker-trait associations (MTAs) were identified by considering associations with a -log10 (P value) ≥3.0. Significant MTAs for eyespot-score BLUEs were found on chromosomes 1D, 2A, 2D, 3D, 5A, 5D, 6A, 7A and 7D for the SSR markers and chromosomes 1B, 2A, 2B, 2D, 3B and 7D for the SNP markers. For 18 varieties (10.5%), a highly resistant phenotype was detected that was linked to the presence of the resistance gene Pch1 on chromosome 7D. The identification of genotypes with recombination events in the introgressed genomic segment from Triticum ventricosum harboring the Pch1 resistance gene on chromosome 7DL allowed the fine-mapping of this gene using additional SNP markers and a potential candidate gene Traes_7DL_973A33763 coding for a CC-NBS-LRR class protein was identified.

Validating the prediction accuracies of marker-assisted and genomic selection of Fusarium head blight resistance in wheat using an independent sample.

KEY MESSAGE: Compared with independent validation, cross-validation simultaneously sampling genotypes and environments provided similar estimates of accuracy for genomic selection, but inflated estimates for marker-assisted selection. Estimates of prediction accuracy of marker-assisted (MAS) and genomic selection (GS) require validations. The main goal of our study was to compare the prediction accuracies of MAS and GS validated in an independent sample with results obtained from fivefold cross-validation using genomic and phenotypic data for Fusarium head blight resistance in wheat. In addition, the applicability of the reliability criterion, a concept originally developed in the context of classic animal breeding and GS, was explored for MAS. We observed that prediction accuracies of MAS were overestimated by 127% using cross-validation sampling genotype and environments in contrast to independent validation. In contrast, prediction accuracies of GS determined in independent samples are similar to those estimated with cross-validation sampling genotype and environments. This can be explained by small population differentiation between the training and validation sets in our study. For European wheat breeding, which is so far characterized by a slow temporal dynamic in allele frequencies, this assumption seems to be realistic. Thus, GS models used to improve European wheat populations are expected to possess a long-lasting validity. Since quantitative trait loci information can be exploited more precisely if the predicted genotype is more related to the training population, the reliability criterion is also a valuable tool to judge the level of prediction accuracy of individual genotypes in MAS.

Exploiting the Repetitive Fraction of the Wheat Genome for High-Throughput Single-Nucleotide Polymorphism Discovery and Genotyping.

Transposable elements (TEs) account for more than 80% of the wheat genome. Although they represent a major obstacle for genomic studies, TEs are also a source of polymorphism and consequently of molecular markers such as insertion site-based polymorphism (ISBP) markers. Insertion site-based polymorphisms have been found to be a great source of genome-specific single-nucleotide polymorphism (SNPs) in the hexaploid wheat ( L.) genome. Here, we report on the development of a high-throughput SNP discovery approach based on sequence capture of ISBP markers. By applying this approach to the reference sequence of chromosome 3B from hexaploid wheat, we designed 39,077 SNPs that are evenly distributed along the chromosome. We demonstrate that these SNPs can be efficiently scored with the KASPar (Kompetitive allele-specific polymerase chain reaction) genotyping technology. Finally, through genetic diversity and genome-wide association studies, we also demonstrate that ISBP-derived SNPs can be used in marker-assisted breeding programs.

Analysis of main effect QTL for thousand grain weight in European winter wheat (Triticum aestivum L.) by genome-wide association mapping.

Grain weight, an essential yield component, is under strong genetic control and at the same time markedly influenced by the environment. Genetic analysis of the thousand grain weight (TGW) by genome-wide association study (GWAS) was performed with a panel of 358 European winter wheat (Triticum aestivum L.) varieties and 14 spring wheat varieties using phenotypic data of field tests in eight environments. Wide phenotypic variations were indicated for the TGW with BLUEs (best linear unbiased estimations) values ranging from 35.9 to 58.2 g with a mean value of 45.4 g and a heritability of H(2) = 0.89. A total of 12 candidate genes for plant height, photoperiodism and grain weight were genotyped on all varieties. Only three candidates, the photoperiodism gene Ppd-D1, dwarfing gene Rht-B1and the TaGW-6A gene were significant explaining up to 14.4, 2.3, and 3.4% of phenotypic variation, respectively. For a comprehensive genome-wide analysis of TGW-QTL genotyping data from 732 microsatellite markers and a set of 7769 mapped SNP-markers genotyped with the 90k iSELECT array were analyzed. In total, 342 significant (-log10 (P-value) ≥ 3.0) marker trait associations (MTAs) were detected for SSR-markers and 1195 MTAs (-log10(P-value) ≥ 3.0) for SNP-markers in all single environments plus the BLUEs. After Bonferroni correction, 28 MTAs remained significant for SSR-markers (-log10 (P-value) ≥ 4.82) and 58 MTAs for SNP-markers (-log10 (P-value) ≥ 5.89). Apart from chromosomes 4B and 6B for SSR-markers and chromosomes 4D and 5D for SNP-markers, MTAs were detected on all chromosomes. The highest number of significant SNP-markers was found on chromosomes 3B and 1B, while for the SSRs most markers were significant on chromosomes 6D and 3D. Overall, TGW was determined by many markers with small effects. Only three SNP-markers had R(2) values above 6%.

Whole genome association mapping of plant height in winter wheat (Triticum aestivum L.).

The genetic architecture of plant height was investigated in a set of 358 recent European winter wheat varieties plus 14 spring wheat varieties based on field data in eight environments. Genotyping of diagnostic markers revealed the Rht-D1b mutant allele in 58% of the investigated varieties, while the Rht-B1b mutant was only present in 7% of the varieties. Rht-D1 was significantly associated with plant height by using a mixed linear model and employing a kinship matrix to correct for population stratification. Further genotyping data included 732 microsatellite markers, resulting in 770 loci, of which 635 markers were placed on the ITMI map plus a set of 7769 mapped SNP markers genotyped with the 90 k iSELECT chip. When Bonferroni correction was applied, a total of 153 significant marker-trait associations (MTAs) were observed for plant height and the SSR markers (-log10 (P-value) ≥ 4.82) and 280 (-log10 (P-value) ≥ 5.89) for the SNPs. Linear regression between the most effective markers and the BLUEs for plant height indicated additive effects for the MTAs of different chromosomal regions. Analysis of syntenic regions in the rice genome revealed closely linked rice genes related to gibberellin acid (GA) metabolism and perception, i.e. GA20 and GA2 oxidases orthologous to wheat chromosomes 1A, 2A, 3A, 3B, 5B, 5D and 7B, ent-kaurenoic acid oxidase orthologous to wheat chromosome 7A, ent-kaurene synthase on wheat chromosome 2B, as well as GA-receptors like DELLA genes orthologous to wheat chromosomes 4B, 4D and 7A and genes of the GID family orthologous to chromosomes 2B and 5B. The data indicated that besides the widely used GA-insensitive dwarfing genes Rht-B1 and Rht-D1 there is a wide spectrum of loci available that could be used for modulating plant height in variety development.

Genetic architecture of main effect QTL for heading date in European winter wheat.

A genome-wide association study (GWAS) for heading date (HD) was performed with a panel of 358 European winter wheat (Triticum aestivum L.) varieties and 14 spring wheat varieties through the phenotypic evaluation of HD in field tests in eight environments. Genotyping data consisted of 770 mapped microsatellite loci and 7934 mapped SNP markers derived from the 90K iSelect wheat chip. Best linear unbiased estimations (BLUEs) were calculated across all trials and ranged from 142.5 to 159.6 days after the 1st of January with an average value of 151.4 days. Considering only associations with a -log10 (P-value) ≥ 3.0, a total of 340 SSR and 2983 SNP marker-trait associations (MTAs) were detected. After Bonferroni correction for multiple testing, a total of 72 SSR and 438 SNP marker-trait associations remained significant. Highly significant MTAs were detected for the photoperiodism gene Ppd-D1, which was genotyped in all varieties. Consistent associations were found on all chromosomes with the highest number of MTAs on chromosome 5B. Linear regression showed a clear dependence of the HD score BLUEs on the number of favorable alleles (decreasing HD) and unfavorable alleles (increasing HD) per variety meaning that genotypes with a higher number of favorable or a low number of unfavorable alleles showed lower HD and therefore flowered earlier. For the vernalization gene Vrn-A2 co-locating MTAs on chromosome 5A, as well as for the photoperiodism genes Ppd-A1 and Ppd-B1 on chromosomes 2A and 2B were detected. After the construction of an integrated map of the SSR and SNP markers and by exploiting the synteny to sequenced species, such as rice and Brachypodium distachyon, we were able to demonstrate that a marker locus on wheat chromosome 5BL with homology to the rice photoperiodism gene Hd6 played a significant role in the determination of the heading date in wheat.

Whole genome association mapping of Fusarium head blight resistance in European winter wheat (Triticum aestivum L.).

A total of 358 recent European winter wheat varieties plus 14 spring wheat varieties were evaluated for resistance to Fusarium head blight (FHB) caused by Fusarium graminearum and Fusarium culmorum in four separate environments. The FHB scores based on FHB incidence (Type I resistance)×FHB severity (Type II resistance) indicated a wide phenotypic variation of the varieties with BLUE (best linear unbiased estimation) values ranging from 0.07 to 33.67. Genotyping with 732 microsatellite markers resulted in 782 loci of which 620 were placed on the ITMI map. The resulting average marker distance of 6.8 cM allowed genome wide association mapping employing a mixed model. Though no clear population structure was discovered, a kinship matrix was used for stratification. A total of 794 significant (-log10(p)-value≥3.0) associations between SSR-loci and environment-specific FHB scores or BLUE values were detected, which included 323 SSR alleles. For FHB incidence and FHB severity a total of 861 and 877 individual marker-trait associations (MTA) were detected, respectively. Associations for both traits co-located with FHB score in most cases. Consistent associations detected in three or more environments were found on all chromosomes except chromosome 6B, and with the highest number of MTA on chromosome 5B. The dependence of the number of favourable and unfavourable alleles within a variety to the respective FHB scores indicated an additive effect of favourable and unfavourable alleles, i.e. genotypes with more favourable or less unfavourable alleles tended to show greater resistance to FHB. Assessment of a marker specific for the dwarfing gene Rht-D1 resulted in strong effects. The results provide a prerequisite for designing genome wide breeding strategies for FHB resistance.

In search of a maize ideotype for cell wall enzymatic degradability using histological and biochemical lignin characterization.

Grass cell wall degradability is conventionally related to the lignin content and to the ferulic-mediated cross-linking of lignins to polysaccharides. To better understand the variations in degradability, 22 maize inbred lines were subjected to image analyses of Fasga- and Mäule-stained stem sections and to chemical analyses of lignins and p-hydroxycinnamic acids. For the first time, the nearness of biochemical and histological estimates of lignin levels was established. Combination of histological and biochemical traits could explain 89% of the variations for cell wall degradability and define a maize ideotype for cell wall degradability. In addition to a reduced lignin level, such an ideotype would contain lignins richer in syringyl than in guaiacyl units and preferentially localized in the cortical region rather than in the pith. Such enrichment in syringyl units would favor wall degradability in grasses, contrary to dicots, and could be related to the fact that grass syringyl units are noticeably p-coumaroylated. This might affect the interaction capabilities of lignins and polysaccharides.

Genetic and molecular basis of grass cell wall biosynthesis and degradability. II. Lessons from brown-midrib mutants.

The brown-midrib mutants of maize have a reddish-brown pigmentation of the leaf midrib and stalk pith, associated with lignified tissues. These mutants progressively became models for lignification genetics and biochemical studies in maize and grasses. Comparisons at silage maturity of bm1, bm2, bm3, bm4 plants highlighted their reduced lignin, but also illustrated the biochemical specificities of each mutant in p-coumarate, ferulate ester and etherified ferulate content, or syringyl/guaiacyl monomer ratio after thioacidolysis. Based on the current knowledge of the lignin pathway, and based on presently developed data and discussions, C3H and CCoAOMT activities are probably major hubs in controlling cell-wall lignification (and digestibility). It is also likely that ferulates arise via the CCoAOMT pathway.