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

doi:10.21203/rs.2.11611/v1

Download PDF

Research article

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

https://doi.org/10.21203/rs.2.11611/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 08 Nov, 2019

Read the published version in BMC Plant Biology →

You are reading this older preprint version

Read the latest preprint version →

Background

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.

Results

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 determined according to International Wheat Genome Sequence Consortium (IWGSC) RefSeq v1.0.

Conclusions

Overall, we achieved a rapid QTL mapping and reasonable predication of causal genes for Dek, laying foundation to conduct fine mapping and dissect the regulatory mechanism underlying Dek trait in wheat.

Plant Molecular Biology and Genetics

Defective kernel

Grain filling

Rapid QTL mapping

SNP array

Triticum aestivem L.

Wheat (Triticum aestivum L.) supplies staple food for more than 35% of population worldwide (http://www.wheatinitiative.org/sites/default/files/attached_file/wheat_facts.pdf). 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-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-10]. With rapid advancements in sequencing technology, high-throughput sequencing or chip-based SNP genotyping platforms have been used for QTL mapping [11-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.

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 F₂ plants and their corresponding F_2:3 lines from the cross between BL31 and BL33 were created for QTL mapping.

Field trials and phenotype evaluation

The F₂ population and parents were planted at Gaoyi (37°33′N, 114°26′E) in Hebei province during the 2016 - 2017 cropping season. The F_2: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 F₂ plants was scored using the percentages of defective kernels (defined as Dek rate) in single plants. For F_2: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 [60]. The Wheat660K SNP array containing 630,517 SNPs, a high-throughput genotyping tool developed by Institute of Crop Sciences, Chinese Academy of Agricultural Sciences (https://wheat.pw.usda.gov/ggpages/topics/Wheat660_SNP_array_developed_by_CAAS.pdf), was used to test two isogenic lines BL31 and BL33 at CapitalBio Corporation (Beijing, China; http://www.capitalbio.com/). The frequency distribution of polymorphic SNPs in each chromosome was analyzed using Microsoft Excel 2016 (https://www.microsoft.com/zh-cn/). 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; http://www.ht-biotech.net/), 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; http://www.genecompany.com/). 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 250V using the JY 600HE Electrophoretic Apparatus (Beijing Junyi Electrophoresis Co., Ltd., Beijing; http://www.bjjunyi.com/) or 2% agarose gels at 220V using the DYY-7C Electrophoresis Apparatus (Beijing Liuyi Biotechnology Co., Ltd.; http://www.ly.com.cn/).

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 [61], 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 (https://www.lynnon.com/index.html) for designing chromosome-specific CAPS markers [62]. 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 (http://www.sas.com) with PROC MIXED procedure. Broad-sense heritability (h_B²) was calculated using the following formula [63]:

h_B² =σ²_g / (σ²_g +σ²_ge/ r +σ²_ε/ re)

in which, σ²_g, σ²_ge and σ²_ε 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 (http://www.kyazma.nl) was used for linkage map construction [64] and genetic distances between markers were estimated based on the Kosambi mapping function [65]. QTL analysis was conducted by inclusive composite interval mapping (ICIM) based on the Dek rate using IciMapping 4.0 software (http://www.isbreeding.net) [66, 67]. 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

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 using Plant Sucrose Assay kit (Beijing Solarbio Science & Technology Co., Ltd., Beijing, http://www.solarbio.net/). Mature grains were used to measure starch contents with Starch Assay Kit (Beijing Solarbio Science & Technology Co., Ltd., Beijing, http://www.solarbio.net/).

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., http://www.tiangen.com/). Reverse transcription was performed using the reverse transcriptase M-MLV (RNase H^-) (Takara Biomedical Technology Co., Ltd., http://www.takara.com). The iTaq^TM 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 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 [68]. 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).

Phenotypic analysis

The kernels of BL31 were normal and plump, whereas the Dek rate of BL33 is 100% in all environments (Additional file 1: Figure S1). The Dek rates in F₂ and F_2:3 populations were highly correlated among three environments with correlation coefﬁcients 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) and the broad-sense heritability of Dek rate in F_2:3 across two environments was 0.95 (Table 2), 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].

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: Figure S2). 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: Table S1) [17-24]. 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 F₂ population, and two linkage groups were constructed in chromosomes 3B and 4A, including 6 and 9 SSR markers, respectively (Fig. 2).

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 (http://plants.ensembl.org/Triticum_aestivum/Info/Index) [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.

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: Table S2; Figure 2; Additional file 5: Figure S3; Additional file 6: Figure S4 a-e). 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: Table S3). 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: Figure S5). 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-7]. Our biochemical analysis also showed that the mature grain starch contents of BL33 were significantly lower than those of BL31 (Additional file 9: Figure S6). In contrast, BL33 had higher sucrose contents and slower sucrose reduction across the filling stages than BL31 (Additional file 10: Figure S7). 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 (Table S4). 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: Table S4).

To further identify candidate genes for the Dek QTL, the spatio-temporal expression patterns of the nine genes were analyzed using the WheatExp (https://wheat.pw.usda.gov/WheatExp/). Of these, four genes, TraesCS3B01G025200, TraesCS3B01G028700, TraesCS3B01G039800, and TraesCS4A01G446700, had an observable expression during grain development (Additional file 12: Table S5) [25]. 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: Figure S8). 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 [26-28]. Consequently, a few sequencing-based and chip-based SNP genotyping platforms have been developed and used for genetic analyses in wheat [29-31]. Some studies have achieved fine mapping of QTL using SNP chips [32-35]. 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 [36, 37]. 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).

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 [38]. QTgw.nfcri-3B for TGW was also linked to Xgwm533 [38]. 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. [39]. Additionally, QMswscf.acs-3B.1 for water soluble carbohydrates (WSC) of main stem at the early anthesis stage, were linked to Xcfd79 [40].

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 [41-48]. QGsc4A, a QTL for grain starch contents, is closely linked to Wx-B1 [43]. Qftsc4A for total grain starch content and Qfams4A for grain amylose content were also linked to Wx-B1 [47]. To identify the genetic relationship between QDek.caas-4AL and Wx-B1, the gene-specific functional marker developed by Nakamura et al. [49] was used to genotype the parents, BL31 and BL33. The result showed that no polymorphism was detected between them (Additional file 14: Figure 9), 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 [50, 51]. β-fructofuranosidase is a key enzyme for synthesis, metabolism and transportation of sucrose and starch, and can irreversibly catalyze sucrose hydrolysis to glucose and fructose [52, 53]. 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 [54]. Neoxanthin is necessary for maintaining the integrity of the light harvesting chlorophyll and photosynthetic function of leaves, involving the metabolism of photosynthetic product [55, 56]. TaSus-4A may be the causal gene of QDek.caas-4AL. Sus deficiency led to the shrunken-1 mutations in maize [57]. 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 [58, 59].

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.

In the study, we detected three main QTLs controlled the DEK, and single QTL could explain more than 10% phenotypic variation. Four genes, TaFba-3B, TaBff-3B, TaABA4-3B, and TaSus-4A, were believed as the candidate genes and had the feature of time-specific expression.

AGPP: adenosine diphosphate glucose pyrophosphorylase; ANOVA: analysis of variance; CAPS: cleaved amplified polymorphic sequence; CWI: cell wall invertase; DAF: days after flowering; Dek: defective kernel; Fba: fructose-bisphosphate aldolase; GFD: grain filling duration; GFR: grain filling rate; GWAS: genome-wide association study; ICIM: inclusive composite interval mapping; IWGSC: International Wheat Genome Sequence Consortium; KASP: kompetitive allele specific polymerase chain reaction; LOD: logarithm of odds; PVE: phenotype variation explained; qPCR: quantitative real-time PCR; SNP: single nucleotide polymorphism; SSR: simple sequence repeat; Sus: sucrose synthase; TGW: thousand grain weight; WSC: water soluble carbohydrates

Authors' contributions

HZH and XXC conceived the study. DJY provided test material for the research. ZY and CSH designed the experiments. FC performed the experiments and wrote the paper. TXL, FLP, XDA, WY and XXT offered help for field trials and data analysis. All authors discussed the results and commented on the manuscript. All authors read and approved the final manuscript.

Funding

This work was funded by the National Key Research and Development Programs of China (2016YFD0100502, 2016YFD0101802, 2016YFE0108600), and CAAS Science and Technology Innovation Program. There is on role of the funding bodies in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Availability of data and materials

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

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Author details

1 Institute of Crop Sciences, National Wheat Improvement Center, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100081, China

2 Wheat Research Institute, Gansu Academy of Agricultural Sciences, Lanzhou 730070, China

3 International Maize and Wheat Improvement Center (CIMMYT), c/o CAAS, 12 Zhongguancun South Street, Beijing 100081, China

*Authors for correspondence: Shuanghe Cao: [email protected]; Yan Zhang: [email protected]

Cheng SH, Yang SM, Zhang BQ, Ji KZ, Zhao BH, Gao DR. Kernel plumpness and filling index in wheat (Triticum aestivum). Jiangsu J Agric Sci. 1993;9:7-10.
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.
Neuffer MG, Sheridan WF. Defective kernel mutants of maize. I. genetic and lethality studies. Genetics. 1980;95:929-44.
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.
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.
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.
Li XJ, Wei G, Sun SL, Chen ZL, C. Jing, W. Song, et al. Defective Kernel 39 encodes a PPR protein required for seed development in maize. J Integr Plant Biol. 2018;60:45-64.
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.
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.
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.
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.
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.
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.
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.
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.
Cao SH, Guo XL, Liu DC, Zhang XQ, Zhang AM. Preliminary gene-mapping of photoperiod-temperature sensitive genic male sterility in wheat (Triticum aestivum ). Acta Genet Sin. 2004;31:293-8.
Song QJ, Fickus EW, Cregan PB. Characterization of trinucleotide SSR motifs in wheat. Theor Appl Genet. 2002;104:286–93.
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. 2002a;105:736-44.
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.
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. 2002b;104:1164-72.
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.
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 ). Funct Integr Genomics. 2004;4:12-25.
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.
Somers DJ, Isaac P, Edwards K. A high-density microsatellite consensus map for bread wheat (Triticum aestivum ). Theor Appl Genet. 2004;109:1105-14.
Zadoks JC, Chang TT, Konzak A decimal code for the growth stages of cereals. Weed Res. 1974;14:415-21.
IWGSC, . A chromosome-based draft sequence of the hexaploid bread wheat (Triticum aestivum) genome. Science. 2014;345:1251788.
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-97.
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.
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.
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.
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-
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.
Cui F, Zhang N, Fan XL, Zhang W, Zhao CH, Yang LJ, R. Pan, C. Mei, H. Jie, X. Zhao, Utilization of a Wheat660K SNP array-derived high-density genetic map for high-resolution mapping of a major QTL for kernel number, Sci. Rep. 7 (2017) 3788.
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.
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). Theor Appl Genet. 2016;129:1-12.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Tian B, Deng ZY, Xie QG, J. Tian JC. Genetic dissection of the developmental behaviour of total starch content and its components in wheat grain. Crop Pasture Sci. 2015;66:445-55.
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.
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.
Nakamura T, Vrinten P, Saito M, Konda M. Rapid classification of partial waxy wheats using PCR-based markers. Genome. 2002;45:1150-6.
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.
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.
Wang LJ, Research progress of sucrose invertase in higher plants. J Anhui Agric Sci. 2014;42:8108-11.
Liao JH, Li JM, Wu Y. Research progress in starch synthase of wheat. J Guizhou Agric Sci. 2014;42:18-22.
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.
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.
Paulsen H. Chlorophyll a/b‐binding proteins. Photochem Photobiol. 1995;62:367-82.
Chourey PS, Nelson OE. The enzymatic deficiency conditioned by the shrunken-1 mutations in maize. Biochem Genet. 1976;14:1041-55.
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.
Jiang QY, Hou J, Hao CY, Wang LF, Ge HF, Dong YS, Zhang XY. The wheat ( aestivum) sucrose synthase 2 gene (TaSus2) active in endosperm development is associated with yield traits. Funct Integr Genomics. 2011;11:49-61.
Guillemaut P, Laurence MD. Isolation of plant DNA: A fast, inexpensive, and reliable method. Plant Mol Biol Report. 1992;10:60-5.
Shahinnia F, Sayed-Tabatabaei BE. Conversion of barley SNPs into PCR-based markers using dCAPS method. Genet Mol Biol. 2009;32:564-7.
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.
Nyquist WE, Baker RJ. Estimation of heritability and prediction of selection response in plant populations. Crit Rev Plant Sci. 1991;10:235-322.
Stam P. Construction of integrated genetic linkage maps by means of a new computer package: Join Map. Plant J. 1993;3:739-44.
Kosambi DD. The estimation of map distances from recombination values. Ann Eugen. 1943;12:172-5.
Li HH, Ye GY, Wang JK. A modified algorithm for the improvement of composite interval mapping. Genetics. 2007;175:361-74.
Wang J. Inclusive composite interval mapping of quantitative trait genes. Acta Agron Sin. 2009;35:239-45.
Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative C_T Nat Protoc. 2008;3:1101-8.

Table 1 Statistical analysis of Dek rates for F₂ and F_2:3 populations across three environments

Environment	Range	Mean ± SD^a	Correlation coefficient
Environment	Range	Mean ± SD^a	Gaoyi 2018	Xinxiang 2018
Gaoyi 2017^b	0~1.00	0.33 ± 0.32	0.70^**	0.73^**
Gaoyi 2018^c	0~1.00	0.24 ± 0.32		0.92^**
Xinxiang 2018^d	0~1.00	0.22 ± 0.29

^a SD, standard deviation.

^bGaoyi 2017, Gaoyi, 2016-2017 cropping season.

^cGaoyi 2018, Gaoyi, 2017-2018 cropping season.

^dXinxiang 2018, Xinxiang, 2017-2018 cropping season.

^** indicates significance levels at P < 0.01.

Table 2 Analysis of variance of Dek rates for F_2:3 lines derived from BL33/BL31

Source of variation	Degrees of freedom	Sum of square	Mean square	F value
Environment	1	0.097	0.097	5.64^*
Genotype	145	51.76	0.36	20.68^***
Genotype × Environment	145	2.21	0.015	0.88
Error	290	5.01	0.017

^* indicates significance levels at P < 0.05.

^*** indicates significance levels at P < 0.001.

Table 3 QTL for Dek trait detected in the F₂and F_2:3 populations across different environments

Environment	QTL	Flanking marker	Physical interval (Mb)^a	LOD	Add^b	PVE (%)^c
Gaoyi 2017^d	QDek.caas-3BS.2	Xcfd79-AX-109516011	18.84-20.00	2.84	0.13	16.61
	QDek.caas-4AL	Xwmc500-AX-109301653	704.34-705.47	3.98	0.06	28.19
Xinxiang 2018^e	QDek.caas-3BS.1	Xbarc133-Xcfd79	7.60-18.84	4.00	0.14	18.17
	QDek.caas-4AL	Xwmc500-AX-109301653	704.34-705.47	4.01	0.12	22.60
Gaoyi 2018^f	QDek.caas-3BS.1	Xbarc133-Xcfd79	7.60-18.84	3.26	0.14	16.39
	QDek.caas-3BS.2	Xcfd79-AX-109516011	18.84-20.00	3.25	0.17	21.83
	QDek.caas-4AL	Xwmc500-AX-109301653	704.34-705.47	3.03	0.09	19.08
Average	QDek.caas-3BS.1	Xbarc133-Xcfd79	7.60-18.84	2.98	0.11	14.78
	QDek.caas-3BS.2	Xcfd79-AX-109516011	18.84-20.00	3.49	0.13	21.80
	QDek.caas-4AL	Xwmc500-AX-109301653	704.34-705.47	3.95	0.09	23.38

^a Physical interval defined by physical positions of flanking markers according to IWGSC RefSeq v1.0.

^b Add, additive effect; positive values indicate an increasing effect from BL33.

^cPVE, phenotypic variance explained.

^dGaoyi 2017, Gaoyi, 2016-2017 cropping season.

^eGaoyi 2018, Gaoyi, 2017-2018 cropping season.

^fXinxiang 2018, Xinxiang, 2017-2018 cropping season.

Table 4 Candidate genes of Dek QTLs

QTL	Temporary name	Gene ID	Candidate gene
QDek.caas-3BS.1	TaFba-3B	TraesCS3B01G025200	Fructose-bisphosphate aldolase
	TaBff-3B	TraesCS3B01G028700	β-fructofuranosidase
QDek.caas-3BS.2	TaABA4-3B	TraesCS3B01G039800	Abscisic acid-deficient 4
QDek.caas-4AL	TaSus-4A	TraesCS4A01G446700	Sucrose synthase

Download PDF

Journal Publication

published 08 Nov, 2019

Read the published version in BMC Plant Biology →

Editorial decision: Minor revision
01 Oct, 2019
Review #2 received at journal
16 Sep, 2019
Review #1 received at journal
16 Sep, 2019
Reviewer #2 agreed at journal
20 Aug, 2019
Reviewers invited by journal
06 Aug, 2019
Reviewer #1 agreed at journal
06 Aug, 2019
Submission checks completed at journal
16 Jul, 2019
Editor invited by journal
11 Jul, 2019
Editor assigned by journal
11 Jul, 2019
First submitted to journal
08 Jul, 2019

You are reading this older preprint version

Read the latest preprint version →

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

Status:

Journal Publication

Version 1

Abstract

Figures

Background

Methods

Results

Discussion

Conclusions

Abbreviations

Declarations

References

Tables

Supplementary Files

Status:

Journal Publication

Version 1