Skip to main content

Rapid identification and characterization of genetic loci for defective kernel in bread wheat



Wheat is a momentous crop and feeds billions of people in the world. The improvement of wheat yield is very important to ensure world food security. Normal development of grain is the essential guarantee for wheat yield formation. The genetic study of grain phenotype and identification of key genes for grain filling are of great significance upon dissecting the molecular mechanism of wheat grain morphogenesis and yield potential.


Here we identified a pair of defective kernel (Dek) isogenic lines, BL31 and BL33, with plump and shrunken mature grains, respectively, and constructed a genetic population from the BL31/BL33 cross. Ten chromosomes had higher frequency of polymorphic single nucleotide polymorphism (SNP) markers between BL31 and BL33 using Wheat660K chip. Totally 783 simple sequence repeat (SSR) markers were chosen from the above chromosomes and 15 of these were integrated into two linkage groups using the genetic population. Genetic mapping identified three QTL, QDek.caas-3BS.1, QDek.caas-3BS.2 and QDek.caas-4AL, explaining 14.78–18.17%, 16.61–21.83% and 19.08–28.19% of phenotypic variances, respectively. Additionally, five polymorphic SNPs from Wheat660K were successfully converted into cleaved amplified polymorphic sequence (CAPS) markers and enriched the target regions of the above QTL. Biochemical analyses revealed that BL33 has significantly higher grain sucrose contents at filling stages and lower mature grain starch contents than BL31, indicating that the Dek QTL may be involved in carbohydrate metabolism. As such, the candidate genes for each QTL were predicated according to International Wheat Genome Sequence Consortium (IWGSC) RefSeq v1.0.


Three major QTL for Dek were identified and their causal genes were predicted, laying a foundation to conduct fine mapping and dissect the regulatory mechanism underlying Dek trait in wheat.


Wheat (Triticum aestivum L.) supplies staple food for more than 35% of population worldwide ( The increase of population, decline of cultivated land and frequent occurrence of extreme climate posed large challenges to meet the consumption of wheat. Identification and mining of genetic loci for grain yield is important to improve the yield potential in wheat. Yield is determined by multiple components, such as grain weight and size. Grain plumpness is an important parameter affecting grain weight and morphological quality [1]. Many genetic loci controlling grain weight have been identified in previous studies [2], whereas genetic analyses on grain plumpness associated with grain filling is rarely reported in wheat. The defective kernel (Dek) mutants in gramineous crops are directly reflected in grain plumpness and usually characterized by shriveled kernel and low grain weight. Thus, identification of genetic loci for Dek is important to improve yield potential and appearance quality. Quite a few Dek mutants were identified in maize [3]. Of them, the genes of dek1, dek35, dek38 and dek39 have been isolated, and all involved in the synthesis and accumulation of starches and/or proteins in the endosperm [4,5,6,7]. Many Dek mutants have been generated in wheat, but their genetic information remains largely unknown. SSR markers are characterized with high polymorphisms, good repeatability, co-dominance, and distribution throughout whole genome, and favored for QTL mapping and molecular breeding in wheat [8,9,10]. With rapid advancements in sequencing technology, high-throughput sequencing or chip-based SNP genotyping platforms have been used for QTL mapping [11,12,13,14]. Additionally, the wheat reference genome was assembled and released recently [15]. The availability of high-quality wheat genome data and high-throughput SNP genotyping systems is very important to improve wheat genetics and breeding. We identified a Dek line in a spontaneous variation population and found that its kernels are severely shrunken and sunken. The objective of this study was to rapidly achieve QTL mapping and causal gene prediction for Dek through multiple genotyping systems, wheat reference genome and biochemical analysis.


Phenotypic analysis

The kernels of BL31 were normal and plump, whereas the Dek rate of BL33 is 100% in all environments (Additional file 1). The Dek rates in F2 and F2:3 populations were highly correlated among three environments with correlation coefficients ranging from 0.70 to 0.92 (Table 1). Variance analyses showed that both genotypes (P < 0.001) and environments (P < 0.05) had significant effects on the target trait (Table 2). The genotype × environment interactions did not reach a statistical significance (P > 0.05) (Table 2) and the broad-sense heritability of Dek rate in F2:3 across two environments was 0.95, denoting that the Dek rate was mainly determined by genotypes. The Dek rates in the populations showed continuous distribution (Fig. 1, a-d), indicating that Dek is controlled by multiple genetic loci. Additionally, the frequency graph underpinning the phenotypes is abnormal distribution with the “valley - peak” pattern, suggesting the presence of major QTL controlling Dek in the populations [16].

Table 1 Statistical analysis of Dek rates for F2 and F2:3 populations across three environments
Table 2 Analysis of variance of Dek rates for F2:3 lines derived from BL33/BL31
Fig. 1
figure 1

The frequency distributions of Dek rates in Gaoyi 2016–2017 cropping season (a), Gaoyi 2017–2018 cropping season (b), Xinxiang 2017–2018 cropping season (c), and the average data across three environments (d)

Chromosomal locations of Dek loci and construction of linkage maps

Theoretically, the more polymorphic markers between isogenic lines one chromosome contains, the higher probability that it harbors QTL for target traits. To rapidly identify the chromosomal locations of Dek QTL, we genotyped the isogenic lines, BL31 and BL33, with contrasting grain phenotypes using the Wheat660K SNP array. Ten chromosomes, 1A, 1B, 2A, 3A, 3B, 4A, 4B, 5A, 6B and 7B, were considered to be potential target chromosomes with Dek QTL based on the frequency distributions of polymorphic SNPs (Additional file 2). To further verify the chromosomal location of Dek QTL, SSR markers from the above target chromosomes were selected to construct genetic map. Totally 783 SSR markers on these chromosomes were selected to construct linkage map for QTL analysis (Additional file 3) [8, 17,18,19,20,21,22,23]. Of them, 118 SSR markers displayed polymorphisms between two parents, while only 35 SSRs in chromosomes 3B (13), 4A (19) and 5A (3), showed consistent polymorphisms between two groups of lines with contrasting phenotypes, i.e. 40 defective kernel and 40 normal kernel lines, respectively. Subsequently, the polymorphic SSR markers were used to genotype the entire F2 population, and two linkage groups were constructed in chromosomes 3B and 4A, including 6 and 9 SSR markers, respectively (Fig. 2).

Fig. 2
figure 2

LOD contours for Dek QTL. The genetic regions of QDek.caas-3BS.1, QDek.caas-3BS.2 and QDek.caas-4AL are highlighted in green, blue and yellow on chromosomes. LOD contours in Gaoyi 2017, Xinxiang 2018, Gaoyi 2018, and the average are indicated with red, green, wathet, and deongaree lines, respectively. LOD threshold of 2.5 is indicated by the dotted vertical line

QTL mapping

Based on the resultant linkage groups above and Dek rates in the genetic populations, we performed QTL mapping for Dek traits; three QTL were detected on chromosomes 3B and 4A, designated as QDek.caas-3BS.1, QDek.caas-3BS.2 and QDek.caas-4AL, respectively (Table 3; Fig. 2). QDek.caas-3BS.1 was flanked by Xbarc133 and Xcfd79 spanning a genetic interval of 6.77 cM; QDek.caas-3BS.2 was located between Xcfd79 and Xwmc808 in an interval of 35.03 cM; QDek.caas-4AL was mapped between Xwmc500 and Xgpw3079 in a genetic interval of 15.17 cM (Fig. 2). Based on the preliminary mapping above, the physical regions of QDek.caas-3BS.1, QDek.caas-3BS.2 and QDek.caas-4AL can be defined by their flanking markers according to IWGSC reference genome RefSeq v1.0 ( [15]. The physical intervals of QDek.caas-3BS.1, QDek.caas-3BS.2 and QDek.caas-4AL are 7.60–18.84 Mb, 18.84–32.75 Mb in chromosomes 3BS and 704.34–710.05 Mb in chromosome 4AL. Genetic analyses showed that QDek.caas-3BS.1, QDek.caas-3BS.2 and QDek.caas-4AL could explain 14.78 to 18.17%, 16.61 to 21.83% and 19.08 to 28.19% of the phenotypic variances, respectively (Table 3). Each of these QTL determined more than 10% phenotypic variation of Dek in each environment, indicating that they are major QTL, consistent with the “peak - valley” distribution pattern of Dek rate in genetic populations. Additionally, all the alleles governing the formation of Dek came from BL33.

Table 3 QTL for Dek trait detected in the F2 and F2:3 populations across different environments

Marker enrichment in the target regions of QTL

To further narrow down the genetic intervals of target QTL, the polymorphic SNPs identified from the Wheat660K chips were selected according to the above physical regions defined by QTL preliminary mapping. Totally 60 polymorphic SNPs were used to develop site-specific CAPS markers and five of these, AX-109027972, AX-109516011, AX-110042240, AX-109301653 and AX-108996126, were successfully mapped in the linkage groups (Additional file 4; Additional file 5; Additional file 6). Subsequently, QDek.caas-3BS.2 and QDek.caas-4AL were narrowed to the genetic intervals of 14.01 cM and 13.18 cM, respectively (Table 3; Fig. 2). Accordingly, the physical intervals of QDek.caas-3BS.2 and QDek.caas-4AL are reduced to 1.16 Mb and 1.13 Mb, respectively (Table 3). Although the polymorphic loci AX-108996126 and AX-110042240 did not further diminish the genetic interval of QDek.caas-3BS.1 and QDek.caas-4AL, they saturated the target regions (Fig. 2). Dek trait displays abnormal development of grain and thus it is usually considered to affect grain weight. To verify this case, we analyzed the effect of the resultant Dek QTL on thousand grain weight (TGW). Compared with BL31, the TGW of BL33 significantly decreased 24.12–28.81% (Additional file 7). QTL analyses based on the linkage maps above showed that only one QTL for TGW, designated as QTGW.caas-3BS, was detected between Xcfd79 and AX-109027972, explaining 19.18 to 23.94% of phenotypic variances (Additional file 8). QTGW.caas-3BS overlapped QDek.caas-3BS.2 on the linkage map, suggesting that the Dek QTL QDek.caas-3BS.2 may affect TGW. No significant effect on TGW was detected in the genetic intervals of the other two QTL, QDek.caas-3BS.1 and QDek.caas-4AL.

Predication of candidate genes for the Dek QTL

The cloned dek genes in maize are involved in the grain growth and development, especially the synthesis and storage of starches and/or proteins in the endosperm [4,5,6,7]. Our biochemical analysis also showed that the mature grain starch contents of BL33 were significantly lower than those of BL31 (Additional file 9). In contrast, BL33 had higher sucrose contents and slower sucrose reduction across the filling stages than BL31 (Additional file 10). Accordingly, we assumed that the causal genes for the three QTL are involved in carbohydrate metabolism. In this study, three major QTL was mapped into the genetic intervals of less than 15 cM. The physical region of QDek.caas-3BS.1 is 11.24 Mb, whereas those of QDek.caas-3BS.2 and QDeg.caas-4AL are only 1.16 and 1.13 Mb, respectively, according to IWGSC RefSeq v1.0 (Table 3). Totally 233 genes were present in the physical intervals of the three QTL (Additional file 11). We found nine genes (Gene ID: TraesCS3B01G025200, TraesCS3B01G028100, TraesCS3B01G028200, TraesCS3B01G028300, TraesCS3B01G028500, TraesCS3B01G028700, TraesCS3B01G039100, TraesCS3B01G039800, and TraesCS4A01G446700) involving carbohydrate metabolism or grain development in or near the target regions of the Dek QTLs (Additional file 11). To further identify candidate genes for the Dek QTL, the spatio-temporal expression patterns of the nine genes were analyzed using the WheatExp ( Of these, four genes, TraesCS3B01G025200, TraesCS3B01G028700, TraesCS3B01G039800, and TraesCS4A01G446700, had an observable expression during grain development (Additional file 12) [24]. Based on functional annotation in IWGSC RefSeq v1.0, the four genes encoded fructose-bisphosphate aldolase (Fba), β-fructofuranosidase, abscisic acid-deficient 4 (ABA4), and sucrose synthase (Sus), respectively, tentatively designated as TaFba-3B, TaBff-3B, TaABA4-3B, and TaSus-4A. Their expression patterns were confirmed in the immature grains harvested at three different stages (5, 15, and 25 DAF) using qPCR (Additional file 13). The results showed that all of the above genes expressed in grain and had the lowest expression level at 5 DAF. TaBff-3B and TaSus-4A was mainly detected at 15 DAF, while TaFba-3B and TaABA4-3B had the highest expression activity in 25-DAF grains.


A rapid and cost-effective QTL mapping strategy through multiple genotyping systems in combination with wheat reference genome

With rapid advances in sequencing technology, more and more wheat varieties were sequenced and a large number of SNPs have been identified [25,26,27]. Consequently, a few sequencing-based and chip-based SNP genotyping platforms have been developed and used for genetic analyses in wheat [28,29,30]. Some studies have achieved fine mapping of QTL using SNP chips [31,32,33,34]. The chip-based SNP arrays are high-throughput and effective for genetic analyses, but the chips are very expensive. In this study, to rapidly and cost-effectively map QTL for Dek, we only genotyped the isogenic parental lines with contrasting phenotypes using the Wheat660K chip. Firstly, the target chromosomes potentially harboring the QTL for Dek were determined according to the frequency distributions of the polymorphic SNPs. Secondly, it is necessary to construct linkage map for QTL analysis through genotyping genetic populations. However, it is well-known that successful conversion of SNP to the flexible markers, such as KASP and CAPS, is laborious and time-consuming, especially in the hexaploid wheat with three homoeologous subgenomes. Considering high polymorphisms, good repeatability, co-dominance, and distribution throughout whole genome of SSR markers, we harvested the handy SSR markers in the candidate chromosomes above and rapidly constructed linkage maps [8, 35]. Thirdly, we performed genetic mapping for Dek QTL and defined their physical intervals according to IWGSC RefSeq v1.0. Finally, we successfully converted some of the polymorphic SNPs in the physical target regions to CAPS markers and narrowed down the genetic regions of the QTL. In all, we set up a rapid and efficient strategy for QTL mapping through combining high-throughput SNP arrays with cost-effective and flexible SSR markers as well as wheat reference genome (Fig. 3). Of course, a genome-wide linkage map is necessary to find the QTL in all chromosomes. In the future, we intend to construct a recombinant inbred line (RIL) population and achieve genome-wide identification for Dek QTL.

Fig. 3
figure 3

The workflow of rapid and cost-effective QTL mapping for Dek trait. The production of our Fig. 3 referred to the frame of Fig. 1 in Ramirez-Gonzalez et al. [36]

Comparisons of the Dek QTL with previously reported Dek-related genetic loci

The comparisons of genetic loci for Dek-related traits can not only improve the understanding of genetic effects of the resultant QTL in the present study but also supply with genetic information for QTL fine mapping. QDek.caas-3BS.1 is linked with Xbarc133 (Table 3; Fig. 2). Two QTL, QGfrmax.nfcri-3B and QGfd.nfcri-3B, controlling grain filling rate (GFR) and grain filling duration (GFD), respectively, were linked to Xgwm533 with a physical distance of 1 Mb from Xbarc133 [37]. QTgw.nfcri-3B for TGW was also linked to Xgwm533 [37]. QDek.caas-3BS.1 and QDek.caas-3BS.2 were identified in adjacent regions and both are linked to Xcfd79 (Table 3; Fig. 2). A QTL for TGW linked with Xcfd79 was detected by Groos et al. [38]. Additionally, QMswscf.acs-3B.1 for water soluble carbohydrates (WSC) of main stem at the early anthesis stage, were linked to Xcfd79 [39]. QDek.caas-4AL is flanked by Xwmc500 and AX-109301653 (Table 3). The physical distance between QDek.caas-4AL and Wx-B1 is 16 Mb approximately. So far, many QTL for the contents of starch, protein and amylose and amylopectin, and grain yield have been reported on chromosome 4AL close to Wx-B1 [40,41,42,43,44,45,46,47]. QGsc4A, a QTL for grain starch contents, is closely linked to Wx-B1 [42]. Qftsc4A for total grain starch content and Qfams4A for grain amylose content were also linked to Wx-B1 [46]. To identify the genetic relationship between QDek.caas-4AL and Wx-B1, the gene-specific functional marker developed by Nakamura et al. [48] was used to genotype the parents, BL31 and BL33. The result showed that no polymorphism was detected between them (Additional file 14), suggesting that Wx-B1 had little effect on Dek trait.

Candidate genes analysis

Two candidate genes, TaFba-3B and TaBff-3B, were detected in QDek.caas-3BS.1 region, and encoded Fba and β-fructofuranosidases, respectively (Table 4). Fba is one of the nonregulated enzymes comprising the Calvin cycle and controlling carbon fixation and starch synthesis in photosynthetic organs, which can catalyze reversible aldol condensation reaction of dihydroxy acetone phosphate and 3-phosphate glyceraldehyde yielding fructose-1,6-diphosphate [49, 50]. β-fructofuranosidase is a key enzyme for synthesis, metabolism and transportation of sucrose and starch, and can irreversibly catalyze sucrose hydrolysis to glucose and fructose [51, 52]. TaFba-3B and TaBff-3B probably regulate carbohydrate metabolism and consequently affect Dek phenotype. TaABA4-3B was predicted as a candidate gene for QDek.caas-3BS.2. ABA4 functions in neoxanthin synthesis [53]. Neoxanthin is necessary for maintaining the integrity of the light harvesting chlorophyll and photosynthetic function of leaves, involving the metabolism of photosynthetic product [54, 55]. TaSus-4A may be the causal gene of QDek.caas-4AL. Sus deficiency led to the shrunken-1 mutations in maize [56]. Association analyses showed that the two sucrose synthase genes, TaSus1 and TaSus2 on wheat homoeologous chromosome groups 7 and 2, respectively, had significant effects on TGW [57, 58]. In near future, we will identify polymorphic sites of these candidate genes to conduct fine mapping for the target QTL and further determine the causal genes for Dek through wheat transformation experiments.

Table 4 Candidate genes of Dek QTLs


In the study, we detected three DEK QTLs, QDek.caas-3BS.1, QDek.caas-3BS.2 and QDek.caas-4AL, each of which could explain more than 10% phenotypic variation. Additionally, their candidate genes, TaFba-3B, TaBff-3B, TaABA4-3B, and TaSus-4A, were predicated according to biochemical analysis and IWGSC RefSeq v1.0 annotation. These findings provide important information for QTL fine mapping and molecular mechanistic dissection underlying Dek in wheat.


Plant materials

A number of Dek lines were generated in Wheat Research Institute, Gansu Academy of Agricultural Science (GAAS). Through multi-generation self-crossing, we identified a pair of isogenic lines, BL31 and BL33, with plump and shrunken mature grains, respectively. Genetic populations comprising 146 F2 plants and their corresponding F2:3 lines from the cross between BL31 and BL33 were created for QTL mapping.

Field trials and phenotype evaluation

The F2 population and parents were planted at Gaoyi (37°33′N, 114°26′E) in Hebei province during the 2016–2017 cropping season. The F2:3 lines were sown at Xinxiang (34°53′N, 113°23′E) in Henan province and Gaoyi during the 2017–2018 cropping season. Field trials were arranged in randomized complete blocks with two replicates at each location. Each plot comprised a single 2 m row with 20 plants and spaced 30 cm apart. Grains were harvested in individual plants or lines. Grain phenotype of F2 plants was scored using the percentages of defective kernels (defined as Dek rate) in single plants. For F2:3 lines, ten representative plants in each line were selected from the middle part of each plot to figure out Dek rate, then the average Dek rate of each plot was used for subsequent statistical and genetic analysis.

Preliminary mapping of QTL on target chromosomes using SNP chip

High-quality genomic DNA was extracted from fresh leaves following Guillemaut and Laurence [59]. The Wheat660K SNP array containing 630,517 SNPs, a high-throughput genotyping tool developed by Institute of Crop Sciences, Chinese Academy of Agricultural Sciences (, was used to test two isogenic lines BL31 and BL33 at CapitalBio Corporation (Beijing, China; The frequency distribution of polymorphic SNPs in each chromosome was analyzed using Microsoft Excel 2016 ( The target chromosomes with potential QTL were determined based on the frequency distribution of polymorphic SNPs in each chromosome.

Genotyping by SSR markers

Polymerase chain reaction (PCR) for SSR markers was performed in a reaction mixture of 15 μl containing 7.5 μl of 2 × Taq PCR Mix (Beijing HT-biotech Co., Ltd., Beijing;, 5 pmol of each primer and 50–100 ng of genomic DNA using the ABI Applied Biosystems Veriti 96 Well Thermal Cycler (Gene Co., Ltd., Shanghai; PCR amplification was performed at 95 °C for 5 min, followed by 35 cycles (94 °C for 30s, 55–60 °C for 30s and 72 °C for 30s) and a final extension at 72 °C for 7 min. The PCR products were separated in 8% polyacrylamide gels at 250 V using the JY 600HE Electrophoretic Apparatus (Beijing Junyi Electrophoresis Co., Ltd., Beijing; or 2% agarose gels at 220 V using the DYY-7C Electrophoresis Apparatus (Beijing Liuyi Biotechnology Co., Ltd.;

Converting chip-based polymorphic SNPs into CAPS markers

To enrich the target genetic region, polymorphic SNPs between BL31 and BL33, are needed to convert into flexible markers, such as kompetitive allele specific PCR (KASP) and CAPS markers. CAPS technology is an efficient tool to genotype SNPs [60], combining PCR with restriction enzymes digestion to differentiate the target alleles with polymorphic SNPs. The polymorphic SNPs between parents in the target region were identified by alignments with DNAMAN software ( for designing chromosome-specific CAPS markers [61]. The primers and restriction enzymes used for CAPS markers were listed in Additional file 4.

Statistical analysis

The analysis of variance (ANOVA) was conducted using the SAS 9.2 software ( with PROC MIXED procedure. Broad-sense heritability (hB2) was calculated using the following formula [62]: hB2 = σ2g / (σ2g + σ2ge / r + σ2ε / re).

in which, σ2g, σ2ge and σ2ε were estimates of genotypes, genotype × environment interaction and residual error variances, respectively; r and e were replicates per environment and the number of environments, respectively. The correlation coefficients of traits among three different environments were computed using the PROC CORR procedure with the SAS 9.2 software.

Linkage map construction and QTL analysis

Joinmap 4.0 software ( was used for linkage map construction [63] and genetic distances between markers were estimated based on the Kosambi mapping function [64]. QTL analysis was conducted by inclusive composite interval mapping (ICIM) based on the Dek rate using IciMapping 4.0 software ( [65, 66]. The logarithm of odds (LOD) threshold score of 2.5 was determined for declaring significant QTL based on 1000 permutations with a type I error of P < 0.05. The phenotypic variance explained (PVE) were calculated through stepwise regression to estimate individual QTL effects.

Assays for grain starch and sucrose contents

Six mature grains were grinded to fine powder in a mortar. The resultant powder (0.1 g) was used to measure starch contents with Starch Assay Kit (Beijing Solarbio Science & Technology Co., Ltd., Beijing, Then, the following formula was used to convert mg/g to mg/grain for scoring the grain starch content: Grain starch content (mg/grain) = starch contents (mg/g) * TGW / 1000.

TGW data was shown in Additional file 7. Six grains of two parents were collected at 5, 10, 15, 20, 25, and 30 days after flowering (DAF) to validate the dynamic trend of grain sucrose contents during the filling stage. The dried grain were grinded to fine powder in mortar, and 0.1 g powder was used to measure grain sucrose contents using Plant Sucrose Assay kit (Beijing Solarbio Science & Technology Co., Ltd., Beijing,

Quantitative real-time PCR

Total RNA was extracted from fresh grain 5, 15, and 25 days after flowering (DAF) using RNAprep pure Plant Kit (Tiangen Biotech Co., Ltd., Reverse transcription was performed using the reverse transcriptase M-MLV (RNase H) (Takara Biomedical Technology Co., Ltd., The iTaq™ Universal SYBR Green Supermix on CFX Real-Time System (Bio-Rad) was used to perform quantitative real-time PCR (qPCR). The profile of qPCR was as follows: 95 °C for 2 min, followed by 40 cycles of 95 °C for 5 s and 58 °C for 20 s. A final dissociation stage was run to generate a melting curve for judgement of amplification specificity. Quantification of qPCR was determined using the 2-ΔΔCt method [67]. Elongation factor 1 α (EF1α) gene was used as the internal control to calibrate the expression level of genes of interest. An additional file shows the primers used in qPCR (see Additional file 15).

Availability of data and materials

The data sets supporting the results of this article are included within the article and its additional files.



adenosine diphosphate glucose pyrophosphorylase


analysis of variance


cleaved amplified polymorphic sequence


cell wall invertase


days after flowering


defective kernel


fructose-bisphosphate aldolase


grain filling duration


grain filling rate


genome-wide association study


inclusive composite interval mapping


International Wheat Genome Sequence Consortium


kompetitive allele specific polymerase chain reaction


logarithm of odds


phenotype variation explained


quantitative real-time PCR


single nucleotide polymorphism


simple sequence repeat


sucrose synthase


thousand grain weight


water soluble carbohydrates


  1. Cheng SH, Yang SM, Zhang BQ, Ji KZ, Zhao BH, Gao DR. Kernel plumpness and filling index in wheat (Triticum aestivum L.). Jiangsu J Agric Sci. 1993;9:7–10.

    Google Scholar 

  2. 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 

  3. Neuffer MG, Sheridan WF. Defective kernel mutants of maize. I. Genetic and lethality studies. Genetics. 1980;95:929–44.

    CAS  PubMed  PubMed Central  Google Scholar 

  4. Lid SE, Gruis D, Jung R, Lorentzen JA, Ananiev E, Chamberlin M, et al. The defective kernel 1 (dek1) gene required for aleurone cell development in the endosperm of maize grains encodes a membrane protein of the calpain gene superfamily. Proc Natl Acad Sci. 2002;99:5460–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Chen XZ, Feng F, Qi WW, Xu LM, Yao DS, Wang Q, et al. Dek35 encodes a PPR protein that affects cis-splicing of mitochondrial nad4 intron 1 and seed development in maize. Mol Plant. 2017;10:427–41.

    Article  CAS  PubMed  Google Scholar 

  6. Garcia N, Li YB, Dooner HK, Messing J. Maize defective kernel mutant generated by insertion of a ds element in a gene encoding a highly conserved TTI2 cochaperone. Proc Natl Acad Sci. 2017;114:5165–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Li XJ, Wei G, Sun SL, Chen ZL, Jing C, Song W, et al. Defective Kernel 39 encodes a PPR protein required for seed development in maize. J Integr Plant Biol. 2018;60:45–64.

    Article  CAS  PubMed  Google Scholar 

  8. Röder MS, Korzun V, Wendehake K, Plaschke J, Tixier MH, Leroy P, et al. A microsatellite map of wheat. Genetics. 1998;149:2007–23.

    PubMed  PubMed Central  Google Scholar 

  9. Wang Y, Yang C, Jin QJ, Zhou DJ, Wang SS, Yu YJ, et al. Genome-wide distribution comparative and composition analysis of the SSRs in Poaceae. BMC Genet. 2015;16:18.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. UPOV-BMT, BMT/36/10 Progress report of the 36th session of the technical committee, the technical working parties and working group on biochemical and molecular techniques and DNA-profiling in particular, Geneva. 2002.

  11. Sukumaran S, Dreisigacker S, Lopes M, Chavez P, Reynolds MP. Genome-wide association study for grain yield and related traits in an elite spring wheat population grown in temperate irrigated environments. Theor Appl Genet. 2015;128:353–63.

    Article  CAS  PubMed  Google Scholar 

  12. Gao LL, Turner MK, Chao S, Kolmer J, Anderson JA. Genome wide association study of seedling and adult plant leaf rust resistance in elite spring wheat breeding lines. PLoS One. 2016;11:e0148671.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  13. Kassa MT, You FM, Hiebert CW, Pozniak CJ, Fobert PR, Sharpe AG. Highly predictive SNP markers for efficient selection of the wheat leaf rust resistance gene Lr16. BMC Plant Biol. 2017;17:45.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  14. Mourad AMI, Sallam A, Belamkar V, Wegulo S, Bowden R, Jin Y. Genome-wide association study for identification and validation of novel SNP markers for Sr6 stem rust resistance gene in bread wheat. Front Plant Sci. 2018;9:380.

    Article  PubMed  PubMed Central  Google Scholar 

  15. IWGSC, Appels R, Eversole K, Feuillet C, Keller B, Rogers J, et al. Shifting the limits in wheat research and breeding using a fully annotated reference genome. Science. 2018;361.

  16. Cao SH, Guo XL, Liu DC, Zhang XQ, Zhang AM. Preliminary gene-mapping of photoperiod-temperature sensitive genic male sterility in wheat (Triticum aestivum L.). Acta Genet Sin. 2004;31:293–8.

    CAS  PubMed  Google Scholar 

  17. Song QJ, Fickus EW, Cregan PB. Characterization of trinucleotide SSR motifs in wheat. Theor Appl Genet. 2002;104:286–93.

    Article  CAS  PubMed  Google Scholar 

  18. Guyomarc'h H, Sourdille P, Edwards K, Bernard M. Studies of the transferability of microsatellites derived from Triticum tauschii to hexaploid wheat and to diploid related species using amplification, hybridization and sequence comparisons. Theor Appl Genet. 2002;105:736–44.

    Article  CAS  PubMed  Google Scholar 

  19. Sourdille P, Tavaud M, Charmet G, Bernard M. Transferability of wheat microsatellites to diploid Triticeae species carrying the a, B and D genomes. Theor Appl Genet. 2001;103:346–52.

    Article  CAS  Google Scholar 

  20. Guyomarc'h H, Sourdille P, Charmet G, Edwards K, Bernard M. Characterisation of polymorphic microsatellite markers from Aegilops tauschii and transferability to the D-genome of bread wheat. Theor Appl Genet. 2002;104:1164–72.

    Article  CAS  PubMed  Google Scholar 

  21. Pestsova E, Ganal EW, Röder MS. Isolation and mapping of microsatellite markers specific for the D genome of bread wheat. Genome. 2000;43:689–97.

    Article  CAS  PubMed  Google Scholar 

  22. Sourdille P, Singh S, Cadalen T, Brown-Guedira GL, Gay G, Qi LL, et al. Microsatellite-based deletion bin system for the establishment of genetic-physical map relationships in wheat (Triticum aestivum L.). Funct Integr Genomics. 2004;4:12–25.

    Article  CAS  PubMed  Google Scholar 

  23. Somers DJ, Isaac P, Edwards K. A high-density microsatellite consensus map for bread wheat (Triticum aestivum L.). Theor Appl Genet. 2004;109:1105–14.

    Article  CAS  PubMed  Google Scholar 

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

    Article  Google Scholar 

  25. IWGSC. A chromosome-based draft sequence of the hexaploid bread wheat (Triticum aestivum) genome. Science. 2014;345:1251788.

    Article  CAS  Google Scholar 

  26. Avni R, Nave M, Barad O, Baruch K, Twardziok SO, Gundlach H, et al. Wild emmer genome architecture and diversity elucidate wheat evolution and domestication. Science. 2017;357:93–7.

    Article  CAS  PubMed  Google Scholar 

  27. Luo MC, Gu YQ, Puiu D, Wang H, Twardziok SO, Deal KR, et al. Genome sequence of the progenitor of the wheat D genome Aegilops tauschii. Nature. 2017;551:498–502.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Wang SC, Wong D, Forrest K, Allen A, Chao S, Huang BE. Characterization of polyploid wheat genomic diversity using a high-density 90,000 single nucleotide polymorphism array. Plant Biotechnol. 2014;12:787–96.

    Article  CAS  Google Scholar 

  29. Winfield MO, Allen AM, Burridge AJ, Barker GLA, Benbow HR, Wilkinson PA, et al. High-density SNP genotyping array for hexaploid wheat and its secondary and tertiary gene pool. Plant Biotechnol. 2016;14:1195–206.

    Article  CAS  Google Scholar 

  30. Allen AM, Winfield MO, Burridge AJ, Downie RC, Benbow HR, Barker GL, et al. Characterization of a wheat breeders’ Array suitable for high-throughput SNP genotyping of global accessions of hexaploid bread wheat (Triticum aestivum). Plant Biotechnol J. 2017;15:390–401.

    Article  CAS  PubMed  Google Scholar 

  31. Jin H, Wen WE, Liu JD, Zhai SN, Zhang Y, Yan J, et al. Genome-wide QTL mapping for wheat processing quality parameters in a Gaocheng 8901/Zhoumai 16 recombinant inbred line population. Front Plant Sci. 2016;7:1032S.

    Google Scholar 

  32. Cui F, Zhang N, Fan XL, Zhang W, Zhao CH, Yang LJ, Pan R, Mei C, Jie H, Zhao X. Utilization of a Wheat660K SNP array-derived high-density genetic map for high-resolution mapping of a major QTL for kernel number. Sci Rep. 2017;7:3788.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  33. Tian XL, Wen WE, Xie L, Fu LP, Xu DA, Fu C, et al. Molecular mapping of reduced plant height gene Rht24 in bread wheat. Front Plant Sci. 2017;8:1379.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Liu JD, He ZH, Wu L, Bai B, Wen WE, Xie CJ, et al. Genome-wide linkage mapping of QTL for black point reaction in bread wheat (Triticum aestivum L.). Theor Appl Genet. 2016;129:1–12.

    Article  CAS  Google Scholar 

  35. Wang LX, Qiu J, Chang LF, Liu LH, Li HB, Pang BS, et al. Assessment of wheat variety distinctness using SSR markers. J Integr Agric. 2015;14:1923–35.

    Article  Google Scholar 

  36. Ramirez-Gonzalez RH, Segovia V, Bird N, Fenwick P, Holdgate S, Berry S, et al. RNA-Seq bulked segregant analysis enables the identification of high-resolution genetic markers for breeding in hexaploid wheat. Plant Biotechnol J. 2015;13:613–24.

    Article  CAS  PubMed  Google Scholar 

  37. Wang RX, Zhang XY, Wu L, Wang R, Hai L, Yan CS, et al. QTL mapping for grain filling rate and thousand-grain weight in different ecological environments in wheat. Acta Agron Sin. 2008;34:1750–6.

    Article  CAS  Google Scholar 

  38. Groos C, Robert N, Bervas E, Charmet G. Genetic analysis of grain protein-content, grain yield and thousand-kernel weight in bread wheat. Theor Appl Genet. 2003;106:1032–40.

    Article  CAS  PubMed  Google Scholar 

  39. Hu YJ, Li MF, Yang DL, Liu Y, Cheng HB, Chang L, et al. Integration and meta-analysis of quantitative trait loci for water-soluble carbohydrate content in wheat grain. J Triticeae Crops. 2016;36:989–95.

    CAS  Google Scholar 

  40. Araki E, Miura H, Sawada S. Identification of genetic loci affecting amylose content and agronomic traits on chromosome 4A of wheat. Theor Appl Genet. 1999;98:977–84.

    Article  CAS  Google Scholar 

  41. Mccartney CA, Somers DJ, Lukow O, Ames N, Noll J, Cloutier S, et al. QTL analysis of quality traits in the spring wheat cross RL4452 × 'AC Domain'. Plant Breed. 2006;125:565–75.

    Article  CAS  Google Scholar 

  42. Tian B, Liu B, Zhu ZL, Xie QG, Tian JC. Conditional and unconditional QTL mapping of grain starch accumulation in wheat. Sci Agric Sin. 2011;44:4551–9.

    Google Scholar 

  43. Wang JS, Liu WH, Wang H, Li LH, Wu J, Yang XM, et al. QTL mapping of yield-related traits in the wheat germplasm 3228. Euphytica. 2011;177:277–92.

    Article  Google Scholar 

  44. Cui F, Zhao CH, Ding AM, Li J, Wang L, Li XF, et al. Construction of an integrative linkage map and QTL mapping of grain yield-related traits using three related wheat RIL populations. Theor Appl Genet. 2014;127:659–75.

    Article  PubMed  Google Scholar 

  45. Tian B, Deng ZY, Xie QG, JC JT. Genetic dissection of the developmental behaviour of total starch content and its components in wheat grain. Crop Pasture Sci. 2015;66:445–55.

    Article  CAS  Google Scholar 

  46. Deng ZY, Hu SN, Chen F, Li WJ, Chen JS, Sun CL, et al. Genetic dissection of interaction between wheat protein and starch using three mapping populations. Mol Breed. 2015;35:12.

    Article  CAS  Google Scholar 

  47. Deng ZY, Fang WQ, Guo X, Zhao XX, Guo H, Hu SN, et al. Genetic dissection of interactions between wheat flour starch and its components in two populations using two QTL mapping methods. Mol Breed. 2018;38:41.

    Article  CAS  Google Scholar 

  48. Nakamura T, Vrinten P, Saito M, Konda M. Rapid classification of partial waxy wheats using PCR-based markers. Genome. 2002;45:1150–6.

    Article  CAS  PubMed  Google Scholar 

  49. Uematsu K, Suzuki N, Iwamae T, Inui M, Yukawa H. Increased fructose 1, 6-bisphosphate aldolase in plastids enhances growth and photosynthesis of tobacco plants. J Exp Bot. 2012;68:3001–9.

    Article  CAS  Google Scholar 

  50. Iwaki T, Wadano A, Yokota A, Himeno M. Aldolase-an important enzyme in controlling the ribulose 1, 5-bisphosphate regeneration rate in photosynthesis. Plant Cell Physiol. 1991;32:1083–91.

    Article  CAS  Google Scholar 

  51. Wang LJ. Research progress of sucrose invertase in higher plants. J Anhui Agric Sci. 2014;42:8108–11.

    CAS  Google Scholar 

  52. Liao JH, Li JM, Wu Y. Research progress in starch synthase of wheat. J Guizhou Agric Sci. 2014;42:18–22.

    CAS  Google Scholar 

  53. North HM, De Almeida A, Boutin JP, Frey A, To A, Botran L, et al. The Arabidopsis ABA-deficient mutant aba4 demonstrates that the major route for stress-induced ABA accumulation is via neoxanthin isomers. Plant J. 2010;50:810–24.

    Article  CAS  Google Scholar 

  54. Wang Q, Wen XG, Lu CM, Zhang QD. Photosynthetic functions of different senescing leaves in the canopy of super high-yield rice ‘Hua-An3’. Acta Phytoecol Sin. 2004;28:39–46.

    CAS  Google Scholar 

  55. Paulsen H. Chlorophyll a/b-binding proteins. Photochem Photobiol. 1995;62:367–82.

    Article  CAS  Google Scholar 

  56. Chourey PS, Nelson OE. The enzymatic deficiency conditioned by the shrunken-1 mutations in maize. Biochem Genet. 1976;14:1041–55.

    Article  CAS  PubMed  Google Scholar 

  57. Hou J, Jiang QY, Hao CY, Wang YQ, Zhang HN, Zhang XY. Global selection on sucrose synthase haplotypes during a century of wheat breeding. Plant Physiol. 2014;164:1918–29.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Jiang QY, Hou J, Hao CY, Wang LF, Ge HF, Dong YS, Zhang XY. The wheat (T. aestivum) sucrose synthase 2 gene (TaSus2) active in endosperm development is associated with yield traits. Funct Integr Genomics. 2011;11:49–61.

    Article  CAS  PubMed  Google Scholar 

  59. Guillemaut P, Laurence MD. Isolation of plant DNA: a fast, inexpensive, and reliable method. Plant Mol Biol Report. 1992;10:60–5.

    Article  CAS  Google Scholar 

  60. Shahinnia F, Sayed-Tabatabaei BE. Conversion of barley SNPs into PCR-based markers using dCAPS method. Genet Mol Biol. 2009;32:564–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Thiel T, Kota R, Grosse I, Stein N, Graner A. SNP2CAPS: a SNP and INDEL analysis tool for CAPS marker development. Nucleic Acids Res. 2004;32:e5.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Nyquist WE, Baker RJ. Estimation of heritability and prediction of selection response in plant populations. Crit Rev Plant Sci. 1991;10:235–322.

    Article  Google Scholar 

  63. Stam P. Construction of integrated genetic linkage maps by means of a new computer package: join map. Plant J. 1993;3:739–44.

    Article  CAS  Google Scholar 

  64. Kosambi DD. The estimation of map distances from recombination values. Ann Eugenics. 1943;12:172–5.

    Article  Google Scholar 

  65. Li HH, Ye GY, Wang JK. A modified algorithm for the improvement of composite interval mapping. Genetics. 2007;175:361–74.

    Article  PubMed  PubMed Central  Google Scholar 

  66. Wang J. Inclusive composite interval mapping of quantitative trait genes. Acta Agron Sin. 2009;35:239–45.

    Article  CAS  Google Scholar 

  67. Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative CT method. Nat Protoc. 2008;3:1101–8.

    Article  CAS  PubMed  Google Scholar 

Download references


The authors are grateful to Prof. R. A. McIntosh, at Plant Breeding Institute, University of Sydney, for critical review of this manuscript.


This work was funded by the National Key Research and Development Programs of China (2016YFD0100502), National Natural Science Foundation of China (31571663) and CAAS Science and Technology Innovation Program. The planting and phenotyping of materials were funded by National Natural Science Foundation of China (31571663) and CAAS Science and Technology Innovation Program, whereas the genotyping, data analysis and manuscript writing were funded by National Key Research and Development Programs of China (2016YFD0100502).

Author information

Authors and Affiliations



FC, TXL, FLP, XDA, WY and XXT performed the experiments and data analysis. DJY created and provided DEK lines for the research. CSH and ZY designed the experiments. FC and CSH wrote the draft. XXC and HZH revised the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Yan Zhang or Shuanghe Cao.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary information

Additional file 1: Figure S1.

The phenotypes of BL31 (a) and BL33 (b) grains.

Additional file 2: Figure S2.

The frequency distributions of polymorphic SNPs on 21 wheat chromosomes. The chromosomes with potential Dek QTL are highlighted in orange. The frequency threshold of 0.08 is indicated by a horizontal line. Distribution frequency is the percentage of polymorphic SNPs in the total ones on each chromosome.

Additional file 3: Table S1.

SSR markers used in this study.

Additional file 4: Table S2.

SNPs used to develop CAPS markers.

Additional file 5: Figure S3.

Design of CAPS marker for AX-109027972 (a), AX-108996126 (b), AX-109301653 (c), AX-109516011 (d) and AX-110042240 (e). The arrows and boxes show the target sites and the positions of restriction sites for the CAPS markers, respectively. The underlined forward and reverse primers for CASP markers are labeled at specific target sites.

Additional file 6: Figure S4.

PCR profile of molecular markers to enrich the target regions of QTL. (a) AX-109027972. (b) AX-108996126. (c) AX-109301653. (d) AX-109516011. (e) AX-110042240. M, Marker (DL2000DNA ladder in a, b, c and e; 20 bp DNA ladder in d; Takara Bio Company). Target fragments are shown with red arrows. Lanes 1 and 2 indicate BL33 and BL31, respectively; lane 3 shows F1 plants; lanes 4–10 are F2 individual plants.

Additional file 7: Table S3.

Comparison of TGW between BL31 and BL33.

Additional file 8: Figure S5.

LOD contours obtained by inclusive composite interval mapping (ICIM) analysis for TGW QTL on chromosomes 3B in F2 and F2:3 populations derived from BL33/BL31. The genetic region of QTGW.caas-3BS is highlighted in red. LOD contour in Gaoyi 2016–2017 cropping season and Gaoyi 2017–2018 cropping season indicated with red and green lines, respectively. LOD threshold of 2.5 is indicated by the dotted vertical line.

Additional file 9: Figure S6.

The comparison of grain starch contents in mature grains between BL31 and BL33.

Additional file 10: Figure S7.

The change of grain sucrose contents after flowering.

Additional file 11: Table S4.

Annotated genes in the physical intervals of QTL.

Additional file 12: Table S5.

Transcriptional patterns of candidate genes in the developing grains.

Additional file 13: Figure S8.

Transcriptional patterns of four candidate genes in the developing grains.

Additional file 14: Figure S9.

PCR profile of the Wx-B1 specific marker. M, Marker (DL2000DNA ladder; Takara Bio Company). Target fragments are shown with red arrows. Lane 1 indicates Chinese Spring; lane 2 indicates Guandong107; lane 3 and 4 show BL33 and BL31, respectively.

Additional file 15: Table S6.

Primers used in quantitative real-time PCR.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Fu, C., Du, J., Tian, X. et al. Rapid identification and characterization of genetic loci for defective kernel in bread wheat. BMC Plant Biol 19, 483 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: