Genome-wide association analyses of 54 traits identified multiple loci for the determination of floret fertility in wheat

Genetic response of a perennial grass to warm and wet environments interacts and is associated with trait means as well as plasticity

The potential for rapid evolution is an important mechanism allowing species to adapt to changing climatic conditions. Although such potential has been largely studied in various short-lived organisms, to what extent we can observe similar patterns in long-lived plant species, which often dominate natural systems, is largely unexplored. We explored the potential for rapid evolution in Festuca rubra, a long-lived grass with extensive clonal growth dominating in alpine grasslands. We used a field sowing experiment simulating expected climate change in our model region. Specifically, we exposed seeds from five independent seed sources to novel climatic conditions by shifting them along a natural climatic grid and explored the genetic profiles of established seedlings after 3 years. Data on genetic profiles of plants selected under different novel conditions indicate that different climate shifts select significantly different pools of genotypes from common seed pools. Increasing soil moisture was more important than increasing temperature or the interaction of the two climatic factors in selecting pressure. This can indicate negative genetic interaction in response to the combined effects or that the effects of different climates are interactive rather than additive. The selected alleles were found in genomic regions, likely affecting the function of specific genes or their expression. Many of these were also linked to morphological traits (mainly to trait plasticity), suggesting these changes may have a consequence on plant performance. Overall, these data indicate that even long-lived plant species may experience strong selection by climate, and their populations thus have the potential to rapidly adapt to these novel conditions.

Dissecting the molecular basis of spike traits by integrating gene regulatory networks and genetic variation in wheat

Spike architecture influences both grain weight and grain number per spike, which are the two major components of grain yield in bread wheat (Triticum aestivum L.). However, the complex wheat genome and the influence of various environmental factors pose challenges in mapping the causal genes that affect spike traits. Here, we systematically identified genes involved in spike trait formation by integrating information on genomic variation and gene regulatory networks controlling young spike development in wheat. We identified 170 loci that are responsible for variations in spike length, spikelet number per spike, and grain number per spike through genome-wide association study and meta-QTL analyses. We constructed gene regulatory networks for young inflorescences at the double ridge stage and the floret primordium stage, in which the spikelet meristem and the floret meristem are predominant, respectively, by integrating transcriptome, histone modification, chromatin accessibility, eQTL, and protein-protein interactome data. From these networks, we identified 169 hub genes located in 76 of the 170 QTL regions whose polymorphisms are significantly associated with variation in spike traits. The functions of TaZF-B1, VRT-B2, and TaSPL15-A/D in establishment of wheat spike architecture were verified. This study provides valuable molecular resources for understanding spike traits and demonstrates that combining genetic analysis and developmental regulatory networks is a robust approach for dissection of complex traits.

Novel PHOTOPERIOD-1 gene variants associate with yield-related and root-angle traits in European bread wheat

KEY MESSAGE

PHOTOPERIOD-1 homoeologous gene copies play a pivotal role in regulation of flowering time in wheat. Here, we show that their influence also extends to spike and shoot architecture and even impacts root development. The sequence diversity of three homoeologous copies of the PHOTOPERIOD-1 gene in European winter wheat was analyzed by Oxford Nanopore amplicon-based multiplex sequencing and molecular markers in a panel of 194 cultivars representing breeding progress over the past 5 decades. A strong, consistent association with an average 8% increase in grain yield was observed for the PpdA1-Hap1 haplotype across multiple environments. This haplotype was found to be linked in 51% of cultivars to the 2NS/2AS translocation, originally introduced from Aegilops ventricosa, which leads to an overestimation of its effect. However, even in cultivars without the 2NS/2AS translocation, PpdA1-Hap1 was significantly associated with increased grain yield, kernel per spike and kernel per m2 under optimal growth conditions, conferring a 4% yield advantage compared to haplotype PpdA1-Hap4. In contrast to Ppd-B1 and Ppd-D1, the Ppd-A1 gene exhibits novel structural variations and a high number of SNPs, highlighting the evolutionary changes that have occurred in this region over the course of wheat breeding history. Additionally, cultivars carrying the photoperiod-insensitive Ppd-D1a allele not only exhibit earlier heading, but also deeper roots compared to those with photoperiod-sensitive alleles under German conditions. PCR and KASP assays have been developed that can be effectively employed in marker-assisted breeding programs to introduce these favorable haplotypes.

Mining of QTLs for Spring Bread Wheat Spike Productivity by Comparing Spring Wheat Cultivars Released in Different Decades of the Last Century

Genome-wide association studies (GWAS) are among the genetic tools for the mining of genomic loci associated with useful agronomic traits. The study enabled us to find new genetic markers associated with grain yield as well as quality. The sample under study consisted of spring wheat cultivars developed in different decades of the last century. A panel of 186 accessions was evaluated at VIR's experiment station in Pushkin across a 3-year period of field trials. In total, 24 SNPs associated with six productivity characteristics were revealed. Along with detecting significant markers for each year of the field study, meta-analyses were conducted. Loci associated with useful yield-related agronomic characteristics were detected on chromosomes 4A, 5A, 6A, 6B, and 7B. In addition to previously described regions, novel loci associated with grain yield and quality were identified during the study. We presume that the utilization of contrast cultivars which originated in different breeding periods allowed us to identify new markers associated with useful agronomic characteristics.

A co-located QTL for seven spike architecture-related traits shows promising breeding use potential in common wheat (Triticum aestivum L.)

KEY MESSAGE

A co-located novel QTL for TFS, FPs, FMs, FFS, FFPs, KWS, and KWPs with potential of improving wheat yield was identified and validated. Spike-related traits, including fertile florets per spike (FFS), kernel weight per spike (KWS), total florets per spike (TFS), florets per spikelet (FPs), florets in the middle spikelet (FMs), fertile florets per spikelet (FFPs), and kernel weight per spikelet (KWPs), are key traits in improving wheat yield. In the present study, quantitative trait loci (QTL) for these traits evaluated under various environments were detected in a recombinant inbred line population (msf/Chuannong 16) mainly genotyped using the 16 K SNP array. Ultimately, we identified 60 QTL, but only QFFS.sau-MC-1A for FFS was a major and stably expressed QTL. It was located on chromosome arm 1AS, where loci for TFS, FPs, FMs, FFS, FFPs, KWS, and KWPs were also simultaneously co-mapped. The effect of QFFS.sau-MC-1A was further validated in three independent segregating populations using a Kompetitive Allele-Specific PCR marker. For the co-located QTL, QFFS.sau-MC-1A, the presence of a positive allele from msf was associate with increases for all traits: + 12.29% TFS, + 10.15% FPs, + 13.97% FMs, + 17.12% FFS, + 14.75% FFPs, + 22.17% KWS, and + 19.42% KWPs. Furthermore, pleiotropy analysis showed that the positive allele at QFFS.sau-MC-1A simultaneously increased the spike length, spikelet number per spike, and thousand-kernel weight. QFFS.sau-MC-1A represents a novel QTL for marker-assisted selection with the potential for improving wheat yield. Four genes, TraesCS1A03G0012700, TraesCS1A03G0015700, TraesCS1A03G0016000, and TraesCS1A03G0016300, which may affect spike development, were predicted in the physical interval harboring QFFS.sau-MC-1A. Our results will help in further fine mapping QFFS.sau-MC-1A and be useful for improving wheat yield.

Grain yield trade-offs in spike-branching wheat can be mitigated by elite alleles affecting sink capacity and post-anthesis source activity

