Skip to main content

Alleles on locus chromosome 4B from different parents confer tiller number and the yield-associated traits in wheat

Abstract

Pleiotropy is frequently detected in agronomic traits of wheat (Triticum aestivum). A locus on chromosome 4B, QTn/Ptn/Sl/Sns/Al/Tgw/Gl/Gw.caas-4B, proved to show pleiotropic effects on tiller, spike, and grain traits using a recombinant inbred line (RIL) population of Qingxinmai × 041133. The allele from Qingxinmai increased tiller numbers, and the allele from line 041133 produced better performances of spike traits and grain traits. Another 52 QTL for the eight traits investigated were detected on 18 chromosomes, except for chromosomes 5D, 6D, and 7B. Several genes in the genomic interval of the locus on chromosome 4B were differentially expressed in crown and inflorescence samples between Qingxinmai and line 041133. The development of the KASP marker specific for the locus on chromosome 4B is useful for molecular marker-assisted selection in wheat breeding.

Peer Review reports

Introduction

Wheat (Triticum aestivum L.) is a long-historically grown crop in more than 40 countries. It feeds about 30% of the world’s population. Despite China, as a leading wheat grower, producing approximately 17% of the global wheat, the demand for wheat is increasing due to the population growth and rapid urbanization [1,2,3]. Increase in wheat yield per annum from the 1920s was estimated to be 1.29% and 1.50% for north and south winter wheat in China [4]. Nevertheless, this yield increment does not meet the future demand for wheat [5]. Further increase in grain yields is the highest priority in most wheat breeding programs throughout the country.

Wheat yield is determined by many agronomic traits, such as the number of spikes per unit area, grain numbers per spike (GNS), and thousand-grain weight (TGW) [6]. Because direct selection of grain yield is difficult, the improvement of yield-associated traits is often conducted instead. In fact, GNS and TGW have been increased in wheat cultivars released in the past several decades in China [4]. Many genetic loci governing grain weight were identified on different wheat chromosomes [7]. Several studies reported quantitative trait loci (QTL) for grain weight on the short arm of chromosome 4B. QTKW.caas-4BS for grain weight was mapped in a 483 kb genomic interval in Doumai, which harbors three genes ZnF, EamA, and Rht-B1 [8]. The functions of Rht-B1 (encoding gibberellin signaling repressor) and ZnF-B (encoding a RING-type E3 ligase) on grain yield were determined using natural deletion of a haploblock ( 500 kb) in wheat cultivar Heng597 [9]. Grain weight is associated with grain size measured by grain length (GL) and grain width (GW) [10,11,12,13]. Chen et al. (2020) [14] identified 30 stable QTL including QTgw.cau-7D and QGw.cau-7D for grain size on chromosome 7D.

Tiller is a branching phenomenon in monocotyledons [15, 16], which determines the number of spikes. A low-tillering gene, tin1, was localized on the short arm of wheat chromosome 1A [17, 18]. Several genes or QTL for increasing or decreasing tiller numbers were identified. Genes tin3 and TaD27 inhibit tiller formation [19, 20]. QMTN.sicau-4D and QPtn.sau-4B increase tiller numbers [21, 22]. QTa.sau-2B-769 for tiller angle was detected on the chromosome arm 2BL. TraesCS2B01G583800, a gene regulating leaf angle, was most likely the candidate gene for this QTL [23].

A wheat spike is composed of spikelets borne on its rachis. Spike length (SL) is closely correlated with spikelet number (SN). Both traits are influenced by interaction between genetic and environmental factors. Li et al. [24] mapped two major QTL for SL, QSc/Sl.cib-5A and QSc/Sl.cib-6 A, explaining 7.13–33.6% of the phenotypic variations. TaAPO-A1 confers total SN in European winter wheat cultivars [25]. QSns.sau-2D was detected on the chromosome arm 2DS and can explain 10.16–45.68% of phenotypic variations [26]. An important pleiotropic QTL, Q.SpnN/SpkLng/PH/SPP.3A, was identified on chromosome 3A. It was found to be associated with spikelet number per spike, spike length, plant height, and spikes per spike [27].

Some wheat cultivars have awns in glumes. Awn is an important organ for respiration and photosynthesis of spikes. It not only plays a role in protection and transmission of seeds, but also impacts on grain yield. Assimilates synthesized by wheat awns are transported to grains nearby [28]. It is believed that B1, B2, and HD are the main genes inhibiting the awn development. Locus B1 (tipped 1) on chromosome 5AL regulates the phenotype of tip awns, and B2 (tipped 2) on chromosome 6BL shortens the awn length [29, 30]. Hd, a dominant allele of knotted 1 (kn-1) on the chromosome arm 4AS, regulates the hook-awn phenotype of wheat [31].

Wheat landrace Qingxinmai is featured by a large number of tillers. It can develop 30–40 spikes per plant. Breeding line 041133 develops large spikes and grains. Incorporation of favorable alleles of loci for tiller number, spike, and grain traits from Qingxinmai and line 041133 is an option to increase yield of wheat. Toward this end, dissection of QTL for the traits of interest is a prerequisite. A recombinant inbred line (RIL) population was developed from a cross between these genotypes. The purpose of this study was to unravel the genetic control of the agronomic traits with the aid of bulk segregant analysis-RNA-Seq (BSR-Seq), bulked segregant exome capture sequencing (BSE-Seq), and a wheat 16 K genotyping by target sequencing (GBTS) single nucleotide polymorphism (SNP) array.

Results

Phenotypic performances

Qingxinmai and line 041133 differed significantly in tiller traits, total tiller number (TN) and productive tiller number (PTN), spike traits, spike length (SL), spikelet number per spike (SNS), and awn length (AL) (including top, central, and bottom positions on spikes), and grain traits, thousand-grain weight (TGW), grain length (GL), and grain width (GW) (P < 0.05). Qingxinmai had more tillers and longer awns, while line 041133 had longer spikes and larger kernels in different field trials (Fig. 1; Table 1). The broad-sense heritability (H2) for these traits ranged from 0.40 to 0.99.

Fig. 1
figure 1

Plants (a), inflorescences (b), and grains (c) of Qingxinmai (left) and 041133 (right)

Table 1 Phenotypic variation of total tiller number (TN), productive tiller number (PTN), awn length (AL), spike length (SL), spikelet number per spike (SNS), grain width (GW), grain length (GL), and thousand-grain weight (TGW) of the parents and the RIL population in different environments

Variations were observed in traits of the Qingxinmai × 041133 RIL population investigated. The coefficients of variation (CV) of the tiller and spike traits (i.e., TN, PTN, SL, SNS, and AL) was greater than that of the grain traits (i.e., TGW, GL, and GW) (Table 1). A nearly normal frequency distribution of TN, PTN, SNS, SL, TGW, GL, and GW for the RILs was observed in each environment and the best linear unbiased estimate (BLUE) datasets with the absolute values of Skewness and Kurtosis coefficients approaching 0, except for TN and PTN obtained at 2020CP. (Table 1, Figure S1), suggesting the quantitative inheritance controlled by multiple loci. Awn lengths measured in the three positions on spikes showed a bimodal frequency distribution with the Kurtosis coefficients > 1, indicating the presence of major gene for these traits.

Correlation analysis demonstrated that the grain traits (TGW, GL, and GW) and the spike traits (SL and SNS) were positively correlated (r = 0.20–0.77, P < 0.01), but they were negatively correlated with TN and PTN (r = 0.14–0.55, P < 0.05) (Figure S2). TN and PTN were positively correlated with each other (r = 0.57, P < 0.01). Lengths of awn on the top, centre, and bottom of spikes were correlated (P < 0.01), and they were also correlated with GL, GW and SNS (P < 0.05). Except for TN, the other 7 traits were significantly correlated among different environments (P < 0.01) (Table S1). A significant interaction was observed between genotypes and environments for TN and PTN, while no such interaction was detected for the other traits (Table S2).

BSE-Seq analysis

The statistical parameters from the BSE-Seq analysis are summarized in Table S3. A total of 10,774 SNP variants were detected between the high- and low-TGW DNA pools (Bulk-HTGW and Bulk-LTGW). The most abundant enrichment of the TGW-associated SNPs and InDels was observed in the genomic intervals of 28.60–206.77 Mb and 342.32–621.14 Mb on chromosome 4B (8075) and 33.00–87.24 Mb on chromosome 7 A (2699) of Chinese Spring (CS) reference genome sequence RefSeq v1.0 (Fig. 2c, Table S3).

Fig. 2
figure 2

Distribution of single nucleotide polymorphisms (SNPs) between the phenotypically contrasting bulks of the Qingxinmai × 041133 RIL population. (a) BSR-Seq analysis of the high- and low-tiller number RNA pools (Bulk-HTN and Bulk-LTN); (b) BSR-Seq analysis of the long- and short-spike RNA pools (Bulk-LS and Bulk-SS); and (c) exon trapping analysis of the high- and low-thousand-grain weight pools (Bulk-HTGW and Bulk-LTGW)

BSR-Seq analysis of crowns and inflorescences

The statistical parameters of the RNA sequencing and comparisons between the RNA pools from contrasting tiller numbers (Bulk-HTN vs. Bulk-LTN) and between the RNA pools from contrasting spike lengths (Bulk-LS vs. Bulk-SS) are shown in Table S3. Raw reads of these RNA samples ranged from 135,474,852 to 156,857,352, and clean reads after trimming ranged from 129,072,674 to 149,787,196. The proportions of clean reads that were mapped to the CS reference genome RefSeq v1.0 were over 95%. The BSR-Seq analysis identified 1061 and 1508 SNP variants by comparing crown (Bulk-HTN vs. Bulk-LTN) and inflorescence (Bulk-LS vs. Bulk-SS) samples, respectively. Most SNP variants from tiller samples (Bulk-HTN vs. Bulk-LTN) were enriched on chromosome 4B (888), 6A (128), and 7A (45). SNPs from inflorescence samples (Bulk-LS vs. Bulk-SS) were mostly anchored on chromosomes 4B (1342), 4D (67), and 6A (99) (Fig. 2, Table S3). Most of the 464 common SNP variants detected in the tiller and inflorescence samples were located on chromosome 4B.

Construction of the genetic linkage map

The 16 K GBTS SNP array applied for the RIL population of Qingxinmai × 041133 generated 37,699 original SNPs. After removing those with the coverage depth < 5×, 14,868 SNPs were retained. Polymorphic SNPs between parents Qingxinmai and line 041133 were 4,939 in 2,430 loci. A total of 2,398 loci with clear positions on the reference genome were used to construct a genetic linkage map (3113.1 cM). Thirty linkage groups were established. Most chromosomes had single linkage groups, except for two for chromosomes 1D, 3D, 5D, 6B, and 7D and three for chromosomes 6A and 6D. The marker density of this genetic linkage map was 3.27 cM per locus and 1.30 cM per marker. Subgenomes A, B, and D consisted of 344, 455, and 154 loci with the map distances of 1179.2, 1282.2, and 651.7 cM, respectively (Table S4). The genetic linkage maps constructed ranged from 50.9 cM (6D) to 197.3 cM (7A). Each chromosome contained 12 (4D) to 232 (3B) SNP markers. The average distance between adjacent SNPs and bin markers was in a range of 0.79 (3A) to 8.13 cM (3D) and 2.12 (2B) to 8.81 cM (3D). The marker orders on most chromosomes were generally consistent with those in the CS reference genome sequence RefSeq v1.0 (Figure S3).

The Chi-squared test revealed a genetic distortion in 377 polymorphic molecular markers (P < 0.05). Among them, 219 (58.09%) and 158 (41.91%) markers were biased towards line 041133 and Qingxinmai, respectively. Fifteen segregation distortion regions (SDR, ≥3SD loci) were detected in the RIL population. Eleven SDR originated from Qingxinmai and 4 from line 041133 (Table S5).

QTL mapping

Sixty QTL for the 8 traits investigated were detected on 18 chromosomes, except for 5D, 6D, and 7B. Nine QTL were detected in multiple environments and the BLUE data sets, explaining phenotypic variations ranging from 3.43 to 65.34%. Fifteen QTL were detected in one or two environments and the BLUE data sets, explaining phenotypic variations ranging from 1.53 to 15.27%. The resting 36 QTL were observed in single environments, explaining 3.15–19.41% of the phenotypic variations (Table 2, Figure S4-S9).

Table 2 Quantitative trait loci (QTL) for total tiller number (TN), productive tiller number (PTN), spikelet number per spike (SNS), spike length (SL), awn length (AL), thousand-grain weight (TGW), grain length (GL), and grain width (GW) identified from different environments in the Qingxinmai × 041133 RIL population

Tiller traits

Twelve QTL conferring TN and PTN were identified on chromosomes 1A, 2D, 3B, 4B, 5A, 5B, and 7A (Figure S4-S9). A major locus, QTn.caas-4B.1/QPtn.caas-4B, was localized in the same genetic interval on chromosome 4B. It explained 13.75–32.41% of the phenotypic variations for TN and PTN with logarithm of odds (LOD) values of 7.26–20.49. The positive alleles of this locus for more tillers were contributed by Qingxinmai. Another minor effective QTL, QPtn.caas-2D for PTN, contributed by line 041133, was observed in one environment (2020XX) and the BLUE datasets. It explained 4.00% and 4.11% of the phenotypic variations. The remaining 9 loci were detected in single environments and explained 3.98–11.73% of the phenotypic variations. Most of them were contributed by Qingxinmai, except for two loci (QTn.caas-1A and QPtn.caas-5B) by line 041133 (Table 2).

Spike and awn traits

Nineteen QTL for SL, SNS, and AL were identified on chromosomes 1A, 1B, 1D, 2A, 2D, 3A, 4A, 4B, 5A, 6A, and 7D (Figure S4-S9). A major locus, QSl.caas-4B/QSns.caas-4B, was detected on chromosome 4B in four environments and the BLUE datasets. This QTL was contributed by line 041133 and explained the phenotypic variations of SL and SNS by 6.7–36.07% (LOD = 3.46–21.89). QSl.caas-2A, QSl.caas-2D, QSl.caas-4A, QSl.caas-6A for SL, and QSns.caas-2A, QSns.caas-2D.1, and QSns.caas-2D.2 for SNS were detected on different chromosomes in one or two environments and the BLUE datasets. They explained 3.26–8.43% of the phenotypic variations (LOD = 3.04 to 6.65). The remaining seven QTL for spike and awn traits were observed in single environments, accounting for 4.01–12.34% of the phenotypic variations (Table 2).

QAl.caas-4B for AL detected in 2021CP and the BLUE dataset was located in the same genetic interval as QSl.caas-4B/QSns.caas-4B. This QTL was contributed by 041133. But ALs on the top, central, and bottom of spikes appeared to be controlled by a major locus QAl.caas-5A. It was detected in all the three environments and the BLUE dataset, and explained the phenotypic variations of 59.96–65.34% (LOD = 56.12–71.09). The physical location of this locus was observed in a genomic interval of 688.17–697.64 Mb (Table 2). The additive effect of QAl.caas-5A was provided by Qingxinmai. The SNP sequence of the closest molecular marker 5A_688174490 obtained by the 16 K GBTS SNP array was converted to a Kompetitive allele-specific PCR (KASP) marker KASP_ 5A_ 688174490 (Fig. 3, Table S6). This KASP marker proved to be tightly linked to QAl.caas-5A by genotyping the entire RIL population (Figure S10a).

Fig. 3
figure 3

Genetic linkage map (a), QTL analysis method (b), and effects (c) of locus QAL.caas-5A on awn lengths. The BLUE values for awn lengths of the Qingxinmai × 041133 RILs were grouped based on the genotypes of the locus-specific KASP marker KASP_5A_ 688174490. **, P < 0.01

Grain traits

Twenty-nine QTL for TGW, GW, and GL were detected on chromosomes 1A, 1B, 2A, 2B, 3A, 3B, 3D, 4A, 4B, 4D, 5A, 6A, 6B, 7A, and 7D (Figure S4-S9). Nine loci on chromosomes 1B, 2A, 2B, 3B, 3D, and 4B were identified in multiple environments. The major QTL, QTgw/Gl/Gw.caas-4B for the three grain traits, were detected in all environments and the BLUE datasets. This locus was contributed by line 041133 and explained 6.47-25.30% of the phenotypic variations (LOD = 3.62–14.77). Another five loci for grain weight, QGl.caas-1B, QGl.caas-2B, QGl.caas-4B.3QGw.caas-2A, and QGw.caas-3D, were identified in two environments or the BLUE dataset, explaining 4.72–14.64% of the phenotypic variations. The favorable alleles of these QTL were also contributed by line 041133, except for QGl.caas-1B and QGl.caas-2B by Qingxinmai. QGl.caas-3B for GL on chromosome 3B were detected in two environments at sites 2021BJ and 2021QS and the BLUE datasets, explaining 5.26–11.37% of the phenotypic variations. The favorable allele of this locus was contributed by Qingxinmai. The remaining 20 minor QTL for the grain traits were detected in single environments, accounting for 3.15–19.41% of the phenotypic variations. Line 041133 contributed 16 positive alleles and Qingxinmai contributed the other 4 positive alleles (Table 2).

Analysis of additive effects of the major QTL

Since the pleiotropic QTL for PTN, SL, SNS, TGW, GW, and GL were detected on chromosome 4B, we analyzed the additive effects of this locus on the corresponding traits using the BLUE datasets of the mapping population. Lines with the favorable alleles of QPtn/Sl/Sns/Tgw/Gl/Gw.caas-4B only increased PTN, TGW, GL, and GW by 17.03%, 6.85%, 2.80%, and 4.23% over those without the alleles (Fig. 4). The addition of alleles at the minor loci QPtn.caas-2D, QTgw.caas-3D, QGw.caas-3D, and QGl.caas-3B further enhanced the performances of those traits by 20.79%, 11.26%, 4.45%, and 5.68%, respectively (Fig. 4). Lines with the favorable allele for SL at QSl.caas-4B only increased SL by 9.84% relative to lines without the favorable allele. The trait performances of SL and SNS appeared to be associated with the number of positive alleles. More favorable alleles increased the trait values by different magnitudes (Figure S11).

Fig. 4
figure 4

Additive effects of the QTL detected on productive tiller number (PTN) (a), thousand-grain weight (TGW) (b), grain width (GW) (c), and grain length (GL) (d) using the BLUE datasets of the Qingxinmai × 041133 RIL population. + and −: presence and absence of the favorable alleles of the target QTL based on genotypes of the flanking markers of the corresponding QTL. **, P < 0.01. ns: no significant difference

Development of the KASP marker specific for QTn/Ptn/Sl/Sns/Tgw/Gl/Gw.caas-4B

Locus QTn/Ptn/Sl/Sns/Tgw/Gl/Gw.caas-4B for the tiller, spike, and grain traits was anchored in the physical interval of 22.92–32.17 Mb on chromosome 4B of the CS reference genome sequence RefSeq v1.0. A KASP marker, KASP_4B_32174878, was converted from the SNP locus 4B_32174878 (Table S6, Figure S10b). Two previously developed QTkw.caas-4BS-linked markers were also used to genotype Qingxinmai and line 041133. Marker AX_89323611 exhibited monomorphism and the other marker AX_110956837 was polymorphic. The two polymorphic markers KASP_4B_32174878 and AX_110956837 placed locus QTn/Ptn/Sl/Sns/Tgw/Gl/Gw.caas-4B in a physical interval of 3.23 Mb (28.95–32.17 Mb) in the CS reference genome RefSeq v1.0 (Fig. 5).

Fig. 5
figure 5

a: Genetic linkage map (a), QTL analysis method (b), and effects (c) of the pleiotropic locus QTn/Ptn/Sl/Sns/Tgw/Gl/Gw.caas-4B on total tiller number (TN), productive tiller number (PTN), spike length (SL), spikelet number per spike (SNS), thousand-grain weight (TGW), grain length (GL), and grain width (GW). The BLUE datasets of the Qingxinmai × 041133 RILs were grouped based on the genotypes of the locus-specific KASP marker KASP_4B_32174878. **, P < 0.01

Analysis of the annotated genes in the target genomic interval of QTn/Ptn/Sl/Sns/Tgw/Gl/Gw.caas-4B

We analyzed the genes in the genomic region where the QTL on chromosome 4B resides using the in silico expression method and the RNA-Seq data generated from the BSR-Seq analysis with crowns and inflorescences. Twenty-two high confidential genes were annotated in the genomic interval of the pleiotropic QTL on chromosome 4B of the Chinese Spring reference genome RefSeq v1.0 (Table S7). The in silico expression of these annotated genes was analyzed in the Hexaploid Wheat Expression Database (IWGSC Annotation v1.1) assembled in the Triticeae Multi-Omics Center (http://202.194.139.32/). Five genes TraesCS4B01G042300 (pleckstrin homology domain), TraesCS4B01G042900 (ZnF), TraesCS4B01G043100 (Rht-B1b), TraesCS4B01G043400 (phytanoyl-CoA dioxygenase), and TraesCS4B01G044300 (microsomal glutathione S-transferase 3) were expressed in spikes and grains (Figure S12a). TraesCS4B01G042900, TraesCS4B01G043100, and TraesCS4B01G043300 were also differentially expressed between Qingxinmai and line 041133 in the BSR-Seq analysis with the crown and inflorescence RNA samples (Figure S12b).

The expression of these genes in crowns and inflorescences of Qingxinmai and line 041133 were further determined by qPCR (Figure S13). Gene TraesCS4B01G042900 was differentially expressed in crowns but not in inflorescences, TraesCS4B01G043100 (Rht-B1b) in both tissues, and TraesCS4B01G043300 in inflorescences but not in crowns. The expression of the three genes was higher in Qingxinmai than in line 041133, except TraesCS4B01G043100 in crowns.

Discussion

We detected a pleiotropic locus on chromosome 4B using a RIL population derived from a wheat landrace Qingxinmai and a breeding line 041133. We further determined that the alleles of this locus from the two parents confer different traits. The allele from Qingxinmai was responsible for increasing TN and PTN. The allele from line 041133 increased spike traits (SL and SNS) and grain traits (TGW, GL, and GW), and even had a minor effect on increasing AL. The effects of this locus on the traits investigated were enhanced by several minor effective QTL, such as QPtn.caas-2D, QTgw.caas-3D, QGw.caas-3D, and QGl.caas-3B.

Many QTL for various wheat plant growth and yield-related traits have been characterized with different mapping populations and various types of molecular markers [7, 32]. Wheat chromosome 4B was associated with tiller number [22, 33, 34], spike length and spikelet number [35,36,37], and grain weight in separate studies [8, 9, 36]. Yet some of those loci may not be localized in the same genomic regions of chromosome 4B. A QTL for tiller number was localized at 482.82 Mb [33]. Four loci for PTN were located on the genomic regions of 75.74-640.97 Mb [34, 38]. Liu et al. [38] and Deng et al. [39] identified a major QTL for PTN at 256.31 Mb. The allele of QTn/Ptn.caas-4B from Qingxinmai was anchored in the 28.95–32.17 Mb genomic region, which appears to be different from the above loci for tiller number. But the genomic location of locus QTn/Ptn.caas-4B overlaps with QPtn.sau-4B (28,941,377 and 32,167,076 bp) on chromosome 4B [22]. QPtn.sau-4B was detected in Chuannong 16, a spring wheat cultivar developed in Sichuan province, which is unrelated to Qingxinmai.

Wheat spike traits are also associated with locus on chromosome 4BS. QSL.caas-4BS was mapped to a physical position of 25.80–46.60 Mb in the Linmai 2 × Zhong 892 RIL population [36]. The association of chromosome 4B and grain weight was reported in a RIL population of Doumai/Shi 4185 [36]. In that study, one of the 11 QTL for grain weight QTkw.caas-4BS explained a high range of phenotypic variation (12.1–45.6%). That locus was located in a genomic region of 25.80–46.60 Mb. A 483-kb deletion in this region in Doumai contains genes ZnF, EamA, and Rht-B1 [8]. In a most recent study, Song et al. [9] reported that wheat cultivar Heng 597 possessed a locus QTgw.cau-4B for grain weight. The deletion of approximately 500 kb fragment, also carrying these three genes, increased grain weight. The knockdown of Rht-B1b in Fielder increases plant height, spike length, and grain weight. But deletion of ZnF-B led to a slight reduction in grain size and plant height with no change in spike length compared to the wild-type Fielder. There is no large fragment deletion in the same genomic interval in Qingxinmai and line 041133 as in Doumai and Heng 597. Further study is needed to characterize the genes associated with the pleiotropic locus on chromosome 4B in the current study.

Selection for awns with minimal extension, colloquially referred to as ‘awnletted’, has been dominated by genes Tipped1 (B1), Tipped2 (B2), and Hooded (Hd). These loci inhibit awn development of wheat [40, 41]. We detected a major effective locus QAl.caas-5A for awn inhibition in the tip-awned line 041133. The target interval of this locus was located on chromosome 5A (688.17–697.64 Mb), which overlaps gene B1 [42, 43]. B1 inhibits awn elongation by impeding cytokinin- and auxin-promoted cell division and directly repressing the expression of TaRAE2 and TaLks2 [44]. Based on the genomic locations, QAl.caas-5A is most likely identical to gene B1.

Segregation distortion of molecular markers is present in different genetic populations. Paillard et al. [45] detected 17% segregation distortion of RFLP and SSR markers in the RILs from cross between Arina and Forno wheats. We detected two large SDRs, one on chromosome 2B from Qingxinmai, and the other on chromosome 6B from line 041133. Two SDRs, SDR-4B.2 (26.9–30.8 cM) and SDR-4B.3 (30.8–34.4 cM) from Qingxinmai, were associated with the pleiotropic locus QPtn.caas-4B for PTN. The SDRs may arise from chromosome recombination, gametophyte lethal genes, and segregation distortion factors [46, 47].

Materials and methods

Plant materials

Qingxinmai, a wheat landrace from Xinjiang, China, is characterized as long awn, slender spike, grain, and culm, and plenty of tillers (Fig. 1a). Line 041133 (pedigree: Jining 13/Tongmai 2) was developed in Qinghai province, with characteristics of tip-awns, thick spike, strong culm, fewer tillers, and larger grains (Fig. 1b). A RIL population consisting of 228 F2:9 lines was developed by a single seed decent method from cross Qingxinmai × 041133 to be used as the mapping population.

Phenotype assessments

During the 2019–2020 and 2020–2021 wheat growing seasons, field plots were set for assessing traits of the mapping population and the parents in the experimental farms of Institute of Crop Sciences, Chinese Academy of Agricultural Sciences in Beijing (2021BJ, 116.33°E, 39.96°N) and Changping, Beijing (2020CP, 2021CP, 116.26°E, 40.17°N), Zhaoxian, Hebei province (2020ZX, 114.78°E, 37.75°N), and Xinxiang, Henan province (2020XX, 113.98°E, 35.32°N), as well as a farm of Gansu Academy of Agricultural Sciences in Qingshui, Gansu province (2021QS, 105.80°E, 34.60°N). About 40 seeds of each line were planted in a one-row plot 2.0 m in length and a row spacing of 30 cm. A randomized complete block design with two replicates was used to arrange the RILs and the parents in each site. Field managements were performed according to the local practices for wheat production. Tiller traits, including TN at the tillering stage [Zadocks growth stage (GS) 31] [48] and PTN at the late milk stage (GS 77) were enumerated in 10 plants from each plot. At maturity, ten plants were randomly harvested from each plot to measure SL from the base of the rachis to the tip of terminal spikelet excluding awns and enumerate SNS. Length of awns at the top, central, and bottom of five spikes were measured. The phenotypic data of TGW, GL, and GW were measured by the Wanshen SC-G Automatic Seed Testing Analysis and Thousand Grain Weight Software (WSeen Inc., Hangzhou, China).

BSE-Seq analysis

Genomic DNA was isolated from grains with a cetyltrimethylammonium bromide method [49]. Bulked DNA pools (Bulk-HTGW and Bulk-LTGW) were constructed by separately mixing equal amounts of DNA samples from 40 high-TGW (36.5–49.9 g) and 40 low-TGW (20.4–29.0 g) RILs. These DNA bulks, together with the parents, were subjected to exome capture sequencing on the WheatPanExomeV2 platform at Chengdu Teuni Technology (Chengdu, China). Uncaptured DNA fragments were removed, and the enriched exons were amplified by PCR. High-throughput DNA sequencing was performed on the Illumina platform (Illumina Inc., San Diego, CA, USA). SNPs obtained by mutation detection were filtered with the criteria of allele frequency < 0.3 or > 0.7 using the SNP-index algorithm to determine the genotype frequency of the extreme bulks [50]. Significantly different SNP sites between the contrasting DNA bulks were statistically screened. Euclidean distance (ED) values of SNPs between the two DNA bulks were calculated using each allele depth with the quantile method. The ED values of SNP exceeding 99% was selected as the filtering threshold [51].

BSR-Seq analysis

Crown and inflorescence were sampled at GS 24 and GS 31 for RNA sequencing. Thirty phenotypically contrasting RILs each were chosen based on phenotypes of tiller and spike traits to construct high- and low-tiller number (Bulk-HTN and Bulk-LTN) and long- and short-spike (Bulk-LS and Bulk-SS) pools. A BSR-Seq analysis was performed following a previously described pipeline [52]. In brief, RNA purified from the bulked samples with an Illumina TruSeq RNA sample preparation kit was sequenced on an Illumina Hiseq 4000 platform (Illumina Inc., San Diego, CA, USA). Adapter and low-quality sequences were truncated using Trimmomatic v0.36 [53]. High-quality reads were aligned against the Chinese Spring (CS) reference genome sequence RefSeq v1.0 (http://wheat-urgi.versailles.infa.fr) with the aid of STARv2.5.1b software [54]. Single nucleotide polymorphism (SNP) variants [P < 1e-8 for the Fisher’s Exact Test (FET) and the allele frequency difference (AFD) > 0.6] were identified from confident alignments using the “Haplotype Caller” module assembled in software GATK v3.6 [55].

QTL mapping

The RILs and their parents were genotyped with genomic DNA samples by the wheat 16 K GBTS SNP array (MolBreeding Biotechnology Co. Ltd., Shijiazhuang, China) (http://www.molbreeding.com). Raw reads generated were processed with fastp v0.20.0 [56] and then aligned to the CS reference genome RefSeq v1.0 with Burrows-Wheeler aligner [57]. High-quality SNPs were obtained and filtered by GATK v3.5 [58], and SNP variants with the read depth < 5 were excluded from further analysis. SNP loci were classified as heterozygous genotypes when the SNP variation frequency ranged from 0.2 to 0.8, and the remaining genotypes were homozygous. Polymorphic SNPs between Qingxinmai and line 041133 were extracted for further analysis. A genetic linkage map was constructed for QTL calling using IciMapping 4.2 [59]. Only one marker was selected as a delegate from each bin to construct the linkage map. QTL for the traits were detected by IciMapping 4.2 with the inclusive composite interval mapping (ICIM) method. A test of 1,000 permutations was used to identify the logarithm of odds (LOD) threshold (> 3.0) that corresponded to a genome-wide false discovery rate of 5% (P < 0.05).

KASP marker development and QTL validation

DNA sequences of selected SNPs were used to develop KASP markers with the help of the Triticeae Multi-omics Center (http://202.194.139.32). Primers were designed using PolyMarker (http://polymarker.tgac.ac.uk/). The probe sequences for the FAM and the HEX signals were separately added to the primers specific for the two parental genotypes. A thermal Cycler (C1000 Touch™, Bio-Rad, Foster City, CA, USA) was used to perform KASP assays. The reaction mixture (10 µl) was prepared by mixing 5 µl of 2× master mix (Wuhan Gentides Biotech Co., Ltd., Wuhan, China), 0.2 µl of primer mix, 3 µl of ddH2O, and 2 µl of DNA template (50–150 ng/µl). The thermal cycling profile included 94 °C for 15 min hot-start activation, a touchdown phase of 10 cycles (94 °C for 20 s, touchdown at 61 °C initially and then decreased by 0.6 °C per cycle for 60 s), and 26 cycles of regular PCR (94 °C for 20 s, 55 °C for 1 min). The following cycling and resting steps set at 94 °C for 20 s and 57 °C for 60 s (3–10 cycles per step) were performed if signals were poorly clustered. End-point fluorescence data were screened using the microplate reader FLUOstar Omega SNP (BMG Labtech, Durham, NC, USA) and analyzed by the Klustering Caller Software (http://www.lgcgroup.com/). The KASP markers polymorphic between the two parents were used to genotype the entire RIL population.

Statistical analysis

The BLUE value was calculated with the Aov [analysis of variance (ANOVA) of multi-environmental trials] function in the QTL IciMapping 4.2 [59] to be used for combined QTL detection, correlation, and normal distribution analyses. The H2 for each trait were was also analyzed through QTL IciMapping 4.2. Phenotypic correlation was computed from the BLUE value of each line in SPSS v. 20.0 for Windows (IBM SPSS, Armonk, NY, USA). The Student’s t-test was used to evaluate the significance of differences in SPSS at P < 0.05 or P < 0.01.

Expression of candidate genes in the target genomic region

Information on gene annotation in target genomic intervals was obtained with JBrowse in the Triticeae Multi-Omics Center (http://202.194.139.32). Gene codes annotated were used to predict the gene expression levels in the GeneExpression toolbar in the Triticeae Multi-Omics Center (http://202.194.139.32). Significantly differential expression was defined as an absolute log2 value (fold change) > 1 at P < 0.01.

Total RNA was isolated from crowns (sampled at GS 31) and inflorescences (sampled at GS 41) of Qingxinmai and line 041133 by a FastPure Universal Plant Total RNA Isolation Kit (Vazyme Biotech Co., Ltd., Nanjing, China). The first-strand cDNA was synthesized using a PrimeScript RT Reagent Kit with gDNA Eraser (https://www.takarabiomed.com.cn/). Gene-specific quantitative real-time PCR (qPCR) primer pairs for candidate genes were designed according to the gene annotations from the CS reference genome RefSeq v1.0 [60]. qPCR assays were performed on a BioRad CFX system with the Taq Pro Universal SYBR qPCR Master Mix (Vazyme Biotech Co., Ltd., Nanjing, China). Wheat Actin gene was amplified as the reference gene. Relative expression was determined using the 2−ΔΔCT method [61]. Three biological replicates were taken from crowns and inflorescences. Three replicates for each RNA sample were run as technical replicates.

Data availability

The BSR-Seq and BSE-Seq data presented in this study are deposited in the National Genomics Data Center repository (https://ngdc.cncb.ac.cn/gsa/browse/CRA013836 and https://bigd.big.ac.cn/gsa/browse/CRA013931), accession number: CRA013836 and CRA013931.

Abbreviations

AFD:

The allele frequency difference

AL:

Awn length

Aov:

Analysis of variance (ANOVA) of multi-environmental trials

B1 :

Tipped1

B2 :

Tipped2

BJ:

Beijing

BLUE:

The best linear unbiased estimate

BSE-Seq:

Bulked segregant exome capture sequencing

BSR-Seq:

Bulked segregant RNA sequencing

Bulk-HTGW:

the high-TGW DNA pool

Bulk-HTN:

The high tiller numbers RNA pool

Bulk-LS:

The long-SL RNA pool

Bulk-LTGW:

The low-TGW DNA pool

Bulk-LTN:

The low tiller numbers RNA pool

Bulk-SS:

The short-SL RNA pool

CP:

Changping

CS:

The Chinese Spring

CTAB:

Cetyltrimethylammonium bromide

CV:

The coefficients of variation

ED:

Euclidean Distance

FET:

The Fisher’s Exact Test

GBTS:

Genotyping by target sequencing

GL:

Grain length

GNS:

Grain numbers per spike

GW:

Grain width

H 2 :

The broad-sense heritability

Hd :

Hooded

ICIM:

The inclusive composite interval mapping

KASP:

Kompetitive allele-specific PCR

LOD:

The logarithm of odds

PTN:

Productive tiller number

qPCR:

Quantitative real-time PCR

QS:

Qingshui

QTL:

Quantitative trait loci

RIL:

Recombinant inbred line, SDR: The segregation distortion regions

SL:

Spike length

SN:

Spikelet number

SNP:

Single nucleotide polymorphism

SNS:

Spikelet number per spike

TGW:

Thousand-grain weight

TN:

Tiller number

XX:

Xinxiang

ZX:

Zhaoxian

References

  1. FAO. Food Outlook: Biannual Report on Global Food Markets. Rome. 2018; pp. 67–74.

  2. Yang XJ. China’s rapid urbanization. Science. 2013;342:310.

    Article  CAS  PubMed  Google Scholar 

  3. Huang JK, Wei W, Cui Q, Xie W. The prospects for China’s food security and imports: will China starve the world via imports? J Integr Agric. 2017;16:2933–44.

    Article  Google Scholar 

  4. Qin XL, Liu FXC, Yu H, Bao BG, Tian SQ, Liao YC, et al. Wheat yield improvements in China: past trends and future directions. Field Crops Res. 2015;177:117–24.

    Article  Google Scholar 

  5. Li YX, Zhang WF, Ma L, Wu L, Shen JB, Davies WJ, et al. An analysis of China’s grain production: looking back and looking forward. Food Energy Sect. 2014;3:19–32.

    Article  Google Scholar 

  6. Yang L, Shao H, Wu QX, Yu J, Ran CF, Li LQ, et al. QTLs mapping and epistasis analysis for the number of tillers and spike number per plant in wheat. J Triticeae Crops. 2013;33:875–82.

    CAS  Google Scholar 

  7. Miao YP, Jing FL, Ma JF, Liu Y, Zhang PP, Chen T, et al. Major genomic regions for wheat grain weight as revealed by QTL linkage mapping and meta-analysis. Front Plant Sci. 2022;13:802310.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Xu DG, Wen WE, Fu LP, Li FJ, Li JH, Xie L, et al. Genetic dissection of a major QTL for kernel weight spanning the Rht–B1 locus in bread wheat. Theor Appl Genet. 2019;132:3191–200.

    Article  CAS  PubMed  Google Scholar 

  9. Song L, Liu J, Cao BL, Liu B, Zhang XP, Chen ZY, et al. Reducing brassinosteroid signalling enhances grain yield in semi-dwarf wheat. Nature. 2023;617:118–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Kumar A, Mantovani EE, Seetan R, Soltani A, Echeverry-Solarte M, Jain S, et al. Dissection of genetic factors underlying wheat Kernel shape and size in an Elite × Nonadapted cross using a high-density SNP linkage map. Plant Genome. 2016;9:0081.

    Article  Google Scholar 

  11. Wang XQ, Dong LH, Hu JM, Pang YL, Hu LQ, Xiao GL, et al. Dissecting genetic loci affecting grain morphological traits to improve grain weight via nested association mapping. Theor Appl Genet. 2019;132:3115–28.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Guan PF, Shen XY, Mu Q, Wang YF, Wang XB, Chen YM, et al. Dissection and validation of a QTL cluster linked to Rht-B1 locus controlling grain weight in common wheat (Triticum aestivum L.) using near-isogenic lines. Theor Appl Genet. 2020;133:2639–53.

    Article  CAS  PubMed  Google Scholar 

  13. Dholakia BB, Ammiraju JSS, Singh H, Lagu MD, Röder MS, Rao VS, et al. Molecular marker analysis of kernel size and shape in bread wheat. Plant Breed. 2003;122:392–5.

    Article  CAS  Google Scholar 

  14. Chen ZY, Cheng XJ, Chai LL, Wang ZH, Bian RL, Li J, et al. Dissection of genetic factors underlying grain size and fine mapping of QTgw.cau-7D in common wheat (Triticum aestivum L). Theor Appl Genet. 2020;133:149–62.

    Article  CAS  PubMed  Google Scholar 

  15. Kebrom TH, Spielmeyer W, Finnegan EJ. Grasses provide new insights into regulation of shoot branching. Trends Plant Sci. 2013;18:41–8.

    Article  CAS  PubMed  Google Scholar 

  16. Hussien A, Tavakol E, Horner DS, Muñoz-Amatriaín M, Muehlbauer GJ, Rossini L. Genetics of tillering in rice and barley. Plant Genome. 2014;7:1–20.

    Article  Google Scholar 

  17. Richards R. A tiller inhibitor gene in wheat and its effect on plant growth. Austr J Agric Res. 1998;39:749–57.

    Article  Google Scholar 

  18. Spielmeyer W, Richards R. Comparative mapping of wheat chromosome 1AS which contains the tiller inhibition gene (tin) with rice chromosome 5S. Theor Appl Genet. 2004;109:1303–10.

    Article  CAS  PubMed  Google Scholar 

  19. Ahmed HI, Heuberger M, Schoen A, Koo D, Quiroz-Chavez J, Adhikari L, et al. Einkorn genomics sheds light on history of the oldest domesticated wheat. Nature. 2023;620:830–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Zhao B, Wu TT, Ma SS, Jiang DJ, Bie XM, Sui N, et al. TaD27-B gene controls the tiller number in hexaploid wheat. Plant Biotechnol J. 2019;18:513–25.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Hu YS, Ren TH, Li Z, Tang YZ, Ren ZL, Yan BJ. Molecular mapping and genetic analysis of a QTL controlling spike formation rate and tiller number in wheat. Gene. 2017;634:15–21.

    Article  CAS  PubMed  Google Scholar 

  22. Liu JJ, Luo W, Qin NN, Ding PY, Zhang H, Yang CC, et al. A 55K SNP array-based genetic map and its utilization in QTL mapping for productive tiller number in common wheat. Theor Appl Genet. 2018;131:2439–50.

    Article  CAS  PubMed  Google Scholar 

  23. Liu JJ, Zhou JG, Tang HP, Tu Y, Mu Y, Gou LL, et al. A major vernalization-independent QTL for tiller angle on chromosome arm 2BL in bread wheat. Crop J. 2022;10:185–93.

    Article  Google Scholar 

  24. Li T, Deng GB, Su Y, Yang Z, Tang YY, Wang JH, et al. 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. Theor Appl Genet. 2021;134:3625–41.

    Article  CAS  PubMed  Google Scholar 

  25. Muqaddasi QH, Brassac J, Koppolu R, Plieske J, Ganal MW, Röder MS. TaAPO-A1, an ortholog of rice Aberrant Panicle Organization 1, is associated with total spikelet number per spike in elite European hexaploid winter wheat (Triticum aestivum L.) varieties. Sci Rep. 2019;9:13853.

  26. Ma J, Ding PY, Liu JJ, Li T, Zou YY, Habib A, et al. Identification and validation of a major and stably expressed QTL for spikelet number per spike in bread wheat. Theor Appl Genet. 2019;132:3155–67.

    Article  CAS  PubMed  Google Scholar 

  27. Saripalli G, Adhikari L, Amos C, Kibriya A, Ahmed HI, Heuberger M, et al. Integration of genetic and genomics resources in einkorn wheat enables precision mapping of important traits. Commun Biol. 2023;6:835.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Blum A. Photosynthesis and transpiration in leaves and ears of wheat and barley varieties. J Exp Bot. 1985;164:432–40.

    Article  Google Scholar 

  29. Sears ER. The aneuploids of common wheat. University of Missouri Agricultural Exp Stn Res Bull. 1954;572:1–58.

    Google Scholar 

  30. Li H, Han YC, Guo XX, Xue F, Wang CY, Ji WQ. Genetic effect of locus B2 inhibiting awning in double-ditelosomic 6B of Triticum durum DR147. Genet Resour Crop Evol. 2014;62:407–18.

    Article  Google Scholar 

  31. Takumi S, Kosugi T, Murai K, Mori N, Nakamura C. Molecular cloning of three homoeologous cDNAs encoding orthologs of the maize KNOTTED1 homeobox protein from young spikes of hexaploid wheat. Gene. 2000;249:171–81.

    Article  CAS  PubMed  Google Scholar 

  32. Tshikunde NM, Mashilo J, Shimelis H, Odindo A. Agronomic and physiological traits, and associated quantitative trait loci (QTL) affecting yield response in wheat (Triticum aestivum L.): a review. Front Plant Sci. 2019;10:1428.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Naruoka Y, Talbert LE, Lanning SP, Blake NK, Martin JM, Sherman JD. Identification of quantitative trait loci for productive tiller number and its relationship to agronomic traits in spring wheat. Theor Appl Genet. 2011;123:1043–53.

    Article  CAS  PubMed  Google Scholar 

  34. Liu G, Jia LJ, Lu LH, Qin DD, Zhang JP, Guan F, et al. Mapping QTLs of yield-related traits using RIL population derived from common wheat and tibetan semi-wild wheat. Theor Appl Genet. 2014;127:2415–32.

    Article  CAS  PubMed  Google Scholar 

  35. Heidari B, Sayed-Tabatabaei BE, Saeidi G, Kearsey M, Suenaga K. Mapping QTL for grain yield, yield components, and spike features in a doubled haploid population of bread wheat. Genome. 2011;54:517–27.

    Article  PubMed  Google Scholar 

  36. Li FJ, Wen WE, He ZH, Liu JD, Jin H, Cao SH, et al. Genome–wide linkage mapping of yield–related traits in three Chinese bread wheat populations using high-density SNP markers. Theor Appl Genet. 2018;131:1903–24.

    Article  PubMed  Google Scholar 

  37. Fan XL, Cui F, Ji J, Zhang W, Zhao XQ, Liu JJ, et al. Dissection of pleiotropic QTL regions controlling wheat spike characteristics under different nitrogen treatments using traditional and conditional QTL mapping. Front Plant Sci. 2019;10:187.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Liu SB, Zhou RG, Dong YC, Li P, Jia JZ. Development, utilization of introgression lines using a synthetic wheat as donor. Theor Appl Genet. 2006;112:1360–73.

    Article  CAS  PubMed  Google Scholar 

  39. Deng SM, Wu XR, Wu YY, Zhou RH, Wang HG, Jia JZ, et al. Characterization and precise mapping of a QTL increasing spike number with pleiotropic effects in wheat. Theor Appl Genet. 2011;122:281–9.

    Article  PubMed  Google Scholar 

  40. McIntosh RA, Hart GE, Devos KM, Rogers J, Gale MD. Catalogue of gene symbols for wheat. Wheat Inf Ser. 1998;86:54–91.

    Google Scholar 

  41. Kosuge K, Watanabe N, Kuboyama T, Melnik VM, Yanchenko VI, Rosova MA, et al. Cytological and microsatellite mapping of mutant genes for spherical grain and compact spikes in durum wheat. Euphytica. 2008;159:289–96.

    Article  CAS  Google Scholar 

  42. DeWitt N, Guedira M, Lauer E, Sarinelli M, Tyagi P, Fu DL, et al. Sequence based mapping identifies a candidate transcription repressor underlying awn suppression at the B1 locus in wheat. New Phytol. 2020;225:326–39.

    Article  CAS  PubMed  Google Scholar 

  43. Huang DQ, Zheng Q, Melchkart T, Bekkaoui Y, Konkin DJF, Kagale S, et al. Dominant inhibition of awn development by a putative zinc-finger transcriptional repressor expressed at the B1 locus in wheat. New Phytol. 2020;225:340–55.

    Article  CAS  PubMed  Google Scholar 

  44. Ke WS, Xing JW, Chen ZY, Zhao YD, Xu WY, Tian LL, et al. The TaTCP4/10–B1 cascade regulates awn elongation in wheat (Triticum aestivum L). Plant Commun. 2023;4:100590.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Paillard S, Schnurbusch T, Winzeler M, Messmer M, Sourdille P, Abderhalden O, et al. An integrative genetic linkage map of winter wheat (Triticum aestivum L). Theor Appl Genet. 2003;107:1235–42.

    Article  CAS  PubMed  Google Scholar 

  46. Nasuda S, Friebe B, Gill BS. Gametocidal genes induce chromosome breakage in the interphase prior to the first mitotic cell division of the male gametophyte in wheat. Genetics. 1998;149:1115–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Marais GF, Marais AS, Groenewald JZ. Evalyation and reduction of Lr19-49, a recombined from of the Lr19 translocation of wheat. Euphytica. 2001;121:289–95.

    Article  CAS  Google Scholar 

  48. Zadoks JC, Chang TT, Konzak CF. A decimal code for the growth stages of cereals. Weed Res. 1974;14:415–21.

    Article  Google Scholar 

  49. Chatterjee A, Moulik S, Majhi P, Sanyal S. Studies on surfactant biopolymer interaction. I. Microcalorimetric investigation on the interaction of cetyltrimethylammonium bromide (CTAB) and sodium dodecylsulfate (SDS) with gelatin (gn), lysozyme (lz) and deoxyribonucleic acid (DNA). Biophys Chem. 2002;98:313–27.

    Article  CAS  PubMed  Google Scholar 

  50. Wang QG, Xiong HC, Guo HJ, Zhao LS, Xie YD, Gu JY, et al. Genetic analysis and mapping of dwarf gene without yield penalty in a γ-ray-induced wheat mutant. Front Plant Sci. 2023;14:1133024.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Hill J, Demarest B, Bisgrove B, Gorsi B, Su Y, Yost H. MMAPPR: mutation mapping analysis pipeline for pooled RNA-seq. Genome Res. 2013;23:687–97.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Xie JZ, Guo GH, Wang Y, Hu TZ, Wang LL, Li JT, et al. A rare single nucleotide variant in Pm5e confers powdery mildew resistance in common wheat. New Phytol. 2020;228:1011–26.

    Article  CAS  PubMed  Google Scholar 

  53. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29:15–21.

    Article  CAS  PubMed  Google Scholar 

  55. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The genome analysis Toolkit: a MapReduce framework for analyzing nextgeneration DNA sequencing data. Genome Res. 2010;20:1297–303.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Chen SF, Zhou YQ, Chen YR, Gu J. Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018;34:i884–90.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25:1754–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. DePristo MA, Banks E, Popli R, Garimella KV, Maguire JR, Hartl C, et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011;43:491–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Meng L, Li HH, Zhang LY, Wang JK. QTL IciMapping: integrated software for genetic linkage map construction and quantitative trait locus mapping in biparental populations. Crop J. 2015;3:269–83.

    Article  Google Scholar 

  60. International Wheat Genome Sequencing Consortium (IWGSC). Shifting the limits in wheat research and breeding using a fully annotated reference genome. Science. 2018;361:eaar7191.

    Article  Google Scholar 

  61. Livak JK, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2–∆∆CT method. Methods. 2001;25:402–8.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The financial support from the National Key Research and Development Program of China (2023YFD1200403) is gratefully acknowledged.

Funding

The financial support from the National Key Research and Development Program of China (2023YFD1200403) is gratefully acknowledged.

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization, H.J.L. Z.Y.L., and Y.J.Z; methodology, Y.H.L., J.H.H., D.Q. and Q.H.W.; software, Y.H.L., Q.H.W., Y.F.Q., and J.H.H; data curation, Y.H.L., J.H.H., J.Y.D, D.Q., Y.F.Q., H.L.L, L.M., L.H., Y.Z., H.J.Z., L.Y., and H.W.L.; writing-review and editing, H.J.L., Y.H.L., Z.Y.L., and Y.J.Z.; funding acquisition and supervision, H.J.L., Z.Y.L, and Y.J.Z. All authors have read and agreed to the published version of the manuscript.

Corresponding authors

Correspondence to Zhiyong Liu, Yijun Zhou or Hongjie Li.

Ethics declarations

Ethics approval and consent to participate

All experimental research methods in this study were performed in accordance with the relevant guidelines and regulations.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, Y., Hu, J., Qu, Y. et al. Alleles on locus chromosome 4B from different parents confer tiller number and the yield-associated traits in wheat. BMC Plant Biol 24, 454 (2024). https://doi.org/10.1186/s12870-024-05079-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-024-05079-4

Keywords