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.

Wheat Grains as a Sustainable Source of Protein for Health

Protein deficiency is recognized among the major global health issues with an underestimation of its importance. Genetic biofortification is a cost-effective and sustainable strategy to overcome global protein malnutrition. This study was designed to focus on protein-dense grains of wheat (Triticum aestivum L.) and identify the genes governing grain protein content (GPC) that improve end-use quality and in turn human health. Genome-wide association was applied using the 90k iSELECT Infinium and 35k Affymetrix arrays with GPC quantified by using a proteomic-based technique in 369 wheat genotypes over three field-year trials. The results showed significant natural variation among bread wheat genotypes that led to detecting 54 significant quantitative trait nucleotides (QTNs) surpassing the false discovery rate (FDR) threshold. These QTNs showed contrasting effects on GPC ranging from -0.50 to +0.54% that can be used for protein content improvement. Further bioinformatics analyses reported that these QTNs are genomically linked with 35 candidate genes showing high expression during grain development. The putative candidate genes have functions in the binding, remobilization, or transport of protein. For instance, the promising QTN AX-94727470 on chromosome 6B increases GPC by +0.47% and is physically located inside the gene TraesCS6B02G384500 annotated as Trehalose 6-phosphate phosphatase (T6P), which can be employed to improve grain protein quality. Our findings are valuable for the enhancement of protein content and end-use quality in one of the major daily food resources that ultimately improve human nutrition.

Identification of genomic regions associated with cereal cyst nematode (Heterodera avenae Woll.) resistance in spring and winter wheat

Cereal cyst nematode (CCN) is a major threat to cereal crop production globally including wheat (Triticum aestivum L.). In the present study, single-locus and multi-locus models of Genome-Wide Association Study (GWAS) were used to find marker trait associations (MTAs) against CCN (Heterodera avenae) in wheat. In total, 180 wheat accessions (100 spring and 80 winter types) were screened against H. avenae in two independent years (2018/2019 "Environment 1" and 2019/2020 "Environment 2") under controlled conditions. A set of 12,908 SNP markers were used to perform the GWAS. Altogether, 11 significant MTAs, with threshold value of -log10 (p-values) ≥ 3.0, were detected using 180 wheat accessions under combined environment (CE). A novel MTA (wsnp_Ex_c53387_56641291) was detected under all environments (E1, E2 and CE) and considered to be stable MTA. Among the identified 11 MTAs, eight were novel and three were co-localized with previously known genes/QTLs/MTAs. In total, 13 putative candidate genes showing differential expression in roots, and known to be involved in plant defense mechanisms were reported. These MTAs could help us to identify resistance alleles from new sources, which could be used to identify wheat varieties with enhanced CCN resistance.

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.

Measurement(s)	agronomic, quality and disease traits
Technology Type(s)	manual measurement in the field
Sample Characteristic - Organism	Triticum aestivum L.

Evaluation of the genetic structure of Bromus inermis populations from chemically and radioactively polluted areas using microsatellite markers from closely related species

Hypothesis The ecotoxicological and radiobiological effects can be manifested in a decrease in genetic diversity with an increase in toxic and radiation load, in an increase in the frequencies of rare and/or unique (private) alleles in impact samples, and in a decrease in the differentiation of B. inermis populations within each pollution area.Materials and methods We have selected a collection of primers for Bromus inermis, consisting of 21 microsatellite (SSR) loci from B. sterilis, B. tectorum and Triticum aestivum. The level of toxic load (chemically polluted area) was 4-19 conventional units, and the absorbed dose rate (the Kyshtym accident area) varied from 0.153 to 21.5 μGy h^-1, which is up to two orders higher than the natural background radiation level (≈ 0.1 μGy h^-1).Results Only eight of 21 (38%) of SSR primers showed good transferability and were used for B. inermis population studies from areas of technogenic pollution (heavy metals and radionuclides). We revealed 42 alleles at eight loci, and the number of alleles per locus varied from one to 13 in B. inermis populations. The percentage of polymorphic loci in B. inermis populations was 48.44%, the polymorphism information content (PIC) value was 0.556, and Shannon information index was 0.69 ± 0.3. A total of 22 rare, 14 private and 9 both rare and private alleles were reported for all B. inermis populations. There were no correlations between geographic and genetic distances. Only 6.8% of the genetic variability was distributed among B. inermis populations.Conclusion There was no decrease in genetic diversity ("genetic erosion") found in B. inermis populations growing for a long time under anthropogenic stress. No significant differences in the number of rare and private alleles in the background and impact populations of B. inermis were found. The smooth brome is characterized by low differentiation of the populations. Possible reasons for this phenomenon are discussed.

A mechanistic view on lodging resistance in rye and wheat: a multiscale comparative study

The development of crop varieties that are resistant to lodging is a top priority for breeding programmes. Herein, we characterize the rye mutant ´Stabilstroh' ('stable straw') possessing an exceptional combination of high lodging resistance, tall posture and high biomass production. Nuclear magnetic resonance imaging displayed the 3-dimensional assembly of vascular bundles in stem. A higher number of vascular bundles and a higher degree of their incline were the features of lodging-resistant versus lodging-prone lines. Histology and electron microscopy revealed that stems are fortified by a higher proportion of sclerenchyma and thickened cell walls, as well as some epidermal invaginations. Biochemical analysis using Fourier-transform infrared spectroscopy and inductively coupled plasma-optical emission spectrometry further identified elevated levels of lignin, xylan, zinc and silicon as features associated with high lodging resistance. Combined effects of above features caused superior culm stability. A simplistic mathematical model showed how mechanical forces distribute within the stem under stress. Main traits of the lodging-resistant parental line were heritable and could be traced back to the genetic structure of the mutant. Evaluation of lodging-resistant wheat 'Babax' ('Baviacora') versus contrasting, lodging-prone, genotype ´Pastor´ agreed with above findings on rye. Our findings on mechanical stability and extraordinary culm properties may be important for breeders for the improvement of lodging resistance of tall posture cereal crops.

