Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Verification of QTL for Grain Starch Content and Its Genetic Correlation with Oil Content Using Two Connected RIL Populations in High-Oil Maize

  • Guohu Yang,

    Affiliation College of Agriculture, Henan Agricultural University, Key Laboratory of Physiological Ecology and Genetic Improvement of Food Crops in Henan Province, Zhengzhou, People’s Republic of China

  • Yongbin Dong,

    Affiliation College of Agriculture, Henan Agricultural University, Key Laboratory of Physiological Ecology and Genetic Improvement of Food Crops in Henan Province, Zhengzhou, People’s Republic of China

  • Yuling Li ,

    yuling_li@126.com

    Affiliation College of Agriculture, Henan Agricultural University, Key Laboratory of Physiological Ecology and Genetic Improvement of Food Crops in Henan Province, Zhengzhou, People’s Republic of China

  • Qilei Wang,

    Affiliation College of Agriculture, Henan Agricultural University, Key Laboratory of Physiological Ecology and Genetic Improvement of Food Crops in Henan Province, Zhengzhou, People’s Republic of China

  • Qingling Shi,

    Affiliation College of Agriculture, Henan Agricultural University, Key Laboratory of Physiological Ecology and Genetic Improvement of Food Crops in Henan Province, Zhengzhou, People’s Republic of China

  • Qiang Zhou

    Affiliation College of Agriculture, Henan Agricultural University, Key Laboratory of Physiological Ecology and Genetic Improvement of Food Crops in Henan Province, Zhengzhou, People’s Republic of China

Abstract

Grain oil content is negatively correlated with starch content in maize in general. In this study, 282 and 263 recombinant inbred lines (RIL) developed from two crosses between one high-oil maize inbred and two normal dent maize inbreds were evaluated for grain starch content and its correlation with oil content under four environments. Single-trait QTL for starch content in single-population and joint-population analysis, and multiple-trait QTL for both starch and oil content were detected, and compared with the result obtained in the two related F2∶3 populations. Totally, 20 single-population QTL for grain starch content were detected. No QTL was simultaneously detected across all ten cases. QTL at bins 5.03 and 9.03 were all detected in both populations and in 4 and 5 cases, respectively. Only 2 of the 16 joint-population QTL had significant effects in both populations. Three single-population QTL and 8 joint-population QTL at bins 1.03, 1.04–1.05, 3.05, 8.04–8.05, 9.03, and 9.05 could be considered as fine-mapped. Common QTL across F2∶3 and RIL generations were observed at bins 5.04, 8.04 and 8.05 in population 1 (Pop.1), and at bin 5.03 in population 2 (Pop.2). QTL at bins 3.02–3.03, 3.05, 8.04–8.05 and 9.03 should be focused in high-starch maize breeding. In multiple-trait QTL analysis, 17 starch-oil QTL were detected, 10 in Pop.1 and 7 in Pop.2. And 22 single-trait QTL failed to show significance in multiple-trait analysis, 13 QTL for starch content and 9 QTL for oil content. However, QTL at bins 1.03, 6.03–6.04 and 8.03–8.04 might increase grain starch content and/or grain oil content without reduction in another trait. Further research should be conducted to validate the effect of these QTL in the simultaneous improvement of grain starch and oil content in maize.

Introduction

Maize is widely used as an important renewable resource for industrial materials, biodiesel production, and dietary consumption by humans and animals all over the world [1], [2], [3], [4]. Since starch is rich in caloric and oil is high in energy content, both starch and oil content in maize grain have been improved for a long time [5]. After Hopkins (1899) initiated the selection of high and low grain oil content in maize [6], several research went along with the improvement of grain chemical composition using different genetic backgrounds [5], [7], [8], [9], [10]. On the whole, two major practical objectives in maize breeding and genetic research are how to improve grain quality for normal maize and how to improve grain yield for quality maize.

Several early research demonstrated that the negative correlation between grain starch and oil content was resulted in the reduction of starch content and grain yield along with the increase of oil content [4], [11], [12], [13]. Wassom et al. [9] found that starch content was positively correlated with kernel mass in BC1S1s (rp = 0.67**) and with yield in TCs (rp = 0.59**), while oil content was negatively correlated with kernel mass and starch content in BC1S1s (rp = –0.29**, –0.75**) and with yield in TCs (rp = –0.30**, –0.66**) by analyzing kernel traits in a backcross and testcross populations derived from Illinois High Oil (IHO) and the recurrent parent B73. Further negative correlations between starch and oil content at different levels were reflected in populations for QTL mapping using different genetic backgrounds, such as Illinois High Protein (IHP) and Illinois Low Protein (ILP) [5], [14], Illinois High Oil (IHO) and Illinois Low Oil (ILO) [9], [15], [16], [17], Beijing high-oil (BHO) [18], [19], tropical [8], Alexho Single-Kernel (ASK) high-oil [20], [21], and popcorn [13]. Great differences in QTL numbers, locations and effects for the same trait were also observed among those studies. Till now, only one high-oil QTL (qHO6) was applied to marker based selection, refined and cloned [20]. And the increase in seed oil with no change in seed weight was observed in its transgenic lines. Of course, the molecular basis of QTL for oil and starch content and their relationship should be extensively revealed in further research.

In our previous study, two connected F2∶3 populations have been used to identify QTL for grain oil and starch content and their associations [10]. Two connected RIL populations developed from the same two crosses as the F2∶3 populations have been derived and used to detect QTL for ear kernel traits and grain oil content [22]. In this study, the two connected RIL populations were used to detect and compare single-trait QTL for grain starch content in single-population and joint-population analysis, and multiple-trait QTL for grain starch and oil content under 4 environments as in the F2∶3 generations. Our first objective was to analyze the influence of genetic backgrounds, environments and generations on QTL detection for grain starch content. The second objective was to further reveal the genetic mechanism in controlling grain starch content and its correlation with oil content in high-oil maize. The third objective was to find consistent QTL in further research in QTL cloning, marker assisted breeding for starch improvement and the simultaneous improvement of starch and oil content in high-oil maize.

Materials and Methods

Population Development and Field Experiment

The development of the two RIL populations and the field experiment had been described in our previous research [22]. Briefly, one high-oil maize inbred GY220 was crossed with two normal corn inbreds 8984 and 8622 to make two connected RIL populations using single-seed descent method. GY220 was derived from the cycle 27 Alexander high-oil maize population and was selected and provided by China Agricultural University, which belonged to the Lancaster heterotic group [21], [23]. The two normal dent corn inbreds were developed in our laboratory and belonged to the Chinese Reid heterotic group. The RIL populations derived from 8984 × GY220 and 8622 × GY220 were 282 and 263 lines, and referred as population 1 (Pop.1) and population 2 (Pop.2), respectively.

The two RIL populations, along with their two respective parent lines, were planted in two adjacent trails under four environments, at three locations (Xuchang, Xinxiang, and Zhengzhou) in summer sowings (12 June) in Henan, and at Yinchuan in spring sowing (15 April), in China, in 2009. The α-design was used with one-row plots and two replications. Each row was 4 m long with 0.67 m between rows. Plots were planted by hand at a density of 60,000 plants ha−1. Standard cultivation management practices were used at each location [22].

Trait Evaluation

As described by Yang et al. [22], three plants were sib-pollinated within each plot by hand to avoid xenia effect. After maturity, ears of the three plants were harvested and naturally dried. Grain starch and oil content were measured on grain samples mixed within each plot with a MATRIX-1 NIR Spectroscope (Bruker, Corporation, Germany) according to the method described by Dudley and Lambert [24].

Phenotypic Data Analysis

Using the statistical software package SPSS 12.0, combined analysis of variance and correlation coefficients between starch and oil content were carried out following the standard procedures of a mixed model with a random genetic effect, and fixed environment and replicate effects according to the Henderson III method [25]. Since the genotype × environment interactions were significant for grain starch content in the two RIL populations, data under each environment were analyzed individually in further analysis. For comparison, combined analysis was conducted using simple average phenotypic data across four environments as in our previous articles [22]. Broad-sense heritabilities (H2B) for the two connected RIL populations on an entry mean basis were calculated by dividing the genotypic variance by the phenotypic variance [26]. Confidence intervals on heritability estimates were determined according to Knapp et al. [27].

QTL Analysis within Single Population

The integrated genetic map for the two populations included 313 SSR markers, and was 2349.4 cM long with an average interval of 7.50 cM [22]. It was obtained through BioMercator 2.1 software [28].

The composite interval mapping (CIM) [29] with Model 6 of the Zmapqtl procedure in QTL Cartographer Version 2.5 was used in QTL mapping for grain starch content under each environment and combined across four environments in both populations. Thresholds for logarithm of odds (LOD) scores to identify QTL were estimated by permutation tests with a minimum of 1,000 replicates [30]. QTL positions were assigned to relevant regions at the point of a maximum LOD score. Confidence intervals were calculated by subtracting one LOD unit on each side of the maximum LOD position [31]. Multiple interval mapping (MIM) in QTL Cartographer Version 2.5 was used to analyze interactions between detected QTL for grain starch content [32], [33]. Significant thresholds were identified by the quick method for computing approximate thresholds for QTL detection [34].

To further test the genetic relationships between grain oil and starch content, a multiple-trait version of the composite interval mapping was used detect multiple-trait QTL [35] with QTL Cartographer Version 2.5 [33]. A significance threshold was also identified by the quick method for computing approximate thresholds for quantitative trait loci detection [34].

Joint-population QTL Analysis for the Two Population

To combine the two populations, a joint-population analysis was performed using the joint inclusive composite interval mapping (JICIM) method [36]. QTL mapping for the nested association mapping (NAM) design in QTL IciMapping Version 3.2 (available from www.isbreeding.net) was used to detect joint-population QTL for grain starch content under each environment and combined across four environments. Empirical LOD score of 2.5 were used to declare the existence of QTL. The additive effects, positions, LOD values, and the phenotypic variance explained (R2) for each detected QTL were obtained.

