Genome-Wide Association Study Uncovers Novel Genomic Regions Associated With Coleoptile Length in Hard Winter Wheat

Characterization of flag leaf morphology identifies a major genomic region controlling flag leaf angle in the US winter wheat (Triticum aestivum L.)

KEY MESSAGE

Multi-environmental characterization of flag leaf morphology traits in the US winter wheat revealed nine stable genomic regions for different flag leaf-related traits including a major region governing flag leaf angle. Flag leaf in wheat is the primary contributor to accumulating photosynthetic assimilates. Flag leaf morphology (FLM) traits determine the overall canopy structure and capacity to intercept the light, thus influencing photosynthetic efficiency. Hence, understanding the genetic control of these traits could be useful for breeding desirable ideotypes in wheat. We used a panel of 272 accessions from the hard winter wheat (HWW) region of the USA to investigate the genetic architecture of five FLM traits including flag leaf length (FLL), width (FLW), angle (FLANG), length-width ratio, and area using multilocation field experiments. Multi-environment GWAS using 14,537 single-nucleotide polymorphisms identified 36 marker-trait associations for different traits, with nine being stable across environments. A novel and major stable region for FLANG (qFLANG.1A) was identified on chromosome 1A accounting for 9-13% variation. Analysis of spatial distribution for qFLANG.1A in a set of 2354 breeding lines from the HWW region showed a higher frequency of allele associated with narrow leaf angle. A KASP assay was developed for allelic discrimination of qFLANG.1A and was used for its independent validation in a diverse set of spring wheat accessions. Furthermore, candidate gene analysis for two regions associated with FLANG identified seven putative genes of interest for each of the two regions. The present study enhances our understanding of the genetic control of FLM in wheat, particularly FLANG, and these results will be useful for dissecting the genes underlying canopy architecture in wheat facilitating the development of climate-resilient wheat varieties.

Genome-wide association mapping of genomic regions associated with drought stress tolerance at seedling and reproductive stages in bread wheat

Understanding the genetic architecture of drought stress tolerance in bread wheat at seedling and reproductive stages is crucial for developing drought-tolerant varieties. In the present study, 192 diverse wheat genotypes, a subset from the Wheat Associated Mapping Initiative (WAMI) panel, were evaluated at the seedling stage in a hydroponics system for chlorophyll content (CL), shoot length (SLT), shoot weight (SWT), root length (RLT), and root weight (RWT) under both drought and optimum conditions. Following that, a genome-wide association study (GWAS) was carried out using the phenotypic data recorded during the hydroponics experiment as well as data available from previously conducted multi-location field trials under optimal and drought stress conditions. The panel had previously been genotyped using the Infinium iSelect 90K SNP array with 26,814 polymorphic markers. Using single as well as multi-locus models, GWAS identified 94 significant marker-trait associations (MTAs) or SNPs associated with traits recorded at the seedling stage and 451 for traits recorded at the reproductive stage. The significant SNPs included several novel, significant, and promising MTAs for different traits. The average LD decay distance for the whole genome was approximately 0.48 Mbp, ranging from 0.07 Mbp (chromosome 6D) to 4.14 Mbp (chromosome 2A). Furthermore, several promising SNPs revealed significant differences among haplotypes for traits such as RLT, RWT, SLT, SWT, and GY under drought stress. Functional annotation and in silico expression analysis revealed important putative candidate genes underlying the identified stable genomic regions such as protein kinases, O-methyltransferases, GroES-like superfamily proteins, NAD-dependent dehydratases, etc. The findings of the present study may be useful for improving yield potential, and stability under drought stress conditions.

Impact of "Green Revolution" gene Rht-B1b on coleoptile length of wheat

Wheat coleoptile is a sheath-like structure that helps to deliver the first leaf from embryo to the soil surface. Here, a RIL population consisting of 245 lines derived from Zhou 8425B × Chinese Spring cross was genotyped by the high-density Illumina iSelect 90K assay for coleoptile length (CL) QTL mapping. Three QTL for CL were mapped on chromosomes 2BL, 4BS and 4DS. Of them, two major QTL QCL.qau-4BS and QCL.qau-4DS were detected, which could explain 9.1%-22.2% of the phenotypic variances across environments on Rht-B1 and Rht-D1 loci, respectively. Several studies have reported that Rht-B1b may reduce the length of wheat CL but no study has been carried out at molecular level. In order to verify that the Rht-B1 gene is the functional gene for the 4B QTL, an overexpression line Rht-B1b-OE and a CRISPR/SpCas9 line Rht-B1b-KO were studied. The results showed that Rht-B1b overexpression could reduce the CL, while loss-of-function of Rht-B1b would increase the CL relative to that of the null transgenic plants (TNL). To dissect the underlying regulatory mechanism of Rht-B1b on CL, comparative RNA-Seq was conducted between Rht-B1b-OE and TNL. Transcriptome profiles revealed a few key pathways involving the function of Rht-B1b in coleoptile development, including phytohormones, circadian rhythm and starch and sucrose metabolism. Our findings may facilitate wheat breeding for longer coleoptiles to improve seedling early vigor for better penetration through the soil crust in arid regions.