Introducing variations in inflorescence architecture, such as the 'Miracle-Wheat' (Triticum turgidum convar. compositum (L.f.) Filat.) with a branching spike, has relevance for enhancing wheat grain yield. However, in the spike-branching genotypes, the increase in spikelet number is generally not translated into grain yield advantage because of reduced grains per spikelet and grain weight. Here, we investigated if such trade-offs might be a function of source-sink strength by using 385 recombinant inbred lines developed by intercrossing the spike-branching landrace TRI 984 and CIRNO C2008, an elite durum (T. durum L.) cultivar; they were genotyped using the 25K array. Various plant and spike architectural traits, including flag leaf, peduncle, and spike senescence rate, were phenotyped under field conditions for 2 consecutive years. On chromosome 5AL, we found a new modifier QTL for spike branching, branched headt3 (bht-A3), which was epistatic to the previously known bht-A1 locus. Besides, bht-A3 was associated with more grains per spikelet and a delay in flag leaf senescence rate. Importantly, favourable alleles, viz. bht-A3 and grain protein content (gpc-B1) that delayed senescence, are required to improve grain number and grain weight in the spike-branching genotypes. In summary, achieving a balanced source-sink relationship might minimize grain yield trade-offs in Miracle-Wheat.

Ameliorating the effects of multiple stresses on agronomic traits in crops: modern biotechnological and omics approaches

While global climate change poses a significant environmental threat to agriculture, the increasing population is another big challenge to food security. To address this, developing crop varieties with increased productivity and tolerance to biotic and abiotic stresses is crucial. Breeders must identify traits to ensure higher and consistent yields under inconsistent environmental challenges, possess resilience against emerging biotic and abiotic stresses and satisfy customer demands for safer and more nutritious meals. With the advent of omics-based technologies, molecular tools are now integrated with breeding to understand the molecular genetics of genotype-based traits and develop better climate-smart crops. The rapid development of omics technologies offers an opportunity to generate novel datasets for crop species. Identifying genes and pathways responsible for significant agronomic traits has been made possible by integrating omics data with genetic and phenotypic information. This paper discusses the importance and use of omics-based strategies, including genomics, transcriptomics, proteomics and phenomics, for agricultural and horticultural crop improvement, which aligns with developing better adaptability in these crop species to the changing climate conditions.

Image facilitated assessment of intra-spike variation in grain size in wheat under high temperature and drought stress

In wheat (Triticum aestivum L.), the grain size varies according to position within the spike. Exposure to drought and high temperature stress during grain development in wheat reduces grain size, and this reduction also varies across the length of the spike. We developed the phenomics approach involving image-based tools to assess the intra-spike variation in grain size. The grains were arranged corresponding to the spikelet position and the camera of smart phone was used to acquire 333 images. The open-source software ImageJ was used to analyze features of each grain and the image-derived parameters were used to calculate intra-spike variation as standard deviation (ISVAD). The effect of genotype and environment were highly significant on the ISVAD of grain area. Sunstar and Raj 4079 contrasted in the ISVAD of grain area under late sown environment, and RNA sequencing of the spike was done at 25 days after anthesis. The genes for carbohydrate transport and stress response were upregulated in Sunstar as compared to Raj 4079, suggesting that these play a role in intra-spike assimilate distribution. The phenomics method developed may be useful for grain phenotyping and identifying germplasm with low intra-spike variation in grain size for their further validation as parental material in breeding.

High-affinity iron uptake is required for optimal Epichloë festucae colonization of Lolium perenne and seed transmission

Epichloë festucae uses a siderophore-mediated system to acquire iron, which is important to maintain endophyte-grass symbioses. Here we investigate the roles of the alternative iron acquisition system, reductive iron assimilation (RIA), via disruption of the fetC gene, which encodes a multicopper ferroxidase, either alone (i.e., ΔfetC) or in combination with disruption of the gene sidA, which encodes a siderophore biosynthesis enzyme (i.e., ΔfetC/ΔsidA). The phenotypic characteristics of these mutants were compared to ΔsidA and wild-type (WT) strains during growth under axenic culture conditions (in culture) and in symbiosis with the host grass, perennial ryegrass (in planta). Under iron deficiency, the colony growth rate of ΔfetC was slightly slower than that of WT, while the growth of ΔsidA and ΔfetC/ΔsidA mutants was severely suppressed. Siderophore analyses indicated that ΔfetC mutants hyperaccumulate ferriepichloënin A (FEA) at low iron concentrations and ferricrocin and FEA at higher iron concentrations. When compared to WT, all mutant strains displayed hyperbranching hyphal structures and a reduced ratio of Epichloë DNA to total DNA in planta. Furthermore, host colonization and vertical transmission through infection of the host seed were significantly reduced in the ΔfetC/ΔsidA mutants, confirming that high-affinity iron uptake is a critical process for Epichloë transmission. Thus, RIA and siderophore iron uptake are complementary systems required for the maintenance of iron metabolism, fungal growth, and symbiosis between E. festucae and perennial ryegrass.

Characterization and fine mapping analysis of a major stable QTL qKnps-4A for kernel number per spike in wheat

KEY MESSAGE

A major stable QTL for kernel number per spike was narrowed down to a 2.19-Mb region containing two potential candidate genes, and its effects on yield-related traits were characterized. Kernel number per spike (KNPS) in wheat is a key yield component. Dissection and characterization of major stable quantitative trait loci (QTLs) for KNPS would be of considerable value for the genetic improvement of yield potential using molecular breeding technology. We had previously reported a major stable QTL controlling KNPS, qKnps-4A. In the current study, primary fine-mapping analysis, based on the primary mapping population, located qKnps-4A to an interval of approximately 6.8-Mb from 649.0 to 655.8 Mb on chromosome 4A refering to 'Kenong 9204' genome. Further fine-mapping analysis based on a secondary mapping population narrowed qKnps-4A to an approximately 2.19-Mb interval from 653.72 to 655.91 Mb. Transcriptome sequencing, gene function annotation analysis and homologous gene related reports showed that TraesKN4A01HG38570 and TraesKN4A01HG38590 were most likely to be candidate genes of qKnps-4A. Phenotypic analysis based on paired near-isogenic lines in the target region showed that qKnps-4A increased KNPS mainly by increasing the number of central florets per spike. We also evaluated the effects of qKnps-4A on other yield-related traits. Moreover, we dissected the QTL cluster of qKnps-4A and qTkw-4A and proved that the phenotypic effects were probably due to close linkage of two or more genes rather than pleiotropic effects of a single gene. This study provides molecular marker resource for wheat molecular breeding designed to improve yield potential, and lay the foundation for gene functional analysis of qKnps-4A.

A high-resolution genotype-phenotype map identifies the TaSPL17 controlling grain number and size in wheat

BACKGROUND

Large-scale genotype-phenotype association studies of crop germplasm are important for identifying alleles associated with favorable traits. The limited number of single-nucleotide polymorphisms (SNPs) in most wheat genome-wide association studies (GWASs) restricts their power to detect marker-trait associations. Additionally, only a few genes regulating grain number per spikelet have been reported due to sensitivity of this trait to variable environments.

RESULTS

We perform a large-scale GWAS using approximately 40 million filtered SNPs for 27 spike morphology traits. We detect 132,086 significant marker-trait associations and the associated SNP markers are located within 590 associated peaks. We detect additional and stronger peaks by dividing spike morphology into sub-traits relative to GWAS results of spike morphology traits. We propose that the genetic dissection of spike morphology is a powerful strategy to detect signals for grain yield traits in wheat. The GWAS results reveal that TaSPL17 positively controls grain size and number by regulating spikelet and floret meristem development, which in turn leads to enhanced grain yield per plant. The haplotypes at TaSPL17 indicate geographical differentiation, domestication effects, and breeding selection.

CONCLUSION

Our study provides valuable resources for genetic improvement of spike morphology and a fast-forward genetic solution for candidate gene detection and cloning in wheat.

Genomics for Yield and Yield Components in Durum Wheat

In recent years, many efforts have been conducted to dissect the genetic basis of yield and yield components in durum wheat thanks to linkage mapping and genome-wide association studies. In this review, starting from the analysis of the genetic bases that regulate the expression of yield for developing new durum wheat varieties, we have highlighted how, currently, the reductionist approach, i.e., dissecting the yield into its individual components, does not seem capable of ensuring significant yield increases due to diminishing resources, land loss, and ongoing climate change. However, despite the identification of genes and/or chromosomal regions, controlling the grain yield in durum wheat is still a challenge, mainly due to the polyploidy level of this species. In the review, we underline that the next-generation sequencing (NGS) technologies coupled with improved wheat genome assembly and high-throughput genotyping platforms, as well as genome editing technology, will revolutionize plant breeding by providing a great opportunity to capture genetic variation that can be used in breeding programs. To date, genomic selection provides a valuable tool for modeling optimal allelic combinations across the whole genome that maximize the phenotypic potential of an individual under a given environment.

Genetic basis of geographical differentiation and breeding selection for wheat plant architecture traits

BACKGROUND

Plant architecture associated with increased grain yield and adaptation to the local environments is selected during wheat (Triticum aestivum) breeding. The internode length of individual stems and tiller length of individual plants are important for the determination of plant architecture. However, few studies have explored the genetic basis of these traits.

RESULTS

Here, we conduct a genome-wide association study (GWAS) to dissect the genetic basis of geographical differentiation of these traits in 306 worldwide wheat accessions including both landraces and traditional varieties. We determine the changes of haplotypes for the associated genomic regions in frequency in 831 wheat accessions that are either introduced from other countries or developed in China from last two decades. We identify 83 loci that are associated with one trait, while the remaining 247 loci are pleiotropic. We also find 163 associated loci are under strong selective sweep. GWAS results demonstrate independent regulation of internode length of individual stems and consistent regulation of tiller length of individual plants. This makes it possible to obtain ideal haplotype combinations of the length of four internodes. We also find that the geographical distribution of the haplotypes explains the observed differences in internode length among the worldwide wheat accessions.

CONCLUSION

This study provides insights into the genetic basis of plant architecture. It will facilitate gene functional analysis and molecular design of plant architecture for breeding.

The pathway of melatonin biosynthesis in common wheat (Triticum aestivum)

Melatonin (Mel) is a multifunctional biomolecule found in both animals and plants. In plants, the biosynthesis of Mel from tryptophan (Trp) has been delineated to comprise of four consecutive reactions. However, while the genes encoding these enzymes in rice are well characterized no systematic evaluation of the overall pathway has, as yet, been published for wheat. In the current study, the relative contents of six Mel-pathway-intermediates including Trp, tryptamine (Trm), serotonin (Ser), 5-methoxy tryptamine (5M-Trm), N-acetyl serotonin (NAS) and Mel, were determined in 24 independent tissues spanning the lifetime of wheat. These studies indicated that Trp was the most abundant among the six metabolites, followed by Trm and Ser. Next, the candidate genes expressing key enzymes involved in the Mel pathway were explored by means of metabolite-based genome-wide association study (mGWAS), wherein two TDC genes, a T5H gene and one SNAT gene were identified as being important for the accumulation of Mel pathway metabolites. Moreover, a 463-bp insertion within the T5H gene was discovered that may be responsible for variation in Ser content. Finally, a ASMT gene was found via sequence alignment against its rice homolog. Validations of these candidate genes were performed by in vitro enzymatic reactions using proteins purified following recombinant expression in Escherichia coli, transient gene expression in tobacco, and transgenic approaches in wheat. Our results thus provide the first comprehensive investigation into the Mel pathway metabolites, and a swift candidate gene identification via forward-genetics strategies, in common wheat.

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.

Dissecting Complex Traits Using Omics Data: A Review on the Linear Mixed Models and Their Application in GWAS

Genome-wide association study (GWAS) is the most popular approach to dissecting complex traits in plants, humans, and animals. Numerous methods and tools have been proposed to discover the causal variants for GWAS data analysis. Among them, linear mixed models (LMMs) are widely used statistical methods for regulating confounding factors, including population structure, resulting in increased computational proficiency and statistical power in GWAS studies. Recently more attention has been paid to pleiotropy, multi-trait, gene-gene interaction, gene-environment interaction, and multi-locus methods with the growing availability of large-scale GWAS data and relevant phenotype samples. In this review, we have demonstrated all possible LMMs-based methods available in the literature for GWAS. We briefly discuss the different LMM methods, software packages, and available open-source applications in GWAS. Then, we include the advantages and weaknesses of the LMMs in GWAS. Finally, we discuss the future perspective and conclusion. The present review paper would be helpful to the researchers for selecting appropriate LMM models and methods quickly for GWAS data analysis and would benefit the scientific society.

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.

Annual Herbaceous Plants Exhibit Altered Morphological Traits in Response to Altered Precipitation and Drought Patterns in Semiarid Sandy Grassland, Northern China

The frequency and intensity of extreme precipitation events and severe drought are predicted to increase in semiarid areas due to global climate change. Plant morphological traits can reflect plant responses to a changing environment, such as altered precipitation or drought patterns. In this study, we examined the response of morphological traits of root, stem, leaf and reproduction meristems of annual herbaceous species to altered precipitation and drought patterns in a semiarid sandy grassland. The study involved a control treatment (100% of background precipitation) and the following six altered precipitation treatments: (1) P(+): precipitation increased by 30%, (2) P(++): precipitation increased by 60%, (3) P(-): precipitation decreased by 30%, (4) P(--): precipitation decreased by 60%, (5) drought 1 (D1): 46-day drought from May 1st to June 15th, and (6) drought 2 (D2): 46-day drought from July 1st to August 15th. P(++) significantly increased root length, flower length-to-width ratio, both P(+) and P(++) significantly increased stem length and flower number in the plant growing seasons, while all of them decreased under P(-) and P(--). The annual herbaceous plants marginally increased the number of second-level stem branches and stem diameter in order to better resist the severe drought stress under P(--). P(+) and P(++) increased the root, stem, leaf, and flower dry weight, with the flower dry weight accounting for a larger proportion than the other aboveground parts. Under D2, the plants used the limited water resources more efficiently by increasing the root-to-shoot ratio compared with P(-), P(--) and D1, which reflects biomass allocation to belowground increased. The linear mixed-effects models and redundancy analysis showed that the root-to-shoot ratio and the dry weight of various plant components were significantly affected by morphological traits and altered precipitation magnitude. Our results showed that the herbaceous species have evolved morphological trait responses that allow them to adapt to climate change. Such differences in morphological traits may ultimately affect the growing patterns of annual herbaceous species, enhancing their drought-tolerant capacity in semiarid sandy grassland during the ongoing climate change.

Spikelet abortion in six-rowed barley is mainly influenced by final spikelet number, with potential spikelet number acting as a suppressor trait

The potential to increase barley grain yield lies in the indeterminate nature of its inflorescence meristem, which produces spikelets, the basic reproductive unit in grasses that are linked to reproductive success. During early reproductive growth, barley spikes pass through the maximum yield potential-a stage after which no new spikelet ridges are produced. Subsequently, spikelet abortion (SA), a phenomenon in which spikelets abort during spike growth, imposes a bottleneck for increasing the grain yield potential. Here, we studied the potential of main culm spikes by counting potential spikelet number (PSN) and final spikelet number (FSN), and computed the corresponding SA (%) in a panel of 417 six-rowed spring barleys. Our phenotypic data analyses showed a significantly large within- and across-years genotypic variation with high broad-sense heritability estimates for all the investigated traits, including SA. Asian accessions displayed the lowest SA, indicating the presence of favourable alleles that may be exploited in breeding programs. A significantly negative Pearson's product-moment correlation was observed between FSN and SA. Our path analysis revealed that PSN and FSN explain 93% of the observed phenotypic variability for SA, with PSN behaving as a suppressor trait that magnifies the effect of FSN. Based on a large set of diverse barley accessions, our results provide a deeper phenotypic understanding of the quantitative genetic nature of SA, its association with traits of high agronomic importance, and a resource for further genetic analyses.

High-LD SNP markers exhibiting pleiotropic effects on salt tolerance at germination and seedlings stages in spring wheat

Salt tolerance at germination and seedling growth stages was investigated. GWAS revealed nine genomic regions with pleiotropic effects on salt tolerance. Salt tolerant genotypes were identified for future breeding program. With 20% of the irrigated land worldwide affected by it, salinity is a serious threat to plant development and crop production. While wheat is the most stable food source worldwide, it has been classified as moderately tolerant to salinity. In several crop plants; such as barley, maize and rice, it has been shown that salinity tolerance at seed germination and seedling establishment is under polygenic control. As yield was the ultimate goal of breeders and geneticists, less attention has been paid to understanding the genetic architecture of salt tolerance at early stages. Thus, the genetic control of salt tolerance at these stages is poorly understood relative to the late stages. In the current study, 176 genotypes of spring wheat were tested for salinity tolerance at seed germination and seedling establishment. Genome-Wide Association Study (GWAS) has been used to identify the genomic regions/genes conferring salt tolerance at seed germination and seedling establishment. Salinity stress negatively impacted all germination and seedling development parameters. A set of 137 SNPs showed significant association with the traits of interest. Across the whole genome, 33 regions showed high linkage disequilibrium (LD). These high LD regions harbored 15 SNPs with pleiotropic effect (i.e. SNPs that control more than one trait). Nine genes belonging to different functional groups were found to be associated with the pleiotropic SNPs. Noteworthy, chromosome 2B harbored the gene TraesCS2B02G135900 that acts as a potassium transporter. Remarkably, one SNP marker, reported in an early study, associated with salt tolerance was validated in this study. Our findings represent potential targets of genetic manipulation to understand and improve salinity tolerance in wheat.

Identification and characterization of QTL for spike morphological traits, plant height and heading date derived from the D genome of natural and resynthetic allohexaploid wheat

QHd.cau-7D.1 for heading date was delimited into the physical interval of approximately 17.38 Mb harboring three CONSTANS-like zinc finger genes. Spike morphological traits, plant height and heading date play important roles in yield improvement of wheat. To reveal the genetic factors that controlling spike morphological traits, plant height and heading date on the D genome, we conducted analysis of quantitative traits locus (QTL) using 198 F_7:8 recombinant inbred lines (RILs) derived from a cross between the common wheat TAA10 and resynthesized allohexaploid wheat XX329 with similar AABB genomes. A total of 23 environmentally stable QTL on the D sub-genome for spike length (SL), fertile spikelet number per spike (FSN), sterile spikelet number per spike (SSN), total spikelet number per spike (TSN), spike compactness (SC), plant height (PHT) and heading date (HD) were detected, among which eight appeared to be novel QTL. Furthermore, QHd.cau-7D.1 and QPht.cau-7D.2 shared identical confidence interval and were delimited into the physical interval of approximately 17.38 Mb with 145 annotated genes, including three CONSTANS-like zinc finger genes (TraesCS7D02G209000, TraesCS7D02G213000 and TraesCS7D02G220300). This study will help elucidate the molecular mechanism of the seven traits (SL, FSN, SSN, TSN, SC, PHT and HD) and provide a potentially valuable resource for genetic improvement.

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.

Raise and characterization of a bread wheat hybrid line (Tulaykovskaya 10 × Saratovskaya 29) with chromosome 6Agi2 introgressed from Thinopyrum intermedium

Wheatgrass Thinopyrum intermedium is a source of agronomically valuable traits for common wheat. Partial wheat–wheatgrass amphidiploids and lines with wheatgrass chromosome substitutions are extensively used as intermediates in breeding programs. Line Agis 1 (6Agi2/6D) is present in the cultivar Tulaykovskaya 10 pedigree. Wheatgrass chromosome 6Agi2 carries multiple resistance to fungal diseases in various ecogeographical zones. In this work, we studied the transfer of chromosome 6Agi2 in hybrid populations Saratovskaya 29 × skaya 10 (S29 × T10) and Tulaykovskaya 10 × Saratovskaya 29 (T10 × S29). Chromosome 6Agi2 was identif ied by PCR with chromosome-specif ic primers and by genomic in situ hybridization (GISH). According to molecular data, 6Agi2 was transmitted to nearly half of the plants tested in the F2 and F3 generations. A new breeding line 49-14 (2n = 42) with chromosome pair 6Agi2 was isolated and characterized in T10 × S29 F5 by GISH. According to the results of our f ield experiment in 2020, the line had high productivity traits. The grain weights per plant (10.04 ± 0.93 g) and the number of grains per plant (259.36 ± 22.49) did not differ signif icantly from the parent varieties. The number of grains per spikelet in the main spike was signif icantly higher than in S29 ( p ≤ 0.001) or T10 ( p ≤ 0.05). Plants were characterized by the ability to set 3.77 ± 0.1 grains per spikelet, and this trait varied among individuals from 2.93 to 4.62. The grain protein content was 17.91 %, and the gluten content, 40.55 %. According to the screening for fungal disease resistance carried out in the f ield in 2018 and 2020, chromosome 6Agi2 makes plants retain immunity to the West Siberian population of brown rust and to dominant races of stem rust. It also provides medium resistant and medium susceptible types of response to yellow rust. The possibility of using lines/varieties of bread wheat with wheatgrass chromosomes 6Agi2 in breeding in order to increase protein content in the grain, to confer resistance to leaf diseases on plants and to create multif lowered forms is discussed.

Strategies of grain number determination differentiate barley row types

Gaining knowledge on fundamental interactions of various yield components is crucial to improve yield potential in small grain cereals. It is well known in barley that increasing grain number greatly improves yield potential; however, the yield components determining grain number and their association in barley row types are less explored. In this study, we assessed different yield components such as potential spikelet number (PSN), spikelet survival (SSL), spikelet number (SN), grain set (GS), and grain survival (GSL), as well as their interactions with grain number by using a selected panel of two- and six-rowed barley types. Also, to analyze the stability of these interactions, we performed the study in the greenhouse and the field. From this study, we found that in two-rowed barley, grain number determination is strongly influenced by PSN rather than SSL and/or GS in both growth conditions. Conversely, in six-rowed barley, grain number is associated with SSL instead of PSN and/or GS. Thus, our study showed that increasing grain number might be possible by augmenting PSN in two-rowed genotypes, while for six-rowed genotypes SSL needs to be improved. We speculate that this disparity of grain number determination in barley row types might be due to the fertility of lateral spikelets. Collectively, this study revealed that grain number in two-rowed barley largely depends on the developmental trait, PSN, while in six-rowed barley, it mainly follows the ability for SSL.

Genome wide association study of plant height and tiller number in hulless barley

Hulless barley (Hordeum vulgare L. var. nudum), also called naked barley, is a unique variety of cultivated barley. The genome-wide specific length amplified fragment sequencing (SLAF-seq) method is a rapid deep sequencing technology that is used for the selection and identification of genetic loci or markers. In this study, we collected 300 hulless barley accessions and used the SLAF-seq method to identify candidate genes involved in plant height (PH) and tiller number (TN). We obtained a total of 1407 M paired-end reads, and 228,227 SLAF tags were developed. After filtering using an integrity threshold of >0.8 and a minor allele frequency of >0.05, 14,504,892 single-nucleotide polymorphisms (SNP) loci were screened out. The remaining SNPs were used for the construction of a neighbour-joining phylogenetic tree, and the three subcluster members showed no obvious differentiation among regional varieties. We used a genome wide association study approach to identify 1006 and 113 SNPs associated with TN and PH, respectively. Based on best linear unbiased predictors (BLUP), 41 and 29 SNPs associated with TN and PH, respectively. Thus, several of genes, including Hd3a and CKX5, may be useful candidates for the future genetic breeding of hulless barley. Taken together, our results provide insight into the molecular mechanisms controlling barley architecture, which is important for breeding and yield.