For each detected QTL, positive and negative additive effects indicated that the allele from the high-oil parent GY220 and the dent corn inbred lines 8984 or 8622 increased the value of the trait, respectively. QTL were named according to “q”+“environment abbreviation”+“trait abbreviation”+“population number”+“–”+“chromosome number”+“–”+“QTL number”. For joint-population QTL, the “population number” was omitted. Multiple-trait QTL for the two traits were named according to “q”+“environment abbreviation”+“abbreviation for the two traits (starch and oil, SO)”+“population number”+“–”+“chromosome number”+“–”+“QTL number”. For example, the front words “qc”, “qx”, “qz”, “qy” and “q” represented QTL detected at Xuchang, Xinxiang, Zhengzhou, Yinchuan and in combined analysis, respectively.

Results

Combined Variance, Heritability and Performance for Grain Starch Content in the Two Connected RIL Populations

The result of combined analysis of variance showed that all variance components were significant for grain starch content in both populations except environment variance (σE2) in Pop.1 (Table 1). The broad sense heritability (hB2) estimates were a little higher in Pop.1 (0.72) than in Pop.2 (0.69). Among the three parent lines, the value for grain starch content for the high-oil maize inbred GY220 was lower than that of the two normal inbreds 8984/8622 under all environments. But the difference between the two normal inbreds was not obvious (Table 1). In the two populations, all values of grain starch content showed normal distribution with a wide range of variation and transgressive segregation exceeding both parent values. The variance coefficients (CV%) were higher in Pop.2 than in Pop.1 under all environments.

thumbnail
Table 1. Combined variance analysis, heritability and phenotypic performance of grain starch content for the two RIL populations under each environment.

https://doi.org/10.1371/journal.pone.0053770.t001

Significantly negative phenotypic and genotypic correlations were observed between grain starch and oil content under four environments and in combined analysis in both connected populations, −0.21**–−0.35** in Pop.1, and −0.26**–−0.40** in Pop.2 (Table 2).

thumbnail
Table 2. Phenotypic and genotypic correlations between grain starch and oil content for the two RIL populations under each environment and in combined analysis.

https://doi.org/10.1371/journal.pone.0053770.t002

QTL Analysis for Grain Starch Content within Single Population

Since the variances of genotype × environment interaction (σGE2) for grain starch content were significant in both populations, QTL mapping was conducted under each environment. For comparison, combined analyses using means across the four environments were also conducted. A total of 20 single-population QTL were detected under four environments and in combined analyses in the two populations (Table 3, Fig. 1), 12 in Pop.1 and 8 in Pop.2. Except no QTL were detected for Pop.2 at Yinchuan, 1–4 QTL were detected in each case.

thumbnail
Figure 1. The integrated genetic map, single-trait QTL for grain starch content in single-population and joint-population analysis, and multiple-trait QTL for starch and oil content on chromosomes.

QTL are showed to the left of chromosome. One-LOD support intervals are indicated by vertical bars, and the maximum LOD peak positions are indicated by horizontal bars.

https://doi.org/10.1371/journal.pone.0053770.g001

thumbnail
Table 3. QTL detected for grain starch content in the two RIL populations under each environment and in combined analysis.

https://doi.org/10.1371/journal.pone.0053770.t003

These QTL were located on chromosomes 1, 3, 5, 8, 9 and 10. The QTL on chromosomes 3, 5 and 9 were all detected in both populations. The QTL on chromosome 9 was located at the same or near bins 9.03–9.04 in 6 cases, 2 cases in Pop.1 and 4 cases in Pop.2. The QTL on chromosome 5 at bins 5.02–5.04/5.03–5.05 were detected in 2 cases for both populations. The QTL on chromosome 3 were detected in 3 cases, 2 cases in Pop.1 and 1 case in Pop.2. But they were located at different bins, at bins 3.07–3.08/3.05 in Pop.1 and at bins 3.02–3.03 in Pop.2. The QTL on chromosomes 1 and 8 were only detected in Pop.1 in 2 and 4 cases, respectively. One QTL on chromosome 10 at bins 10.05–10.06 was detected in Pop.2. Except five QTL at bins 1.03, 3.02–3.03, 8.03–8.04 and 10.05, the positive alleles of other 14 QTL were all contributed by the normal corn parents. The contribution to phenotypic variation for qcSTA2-3-1 was the largest (12.35%), which could be considered as a major QTL. But those for other 19 individual QTL were all smaller than 10% (between 3.24% and 8.03%), they were all minor QTL. The marker intervals of QTL positions were 2.4–33.55 cM, with an average of 12.59 cM. The QTL numbers with marker intervals <5 cM, 5–10 cM, 10–15 cM and >15 cM were 3, 5, 6 and 6, accounting for 15.0%, 25.0%, 30.0% and 30.0%, respectively. Three QTL (qxSTA1-8-1, qySTA1-3-1 and qySTA1-8-1) with intervals 3.9 cM, 2.4 cM and 3.9 cM could be considered as being fine-mapped.