Juvenile Heat Tolerance in Wheat for Attaining Higher Grain Yield by Shifting to Early Sowing in October in South Asia

Farmers in northwestern and central India have been exploring to sow their wheat much earlier (October) than normal (November) to sustain productivity by escaping terminal heat stress and to utilize the available soil moisture after the harvesting of rice crop. However, current popular varieties are poorly adapted to early sowing due to the exposure of juvenile plants to the warmer temperatures in the month of October and early November. Therefore, a study was undertaken to identify wheat genotypes suited to October sowing under warmer temperatures in India. A diverse collection of 3322 bread wheat varieties and elite lines was prepared in CIMMYT, Mexico, and planted in the 3rd week of October during the crop season 2012-2013 in six locations (Ludhiana, Karnal, New Delhi, Indore, Pune and Dharwad) spread over northwestern plains zone (NWPZ) and central and Peninsular zone (CZ and PZ; designated as CPZ) of India. Agronomic traits data from the seedling stage to maturity were recorded. Results indicated substantial diversity for yield and yield-associated traits, with some lines showing indications of higher yields under October sowing. Based on agronomic performance and disease resistance, the top 48 lines (and two local checks) were identified and planted in the next crop season (2013-2014) in a replicated trial in all six locations under October sowing (third week). High yielding lines that could tolerate higher temperature in October sowing were identified for both zones; however, performance for grain yield was more promising in the NWPZ. Hence, a new trial of 30 lines was planted only in NWPZ under October sowing. Lines showing significantly superior yield over the best check and the most popular cultivars in the zone were identified. The study suggested that agronomically superior wheat varieties with early heat tolerance can be obtained that can provide yield up to 8 t/ha by planting in the third to fourth week of October.

Quantitative trait loci for yellow rust resistance in spring wheat doubled haploid populations developed from the German Federal ex situ genebank genetic resources

Novel resistance sources to the pathogen Puccinia striiformis f. sp. tritici, which causes yellow rust (stripe rust), a widespread devastating foliar disease in wheat (Triticum aestivum L.), are in demand. Here, we tested two doubled haploid (DH) spring wheat populations derived from the genetic resources for resistance to yellow rust in field trials in Germany and Egypt. Additionally, we performed tests for all-stage resistance (seedling resistance). We performed linkage mapping based on 15k Infinium SNP chip genotyping data that resulted in 3,567 and 3,457 polymorphic markers for DH Population 1 (103 genotypes) and DH Population 2 (148 genotypes), respectively. In DH Population 1, we identified a major and consistent quantitative trait locus (QTL) on chromosome 1B that explained up to 28 and 39% of the phenotypic variation in the field and seedling tests, respectively. The favorable allele was contributed by the line 'TRI-5645', a landrace from Iran, and is most probably the yellow rust resistance (Yr) gene Yr10. In DH Population 2, the favorable allele of a major QTL on chromosome 6B was contributed by the line 'TRI-5310', representing the variety 'Eureke' from France. This QTL was mainly effective in the German environments and explained up to 36% of the phenotypic variation. In Egypt, however, only a moderate resistance QTL was identified in the field tests and no resistance QTL was observed in the seedling tests. Our results demonstrate the usefulness of genetic resources to identify novel sources of resistance to yellow rust, including the "Warrior" race PstS10.

Toward identification of a putative candidate gene for nutrient mineral accumulation in wheat grains for human nutrition purposes

A multilocus genome-wide association study of a panel of 369 diverse wheat (Triticum aestivum) genotypes was carried out in order to examine the genetic basis of variations in nutrient mineral concentrations in the grains. The panel was grown under field conditions for three consecutive years and the concentrations of Ca, K, Mg, Mn, P, and S were determined. Wide ranges of natural variation were detected among the genotypes. Strong positive correlations were found among the minerals except for K, which showed negative correlation trends with the other minerals. Genetic association analysis detected 86 significant marker-trait associations (MTAs) underlying the natural variations in mineral concentrations in grains. The major MTA was detected on the long arm of chromosome 5A and showed a pleiotropic effect on Ca, K, Mg, Mn, and S. Further significant MTAs were distributed among the whole genome except for chromosomes 3D and 6D. We identified putative candidate genes that are potentially involved in metal uptake, transport, and assimilation, including TraesCS5A02G542600 on chromosome 5A, which was annotated as a Major Facilitator Superfamily transporter and acted on all the minerals except K. TraesCS5A02G542600 was highly expressed in seed coat, and to a lesser extent in the peduncle, awns, and lemma. Our results provide important insights into the genetic basis of enhancement of nutrient mineral concentrations that can help to inform future breeding studies in order to improve human nutrition.

Omics for the Improvement of Abiotic, Biotic, and Agronomic Traits in Major Cereal Crops: Applications, Challenges, and Prospects

Omics technologies, namely genomics, transcriptomics, proteomics, metabolomics, and phenomics, are becoming an integral part of virtually every commercial cereal crop breeding program, as they provide substantial dividends per unit time in both pre-breeding and breeding phases. Continuous advances in omics assure time efficiency and cost benefits to improve cereal crops. This review provides a comprehensive overview of the established omics methods in five major cereals, namely rice, sorghum, maize, barley, and bread wheat. We cover the evolution of technologies in each omics section independently and concentrate on their use to improve economically important agronomic as well as biotic and abiotic stress-related traits. Advancements in the (1) identification, mapping, and sequencing of molecular/structural variants; (2) high-density transcriptomics data to study gene expression patterns; (3) global and targeted proteome profiling to study protein structure and interaction; (4) metabolomic profiling to quantify organ-level, small-density metabolites, and their composition; and (5) high-resolution, high-throughput, image-based phenomics approaches are surveyed in this review.

Omics in Major Cereals: Applications, Challenges, and Prospects.. [PMID: 0 DOI: 10.20944/preprints202104.0531.v1]