Genome-wide association analysis uncovers the genetic architecture of tradeoff between flowering date and yield components in sesame

BACKGROUND

Unrevealing the genetic makeup of crop morpho-agronomic traits is essential for improving yield quality and sustainability. Sesame (Sesamum indicum L.) is one of the oldest oil-crops in the world. Despite its economic and agricultural importance, it is an 'orphan crop-plant' that has undergone limited modern selection, and, as a consequence preserved wide genetic diversity. Here we established a new sesame panel (SCHUJI) that contains 184 genotypes representing wide phenotypic variation and is geographically distributed. We harnessed the natural variation of this panel to perform genome-wide association studies for morpho-agronomic traits under the Mediterranean climate conditions.

RESULTS

Field-based phenotyping of the SCHUJI panel across two seasons exposed wide phenotypic variation for all traits. Using 20,294 single-nucleotide polymorphism markers, we detected 50 genomic signals associated with these traits. Major genomic region on LG2 was associated with flowering date and yield-related traits, exemplified the key role of the flowering date on productivity.

CONCLUSIONS

Our results shed light on the genetic architecture of flowering date and its interaction with yield components in sesame and may serve as a basis for future sesame breeding programs in the Mediterranean basin.

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.

Large-scale integration of meta-QTL and genome-wide association study discovers the genomic regions and candidate genes for yield and yield-related traits in bread wheat

Based on the large-scale integration of meta-QTL and Genome-Wide  Association Study, 76 high-confidence MQTL regions and 237 candidate genes that affected wheat yield and yield-related traits were discovered. Improving yield and yield-related traits are key goals in wheat breeding program. The integration of accumulated wheat genetic resources provides an opportunity to uncover important genomic regions and candidate genes that affect wheat yield. Here, a comprehensive meta-QTL analysis was conducted on 2230 QTL of yield-related traits obtained from 119 QTL studies. These QTL were refined into 145 meta-QTL (MQTL), and 89 MQTL were verified by GWAS with different natural populations. The average confidence interval (CI) of these MQTL was 2.92 times less than that of the initial QTL. Furthermore, 76 core MQTL regions with a physical distance less than 25 Mb were detected. Based on the homology analysis and expression patterns, 237 candidate genes in the MQTL involved in photoperiod response, grain development, multiple plant growth regulator pathways, carbon and nitrogen metabolism and spike and flower organ development were determined. A novel candidate gene TaKAO-4A was confirmed to be significantly associated with grain size, and a CAPS marker was developed based on its dominant haplotype. In summary, this study clarified a method based on the integration of meta-QTL, GWAS and homology comparison to reveal the genomic regions and candidate genes that affect important yield-related traits in wheat. This work will help to lay a foundation for the identification, transfer and aggregation of these important QTL or candidate genes in wheat high-yield breeding.

Long-range assembly of sequences helps to unravel the genome structure and small variation of the wheat-Haynaldia villosa translocated chromosome 6VS.6AL

Genomics studies in wild species of wheat have been limited due to the lack of references; however, new technologies and bioinformatics tools have much potential to promote genomic research. The wheat-Haynaldia villosa translocation line T6VS·6AL has been widely used as a backbone parent of wheat breeding in China. Therefore, revealing the genome structure of translocation chromosome 6VS·6AL will clarify how this chromosome formed and will help to determine how it affects agronomic traits. In this study, chromosome flow sorting, NGS sequencing and Chicago long-range linkage assembly were innovatively used to produce the assembled sequences of 6VS·6AL, and gene prediction and genome structure characterization at the molecular level were effectively performed. The analysis discovered that the short arm of 6VS·6AL was actually composed of a large distal segment of 6VS, a small proximal segment of 6AS and the centromere of 6A, while the collinear region in 6VS corresponding to 230-260 Mb of 6AS-Ta was deleted when the recombination between 6VS and 6AS occurred. In addition to the molecular mechanism of the increased grain weight and enhanced spike length produced by the translocation chromosome, it may be correlated with missing GW2-V and an evolved NRT-V cluster. Moreover, a fine physical bin map of 6VS was constructed by the high-throughput developed 6VS-specific InDel markers and a series of newly identified small fragment translocation lines involving 6VS. This study will provide essential information for mining of new alien genes carried by the 6VS·6AL translocation chromosome.

Genome-wide association study of yield components in spring wheat collection harvested under two water regimes in Northern Kazakhstan

Background

Bread wheat is the most important cereal in Kazakhstan, where it is grown on over 12 million hectares. One of the major constraints affecting wheat grain yield is drought due to the limited water supply. Hence, the development of drought-resistant cultivars is critical for ensuring food security in this country. Therefore, identifying quantitative trait loci (QTLs) associated with drought tolerance as an essential step in modern breeding activities, which rely on a marker-assisted selection approach.

Methods

A collection of 179 spring wheat accessions was tested under irrigated and rainfed conditions in Northern Kazakhstan over three years (2018, 2019, and 2020), during which data was collected on nine traits: heading date (HD), seed maturity date (SMD), plant height (PH), peduncle length (PL), number of productive spikes (NPS), spike length (SL), number of kernels per spike (NKS), thousand kernel weight (TKW), and kernels yield per m² (YM2). The collection was genotyped using a 20,000 (20K) Illumina iSelect SNP array, and 8,662 polymorphic SNP markers were selected for a genome-wide association study (GWAS) to identify QTLs for targeted agronomic traits.

Results

Out of the total of 237 discovered QTLs, 50 were identified as being stable QTLs for irrigated and rainfed conditions in the Akmola region, Northern Kazakhstan; the identified QTLs were associated with all the studied traits except PH. The results indicate that nine QTLs for HD and 11 QTLs for SMD are presumably novel genetic factors identified in the irrigated and rainfed conditions of Northern Kazakhstan. The identified SNP markers of the QTLs for targeted traits in rainfed conditions can be applied to develop new competitive spring wheat cultivars in arid zones using a marker-assisted selection approach.

Mapping QTL for spike fertility and related traits in two doubled haploid wheat (Triticum aestivum L.) populations

BACKGROUND

In breeding programs, the selection of cultivars with the highest yield potential consisted in the selection of the yield per se, which resulted in cultivars with higher grains per spike (GN) and occasionally increased grain weight (GW) (main numerical components of the yield). In this study, quantitative trait loci (QTL) for GW, GN and spike fertility traits related to GN determination were mapped using two doubled haploid (DH) populations (Baguette Premium 11 × BioINTA 2002 and Baguette 19 × BioINTA 2002).

RESULTS

In total 305 QTL were identified for 14 traits, out of which 12 QTL were identified in more than three environments and explained more than 10% of the phenotypic variation in at least one environment. Eight hotspot regions were detected on chromosomes 1A, 2B, 3A, 5A, 5B, 7A and 7B in which at least two major and stable QTL sheared confidence intervals. QTL on two of these regions (R5A.1 and R5A.2) have previously been described, but the other six regions are novel.

CONCLUSIONS

Based on the pleiotropic analysis within a robust physiological model we conclude that two hotspot genomic regions (R5A.1 and R5A.2) together with the QGW.perg-6B are of high relevance to be used in marker assisted selection in order to improve the spike yield potential. All the QTL identified for the spike related traits are the first step to search for their candidate genes, which will allow their better manipulation in the future.