Among detected QTL, 9 pairs of digenic interactions were identified under each environment and in combined analysis in the two populations (Table 4), which included 3 QTL×QTL (additive×additive effects), 4 QTL×genetic background (marker intervals on which none QTL were detected) and 1 genetic background×genetic background interactions. Those interactions were related with all the six chromosomes with QTL distributed. Besides, 5 marker intervals on chromosomes 4 and 7 were also included, on which QTL were failed to be detected. All the values of interaction effects were low, from 0.4%–4.3%. This result suggested that the contributions of digenic interactions to the performance of grain starch content were all minimal in both populations. Obvious influence of environments on epistasis could also be observed. For example, although both QTL on chromosomes 8 and 9 (at the same marker interval umc1688–umc1771) were detected at Xuchang and Xinxiang in Pop.1, different digenic interactions were detected, umc1847–bnlg2244 ×qcSTA1-9-1 at Xuchang vs. qxSTA1-8-1×qxSTA1-9-1 at Xinxiang.

thumbnail
Table 4. Digenic epistatic interactions among detected QTL for grain starch content under each environment and in combined analysis for the two RIL populations.

https://doi.org/10.1371/journal.pone.0053770.t004

Joint-population QTL Analysis for Grain Starch Content Across the Two Connected RIL Populations

Using joint-population analysis across the two connected RIL populations by JICIM, 16 joint-population QTL for grain starch content were detected under each environment and in combined analysis, 3 at Xuchang, 5 at Xinxiang, 3 at Zhengzhou and 5 in combined analysis (Table 5, Fig. 1). No QTL was detected at Yinchuan. The contribution to phenotypic variation for an individual QTL varied between 11.39% and 23.45%, with 2.00%–95.82% in Pop.1 and 1.29%–26.33% in Pop.2. The marker intervals of QTL positions were 3.2–36.4 cM, with an average of 9.59 cM. The QTL numbers with marker intervals <5 cM, 5–10 cM, 10–15 cM and >15 cM were 8, 2, 5 and 1, accounting for 50.0%, 12.5%, 31.25% and 6.25%, respectively. Eight QTL (qcSTA8-1, qcSTA9-1, qxSTA1-1, qxSTA9-1, qzSTA8-1, qSTA1-1, qSTA8-1 and qSTA9-2) at bins 1.03, 1.04–1.05, 8.06, 9.03 and 9.05 with intervals 3.6 cM, 4.1 cM, 4.8 cM, 3.9 cM and 3.2 cM could be considered as being fine-mapped.

thumbnail
Table 5. QTL for grain starch content detected by joint-population analysis for the two population using JICIM.

https://doi.org/10.1371/journal.pone.0053770.t005

In comparison, 12 of 16 joint-population QTL were detected in the single-population analysis (Table 3). Although single-population QTL qSTA2-10-1 failed to be detected herein, 4 joint-population QTL (qxSTA1-1, qxSTA7-1, qSTA8-1 and qSTA9-1) could be considered as additional QTL detected in joint-population analysis. The LOD values and the contribution to phenotypic variations for most joint-population QTL were much higher than those in the single-population mapping, while their interval distances were much smaller. This reflected the higher mapping power of joint-population analysis. Ten of the 16 joint-population QTL with significant effects in only one population could be considered as rare QTL. However, 2 joint-population QTL (qcSTA9-1 and qzSTA5-1) had significant genetic effects in both populations.

Multiple-trait QTL Analysis for Grain Starch Content with Grain Oil Content in the Two Connected RIL Populations

To further analyze the genetic correlations between starch and oil content, multiple-trait analysis for grain starch and oil content was conducted in both populations (Table 6; Fig. 1). Totally, 17 starch-oil QTL (qSO) were detected, including 3 at Xuchuang, 3 at Xinxiang, 3 at Zhengzhou and 8 in combined analysis. No starch-oil QTL were detected at Yinchuan in both populations. Ten and 7 starch-oil QTL were detected in Pop.1 and Pop.2, respectively. These QTL were located on chromosomes 1, 3, 4, 5 and 8 at bins 1.03, 3.05/3.01–3.03, 4.05–4.06, 5.05–5.06/5.06–5.07, 8.04–8.05, respectively.

thumbnail
Table 6. Multiple-trait QTL analysis of grain starch content with oil content in the two RIL populations under each environment and in combined analysis.

https://doi.org/10.1371/journal.pone.0053770.t006