Omics technologies, viz., genomics, transcriptomics, proteomics, metabolomics, and phenomics, are becoming an integral part of virtually every commercial cereal breeding program because they provide substantial dividends per unit time in both pre-breeding and breeding phases. Continuous advances in cereal-omics promise—in combination with time efficiency—the cost benefits. In this review, we provide a comprehensive overview of the established cereal-omics methods in five major cereals, viz., rice, sorghum, maize, barley, and bread wheat. We cover the evolution of technologies in each omics section independently and concentrate on their use to improve economically important agronomic as well as biotic and abiotic stress-related traits. Advancements in the (1) identification, mapping, and sequencing of molecular/structural variants, (2) high-density transcriptomics data to study gene expression patterns, (3) global and targeted proteome profiling to study protein structure and interaction, (4) metabolomic profiling to quantify organ level small-density metabolites and their composition, and (5) high-resolution high-throughput image-based phenomics approaches are surveyed in this review.

Genome-Wide Association Studies and Prediction of Tan Spot (Pyrenophora tritici-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 simple sequence repeats (SSRs) plus 35k and 90k single nucleotide polymorphism (SNP) marker platforms. In our phenotypic data analyses, Ptr infection showed a significant genotypic variance and a significant negative correlation with plant height. Genome-wide association studies revealed a highly quantitative nature of Ptr infection and identified two quantitative trait loci (QTL), viz., QTs.ipk-7A and QTs.ipk-7B, which imparted 21.23 and 5.84% of the genotypic variance, respectively. Besides, the Rht-D1 gene showed a strong allelic influence on the infection scores. Due to the complex genetic nature of the Ptr infection, the potential of genome-wide prediction (GP) was assessed via three different genetic models on individual and combined marker platforms. The GP results indicated that the marker density and marker platforms do not considerably impact prediction accuracy (~40–42%) and that higher-order epistatic interactions may not be highly pervasive. Our results provide a further understanding of Ptr-infection’s genetic nature, serve as a resource for marker-assisted breeding, and highlight the potential of genome-wide selection for improved Ptr resistance.

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 (Triticumaestivum 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 FLOWERINGLOCUST (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.

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.

Tight DNA-protein complexes isolated from barley seedlings are rich in potential guanine quadruplex sequences

Background

The concept of chromatin domains attached to the nuclear matrix is being revisited, with nucleus described as a set of topologically associating domains. The significance of the tightly bound to DNA proteins (TBP), a protein group that remains attached to DNA after its deproteinization should be also revisited, as the existence of these interactions is in good agreement with the concept of the topologically associating domain. The work aimed to characterize the DNA component of TBP isolated from barley seedlings.

Methods

The tight DNA-protein complexes from the first leaves, coleoptiles, and roots of barley seedlings were isolated by purification with chromatography on nitrocellulose or exhaustive digestion of DNA with DNase I. Cloning and transformation were performed using pMOSBBlue Blunt Ended Cloning Kit. Inserts were amplified by PCR, and sequencing was performed on the MegaBace 1000 Sequencing System. The BLAST search was performed using sequence databases at NCBI, CR-EST, and TREP and Ensembl Plants databases. Comparison to MAR/SAR sequences was performed using http://smartdb.bioinf.med.uni-goettingen.de/cgi-bin/SMARtDB/smar.cgi database. The prediction of G quadruplexes (GQ) was performed with the aid of R-studio library pqsfinder. CD spectra were recorded on a Chirascan CS/3D spectrometer.

Results

Although the barley genome is AT-rich (43% of GC pairs), most DNA fragments associated with TBP were GC-rich (up to 70% in some fractions). Both fractionation procedures yielded a high proportion of CT-motif sequences presented predominantly by the 16-bp CC(TCTCCC)₂ TC fragment present in clones derived from the TBP-bound DNA and absent in free DNA. BLAST analysis revealed alignment with different barley repeats. Some clones, however, aligned with both nuclear and chloroplast structural genes. Alignments with MAR/SAR motifs were very few. The analysis produced by the pqsfinder program revealed numerous potential quadruplex-forming sites in the TBP-bound sequences. A set of oligonucleotides containing sites of possible GQs were designed and ordered. Three of them represented the minus strand of the CT-repeat. Two were derived from sequences of two clones of nitrocellulose retained fraction from leaves and contained GC-rich motifs different from the CT motif. Circular dichroism spectroscopy revealed profound changes in spectra when oligonucleotides were incubated with 100 mM KCl. There was either an increase of positive band in the area of 260 nm or the formation of a positive band at 290 nm. In the former case, changes are typical for parallel G-quadruplexes and, in the latter, 3 + 1 structures.

Discussion

The G-quadruplexes anchor proteins are probably involved in the maintenance of the topologically associated domain structure.

Suitability of Single-Nucleotide Polymorphism Arrays Versus Genotyping-By-Sequencing for Genebank Genomics in Wheat

Genebank genomics promises to unlock valuable diversity for plant breeding but first, one key question is which marker system is most suitable to fingerprint entire genebank collections. Using wheat as model species, we tested for the presence of an ascertainment bias and investigated its impact on estimates of genetic diversity and prediction ability obtained using three marker platforms: simple sequence repeat (SSR), genotyping-by-sequencing (GBS), and array-based SNP markers. We used a panel of 378 winter wheat genotypes including 190 elite lines and 188 plant genetic resources (PGR), which were phenotyped in multi-environmental trials for grain yield and plant height. We observed an ascertainment bias for the array-based SNP markers, which led to an underestimation of the molecular diversity within the population of PGR. In contrast, the marker system played only a minor role for the overall picture of the population structure and precision of genome-wide predictions. Interestingly, we found that rare markers contributed substantially to the prediction ability. This combined with the expectation that valuable novel diversity is most likely rare suggests that markers with minor allele frequency deserve careful consideration in the design of a pre-breeding program.

Automated Spike Detection in Diverse European Wheat Plants Using Textural Features and the Frangi Filter in 2D Greenhouse Images

Spike is one of the crop yield organs in wheat plants. Determination of the phenological stages, including heading time point (HTP), and area of spike from non-invasive phenotyping images provides the necessary information for the inference of growth-related traits. The algorithm previously developed by Qiongyan et al. for spike detection in 2-D images turns out to be less accurate when applied to the European cultivars that produce many more leaves. Therefore, we here present an improved and extended method where (i) wavelet amplitude is used as an input to the Laws texture energy-based neural network instead of original grayscale images and (ii) non-spike structures (e.g., leaves) are subsequently suppressed by combining the result of the neural network prediction with a Frangi-filtered image. Using this two-step approach, a 98.6% overall accuracy of neural network segmentation based on direct comparison with ground-truth data could be achieved. Moreover, the comparative error rate in spike HTP detection and growth correlation among the ground truth, the algorithm developed by Qiongyan et al., and the proposed algorithm are discussed in this paper. The proposed algorithm was also capable of significantly reducing the error rate of the HTP detection by 75% and improving the accuracy of spike area estimation by 50% in comparison with the Qionagyan et al. method. With these algorithmic improvements, HTP detection on a diverse set of 369 plants was performed in a high-throughput manner. This analysis demonstrated that the HTP of 104 plants (comprises of 57 genotypes) with lower biomass and tillering range (e.g., earlier-heading types) were correctly determined. However, fine-tuning or extension of the developed method is required for high biomass plants where spike emerges within green bushes. In conclusion, our proposed method allows significantly more reliable results for HTP detection and spike growth analysis to be achieved in application to European cultivars with earlier-heading types.

Identification of consistent QTL with large effect on anther extrusion in doubled haploid populations developed from spring wheat accessions in German Federal ex situ Genebank

Novel large-effect consistent QTL for anther extrusion (AE) to improve cross-pollination were mapped in doubled haploid populations derived from IPK gene bank spring wheat accessions. TaAP2-D, an ortholog of Cleistogamy1 in barley, is a likely candidate gene for AE in wheat. To establish a robust hybrid wheat breeding system, male lines harboring alleles that promote outcrossing should be developed. In this study, we developed two doubled haploid (DH) populations of hexaploid spring wheat (Triticum aestivum L.) by crossing accessions taken from IPK gene bank. In both populations, the phenotypic data of anther extrusion (AE) based on three years of field trials showed a wide variation and approximated a normal distribution. Both populations were genotyped with a 15 k Infinium single nucleotide polymorphism (SNP) array resulting in 3567 and 3457 polymorphic SNP markers for DH population-1 and DH population-2, respectively. Composite interval mapping identified quantitative trait loci (QTL) on chromosomes 1D, 2D, 4A, 4B, 5A, 5D, 6A, and 6B; with consistent QTL (that are identified in all the years) on chromosome 4A in DH population-1, and on chromosomes 2D and 6B in DH population-2. The consistent QTL explained 17.2%, 32.9%, and 12.3% of the phenotypic variances, respectively. Genic scan of the chromosome 2D-QTL showed that the wheat gene TaAP2-D, an ortholog of Cleistogamy1 which promotes AE via swelling of the lodicules in barley, lies within the QTL region. A diagnostic marker was developed for TaAP2-D that showed co-segregation with the AE phenotype. This study shows the use of gene bank diversity reservoir to find alleles which are otherwise difficult to detect in elite populations. The identification of large-effect consistent QTL for AE is expected to help form efficient male parental lines suitable for hybrid wheat seed production and serve as a source for map-based cloning.

Genetic analysis of drought response of wheat following either chemical desiccation or the use of a rain-out shelter

Simulating drought stress during the breeding process has been proposed as a way to select varieties under naturally non-stressful conditions. The aim of the study was to characterise the genetic basis of the response of 111 spring wheat (Triticum aestivum L.) varieties and landraces to chemical desiccation and to rain-out shelter drought. The effect of the rain-out shelter was a 15% reduction in plant height, spike length and thousand seed weight (TSW); in contrast, the desiccant treatment induced a 15% reduction in seed number, a 35-72% loss in TSW and a reduction in subsequent germination of 12%. A genome-wide association analysis revealed 263 significant marker-trait associations (MTAs), of which 246 involved days to anthesis, plant height, spike length, number of spikelets, seed number, TSW and germination from the non-treated plants. Only four and five MTAs involved TSW from plants grown under the rain-out shelter and the chemical desiccation, respectively, and harboured the Sugar-Dependent6 gene. Seven MTAs involved seed number for chemical desiccated plants. Both, chemical desiccation and rain-out shelter drought identified same tolerant genotypes. Concluding, both approaches are suitable to simulate different drought scenarios. However, there was a strong environmental impact for chemical desiccation which may increase the complexity of this tolerance mechanism.

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.

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.

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.

Genome-wide association mapping in bread wheat subjected to independent and combined high temperature and drought stress

Most investigations to date aiming to identify the genetic basis of the stress response of wheat (Triticum aestivum L.) have focused on the response to single stress agents such as high temperature or drought, even though in the natural situation, these stresses tend often to occur together. Here, a panel of 108 spring type bread wheat cultivars was phenotyped for 15 yield and yield related traits for two years (2014/15 and 2015/16) under non-stressed conditions, under high temperature stress, under drought and under a combined high temperature and drought regime. The mean loss in grain yield caused by all stress regimes was 51.33%. Analysis of variance (ANOVA) of yield trait showed significant differences among genotypes environments and their interactions (G×E). All the studied traits had higher heritability values which ranged from 0.35-0.94 under [C], 0.58-0.95 under [D], 0.62-0.93 under [H] and 0.60-0.95 under [HD]. GWAS was performed by using 9,646 informative SNP markers and based on these polymorphic SNPs population structure analysis divided whole germplasm into five major sub-populations. Mixed model association analysis detected 503 marker-trait associations (MTAs) at P ≤0.001 while 329 MTAs crossed FDR ≤ 0.05 for all traits with phenotypic variances (R2) ranged from 24.83% to 12.51%. Seven new pleiotropic SNPs on chromosome 7D and IAAV8258 (86.91cM) and wsnp_Ex_c7168_12311649 (57.93cM) on chromosome 5A were most stable association in present study. Furthermore, candidate genes Psy and Sr25 (TG0040) were also significant in present study, these genes were previously mapped on 7A and 7D. The region on 7D was assiociated with 7DL.7Ag translocation from Lophopyrum carring rust resistance Yr16 and many other genes. Similarly region on chromosome 7A which was associated with Psy gene was linked with grain yellow pigment content QTLs. Favourable alleles controlling grain yield were of vital importance and incorporation of these alleles after validation through marker assisted selection and fine mapping could be helpful in wheat yield improvement stress and non-stress conditions.

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.

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.

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.

Genetic and physical mapping of anther extrusion in elite European winter wheat

The production and cultivation of hybrid wheat is a possible strategy to close the yield gap in wheat. Efficient hybrid wheat seed production largely depends on high rates of cross-pollination which can be ensured through high anther extrusion (AE) by male parental lines. Here, we report the AE capacity and elucidate its genetics in 514 elite European winter wheat varieties via genome-wide association studies (GWAS). We observed a wide range of variation among genotypes and a high heritability (0.80) for AE. The whole panel was genotyped with the 35k Affymetrix and 90k iSELECT single nucleotide polymorphism (SNP) arrays plus Ppd-D1, Rht-B1 and Rht-D1 candidate markers. GWAS revealed 51 marker-trait associations (MTAs) on chromosomes 1A, 1B, 2A, 4D and 5B, with Rht-D1 (4D) being the most significant marker. Division of whole panel according to the Rht-D1 genotype resulted in 212 and 294 varieties harboring Rht-D1a and Rht-D1b allele, respectively. The presence of Rht-D1a compared to Rht-D1b (mutant) allele had a large phenotypic influence on AE resulting in its ~17% increase. GWAS performed on the sub-panels detected novel MTAs on chromosomes 2D, 3B and 6A with increased phenotypic variance imparted by individual markers. Our study shows that AE is a highly quantitative trait and wild type Rht-D1a allele greatly improves AE. Moreover, demarcating the quantitative trait loci regions based on intra-chromosomal linkage disequilibrium revealed AE’s candidate genes/genomic regions. Understanding the genetics of AE in elite European wheat and utilizing the linked markers in breeding programs can help to enhance cross-pollination for better exploitation of heterosis.

Genome-Wide Association Study of Calcium Accumulation in Grains of European Wheat Cultivars

Mineral concentrations in cereals are important for human health, especially for people who depend mainly on consuming cereal diet. In this study, we carried out a genome-wide association study (GWAS) of calcium concentrations in wheat (Triticum aestivum L.) grains using a European wheat diversity panel of 353 varieties [339 winter wheat (WW) plus 14 of spring wheat (SW)] and phenotypic data based on two field seasons. High genotyping densities of single-nucleotide polymorphism (SNP) markers were obtained from the application of the 90k iSELECT ILLUMINA chip and a 35k Affymetrix chip. Inductively coupled plasma optical emission spectrometry (ICP-OES) was used to measure the calcium concentrations of the wheat grains. Best linear unbiased estimates (BLUEs) for calcium were calculated across the seasons and ranged from 288.20 to 647.50 among the varieties (μg g-1 DW) with a mean equaling 438.102 (μg g-1 DW), and the heritability was 0.73. A total of 485 SNP marker-trait associations (MTAs) were detected in data obtained from grains cultivated in both of the two seasons and BLUE values by considering associations with a -log10 (P-value) ≥3.0. Among these SNP markers, we detected 276 markers with a positive allele effect and 209 markers with a negative allele effect. These MTAs were found on all chromosomes except chromosomes 3D, 4B, and 4D. The most significant association was located on chromosome 5A (114.5 cM) and was linked to a gene encoding cation/sugar symporter activity as a potential candidate gene. Additionally, a number of candidate genes for the uptake or transport of calcium were located near significantly associated SNPs. This analysis highlights a number of genomic regions and candidate genes for further analysis as well as the challenges faced when mapping environmentally variable traits in genetically highly diverse variety panels. The research demonstrates the feasibility of the GWAS approach for illuminating the genetic architecture of calcium-concentration in wheat grains and for identifying putative candidate genes underlying this trait.

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 analyses of 54 traits identified multiple loci for the determination of floret fertility in wheat

Increasing grain yield is still the main target of wheat breeding; yet today's wheat plants utilize less than half of their yield potential. Owing to the difficulty of determining grain yield potential in a large population, few genetic factors regulating floret fertility (i.e. the difference between grain yield potential and grain number) have been reported to date. In this study, we conducted a genome-wide association study (GWAS) by quantifying 54 traits (16 floret fertility traits and 38 traits for assimilate partitioning and spike morphology) in 210 European winter wheat accessions. The results of this GWAS experiment suggested potential associations between floret fertility, assimilate partitioning and spike morphology revealed by shared quantitative trait loci (QTLs). Several candidate genes involved in carbohydrate metabolism, phytohormones or floral development colocalized with such QTLs, thereby providing potential targets for selection. Based on our GWAS results we propose a genetic network underlying floret fertility and related traits, nominating determinants for improved yield performance.

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

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.

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

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 ² = 0.91. A total of 108 SSR and 235 SNP marker-trait associations (MTAs) were identified by considering associations with a -log₁₀ (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.

Genetic Architecture of Anther Extrusion in Spring and Winter Wheat

Hybrid wheat breeding is gaining prominence worldwide because it ensures higher and more static yield than conventionally bred varieties. The cleistogamous floral architecture of wheat (Triticum aestivum L.) impedes anthers inside the floret, making it largely an inbreeder. For hybrid seed production, high anther extrusion is needed to promote cross pollination and to ensure a high level of pollen availability for the seed plant. This study, therefore, aimed at the genetic dissection of anther extrusion (AE) in panels of spring (SP), and winter wheat (WP) accessions by genome wide association studies (GWAS). We performed GWAS to identify the SNP markers potentially linked with AE in each panel separately. Phenotypic data were collected for 3 years for each panel. The average levels of Pearson's correlation (r) among all years and their best linear unbiased estimates (BLUEs) within both panels were high (r(SP) = 0.75, P < 0.0001;r(WP) = 0.72, P < 0.0001). Genotypic data (with minimum of 0.05 minor allele frequency applied) included 12,066 and 12,191 SNP markers for SP and WP, respectively. Both genotypes and environment influenced the magnitude of AE. In total, 23 significant (|log₁₀(P)| > 3.0) marker trait associations (MTAs) were detected (SP = 11; WP = 12). Anther extrusion behaved as a complex trait with significant markers having either favorable or unfavorable additive effects and imparting minor to moderate levels of phenotypic variance (R²(SP) = 9.75-14.24%; R² (WP) = 9.44-16.98%). All mapped significant markers as well as the markers within their significant linkage disequilibrium (r² ≥ 0.30) regions were blasted against wheat genome assembly (IWGSC1+popseq) to find the corresponding genes and their high confidence descriptions were retrieved. These genes and their orthologs in Hordeum vulgare, Brachypodium distachyon, Oryza sativa, and Sorghum bicolor revealed syntenic genomic regions potentially involved in flowering-related traits. Moreover, the expression data of these genes suggested potential candidates for AE. Our results suggest that the use of significant markers can help to introduce AE in high yielding varieties to increase cross fertilization rates and improve hybrid-seed production in wheat.

Genome-Wide Association Mapping of Anther Extrusion in Hexaploid Spring Wheat

In a number of crop species hybrids are able to outperform line varieties. The anthers of the autogamous bread wheat plant are normally extruded post anthesis, a trait which is unfavourable for the production of F1 hybrid grain. Higher anther extrusion (AE) promotes cross fertilization for more efficient hybrid seed production. Therefore, this study aimed at the genetic dissection of AE by genome wide association studies (GWAS) and determination of the main effect QTL. We applied GWAS approach to identify DArT markers potentially linked to AE to unfold its genetic basis in a panel of spring wheat accessions. Phenotypic data were collected for three years and best linear unbiased estimate (BLUE) values were calculated across all years. The extent of the AE correlation between growing years and BLUE values ranged from r = +0.56 (2013 vs 2015) to 0.91 (2014 vs BLUE values). The broad sense heritability was 0.84 across all years. Six accessions displayed stable AE >80% across all the years. Genotyping data included 2,575 DArT markers (with minimum of 0.05 minor allele frequency applied). AE was influenced both by genotype and by the growing environment. In all, 131 significant marker trait associations (MTAs) (|log₁₀ (P)| >FDR) were established for AE. AE behaved as a quantitative trait, with five consistently significant markers (significant across at least two years with a significant BLUE value) contributing a minor to modest proportion (4.29% to 8.61%) of the phenotypic variance and affecting the trait either positively or negatively. For this reason, there is potential for breeding for improved AE by gene pyramiding. The consistently significant markers linked to AE could be helpful for marker assisted selection to transfer AE to high yielding varieties allowing to promote the exploitation of hybrid-heterosis in the key crop wheat.

Are PECTIN ESTERASE INHIBITOR Genes Involved in Mediating Resistance to Rhynchosporium commune in Barley?

A family of putative PECTIN ESTERASE INHIBITOR (PEI) genes, which were detected in the genomic region co-segregating with the resistance gene Rrs2 against scald caused by Rhynchosporium commune in barley, were characterized and tested for their possible involvement in mediating resistance to the pathogen by complementation and overexpression analysis. The sequences of the respective genes were derived from two BAC contigs originating from the susceptible cultivar ‘Morex’. For the genes HvPEI2, HvPEI3, HvPEI4 and HvPEI6, specific haplotypes for 18 resistant and 23 susceptible cultivars were detected after PCR-amplification and haplotype-specific CAPS-markers were developed. None of the tested candidate genes HvPEI2, HvPEI3 and HvPEI4 alone conferred a high resistance level in transgenic over-expression plants, though an improvement of the resistance level was observed especially with OE-lines for gene HvPEI4. These results do not confirm but also do not exclude an involvement of the PEI gene family in the response to the pathogen. A candidate for the resistance gene Rrs2 could not be identified yet. It is possible that Rrs2 is a PEI gene or another type of gene which has not been detected in the susceptible cultivar ‘Morex’ or the full resistance reaction requires the presence of several PEI genes.

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.

Potential and limits to unravel the genetic architecture and predict the variation of Fusarium head blight resistance in European winter wheat (Triticum aestivum L.)

Genome-wide mapping approaches in diverse populations are powerful tools to unravel the genetic architecture of complex traits. The main goals of our study were to investigate the potential and limits to unravel the genetic architecture and to identify the factors determining the accuracy of prediction of the genotypic variation of Fusarium head blight (FHB) resistance in wheat (Triticum aestivum L.) based on data collected with a diverse panel of 372 European varieties. The wheat lines were phenotyped in multi-location field trials for FHB resistance and genotyped with 782 simple sequence repeat (SSR) markers, and 9k and 90k single-nucleotide polymorphism (SNP) arrays. We applied genome-wide association mapping in combination with fivefold cross-validations and observed surprisingly high accuracies of prediction for marker-assisted selection based on the detected quantitative trait loci (QTLs). Using a random sample of markers not selected for marker-trait associations revealed only a slight decrease in prediction accuracy compared with marker-based selection exploiting the QTL information. The same picture was confirmed in a simulation study, suggesting that relatedness is a main driver of the accuracy of prediction in marker-assisted selection of FHB resistance. When the accuracy of prediction of three genomic selection models was contrasted for the three marker data sets, no significant differences in accuracies among marker platforms and genomic selection models were observed. Marker density impacted the accuracy of prediction only marginally. Consequently, genomic selection of FHB resistance can be implemented most cost-efficiently based on low- to medium-density SNP arrays.

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.

Development of SNP markers for genes of the phenylpropanoid pathway and their association to kernel and malting traits in barley

BACKGROUND

Flavonoids are an important class of secondary compounds in angiosperms. Next to certain biological functions in plants, they play a role in the brewing process and have an effect on taste, color and aroma of beer. The aim of this study was to reveal the haplotype diversity of candidate genes involved in the phenylpropanoid biosynthesis pathway in cultivated barley varieties (Hordeum vulgare L.) and to determine associations to kernel and malting quality parameters.

RESULTS

Five genes encoding phenylalanine ammonia-lyase (PAL), cinnamate 4-hydroxylase (C4H), chalcone synthase (CHS), flavanone 3-hydroxylase (F3H) and dihydroflavonol reductase (DFR) of the phenylpropanoid biosynthesis pathway were partially resequenced in 16 diverse barley reference genotypes. Their localization in the barley genome, their genetic structure, and their genetic variation e.g. single nucleotide polymorphism (SNP) and Insertion/Deletion (InDel) patterns were revealed. In total, 130 SNPs and seven InDels were detected. Of these, 21 polymorphisms were converted into high-throughput pyrosequencing markers. The resulting SNP and haplotype patterns were used to calculate associations with kernel and malting quality parameters.

CONCLUSIONS

SNP patterns were found to be highly variable for the investigated genes. The developed high-throughput markers are applicable for assessing the genetic variability and for the determination of haplotype patterns in a set of barley accessions. The candidate genes PAL, C4H and F3H were shown to be associated to several malting properties like glassiness (PAL), viscosity (C4H) or to final attenuation (F3H).

Nitrogen-metabolism related genes in barley - haplotype diversity, linkage mapping and associations with malting and kernel quality parameters

Background

Several studies report about intra-specific trait variation of nitrogen-metabolism related traits, such as N(itrogen)-use efficiency, protein content, N-storage and remobilization in barley and related grass species. The goal of this study was to assess the intra-specific genetic diversity present in primary N-metabolism genes of barley and to investigate the associations of the detected haplotype diversity with malting and kernel quality related traits.

Results

Partial sequences of five genes related to N-metabolism in barley (Hordeum vulgare L.) were obtained, i.e. nitrate reductase 1, glutamine synthetase 2, ferredoxin-dependent glutamate synthase, aspartate aminotransferase and asparaginase. Two to five haplotypes in each gene were discovered in a set of 190 various varieties. The development of 33 SNP markers allowed the genotyping of all these barley varieties consisting of spring and winter types. Furthermore, these markers could be mapped in several doubled haploid populations. Cluster analysis based on haplotypes revealed a more uniform pattern of the spring barleys as compared to the winter barleys. Based on linear model approaches associations to several malting and kernel quality traits including soluble N and protein were identified.

Conclusions

A study was conducted to investigate the presence of sequence variation of several genes related to the primary N-metabolism in barley. The detected diversity could be related to particular phenotypic traits. Specific differences between spring and winter barleys most likely reflect different breeding aims. The developed markers can be used as tool for further genetic studies and marker-assisted selection during breeding of barley.

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.

QTL mapping of terminal heat tolerance in hexaploid wheat (T. aestivum L.)

High temperature (>30 °C) at the time of grain filling is one of the major causes of yield reduction in wheat in many parts of the world, especially in tropical countries. To identify quantitative trait loci (QTL) for heat tolerance under terminal heat stress, a set of 148 recombinant inbred lines was developed by crossing a heat-tolerant hexaploid wheat (Triticum aestivum L.) cultivar (NW1014) and a heat-susceptible (HUW468) cultivar. The F(5), F(6), and F(7) generations were evaluated in two different sowing dates under field conditions for 2 years. Using the trait values from controlled and stressed trials, four different traits (1) heat susceptibility index (HSI) of thousand grain weight (HSITGW); (2) HSI of grain fill duration (HSIGFD); (3) HSI of grain yield (HSIYLD); and (4) canopy temperature depression (CTD) were used to determine heat tolerance. Days to maturity was also investigated. A linkage map comprising 160 simple sequence repeat markers was prepared covering the whole genome of wheat. Using composite interval mapping, significant genomic regions on 2B, 7B and 7D were found to be associated with heat tolerance. Of these, two (2B and 7B) were co-localized QTL and explained more than 15 % phenotypic variation for HSITGW, HSIGFD and CTD. In pooled analysis over three trials, QTL explained phenotypic variation ranging from 9.78 to 20.34 %. No QTL × trial interaction was detected for the identified QTL. The three major QTL obtained can be used in marker-assisted selection for heat stress in wheat.

Identification and characterization of a novel powdery mildew resistance gene PmG3M derived from wild emmer wheat, Triticum dicoccoides

Powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt) is one of the most important wheat diseases worldwide. Wild emmer wheat, Triticum turgidum ssp. dicoccoides, the tetraploid ancestor (AABB) of domesticated bread and durum wheat, harbors many important alleles for resistance to various diseases, including powdery mildew. In the current study, two tetraploid wheat mapping populations, derived from a cross between durum wheat (cv. Langdon) and wild emmer wheat (accession G-305-3M), were used to identify and map a novel powdery mildew resistance gene. Wild emmer accession G-305-3M was resistant to all 47 Bgt isolates tested, from Israel and Switzerland. Segregation ratios of F(2) progenies and F(6) recombinant inbred line (RIL) mapping populations, in their reactions to inoculation with Bgt, revealed a Mendelian pattern (3:1 and 1:1, respectively), indicating the role of a single dominant gene derived from T. dicoccoides accession G-305-3M. This gene, temporarily designated PmG3M, was mapped on chromosome 6BL and physically assigned to chromosome deletion bin 6BL-0.70-1.00. The F(2) mapping population was used to construct a genetic map of the PmG3M gene region consisted of six simple sequence repeats (SSR), 11 resistance gene analog (RGA), and two target region amplification polymorphism (TRAP) markers. A second map, constructed based on the F(6) RIL population, using a set of skeleton SSR markers, confirmed the order of loci and distances obtained for the F(2) population. The discovery and mapping of this novel powdery mildew resistance gene emphasize the importance of the wild emmer wheat gene pool as a source for crop improvement.

High-density genetic and physical bin mapping of wheat chromosome 1D reveals that the powdery mildew resistance gene Pm24 is located in a highly recombinogenic region

Genetic maps of wheat chromosome 1D consisting of 57 microsatellite marker loci were constructed using Chinese Spring (CS) × Chiyacao F(2) and the International Triticeae Mapping Initiative (ITMI) recombinant inbred lines (RILs) mapping populations. Marker order was consistent, but genetic distances of neighboring markers were different in two populations. Physical bin map of 57 microsatellite marker loci was generated by means of 10 CS 1D deletion lines. The physical bin mapping indicated that microsatellite marker loci were not randomly distributed on chromosome 1D. Nineteen of the 24 (79.2%) microsatellite markers were mapped in the distal 30% genomic region of 1DS, whereas 25 of the 33 (75.8%) markers were assigned to the distal 59% region of 1DL. The powdery mildew resistance gene Pm24, originating from the Chinese wheat landrace Chiyacao, was previously mapped in the vicinity of the centromere on the short arm of chromosome 1D. A high density genetic map of chromosome 1D was constructed, consisting of 36 markers and Pm24, with a total map length of 292.7 cM. Twelve marker loci were found to be closely linked to Pm24. Pm24 was flanked by Xgwm789 (Xgwm603) and Xbarc229 with genetic distances of 2.4 and 3.6 cM, respectively, whereas a microsatellite marker Xgwm1291 co-segregated with Pm24. The microsatellite marker Xgwm1291 was assigned to the bin 1DS5-0.70-1.00 of the chromosome arm 1DS. It could be concluded that Pm24 is located in the '1S0.8 gene-rich region', a highly recombinogenic region of wheat. The results presented here would provide a start point for the map-based cloning of Pm24.

Haplotyping, linkage mapping and expression analysis of barley genes regulated by terminal drought stress influencing seed quality

BACKGROUND

The increasingly narrow genetic background characteristic of modern crop germplasm presents a challenge for the breeding of cultivars that require adaptation to the anticipated change in climate. Thus, high priority research aims at the identification of relevant allelic variation present both in the crop itself as well as in its progenitors. This study is based on the characterization of genetic variation in barley, with a view to enhancing its response to terminal drought stress.

RESULTS

The expression patterns of drought regulated genes were monitored during plant ontogeny, mapped and the location of these genes was incorporated into a comprehensive barley SNP linkage map. Haplotypes within a set of 17 starch biosynthesis/degradation genes were defined, and a particularly high level of haplotype variation was uncovered in the genes encoding sucrose synthase (types I and II) and starch synthase. The ability of a panel of 50 barley accessions to maintain grain starch content under terminal drought conditions was explored.

CONCLUSION

The linkage/expression map is an informative resource in the context of characterizing the response of barley to drought stress. The high level of haplotype variation among starch biosynthesis/degradation genes in the progenitors of cultivated barley shows that domestication and breeding have greatly eroded their allelic diversity in current elite cultivars. Prospective association analysis based on core drought-regulated genes may simplify the process of identifying favourable alleles, and help to understand the genetic basis of the response to terminal drought.

Assessment of genetic diversity among Syrian durum (Triticum turgidum ssp. durum) and bread wheat (Triticum aestivum L.) using SSR markers

Genetic diversity among 49 wheat varieties (37 durum and 12 bread wheat) was assayed using 32 microsatellites representing 34 loci covering almost the whole wheat genome. The polymorphic information content (PIC) across the tested loci ranged from 0 to 0.88 with average values of 0.57 and 0.65 for durum and bread wheat respectively. B genome had the highest mean number of alleles (10.91) followed by A genome (8.3) whereas D genome had the lowest number (4.73). The correlation between PIC and allele number was significant in all genome groups accounting for 0.87, 074 and 0.84 for A, B and D genomes respectively, and over all genomes, the correlation was higher in tetraploid (0.8) than in hexaploid wheat varieties (0.5). The cluster analysis discriminated all varieties and clearly divided the two ploidy levels into two separate clusters that reflect the differences in genetic diversity within each cluster. This study demonstrates that microsatellites markers have unique advantages compared to other molecular and biochemical fingerprinting techniques in revealing the genetic diversity in Syrian wheat varieties that is crucial for wheat improvement.

Delineating the structural, functional and evolutionary relationships of sucrose phosphate synthase gene family II in wheat and related grasses

BACKGROUND

Sucrose phosphate synthase (SPS) is an important component of the plant sucrose biosynthesis pathway. In the monocotyledonous Poaceae, five SPS genes have been identified. Here we present a detailed analysis of the wheat SPSII family in wheat. A set of homoeologue-specific primers was developed in order to permit both the detection of sequence variation, and the dissection of the individual contribution of each homoeologue to the global expression of SPSII.

RESULTS

The expression in bread wheat over the course of development of various sucrose biosynthesis genes monitored on an Affymetrix array showed that the SPS genes were regulated over time and space. SPSII homoeologue-specific assays were used to show that the three homoeologues contributed differentially to the global expression of SPSII. Genetic mapping placed the set of homoeoloci on the short arms of the homoeologous group 3 chromosomes. A resequencing of the A and B genome copies allowed the detection of four haplotypes at each locus. The 3B copy includes an unspliced intron. A comparison of the sequences of the wheat SPSII orthologues present in the diploid progenitors einkorn, goatgrass and Triticum speltoides, as well as in the more distantly related species barley, rice, sorghum and purple false brome demonstrated that intronic sequence was less well conserved than exonic. Comparative sequence and phylogenetic analysis of SPSII gene showed that false purple brome was more similar to Triticeae than to rice. Wheat - rice synteny was found to be perturbed at the SPS region.

CONCLUSION

The homoeologue-specific assays will be suitable to derive associations between SPS functionality and key phenotypic traits. The amplicon sequences derived from the homoeologue-specific primers are informative regarding the evolution of SPSII in a polyploid context.