Genome-Wide Association Analyses Identify QTL Hotspots for Yield and Component Traits in Durum Wheat Grown under Yield Potential, Drought, and Heat Stress Environments

Genomic analysis for heat and combined heat-drought resilience in bread wheat under field conditions

KEY MESSAGE

GWAS on a bread wheat panel with high D genome diversity identified novel alleles and QTLs associated with resilience to combined heat and drought stress under natural field conditions. As heat (H) and drought stresses occur concurrently under field conditions, studying them separately offers limited opportunities for wheat improvement. Here, a wheat diversity panel containing Aegilops tauschii introgressions was evaluated under H and combined heat-drought (HD) stresses to identify quantitative trait loci (QTLs) associated with resilience to the stresses, and to assess the practicability of harnessing Ae. tauschii diversity for breeding for combined stress resilience. Using genome-wide analysis, we identified alleles and QTLs on chromosomes 3D, 5D, and 7A controlling grain yield (GY), kernel number per spike, and thousand-kernel weight, and on 3D (521-549 Mbp) controlling GY alone. A strong marker-trait association (MTA) for GY stability on chromosome 3D (508.3 Mbp) explained 20.3% of the variation. Leaf traits-canopy temperature, vegetation index, and carbon isotope composition-were controlled by five QTLs on 2D (23-96, 511-554, and 606-614 Mbp), 3D (155-171 Mbp), and 5D (407-413 Mbp); some of them were pleiotropic for GY and yield-related traits. Further analysis revealed candidate genes, including GA20ox, regulating GY stability, and CaaX prenyl protease 2, regulating canopy temperature at the flowering stage, under H and HD stresses. As genome-wide association studies under HD in field conditions are scarce, our results provide genomic landmarks for wheat breeding to improve adaptation to H and HD conditions under climate change.

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.

Genome-Wide Association Study of Root Mealiness and Other Texture-Associated Traits in Cassava

Cassava breeders have made significant progress in developing new genotypes with improved agronomic characteristics such as improved root yield and resistance against biotic and abiotic stresses. However, these new and improved cassava (Manihot esculenta Crantz) varieties in cultivation in Nigeria have undergone little or no improvement in their culinary qualities; hence, there is a paucity of genetic information regarding the texture of boiled cassava, particularly with respect to its mealiness, the principal sensory quality attribute of boiled cassava roots. The current study aimed at identifying genomic regions and polymorphisms associated with natural variation for root mealiness and other texture-related attributes of boiled cassava roots, which includes fibre, adhesiveness (ADH), taste, aroma, colour, and firmness. We performed a genome-wide association (GWAS) analysis using phenotypic data from a panel of 142 accessions obtained from the National Root Crops Research Institute (NRCRI), Umudike, Nigeria, and a set of 59,792 high-quality single nucleotide polymorphisms (SNPs) distributed across the cassava genome. Through genome-wide association mapping, we identified 80 SNPs that were significantly associated with root mealiness, fibre, adhesiveness, taste, aroma, colour and firmness on chromosomes 1, 4, 5, 6, 10, 13, 17 and 18. We also identified relevant candidate genes that are co-located with peak SNPs linked to these traits in M. esculenta. A survey of the cassava reference genome v6.1 positioned the SNPs on chromosome 13 in the vicinity of Manes.13G026900, a gene recognized as being responsible for cell adhesion and for the mealiness or crispness of vegetables and fruits, and also known to play an important role in cooked potato texture. This study provides the first insights into understanding the underlying genetic basis of boiled cassava root texture. After validation, the markers and candidate genes identified in this novel work could provide important genomic resources for use in marker-assisted selection (MAS) and genomic selection (GS) to accelerate genetic improvement of root mealiness and other culinary qualities in cassava breeding programmes in West Africa, especially in Nigeria, where the consumption of boiled and pounded cassava is low.

Functional phenomics and genetics of the root economics space in winter wheat using high-throughput phenotyping of respiration and architecture

The root economics space is a useful framework for plant ecology but is rarely considered for crop ecophysiology. In order to understand root trait integration in winter wheat, we combined functional phenomics with trait economic theory, utilizing genetic variation, high-throughput phenotyping, and multivariate analyses. We phenotyped a diversity panel of 276 genotypes for root respiration and architectural traits using a novel high-throughput method for CO2 flux and the open-source software RhizoVision Explorer to analyze scanned images. We uncovered substantial variation in specific root respiration (SRR) and specific root length (SRL), which were primary indicators of root metabolic and structural costs. Multiple linear regression analysis indicated that lateral root tips had the greatest SRR, and the residuals from this model were used as a new trait. Specific root respiration was negatively correlated with plant mass. Network analysis, using a Gaussian graphical model, identified root weight, SRL, diameter, and SRR as hub traits. Univariate and multivariate genetic analyses identified genetic regions associated with SRR, SRL, and root branching frequency, and proposed gene candidates. Combining functional phenomics and root economics is a promising approach to improving our understanding of crop ecophysiology. We identified root traits and genomic regions that could be harnessed to breed more efficient crops for sustainable agroecosystems.

Functional phenomics and genetics of the root economics space in winter wheat using high-throughput phenotyping of respiration and architecture

The root economics space is a useful framework for plant ecology but is rarely considered for crop ecophysiology. In order to understand root trait integration in winter wheat, we combined functional phenomics with trait economic theory, utilizing genetic variation, high-throughput phenotyping, and multivariate analyses. We phenotyped a diversity panel of 276 genotypes for root respiration and architectural traits using a novel high-throughput method for CO₂ flux and the open-source software RhizoVision Explorer to analyze scanned images. We uncovered substantial variation in specific root respiration (SRR) and specific root length (SRL), which were primary indicators of root metabolic and structural costs. Multiple linear regression analysis indicated that lateral root tips had the greatest SRR, and the residuals from this model were used as a new trait. Specific root respiration was negatively correlated with plant mass. Network analysis, using a Gaussian graphical model, identified root weight, SRL, diameter, and SRR as hub traits. Univariate and multivariate genetic analyses identified genetic regions associated with SRR, SRL, and root branching frequency, and proposed gene candidates. Combining functional phenomics and root economics is a promising approach to improving our understanding of crop ecophysiology. We identified root traits and genomic regions that could be harnessed to breed more efficient crops for sustainable agroecosystems.

Identification of Quantitative Trait Loci Hotspots Affecting Agronomic Traits and High-Throughput Vegetation Indices in Rainfed Wheat

Understanding the genetic basis of agronomic traits is essential for wheat breeding programs to develop new cultivars with enhanced grain yield under climate change conditions. The use of high-throughput phenotyping (HTP) technologies for the assessment of agronomic performance through drought-adaptive traits opens new possibilities in plant breeding. HTP together with a genome-wide association study (GWAS) mapping approach can be a useful method to dissect the genetic control of complex traits in wheat to enhance grain yield under drought stress. This study aimed to identify molecular markers associated with agronomic and remotely sensed vegetation index (VI)-related traits under rainfed conditions in bread wheat and to use an in silico candidate gene (CG) approach to search for upregulated CGs under abiotic stress. The plant material consisted of 170 landraces and 184 modern cultivars from the Mediterranean basin. The collection was phenotyped for agronomic and VI traits derived from multispectral images over 3 and 2 years, respectively. The GWAS identified 2,579 marker-trait associations (MTAs). The quantitative trait loci (QTL) overview index statistic detected 11 QTL hotspots involving more than one trait in at least 2 years. A CG analysis detected 12 CGs upregulated under abiotic stress in six QTL hotspots and 46 downregulated CGs in 10 QTL hotspots. The current study highlights the utility of VI to identify chromosome regions that contribute to yield and drought tolerance under rainfed Mediterranean conditions.

Dissection of quantitative trait loci in the Lachancea waltii yeast species highlights major hotspots

Dissecting the genetic basis of complex trait remains a real challenge. The budding yeast Saccharomyces cerevisiae has become a model organism for studying quantitative traits, successfully increasing our knowledge in many aspects. However, the exploration of the genotype-phenotype relationship in non-model yeast species could provide a deeper insight into the genetic basis of complex traits. Here, we have studied this relationship in the Lachancea waltii species which diverged from the S. cerevisiae lineage prior to the whole-genome duplication. By performing linkage mapping analyses in this species, we identified 86 quantitative trait loci (QTL) impacting the growth in a large number of conditions. The distribution of these loci across the genome has revealed two major QTL hotspots. A first hotspot corresponds to a general growth QTL, impacting a wide range of conditions. By contrast, the second hotspot highlighted a trade-off with a disadvantageous allele for drug-free conditions which proved to be advantageous in the presence of several drugs. Finally, a comparison of the detected QTL in L. waltii with those which had been previously identified for the same trait in a closely related species, Lachancea kluyveri was performed. This analysis clearly showed the absence of shared QTL across these species. Altogether, our results represent a first step toward the exploration of the genetic architecture of quantitative trait across different yeast species.