In our previous research, 12 and 14 QTL for grain oil content were detected in the two populations, respectively [22]. They were located on chromosomes 1, 3, 4, 5, 6, 8 and 10, and related with 14 bin loci (1.03, 1.03–1.04, 1.05–1.06, 1.08, 3.01–3.03, 3.03–3.04, 3.06–3.08, 4.05–4.06, 5.05–5.06, 6.03–6.04, 6.04, 8.03, 8.06 and 10.04–10.05). Compared with single-trait QTL for both traits, no additional QTL was detected, which indicated QTL only detected in multiple -trait analysis. The graph of LOD curve peaks for grain starch and oil content changed simultaneously and in the same direction at 5 marker intervals in 8 cases, including phi053–umc1174 on chromosome 3 at Zhengzhou, umc1548–umc1329 on chromosome 4 in combined analysis in Pop.1, and umc2316–umc1979 on chromosome 6 at Xinxiang, at Zhengzhou and in combined analysis, phi100175–bnlg1863 on chromosome 8 at Zhengzhou and in combined analysis, and umc1601–bnlg1615 on chromosome 1 in combined analysis in Pop.2. These data suggested that there existed pleiotropic QTL controlling both traits simultaneously. Besides, the peaks of LOD curve graph for both traits changed in the same close direction at 4 marker intervals in 7 cases, including umc2026–umc1019 and umc2305–bnlg1306 on chromosome 5 at Xuchang and in combined analysis, and umc1044–phi109275 on chromosome 1 at Xinxiang and in combined analysis in Pop.1, and phi031–bnlg1617 on chromosome 6 in Pop.2. These results suggested that there were tightly linked QTL controlling grain starch and oil content in these marker intervals.

However, 22 single-trait QTL failed to show significance in multiple-trait analysis, including 6 QTL for starch content at bins 3.07–3.08, 5.02–5.04, 8.06 and 9.03–9.04, and 3 QTL for oil content at bins 6.04 and 10.04–10.05 in Pop.1, and 7 QTL for starch content at bins 3.07–3.08, 5.03/5.03–5.05, 9.03/9.05 and 10.05–10.06, and 6 QTL for oil content at bins 1.03–1.04/1.05–1.06/1.08, 3.06/3.06–3.08 and 8.06 in Pop.2 (Table 3, Table 6). Those QTL might have effects with opposite directions for both traits. QTL at bins 1.03, 3.05, 8.04-8.05 in Pop.1, and at bin 3.03 in Pop.2, were all detected both in single-trait QTL mapping for starch content and in multiple-trait QTL mapping for both traits. These QTL for starch content might not be influenced by oil content.

Discussion

Comparison of QTL Detected for Grain Starch Content in the Two Connected RIL Populations and in Previous Studies

In this study, 20 single-population and 16 joint-population QTL for grain starch content were identified for the two connected RIL populations under four environments and in combined analysis. In single-population analysis, no QTL showed consistency across all cases. QTL at the same bins 5.03 and 9.03 were detected in both populations in 4 and 5 cases, respectively. Three QTL on chromosome 3 were located at different bins in both populations, at bins 3.05 and 3.07–3.08 in Pop.1 and at bins 3.02–3.03 in Pop.2. The QTL on chromosomes 1 and 8 in Pop.1, and the QTL on chromosome 10 in Pop.2 could be considered as population-specific QTL. Joint-population analysis showed that only 2 QTL (qcSTA9-1 and qzSTA5-1) had significant effects in both populations. Ten of the 16 joint-population QTL had significant effects in only one population, which could be considered as rare QTL.

In our previous study, two F2∶3 populations developed from the same two crosses as in this study had been used to map QTL for grain starch content [10]. In comparison, QTL linked with the same markers at bins 5.04 (umc1162), 8.04 (bnlg2046) and 8.05 (umc1562) in Pop.1, and at bin 5.03 (umc1389) in Pop.2 were observed across both F2∶3 and RIL generations. However, QTL showing generation-specific were observed in both populations, which included QTL on chromosomes 6 and 10 in F2∶3, and on chromosomes 1, 3 and 9 in RIL for Pop.1, and on chromosomes 1, 2, 4 and 6 in the F2∶3 generation, and on chromosomes 3, 9 and 10 in the RIL generations for Pop.2. Great influences of genetic backgrounds (different populations and different generations derived from the same cross) and environments on QTL detection for grain starch content had been commonly reported in previous research [10], [19], [20]. Since the field experiments for the F2∶3 and RIL populations were conducted in different years and at different locations, the discrepancies of QTL detected in the two generations could be brought by different genetic backgrounds, environments and their interactions. For the dramatically smaller percentage of variation explained by each QTL at chromosomes 5 and 8 herein, another reason might be the dominance effects at different levels (partial dominance, dominance, and over-dominance) for most QTL detected in the F2∶3 generation. Therefore, it was necessary to conduct QTL mapping using various genetic backgrounds under diverse environments thoroughly. Only such consistent QTL would have broad values in theoretical research and practical breeding.

Simultaneously considering the results in other previous research [5], [9], [13], [19], [21], [37], QTL detected herein at all bin loci except bin 8.06 have been detected in 1–4 previous reports (Table 7). qcSTA2-3-1 located at bins 3.02–3.03 with the highest contribution to phenotypic variation (12.35%) has also been identified by Liu et al. [13] and Wassom et al. [9]. The QTL at bin 5.03 with consistency across populations (Pop.1 and Pop.2), generations (RIL and F2∶3) and environments has also been identified by Goldman et al. [5], Zhang et al. [19], Wassom et al. [9] and Wang et al. [10]. The QTL at bins 8.04-8.05 with consistency across generations and environments in Pop.1 has also been identified by Zhang et al. [19], Wassom et al. [9] and Wang et al. [10]. The QTL at bin 9.03 with consistency across populations and environments has also been identified by Goldman et al. [5] and Wassom et al. [9].

thumbnail
Table 7. Comparison of QTL for grain starch content detected in the two RIL populations and in previous studies.

https://doi.org/10.1371/journal.pone.0053770.t007

The biosynthesis process of starch is very complex and is regulated by many genes. QTL co-located with genes encoding for related functions could be considered as candidate genes [40], [41]. ADP-Glc pyrophosphorylase, starch synthase (granule-bound starch synthase and soluble starch synthases), starch branching enzyme and starch de-branching enzyme are four kinds of important enzymes in starch synthesis and accumulation [42], [43], [44]. According to the public genetic linkage map (www.maizegdb.org), several related functional genes were located in the marker intervals of QTL on chromosomes 1, 3, 5 and 8, such as collapsed kernel3 (p3), defective kernel1(dek1), dek32, pyruvate decarboxylase3 (pdc3), soft endosperm3(sen3), NADH ubiquinone oxidoreductase1(nad1) on chromosome 1, beta-D-glucosidase (glu1), NADP malic enzyme3 (me3), 6-phosphogluconate dehydrogenase (pdh1), and sen1 on chromosome 3, acetolactate synthase 2 (als2), dek27 and sucrose export defective1 (sxd1) on chromosome 5, and sugar transport 1 (stp1), starch branching enzyme 3 (sbe3), and sucrose phosphate synthase1 (sps1) on chromosome 8. Further studies should be concentrated on the chromosome intervals of QTL at bins 3.02–3.03, 3.05, 5.03, 8.04–8.05 and 9.03 in near isogenic line (NIL) construction, fine-mapping, and marker-assisted selection (MAS) in starch improvement. In fact, three QTL (qxSTA1-8-1, qySTA1-3-1 and qySTA1-8-1) in the single-population analysis and eight QTL (qcSTA8-1, qcSTA9-1, qxSTA1-1, qxSTA9-1, qzSTA8-1, qaSTA1-1, qaSTA8-1 and qaSTA9-1) in the joint-population analysis at bins 1.03, 1.04–1.05, 3.05, 8.04–8.05, 9.03, and 9.05 could be considered as being fine-mapped. MAS has been used to select their NILs in our present research.

Comparison of QTL Detected for Starch and Oil Content and their Correlations

In maize kernel, 98% of starch is stored in endosperm and 85% of oil is located in embryo [38]. Selection for oil increase in kernel was accompanied by decrease in kernel weight in IHO [9], [12], a synthetic population [39], and normal corn germplasm [45]. Negative relationship between grain oil content and starch content was commonly reported in research using Illinois strains [4,9,46], BHO [19], ASK high-oil [10], and popcorn germplasm [13]. In the present study, significant negative phenotypic and genotypic correlations between grain starch and oil content were observed for the two RIL populations under each environment and in combined analysis. Comparing with the single-trait QTL detected for starch and oil content, QTL for both traits were all detected on chromosomes 1, 3 and 8 in Pop.1, and on chromosome 3 in Pop.2. Only the QTL at bin 1.03 for both traits was located in the same marker interval (umc1044–phi109275), and with the favorable alleles contributed by the high-oil parent GY220. Specific QTL were observed for grain oil content on chromosome 6 and for grain starch content on chromosome 9. However, tightly linked and/or pleiotropic QTL for both traits were reflected in multiple-trait analysis at bins 1.03, 3.05, 4.05–4.06, 5.05–5.06 and 5.06–5.07 in Pop.1, and at bins 1.05–1.06, 6.03–6.04, 6.04–6.05 and 8.03 in Pop.2. In addition, most single-trait QTL for oil and starch content failed to be detected in multiple-trait analysis. Similar results have been observed by Berke and Rocheford [15], Zhang et al. [19] and Wang et al. [10]. Obviously, the relationship of grain starch content and oil content was very complicated [10], [18]. This reflected that tightly linked and/or pleiotropic QTL with opposite effects might exist at those chromosome regions. It was not easy to improve grain starch and oil content simultaneously.

