Genetic Dissection of Three Major Quantitative Trait Loci for Spike Compactness and Length in Bread Wheat (Triticum aestivum L.)

Discovering novel genomic regions explaining adaptation of bread wheat to conservation agriculture through GWAS

To sustainably increase wheat yield to meet the growing world population's food demand in the face of climate change, Conservation Agriculture (CA) is a promising approach. Still, there is a lack of genomic studies investigating the genetic basis of crop adaptation to CA. To dissect the genetic architecture of 19 morpho-physiological traits that could be involved in the enhanced adaptation and performance of genotypes under CA, we performed GWAS to identify MTAs under four contrasting production regimes viz., conventional tillage timely sown (CTTS), conservation agriculture timely sown (CATS), conventional tillage late sown (CTLS) and conservation agriculture late sown (CALS) using an association panel of 183 advanced wheat breeding lines along with 5 checks. Traits like Phi2 (Quantum yield of photosystem II; CATS:0.37, CALS: 0.31), RC (Relative chlorophyll content; CATS:55.51, CALS: 54.47) and PS1 (Active photosystem I centers; CATS:2.45, CALS: 2.23) have higher mean values in CA compared to CT under both sowing times. GWAS identified 80 MTAs for the studied traits across four production environments. The phenotypic variation explained (PVE) by these QTNs ranged from 2.15 to 40.22%. Gene annotation provided highly informative SNPs associated with Phi2, NPQ (Quantum yield of non-photochemical quenching), PS1, and RC which were linked with genes that play crucial roles in the physiological adaptation under both CA and CT. A highly significant SNP AX94651261 (9.43% PVE) was identified to be associated with Phi2, while two SNP markers AX94730536 (30.90% PVE) and AX94683305 (16.99% PVE) were associated with NPQ. Identified QTNs upon validation can be used in marker-assisted breeding programs to develop CA adaptive genotypes.

Identification, validation and candidate gene analysis of major QTL for Supernumerary spikelets in wheat

BACKGROUND

The number of spikelets per spike is a key trait that affects the yield of bread wheat (Triticum aestivum L.). Identification of the QTL for spikelets per spike and its genetic effects that could be used in molecular assistant breeding in the future.

RESULTS

In this study, four recombinant inbred line (RIL) populations were generated and used, having YuPi branching wheat (YP), with Supernumerary Spikelets (SS) phenotype, as a common parent. QTL (QSS.sicau-2 A and QSS.sicau-2D) related to SS trait were mapped on chromosomes 2 A and 2D through bulked segregant exome sequencing (BSE-Seq). Fourteen molecular markers were further developed within the localization interval, and QSS.sicau-2 A was narrowed to 3.0 cM covering 7.6 Mb physical region of the reference genome, explaining 13.7 - 15.9% the phenotypic variance. Similarly, the QSS.sicau-2D was narrowed to 1.8 cM covering 2.4 Mb physical region of the reference genome, and it explained 27.4 - 32.9% the phenotypic variance. These two QTL were validated in three different genetic backgrounds using the linked markers. QSS.sicau-2 A was identified as WFZP-A, and QSS.sicau-2D was identified a novel locus, different to the previously identified WFZP-D. Based on the gene expression patterns, gene annotation and sequence analysis, TraesCS2D03G0260700 was predicted to be a potential candidate gene for QSS.sicau-2D.

CONCLUSION

Two significant QTL for SS, namely QSS.sicau-2 A and QSS.sicau-2D were identified in multiple environments were identified and their effect in diverse genetic populations was assessed. QSS.sicau-2D is a novel QTL associated with the SS trait, with TraesCS2D03G0260700 predicted as its candidate gene.

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.

Dissection of a novel major stable QTL on chromosome 7D for grain hardness and its breeding value estimation in bread wheat

Grain hardness (Gh) is important for wheat processing and end-product quality. Puroindolines polymorphism explains over 60% of Gh variation and the novel genetic factors remain to be exploited. In this study, a total of 153 quantitative trait loci (QTLs), clustered into 12 genomic intervals (C1-C12), for 13 quality-related traits were identified using a recombinant inbred line population derived from the cross of Zhongkemai138 (ZKM138) and Chuanmai44 (CM44). Among them, C7 (harboring eight QTLs for different quality-related traits) and C8 (mainly harboring QGh.cib-5D.1 for Gh) were attributed to the famous genes, Rht-D1 and Pina, respectively, indicating that the correlation of involved traits was supported by the pleotropic or linked genes. Notably, a novel major stable QTL for Gh was detected in the C12, QGh.cib-7D, with ZKM138-derived allele increasing grain hardness, which was simultaneously mapped by the BSE-Seq method. The geographic pattern and transmissibility of this locus revealed that the increasing-Gh allele is highly frequently present in 85.79% of 373 worldwide wheat varieties and presented 99.31% transmissibility in 144 ZKM138-derivatives, indicating the non-negative effect on yield performance and that its indirect passive selection has happened during the actual breeding process. Thus, the contribution of this new Gh-related locus was highlighted in consideration of improving the efficiency and accuracy of the soft/hard material selection in the molecular marker-assisted process. Further, TraesCS7D02G099400, TraesCS7D02G098000, and TraesCS7D02G099500 were initially deduced to be the most potential candidate genes of QGh.cib-7D. Collectively, this study provided valuable information of elucidating the genetic architecture of Gh for wheat quality improvement.

Genetic dissection of major QTL for grain number per spike on chromosomes 5A and 6A in bread wheat (Triticum aestivum L.)

Grain number per spike (GNS) is a crucial component of grain yield and plays a significant role in improving wheat yield. To identify quantitative trait loci (QTL) associated with GNS, a recombinant inbred line (RIL) population derived from the cross of Zhongkemai 13F10 and Chuanmai 42 was employed to conduct QTL mapping across eight environments. Based on the bulked segregant exome sequencing (BSE-Seq), genomic regions associated with GNS were detected on chromosomes 5A and 6A. According to the constructed genetic maps, two major QTL QGns.cib-5A (LOD = 4.35-8.16, PVE = 8.46-14.43%) and QGns.cib-6A (LOD = 3.82-30.80, PVE = 5.44-12.38%) were detected in five and four environments, respectively. QGns.cib-6A is a QTL cluster for other seven yield-related traits. QGns.cib-5A and QGns.cib-6A were further validated using linked Kompetitive Allele Specific PCR (KASP) markers in different genetic backgrounds. QGns.cib-5A exhibited pleiotropic effects on productive tiller number (PTN), spike length (SL), fertile spikelet number per spike (FSN), and ratio of grain length to grain width (GL/GW) but did not significantly affect thousand grain weight (TGW). Haplotype analysis revealed that QGns.cib-5A and QGns.cib-6A were the targets of artificial selection during wheat improvement. Candidate genes for QGns.cib-5A and QGns.cib-6A were predicted by analyzing gene annotation, spatiotemporal expression patterns, and orthologous and sequence differences. These findings will be valuable for fine mapping and map-based cloning of genes underlying QGns.cib-5A and QGns.cib-6A.

Identification of a novel SNP in the miR172 binding site of Q homoeolog AP2L-D5 is associated with spike compactness and agronomic traits in wheat (Triticum aestivum L.)

KEY MESSAGE

This study found that the compact spike locus of ANK-15 is on chromosome 5D instead of 2B. We have identified a new allele of AP2L-D5 as the candidate causal polymorphism. Spike architecture is a key determinant of wheat yield, a crop which supports much of the human diet but whose yield gains are stagnating. Spike architecture mutants offer opportunities to identify genetic factors contributing to inflorescence development. Here, we investigate the locus underlying the compact spike phenotype of mutant line ANK-15 by conducting mRNA-sequencing and genetic mapping using ANK-15 and its non-compact spike near-isogenic line Novosibirskaya 67 (N67). Previous literature has placed the compact spike locus of ANK-15 to chromosome 2B. However, based on the single nucleotide polymorphisms (SNPs) identified using mRNA-seq data, we were unable to detect polymorphisms between N67 and ANK-15 in the putative chromosome 2B region. We performed differential expression analysis of developing rachis and found that AP2L-D5, the D homoeolog of the domestication Q gene, is upregulated in ANK-15 in comparison to N67. ANK-15 carries a SNP in the microRNA172 binding site of AP2L-D5, which is predicted to lead to higher expression of AP2L-D5 due to decreased miRNA172-mediated degradation. Furthermore, we performed genetic mapping using an ANK-15 × N67 F2 population and found a single quantitative trait locus on chromosome 5D coinciding with the position of AP2L-D5. This result suggests that AP2L-D5 is likely the underlying causal gene for the compact spike phenotype in ANK-15. We performed a field trial to investigate the effect of the AP2L-D5 allele on agronomic traits and found that the AP2L-D5 allele from ANK-15 is associated with a significant reduction in height, increased thousand grain weight (TGW), and increased grain width.

Physical mapping of the wheat genes in low-recombination regions: radiation hybrid mapping of the C-locus

KEY MESSAGE

This work reports the physical mapping of an important gene affecting spike compactness located in a low-recombination region of hexaploid wheat. This work paves the way for the eventual isolation and characterization of the factor involved but also opens up possibilities to use this approach to precisely map other wheat genes located on proximal parts of wheat chromosomes that show highly reduced recombination. Mapping wheat genes, in the centromeric and pericentromeric regions (~ 2/3rd of a given chromosome), poses a formidable challenge due to highly suppressed recombination. Using an example of compact spike locus (C-locus), this study provides an approach to precisely map wheat genes in the pericentromeric and centromeric regions that house ~ 30% of wheat genes. In club-wheat, spike compactness is controlled by the dominant C-locus, but previous efforts have failed to localize it, on a particular arm of chromosome 2D. We integrated radiation hybrid (RH) and high-resolution genetic mapping to locate C-locus on the short arm of chromosome 2D. Flanking markers of the C-locus span a physical distance of 11.0 Mb (231.0-242 Mb interval) and contain only 11 high-confidence annotated genes. This work demonstrates the value of this integrated strategy in mapping dominant genes in the low-recombination regions of the wheat genome. A comparison of the mapping resolutions of the RH and genetic maps using common anchored markers indicated that the RH map provides ~ 9 times better resolution that the genetic map even with much smaller population size. This study provides a broadly applicable approach to fine map wheat genes in regions of suppressed recombination.