Mapping and characterization of major QTL for spike traits in common wheat

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.

Characterization of a wheat stable QTL for spike length and its genetic effects on yield-related traits

Spike length (SL) is one of the most important agronomic traits affecting yield potential and stability in wheat. In this study, a major stable quantitative trait locus (QTL) for SL, i.e., qSl-2B, was detected in multiple environments in a recombinant inbred line (RIL) mapping population, KJ-RILs, derived from a cross between Kenong 9204 (KN9204) and Jing 411 (J411). The qSl-2B QTL was mapped to the 60.06-73.06 Mb region on chromosome 2B and could be identified in multiple mapping populations. An InDel molecular marker in the target region was developed based on a sequence analysis of the two parents. To further clarify the breeding use potential of qSl-2B, we analyzed its genetic effects and breeding selection effect using both the KJ-RIL population and a natural mapping population, which consisted of 316 breeding varieties/advanced lines. The results showed that the qSl-2B alleles from KN9204 showed inconsistent genetic effects on SL in the two mapping populations. Moreover, in the KJ-RILs population, the additive effects analysis of qSl-2B showed that additive effect was higher when both qSl-2D and qSl-5A harbor negative alleles under LN and HN. In China, a moderate selection utilization rate for qSl-2B was found in the Huanghuai winter wheat area and the selective utilization rate for qSl-2B continues to increase. The above findings provided a foundation for the genetic improvement of wheat SL in the future via molecular breeding strategies.

Identification and Validation of a Stable Major-Effect Quantitative Trait Locus for Kernel Number per Spike on Chromosome 2D in Wheat (Triticum aestivum L.)

A recombinant inbred line population including 371 lines was developed by a high kernel number per spike (KNPS) genotype T1208 and a low KNPS genotype Chuannong18 (CN18). A genetic linkage map consisting of 11,583 markers was constructed by the Wheat55K SNP Array. The quantitative trait loci (QTLs) related to KNPS were detected in three years. Eight, twenty-seven, and four QTLs were identified using the ICIM-BIP, ICIM-MET, and ICIM-EPI methods, respectively. One QTL, QKnps.sau-2D.1, which was mapped on chromosome 2D, can explain 18.10% of the phenotypic variation (PVE) on average and be considered a major and stable QTL for KNPS. This QTL was located in a 0.89 Mb interval on chromosome 2D and flanked by the markers AX-109283238 and AX-111606890. Moreover, KASP-AX-111462389, a Kompetitive Allele-Specific PCR (KASP) marker which closely linked to QKnps.sau-2D.1, was designed. The genetic effect of QKnps.sau-2D.1 on KNPS was successfully confirmed in two RIL populations. The results also showed that the significant increase of KNPS and 1000-kernel weight (TKW) was caused by QKnps.sau-2D.1 overcoming the disadvantage due to the decrease of spike number (SN) and finally lead to a significant increase of grain yield. In addition, within the interval in which QKnps.sau-2D.1 is located in Chinese Spring reference genomes, only fifteen genes were found, and two genes that might associate with KNPS were identified. QKnps.sau-2D.1 may provide a new resource for the high-yield breeding of wheat in the future.

Identification and validation of new quantitative trait loci for spike-related traits in two RIL populations

Wheat (Triticum aestivum L.) is one of the most important cereal crops for ensuring food security worldwide. Identification of major quantitative trait loci (QTL) for spike-related traits is important for improvement of yield potential in wheat breeding. In this study, by using the wheat 55K single nucleotide polymorphism (SNP) array and diversity array technology (DArT), two recombinant inbred line populations derived from crosses avocet/chilero and avocet/huites were used to map QTL for kernel number per spike (KNS), total spikelet number per spike (TSS), fertile spikelet number per spike (FSS), and spike compactness (SC). Forty-two QTLs were identified on chromosomes 2A (4), 2B (3), 3A (2), 3B (7), 5A (11), 6A (4), 6B, and 7A (10), explaining 3.13-21.80% of the phenotypic variances. Twelve QTLs were detected in multi-environments on chromosomes 2A, 3B (2), 5A (4), 6A (3), 6B, and 7A, while four QTL clusters were detected on chromosomes 3A, 3B, 5A, and 7A. Two stable and new QTL clusters, QKns/Tss/Fss/SC.haust-5A and QKns/Tss/Fss.haust-7A, were detected in the physical intervals of 547.49-590.46 Mb and 511.54-516.15 Mb, accounting for 7.53-14.78% and 7.01-20.66% of the phenotypic variances, respectively. High-confidence annotated genes for QKns/Tss/Fss/SC.haust-5A and QKns/Tss/Fss.haust-7A were more highly expressed in spike development. The results provide new QTL and molecular markers for marker-assisted breeding in wheat.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01401-4.

Identification and validation of a major quantitative trait locus for spike length and compactness in the wheat (Triticum aestivum L.) line Chuanyu12D7