However, QTL for grain starch and oil content were all detected with all the favorable alleles contributed by the high-oil parent GY220 at bins 1.03 and 8.03–8.04. There might exist tightly linked and/or pleiotropic QTL with the same direction effect for both traits. QTL located at bin 6.03–6.04 for grain oil content might increase grain oil content without reduction in grain starch content. In fact, the high-oil QTL (qHO6) cloned by Zheng et al. [20] could increase grain oil content with no effect on seed weight, which was located at bin 6.04. Some candidate genes encoding the enzymes involved in the synthesis and modification pathways of grain oil were located in the marker intervals of QTL on chromosomes 1 and 6, such as fatty acid desaturase8 (fad8) and oleic acid content1 (olc1) on chromosome 1, and linoleic-acid content (ln1) on chromosome 6. Zheng et al. [20] considered that ln1and qHO6 were likely the same genes. Further research should be conducted to prove the effect of these QTL in the simultaneous improvement of grain starch and oil content in maize. Certainly, what is the real situation in our populations should be proved in further research.

Acknowledgments

We greatly thank China Agricultural University for providing us the high-oil maize inbred line GY220. We are grateful to Jiankang Wang and Huihui Li in helping us for joint-population QTL analysis.

Author Contributions

Conceived and designed the experiments: YL. Performed the experiments: GY QW QS QZ. Analyzed the data: GY YD. Wrote the paper: GY YL.