Meta-analysis of genome-wide association studies reveal common loci controlling agronomic and quality traits in a wide range of normal and heat stressed environments

Several stable QTL were detected using metaGWAS analysis for different agronomic and quality traits under 26 normal and heat stressed environments. Heat stress, exacerbated by global warming, has a negative influence on wheat production worldwide and climate resilient cultivars can help mitigate these impacts. Selection decisions should therefore depend on multi-environment experiments representing a range of temperatures at critical stages of development. Here, we applied a meta-genome wide association analysis (metaGWAS) approach to detect stable QTL with significant effects across multiple environments. The metaGWAS was applied to 11 traits scored in 26 trials that were sown at optimal or late times of sowing (TOS1 and TOS2, respectively) at five locations. A total of 2571 unique wheat genotypes (13,959 genotypes across all environments) were included and the analysis conducted on TOS1, TOS2 and both times of sowing combined (TOS1&2). The germplasm was genotyped using a 90 k Infinium chip and imputed to exome sequence level, resulting in 341,195 single nucleotide polymorphisms (SNPs). The average accuracy across all imputed SNPs was high (92.4%). The three metaGWAS analyses revealed 107 QTL for the 11 traits, of which 16 were detected in all three analyses and 23 were detected in TOS1&2 only. The remaining QTL were detected in either TOS1 or TOS2 with or without TOS1&2, reflecting the complex interactions between the environments and the detected QTL. Eight QTL were associated with grain yield and seven with multiple traits. The identified QTL provide an important resource for gene enrichment and fine mapping to further understand the mechanisms of gene × environment interaction under both heat stressed and unstressed conditions.

Dissection of Genetic Basis Underpinning Kernel Weight-Related Traits in Common Wheat

Genetic dissection kernel weight-related traits is of great significance for improving wheat yield potential. As one of the three major yield components of wheat, thousand kernel weight (TKW) was mainly affected by grain length (GL) and grain width (GW). To uncover the key loci for these traits, we carried out a quantitative trait loci (QTL) analysis of an F₆ recombinant inbred lines (RILs) population derived from a cross of Henong 5290 (small grain) and 06Dn23 (big grain) with a 50 K single nucleotide polymorphism (SNP) array. A total of 17 stable and big effect QTL, including 5 for TKW, 8 for GL and 4 for GW, were detected on the chromosomes 1B, 2A, 2B, 2D, 4B, 5A, 6A and 6D, respectively. Among these, there were two co-located loci for three traits that were mapped on the chromosome 4BS and 6AL. The QTL on 6AL was the most stable locus and explained 15.4–24.8%, 4.1–8.8% and 15.7–24.4% of TKW, GW and GL variance, respectively. In addition, two more major QTL of GL were located on chromosome arm 2BL and 2DL, accounting for 9.7–17.8% and 13.6–19.8% of phenotypic variance, respectively. In this study, we found one novel co-located QTL associated with GL and TKW in 2DL, QGl.haaf-2DL.2/QTkw.haaf-2DL.2, which could explain 13.6–19.8% and 9.8–10.7% phenotypic variance, respectively. Genetic regions and linked markers of these stable QTL will help to further refine mapping of the corresponding loci and marker-assisted selection (MAS) breeding for wheat grain yield potential improvement.

Genome-wide association study of multiple traits linked to heat tolerance in emmer-derived hexaploid wheat genotypes

Heat stress tolerance in plants is a complex trait controlled by multiple genes of minor effect which are influenced by the environment and this makes breeding and selection complicated. Emmer wheat (Triticum dicoccon Schrank) carries valuable diversity that can be used to improve the heat tolerance of modern bread wheat. A diverse set of emmer-based genotypes was developed by crossing emmer wheat with hexaploid wheat. These materials, along with their hexaploid recurrent parents and commercial cultivars, were evaluated at optimum (E1) and heat stressed (E2) sowing times in the field for three consecutive years (2014-2016). The material was genotyped using the Infinium iSelect SNP 90K SNP Assay. The phenotypic data were combined across years within each sowing time and best linear unbiased estimators calculated for each genotype in each environment. These estimates were used for GWAS analysis. Significant phenotypic and genotypic variation was observed for all traits. A total of 125 and 142 marker-trait associations (MTAs) were identified in E1 and E2, respectively. The highest number of MTAs were observed on the A genome (106), followed by the B (105) and D (56) genomes. MTAs with pleiotropic effects within and across the environments were observed. Many of the MTAs found were reported previously for various traits, and a few significant MTAs under heat stress were new and linked to emmer genome. Genomic regions identified on chromosomes 2B and 3A had a significant positive impact on grain yield under stress with a 7% allelic effect. Genomic regions on chromosomes 1A and 4B contributed 11% and 9% of the variation for thousand kernel weight (TKW) under heat stress respectively. Following fine mapping, these regions could be used for marker-assisted selection to improve heat tolerance in wheat.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01222-3.

Labelling Selective Sweeps Used in Durum Wheat Breeding from a Diverse and Structured Panel of Landraces and Cultivars

Simple Summary

Evaluation of the genetic diversity of a crop species is a critical step for breeding. Landraces are essential to avoid genetic erosion, and Mediterranean landraces are an important group of genetic resources due to their high genetic variability, adaptation to local conditions in rainfed environments, and their resilience to pests and pathogens. This study uses a genome-wide association approach employing eigenvectors to identify selective sweeps among Mediterranean durum wheat landraces and a world panel of modern cultivars.

Abstract

A panel of 387 durum wheat genotypes including Mediterranean landraces and modern cultivars was characterized with 46,161 diversity arrays technology (DArTseq) markers. Analysis of population structure uncovered the existence of five subpopulations (SP) related to the pattern of migration of durum wheat from the domestication area to the west of the Mediterranean basin (SPs 1, 2, and 3) and further improved germplasm (SPs 4 and 5). The total genetic diversity (H_T) was 0.40 with a genetic differentiation (G_ST) of 0.08 and a mean gene flow among SPs of 6.02. The lowest gene flow was detected between SP 1 (presumably the ancient genetic pool of the panel) and SPs 4 and 5. However, gene flow from SP 2 to modern cultivars was much higher. The highest gene flow was detected between SP 3 (western Mediterranean germplasm) and SP 5 (North American and European cultivars). A genome wide association study (GWAS) approach using the top ten eigenvectors as phenotypic data revealed the presence of 89 selective sweeps, represented as quantitative trait loci (QTL) hotspots, widely distributed across the durum wheat genome. A principal component analysis (PCoA) using 147 markers with −log₁₀p > 5 identified three regions located on chromosomes 2A, 2B and 3A as the main drivers for differentiation of Mediterranean landraces. Gene flow between SPs offers clues regarding the putative use of Mediterranean old durum germplasm by the breeding programs represented in the structure analysis. EigenGWAS identified selective sweeps among landraces and modern cultivars. The analysis of the corresponding genomic regions in the ‘Zavitan’, ‘Svevo’ and ‘Chinese Spring’ genomes discovered the presence of important functional genes including Ppd, Vrn, Rht, and gene models involved in important biological processes including LRR-RLK, MADS-box, NAC, and F-box.

QTL mapping of PEG-induced drought tolerance at the early seedling stage in sesame using whole genome re-sequencing

Improvement in sesame drought tolerance at seedling stage is important for yield stability. Genetic approaches combing with conventional breeding is the most effective way to develop drought-tolerant cultivars. In this study, three traits and their relative values, including seedling weight (SW), shoot length (SL) and root length (RL), were evaluated under control and osmotic conditions in a recombinant inbred line (RIL) population derived from cross of Zhushanbai and Jinhuangma. Significant variation and high broad sense heritability were observed for all traits except SW under stress condition in the population. With this population, a high-density linkage map with 1354 bin markers was constructed through whole genome re-sequencing (WGS) strategy. Quantitative trait loci (QTL) mapping was performed for all the traits. A total of 34 QTLs were detected on 10 chromosomes. Among them, 13 stable QTLs were revealed in two independent experiments, eight of them were associated with traits under water stress condition. One region on chromosome 12 related to RL under osmotic condition and relative RL had the highest LOD value and explained the largest phenotypic variation among all the QTLs detected under water stress condition. These findings will provide new genetic resources for molecular improvement of drought tolerance and candidate gene identification in sesame.

