ORIGINAL RESEARCH article

Front. Plant Sci., 06 October 2022

Sec. Plant Breeding

Volume 13 - 2022 | https://doi.org/10.3389/fpls.2022.950664

Genetic dissection of QTLs for starch content in four maize DH populations

  • 1. Quality and Safety Institute of Agricultural Products, Heilongjiang Academy of Agricultural Sciences, Harbin, China

  • 2. Institute of Advanced Agricultural Technology, Qilu Normal University, Jinan, China

  • 3. Maize Yufeng Biotechnology LLC, Beijing, China

Article metrics

View details

7

Citations

2k

Views

868

Downloads

Abstract

Starch is the principal carbohydrate source in maize kernels. Understanding the genetic basis of starch content (SC) benefits greatly in improving maize yield and optimizing end-use quality. Here, four double haploid (DH) populations were generated and were used to identify quantitative trait loci (QTLs) associated with SC. The phenotype of SC exhibited continuous and approximate normal distribution in each population. A total of 13 QTLs for SC in maize kernels was detected in a range of 3.65–16.18% of phenotypic variation explained (PVE). Among those, only some partly overlapped with QTLs previously known to be related to SC. Meanwhile, 12 genes involved in starch synthesis and metabolism located within QTLs were identified in this study. These QTLs will lay the foundation to explore candidate genes regulating SC in maize kernel and facilitate the application of molecular marker-assisted selection for a breeding program to cultivate maize varieties with a deal of grain quality.

Introduction

Maize (Zea mays L.) is one of the important crops in the world. Starch is the main component of maize kernel, which accounts for about 70% of kernel weight and is the main energy source to supply adequate food to humans and feed animals and is a source of demand in bio-ethanol production and other industrial applications (Nelson and Pan, 1995; Balter, 2007). Therefore, understanding processes related to starch biosynthesis in maize kernel is of fundamental biological interest, and such knowledge also benefits agricultural applications by providing means to improve grain yield and quality in seed production (Comparot-Moss and Denyer, 2009; Jeon et al., 2010).

To elucidate genetic variation in starch biosynthesis and regulation, over the past few years, numerous quantitative trait locus (QTL) studies have been carried out using different mapping methods and populations, and many QTLs associated with starch content (SC) in maize kernel have been revealed (Liu et al., 2008; Wassom et al., 2008; Zhang et al., 2008; Li et al., 2009; Wang Y. Z. et al., 2010; Cook et al., 2012; Guo et al., 2013; Nancy et al., 2014; Yang et al., 2014; Dong et al., 2015; Alves et al., 2019). For instance, six QTLs were identified in a recombinant inbred line (RIL) population, and four genes were considered candidate genes likely acting as direct regulators of starch biosynthesis (Wang et al., 2015). A genome-wide association study (GWAS) revealed 27 QTLs involving 39 candidate genes that are associated with amylose content, including transcription factors, glycosyltransferases, glycosidases, and hydrolases (Li et al., 2018). Eight QTLs were revealed in an intermated B73×Mo17 population, and the gene encoding GLABRA2 expression modulator was nominated as the candidate gene of a major QTL by GWAS and the gene co-expression network analysis (Lin et al., 2019). Recently, in a total of 50, 18 novel QTLs were identified by integrating single linkage mapping, joint linkage mapping, and GWAS in a multi-parent population containing six RIL populations (Hu et al., 2021). Thus, although common QTLs associated with SC has been reported in different studies, unique QTLs were always specialized in the context of distinct populations and parents, which prompt us to conduct further investigations using relevant germplasm to extend our knowledge on the genetic basis of SC.

Distinct strategies of mapping populations were featured with strengths and limitations, which caused great impacts on the outputs of QTL (Odell et al., 2022). In particular, different types of populations tend to vary with two main characteristics: (1) their ability to capture genetic diversity and (2) their power to detect QTL of small effect (Odell et al., 2022). Double haploid (DH) segregating populations have been commonly used in QTL analysis for several particular advantages (Chaikam et al., 2019). Complete homozygosity of DH lines allows more accurate phenotyping over multiple locations and years compared to families in early selfing generations (Foiada et al., 2015; Yan et al., 2017). In addition, DH populations enable the removal of any residual heterozygosity, ensuring replicates are genetically identical (Odell et al., 2022). Moreover, the relatively high genetic variance in DH lines increases selection response by stabilizing the heritability of various traits during perse and test cross evaluation (Bordes et al., 2006; Gallais and Bordes, 2007; Mayor and Bernardo, 2009). In this study, we utilized four DH populations derived from the practical breeding program to further dissect the genetic basis and QTLs controlling the phenotypic variation of SC in maize kernels. We intended to mine novel alleles and genes from DH populations to improve the starch content of advanced maize germplasms in-house.

Materials and methods

Plant materials

Four DH populations (SC1, SC2, SC3, and SC4) including 345, 275, 134, and 258 lines, respectively, were developed from biparental crosses of eight inbred lines exhibiting the variation in SC (Table 1). These parental lines belonged to elite inbred lines that were widely used by the breeding program at Maize Yufeng Biotechnology LLC to optimize grain nutritional quality. Parents of SC1 and SC2 belonged to Reid Yellow Dent germplasm, and parents of SC3 and SC4 belonged to maize Lancaste germplasm (Table 1). Plants were grown in Liaoning province, China (LN, 40°‘82′N, 123°56′E) with three replication blocks in 2021. Each line was grown in a single-row plot with a row length of 150 and 60 cm between rows under natural field conditions. All plants in each row were self-pollinated and harvested after maturity. Kernels from the middle part of three well-grown ears were bulked for the measurement of starch. We declared that all the collections of plant and seed specimens related to this study were performed in accordance with the relevant guidelines and regulations by the Ministry of Agriculture (MOA) of the People's Republic of China.