Spike length (SL) and spike compactness (SC) are crucial traits related to wheat (Triticum aestivum L.) yield potential. In this study, a backcrossed inbred lines (BILs) population segregating for SL/SC was developed by using a commercial variety chuanyu25 as recurrent parent and a backbone parent Chuanyu12D7. Bulked segregant analysis (BSA) combined with the Wheat 660K SNP array was performed to conduct quantitative trait locus (QTL) mapping. A major and stable SL/SC QTL (designated as QSl/Sc.cib-2D.1) was identified on chromosome 2DS, explaining 45.63-59.72% of the phenotypic variation. QSl/Sc.cib-2D.1 was mapped to a 102.29-Kb interval by flanking SNPs AX-110276364 and AX-111593853 using a BC₄F_2:3 population. Since QSl/Sc.cib-2D.1 is linked to the Rht8 gene, their additive effects on plant type and spike type were analysed. Remarkably, the superior allele of QSl/Sc.cib-2D.1 combined with Rht8 can increase SL and TGW, and decrese SC without any apparent trade-offs in other yield-related traits. In addition, the closely linked kompetitive allele-specific PCR (KASP) markers of this locus were developed for marker-assisted selection (MAS) breeding. Four genes within the physical interval were considered as potential candidates based on expression patterns as well as orthologous gene functions. These results laid the foundation for map-based cloning of the gene(s) underlying QSl/Sc.cib-2D.1 and its potential application in wheat ideotype breeding.

Construction of a high-density genetic map and localization of grazing-tolerant QTLs in Medicago falcata L

BACKGROUND

Using genomic DNA from 79 F1 plants resulted from a crossing between parents with strong and weak grazing tolerance in Medicago falcata L., we generated an EcoRI restriction site-associated DNA (RAD) sequencing library. After sequencing and assembly, a high-density genetic map with high-quality SNP markers was constructed, with a total length of 1312.238 cM and an average density of 0.844 SNP/cM.

METHODS

The phenotypic traits of 79 F1 families were observed and the QTLS of 6 traits were analyzed by interval mapping.

RESULTS

Sixty three QTLs were identified for seven traits with LOD values from 3 to 6 and the contribution rates from 15% to 30%. Among the 63 QTLs, 17 were for natural shoot height, 12 for rhizome Length, 10 for Shoot canopy diameter, 9 for Basal plant diameter, 6 for stem number, 5 for absolute shoot height, and 4 for rhizome width. These QTLs were concentrated on LG2, LG4, LG5, LG7, and LG8. LG6 had only 6 QTLs. According to the results of QTL mapping, comparison of reference genomes, and functional annotation, 10 candidate genes that may be related to grazing tolerance were screened. qRT-PCR analysis showed that two candidate genes (LOC11412291 and LOC11440209) may be the key genes related to grazing tolerance of M. falcata.

CONCLUSION

The identified trait-associated QTLs and candidate genes in this study will provide a solid foundation for future molecular breeding for enhanced grazing-tolerance in M. falcata.

RAD-Seq-Based High-Density Linkage Maps Construction and Quantitative Trait Loci Mapping of Flowering Time Trait in Alfalfa (Medicago sativa L.)

Alfalfa (Medicago sativa L.) is a perennial forage crop known as the "Queen of Forages." To dissect the genetic mechanism of flowering time (FT) in alfalfa, high-density linkage maps were constructed for both parents of an F1 mapping population derived from a cross between Cangzhou (P1) and ZhongmuNO.1 (P2), consisting of 150 progenies. The FT showed a transgressive segregation pattern in the mapping population. A total of 13,773 single-nucleotide polymorphism markers was obtained by using restriction-site associated DNA sequencing and distributed on 64 linkage groups, with a total length of 3,780.49 and 4,113.45 cM and an average marker interval of 0.58 and 0.59 cM for P1 and P2 parent, respectively. Quantitative trait loci (QTL) analyses were performed using the least square means of each year as well as the best linear unbiased prediction values across 4 years. Sixteen QTLs for FT were detected for P1 and 22 QTLs for P2, accounting for 1.40-16.04% of FT variation. RNA-Seq analysis at three flowering stages identified 5,039, 7,058, and 7,996 genes that were differentially expressed between two parents, respectively. Based on QTL mapping, DEGs analysis, and functional annotation, seven candidate genes associated with flowering time were finally detected. This study discovered QTLs and candidate genes for alfalfa FT, making it a useful resource for breeding studies on this essential crop.

Identification and validation of QTLs for kernel number per spike and spike length in two founder genotypes of wheat

BACKGROUND

Kernel number per spike (KNS) and spike length (SL) are important spike-related traits in wheat variety improvement. Discovering genetic loci controlling these traits is necessary to elucidate the genetic basis of wheat yield traits and is very important for marker-assisted selection breeding.

RESULTS

In the present study, we used a recombinant inbred line population with 248 lines derived from the two founder genotypes of wheat, Bima4 and BainongAK58, to construct a high-density genetic map using wheat 55 K genotyping assay. The final genetic linkage map consists of 2356 bin markers (14,812 SNPs) representing all 21 wheat chromosomes, and the entire map spanned 4141.24 cM. A total of 7 and 18 QTLs were identified for KNS and SL, respectively, and they were distributed on 11 chromosomes. The allele effects of the flanking markers for 12 stable QTLs, including four QTLs for KNS and eight QTLs for SL, were estimated based on phenotyping data collected from 15 environments in a diverse wheat panel including 384 elite cultivars and breeding lines. The positive alleles at seven loci, namely, QKns.his-7D2-1, QKns.his-7D2-2, QSl.his-4A-1, QSl.his-5D1, QSl.his-4D2-2, QSl.his-5B and QSl.his-5A-2, significantly increased KNS or SL in the diverse panel, suggesting they are more universal in their effects and are valuable for gene pyramiding in breeding programs. The transmission of Bima4 allele indicated that the favorite alleles at five loci (QKns.his-7D2-1, QSl.his-5A-2, QSl.his-2D1-1, QSl.his-3A-2 and QSl.his-3B) showed a relatively high frequency or an upward trend following the continuity of generations, suggesting that they underwent rigorous selection during breeding. At two loci (QKns.his-7D2-1 and QSl.his-5A-2) that the positive effects of the Bima4 alleles have been validated in the diverse panel, two and one kompetitive allele-specific PCR (KASP) markers were further developed, respectively, and they are valuable for marker-assisted selection breeding.

CONCLUSION

Important chromosome regions controlling KNS and SL were identified in the founder parents. Our results are useful for knowing the molecular mechanisms of founder parents and future molecular breeding in wheat.

Spike Density Quantitative Trait Loci Detection and Analysis in Tetraploid and Hexaploid Wheat Recombinant Inbred Line Populations

Spike density (SD) is an agronomically important character in wheat. In addition, an optimized spike structure is a key basis for high yields. Identification of quantitative trait loci (QTL) for SD has provided a genetic basis for constructing ideal spike morphologies in wheat. In this study, two recombinant inbred line (RIL) populations (tetraploid RIL AM and hexaploid RIL 20828/SY95-71 (2SY)) previously genotyped using the wheat55K SNP array were used to identify SD QTL. A total of 18 QTL were detected, and three were major and one was stably expressed (QSd.sau-2SY-7A.2, QSd.sau-AM-5A.2, QSd.sau-AM-7B, and QSd.sau-2SY-2D). They can explain up to 23.14, 19.97, 12.00, and 9.44% of phenotypic variation, respectively. QTL × environment and epistatic interactions for SD were further analyzed. In addition, pyramiding analysis further revealed that there were additive effects between QSd.sau-2SY-2D and QSd.sau-2SY-7A.2 in 2SY, and QSd.sau-AM-5A.2 and QSd.sau-AM-7B in AM. Pearson's correlation between SD and other agronomic traits, and effects of major or stable QTL on yield related traits indicated SD significantly impacted spike length (SL), spikelet number per spike (SNS) and kernel length (KL). Several genes related to spike development within the physical intervals of major or stable QTL were predicted and discussed. Collectively, our research identified QTL with potential applications for modern wheat breeding and broadening the genetic basis of SD.

Identification and validation of two major QTLs for spike compactness and length in bread wheat (Triticum aestivum L.) showing pleiotropic effects on yield-related traits

Two major and stable QTLs for spike compactness and length were detected and validated in multiple genetic backgrounds and environments, and their pleiotropic effects on yield-related traits were analyzed. Spike compactness (SC) and length (SL) are greatly associated with wheat (Triticum aestivum L.) grain yield. To detect quantitative trait loci (QTL) associated with SC and SL, two biparental populations derived from crosses of Chuanmai42/Kechengmai1 and Chuanmai42/Chuannong16 were employed to perform QTL mapping in five environments. A total of 34 QTLs were identified, in which six major QTLs were repeatedly detected in more than four environments and the best linear unbiased prediction datasets, explaining 7.13-33.6% of phenotypic variation. These major QTLs were co-located in two genomic regions on chromosome 5A and 6A, namely QSc/Sl.cib-5A and QSc/Sl.cib-6A, respectively. By developing kompetitive allele-specific PCR (KASP) markers that linked to them, the two loci were validated in different genetic backgrounds, and their interactions were also analyzed. Comparison analysis showed that QSc/Sl.cib-5A was not Vrn-A1 and Q, and QSc/Sl.cib-6A was likely a new locus for SC and SL. Both QSc/Sl.cib-5A and QSc/Sl.cib-6A had pleiotropic effects on other yield-related traits including plant height, thousand grain weight and grain length. Therefore, the two loci combined with the developed KASP markers might be potentially applicable in wheat breeding. Furthermore, based on the spatiotemporal expression patterns, gene annotation, orthologous search and sequence differences, TraesCS5A01G301400 and TraesCS6A01G090300 were considered as potential candidates for QSc/Sl.cib-5A and QSc/Sl.cib-6A, respectively. These results provided valuable information for fine mapping and cloning of the two loci in the future.

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

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