Multi-locus genome-wide association mapping for spike-related traits in bread wheat (Triticum aestivum L.)

Identification of molecular markers and candidate regions associated with grain number per spike in Pubing3228 using SLAF-BSA

Grain number per spike, a pivotal agronomic trait dictating wheat yield, lacks a comprehensive understanding of its underlying mechanism in Pubing3228, despite the identification of certain pertinent genes. Thus, our investigation sought to ascertain molecular markers and candidate regions associated with grain number per spike through a high-density genetic mapping approach that amalgamates site-specific amplified fragment sequencing (SLAF-seq) and bulked segregation analysis (BSA). To facilitate this, we conducted a comparative analysis of two wheat germplasms, Pubing3228 and Jing4839, known to exhibit marked discrepancies in spike shape. By leveraging this methodology, we successfully procured 2,810,474 SLAF tags, subsequently resulting in the identification of 187,489 single nucleotide polymorphisms (SNPs) between the parental strains. We subsequently employed the SNP-index association algorithm alongside the extended distribution (ED) association algorithm to detect regions associated with the trait. The former algorithm identified 24 trait-associated regions, whereas the latter yielded 70. Remarkably, the intersection of these two algorithms led to the identification of 25 trait-associated regions. Amongst these regions, we identified 399 annotated genes, including three genes harboring non-synonymous mutant SNP loci. Notably, the APETALA2 (AP2) transcription factor families, which exhibited a strong correlation with spike type, were also annotated. Given these findings, it is plausible to hypothesize that these genes play a critical role in determining spike shape. In summation, our study contributes significant insights into the genetic foundation of grain number per spike. The molecular markers and candidate regions we have identified can be readily employed for marker-assisted breeding endeavors, ultimately leading to the development of novel wheat cultivars possessing enhanced yield potential. Furthermore, conducting further functional analyses on the identified genes will undoubtedly facilitate a comprehensive elucidation of the underlying mechanisms governing spike development in wheat.

Association mapping reveals novel genes and genomic regions controlling grain size architecture in mini core accessions of Indian National Genebank wheat germplasm collection

Wheat (Triticum aestivum L.) is a staple food crop for the global human population, and thus wheat breeders are consistently working to enhance its yield worldwide. In this study, we utilized a sub-set of Indian wheat mini core germplasm to underpin the genetic architecture for seed shape-associated traits. The wheat mini core subset (125 accessions) was genotyped using 35K SNP array and evaluated for grain shape traits such as grain length (GL), grain width (GW), grain length, width ratio (GLWR), and thousand grain weight (TGW) across the seven different environments (E₁, E₂, E₃, E₄, E₅, E₅, E₆, and E₇). Marker-trait associations were determined using a multi-locus random-SNP-effect Mixed Linear Model (mrMLM) program. A total of 160 non-redundant quantitative trait nucleotides (QTNs) were identified for four grain shape traits using two or more GWAS models. Among these 160 QTNs, 27, 36, 38, and 35 QTNs were associated for GL, GW, GLWR, and TGW respectively while 24 QTNs were associated with more than one trait. Of these 160 QTNs, 73 were detected in two or more environments and were considered reliable QTLs for the respective traits. A total of 135 associated QTNs were annotated and located within the genes, including ABC transporter, Cytochrome450, Thioredoxin_M-type, and hypothetical proteins. Furthermore, the expression pattern of annotated QTNs demonstrated that only 122 were differentially expressed, suggesting these could potentially be related to seed development. The genomic regions/candidate genes for grain size traits identified in the present study represent valuable genomic resources that can potentially be utilized in the markers-assisted breeding programs to develop high-yielding varieties.

Identification of novel putative alleles related to important agronomic traits of wheat using robust strategies in GWAS

Principal component analysis (PCA) is widely used in various genetics studies. In this study, the role of classical PCA (cPCA) and robust PCA (rPCA) was evaluated explicitly in genome-wide association studies (GWAS). We evaluated 294 wheat genotypes under well-watered and rain-fed, focusing on spike traits. First, we showed that some phenotypic and genotypic observations could be outliers based on cPCA and different rPCA algorithms (Proj, Grid, Hubert, and Locantore). Hubert's method provided a better approach to identifying outliers, which helped to understand the nature of these samples. These outliers led to the deviation of the heritability of traits from the actual value. Then, we performed GWAS with 36,000 single nucleotide polymorphisms (SNPs) based on the traditional approach and two robust strategies. In the conventional approach and using the first three components of cPCA as population structure, 184 and 139 marker-trait associations (MTAs) were identified for five traits in well-watered and rain-fed environments, respectively. In the first robust strategy and when rPCA was used as population structure in GWAS, we observed that the Hubert and Grid methods identified new MTAs, especially for yield and spike weight on chromosomes 7A and 6B. In the second strategy, we followed the classical and robust principal component-based GWAS, where the first two PCs obtained from phenotypic variables were used instead of traits. In the recent strategy, despite the similarity between the methods, some new MTAs were identified that can be considered pleiotropic. Hubert's method provided a better linear combination of traits because it had the most MTAs in common with the traditional approach. Newly identified SNPs, including rs19833 (5B) and rs48316 (2B), were annotated with important genes with vital biological processes and molecular functions. The approaches presented in this study can reduce the misleading GWAS results caused by the adverse effect of outlier observations.

Genetics of spot blotch resistance in bread wheat (Triticum aestivum L.) using five models for GWAS

Genetic architecture of resistance to spot blotch in wheat was examined using a Genome-Wide Association Study (GWAS) involving an association panel comprising 303 diverse genotypes. The association panel was evaluated at two different locations in India including Banaras Hindu University (BHU), Varanasi (Uttar Pradesh), and Borlaug Institute for South Asia (BISA), Pusa, Samastipur (Bihar) for two consecutive years (2017-2018 and 2018-2019), thus making four environments (E1, BHU 2017-18; E2, BHU 2018-19; E3, PUSA, 2017-18; E4, PUSA, 2018-19). The panel was genotyped for 12,196 SNPs based on DArT-seq (outsourced to DArT Ltd by CIMMYT); these SNPs included 5,400 SNPs, which could not be assigned to individual chromosomes and were therefore, described as unassigned by the vendor. Phenotypic data was recorded on the following three disease-related traits: (i) Area Under Disease Progress Curve (AUDPC), (ii) Incubation Period (IP), and (iii) Lesion Number (LN). GWAS was conducted using each of five different models, which included two single-locus models (CMLM and SUPER) and three multi-locus models (MLMM, FarmCPU, and BLINK). This exercise gave 306 MTAs, but only 89 MTAs (33 for AUDPC, 30 for IP and 26 for LN) including a solitary MTA detected using all the five models and 88 identified using four of the five models (barring SUPER) were considered to be important. These were used for further analysis, which included identification of candidate genes (CGs) and their annotation. A majority of these MTAs were novel. Only 70 of the 89 MTAs were assigned to individual chromosomes; the remaining 19 MTAs belonged to unassigned SNPs, for which chromosomes were not known. Seven MTAs were selected on the basis of minimum P value, number of models, number of environments and location on chromosomes with respect to QTLs reported earlier. These 7 MTAs, which included five main effect MTAs and two for epistatic interactions, were considered to be important for marker-assisted selection (MAS). The present study thus improved our understanding of the genetics of resistance against spot blotch in wheat and provided seven MTAs, which may be used for MAS after due validation.

Physiological implications of SWEETs in plants and their potential applications in improving source-sink relationships for enhanced yield

The sugars will eventually be exported transporters (SWEET) family of transporters in plants is identified as a novel class of sugar carriers capable of transporting sugars, sugar alcohols and hormones. Functioning in intercellular sugar transport, SWEETs influence a wide range of physiologically important processes. SWEETs regulate the development of sink organs by providing nutritional support from source leaves, responses to abiotic stresses by maintaining intracellular sugar concentrations, and host-pathogen interactions through the modulation of apoplastic sugar levels. Many bacterial and fungal pathogens activate the expression of SWEET genes in species such as rice and Arabidopsis to gain access to the nutrients that support virulence. The genetic manipulation of SWEETs has led to the generation of bacterial blight (BB)-resistant rice varieties. Similarly, while the overexpression of the SWEETs involved in sucrose export from leaves and pathogenesis led to growth retardation and yield penalties, plants overexpressing SWEETs show improved disease resistance. Such findings demonstrate the complex functions of SWEETs in growth and stress tolerance. Here, we review the importance of SWEETs in plant-pathogen and source-sink interactions and abiotic stress resistance. We highlight the possible applications of SWEETs in crop improvement programmes aimed at improving sink and source strengths important for enhancing the sustainability of yield. We discuss how the adverse effects of the overexpression of SWEETs on plant growth may be overcome.

Genetics of the Inverse Relationship between Grain Yield and Grain Protein Content in Common Wheat

Grain protein content (GPC) is one of the most important criteria to determine the quality of common wheat (Triticum aestivum). One of the major obstacles for bread wheat production is the negative correlation between GPC and grain yield (GY). Previous studies demonstrated that the deviation from this inverse relationship is highly heritable. However, little is known about the genetics controlling these deviations in common wheat. To fill this gap, we performed quantitative trait locus (QTL) analysis for GY, GPC, and four derived GY-GPC indices using an eight-way multiparent advanced generation intercross population comprising 394 lines. Interval mapping was conducted using phenotypic data from up to nine environments and genotypic data from a 20k single-nucleotide polymorphism array. The four indices were highly heritable (0.76–0.88) and showed distinct correlations to GY and GPC. Interval mapping revealed that GY, GPC, and GY-GPC indices were controlled by 6, 12, and 12 unique QTL, of which each explained only a small amount of phenotypic variance (R² ≤ 10%). Ten of the 12 index QTL were independent of loci affecting GY and GPC. QTL regions harboured several candidate genes, including Rht-1, WAPO-A1, TaTEF-7A, and NRT2.6-7A. The study confirmed the usefulness of indices to mitigate the inverse GY-GPC relationship in breeding, though the selection method should reflect their polygenic inheritance.

Identification and Validation of Quantitative Trait Loci Mapping for Spike-Layer Uniformity in Wheat

Spike-layer uniformity (SLU), the consistency of the spike distribution in the vertical space, is an important trait. It directly affects the yield potential and appearance. Revealing the genetic basis of SLU will provide new insights into wheat improvement. To map the SLU-related quantitative trait loci (QTL), 300 recombinant inbred lines (RILs) that were derived from a cross between H461 and Chinese Spring were used in this study. The RILs and parents were tested in fields from two continuous years from two different pilots. Phenotypic analysis showed that H461 was more consistent in the vertical spatial distribution of the spike layer than in Chinese Spring. Based on inclusive composite interval mapping, four QTL were identified for SLU. There were two major QTL on chromosomes 2BL and 2DL and two minor QTL on chromosomes 1BS and 2BL that were identified. The additive effects of QSlu.sicau-1B, Qslu.sicau-2B-2, and QSlu.sicau-2D were all from the parent, H461. The major QTL, QSlu.sicau-2B-2 and QSlu.sicau-2D, were detected in each of the conducted trials. Based on the best linear unbiased prediction values, the two loci explained 23.97% and 15.98% of the phenotypic variation, respectively. Compared with previous studies, the two major loci were potentially novel and the two minor loci were overlapped. Based on the kompetitive allele-specific PCR (KASP) marker, the genetic effects for QSlu.sicau-2B-2 were validated in an additional RIL population. The genetic effects ranged from 26.65% to 32.56%, with an average value of 30.40%. In addition, QSlu.sicau-2B-2 showed a significant (p < 0.01) and positive influence on the spike length, spikelet number, and thousand kernel weight. The identified QTL and the developed KASP marker will be helpful for fine-mapping these loci, finally contributing to wheat breeding programs in a marker-assisted selection way.

Association analysis for agronomic traits in wheat under terminal heat stress

•

Terminal heat stress leads to irreversible damage in wheat.

•

Marker assisted selection and gene pyramiding for portrayal of heat tolerance.

•

Allelic frequency and polymorphic information showed significant variability.

•

Markers xcfa2147 and xwmc671 could be potentail for heat stress tolerance.

Terminal heat stress causes irreversible damage to wheat crop productivity. It reduces the vegetative growth and flowering period that consequently declines the efficiency to capture available stem reserves (carbohydrates) in grains. Markers associated with thermotolerant traits ease in marker assisted selection (MAS) for crop improvement. It identifies the genomic regions associated with thermotolerant traits in wheat, but the scarcity of markers is the major hindrance in crop improvement. Therefore, 158 wheat genotypes were subjected to genotyping with 165 simple sequence repeat markers dispersed on three genomes (A, B and D). Allelic frequency and polymorphic information content values were highest on genome A (5.34 (14% greater than the lowest value at genome D) and 0.715 (3% greater than the lowest value at genome D)), chromosome 4 (5.40 (16% greater than the lowest value at chromosome 2) and 0.725 (5% greater than the lowest value at chromosome 6)) and marker xgwm44 (13.0 (84% greater than the lowest value at marker xbarc148) and 0.916 (46% greater than the lowest value at marker xbarc148)). Bayesian based population structure discriminated the wheat genotypes into seven groups based on genetic similarity indicating their ancestral origin and geographical ecotype. Linkage disequilibrium pattern had highest significant (P < 0.001) linked loci pairs 732 on genome A at r² > 0.1 whereas, 58 on genome B at r² > 0.5. Linkage disequilibrium decay (P < 0.01 and r² > 0.1) had larger LD block (5–10 cM) on genome A. Highly significant MTAs (P < 0.000061) under heat stress conditions were identified for flag leaf area (xwmc336), spikelet per spike (xwmc553), grains per spike (cxfa2147, xwmc418 and xwmc121), biomass (xbarc7) and grain yield (xcfa2147 and xwmc671). The identified markers in this study could facilitate in MAS and gene pyramiding against heat stress in wheat.

Multi-Locus GWAS for Grain Weight-Related Traits Under Rain-Fed Conditions in Common Wheat (Triticum aestivum L.)

In wheat, a multi-locus genome-wide association study (ML-GWAS) was conducted for the four grain weight-related traits (days to anthesis, grain filling duration, grain number per ear, and grain weight per ear) using data recorded under irrigated (IR) and rain-fed (RF) conditions. Seven stress-related indices were estimated for these four traits: (i) drought resistance index (DI), (ii) geometric mean productivity (GMP), (iii) mean productivity index (MPI), (iv) relative drought index (RDI), (v) stress tolerance index (STI), (vi) yield index, and (vii) yield stability index (YSI). The association panel consisted of a core collection of 320 spring wheat accessions representing 28 countries. The panel was genotyped using 9,627 single nucleotide polymorphisms (SNPs). The genome-wide association (GWA) analysis provided 30 significant marker-trait associations (MTAs), distributed as follows: (i) IR (15 MTAs), (ii) RF (14 MTAs), and (iii) IR+RF (1 MTA). In addition, 153 MTAs were available for the seven stress-related indices. Five MTAs co-localized with previously reported QTLs/MTAs. Candidate genes (CGs) associated with different MTAs were also worked out. Gene ontology (GO) analysis and expression analysis together allowed the selection of the two CGs, which may be involved in response to drought stress. These two CGs included: TraesCS1A02G331000 encoding RNA helicase and TraesCS4B02G051200 encoding microtubule-associated protein 65. The results supplemented the current knowledge on genetics for drought tolerance in wheat. The results may also be used for future wheat breeding programs to develop drought-tolerant wheat cultivars.