Table 1

Trait aPopulations
SC1SC2SC3SC4
Parents
means ± SD (%)AJ51700172.10 ± 0.47AJ51700375.10 ± 0.41KB51900173.91 ± 0.99KB51900372.99 ± 0.46
AJ51700273.67 ± 0.35AJ51700470.45 ± 0.51KB51900275.80 ± 0.34KB51900471.90 ± 0.89
p-valueb0.058ns0.009**0.0194*0.1312ns
DHs
Size345275134258
means ± SD (%)71.18 ± 1.3171.55 ± 2.3473.85 ± 1.6572.90 ± 1.98
Range (%)66.38–74.7661.22–75.4667.40–77.7164.87–76.97
1.454.412.113.25
1.091.351.291.31
1.533.681.343.26
h2 (%)f73.9078.2082.5274.94

Phenotypic performance, variance, and broad-sense heritability of starch content (SC) in the four double-haploid (DH) populations.

a

SC;

b

P value based on a t-test evaluating two parental lines;

c

genetic variance;

d

environmental variance;

e

residual variance;

f

broad-sense heritability (h2);

*

p ≤ 0.05,

**

p ≤ 0.01, ns no significant different.

Starch content measurement and phenotypic data analysis

The starch content in maize kernels was determined by using a near infrared reflectance (NIR) spectrometer (DA 7250, Perten Instruments Inc., Sweden). The reflectance spectra were collected in a range of 400 to 2500 nm with 10-nm intervals in the NIR region. Each sample was bulked with at least 50 kernels and scanned three times.

R Version 4.0.1 (www.R-project.org) was used to perform all statistical analyses as previously described (Zhang et al., 2021). The variances of SC were estimated by using the R ‘AOV' function. The model for the variance analysis was as follows: y = μ + αg + βe + ε, where αg is the effect of the gth line, βe is the effect of the eth environment, and ε is the error. The effects in the model were defined by random. The broad-sense heritability as h2 = /( + /e) (Knapp et al., 1985) was calculated by using these variance components, where is the genetic variance, is the residual error, and e is the number of environments. The best linear unbiased predictor (BLUP) value of each line was calculated to eliminate the influence of environmental effects by using a linear mixed model. The model was yij = μ + ei + fj + εij, where yij is the phenotypic value of individual j in environment i, μ is the grand mean, ei is the effect of different environments, fj is the genetic effect, and εij is the random error. Both genotype and environment were considered as random effects in the R function “LME4”.

Genotyping and constructing a genetic linkage map