Genome-wide association of rice response to blast fungus identifies loci for robust resistance under high nitrogen

BACKGROUND

Nitrogen fertilization is known to increase disease susceptibility, a phenomenon called Nitrogen-Induced Susceptibility (NIS). In rice, this phenomenon has been observed in infections with the blast fungus Magnaporthe oryzae. A previous classical genetic study revealed a locus (NIS1) that enhances susceptibility to rice blast under high nitrogen fertilization. In order to further address the underlying genetics of plasticity in susceptibility to rice blast after fertilization, we analyzed NIS under greenhouse-controlled conditions in a panel of 139 temperate japonica rice strains. A genome-wide association analysis was conducted to identify loci potentially involved in NIS by comparing susceptibility loci identified under high and low nitrogen conditions, an approach allowing for the identification of loci validated across different nitrogen environments. We also used a novel NIS Index to identify loci potentially contributing to plasticity in susceptibility under different nitrogen fertilization regimes.

RESULTS

A global NIS effect was observed in the population, with the density of lesions increasing by 8%, on average, under high nitrogen fertilization. Three new QTL, other than NIS1, were identified. A rare allele of the RRobN1 locus on chromosome 6 provides robust resistance in high and low nitrogen environments. A frequent allele of the NIS2 locus, on chromosome 5, exacerbates blast susceptibility under the high nitrogen condition. Finally, an allele of NIS3, on chromosome 10, buffers the increase of susceptibility arising from nitrogen fertilization but increases global levels of susceptibility. This allele is almost fixed in temperate japonicas, as a probable consequence of genetic hitchhiking with a locus involved in cold stress adaptation.

CONCLUSIONS

Our results extend to an entire rice subspecies the initial finding that nitrogen increases rice blast susceptibility. We demonstrate the usefulness of estimating plasticity for the identification of novel loci involved in the response of rice to the blast fungus under different nitrogen regimes.

Genome-wide association study of yield and related traits in common wheat under salt-stress conditions

BACKGROUND

Soil salinization is a major threat to wheat production. It is essential to understand the genetic basis of salt tolerance for breeding and selecting new salt-tolerant cultivars that have the potential to increase wheat yield.

RESULT

In this study, a panel of 191 wheat accessions was subjected to genome wide association study (GWAS) to identify SNP markers linked with adult-stage characters. The population was genotyped by Wheat660K SNP array and eight phenotype traits were investigated under low and high salinity environments for three consecutive years. A total of 389 SNPs representing 11 QTLs were significantly associated with plant height, spike number, spike length, grain number, thousand kernels weight, yield and biological mass under different salt treatments, with the phenotypic explanation rate (R²) ranging from 9.14 to 50.45%. Of these, repetitive and pleiotropic loci on chromosomes 4A, 5A, 5B and 7A were significantly linked to yield and yield related traits such as thousand kernels weight, spike number, spike length, grain number and so on under low salinity conditions. Spike length-related loci were mainly located on chromosomes 1B, 3B, 5B and 7A under different salt treatments. Two loci on chromosome 4D and 5A were related with plant height in low and high salinity environment, respectively. Three salt-tolerant related loci were confirmed to be important in two bi-parental populations. Distribution of favorable haplotypes indicated that superior haplotypes of pleiotropic loci on group-5 chromosomes were strongly selected and had potential for increasing wheat salt tolerance. A total of 14 KASP markers were developed for nine loci associating with yield and related traits to improve the selection efficiency of wheat salt-tolerance breeding.

CONCLUSION

Utilizing a Wheat660K SNPs chip, QTLs for yield and its related traits were detected under salt treatment in a natural wheat population. Important salt-tolerant related loci were validated in RIL and DH populations. This study provided reliable molecular markers that could be crucial for marker-assisted selection in wheat salt tolerance breeding programs.

Breeding Driven Enrichment of Genetic Variation for Key Yield Components and Grain Starch Content Under Drought Stress in Winter Wheat

Drought is one of the major abiotic stress factors limiting wheat production worldwide, thus threatening food security. The dissection of the genetic footprint of drought stress response offers strong opportunities toward understanding and improving drought tolerance (DT) in wheat. In this study, we investigated the genotypic variability for drought response among 200 diverse wheat cultivars (genotypes) using agronomic, developmental, and grain quality traits (GQT), and conducted genome-wide association studies (GWAS) to uncover the genetic architectures of these important traits. Results indicated significant effects of genotype, water regime and their interactions for all agronomic traits. Grain yield (GY) was the most drought-responsive trait and was highly correlated with kernels number per meter square (KN). Genome-wide association studies revealed 17 and 20 QTL regions under rainfed and drought conditions, respectively, and identified one LD block on chromosome 3A and two others on 5D associated with breeding progress (BP). The major haplotypes of these LD blocks have been positively selected through breeding and are associated with higher starch accumulation and GY under drought conditions. Upon validation, the identified QTL regions caring favorable alleles for high starch and yield will shed light on mechanisms of tolerance to drought and can be used to develop drought resistant cultivars.

Genome-wide association study identifies QTL for thousand grain weight in winter wheat under normal- and late-sown stressed environments

GWAS identified stable loci for TGW and stress tolerance in winter wheat based on two sowing conditions, which will provide opportunities for developing new cultivars with high yield and yield stability. Wheat is an important food crop widely cultivated in the world. Breeding new varieties with high yields and superior adaptability is the main goal of modern wheat breeding program. In order to determine the marker-trait associations (MATs), a set of 688 diverse winter wheat accessions were subjected to genome-wide association study (GWAS) using the wheat 90K array. Field trials under normal-sown (NS) and late-sown (LS) conditions were conducted for thousand grain weight (TGW) and stress susceptibility index (SSI) at three different sites across two consecutive years. A total of 179 (NS) and 158 (LS) MATs corresponded with TGW; of these, 16 and 6 stable MATs for TGW_NS and TGW_LS were identified on chromosomes 1B, 2B, 3A, 3B, 5A, 5B, 5D, 6B, and 7D across at least three environments. Notably, a QTL hot spot controlling TGW under NS and LS conditions was found on chromosome 5A (140-142 cM). Moreover, 8 of 228 stable MATs on chromosomes 4B, 5A, and 5D for SSI were detected. A haplotype block associated with TGW and SSI was located on chromosome 5A at 91 cM, nearby the vernalization gene VRN-A1. Additionally, analysis of wheat varieties from the different eras revealed that the grain weight and stress tolerance are not improved concurrently. Overall, our results provide promising alleles controlling grain weight and stress tolerance (particularly for thermotolerance) for wheat breeders, which can be used in marker-assisted selection for improving grain yield and yield stability in wheat.

Utilization of genetic diversity and marker-trait to improve drought tolerance in rice (Oryza sativa L.)

Drought stress is one of the main problems for the rice crop, as it reduces the production and productivity of the grain yield significantly. In Egypt, many restrictions were made on the cultivation of rice due to its high-water demand. Producing promising drought-tolerant rice cultivars combined with high yielding is one of the main targets for rice breeders. A set of 22 highly diverse rice genotypes were evaluated under normal and drought conditions. Morphological, physiological, and yield traits were recorded on each genotype. Drought susceptibility index (DSI) was estimated for six yield traits to identify the most drought-tolerant rice genotypes. A high genetic variation was found among genotypes tested in the experiment. Under normal conditions, the highest phenotypic correlation was found between grain yield (GY) and sterility percentage (SP) (- 0.73**), while it was among GY and chlorophyll content (CC) (0.82**) under drought conditions. To identify quantitative trait loci (QTL) controlling yielding traits under drought and normal, a single marker analysis was performed between all yield traits under both conditions and a set of 106 simple sequence repeat (SSR) marker alleles. The genetic association analysis revealed 14 and 17 QTL under drought and normal conditions, respectively. The most drought-tolerant genotypes were selected based on phenotypic traits, the number of QTL in each selected genotype, and the level of genetic diversity existed among the genotypes. As a result, five genotypes (Giza 178, IET1444, GZ1368-S-5-4, Nahda, Giza 14) were identified as the most promising drought-tolerant rice genotypes. Eight QTL controlling drought tolerance were identified in Giza 178, Nahda, and GZ1368-S-5-4, while four QTL were found in IET1444. The number of different QTLs were estimated among the five selected genotypes. Giza 178 and GZ1368-S-5-4 shared the same QTLs. Seven different QTLs were found among Nahda, IET1444, GZ1368-S-5-4, and Giza 14. Combining information from phenotypic traits, genetic diversity analysis, and QTL analysis was very useful in identifying the true drought-tolerant rice genotypes that can be used for crossing in the future breeding program.