Genome-wide association analysis of spike and kernel traits in the U.S. hard winter wheat

A better understanding of the genetic control of spike and kernel traits that have higher heritability can help in the development of high-yielding wheat varieties. Here, we identified the marker-trait associations (MTAs) for various spike- and kernel-related traits in winter wheat (Triticum aestivum L.) through genome-wide association studies (GWAS). An association mapping panel comprising 297 hard winter wheat accessions from the U.S. Great Plains was evaluated for eight spike- and kernel-related traits in three different environments. A GWAS using 15,590 single-nucleotide polymorphisms (SNPs) identified a total of 53 MTAs for seven spike- and kernel-related traits, where the highest number of MTAs were identified for spike length (16) followed by the number of spikelets per spike (15) and spikelet density (11). Out of 53 MTAs, 14 were considered to represent stable quantitative trait loci (QTL) as they were identified in multiple environments. Five multi-trait MTAs were identified for various traits including the number of spikelets per spike (NSPS), spikelet density (SD), kernel width (KW), and kernel area (KA) that could facilitate the pyramiding of yield-contributing traits. Further, a significant additive effect of accumulated favorable alleles on the phenotype of four spike-related traits suggested that breeding lines and cultivars with a higher number of favorable alleles could be a valuable resource for breeders to improve yield-related traits. This study improves the understanding of the genetic basis of yield-related traits in hard winter wheat and provides reliable molecular markers that will facilitate marker-assisted selection (MAS) in wheat breeding programs.

Genome-wide association study to identify genomic loci associated with early vigor in bread wheat under simulated water deficit complemented with quantitative trait loci meta-analysis

A genome-wide association study (GWAS) was used to identify associated loci with early vigor under simulated water deficit and grain yield under field drought in a diverse collection of Iranian bread wheat landraces. In addition, a meta-quantitative trait loci (MQTL) analysis was used to further expand our approach by retrieving already published quantitative trait loci (QTL) from recombinant inbred lines, double haploids, back-crosses, and F2 mapping populations. In the current study, around 16%, 14%, and 16% of SNPs were in significant linkage disequilibrium (LD) in the A, B, and D genomes, respectively, and varied between 5.44% (4A) and 21.85% (6A). Three main subgroups were identified among the landraces with different degrees of admixture, and population structure was further explored through principal component analysis. Our GWAS identified 54 marker-trait associations (MTAs) that were located across the wheat genome but with the highest number found in the B sub-genome. The gene ontology (GO) analysis of MTAs revealed that around 75% were located within or closed to protein-coding genes. In the MQTL analysis, 23 MQTLs, from a total of 215 QTLs, were identified and successfully projected onto the reference map. MQT-YLD4, MQT-YLD9, MQT-YLD13, MQT-YLD17, MQT-YLD18, MQT-YLD19, and MQTL-RL1 contributed to the highest number of projected QTLs and were therefore regarded as the most reliable and stable QTLs under water deficit conditions. These MQTLs greatly facilitate the identification of putative candidate genes underlying at each MQTL interval due to the reduced confidence of intervals associated with MQTLs. These findings provide important information on the genetic basis of early vigor traits and grain yield under water deficit conditions and set the foundation for future investigations into adaptation to water deficit in bread wheat.

Identification of leaf rust resistance loci in a geographically diverse panel of wheat using genome-wide association analysis

Leaf rust, caused by Puccinia triticina (Pt) is among the most devastating diseases posing a significant threat to global wheat production. The continuously evolving virulent Pt races in North America calls for exploring new sources of leaf rust resistance. A diversity panel of 365 bread wheat accessions selected from a worldwide population of landraces and cultivars was evaluated at the seedling stage against four Pt races (TDBJQ, TBBGS, MNPSD and, TNBJS). A wide distribution of seedling responses against the four Pt races was observed. Majority of the genotypes displayed a susceptible response with only 28 (9.8%), 59 (13.5%), 45 (12.5%), and 29 (8.1%) wheat accessions exhibiting a highly resistant response to TDBJQ, TBBGS, MNPSD and, TNBJS, respectively. Further, we conducted a high-resolution multi-locus genome-wide association study (GWAS) using a set of 302,524 high-quality single nucleotide polymorphisms (SNPs). The GWAS analysis identified 27 marker-trait associations (MTAs) for leaf rust resistance on different wheat chromosomes of which 20 MTAs were found in the vicinity of known Lr genes, MTAs, or quantitative traits loci (QTLs) identified in previous studies. The remaining seven significant MTAs identified represent genomic regions that harbor potentially novel genes for leaf rust resistance. Furthermore, the candidate gene analysis for the significant MTAs identified various genes of interest that may be involved in disease resistance. The identified resistant lines and SNPs linked to the QTLs in this study will serve as valuable resources in wheat rust resistance breeding programs.

Separation of the effects of two reduced height (Rht) genes and genomic background to select for less Fusarium head blight of short-strawed winter wheat (Triticum aestivum L.) varieties

FHB resistance shared pleiotropic loci with plant height and anther retention. Genomic prediction allows to select for genomic background reducing FHB susceptibility in the presence of the dwarfing allele Rht-D1b. With the high interest for semi-dwarf cultivars in wheat, finding locally adapted resistance sources against Fusarium head blight (FHB) and FHB-neutral reduced height (Rht) genes is of utmost relevance. In this study, 401 genotypes of European origin without/with dwarfing alleles of Rht-D1 and/or Rht24 were analysed across five environments on FHB severity and the morphological traits such as plant height (PH), anther retention (AR), number of spikelets per ear, ear length and ear density. Data were analysed by combined correlation and path analyses, association mapping and coupling single- and multi-trait genome-wide association studies (ST-GWAS and MT-GWAS, respectively) and genomic prediction (GP). All FHB data were corrected for flowering date or heading stage. High genotypic correlation (r_g = 0.74) and direct path effect (0.57) were detected between FHB severity and anther retention (AR). Moderate correlation (r_g = - 0.55) was found between FHB severity and plant height (PH) with a high indirect path via AR (- 0.31). Indirect selection for FHB resistance should concentrate on AR and PH. ST-GWAS identified 25 quantitative trait loci (QTL) for FHB severity, PH and AR, while MT-GWAS detected six QTL across chromosomes 2A, 4D, 5A, 6B and 7B conveying pleiotropic effects on the traits. Rht-D1b was associated with high AR and FHB susceptibility. Our study identified a promising positively acting pleiotropic QTL on chromosome 7B which can be utilized to improve FHB resistance while reducing PH and AR. Rht-D1b genotypes having a high resistance genomic background exhibited lower FHB severity and AR. The use of GP for estimating the genomic background was more effective than selection of GWAS-detected markers. We demonstrated that GP has a great potential and should be exploited by selecting for semi-dwarf winter wheat genotypes with higher FHB resistance due to their genomic background.

Consensus genomic regions associated with grain protein content in hexaploid and tetraploid wheat

A meta-analysis of QTLs associated with grain protein content (GPC) was conducted in hexaploid and tetraploid wheat to identify robust and stable meta-QTLs (MQTLs). For this purpose, as many as 459 GPC-related QTLs retrieved from 48 linkage-based QTL mapping studies were projected onto the newly developed wheat consensus map. The analysis resulted in the prediction of 57 MQTLs and 7 QTL hotspots located on all wheat chromosomes (except chromosomes 1D and 4D) and the average confidence interval reduced 2.71-fold in the MQTLs and QTL hotspots compared to the initial QTLs. The physical regions occupied by the MQTLs ranged from 140 bp to 224.02 Mb with an average of 15.2 Mb, whereas the physical regions occupied by QTL hotspots ranged from 1.81 Mb to 36.03 Mb with a mean of 8.82 Mb. Nineteen MQTLs and two QTL hotspots were also found to be co-localized with 45 significant SNPs identified in 16 previously published genome-wide association studies in wheat. Candidate gene (CG) investigation within some selected MQTLs led to the identification of 705 gene models which also included 96 high-confidence CGs showing significant expressions in different grain-related tissues and having probable roles in GPC regulation. These significantly expressed CGs mainly involved the genes/gene families encoding for the following proteins: aminotransferases, early nodulin 93, glutamine synthetases, invertase/pectin methylesterase inhibitors, protein BIG GRAIN 1-like, cytochrome P450, glycosyl transferases, hexokinases, small GTPases, UDP-glucuronosyl/UDP-glucosyltransferases, and EamA, SANT/Myb, GNAT, thioredoxin, phytocyanin, and homeobox domains containing proteins. Further, eight genes including GPC-B1, Glu-B1-1b, Glu-1By9, TaBiP1, GSr, TaNAC019-A, TaNAC019-D, and bZIP-TF SPA already known to be associated with GPC were also detected within some of the MQTL regions confirming the efficacy of MQTLs predicted during the current study.

A GBS-based genome-wide association study reveals the genetic basis of salinity tolerance at the seedling stage in bread wheat (Triticum aestivum L.)

High salinity levels affect 20% of the cultivated area and 9%–34% of the irrigated agricultural land worldwide, ultimately leading to yield losses of crops. The current study evaluated seven salt tolerance-related traits at the seedling stage in a set of 138 pre-breeding lines (PBLs) and identified 63 highly significant marker-trait associations (MTAs) linked to salt tolerance. Different candidate genes were identified in in silico analysis, many of which were involved in various stress conditions in plants, including glycine-rich cell wall structural protein 1-like, metacaspase-1, glyceraldehyde-3-phosphate dehydrogenase GAPA1, and plastidial GAPA1. Some of these genes coded for structural protein and participated in cell wall structure, some were linked to programmed cell death, and others were reported to show abiotic stress response roles in wheat and other plants. In addition, using the Multi-Trait Genotype-Ideotype Distance Index (MGIDI) protocol, the best-performing lines under salt stress were identified. The SNPs identified in this study and the genotypes with favorable alleles provide an excellent source to impart salt tolerance in wheat.

Genome-wide association study of coleoptile length with Shanxi wheat

In arid and semi-arid regions, coleoptile length is a vital agronomic trait for wheat breeding. The coleoptile length determines the maximum depth that seeds can be sown, and it is critical for establishment of the crop. Therefore, identifying loci associated with coleoptile length in wheat is essential. In the present study, 282 accessions from Shanxi Province representing wheat breeding for the Loess Plateau were grown under three experimental conditions to study coleoptile length. The results of phenotypic variation indicated that drought stress and light stress could lead to shortening of coleoptile length. Under drought stress the growth rate of environmentally sensitive cultivars decreased more than insensitive cultivars. The broad-sense heritability (H 2) of BLUP (best linear unbiased prediction) under various conditions showed G × E interaction for coleoptile length but was mainly influenced by heredity. Correlation analysis showed that correlation between plant height-related traits and coleoptile length was significant in modern cultivars whereas it was not significant in landraces. A total of 45 significant marker-trait associations (MTAs) for coleoptile length in the three conditions were identified using the 3VmrMLM (3 Variance-component multi-locus random-SNP-effect Mixed Linear Model) and MLM (mixed linear model). In total, nine stable genetic loci were identified via 3VmrMLM under the three conditions, explaining 2.94-7.79% of phenotypic variation. Five loci on chromosome 2B, 3A, 3B, and 5B have not been reported previously. Six loci had additive effects toward increasing coleoptile length, three of which are novel. Molecular markers for the loci with additive effects on coleoptile length can be used to breed cultivars with long coleoptiles.

Whole-genome analysis of hard winter wheat germplasm identifies genomic regions associated with spike and kernel traits

Genetic dissection of yield component traits including spike and kernel characteristics is essential for the continuous improvement in wheat yield. Genome-wide association studies (GWAS) have been frequently used to identify genetic determinants for spike and kernel-related traits in wheat, though none have been employed in hard winter wheat (HWW) which represents a major class in US wheat acreage. Further, most of these studies relied on assembled diversity panels instead of adapted breeding lines, limiting the transferability of results to practical wheat breeding. Here we assembled a population of advanced/elite breeding lines and well-adapted cultivars and evaluated over four environments for phenotypic analysis of spike and kernel traits. GWAS identified 17 significant multi-environment marker-trait associations (MTAs) for various traits, representing 12 putative quantitative trait loci (QTLs), with five QTLs affecting multiple traits. Four of these QTLs mapped on three chromosomes 1A, 5B, and 7A for spike length, number of spikelets per spike (NSPS), and kernel length are likely novel. Further, a highly significant QTL was detected on chromosome 7AS that has not been previously associated with NSPS and putative candidate genes were identified in this region. The allelic frequencies of important quantitative trait nucleotides (QTNs) were deduced in a larger set of 1,124 accessions which revealed the importance of identified MTAs in the US HWW breeding programs. The results from this study could be directly used by the breeders to select the lines with favorable alleles for making crosses, and reported markers will facilitate marker-assisted selection of stable QTLs for yield components in wheat breeding.

Delineating meta-quantitative trait loci for anthracnose resistance in common bean (Phaseolus vulgaris L.)

Anthracnose, caused by the fungus Colletotrichum lindemuthianum, is one of the devastating disease affecting common bean production and productivity worldwide. Several quantitative trait loci (QTLs) for anthracnose resistance have been identified. In order to make use of these QTLs in common bean breeding programs, a detailed meta-QTL (MQTL) analysis has been conducted. For the MQTL analysis, 92 QTLs related to anthracnose disease reported in 18 different earlier studies involving 16 mapping populations were compiled and projected on to the consensus map. This meta-analysis led to the identification of 11 MQTLs (each involving QTLs from at least two different studies) on 06 bean chromosomes and 10 QTL hotspots each involving multiple QTLs from an individual study on 07 chromosomes. The confidence interval (CI) of the identified MQTLs was found 3.51 times lower than the CI of initial QTLs. Marker-trait associations (MTAs) reported in published genome-wide association studies (GWAS) were used to validate nine of the 11 identified MQTLs, with MQTL4.1 overlapping with as many as 40 MTAs. Functional annotation of the 11 MQTL regions revealed 1,251 genes including several R genes (such as those encoding for NBS-LRR domain-containing proteins, protein kinases, etc.) and other defense related genes. The MQTLs, QTL hotspots and the potential candidate genes identified during the present study will prove useful in common bean marker-assisted breeding programs and in basic studies involving fine mapping and cloning of genomic regions associated with anthracnose resistance in common beans.

Consensus genomic regions associated with multiple abiotic stress tolerance in wheat and implications for wheat breeding

In wheat, a meta-analysis was performed using previously identified QTLs associated with drought stress (DS), heat stress (HS), salinity stress (SS), water-logging stress (WS), pre-harvest sprouting (PHS), and aluminium stress (AS) which predicted a total of 134 meta-QTLs (MQTLs) that involved at least 28 consistent and stable MQTLs conferring tolerance to five or all six abiotic stresses under study. Seventy-six MQTLs out of the 132 physically anchored MQTLs were also verified with genome-wide association studies. Around 43% of MQTLs had genetic and physical confidence intervals of less than 1 cM and 5 Mb, respectively. Consequently, 539 genes were identified in some selected MQTLs providing tolerance to 5 or all 6 abiotic stresses. Comparative analysis of genes underlying MQTLs with four RNA-seq based transcriptomic datasets unravelled a total of 189 differentially expressed genes which also included at least 11 most promising candidate genes common among different datasets. The promoter analysis showed that the promoters of these genes include many stress responsiveness cis-regulatory elements, such as ARE, MBS, TC-rich repeats, As-1 element, STRE, LTR, WRE3, and WUN-motif among others. Further, some MQTLs also overlapped with as many as 34 known abiotic stress tolerance genes. In addition, numerous ortho-MQTLs among the wheat, maize, and rice genomes were discovered. These findings could help with fine mapping and gene cloning, as well as marker-assisted breeding for multiple abiotic stress tolerances in wheat.

Consensus genomic regions associated with multiple abiotic stress tolerance in wheat and implications for wheat breeding

In wheat, a meta-analysis was performed using previously identified QTLs associated with drought stress (DS), heat stress (HS), salinity stress (SS), water-logging stress (WS), pre-harvest sprouting (PHS), and aluminium stress (AS) which predicted a total of 134 meta-QTLs (MQTLs) that involved at least 28 consistent and stable MQTLs conferring tolerance to five or all six abiotic stresses under study. Seventy-six MQTLs out of the 132 physically anchored MQTLs were also verified with genome-wide association studies. Around 43% of MQTLs had genetic and physical confidence intervals of less than 1 cM and 5 Mb, respectively. Consequently, 539 genes were identified in some selected MQTLs providing tolerance to 5 or all 6 abiotic stresses. Comparative analysis of genes underlying MQTLs with four RNA-seq based transcriptomic datasets unravelled a total of 189 differentially expressed genes which also included at least 11 most promising candidate genes common among different datasets. The promoter analysis showed that the promoters of these genes include many stress responsiveness cis-regulatory elements, such as ARE, MBS, TC-rich repeats, As-1 element, STRE, LTR, WRE3, and WUN-motif among others. Further, some MQTLs also overlapped with as many as 34 known abiotic stress tolerance genes. In addition, numerous ortho-MQTLs among the wheat, maize, and rice genomes were discovered. These findings could help with fine mapping and gene cloning, as well as marker-assisted breeding for multiple abiotic stress tolerances in wheat.

Multi-Locus Genome-Wide Association Studies to Characterize Fusarium Head Blight (FHB) Resistance in Hard Winter Wheat

Fusarium head blight (FHB), caused by the fungus Fusarium graminearum Schwabe is an important disease of wheat that causes severe yield losses along with serious quality concerns. Incorporating the host resistance from either wild relatives, landraces, or exotic materials remains challenging and has shown limited success. Therefore, a better understanding of the genetic basis of native FHB resistance in hard winter wheat (HWW) and combining it with major quantitative trait loci (QTLs) can facilitate the development of FHB-resistant cultivars. In this study, we evaluated a set of 257 breeding lines from the South Dakota State University (SDSU) breeding program to uncover the genetic basis of native FHB resistance in the US hard winter wheat. We conducted a multi-locus genome-wide association study (ML-GWAS) with 9,321 high-quality single-nucleotide polymorphisms (SNPs). A total of six distinct marker-trait associations (MTAs) were identified for the FHB disease index (DIS) on five different chromosomes including 2A, 2B, 3B, 4B, and 7A. Further, eight MTAs were identified for Fusarium-damaged kernels (FDK) on six chromosomes including 3B, 5A, 6B, 6D, 7A, and 7B. Out of the 14 significant MTAs, 10 were found in the proximity of previously reported regions for FHB resistance in different wheat classes and were validated in HWW, while four MTAs represent likely novel loci for FHB resistance. Accumulation of favorable alleles of reported MTAs resulted in significantly lower mean DIS and FDK score, demonstrating the additive effect of FHB resistance alleles. Candidate gene analysis for two important MTAs identified several genes with putative proteins of interest; however, further investigation of these regions is needed to identify genes conferring FHB resistance. The current study sheds light on the genetic basis of native FHB resistance in the US HWW germplasm and the resistant lines and MTAs identified in this study will be useful resources for FHB resistance breeding via marker-assisted selection.

Germination of Triticum aestivum L.: Effects of Soil–Seed Interaction on the Growth of Seedlings

Seed size, sowing depth, and seed disinfection can affect seed germination and seedling establishment, which, in turn, can directly affect crop growth and yield. The current study was comprised of two experiments, the first of which was conducted in the laboratory, and a second which was performed under glasshouse conditions. The objective of these experiments was to investigate the effects of seed size, sowing depth, and seed disinfection on seed germination and initial seedling growth of selected wheat (Triticum aestivum L.) cultivars. The treatments in laboratory experiment were arranged in a completely randomized design, which included: (Ι) four wheat cultivars (Pishgam, Haydari, Soissons, and Mihan), (ΙΙ) two seed size classes (x < 2.25 mm, and x > 2.25 mm), and two disinfection treatments (no-disinfection and disinfection), (ΙΙΙ) with five replicates. In addition to the aforementioned treatments, the effect of planting depth (4, 6, and 8 cm) was also investigated in the subsequent glasshouse experiment. The best results were obtained at a sowing depth of 4 cm, in the non-disinfected treatment, using large seeds. In contrast, the lowest percentage and speed of seed germination and vigor index were observed in seeds sown at 8 cm depth, in the disinfected seed treatment, using small seeds. Large seeds contain larger nutrient stores which may improve seed germination indices, which would therefore result in improved percentage and speed of seed germination, followed by faster coleoptile and seedling growth, higher seedling dry weight and seed vigor. These data also illustrated that seed disinfection in the Pishgam and Haydari cultivars had inhibitory effects upon coleoptile growth and seedling length, which could be related to the fungicide’s chemical composition. Unlike other cultivars, disinfection did not show a significant effect on the Soissons cultivar. Based on our data, in order to improve both the speed of wheat seed germination and subsequent plant growth and development; it is necessary to select high-quality, large seeds, planted at a specific planting depth, which have been treated with an effective disinfectant; all of which will be specific for the wheat cultivar in question. Overall, the current study has provided useful information on the effect size seed, sowing depth, and disinfection have upon germination characteristics and seedling growth of wheat cultivars, which can form the basis for future field scale trails.

Comprehensive evaluation of mapping complex traits in wheat using genome-wide association studies

Genome-wide association studies (GWAS) are effectively applied to detect the marker trait associations (MTAs) using whole genome-wide variants for complex quantitative traits in different crop species. GWAS has been applied in wheat for different quality, biotic and abiotic stresses, and agronomic and yield-related traits. Predictions for marker-trait associations are controlled with the development of better statistical models taking population structure and familial relatedness into account. In this review, we have provided a detailed overview of the importance of association mapping, population design, high-throughput genotyping and phenotyping platforms, advancements in statistical models and multiple threshold comparisons, and recent GWA studies conducted in wheat. The information about MTAs utilized for gene characterization and adopted in breeding programs is also provided. In the literature that we surveyed, as many as 86,122 wheat lines have been studied under various GWA studies reporting 46,940 loci. However, further utilization of these is largely limited. The future breakthroughs in area of genomic selection, multi-omics-based approaches, machine, and deep learning models in wheat breeding after exploring the complex genetic structure with the GWAS are also discussed. This is a most comprehensive study of a large number of reports on wheat GWAS and gives a comparison and timeline of technological developments in this area. This will be useful to new researchers or groups who wish to invest in GWAS.

Phenotypic Evaluation and Genetic Analysis of Seedling Emergence in a Global Collection of Wheat Genotypes (Triticum aestivum L.) Under Limited Water Availability

The challenge in establishing an early-sown wheat crop in southern Australia is the need for consistently high seedling emergence when sowing deep in subsoil moisture (>10 cm) or into dry top-soil (4 cm). However, the latter is strongly reliant on a minimum soil water availability to ensure successful seedling emergence. This study aimed to: (1) evaluate 233 Australian and selected international wheat genotypes for consistently high seedling emergence under limited soil water availability when sown in 4 cm of top-soil in field and glasshouse (GH) studies; (2) ascertain genetic loci associated with phenotypic variation using a genome-wide association study (GWAS); and (3) compare across loci for traits controlling coleoptile characteristics, germination, dormancy, and pre-harvest sprouting. Despite significant (P < 0.001) environment and genotype-by-environment interactions within and between field and GH experiments, eight genotypes that included five cultivars, two landraces, and one inbred line had consistently high seedling emergence (mean value > 85%) across nine environments. Moreover, 21 environment-specific quantitative trait loci (QTL) were detected in GWAS analysis on chromosomes 1B, 1D, 2B, 3A, 3B, 4A, 4B, 5B, 5D, and 7D, indicating complex genetic inheritance controlling seedling emergence. We aligned QTL for known traits and individual genes onto the reference genome of wheat and identified 16 QTL for seedling emergence in linkage disequilibrium with coleoptile length, width, and cross-sectional area, pre-harvest sprouting and dormancy, germination, seed longevity, and anthocyanin development. Therefore, it appears that seedling emergence is controlled by multifaceted networks of interrelated genes and traits regulated by different environmental cues.

Genome-wide association analysis permits characterization of Stagonospora nodorum blotch (SNB) resistance in hard winter wheat

Stagonospora nodorum blotch (SNB) is an economically important wheat disease caused by the necrotrophic fungus Parastagonospora nodorum. SNB resistance in wheat is controlled by several quantitative trait loci (QTLs). Thus, identifying novel resistance/susceptibility QTLs is crucial for continuous improvement of the SNB resistance. Here, the hard winter wheat association mapping panel (HWWAMP) comprising accessions from breeding programs in the Great Plains region of the US, was evaluated for SNB resistance and necrotrophic effectors (NEs) sensitivity at the seedling stage. A genome-wide association study (GWAS) was performed to identify single‐nucleotide polymorphism (SNP) markers associated with SNB resistance and effectors sensitivity. We found seven significant associations for SNB resistance/susceptibility distributed over chromosomes 1B, 2AL, 2DS, 4AL, 5BL, 6BS, and 7AL. Two new QTLs for SNB resistance/susceptibility at the seedling stage were identified on chromosomes 6BS and 7AL, whereas five QTLs previously reported in diverse germplasms were validated. Allele stacking analysis at seven QTLs explained the additive and complex nature of SNB resistance. We identified accessions (‘Pioneer-2180’ and ‘Shocker’) with favorable alleles at five of the seven identified loci, exhibiting a high level of resistance against SNB. Further, GWAS for sensitivity to NEs uncovered significant associations for SnToxA and SnTox3, co-locating with previously identified host sensitivity genes (Tsn1 and Snn3). Candidate region analysis for SNB resistance revealed 35 genes of putative interest with plant defense response-related functions. The QTLs identified and validated in this study could be easily employed in breeding programs using the associated markers to enhance the SNB resistance in hard winter wheat.

Position Validation of the Dwarfing Gene Dw6 in Oat (Avena sativa L.) and Its Correlated Effects on Agronomic Traits

An F6 : 8 recombinant inbred line (RIL) population derived from the cross between WAOAT2132 (Dw6) and Caracas along with the two parents were used to evaluate the genetic effects of Dw6 dwarfing gene on plant height and other agronomic traits in oat (Avena sativa L.) across three environments, and develop closely linked markers for marker-assisted selection (MAS) for Dw6. The two parents differed in all investigated agronomic traits except for the number of whorls. The RIL lines showed a bimodal distribution for plant height in all three tested environments, supporting the height of this population was controlled by a single gene. Dw6 significantly reduced plant height (37.66∼44.29%) and panicle length (13.99∼22.10%) but without compromising the coleoptile length which was often positively associated with the reduced stature caused by dwarfing genes. Dw6 has also strong negative effects on hundred kernel weight (14.00∼29.55%), and kernel length (4.21∼9.47%), whereas the effects of Dw6 on the kernel width were not uniform across three environments. By contrast, lines with Dw6 produced more productive tillers (10.11∼10.53%) than lines without Dw6. All these together suggested the potential yield penalty associated with Dw6 might be partially due to the decrease of kernel weight which is attributed largely to the reduction of kernel length. Eighty-one simple sequence repeat (SSR) primer pairs from chromosome 6D were tested, five of them were polymorphic in two parents and in two contrasting bulks, confirming the 6D location of Dw6. By using the five polymorphic markers, Dw6 was mapped to an interval of 1.0 cM flanked by markers SSR83 and SSR120. Caution should be applied in using this information since maker order conflicts were observed. The close linkages of these two markers to Dw6 were further validated in a range of oat lines. The newly developed markers will provide a solid basis for future efforts both in the identification of Dw6 in oat germplasm and in the determination of the nature of the gene through positional cloning.

Insights into the Genetic Architecture of Bran Friability and Water Retention Capacity, Two Important Traits for Whole Grain End-Use Quality in Winter Wheat

Bran friability (particle size distribution after milling) and water retention capacity (WRC) impact wheat bran functionality in whole grain milling and baking applications. The goal of this study was to identify genomic regions and underlying genes that may be responsible for these traits. The Hard Winter Wheat Association Mapping Panel, which comprised 299 lines from breeding programs in the Great Plains region of the US, was used in a genome-wide association study. Bran friability ranged from 34.5% to 65.9% (median, 51.1%) and WRC ranged from 159% to 458% (median, 331%). Two single-nucleotide polymorphisms (SNPs) on chromosome 5D were significantly associated with bran friability, accounting for 11-12% of the phenotypic variation. One of these SNPs was located within the Puroindoline-b gene, which is known for influencing endosperm texture. Two SNPs on chromosome 4A were tentatively associated with WRC, accounting for 4.6% and 4.4% of phenotypic variation. The favorable alleles at the SNP sites were present in only 15% (friability) and 34% (WRC) of lines, indicating a need to develop new germplasm for these whole-grain end-use quality traits. Validation of these findings in independent populations will be useful for breeding winter wheat cultivars with improved functionality for whole grain food applications.

Phenotypic and genetic variation in phosphorus-deficiency-tolerance traits in Chinese wheat landraces

BACKGROUND

Phosphorus deficiency is a major limiting factors for affecting crop production globally. To understand the genetic variation of phosphorus-deficiency-tolerance, a total of 15 seedling traits were evaluated among 707 Chinese wheat landraces under application of phosphorus (AP) and non-application of phosphorus (NP). A total of 18,594 single-nucleotide polymorphisms and 38,678 diversity arrays technology sequencing markers were used to detect marker-trait associations under AP and NP.

RESULTS

Top ten genotypes with extremely tolerance and bottommost ten genotypes with extremely sensitivity were selected from 707 Chinese wheat landraces for future breeding and genetic analysis. A total of 55 significant markers (81 marker-trait associations) for 13 traits by both CMLM and SUPER method. These were distributed on chromosomes 1A, 1B, 2A, 2B, 2D, 3A, 4B, 5A, 5B, 6A, 6B, 6D, 7A and 7B. Considering the linkage disequilibrium decay distance, 25 and 12 quantitative trait loci (QTL) were detected under AP and NP, respectively (9 QTL were specific to NP).

CONCLUSIONS

The extremely tolerant landraces could be used for breeding phosphorus-deficiency-tolerant cultivars. The QTL could be useful in wheat breeding through marker-assisted selection. Our findings provide new insight into the genetic analysis of P-deficiency-tolerance, and will be helpful for breeding P-deficiency-tolerant cultivars.