Exploring the genic resources underlying metabolites through mGWAS and mQTL in wheat: From large-scale gene identification and pathway elucidation to crop improvement

Common wheat (Triticum aestivum L.) is a leading cereal crop, but has lagged behind with respect to the interpretation of the molecular mechanisms of phenotypes compared with other major cereal crops such as rice and maize. The recently available genome sequence of wheat affords the pre-requisite information for efficiently exploiting the potential molecular resources for decoding the genetic architecture of complex traits and identifying valuable breeding targets. Meanwhile, the successful application of metabolomics as an emergent large-scale profiling methodology in several species has demonstrated this approach to be accessible for reaching the above goals. One such productive avenue is combining metabolomics approaches with genetic designs. However, this trial is not as widespread as that for sequencing technologies, especially when the acquisition, understanding, and application of metabolic approaches in wheat populations remain more difficult and even arguably underutilized. In this review, we briefly introduce the techniques used in the acquisition of metabolomics data and their utility in large-scale identification of functional candidate genes. Considerable progress has been made in delivering improved varieties, suggesting that the inclusion of information concerning these metabolites and genes and metabolic pathways enables a more explicit understanding of phenotypic traits and, as such, this procedure could serve as an -omics-informed roadmap for executing similar improvement strategies in wheat and other species.

The physiology and genetics behind fruiting efficiency: a promising spike trait to improve wheat yield potential

Fruiting efficiency (FE, grains per g of spike dry weight at anthesis) was proposed as a promising spike trait to improve wheat yield potential, based on its functional relationship with grain number determination and the evidence of trait variability in elite germplasm. During the last few years, we have witnessed great advances in the understanding of the physiological and genetic basis of this trait. The present review summarizes the recent heritability estimations and the genetic gains obtained when fruiting efficiency was measured at maturity (FEm, grains per g of chaff) and used as selection criterion. In addition, we propose spike ideotypes for contrasting fruiting efficiencies based on the fertile floret efficiency (FFE, fertile florets per g of spike dry weight at anthesis) and grain set (grains per fertile floret), together with other spike fertility-related traits. We also review novel genes and quantitative trait loci available for using marker-assisted selection for fruiting efficiency and other spike fertility traits. The possible trade-off between FE and grain weight and the genes reported to alter this relation are also considered. Finally, we discuss the benefits and future steps towards the use of fruiting efficiency as a selection criterion in breeding programs.

Identification and Validation of a Novel Locus Controlling Spikelet Number in Bread Wheat (Triticum aestivum L.)

Spikelet number is an important target trait for wheat yield improvement. Thus, the identification and verification of novel quantitative trait locus (QTL)/genes controlling spikelet number are essential for dissecting the underlying molecular mechanisms and hence for improving grain yield. In the present study, we constructed a high-density genetic map for the Kechengmai1/Chuanmai42 doubled haploid (DH) population using 13,068 single-nucleotide polymorphism (SNP) markers from the Wheat 55K SNP array. A comparison between the genetic and physical maps indicated high consistence of the marker orders. Based on this genetic map, a total of 27 QTLs associated with total spikelet number per spike (TSN) and fertile spikelet number per spike (FSN) were detected on chromosomes 1B, 1D, 2B, 2D, 3D, 4A, 4D, 5A, 5B, 5D, 6A, 6B, and 7D in five environments. Among them, five QTLs on chromosome 2D, 3D, 5A, and 7D were detected in multiple environments and combined QTL analysis, explaining the phenotypic variance ranging from 3.64% to 23.28%. Particularly, QTsn/Fsn.cib-3D for TSN and FSN [phenotypic variation explained (PVE) = 5.97-23.28%, limit of detection (LOD) = 3.73-18.51] is probably a novel locus and located in a 4.5-cM interval on chromosome arm 3DL flanking by the markers AX-110914105 and AX-109429351. This QTL was further validated in other two populations with different genetic backgrounds using the closely linked Kompetitive Allele-Specific PCR (KASP) marker KASP_AX-110914105. The results indicated that QTsn/Fsn.cib-3D significantly increased the TSN (5.56-7.96%) and FSN (5.13-9.35%), which were significantly correlated with grain number per spike (GNS). We also preliminary analyzed the candidate genes within this locus by sequence similarity, spatial expression patterns, and collinearity analysis. These results provide solid foundation for future fine mapping and cloning of QTsn/Fsn.cib-3D. The developed and validated KASP markers could be utilized in molecular breeding aiming to increase the grain yield in wheat.

Genome-Wide Association Study of Morpho-Physiological Traits in Aegilops tauschii to Broaden Wheat Genetic Diversity

Aegilops tauschii, the D-genome donor of bread wheat, is a storehouse of genetic diversity that can be used for wheat improvement. This species consists of two main lineages (TauL1 and TauL2) and one minor lineage (TauL3). Its morpho-physiological diversity is large, with adaptations to a wide ecological range. Identification of allelic diversity in Ae. tauschii is of utmost importance for efficient breeding and widening of the genetic base of wheat. This study aimed at identifying markers or genes associated with morpho-physiological traits in Ae. tauschii, and at understanding the difference in genetic diversity between the two main lineages. We performed genome-wide association studies of 11 morpho-physiological traits of 343 Ae. tauschii accessions representing the entire range of habitats using 34,829 DArTseq markers. We observed a wide range of morpho-physiological variation among all accessions. We identified 23 marker-trait associations (MTAs) in all accessions, 15 specific to TauL1 and eight specific to TauL2, suggesting independent evolution in each lineage. Some of the MTAs could be novel and have not been reported in bread wheat. The markers or genes identified in this study will help reveal the genes controlling the morpho-physiological traits in Ae. tauschii, and thus in bread wheat even if the plant morphology is different.

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.

Cytological and molecular characterizations of a novel 2A nullisomic line derived from a widely-grown wheat cultivar Zhoumai 18 conferring male sterility

A dwarf, multi-pistil and male sterile dms mutant was previously reported by us. However, the genetic changes in this dms are unclear. To examine the genetic changes, single nucleotide polymorphism (SNP) association, chromosome counting, and high-resolution chromosome fluorescence in situ hybridization (FISH) techniques were employed. By comparing tall plants (T) with dwarf plants (D) in the offspring of dms mutant plants, SNP association analysis indicated that most SNPs were on chromosome 2A. There were three types in offspring of dms plants, with 42, 41 and 40 chromosomes respectively. High-resolution chromosome painting analysis demonstrated that T plants had all 42 wheat chromosomes; the medium plants (M) had 41 chromosomes, lacking one chromosome 2A; while D plants had 40 wheat chromosomes, and lacked both 2A chromosomes. These data demonstrated that dms resulted from a loss of chromosome 2A. We identified 23 genes on chromosome 2A which might be involved in the development of stamens or pollen grains. These results lay a solid foundation for further analysis of the molecular mechanisms of wheat male sterility. Because D plants can be used as a female parent to cross with other wheat genotypes, dms is a unique germplasm for any functional study of chromosome 2A and wheat breeding specifically targeting genes on 2A.

Unlocking the relationships among population structure, plant architecture, growing season, and environmental adaptation in Henan wheat cultivars

BACKGROUND

Ecological environments shape plant architecture and alter the growing season, which provides the basis for wheat genetic improvement. Therefore, understanding the genetic basis of grain yield and yield-related traits in specific ecological environments is important.

RESULTS

A structured panel of 96 elite wheat cultivars grown in the High-yield zone of Henan province in China was genotyped using an Illumina iSelect 90 K SNP assay. Selection pressure derived from ecological environments of mountain front and plain region provided the initial impetus for population divergence. This determined the dominant traits in two subpopulations (spike number and spike percentage were dominance in subpopulation 2:1; thousand-kernel weight, grain filling rate (GFR), maturity date (MD), and fertility period (FP) were dominance in subpopulation 2:2), which was also consistent with their inheritance from the donor parents. Genome wide association studies identified 107 significant SNPs for 12 yield-related traits and 10 regions were pleiotropic to multiple traits. Especially, GY was co-located with MD/FP, GFR and HD at QTL-ple5A, QTL-ple7A.1 and QTL-ple7B.1 region. Further selective sweep analysis revealled that regions under selection were around QTLs for these traits. Especially, grain yield (GY) is positively correlated with MD/FP and they were co-located at the VRN-1A locus. Besides, a selective sweep signal was detected at VRN-1B locus which was only significance to MD/FP.

CONCLUSIONS

The results indicated that extensive differential in allele frequency driven by ecological selection has shaped plant architecture and growing season during yield improvement. The QTLs for yield and yield components detected in this study probably be selectively applied in molecular breeding.

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

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

Appraising the Genetic Architecture of Kernel Traits in Hexaploid Wheat Using GWAS

Kernel morphology is one of the major yield traits of wheat, the genetic architecture of which is always important in crop breeding. In this study, we performed a genome-wide association study (GWAS) to appraise the genetic architecture of the kernel traits of 319 wheat accessions using 22,905 single nucleotide polymorphism (SNP) markers from a wheat 90K SNP array. As a result, 111 and 104 significant SNPs for Kernel traits were detected using four multi-locus GWAS models (mrMLM, FASTmrMLM, FASTmrEMMA, and pLARmEB) and three single-locus models (FarmCPU, MLM, and MLMM), respectively. Among the 111 SNPs detected by the multi-locus models, 24 SNPs were simultaneously detected across multiple models, including seven for kernel length, six for kernel width, six for kernels per spike, and five for thousand kernel weight. Interestingly, the five most stable SNPs (RAC875_29540_391, Kukri_07961_503, tplb0034e07_1581, BS00074341_51, and BobWhite_049_3064) were simultaneously detected by at least three multi-locus models. Integrating these newly developed multi-locus GWAS models to unravel the genetic architecture of kernel traits, the mrMLM approach detected the maximum number of SNPs. Furthermore, a total of 41 putative candidate genes were predicted to likely be involved in the genetic architecture underlining kernel traits. These findings can facilitate a better understanding of the complex genetic mechanisms of kernel traits and may lead to the genetic improvement of grain yield in wheat.

Metabolite-based genome-wide association study enables dissection of the flavonoid decoration pathway of wheat kernels

The marriage of metabolomic approaches with genetic design has proven a powerful tool in dissecting diversity in the metabolome and has additionally enhanced our understanding of complex traits. That said, such studies have rarely been carried out in wheat. In this study, we detected 805 metabolites from wheat kernels and profiled their relative contents among 182 wheat accessions, conducting a metabolite-based genome-wide association study (mGWAS) utilizing 14 646 previously described polymorphic SNP markers. A total of 1098 mGWAS associations were detected with large effects, within which 26 candidate genes were tentatively designated for 42 loci. Enzymatic assay of two candidates indicated they could catalyse glucosylation and subsequent malonylation of various flavonoids and thereby the major flavonoid decoration pathway of wheat kernel was dissected. Moreover, numerous high-confidence genes associated with metabolite contents have been provided, as well as more subdivided metabolite networks which are yet to be explored within our data. These combined efforts presented the first step towards realizing metabolomics-associated breeding of wheat.

Multi-Locus GWAS of Quality Traits in Bread Wheat: Mining More Candidate Genes and Possible Regulatory Network

In wheat breeding, improved quality traits, including grain quality and dough rheological properties, have long been a critical goal. To understand the genetic basis of key quality traits of wheat, two single-locus and five multi-locus GWAS models were performed for six grain quality traits and three dough rheological properties based on 19, 254 SNPs in 267 bread wheat accessions. As a result, 299 quantitative trait nucleotides (QTNs) within 105 regions were identified to be associated with these quality traits in four environments. Of which, 40 core QTN regions were stably detected in at least three environments, 19 of which were novel. Compared with the previous studies, these novel QTN regions explained smaller phenotypic variation, which verified the advantages of the multi-locus GWAS models in detecting important small effect QTNs associated with complex traits. After characterization of the function and expression in-depth, 67 core candidate genes involved in protein/sugar synthesis, histone modification and the regulation of transcription factor were observed to be associated with the formation of grain quality, which showed that multi-level regulations influenced wheat grain quality. Finally, a preliminary network of gene regulation that may affect wheat quality formation was inferred. This study verified the power and reliability of multi-locus GWAS methods in wheat quality trait research, and increased the understanding of wheat quality formation mechanisms. The detected QTN regions and candidate genes in this study could be further used for gene cloning and marker-assisted selection in high-quality breeding of bread wheat.

TaSPL13 regulates inflorescence architecture and development in transgenic wheat (Triticum aestivum L.)

The SQUAMOSA promoter-binding protein-like (SPL) proteins play vital roles in plant growth and development in rice (Oryza sative L.) and Arabidopsis thaliana (L.) Heynh. However, few studies regarding the SPL proteins have been reported in wheat. In this study, 56 TaSPLs were clustered into eight groups according to an OsSPL phylogenetic comparison analysis. The expression patterns of TaSPLs in different tissues were analysed by RNA-seq data, and partial results were confirmed by qRT-PCR. Based on the above results, genes such as TaSPL13 and TaSPL15 may be involved in spike or seed development in wheat. Multiple genes that regulate the inflorescence architecture of rice have been identified. Additionally, studies on the genes associated with spikelet development in wheat have been reported relatively rarely. Here, TaSPL13-2B was transferred into wheat cv. Bobwhite. Compared with the wild type, the transgenic lines showed significant increases in the number of florets and grains per spike, indicating that TaSPL13-2B could influence the floret development of wheat. TaSPL13-2B was transferred into rice cv. Nipponbare, which demonstrated that TaSPL13-2B can modify panicle architecture in transgenic rice, with significant increases in panicle length, the number and length of primary branches, and the number of secondary branches.

The Wheat 660K SNP array demonstrates great potential for marker-assisted selection in polyploid wheat

The rapid development and application of molecular marker assays have facilitated genomic selection and genome-wide linkage and association studies in wheat breeding. Although PCR-based markers (e.g. simple sequence repeats and functional markers) and genotyping by sequencing have contributed greatly to gene discovery and marker-assisted selection, the release of a more accurate and complete bread wheat reference genome has resulted in the design of single-nucleotide polymorphism (SNP) arrays based on different densities or application targets. Here, we evaluated seven types of wheat SNP arrays in terms of their SNP number, distribution, density, associated genes, heterozygosity and application. The results suggested that the Wheat 660K SNP array contained the highest percentage (99.05%) of genome-specific SNPs with reliable physical positions. SNP density analysis indicated that the SNPs were almost evenly distributed across the whole genome. In addition, 229 266 SNPs in the Wheat 660K SNP array were located in 66 834 annotated gene or promoter intervals. The annotated genes revealed by the Wheat 660K SNP array almost covered all genes revealed by the Wheat 35K (97.44%), 55K (99.73%), 90K (86.9%) and 820K (85.3%) SNP arrays. Therefore, the Wheat 660K SNP array could act as a substitute for other 6 arrays and shows promise for a wide range of possible applications. In summary, the Wheat 660K SNP array is reliable and cost-effective and may be the best choice for targeted genotyping and marker-assisted selection in wheat genetic improvement.

Genetics of yield, abiotic stress tolerance and biofortification in wheat (Triticum aestivum L.)

A review of the available literature on genetics of yield and its component traits, tolerance to abiotic stresses and biofortification should prove useful for future research in wheat in the genomics era. The work reviewed in this article mainly covers the available information on genetics of some important quantitative traits including yield and its components, tolerance to abiotic stresses (heat, drought, salinity and pre-harvest sprouting = PHS) and biofortification (Fe/Zn and phytate contents with HarvestPlus Program) in wheat. Major emphasis is laid on the recent literature on QTL interval mapping and genome-wide association studies, giving lists of known QTL and marker-trait associations. Candidate genes for different traits and the cloned and characterized genes for yield traits along with the molecular mechanism are also described. For each trait, an account of the present status of marker-assisted selection has also been included. The details of available results have largely been presented in the form of tables; some of these tables are included as supplementary files.

Genome-wide and SNP network analyses reveal genetic control of spikelet sterility and yield-related traits in wheat

Revealing the genetic factors underlying yield and agronomic traits in wheat are an imperative need for covering the global food demand. Yield boosting requires a deep understanding of the genetic basis of grain yield-related traits (e.g., spikelet fertility and sterility). Here, we have detected much natural variation among ancient hexaploid wheat accessions in twenty-two agronomic traits collected over eight years of field experiments. A genome-wide association study (GWAS) using 15 K single nucleotide polymorphisms (SNPs) was applied to detect the genetic basis of studied traits. Subsequently, the GWAS output was reinforced via other statistical and bioinformatics analyses to detect putative candidate genes. Applying the genome-wide SNP-phenotype network defined the most decisive SNPs underlying the traits. Six pivotal SNPs, co-located physically within the genes encoding enzymes, hormone response, metal ion transport, and response to oxidative stress have been identified. Of these, metal ion transport and Gibberellin 2-oxidases (GA2oxs) genes showed strong involvement in controlling the spikelet sterility, which had not been reported previously in wheat. SNP-gene haplotype analysis confirmed that these SNPs influence spikelet sterility, especially the SNP co-located on the exon of the GA2ox gene. Interestingly, these genes were highly expressed in the grain and spike, demonstrating their pivotal role in controlling the trait. The integrative analysis strategy applied in this study, including GWAS, SNP-phenotype network, SNP-gene haplotype, expression analysis, and genome-wide prediction (GP), empower the identification of functional SNPs and causal genes. GP outputs obtained in this study are encouraging for the implementation of the traits to accelerate yield improvement by making an early prediction of complex yield-related traits in wheat. Our findings demonstrate the usefulness of the ancient wheat material as a valuable resource for yield-boosting. This is the first comprehensive genome-wide analysis for spikelet sterility in wheat, and the results provide insights into yield improvement.

APETALA 2-like genes AP2L2 and Q specify lemma identity and axillary floral meristem development in wheat

The spikelet is the basic unit of the grass inflorescence. In tetraploid (Triticum turgidum) and hexaploid wheat (Triticum aestivum), the spikelet is a short indeterminate branch with two proximal sterile bracts (glumes) followed by a variable number of florets, each including a bract (lemma) with an axillary flower. Varying levels of miR172 and/or its target gene Q (AP2L5) result in gradual transitions of glumes to lemmas, and vice versa. Here, we show that AP2L5 and its related paralog AP2L2 play critical and redundant roles in the specification of axillary floral meristems and lemma identity. AP2L2, also targeted by miR172, displayed similar expression profiles to AP2L5 during spikelet development. Loss-of-function mutants in both homeologs of AP2L2 (henceforth ap2l2) developed normal spikelets, but ap2l2 ap2l5 double mutants generated spikelets with multiple empty bracts before transitioning to florets. The coordinated nature of these changes suggest an early role of these genes in floret development. Moreover, the flowers of ap2l2 ap2l5 mutants showed organ defects in paleas and lodicules, including the homeotic conversion of lodicules into carpels. Mutations in the miR172 target site of AP2L2 were associated with reduced plant height, more compact spikes, promotion of lemma-like characters in glumes and smaller lodicules. Taken together, our results show that the balance in the expression of miR172 and AP2-like genes is crucial for the correct development of spikelets and florets, and that this balance has been altered during the process of wheat and barley (Hordeum vulgare) domestication. The manipulation of this regulatory module provides an opportunity to modify spikelet architecture and improve grain yield.

Haplotype block analysis of an Argentinean hexaploid wheat collection and GWAS for yield components and adaptation

BACKGROUND

Increasing wheat (Triticum aestivum L.) production is required to feed a growing human population. In order to accomplish this task a deeper understanding of the genetic structure of cultivated wheats and the detection of genomic regions significantly associated with the regulation of important agronomic traits are necessary steps. To better understand the genetic basis and relationships of adaptation and yield related traits, we used a collection of 102 Argentinean hexaploid wheat cultivars genotyped with the 35k SNPs array, grown from two to six years in three different locations. Based on SNPs data and gene-related molecular markers, we performed a haplotype block characterization of the germplasm and a genome-wide association study (GWAS).

RESULTS

The genetic structure of the collection revealed four subpopulations, reflecting the origin of the germplasm used by the main breeding programs in Argentina. The haplotype block characterization showed 1268 blocks of different sizes spread along the genome, including highly conserved regions like the 1BS chromosome arm where the 1BL/1RS wheat/rye translocation is located. Based on GWAS we identified ninety-seven chromosome regions associated with heading date, plant height, thousand grain weight, grain number per spike and fruiting efficiency at harvest (FEh). In particular FEh stands out as a promising trait to raise yield potential in Argentinean wheats; we detected fifteen haplotypes/markers associated with increased FEh values, eleven of which showed significant effects in all three evaluated locations. In the case of adaptation, the Ppd-D1 gene is consolidated as the main determinant of the life cycle of Argentinean wheat cultivars.

CONCLUSION

This work reveals the genetic structure of the Argentinean hexaploid wheat germplasm using a wide set of molecular markers anchored to the Ref Seq v1.0. Additionally GWAS detects chromosomal regions (haplotypes) associated with important yield and adaptation components that will allow improvement of these traits through marker-assisted selection.

Genome-wide association study reveals genomic regions controlling root and shoot traits at late growth stages in wheat

BACKGROUND AND AIMS

Root system morphology is important for sustainable agriculture, but the genetic basis of root traits and their relationship to shoot traits remain to be elucidated. The aim of the present study was to dissect the genetic basis of root traits at late growth stages and its implications on shoot traits in wheat.

METHODS

Among 323 wheat accessions, we investigated phenotypic differences in root traits at booting and mid-grain fill stages in PVC tubes, shoot traits including plant height (PH), canopy temperature (CT) and grain yield per plant (YPP) in a field experiment, and performed a genome-wide association study with a Wheat 660K SNP Array.

KEY RESULTS

Deep-rooted accessions had lower CT and higher YPP than those with shallow roots, but no significant relationship was identified between root dry weight and shoot traits. Ninety-three significantly associated loci (SALs) were detected by the mixed linear model, among which three were hub SALs (Co-6A, Co-6B and Co-6D) associated with root depth at both booting and mid-grain fill stages, as well as CT and YPP. Minirhizotron system scanning results suggested that the causal genes in the three SALs may regulate root elongation in the field. The heritable independence between root depth and PH was demonstrated by linkage disequilibrium analysis. The YPP was significantly higher in genotypes which combined favourable marker alleles (FMAs) for root depth and PH, suggesting that a deep root and shorter plant height are suitable traits for pyramiding target alleles by molecular marker-assisted breeding.

CONCLUSIONS

These results uncovered promising genomic regions for functional gene discovery of root traits in the late growth period, enhanced understanding of correlation between root and shoot traits, and will facilitate intensive study on root morphology and breeding through molecular design.