Transcripts of wheat at a target locus on chromosome 6B associated with increased yield, leaf mass and chlorophyll index under combined drought and heat stress

Drought and heat stress constrain wheat (Triticum aestivum L.) yields globally. To identify putative mechanisms and candidate genes associated with combined drought and heat stress tolerance, we developed bread wheat near-isogenic lines (NILs) targeting a quantitative trait locus (QTL) on chromosome 6B which was previously associated with combined drought and heat stress tolerance in a diverse panel of wheats. Genotyping-by-sequencing was used to identify additional regions that segregated in allelic pairs between the recurrent and the introduced exotic parent, genome-wide. NILs were phenotyped in a gravimetric platform with precision irrigation and exposed to either drought or to combined drought and heat stress from three days after anthesis. An increase in grain weight in NILs carrying the exotic allele at 6B locus was associated with thicker, greener leaves, higher photosynthetic capacity and increased water use index after re-watering. RNA sequencing of developing grains at early and later stages of treatment revealed 75 genes that were differentially expressed between NILs across both treatments and timepoints. Differentially expressed genes coincided with the targeted QTL on chromosome 6B and regions of genetic segregation on chromosomes 1B and 7A. Pathway enrichment analysis showed the involvement of these genes in cell and gene regulation, metabolism of amino acids and transport of carbohydrates. The majority of these genes have not been characterized previously under drought or heat stress and they might serve as candidate genes for improved abiotic stress tolerance.

Genome-Wide Association Mapping of Seedling Drought Tolerance in Winter Wheat

In the southern Great Plains of the United States, winter wheat grown for dual-purpose is often planted early, which puts it at risk for drought stress at the seedling stage in the autumn. To map quantitative trait loci (QTL) associated with seedling drought tolerance, a genome-wide association study (GWAS) was performed on a hard winter wheat association mapping panel. Two sets of plants were planted in the greenhouse initially under well-watered conditions. At the five-leaf stage, one set continued to receive the optimum amount of water, whereas watering was withdrawn from the other set (drought stress treatment) for 14 days to mimic drought stress. Large phenotypic variation was observed in leaf chlorophyll content, leaf chlorophyll fluorescence, shoot length, number of leaves per seedling, and seedling recovery. A mixed linear model analysis detected multiple significant QTL associated with seedling drought tolerance-related traits on chromosomes 1B, 2A, 2B, 2D, 3A, 3B, 3D, 4B, 5A, 5B, 6B, and 7B. Among those, 12 stable QTL responding to drought stress for various traits were identified. Shoot length and leaf chlorophyll fluorescence were good indicators in responding to drought stress because most of the drought responding QTL detected using means of these two traits were also detected in at least two experimental repeats. These stable QTL are more valuable for use in marker-assisted selection during wheat breeding. Moreover, different traits were mapped on several common chromosomes, such as 1B, 2B, 3B, and 6B, and two QTL clusters associated with three or more traits were located at 107-130 and 80-83 cM on chromosomes 2B and 6B, respectively. Furthermore, some QTL detected in this study co-localized with previously reported QTL for root and shoot traits at the seedling stage and canopy temperature at the grain-filling stage of wheat. In addition, several of the mapped chromosomes were also associated with drought tolerance during the flowering or grain-filling stage in wheat. Some significant single-nucleotide polymorphisms (SNPs) were aligned to candidate genes playing roles in plant abiotic stress responses. The SNP markers identified in this study will be further validated and used for marker-assisted breeding of seedling drought tolerance during dual-purpose wheat breeding.

Molecular Mapping and Genomics of Grain Yield in Durum Wheat: A Review

Durum wheat is the most relevant cereal for the whole of Mediterranean agriculture, due to its intrinsic adaptation to dryland and semi-arid environments and to its strong historical cultivation tradition. It is not only relevant for the primary production sector, but also for the food industry chains associated with it. In Mediterranean environments, wheat is mostly grown under rainfed conditions and the crop is frequently exposed to environmental stresses, with high temperatures and water scarcity especially during the grain filling period. For these reasons, and due to recurrent disease epidemics, Mediterranean wheat productivity often remains under potential levels. Many studies, using both linkage analysis (LA) and a genome-wide association study (GWAS), have identified the genomic regions controlling the grain yield and the associated markers that can be used for marker-assisted selection (MAS) programs. Here, we have summarized all the current studies identifying quantitative trait loci (QTLs) and/or candidate genes involved in the main traits linked to grain yield: kernel weight, number of kernels per spike and number of spikes per unit area.

Quantitative trait loci for agronomic traits in tetraploid wheat for enhancing grain yield in Kazakhstan environments

Durum wheat (Triticum turgidum L. ssp. durum) is one of the top crops in Kazakhstan, where it is cultivated in different ecological niches, mainly at higher latitudes in the steppe zone of the northern region. Therefore, local breeding programs for durum wheat are primarily focused on selection for high productivity in Northern Kazakhstan based on the introduction of promising foreign germplasm and the adoption of marker-assisted selection. In this study, a world tetraploid wheat collection consisted of 184 primitive and domesticated accessions, which were previously genotyped using 16,425 polymorphic SNP markers, was field-tested in Northern and South-eastern Kazakhstan. The field tests have allowed the identification of 80 durum wheat promising lines in Northern Kazakhstan in comparison with a local standard cultivar. Also, GGE (Genotype and Genotype by Environment) biplot analyses for yield performance revealed that accessions of T. dicoccum, T. carthlicum, and T. turanicum also have potential to improve durum wheat yield in the region. The genome-wide association study (GWAS) has allowed the identification of 83 MTAs (marker-trait associations) for heading date, seed maturation time, plant height, spike length, number of fertile spikes, number of kernels per spike, and thousand kernel weight. The comparison of the 83 identified MTAs with those previously reported in GWAS for durum wheat suggests that 38 MTAs are presumably novel, while the co-localization of a large number of MTAs with those previously published confirms the validity of the results of this study. The MTAs reported herewith will provide the opportunity to implement marker-assisted selection in ongoing durum wheat breeding projects targeting higher productivity in the region.

Genome-Wide Association Studies in Diverse Spring Wheat Panel for Stripe, Stem, and Leaf Rust Resistance

Among several important wheat foliar diseases, Stripe rust (YR), Leaf rust (LR), and Stem rust (SR) have always been an issue of concern to the farmers and wheat breeders. Evolution of virulent pathotypes of these rusts has posed frequent threats to an epidemic. Pyramiding rust-resistant genes are the most economical and environment-friendly approach in postponing this inevitable threat. To achieve durable long term resistance against the three rusts, an attempt in this study was made searching for novel sources of resistant alleles in a panel of 483 spring wheat genotypes. This is a unique and comprehensive study where evaluation of a diverse panel comprising wheat germplasm from various categories and adapted to different wheat agro-climatic zones was challenged with 18 pathotypes of the three rusts with simultaneous screening in field conditions. The panel was genotyped using 35K SNP array and evaluated for each rust at two locations for two consecutive crop seasons. High heritability estimates of disease response were observed between environments for each rust type. A significant effect of population structure in the panel was visible in the disease response. Using a compressed mixed linear model approach, 25 genomic regions were found associated with resistance for at least two rusts. Out of these, seven were associated with all the three rusts on chromosome groups 1 and 6 along with 2B. For resistance against YR, LR, and SR, there were 16, 18, and 27 QTL (quantitative trait loci) identified respectively, associated at least in two out of four environments. Several of these regions got annotated with resistance associated genes viz. NB-LRR, E3-ubiquitin protein ligase, ABC transporter protein, etc. Alien introgressed (on 1B and 3D) and pleiotropic (on 7D) resistance genes were captured in seedling and adult plant disease responses, respectively. The present study demonstrates the use of genome-wide association for identification of a large number of favorable alleles for leaf, stripe, and stem rust resistance for broadening the genetic base. Quick conversion of these QTL into user-friendly markers will accelerate the deployment of these resistance loci in wheat breeding programs.

Genome-wide association study of Striga resistance in early maturing white tropical maize inbred lines

BACKGROUND

Striga hermonthica (Benth.) parasitism militates against increased maize production and productivity in savannas of sub-Saharan Africa (SSA). Identification of Striga resistance genes is important in developing genotypes with durable resistance. So far, there is only one report on the existence of QTL for Striga resistance on chromosome 6 of maize. The objective of this study was to identify genomic regions significantly associated with grain yield and other agronomic traits under artificial Striga field infestation. A panel of 132 early-maturing maize inbreds were phenotyped for key agronomic traits under Striga-infested and Striga-free conditions. The inbred lines were also genotyped using 47,440 DArTseq markers from which 7224 markers were retained for population structure analysis and genome-wide association study (GWAS).

RESULTS

The inbred lines were grouped into two major clusters based on structure analysis as well as the neighbor-joining hierarchical clustering. A total of 24 SNPs significantly associated with grain yield, Striga damage at 8 and 10 weeks after planting (WAP), ears per plant and ear aspect under Striga infestation were detected. Under Striga-free conditions, 11 SNPs significantly associated with grain yield, number of ears per plant and ear aspect were identified. Three markers physically located close to the putative genes GRMZM2G164743 (bin 10.05), GRMZM2G060216 (bin 3.06) and GRMZM2G103085 (bin 5.07) were detected, linked to grain yield, Striga damage at 8 and 10 WAP and number of ears per plant under Striga infestation, explaining 9 to 42% of the phenotypic variance. Furthermore, the S9_154,978,426 locus on chromosome 9 was found at 2.61 Mb close to the ZmCCD1 gene known to be associated with the reduction of strigolactone production in the maize roots.

CONCLUSIONS

Presented in this study is the first report of the identification of significant loci on chromosomes 9 and 10 of maize that are closely linked to ZmCCD1 and amt5 genes, respectively and may be related to plant defense mechanisms against Striga parasitism. After validation, the identified loci could be targets for breeders for marker-assisted selection (MAS) to accelerate genetic enhancement of maize for Striga resistance in the tropics, particularly in SSA, where the parasitic weed is endemic.

Genome-wide association study of agronomic and quality traits in a world collection of the wild wheat relative Triticum urartu

Triticum urartu (2n = 2x = 14, subgenome A^u A^u ), a wild diploid wheat progenitor, features broad allelic diversity for a number of traits of agronomic relevance. A thorough characterization of the diversity of T. urartu natural accessions may provide wheat breeders with new alleles potentially contributing to wheat improvement. In this study, we performed an extensive genotypic and phenotypic characterization of a world collection of 299 T. urartu ex situ accessions, developing 441 327 single nucleotide polymorphisms and recording trait values for agronomic and quality traits. The collection was highly diverse, with broad variation in phenology and plant architecture traits. Seed features were also varied, and analyses of flour quality reported 18 distinct patterns of glutenins, and carotenoid concentrations and sedimentation volumes in some cases surpassing those of cultivated materials. The genome-wide molecular markers developed on the collection were used to conduct a genome-wide association study reporting 25 highly significant quantitative trait nucleotides for the traits under examination, only partially overlapping loci already reported in wheat. Our data show that T. urartu may be considered a valuable allele pool to support the improvement of wheat agronomy and quality.

Genetic dissection of heat-responsive physiological traits to improve adaptation and increase yield potential in soft winter wheat

BACKGROUND

Climate change, including higher temperatures (HT) has a detrimental impact on wheat productivity and modeling studies predict more frequent heat waves in the future. Wheat growth can be impaired by high daytime and nighttime temperature at any developmental stage, especially during the grain filling stage. Leaf chlorophyll content, leaf greenness, cell membrane thermostability, and canopy temperature have been proposed as candidate traits to improve crop adaptation and yield potential of wheat under HT. Nonetheless, a significant gap exists in knowledge of genetic backgrounds associated with these physiological traits. Identifying genetic loci associated with these traits can facilitate physiological breeding for increased yield potential under high temperature stress condition in wheat.

RESULTS

We conducted genome-wide association study (GWAS) on a 236 elite soft wheat association mapping panel using 27,466 high quality single nucleotide polymorphism markers. The panel was phenotyped for three years in two locations where heat shock was common. GWAS identified 500 significant marker-trait associations (MTAs) (p ≤ 9.99 × 10- 4). Ten MTAs with pleiotropic effects detected on chromosomes 1D, 2B, 3A, 3B, 6A, 7B, and 7D are potentially important targets for selection. Five MTAs associated with physiological traits had pleiotropic effects on grain yield and yield-related traits. Seventy-five MTAs were consistently expressed over several environments indicating stability and more than half of these stable MTAs were found in genes encoding different types of proteins associated with heat stress.

CONCLUSIONS

We identified 500 significant MTAs in soft winter wheat under HT stress. We found several stable loci across environments and pleiotropic markers controlling physiological and agronomic traits. After further validation, these MTAs can be used in marker-assisted selection and breeding to develop varieties with high stability for grain yield under high temperature.

Incorporating Genome-Wide Association Mapping Results Into Genomic Prediction Models for Grain Yield and Yield Stability in CIMMYT Spring Bread Wheat

Untangling the genetic architecture of grain yield (GY) and yield stability is an important determining factor to optimize genomics-assisted selection strategies in wheat. We conducted in-depth investigation on the above using a large set of advanced bread wheat lines (4,302), which were genotyped with genotyping-by-sequencing markers and phenotyped under contrasting (irrigated and stress) environments. Haplotypes-based genome-wide-association study (GWAS) identified 58 associations with GY and 15 with superiority index Pi (measure of stability). Sixteen associations with GY were "environment-specific" with two on chromosomes 3B and 6B with the large effects and 8 associations were consistent across environments and trials. For Pi, 8 associations were from chromosomes 4B and 7B, indicating 'hot spot' regions for stability. Epistatic interactions contributed to an additional 5-9% variation on average. We further explored whether integrating consistent and robust associations identified in GWAS as fixed effects in prediction models improves prediction accuracy. For GY, the model accounting for the haplotype-based GWAS loci as fixed effects led to up to 9-10% increase in prediction accuracy, whereas for Pi this approach did not provide any advantage. This is the first report of integrating genetic architecture of GY and yield stability into prediction models in wheat.

Multi-Locus Genome Wide Association Mapping for Yield and Its Contributing Traits in Hexaploid Wheat under Different Water Regimes

Multi-locus genome wide association study was undertaken using a set of 320 diverse spring wheat accessions, which were each genotyped for 9,626 SNPs. The association panel was grown in replicated trials in four environments [two each in irrigated (IR) and rainfed (RF) environments], and phenotypic data were recorded for five traits including days to heading, days to maturity, plant height, thousand grain weight and grain yield. Forty-six significant marker-trait associations (MTAs) were identified for five traits. These included 20 MTAs in IR and 19 MTAs in RF environments; seven additional MTAs were common to both the environments. Five of these MTAs were co-localized with previously known QTL/MTAs and the remaining MTAs were novel and add to the existing knowledge. Three desirable haplotypes for agronomic traits, one for improvement in RF environment and two for improvement in IR environment were identified. Eighteen (18) promising candidate genes (CGs) involved in seven different biological activities were also identified. The expression profiles of four (Trehalose-6-Phosphate, APETALA2/Ethylene-responsive factor, DNA-binding One Zinc Finger and Gibberellin-dioxygenases) of the 18 genes showed that they were induced by drought stress in the wheat seedlings. The MTAs, haplotypes and CG-based markers may be used in marker-assisted breeding for drought tolerance in wheat.

Genome-wide association study of drought tolerance and biomass allocation in wheat