All lines were genotyped with the GenoBaits Maize 1K marker panel containing 4,589 SNP markers that were developed by Mol Breeding Biotechnology Co., Ltd., Shijiazhuang, China (http://www.molbreeding.com/), based on genotyping by target sequencing platform in maize (Guo et al., 2019). The SNP positions in the B73 reference genome Version 3 were converted to those in Version 4 using CrossMap Version 2.0.5 (Zhao et al., 2014). SNPs with minor allele frequency (MAF) < 0.1 or missing rate > 0.6 were filtered out in each population. Finally, high-quality SNPs in each population were then used to construct the genetic linkage maps via the R/qtl package functions est.rf and est.map (Broman et al., 2003) with the kosambi mapping method.

QTL mapping

The quantitative trait loci were analyzed by the composite interval mapping (CIM) method implemented in Windows QTL Cartographer 2.5 (Wang S. et al., 2010). The genome was scanned at every 1.0 cM interval between markers using a 10 cM window size. The forward and backward stepwise regressions with five controlling markers were conducted to control the background from flanking markers. The empirical logarithm of the odds (LOD) threshold was determined by the 1,000 permutations at a significance (p < 0.05) and used to identify the significant QTLs (Churchill and Doerge, 1994). These threshold LOD values ranged from 2.86 to 3.09 in four DH populations, respectively. For clarity, we used 3.00 as the LOD threshold for the four DH populations. With the 1.5-LOD support interval method, the confidence interval for each QTL position was estimated (Lander and Botstein, 1989). The additive × additive epistatic interactions were performed by the “ICIM-EPI” method in IciMapping Version 4.2 (Li et al., 2008).

Results

Phenotypic variation and heritability in kernel SC

Eight inbred lines (SC with a range of 70.45–75.80%) were used to develop four DH populations, which included 134–345 lines, respectively (Table 1). The SC exhibited continuous and approximate normal distribution in each DH population with a range of 61.22–77.71% (Figure 1, Table 1). The analysis of variance (ANOVA) revealed that the genotype variance was greater than the environmental variance in all populations, indicating that phenotypic variations were mainly controlled by genetic factors. Broad-sense heritability estimates were calculated and showed moderate to high heritability for SC in the four DH populations, with a range of 73.90–82.52% (Table 1). Together, these results indicated that the majority of SC variations are controlled by genetic factors and are suitable for further QTL mapping.

Figure 1

Genotyping and genetic linkage map

All DH lines in four populations were genotyped using the GenoBaits Maize 1K marker panel containing 4,589 SNP markers and further refined by eliminating SNPs with MAF < 0.1 or missing rate > 0.6. This resulted in a total of 776, 1,164, 831, and 1,216 polymorphic SNPs, with their precise physical positions mapped on the B73 reference sequence Version 4. In each DH population, the missing rate of SNPs in most lines was <2% (Supplementary Figure S1). Four linkage maps spanned a total of 775.88, 858.86, 627.75, and 773.21 cM genetic distances, respectively (Supplementary Figure S2). The average genetic distance between every two adjacent markers was 1.01, 0.74, 0.76, and 0.64 cM in each DH population, respectively (Supplementary Table S1).

Identification of QTLs for SC in four DH populations

The quantitative trait locus mapping for SC in the four DH populations and the related genetic features are summarized in Table 2. In total, 13 QTLs were identified with a LOD threshold of above 3.00 in the four populations (Table 2, Figure 2) and distributed on chromosomes 1, 3, 4, 5, 7, 8, and 9. The average of QTL genetic intervals was 9.82 cM in a range of 4.95–17.20 cM. The average of QTL physical intervals was 46.51 Mb in a range of 5.60–149.78 Mb. The average of the total PVE explained by all identified QTLs in a population was 26.13 and ranged from 12.80 (SC1) to 41.57% (SC3). This was less than broad-sense heritability (Tables 1, 2), suggesting that only part of the QTLs have been detected in bi-parent populations.

Table 2

PopulationsQTLChr.aG-Peak (cM)bP-Peak (Mb)_V4cG-Range (cM)dP-Range (Mb)_V4eLODPVE%fAdd.gParent h+PVE(%) -ALLi
SC1qSC-1-1356.11199.1650.83–63.40192.97–200.454.575.060.27AJ51700112.80
qSC-1-2460.66229.3357.50–63.48208.9–240.73.023.650.23AJ517001
qSC-1-3814.457.179.87–22.355.86–16.435.176.38−0.30AJ517002
SC2qSC-2-1141.9561.6740.66–45.6160.19–83.653.343.830.41AJ51700331.40
qSC-2-2427.97147.2018.96–31.7318.47–156.858.8110.580.68AJ517003
qSC-2-3589.14197.2784.56–92.14183.23–197.2712.2216.120.84AJ517003
qSC-2-488.7626.155.36–10.6020.78–65.484.985.790.50AJ517003
SC3qSC-3-1325.61136.6022.24–29.2645.76–157.174.278.03−0.40KB51900241.57
qSC-3-2729.11143.8726.37–35.11138.27–143.877.5416.18−0.57KB519002
qSC-3-398.0935.224.85–17.0913.71–35.226.3313.32−0.52KB519002
SC4qSC-4-1132.3044.5229.38–35.8239.4–62.756.088.430.54KB51900318.73
qSC-4-2193.48282.6189.95–107.15279.59–302.145.338.36−0.54KB519004
qSC-4-3539.4766.0230.35–44.9118.85–168.633.815.320.43KB519003

Individual quantitative trait locus (QTL) for SC in the four DH populations.

a

Chromosome;

b

Genetic position in centimorgans (cM) of QTL with the highest LOD;

c

Physical position of QTL based on the B73 reference sequence (V4);

d

Genetic position range in centimorgans (cM) of QTL with the highest LOD;

e

Physical position range of QTL based on the B73 reference sequence (V4);

f

Percentage of the phenotypic variation explained by the additive effect of QTL;

g

Additive effect of QTL;

h

Parental source of QTL;

i

Percentage of the phenotypic variation explained by the additive effect of all QTL.

Figure 2

In SC1, three QTLs (qSC-1-1, qSC-1-2, and qSC-1-3) were detected on chromosomes 3, 4, and 8, respectively. The phenotypic variation explained by these QTLs was 3.65–6.38% and explained by the additive effect of all QTLs was 12.8%. The parental AJ517001 allele at qSC-1-1 and qSC-1-2 increased the SC, while the parental AJ517002 allele increased the SC at qSC-1-3.

In SC2, a total of four QTLs (qSC-2-1, qSC-2-2, qSC-2-3, and qSC-2-4) were identified on chromosomes 1, 4, 5, and 8, respectively. The phenotypic variation of qSC-2-1 and qSC-2-4 was <10%, whereas the phenotypic variation of qSC-2-2 and qSC-2-3 was 10.58 and 16.12%, respectively, suggesting that qSC-2-2 and qSC-2-3 are two major QTLs controlling SC in SC2. The phenotypic variation explained by the additive effect of all QTLs was 31.4% and the parental AJ517003 allele at all QTLs in this population increased the SC.

In SC3, a total of three QTLs (qSC-3-1, qSC-3-2, and qSC-3-3) were detected on chromosomes 3, 7, and 9, respectively. The phenotypic variation explained by each individual QTL was 8.03-16.18%. qSC-3-2 and qSC-3-3 were the major QTLs explaining the phenotypic variation of 16.18 and 13.32%, respectively. The phenotypic variation explained by the additive effect of all QTLs was 41.57%, and the allele from parental KB519002 at these QTLs increased the SC.

In SC4, a total of three QTLs (qSC-4-1, qSC-4-2, and qSC-4-3) were identified. Among those, qSC-4-1 and qSC-4-2 were located on chromosome 1, and qSC-4-3 was located on chromosome 5. The phenotypic variation explained by these QTLs was 5.32–8.43%, and the phenotypic variation explained by the additive effect of all QTLs was 18.73%. Parental KB519003 allele at qSC-4-1 and qSC-4-3 increased the SC, while parental KB519004 increased the SC at qSC-4-2.

Additionally, the additive × additive epistatic interaction analysis was also investigated. However, no epistatic interactions could be detected, indicating that the genetic component of SC in these DH populations is mainly characterized by additive gene actions.

Genetic overlap of QTLs in the four DH populations with other populations

To evaluate genetic overlaps among different mapping populations, the 1.5-LOD support interval of QTLs in the four DH populations and in other populations for the SC previously reported were compared. QTLs with overlapping support intervals were considered common QTLs. Among the four DH populations, a 2.56 Mb overlap was observed between qSC-2-1 and qSC-4-2.

Importantly, a total of 93 overlaps were detected after comparing the four DH populations with results from all types of other populations, including RILs (Cook et al., 2012; Guo et al., 2013; Yang et al., 2013; Nancy et al., 2014; Dong et al., 2015; Lin et al., 2019; Hu et al., 2021), F2:3 families (Liu et al., 2008; Zhang et al., 2008; Li et al., 2009; Wang Y. Z. et al., 2010), and backcross-derived lines (Wassom et al., 2008) (Figure 3, Supplementary Table S2). For instance, qSC-1-1 colocalized with JLM10 (Hu et al., 2021); qSC-1-2 colocalized with qSTA4-3 (Hu et al., 2021); qSC-1-3 colocalized with JLM24 and qSTA8-1 (Hu et al., 2021); qSC-3-1 colocalized with qCT-3-1 (Dong et al., 2015); qSC-4-1 colocalized with qSTA1-1 (Hu et al., 2021) and qzSTA1-1-1 (Yang et al., 2013); and qSC-4-3 colocalized with qCT-5-1 (Dong et al., 2015), stc5 (Guo et al., 2013), qSTA5-3 and qSTA5-5 (Hu et al., 2021), qSTA1-5-1 and qSTA1-5-2 (Wang Y. Z. et al., 2010), and qxSTA2-5-1 (Yang et al., 2013). These results suggest that, although unique QTLs are constantly specialized in each individual population, some genetic loci may have a common effect on SC among different populations.

Figure 3

Discussion

QTL mapping precision

The starch content in maize kernel is a complex quantitative trait. Elucidation of QTLs or genes controlling phenotypic variation of SC could decipher the genetic architecture of starch in maize kernel (Huang et al., 2021). Molecular markers are useful tools to efficiently facilitate the selection process. SNP markers are the most frequent variations in genomes, and the application of SNP markers in plant breeding has guaranteed the precision of QTL mapping and genetic analysis (Bhattramakki et al., 2002; Mammadov et al., 2012; Flutre et al., 2022; Kaur et al., 2022). In the present study, a total of 13 QTLs were distributed on chromosomes 1, 3, 4, 5, 7, 8, and 9. Three QTLs (qSC-2-2, qSC-3-1, and qSC-4-3) spanned a large physical interval (111.41-149.78 Mb), partly due to the limited size of the mapping DH population. The other QTLs spanned physical intervals of <50 Mb, and two spanned <10 Mb.

Genetic basis of SC in our DH populations

The starch content in the four DH populations examined in this study exhibited a broad range of phenotypic variations with normal distribution. The genetic analysis indicated that SC is highly heritable and the heritability is fairly high in all populations, indicating a superior genetic effect on SC in DH populations. In addition, except for environmental variation, none of the QTLs were shared by all DH populations, reflecting the complexity of SC regulation in diverse maize populations. Most of these QTLs had moderate additive effects. The PVE for each QTL ranged from 3.65 (qSC-1-2 in SC1) to 16.18% (qSC-3-2 in SC3). These results suggested that a few large-effect QTLs, together with a large number of minor-effect QTLs, mainly contributed to the genetic component of SC. This is consistent with the quantitative nature of SC, reflecting the complexity of starch biosynthesis and accumulation in maize kernels (Glowinski and Flint-Garcia, 2018).

The results of QTL detection derived from different studies may exhibit consistency to a certain degree across different germplasms/genetic backgrounds and environments. Indeed, QTLs including qSC-2-2, qSC-2-3, qSC-3-2, and qSC-3-3 with the highest contribution to phenotypic variation (10.58–16.18%) displayed a high degree of overlap with regions associated for SC in multiple former studies (Wassom et al., 2008; Zhang et al., 2008; Cook et al., 2012; Yang et al., 2013; Nancy et al., 2014; Dong et al., 2015; Alves et al., 2019; Lin et al., 2019; Hu et al., 2021). Meanwhile, our QTLs with the second largest effect (3.65–8.43%) showed less degree of overlap with other studies (Liu et al., 2008; Li et al., 2009; Wang Y. Z. et al., 2010; Guo et al., 2013). Moreover, qSC-4-2 is a new QTL, which is definitely worth conducting further research on this QTL via near-isogenic lines (NILs), fine mapping, molecular marker-assisted selection (MAS), and ultimate cloning. Furthermore, considering that each of the four DH populations was developed from biparental crosses, the derived QTLs could only explain the variation between the two parents in regulating SC. Therefore, it would be necessary to assess SC in extra germplasms to get a comprehensive insight into the genetic factors regulating the natural variation of SC.

Importance of QTLs relevant to SC in maize genetic and breeding

The starch content in maize kernel is regulated by many genes (Zhong et al., 2020). Outstanding progress has been made in the understanding of the genetic and biochemistry of starch synthesis, which involves the coordinated activities of a series of starch metabolic enzymes, including sucrose synthase (SUS), adenosine 5'diphosphate-glucose pyrophosphorylase (AGPase), starch synthases (SSs), starch branching enzymes (BEs), and starch debranching enzymes (DBEs) (Nelson and Pan, 1995). Candidate genes underlying QTLs associated with SC may be suggested in view of the co-location of QTLs with genes encoding these enzymes (Prioul et al., 1997; Thévenot et al., 2005). Moreover, the co-location analysis could provide information about the functional relationships between gene expression and some QTLs of the starch biosynthesis pathway (Thévenot et al., 2005). Our study involved 471 genes in the glycolysis, sucrose, and starch pathway, of which 123, including 12 well-known genes encoding key enzymes in maize starch regulation of synthesis and metabolism, were located within QTL intervals (Figure 4, Supplementary Table S3).

Figure 4

The starch synthesis pathway in maize endosperm initiates with the cleavage of sucrose into fructose and UDP-glucose, which is catalyzed by SUS (Nelson and Pan, 1995). Sus2 is one of three SUS-encoding genes in maize (Duncan et al., 2006) and is located inside qSC-4-1. Transactivation of Sus1 and Sus2 by Opaque2 is an essential supplement to sucrose synthase-mediated endosperm filling in maize, indicating that Sus2 might have a unique role in cytoplasmic sucrose metabolism (Deng et al., 2020). IVR2 in qSC-4-3 is annotated to be one of the vacuolar invertases in the maize genome (Kim et al., 2000). It is also a key enzyme of carbon metabolism in both source and sink tissues that irreversibly hydrolyze sucrose to fructose and glucose and regulates sugar accumulation in sink organs (Juárez-Colunga et al., 2018). This gene was also identified in a GWAS study and linked with natural variation in the amylose content (Li et al., 2018). Brittle2 (Bt2) in qSC-2-1 encodes the small unit of AGPase, which catalyzes the upstream products converting into adenosine 5'diphosphate-glucose (ADPG), that is, the glucosyl donors of amylose and amylopectin (Preiss et al., 1990) and the first rate-limiting enzyme in the starch biosynthetic pathways (Ma et al., 2021; Finegan et al., 2022). ZMHXK3a in qSC-3-1, PHI1 in qSC-4-2, and TPI3 in qSC-1-3 also play very important roles in glycolysis and starch metabolism. ZMHXK3a encodes a hexokinase, which catalyzes the conversion of glucose to glucose-6-phosphate (Xiao et al., 2016). PHI1 encodes phosphohexose isomerase 1 that converts fructose into glucose-1-phosphate and TPI3 encodes triosephosphate isomerase isozymes representing the cytosol and involved in glycolysis (Wendel et al., 1989). WX1 and SS1 in qSC-3-3 and SS2C in qSC-4-3 belong to starch synthases. WX1 encoding granule-bound starch synthase I (GBSSI) is solely responsible for amylose production (Shure et al., 1983), whereas SS1 (Knight et al., 1998) and SS2C (Yan et al., 2009) encoding soluble starch synthases are responsible for the synthesis of amylopectin. SBE1 in qSC-4-3 encodes an isozyme of SBE that generates amylopectin by cleaving internal amylase a-(1,4) glucosidic bonds and transferring the reducing ends to C6 hydroxyls to form a-(1,6) bonds (Jiang et al., 2013). Besides, SBE1 was shown to be related to amylose content and starch molecular structure (Zhong et al., 2021). SU1 in qSC-2-2 encodes the DBE isoforms ISA1, which is responsible for the ordering of branch linkages and linear chains as depicted in the widely accepted cluster model of the amylopectin structure (James et al., 1995).

Conclusion

In the present study, we interpreted the genetic basis and QTL mapping of SC in four DH populations. A total of 13 QTLs accounted for 12.80–41.57% of the starch variation, with four QTLs explaining that more than 10% of the phenotypic variation were identified. One novel QTL has never been reported in any previous studies. These results indicated that the phenotypic variation in SC is coordinated by large-effect QTLs and minor-effect QTLs. Our results further enhanced the understanding of genetic variation in SC and offered prospective routes to modify SC through molecular marker-assisted selection in the maize breeding program.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 32201798.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Statements

Data availability statement

The original contributions presented in the study are publicly available. This data can be found here: https://doi.org/10.6084/m9.figshare.20188547.v2.

Author contributions

RZ, CZ, and MW conceived and designed the experiments. XZ, CD, and HG performed the research. MW analyzed the data. XZ and MW wrote the manuscript. All authors read and approved the manuscript, contributed to the article, and approved the submitted version.

Acknowledgments

We thank all members of our laboratories for helpful discussion and assistance during this research.

Conflict of interest

Author CZ was employed by Maize Yufeng Biotechnology LLC (Beijing, China). The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.950664/full#supplementary-material

References

  • 1

    AlvesM. L.CarbasB.GasparD.PauloM.BritesC.Mendes-MoreiraP.et al. (2019). Genome-wide association study for kernel composition and flour pasting behavior in wholemeal maize flour. BMC Plant Biol. 19, 123. 10.1186/s12870-019-1729-7

  • 2

    BalterM. (2007). Plant science: starch reveals crop identities. Science.316, 1834. 10.1126/science.316.5833.1834

  • 3

    BhattramakkiD.DolanM.HanafeyM.WinelandR.VaskeD.RegisterJ. C.et al. (2002). Insertion-deletion polymorphisms in 3' regions of maize genes occur frequently and can be used as highly informative genetic markers. Plant Mol. Biol. 48, 539547. 10.1023/a:1014841612043

  • 4

    BordesJ.CharmetG.de VaulxR. D.PollacsekM.BeckertM.GallaisA. (2006). Doubled haploid versus S1 family recurrent selection for testcross performance in a maize population. Theor. Appl. Genet. 112, 10631072. 10.1007/s00122-006-0208-3

  • 5

    BromanK. W.WuH.SenS.ChurchillG. A. (2003). R/qtl: QTL mapping in experimental crosses. Bioinformatics19, 889890. 10.1093/bioinformatics/btg112

  • 6

    ChaikamV.MolenaarW.MelchingerA. E.BoddupalliP. M. (2019). Doubled haploid technology for line development in maize: technical advances and prospects. Theor. Appl. Genet. 132, 32273243. 10.1007/s00122-019-03433-x

  • 7

    ChurchillG. A.DoergeR. W. (1994). Empirical threshold values for quantitative trait mapping. Genetics.138, 963971. 10.1093/genetics/138.3.963

  • 8

    Comparot-MossS.DenyerK. (2009). The evolution of the starch biosynthetic pathway in cereals and other grasses. J. Exp. Bot.60, 24812492. 10.1093/jxb/erp141

  • 9

    CookJ. P.McMullenM. D.HollandJ. B.TianF.BradburyP.Ross-IbarraJ.et al. (2012). Genetic architecture of maize kernel composition in the nested association mapping and inbred association panels. Plant Physiol. 158, 824834. 10.1104/pp.111.185033

  • 10

    DengY. T.WangJ. C.ZhangZ. Y.WuY. R. (2020). Transactivation of Sus1 and Sus2 by Opaque2 is an essential supplement to sucrose synthase-mediated endosperm filling in maize. Plant Biotechnol. J.18, 18971907. 10.1111/pbi.13349

  • 11

    DongY. B.ZhangZ. W.ShiQ. L.WangQ. L.ZhouQ.LiY. L. (2015). QTL identification and meta-analysis for kernel composition traits across three generations in popcorn. Euphytica204, 649660. 10.1007/s10681-015-1360-0

  • 12

    DuncanK. A.HardinS. C.HuberS. C. (2006). The three maize sucrose synthase isoforms differ in distribution, localization, and phosphorylation. Plant Cell Physiol. 47, 959971. 10.1093/pcp/pcj068

  • 13

    FineganC.BoehleinS. K.LeachK. A.MadridG.HannahL. C.KochK. E.et al. (2022). Genetic perturbation of the starch biosynthesis in maize endosperm reveals sugar-responsive gene networks. Front. Plant Sci.12, 800326. 10.3389/fpls.2021.800326

  • 14

    FlutreT.Le CunffL.FodorA.LaunayA.RomieuC.BergerG.et al. (2022). A genome-wide association and prediction study in grapevine deciphers the genetic architecture of multiple traits and identifies genes under many new QTLs. G3 (Bethesda).29, jkac103. 10.1093/g3journal/jkac103

  • 15

    FoiadaF.WestermeierP.KesselB.OuzunovaM.WimmerV.MayerhoferW.et al. (2015). Improving resistance to the European corn borer: a comprehensive study in elite maize using QTL mapping and genome-wide prediction. Theor. Appl. Genet. 128, 875891. 10.1007/s00122-015-2477-1

  • 16

    GallaisA.BordesJ. (2007). The use of doubled haploids in recurrent selection and hybrid development in maize. Crop Sci. 47, S190S201. 10.2135/cropsci2007.04.0019IPBS

  • 17

    GlowinskiA.Flint-GarciaS. (2018). “Germplasm resources for mapping quantitative traits in maize,” inThe Maize Genome. Compendium of Plant Genomes, Tuberosa, R. (ed). New York, NYC: Springer Press. p. 143159.

  • 18

    GuoY. Q.YangX. H.ChanderS.YanJ. B.ZhangJ.SongT. M.et al. (2013). Identification of unconditional and conditional QTL for oil, protein and starch content in maize. Crop J. 1, 3442. 10.1016/j.cj.2013.07.010

  • 19

    GuoZ. F.WangH. W.TaoJ. J.RenY. H.XuC.WuK. S.et al. (2019). Development of multiple snp marker panels affordable to breeders through genotyping by target sequencing (GBTS) in maize. Mol. Breeding.39, 37. 10.1007/s11032-019-0940-4

  • 20

    HuS. T.WangM.ZhangX.ChenW. K.SongX. R.FuX. Y.et al. (2021). Genetic basis of kernel starch content decoded in a maize multi-parent population. Plant Biotechnol. J. 19, 21922205. 10.1111/pbi.13645

  • 21

    HuangL. C.TanH. Y.ZhangC. Q.LiQ. F.LiuQ. Q. (2021). Starch biosynthesis in cereal endosperms: An updated review over the last decade. Plant Commun. 2, 100237. 10.1016/j.xplc.2021.100237

  • 22

    JamesM. G.RobertsonD. S.MyersA. M. (1995). Characterization of the maize gene sugary1, a determinant of starch composition in kernels. Plant Cell.7, 417429. 10.1105/tpc.7.4.417

  • 23

    JeonJ. S.RyooN.HahnT. R.WaliaH.NakamuraY. (2010). Starch biosynthesis in cereal endosperm. Plant Physiol. Biochem.48, 383392. 10.1016/j.plaphy.2010.03.006

  • 24

    JiangL. L.YuX. M.QiX.YuQ.DengS.BaiB.et al. (2013). Multigene engineering of starch biosynthesis in maize endosperm increases the total starch content and the proportion of amylose. Transgenic Res. 22, 11331142. 10.1007/s11248-013-9717-4

  • 25

    Juárez-ColungaS.López-GonzálezC.Morales-ElíasN. C.Massange-SánchezJ. A.TrachselS.TiessenA. (2018). Genome-wide analysis of the invertase gene family from maize. Plant Mol. Biol.97, 385406. 10.1007/s11103-018-0746-5

  • 26

    KaurG.PathakM.SinglaD.ChhabraG.ChhunejaP.Kaur SaraoN. (2022). Quantitative trait loci mapping for earliness, fruit, and seed related traits using high density genotyping-by-sequencing-based genetic map in bitter gourd (Momordica charantia L.). Front. Plant Sci. 12, 799932. 10.3389/fpls.2021.799932

  • 27

    KimJ. Y.Mah,éA.BrangeonJ.PrioulJ. L. (2000). A maize vacuolar invertase, IVR2, is induced by water stress. Organ/tissue specificity and diurnal modulation of expression. Plant Physiol.124, 7184. 10.1104/pp.124.1.71

  • 28

    KnappS. J.StroupW. W.RossW. M. (1985). Exact confidence-intervals for heritability on a progeny mean basis. Crop Sci. 25, 192194. 10.2135/cropsci1985.0011183X002500010046x

  • 29

    KnightM. E.HarnC.LilleyC. E.GuanH.SingletaryG. W.MuForsterC.et al. (1998). Molecular cloning of starch synthase I from maize (W64) endosperm and expression in Escherichia coli. Plant J.14, 613622. 10.1046/j.1365-313x.1998.00150.x

  • 30

    LanderE. S.BotsteinD. (1989). Mapping mendelian factors underlying quantitative traits using RFLP linkage maps. Genetics121, 185199. 10.1093/genetics/121.1.185

  • 31

    LiC. S.HuangY. C.HuangR. D.WuY. R.WangW. Q. (2018). The genetic architecture of amylose biosynthesis in maize kernel. Plant Biotechnol. J.16, 688695. 10.1111/pbi.12821

  • 32

    LiH. H.RibautJ. M.LiZ. L.WangJ. K. (2008). Inclusive composite interval mapping (ICIM) for digenic epistasis of quantitative traits in biparental populations. Theor. Appl. Genet.116, 243260. 10.1007/s00122-007-0663-5

  • 33

    LiY. L.WangY. Z.WeiM. G.LiX. H.FuJ. F. (2009). QTL identification of grain protein concentration and its genetic correlation with starch concentration and grain weight using two populations in maize (Zea mays L.). J. Genet. 88. 6167. 10.1007/s12041-009-0008-z

  • 34

    LinF.ZhouL.HeB.ZhangX. L.DaiH. X.QianY. L.et al. (2019). QTL mapping for maize starch content and candidate gene prediction combined with co-expression network analysis. Theor. Appl. Genet. 132, 19311941. 10.1007/s00122-019-03326-z

  • 35

    LiuY. Y.DongY. B.NiuS. Z.CuiD. Q.WangY. Z.WeiM. G.et al. (2008). QTL identification of kernel composition traits with popcorn using both F2:3 and BC2F2 populations developed from the same cross. J. Cereal Sci. 48, 625631. 10.1016/j.jcs.2008.02.003

  • 36

    MaB.ZhangL.GaoQ. F.WangJ. M.LiX. Y.WangH.et al. (2021). A plasma membrane transporter coordinates phosphate reallocation and grain filling in cereals. Nat. Genet. 53, 906915. 10.1038/s41588-021-00855-6

  • 37

    MammadovJ.AggarwalR.BuyyarapuR.KumpatlaS. (2012). SNP markers and their impact on plant breeding. Int. J. Plant Genomics2012, 728398. 10.1155/2012/728398

  • 38

    MayorP. J.BernardoR. (2009). Genomewide selection and marker-assisted recurrent selection in doubled haploid versus f populations. Crop Sci. 49, 17191725. 10.2135/cropsci2008.10.0587

  • 39

    NancyS. YPineda-HidalgoK. V.Chavez-OntiverosJ.Gutierrez-DoradoR.Reyes-MorenoC.Bello-PérezL. A.et al. (2014). Biochemical characterization of QTLs associated with endosperm modification in quality protein maize. J. Cereal Sci.60, 255263. 10.1016/j.jcs.2014.04.004

  • 40

    NelsonO.PanD. (1995). Starch synthesis in maize endosperms. Annu. Rev. Plant Biol. 46, 475496. 10.1146/annurev.pp.46.060195.002355

  • 41

    OdellS. G.HudsonA. I.PraudS.DubreuilP.TixierM. H.Ross-IbarraJ.et al. (2022). Modeling allelic diversity of multiparent mapping populations affects detection of quantitative trait loci. G3 (Bethesda). 12, jkac011. 10.1093/g3journal/jkac011

  • 42

    PreissJ.DannerS.SummersP. S.MorellM.BartonC. R.YangL.et al. (1990). Molecular characterization of the Brittle-2 gene effect on maize endosperm ADPglucose pyrophosphorylase subunits. Plant Physiol. 92, 881885. 10.1104/pp.92.4.881

  • 43

    PrioulJ. L.QuarrieS.CausseM.de VienneD. (1997). Dissecting complex physiological functions through the use of molecular quantitative genetics. J. Exp. Bot. 48, 11511163. 10.1093/jxb/48.6.1151

  • 44

    ShureM.WesslerS.FedoroffN. (1983). Molecular identification and isolation of the waxy locus in maize. Cell.35, 225233. 10.1016/0092-8674(83)90225-8

  • 45

    ThévenotC.Simond-CôteE.ReyssA.ManicacciD.TrouverieJ.Le GuillouxM.et al. (2005). QTLs for enzyme activities and soluble carbohydrates involved in starch accumulation during grain filling in maize. J. Exp. Bot. 56, 945958. 10.1093/jxb/eri087

  • 46

    WangS.BastenC. J.ZengZ. B. (2010). Windows QTL Cartographer V2.5_011. Raleigh: Dep. Stat. North Carolina State University.

  • 47

    WangT. T.WangM.HuS. T.XiaoY. N.TongH.PanQ. C.et al. (2015). Genetic basis of maize kernel starch content revealed by high-density single nucleotide polymorphism markers in a recombinant inbred line population. BMC Plant Biol. 15, 288. 10.1186/s12870-015-0675-2

  • 48

    WangY. Z.LiJ. Z.LiY. L.WeiM. G.FuJ. F. (2010). QTL detection for grain oil and starch content and their associations in two connected F2:3 populations in high-oil maize. Euphytica.174, 239252. 10.1007/s10681-010-0123-1

  • 49

    WassomJ. J.WongJ. C.MartinezE.KingJ. J.DebaeneJ.HotchkissJ. R.et al. (2008). QTL associated with maize kernel oil, protein, and starch concentrations; kernel mass; and grain yield in illinois high oil × B73 backcross-derived lines. Crop Sci. 48, 243252. 10.2135/cropsci2007.04.0205

  • 50

    WendelJ. F.StuberC. W.GoodmanM. M.BeckettJ. B. (1989). Duplicated plastid and triplicated cytosolic isozymes of triosephosphate isomerase in maize (Zea mays L.). J. Hered. 80, 218228. 10.1093/oxfordjournals.jhered.a110839

  • 51

    XiaoY. N.ThatcherS.WangM.WangT. T.BeattyM.Zastrow-HayesG.et al. (2016). Transcriptome analysis of near-isogenic lines provides molecular insights into starch biosynthesis in maize kernel. J. Integr. Plant Biol.58, 713723. 10.1111/jipb.12455

  • 52

    YanG. J.LiuH.WangH. B.LuZ. Y.WangY. X.MullanD.et al. (2017). Accelerated generation of selfed pure line plants for gene identification and crop breeding. Front. Plant Sci. 8, 1786. 10.3389/fpls.2017.01786

  • 53

    YanH. B.PanX. X.JiangH. W.WuG. J. (2009). Comparison of the starch synthesis genes between maize and rice: copies, chromosome location and expression divergence. Theor. Appl. Genet.119, 815825. 10.1007/s00122-009-1091-5

  • 54

    YangG. H.DongY. B.LiY. L.WangQ. L.ShiQ. L.ZhouQ. (2013). Verification of QTL for grain starch content and its genetic correlation with oil content using two connected RIL populations in high-oil maize. PLoS ONE.8, e53770. 10.1371/journal.pone.0053770

  • 55

    YangN.LuY. L.YangX. H.HuangJ.ZhouY.AliF.et al. (2014). Genome wide association studies using a new nonparametric model reveal the genetic architecture of 17 agronomic traits in an enlarged maize association panel. PLoS Genet. 10, e1004573. 10.1371/journal.pgen.1004573

  • 56

    ZhangJ.LuX. Q.SongX. F.YanJ. B.SongT. M.DaiJ. R.et al. (2008). Mapping quantitative trait loci for oil, starch, and protein concentrations in grain with high-oil maize by SSR markers. Euphytica.162, 335344. 10.1007/s10681-007-9500-9

  • 57

    ZhangX. L.LuM.XiaA. A.XuT.CuiZ. H.ZhangR. Y.et al. (2021). Genetic analysis of three maize husk traits by QTL mapping in a maize-teosinte population. BMC Genomics.22, 386. 10.1186/s12864-021-07723-x

  • 58

    ZhaoH.SunZ. F.WangJ.HuangH. J.KocherJ. P.WangL. G. (2014). CrossMap: a versatile tool for coordinate conversion between genome assemblies. Bioinformatics.30, 10061007. 10.1093/bioinformatics/btt730

  • 59

    ZhongY. Y.LiuL. S.QuJ. Z.LiS. L.BlennowA.SeytahmetovnaS. A.et al. (2020). The relationship between the expression pattern of starch biosynthesis enzymes and molecular structure of high amylose maize starch. Carbohydr. Polym. 247, 116681. 10.1016/j.carbpol.2020.116681

  • 60

    ZhongY. Y.QuJ. Z.BlennowA.LiuX. X.GuoD. W. (2021). Expression pattern of starch biosynthesis genes in relation to the starch molecular structure in high-amylose maize. J. Agric. Food Chem. 69, 28052815. 10.1021/acs.jafc.0c07354

Summary

Keywords

maize, DH, kernel, starch, QTL

Citation

Zhang X, Wang M, Zhang C, Dai C, Guan H and Zhang R (2022) Genetic dissection of QTLs for starch content in four maize DH populations. Front. Plant Sci. 13:950664. doi: 10.3389/fpls.2022.950664

Received

23 May 2022

Accepted

30 June 2022

Published

06 October 2022

Volume

13 - 2022

Edited by

Ting Peng, Henan Agricultural University, China

Reviewed by

Zhanhui Zhang, Henan Agricultural University, China; Yongrui Wu, Institute of Plant Physiology and Ecology, Shanghai Institutes for Biological Sciences (CAS), China

Updates

Copyright

*Correspondence: Ruiying Zhang

†These authors have contributed equally to this work and share first authorship

This article was submitted to Plant Breeding, a section of the journal Frontiers in Plant Science

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics