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

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.

A major QTL simultaneously increases the number of spikelets per spike and thousand-kernel weight in a wheat line

KEY MESSAGE

A novel and stably expressed QTL QSNS.sicau-SSY-7A for spikelet number per spike in wheat without negative effects on thousand-kernel weight was identified and validated in different genetic backgrounds. Spikelet number per spike (SNS) is an important determinant of yield in wheat. In the present study, we combined bulked segregant analysis (BSA) and the wheat 660 K single-nucleotide polymorphism (SNP) array to rapidly identify genomic regions associated with SNS from a recombinant inbred line (RIL) population derived from a cross between the wheat lines S849-8 and SY95-71. A genetic map was constructed using Kompetitive Allele Specific PCR markers in the SNP-enriched region on the long arm of chromosome 7A. A major and stably expressed QTL, QSNS.sicau-SSY-7A, was detected in multiple environments. It was located in a 1.6 cM interval on chromosome arm 7AL flanked by the markers AX-109983514 and AX-109820548. This QTL explained 6.86-15.72% of the phenotypic variance, with LOD values ranging from 3.66 to 8.66. Several genes associated with plant growth and development were identified in the interval where QSNS.sicau-SSY-7A was located on the 'Chinese Spring' wheat and wild emmer reference genomes. Furthermore, the effects of QSNS.sicau-SSY-7A and WHEAT ORTHOLOG OFAPO1(WAPO1) on SNS were analyzed. Interestingly, QSNS.sicau-SSY-7A significantly increased SNS without negative effects on thousand-kernel weight, anthesis date and plant height, demonstrating its great potential for breeding aimed at improving grain yield. Taken together, these results indicate that QSNS.sicau-SSY-7A is a promising locus for yield improvement, and its linkage markers are helpful for fine mapping and molecular breeding.

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.

Major and stably expressed QTL for traits related to the mature wheat embryo independent of kernel size

Two major and stably expressed QTL for traits related to mature wheat embryo independent of kernel size were identified and validated in a natural population that contained 171 Sichuan wheat accessions and 49 Sichuan wheat landraces. As the juvenile of a highly differentiated plant, mature wheat (Triticum aestivum L.) embryos are highly significant to agricultural production. To understand the genetic basis of traits related to wheat embryo size, the embryo of mature kernels in a recombination inbred line that contained 126 lines from four environments was measured. The genetic loci of embryo size, including embryo length (EL), embryo width (EW), embryo area (EA), embryo length/kernel length (EL/KL), embryo width/kernel width (EW/KW), and EL/EW, were identified based on a genetic linkage map constructed based on PCR markers and the Wheat 55 K single nucleotide polymorphism (SNP) array. A total of 50 quantitative trait loci (QTL) for traits related to wheat embryo size were detected. Among them, QEL.sicau-2SY-4A for EL and QEW.sicau-2SY-7B for EW were major and stably expressed and were genetically independent of KL and KW, respectively. Their effects were further verified in a natural population that contained 171 Sichuan wheat accessions and 49 Sichuan wheat landraces. Further analysis showed that TraesCS4A02G343300 and TraesCS7B02G006800 could be candidate genes for QEL.sicau-2SY-4A and QEW.sicau-2SY-7B, respectively. In addition, significant positive correlations between EL and kernel-related traits and the 1,000-grain weight were detected. Collectively, this study broadens our understanding of the genetic basis of wheat embryo size and will be helpful for the further fine-mapping of interesting loci in the future.

Genome-wide association analysis of spike and kernel traits in the U.S. hard winter wheat

A better understanding of the genetic control of spike and kernel traits that have higher heritability can help in the development of high-yielding wheat varieties. Here, we identified the marker-trait associations (MTAs) for various spike- and kernel-related traits in winter wheat (Triticum aestivum L.) through genome-wide association studies (GWAS). An association mapping panel comprising 297 hard winter wheat accessions from the U.S. Great Plains was evaluated for eight spike- and kernel-related traits in three different environments. A GWAS using 15,590 single-nucleotide polymorphisms (SNPs) identified a total of 53 MTAs for seven spike- and kernel-related traits, where the highest number of MTAs were identified for spike length (16) followed by the number of spikelets per spike (15) and spikelet density (11). Out of 53 MTAs, 14 were considered to represent stable quantitative trait loci (QTL) as they were identified in multiple environments. Five multi-trait MTAs were identified for various traits including the number of spikelets per spike (NSPS), spikelet density (SD), kernel width (KW), and kernel area (KA) that could facilitate the pyramiding of yield-contributing traits. Further, a significant additive effect of accumulated favorable alleles on the phenotype of four spike-related traits suggested that breeding lines and cultivars with a higher number of favorable alleles could be a valuable resource for breeders to improve yield-related traits. This study improves the understanding of the genetic basis of yield-related traits in hard winter wheat and provides reliable molecular markers that will facilitate marker-assisted selection (MAS) in wheat breeding programs.

QTL Mapping for Important Agronomic Traits Using a Wheat55K SNP Array-Based Genetic Map in Tetraploid Wheat

Wheat yield is highly correlated with plant height, heading date, spike characteristics, and kernel traits. In this study, we used the wheat55K single nucleotide polymorphism array to genotype a recombinant inbred line population of 165 lines constructed by crossing two tetraploid wheat materials, Icaro and Y4. A genetic linkage map with a total length of 6244.51 cM was constructed, covering 14 chromosomes of tetraploid wheat. QTLs for 12 important agronomic traits, including plant height (PH), heading date (HD), awn color (AC), spike-branching (SB), and related traits of spike and kernel, were mapped in multiple environments, while combined QTL-by-environment interactions and epistatic effects were analyzed for each trait. A total of 52 major or stable QTLs were identified, among which may be some novel loci controlling PH, SB, and kernel length-width ratio (LWR), etc., with LOD values ranging from 2.51 to 54.49, thereby explaining 2.40-66.27% of the phenotypic variation. Based on the 'China Spring' and durum wheat reference genome annotations, candidate genes were predicted for four stable QTLs, QPH.nwafu-2B.2 (165.67-166.99 cM), QAC.nwafu-3A.1 (419.89-420.52 cM), QAC.nwafu-4A.1 (424.31-447.4 cM), and QLWR.nwafu-7A.1 (166.66-175.46 cM). Thirty-one QTL clusters and 44 segregation distortion regions were also detected, and 38 and 18 major or stable QTLs were included in these clusters and segregation distortion regions, respectively. These results provide QTLs with breeding application potential in tetraploid wheat that broadens the genetic basis of important agronomic traits such as PH, HD, AC, SB, etc., and benefits wheat breeding.

Identification and validation of a locus for wheat maximum root length independent of parental reproductive environment

Maximum root length (MRL) plays an important role in the uptake of nutrients and resisting abiotic stresses. Understanding the genetic mechanism of root development is of great significance for genetic improvement of wheat. Previous studies have confirmed that parental reproductive environment (PRE) has a significant impact on growth and development of the next generation in the whole life cycle of a given plant. In this study, a recombinant inbred line population genotyped using the Wheat55K SNP array, was used to map quantitative trait loci (QTL) for wheat seedling MRL based on the harvested seeds from five different PREs. A total of 5 QTL located on chromosomes 3D and 7A were identified. Among them, QMrl.sicau-2SY-3D.2 located in a 4.0 cM interval on chromosome 3D was likely independent of PREs. QMrl.sicau-2SY-7A.2 was detected in two tests and probably influenced by PREs. The effect of QMrl.sicau-2SY-3D.2 was further validated using the tightly linked kompetitive allele specific PCR (KASP) marker, KASP-AX-111589572, in populations with different genetic backgrounds. Lines with a combination of positive alleles from QMrl.sicau-2SY-3D.2 and QMrl.sicau-2SY-7A.2 have significantly longer MRL. Furthermore, four genes (TraesCS3D03G0612000, TraesCS3D03G0608400, TraesCS3D03G0613600, and TraesCS3D03G0602400) mainly expressed in wheat root were predicted to be associated with root growth. Taken together, this study reports on a major QTL independent of PREs and lays a foundation for understanding the regulation mechanism of wheat MRL at the seedling stage.

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

Spike compactness (SC) and length (SL) are the components of spike morphology and are strongly related to grain yield in wheat (Triticum aestivum L.). To investigate quantitative trait loci (QTL) associated with SC and SL, a recombinant inbred lines (RIL) population derived from the cross of Bailangmai (BLM, a Tibet landrace) and Chuanyu 20 (CY20, an improved variety) was employed in six environments. Three genomic regions responsible for SC and SL traits were identified on chromosomes 2A and 2D using bulked segregant exome sequencing (BSE-Seq). By constructing genetic maps, six major QTL were repeatedly detected in more than four environments and the best linear unbiased estimation (BLUE) datasets, explaining 7.00-28.56% of the phenotypic variation and the logarithm of the odd (LOD) score varying from 2.50 to 13.22. They were co-located on three loci, designed as QSc/Sl.cib-2AS, QSc/Sl.cib-2AL, and QSc/Sl.cib-2D, respectively. Based on the flanking markers, their interactions and effects on the corresponding trait and other agronomic traits were also analyzed. Comparison analysis showed that QSc/Sl.cib-2AS and QSc/Sl.cib-2AL were possibly two novel loci for SC and SL. QSc/Sl.cib-2AS and QSc/Sl.cib-2D showed pleiotropic effects on plant height and grain morphology, while QSc/Sl.cib-2AL showed effects on spikelet number per spike (SNS) and grain width (GW). Based on the gene annotation, orthologous search, and spatiotemporal expression patterns of genes, TraesCS2A03G0410600 and TraesCS2A03G0422300 for QSc/Sl.cib-2AS, and TraesCS2D03G1129300 and TraesCS2D03G1131500 for QSc/Sl.cib-2D were considered as potential candidate genes, respectively. These results will be useful for fine mapping and developing new varieties with high yield in the future.