Genome wide association studies (GWAS) are important in discerning the genetic architecture of complex traits such as biomass allocation for improving drought tolerance and carbon sequestration potential of wheat. The objectives of this study were to deduce the population structure and marker-trait association for biomass traits in wheat under drought-stressed and non-stressed conditions. A 100-wheat (Triticum aestivum L.) genotype panel was phenotyped for days to heading (DTH), days to maturity (DTM), shoot biomass (SB), root biomass (RB), root to shoot ratio (RS) and grain yield (GY). The panel was sequenced using 15,600 single nucleotide polymorphism (SNPs) markers and subjected to genetic analysis using the compressed mixed linear model (CMLM) at false discovery rate (FDR < 0.05). Population structure analysis revealed six sub-clusters with high membership ancestry coefficient of ≤0.65 to their assigned sub-clusters. A total of 75 significant marker-trait associations (MTAs) were identified with a linkage disequilibrium threshold of 0.38 at 5cM. Thirty-seven of the MTAs were detected under drought-stressed condition and 48% were on the B genome, where most quantitative trait loci (QTLs) for RB, SB and GY were previously identified. There were seven pleiotropic markers for RB and SB that may facilitate simultaneous selection. Thirty-seven putative candidate genes were mined by gene annotation on the IWGSC RefSeq 1.1. The significant MTAs observed in this study will be useful in devising strategies for marker-assisted breeding for simultaneous improvement of drought tolerance and to enhance C sequestration capacity of wheat.

Understanding the Genetic Basis of Spike Fertility to Improve Grain Number, Harvest Index, and Grain Yield in Wheat Under High Temperature Stress Environments

Moderate heat stress accompanied by short episodes of extreme heat during the post-anthesis stage is common in most US wheat growing areas and causes substantial yield losses. Sink strength (grain number) is a key yield limiting factor in modern wheat varieties. Increasing spike fertility (SF) and improving the partitioning of assimilates can optimize sink strength which is essential to improve wheat yield potential under a hot and humid environment. A genome-wide association study (GWAS) allows identification of novel quantitative trait loci (QTLs) associated with SF and other partitioning traits that can assist in marker assisted breeding. In this study, GWAS was performed on a soft wheat association mapping panel (SWAMP) comprised of 236 elite lines using 27,466 single nucleotide polymorphisms (SNPs). The panel was phenotyped in two heat stress locations over 3 years. GWAS identified 109 significant marker-trait associations (MTAs) (p ≤ 9.99 x 10-5) related to eight phenotypic traits including SF (a major component of grain number) and spike harvest index (SHI, a major component of grain weight). MTAs detected on chromosomes 1B, 3A, 3B, and 5A were associated with multiple traits and are potentially important targets for selection. More than half of the significant MTAs (60 out of 109) were found in genes encoding different types of proteins related to metabolism, disease, and abiotic stress including heat stress. These MTAs could be potential targets for further validation study and may be used in marker-assisted breeding for improving wheat grain yield under post-anthesis heat stress conditions. This is the first study to identify novel QTLs associated with SF and SHI which represent the major components of grain number and grain weight, respectively, in wheat.

Multi-Trait, Multi-Environment Genomic Prediction of Durum Wheat With Genomic Best Linear Unbiased Predictor and Deep Learning Methods

Although durum wheat (Triticum turgidum var. durum Desf.) is a minor cereal crop representing just 5-7% of the world's total wheat crop, it is a staple food in Mediterranean countries, where it is used to produce pasta, couscous, bulgur and bread. In this paper, we cover multi-trait prediction of grain yield (GY), days to heading (DH) and plant height (PH) of 270 durum wheat lines that were evaluated in 43 environments (country-location-year combinations) across a broad range of water regimes in the Mediterranean Basin and other locations. Multi-trait prediction analyses were performed by implementing a multi-trait deep learning model (MTDL) with a feed-forward network topology and a rectified linear unit activation function with a grid search approach for the selection of hyper-parameters. The results of the multi-trait deep learning method were also compared with univariate predictions of the genomic best linear unbiased predictor (GBLUP) method and the univariate counterpart of the multi-trait deep learning method (UDL). All models were implemented with and without the genotype × environment interaction term. We found that the best predictions were observed without the genotype × environment interaction term in the UDL and MTDL methods. However, under the GBLUP method, the best predictions were observed when the genotype × environment interaction term was taken into account. We also found that in general the best predictions were observed under the GBLUP model; however, the predictions of the MTDL were very similar to those of the GBLUP model. This result provides more evidence that the GBLUP model is a powerful approach for genomic prediction, but also that the deep learning method is a practical approach for predicting univariate and multivariate traits in the context of genomic selection.

Genetic dissection of drought and heat-responsive agronomic traits in wheat

High yield and wide adaptation are principal targets of wheat breeding but are hindered by limited knowledge on genetic basis of agronomic traits and abiotic stress tolerances. In this study, 277 wheat accessions were phenotyped across 30 environments with non-stress, drought-stressed, heat-stressed, and drought-heat-stressed treatments and were subjected to genome-wide association study using 395 681 single nucleotide polymorphisms. We detected 295 associated loci including consistent loci for agronomic traits across different treatments and eurytopic loci for multiple abiotic stress tolerances. A total of 22 loci overlapped with quantitative trait loci identified by biparental quantitative trait loci mapping. Six loci were simultaneously associated with agronomic traits and abiotic stress tolerance, four of which fell within selective sweep regions. Selection in Chinese wheat has increased the frequency of superior marker alleles controlling yield-related traits in the four loci during past decades, which conversely diminished favourable genetic variation controlling abiotic stress tolerance in the same loci; two promising candidate paralogous genes colocalized with such loci, thereby providing potential targets for studying the molecular mechanism of stress tolerance-productivity trade-off. These results uncovering promising alleles controlling agronomic traits and/or multiple abiotic stress tolerances, providing insights into heritable covariation between yield and abiotic stress tolerance, will accelerate future efforts for wheat improvement.

A Genome-Wide Association Study of Highly Heritable Agronomic Traits in Durum Wheat

Uncovering the genetic basis of key agronomic traits, and particularly of drought tolerance, addresses an important priority for durum wheat improvement. Here, a genome-wide association study (GWAS) in 493 durum wheat accessions representing a worldwide collection was employed to address the genetic basis of 17 agronomically important traits and a drought wilting score. Using a linear mixed model with 4 inferred subpopulations and a kinship matrix, we identified 90 marker-trait-associations (MTAs) defined by 78 markers. These markers could be merged into 44 genomic loci by linkage disequilibrium (r 2 > 0.2). Based on sequence alignment of the markers to the reference genome of bread wheat, we identified 14 putative candidate genes involved in enzymes, hormone-response, and transcription factors. The GWAS in durum wheat and a previous quantitative trait locus (QTL) analysis in bread wheat identified a consensus QTL locus.4B.1 conferring drought tolerance, which was further scanned for the presence of potential candidate genes. A haplotype analysis of this region revealed that two minor haplotypes were associated with both drought tolerance and reduced plant stature, thought to be the effect of linkage with the semi-dwarfing gene Rht-B1. Haplotype variants in the key chromosome 4B region were informative regarding evolutionary divergence among durum, emmer and bread wheat. Over all, the data are relevant in the context of durum wheat improvement and the isolation of genes underlying variation in some important quantitative traits.

Genetic Dissection of the Seminal Root System Architecture in Mediterranean Durum Wheat Landraces by Genome-Wide Association Study

Roots are crucial for adaptation to drought stress. However, phenotyping root systems is a difficult and time-consuming task due to the special feature of the traits in the process of being analyzed. Correlations between root system architecture (RSA) at the early stages of development and in adult plants have been reported. In this study, the seminal RSA was analysed on a collection of 160 durum wheat landraces from 21 Mediterranean countries and 18 modern cultivars. The landraces showed large variability in RSA, and differences in root traits were found between previously identified genetic subpopulations. Landraces from the eastern Mediterranean region, which is the driest and warmest within the Mediterranean Basin, showed the largest seminal root size in terms of root length, surface, and volume and the widest root angle, whereas landraces from eastern Balkan countries showed the lowest values. Correlations were found between RSA and yield-related traits in a very dry environment. The identification of molecular markers linked to the traits of interest detected 233 marker-trait associations for 10 RSA traits and grouped them in 82 genome regions named marker-train association quantitative trait loci (MTA-QTLs). Our results support the use of ancient local germplasm to widen the genetic background for root traits in breeding programs.

A large nested association mapping population for breeding and quantitative trait locus mapping in Ethiopian durum wheat

The Ethiopian plateau hosts thousands of durum wheat (Triticum turgidum subsp. durum) farmer varieties (FV) with high adaptability and breeding potential. To harness their unique allelic diversity, we produced a large nested association mapping (NAM) population intercrossing fifty Ethiopian FVs with an international elite durum wheat variety (Asassa). The Ethiopian NAM population (EtNAM) is composed of fifty interconnected bi-parental families, totalling 6280 recombinant inbred lines (RILs) that represent both a powerful quantitative trait loci (QTL) mapping tool, and a large pre-breeding panel. Here, we discuss the molecular and phenotypic diversity of the EtNAM founder lines, then we use an array featuring 13 000 single nucleotide polymorphisms (SNPs) to characterize a subset of 1200 EtNAM RILs from 12 families. Finally, we test the usefulness of the population by mapping phenology traits and plant height using a genome wide association (GWA) approach. EtNAM RILs showed high allelic variation and a genetic makeup combining genetic diversity from Ethiopian FVs with the international durum wheat allele pool. EtNAM SNP data were projected on the fully sequenced AB genome of wild emmer wheat, and were used to estimate pairwise linkage disequilibrium (LD) measures that reported an LD decay distance of 7.4 Mb on average, and balanced founder contributions across EtNAM families. GWA analyses identified 11 genomic loci individually affecting up to 3 days in flowering time and more than 1.6 cm in height. We argue that the EtNAM is a powerful tool to support the production of new durum wheat varieties targeting local and global agriculture.

Elucidating the genetic basis of biomass accumulation and radiation use efficiency in spring wheat and its role in yield potential

One of the major challenges for plant scientists is increasing wheat (Triticum aestivum) yield potential (YP). A significant bottleneck for increasing YP is achieving increased biomass through optimization of radiation use efficiency (RUE) along the crop cycle. Exotic material such as landraces and synthetic wheat has been incorporated into breeding programmes in an attempt to alleviate this; however, their contribution to YP is still unclear. To understand the genetic basis of biomass accumulation and RUE, we applied genome-wide association study (GWAS) to a panel of 150 elite spring wheat genotypes including many landrace and synthetically derived lines. The panel was evaluated for 31 traits over 2 years under optimal growing conditions and genotyped using the 35K wheat breeders array. Marker-trait association identified 94 SNPs significantly associated with yield, agronomic and phenology-related traits along with RUE and final biomass (BM_PM) at various growth stages that explained 7%-17% of phenotypic variation. Common SNP markers were identified for grain yield, BM_PM and RUE on chromosomes 5A and 7A. Additionally, landrace and synthetic derivative lines showed higher thousand grain weight (TGW), BM_PM and RUE but lower grain number (GM2) and harvest index (HI). Our work demonstrates the use of exotic material as a valuable resource to increase YP. It also provides markers for use in marker-assisted breeding to systematically increase BM_PM, RUE and TGW and avoid the TGW/GM2 and BM_PM/HI trade-off. Thus, achieving greater genetic gains in elite germplasm while also highlighting genomic regions and candidate genes for further study.

Drought Stress Tolerance in Wheat and Barley: Advances in Physiology, Breeding and Genetics Research

Climate change is a major threat to most of the agricultural crops grown in tropical and sub-tropical areas globally. Drought stress is one of the consequences of climate change that has a negative impact on crop growth and yield. In the past, many simulation models were proposed to predict climate change and drought occurrences, and it is extremely important to improve essential crops to meet the challenges of drought stress which limits crop productivity and production. Wheat and barley are among the most common and widely used crops due to their economic and social values. Many parts of the world depend on these two crops for food and feed, and both crops are vulnerable to drought stress. Improving drought stress tolerance is a very challenging task for wheat and barley researchers and more research is needed to better understand this stress. The progress made in understanding drought tolerance is due to advances in three main research areas: physiology, breeding, and genetic research. The physiology research focused on the physiological and biochemical metabolic pathways that plants use when exposed to drought stress. New wheat and barley genotypes having a high degree of drought tolerance are produced through breeding by making crosses from promising drought-tolerant genotypes and selecting among their progeny. Also, identifying genes contributing to drought tolerance is very important. Previous studies showed that drought tolerance is a polygenic trait and genetic constitution will help to dissect the gene network(s) controlling drought tolerance. This review explores the recent advances in these three research areas to improve drought tolerance in wheat and barley.

Drought Stress Tolerance in Wheat and Barley: Advances in Physiology, Breeding and Genetics Research

Climate change is a major threat to most of the agricultural crops grown in tropical and sub-tropical areas globally. Drought stress is one of the consequences of climate change that has a negative impact on crop growth and yield. In the past, many simulation models were proposed to predict climate change and drought occurrences, and it is extremely important to improve essential crops to meet the challenges of drought stress which limits crop productivity and production. Wheat and barley are among the most common and widely used crops due to their economic and social values. Many parts of the world depend on these two crops for food and feed, and both crops are vulnerable to drought stress. Improving drought stress tolerance is a very challenging task for wheat and barley researchers and more research is needed to better understand this stress. The progress made in understanding drought tolerance is due to advances in three main research areas: physiology, breeding, and genetic research. The physiology research focused on the physiological and biochemical metabolic pathways that plants use when exposed to drought stress. New wheat and barley genotypes having a high degree of drought tolerance are produced through breeding by making crosses from promising drought-tolerant genotypes and selecting among their progeny. Also, identifying genes contributing to drought tolerance is very important. Previous studies showed that drought tolerance is a polygenic trait and genetic constitution will help to dissect the gene network(s) controlling drought tolerance. This review explores the recent advances in these three research areas to improve drought tolerance in wheat and barley.

Uncovering Genomic Regions Associated With 36 Agro-Morphological Traits in Indian Spring Wheat Using GWAS

Wheat genetic improvement by integration of advanced genomic technologies is one way of improving productivity. To facilitate the breeding of economically important traits in wheat, SNP loci and underlying candidate genes associated with the 36 agro-morphological traits were studied in a diverse panel of 404 genotypes. By using Breeders' 35K Axiom array in a comprehensive genome-wide association study covering 4364.79 cM of the wheat genome and applying a compressed mixed linear model, a total of 146 SNPs (-log10 P ≥ 4) were found associated with 23 traits out of 36 traits studied explaining 3.7-47.0% of phenotypic variance. To reveal this a subset of 260 genotypes was characterized phenotypically for six quantitative traits [days to heading (DTH), days to maturity (DTM), plant height (PH), spike length (SL), awn length (Awn_L), and leaf length (Leaf_L)] under five environments. Gene annotations mined ∼38 putative candidate genes which were confirmed using tissue and stage specific gene expression data from RNA Seq. We observed strong co-localized loci for four traits (glume pubescence, SL, PH, and awn color) on chromosome 1B (24.64 cM) annotated five putative candidate genes. This study led to the discovery of hitherto unreported loci for some less explored traits (such as leaf sheath wax, awn attitude, and glume pubescence) besides the refined chromosomal regions of known loci associated with the traits. This study provides valuable information of the genetic loci and their potential genes underlying the traits such as awn characters which are being considered as important contributors toward yield enhancement.

Genome-wide association analyses for yield and yield-related traits in bread wheat (Triticum aestivum L.) under pre-anthesis combined heat and drought stress in field conditions

Understanding the genetic basis of heat and drought stress tolerance in wheat is prerequisite for wheat breeding program. In the present study, a wheat panel comprising of 192 elite bread wheat genotypes was phenotyped in eight environments for yield and related traits in field conditions. Four stress environments were created by implying four different treatments differing in sowing date and water availability, panel was evaluated for two years in field conditions. The panel was genotyped with 15K Illumina chip and 9236 polymorphic markers concentrated on B genome were employed in GWAS analysis. Consistent, fast LD decay was observed on D genome and structure analysis germplasm divided panel into three major populations. GWAS was performed using BLUEs values of combined environment data in R package GAPIT using log10(P) = 3.96 as significance threshold. The significance of association was further checked using FDR<0.05 threshold. The GWAS identified 487 loci associated with the traits and were significant at log10(p) threshold out of these 350 loci were significant at FDR threshold. For two stress indices 108 associations were significant at FDR threshold. Nine genomic regions were shared among all treatment, while multiple pleiotropic regions were present on chromosome 7D followed by unmapped chromosome. The present study validated many marker trait associations for yield and other traits, MTAs significant under combined drought and heat stress were novel. These regions are important and can be used for fine mapping and marker assisted selection to discover new genes responsible for heat and drought tolerance in wheat.

Genetic dissection of agronomic and quality traits based on association mapping and genomic selection approaches in durum wheat grown in Southern Spain

Climatic conditions affect the growth, development and final crop production. As wheat is of paramount importance as a staple crop in the human diet, there is a growing need to study its abiotic stress adaptation through the performance of key breeding traits. New and complementary approaches, such as genome-wide association studies (GWAS) and genomic selection (GS), are used for the dissection of different agronomic traits. The present study focused on the dissection of agronomic and quality traits of interest (initial agronomic score, yield, gluten index, sedimentation index, specific weight, whole grain protein and yellow colour) assessed in a panel of 179 durum wheat lines (Triticum durum Desf.), grown under rainfed conditions in different Mediterranean environments in Southern Spain (Andalusia). The findings show a total of 37 marker-trait associations (MTAs) which affect phenotype expression for three quality traits (specific weight, gluten and sedimentation indexes). MTAs could be mapped on the A and B durum wheat subgenomes (on chromosomes 1A, 1B, 2A, 2B and 3A) through the recently available bread wheat reference assembly (IWGSC RefSeqv1). Two of the MTAs found for quality traits (gluten index and SDS) corresponded to the known Glu-B1 and Glu-A1 loci, for which candidate genes corresponding to high molecular weight glutenin subunits could be located. The GS prediction ability values obtained from the breeding materials analyzed showed promising results for traits as grain protein content, sedimentation and gluten indexes, which can be used in plant breeding programs.

Genome-wide association studies for yield-related traits in soft red winter wheat grown in Virginia

Grain yield is a trait of paramount importance in the breeding of all cereals. In wheat (Triticum aestivum L.), yield has steadily increased since the Green Revolution, though the current rate of increase is not forecasted to keep pace with demand due to growing world population and increasing affluence. While several genome-wide association studies (GWAS) on yield and related component traits have been performed in wheat, the previous lack of a reference genome has made comparisons between studies difficult. In this study, a GWAS for yield and yield-related traits was carried out on a population of 322 soft red winter wheat lines across a total of four rain-fed environments in the state of Virginia using single-nucleotide polymorphism (SNP) marker data generated by a genotyping-by-sequencing (GBS) protocol. Two separate mixed linear models were used to identify significant marker-trait associations (MTAs). The first was a single-locus model utilizing a leave-one-chromosome-out approach to estimating kinship. The second was a sub-setting kinship estimation multi-locus method (FarmCPU). The single-locus model identified nine significant MTAs for various yield-related traits, while the FarmCPU model identified 74 significant MTAs. The availability of the wheat reference genome allowed for the description of MTAs in terms of both genetic and physical positions, and enabled more extensive post-GWAS characterization of significant MTAs. The results indicate a number of promising candidate genes contributing to grain yield, including an ortholog of the rice aberrant panicle organization (APO1) protein and a gibberellin oxidase protein (GA2ox-A1) affecting the trait grains per square meter, an ortholog of the Arabidopsis thaliana mother of flowering time and terminal flowering 1 (MFT) gene affecting the trait seeds per square meter, and a B2 heat stress response protein affecting the trait seeds per head.

Novel Alleles for Combined Drought and Heat Stress Tolerance in Wheat

Drought and heat waves commonly co-occur in many wheat-growing regions causing significant crop losses. The identification of stress associated quantitative trait loci, particularly those for yield, is problematic due to their association with plant phenology and the high genetic × environment interaction. Here we studied a panel of 315 diverse, spring type accessions of bread wheat (Triticum aestivum) in pots in a semi-controlled environment under combined drought and heat stress over 2 years. Importantly, we treated individual plants according to their flowering time. We found 134 out of the 145 identified loci for grain weight that were not associated with either plant phenology or plant height. The majority of loci uncovered here were novel, with favorable alleles widespread in Asian and African landraces providing opportunities for their incorporation into modern varieties through breeding. Using residual heterozygosity in lines from a nested association mapping population, we were able to rapidly develop near-isogenic lines for important target loci. One target locus on chromosome 6A contributed to higher grain weight, harvest index, thousand kernel weight, and grain number under drought and heat stress in field conditions consistent with allelic effects demonstrated in the genome-wide association study.

Identification of Two Novel Wheat Drought Tolerance-Related Proteins by Comparative Proteomic Analysis Combined with Virus-Induced Gene Silencing

Drought is a major adversity that limits crop yields. Further exploration of wheat drought tolerance-related genes is critical for the genetic improvement of drought tolerance in this crop. Here, comparative proteomic analysis of two wheat varieties, XN979 and LA379, with contrasting drought tolerance was conducted to screen for drought tolerance-related proteins/genes. Virus-induced gene silencing (VIGS) technology was used to verify the functions of candidate proteins. A total of 335 differentially abundant proteins (DAPs) were exclusively identified in the drought-tolerant variety XN979. Most DAPs were mainly involved in photosynthesis, carbon fixation, glyoxylate and dicarboxylate metabolism, and several other pathways. Two DAPs (W5DYH0 and W5ERN8), dubbed TaDrSR1 and TaDrSR2, respectively, were selected for further functional analysis using VIGS. The relative electrolyte leakage rate and malonaldehyde content increased significantly, while the relative water content and proline content significantly decreased in the TaDrSR1- and TaDrSR2-knock-down plants compared to that in non-knocked-down plants under drought stress conditions. TaDrSR1- and TaDrSR2-knock-down plants exhibited more severe drooping and wilting phenotypes than non-knocked-down plants under drought stress conditions, suggesting that the former were more sensitive to drought stress. These results indicate that TaDrSR1 and TaDrSR2 potentially play vital roles in conferring drought tolerance in common wheat.

Haplotype Loci Under Selection in Canadian Durum Wheat Germplasm Over 60 Years of Breeding: Association With Grain Yield, Quality Traits, Protein Loss, and Plant Height

Durum wheat was introduced in the southern prairies of western Canada in the late nineteenth century. Breeding efforts have mainly focused on improving quality traits to meet the pasta industry demands. For this study, 192 durum wheat lines were genotyped using the Illumina 90K Infinium iSelect assay, and resulted in a total of 14,324 polymorphic SNPs. Genetic diversity changed over time, declining during the first 20 years of breeding in Canada, then increased in the late 1980s and early 1990s. We scanned the genome for signatures of selection, using the total variance Fst-based outlier detection method (Lositan), the hierarchical island model (Arlequin) and the Bayesian genome scan method (BayeScan). A total of 407 outliers were identified and clustered into 84 LD-based haplotype loci, spanning all 14 chromosomes of the durum wheat genome. The association analysis detected 54 haplotype loci, of which 39% contained markers with a complete reversal of allelic state. This tendency to fixation of favorable alleles corroborates the success of the Canadian durum wheat breeding programs over time. Twenty-one haplotype loci were associated with multiple traits. In particular, hap_4B_1 explained 20.6, 17.9 and 16.6% of the phenotypic variance of pigment loss, pasta b∗ and dough extensibility, respectively. The locus hap_2B_9 explained 15.9 and 17.8% of the variation of protein content and protein loss, respectively. All these pleiotropic haplotype loci offer breeders the unique opportunity for further improving multiple traits, facilitating marker-assisted selection in durum wheat, and could help in identifying genes as functional annotations of the wheat genome become available.

Genomic-enabled Prediction Accuracies Increased by Modeling Genotype × Environment Interaction in Durum Wheat

Genomic prediction studies incorporating genotype × environment (G×E) interaction effects are limited in durum wheat. We tested the genomic-enabled prediction accuracy (PA) of Genomic Best Linear Unbiased Predictor (GBLUP) models-six non-G × E and three G × E models-on three basic cross-validation (CV) schemes- in predicting incomplete field trials (CV2), new lines (CV1), and lines in untested environments (CV0)- in a durum wheat panel grown under yield potential, drought stress, and heat stress conditions. For CV0, three scenarios were considered: (i) leave-one environment out (CV0-Env); (ii) leave one site out (CV0-Site); and (iii) leave 1 yr out (CV0-Year). The reaction norm models with G × E effects showed higher PA than the non-G × E models. Among the CV schemes, CV2 and CV0-Env had higher PA (0.58 each) than the CV1 scheme (0.35). When the average of all the models and CV schemes were considered, among the eight traits- grain yield, thousand grain weight, grain number, days to anthesis, days to maturity, plant height, and normalized difference vegetation index at vegetative (NDVIvg) and grain filling (NDVIllg)-, plant height had the highest PA (0.68) and moderate values were observed for grain yield (0.34). The results indicated that genomic selection models incorporating G × E interaction show great promise for forward prediction and application in durum wheat breeding to increase genetic gains.