References

  1. 1. Adams KL, Jensen AH (1988) High-fat maize in diets for pigs and sows. Anim Feed Sci Technol 17: 201–212.
  2. 2. Goss J, Kerr PS (1992) Challenges and opportunities for identity preserved varieties. In: Proc 47th Annu Corn and Sorghum Res Conf, American Seed Trade Association, Washington, DC, pp82–92.
  3. 3. Han Y, Parsons CM, Alexander DE (1987) Nutritive value of high oil corn for poultry. Poultry Sci 66: 103–111.
  4. 4. Lambert RJ (2001) High-oil corn hybrids. In: Hallauer AR (ed) Specialty corns. CRC, Boca Raton, FL, 131–154.
  5. 5. Goldman IL, Rocheford TR, Dudley JW (1994) Molecular marker associated with maize kernel oil concentration in an Illinois high protein and Illinois low protein cross. Crop Sci 34: 908–915.
  6. 6. Hopkins CG (1899) Improvement in the chemical composition of the corn kernel. Illinois Agr Exp Sta Bul 55: 205–240.5.
  7. 7. Alexander DE, Lambert RJ (1968) Relationship of kernel oil content to yield in maize. Crop Sci 8: 272–274.
  8. 8. Mangolin CA, de Souza Jr CJ, Garcia AAF, Garcia AF, Sibov ST, et al. (2004) Mapping QTLs for kernel oil content in a tropical maize population. Euphytica 137: 251–259.
  9. 9. Wassom JJ, Wong JC, Martinez E, King JJ, DeBaene J, 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: 243–252.
  10. 10. Wang YZ, Li JZ, Li YL, Wei MG, Li XH, et al. (2010) QTL detection for grain oil and starch content and their associations in two connected F2:3 populations in high-oil maize. Euphytica 174: 239–252.
  11. 11. Woodworth CM, Leng ER, Jugenheimer RW (1952) Fifty generations of selection for protein and oil content in corn. Agron J 44: 60–65.
  12. 12. Dudley JW, Lambert RJ, Alexander DE (1974) Seventy generations of selection for oil and protein concentration in the corn kernel. in Dudley JW (Ed) Seventy Generations for Oil and Protein in Maize. ASA, CSSA, Madison, W. I: 181–212.
  13. 13. Liu YY, Dong YB, Niu SZ, Cui DQ, Wang YZ, et al. (2008) QTL identification of kernel quality traits with popcorn using both F2:3 and BC2F2 populations developed from the same cross. J Cereal Sci 48: 625–631.
  14. 14. Dudley JW, Dijkhuizen A, Paul C, Coates ST, Rocheford TR (2004) Effects of random mating on marker-QTL associations in the cross of the Illinois High Protein-Illinois Low Protein maize strains. Crop Sci 44: 1419–1428.
  15. 15. Berke TG, Rocheford TR (1995) Quantitative trait locus for flowering, plant, and ear height, and kernel traits in maize. Crop Sci 35: 1542–1549.
  16. 16. Laurie CC, Chasalow SD, Ledeaux JR, McCarroll R, Bush D, et al. (2004) The genetic architecture of response to long-term artificial selection for oil concentration in the maize kernel. Genetics 168: 2141–2155.
  17. 17. Willmot DB, Dudley JW, Rocheford TR, Bari A (2006) Effect of random mating on marker-QTL associations for grain quality traits in the cross of Illinois High Oil×Illinois Low Oil. Maydica 51: 187–199.
  18. 18. Song XF, Song TM, Dai JR, Rocheford, Li JS (2004) QTL mapping of kernel oil concentration with high-oil maize by SSR markers. Maydica 49: 41–48.
  19. 19. Zhang J, Lu XQ, Song XF, Yan JB, Song TM, et al. (2008) Mapping quantitative trait locus for oil, starch, and protein concentrations in grain with high-oil maize by SSR markers. Euphytica 162: 335–344.
  20. 20. Zheng PZ, Allen WB, Roesler K, Williams ME, Zhang SR, et al. (2008) A phenylalanine in DGAT is a key determinant of oil content and composition in maize. Nat Genet 40: 367–372.
  21. 21. Li YL, Wang YZ, Wei MG, Li XH, Fu JF (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.). Journal of Genetics 88: 61–67.
  22. 22. Yang GH, Li YL, Wang QL, Zhou YG, Zhou Q, et al. (2012) Detection and integration of quantitative trait loci for grain yield components and oil content in two connected recombinant inbred line populations of high-oil maize. Molecular Breeding 29: 313–333.
  23. 23. Jiang HY, Chen SJ, Gao LF, Xing JM, Song TM, et al. (2005) Analysis on heterotic groups and heterotic patterns of high oil corn inbred lines. Acta Agron Sin 31: 361–367.
  24. 24. Dudley JW, Lambert RJ (1992) Ninety generations of selection for oil and protein in maize. Maydica 37: 81–87.
  25. 25. Henderson CR (1953) Estimation of variance and covariance. Biometrics 9: 226–252.
  26. 26. Hallauer AR, Miranda JB (1981) Quantitative genetics in maize breeding. Iowa State University Press, Ames.
  27. 27. Knapp SJ, Stroup WW, Ross WM (1985) Exact confidence intervals for heritability on progeny mean basis. Crop Sci 25: 192–194.
  28. 28. Arcade A, Labourdette A, Falque M, Mangin B, Chardon F, et al. (2004) BioMercator: integrating genetic maps and QTL towards discovery of candidate genes. Bioinformatics 20: 2324–2326.
  29. 29. Zeng ZB (1994) Precision mapping of quantitative trait locus. Genetics 136: 1457–146830.
  30. 30. Churchill GA, Doerge RW (1994) Empirical threshold values for quantitative trait mapping. Genetics 138: 963–971.
  31. 31. Lynch M, Walsh B (1998) Genetics and analysis of quantitative traits. Sinauer Associates Inc, Sunderland, Massachusetts, USA, 448.
  32. 32. Kao CH, Zeng ZB, Robert DT (1999) Multiple interval mapping for quantitative trait loci. Genetics 152: 1203–1216.
  33. 33. Wang S, Basten CJ, Zeng ZB (2006) Windows QTL Cartographer 2.5. Department of Statistics, North Carolina State University, Raleigh, NC. Available: http://statgen.ncsu.edu/qtlcart/WQTLCart.htm.
  34. 34. Piepho HP (2001) A quick method for computing approximate thresholds for quantitative trait loci detection. Genetics 157: 425–432.
  35. 35. Jiang C, Zeng ZB (1995) Multiple trait analysis of genetic mapping for quantitative trait loci. Genetics 140: 1111–1127.
  36. 36. Li HH, Bradbury P, Ersoz E, Buckler ES, Wang JK (2011) Joint QTL linkage mapping for multiple-cross mating design sharing one common parent. PloS ONE 6: e17573.
  37. 37. Austin DF, Lee M (1996) Comparative mapping in F2:3 and F6:7 generations of quantitative trait locus for grain yield and yield composition traits in maize. Theor Appl Genet 92: 817–826.
  38. 38. Watson SA (1987) Structure and composition. In Corn: Chemistry and Technology (Edited by Watson SA and Ramstad PE), 53–82. Americtn Associ. Cereal Chemists. Inc. USA.
  39. 39. Mišević D, Alexander DE (1989) Twenty-four cycles of phenotypic recurrent selection for percent oil in maize. I Perse and test-cross performance. Crop Sci 29: 320–324.
  40. 40. Prioul JL, Quarrie S, Causse M, de Vienne D (1997) Dissecting complex physiological functions through the use of molecular quantitative genetics. J Exp Bot 48: 1151–1163.
  41. 41. Thévenot C, Simond-Côte E, Reyss A, Manicacci D, Trouverie J, et al. (2005) QTLs for enzyme activities and soluble carbohydrates involved in starch accumulation during grain filling in maize. J Exp Bot 56: 945–958.
  42. 42. Prioul JL, Reyss A, Schwebel-Duguè N (1990) Relationships between carbohydrate metabolism in ear and in the adjacent leaf during grain filling in maize genotypes. Plant Physiol Bioch 28: 489–493.
  43. 43. Prioul JL, Jeannette E, Reyss A, Grégory N, Giroux M, et al. (1994) Expression of ADP-glucose pyrophosphorylase in maize (Zea mays L.) grain and source leaf during grain filling. Plant Physiol 104: 179–187.
  44. 44. Slattery CJ, Kavakli IH, Okita TW (2000) Engineering starch for increased quantity. Trends in Plant Sci 5: 291–298.45. Miller RL, Dudley JW, Alexander DE (1981) High intensity selection for per cent oil in corn. Crop Sci 21: 433–437.
  45. 45. Lambert RJ, Alexander DE, Mollring EL, Wiggens B (1997) Selection for increased oil concentration in maize kernels and associated changes in several kernel traits. Maydica 42: 39–43.