- Research article
- Open Access
Ploidy effect and genetic architecture exploration of stalk traits using DH and its corresponding haploid populations in maize
BMC Plant Biologyvolume 16, Article number: 50 (2016)
Doubled haploid (DH) lines produced via in vivo haploid induction have become indispensable in maize research and practical breeding, so it is important to understand traits characteristics in DH and its corresponding haploids which derived from each DH lines. In this study, a DH population derived from Zheng58 × Chang7-2 and a haploid population, were developed, genotyped and evaluated to investigate genetic architecture of eight stalk traits, especially rind penetrometer resistance (RPR) and in vitro dry matter digestion (IVDMD), which affecting maize stalk lodging-resistance and feeding values, respectively.
Phenotypic correlation coefficients ranged from 0.38 to 0.69 between the two populations for eight stalk traits. Heritability values of all stalk traits ranged from 0.49 to 0.81 in the DH population, and 0.58 to 0.89 in the haploid population. Quantitative trait loci (QTL) mapping study showed that a total of 47 QTL for all traits accounting for genetic variations ranging from 1.6 to 36.5 % were detected in two populations. One or more QTL sharing common region for each trait were detected between two different ploidy populations. Potential candidate genes predicated from the four QTL support intervals for RPR and IVDMD were indirectly or directly involved with cellulose and lignin biosynthesis, which participated in cell wall formation. The increased expression levels of lignin and cellulose synthesis key genes in the haploid situation illustrated that dosage compensation may account for genome dosage effect in our study.
The current investigation extended understanding about the genetic basis of stalk traits and correlations between DH and its haploid populations, which showed consistence and difference between them in phenotype, QTL characters, and gene expression. The higher heritabilities and partly higher QTL detection power were presented in haploid population than in DH population. All of which described above could lay a preliminary foundation for genetic architecture study with haploid population and may benefit selection in haploid-stage to reduce cost in DH breeding.
Maize (Zea mays L.) is one of the important grain and feed crops in which the stalk, as one indispensable part of plant morphology, serves as the conductor of transporting water and nutrients. Stalk lodging lead to yield losses estimated to range from 5 to 20 % annually worldwide . Rind penetrometer resistance (RPR), which is one of the reliable indicators of stalk strength, has been widely used to measure stalk strength and improve stalk lodging resistance [2, 3]. Maize is also one of the most important annual forage crops. In vitro dry matter digestion (IVDMD) has been the most useful evaluating indicator for maize forage variety examination in many European countries . Therefore, a further and better understanding of the molecular basis for RPR and IVDMD is crucial for breeding lodging-resistant and highly digestible maize .
The genetic analysis of quantitative traits is difficult and complex in maize, and quantitative traits are affected by key genes and interacting networks of small-effect genes. Therefore, different studies have provided different results including quantitative trait loci (QTL) number, distribution, and genetic effects for one trait [6, 7]. This lack of conformity may also be explained by the many differences in parental materials, segregation-population types, ecological conditions, genetic maps, analytical methods and phenotype evaluation [8, 9]. Moreover, high genome dosage levels have effect on genetic analysis [10, 11].
Due to the advantages of time-saving and high genetic variance, doubled haploid (DH) technology is routinely used in modern maize breeding for production of homozygous parental lines for maize hybrid breeding and constructing DH populations for genetic research [12–15]. Although haploid populations possess the characteristics of genetic homozygosity and have one genome dosage, moderate to strong correlations have been identified between small size DH populations and their haploid version populations for some agronomic traits . Moreover, haploid lines could react more sensitively to biotic and abiotic stresses and, therefore, they would effectively uncover susceptibility to diseases and environmental constraints. In A. thaliana, the utility and power of haploid genetics had been reported. Haploids can provide genetic analysis advantages that are not available in diploids, such as specifically pyramiding multiple mutant combinations, forward mutagenesis screens and swapping of nuclear and cytoplasmic genomes . In yeast, haploid screens represent an ideal platform for negative selection since a certain genetic lesion set by mutagenesis will exert equal effects in all cells . In this regard, the haploid lines may also be interesting in the genetic architecture exploration of maize quantitative traits.
Different segregating populations have been used in linkage analysis or genome-wide association study of RPR, and the genome set number of all these populations was two. The results suggested the genetic complexity of RPR. Flint Garcia et al.  first detected 35 RPR QTL in four F2:3 populations, which accounted for more than 33 % of the total phenotype variation. Hu et al.  detected 9 QTL in a RIL population developed from the cross of B73 × Ce3005, which could explain 1.15–12.43 % of the phenotypic variation. Li et al.  narrowed the QTL interval which had the largest effect among the 7 QTL of RPR detected in two RIL populations by the method of haplotype analysis. Peiffer et al.  reported that 18 family-nested QTL and 141 significant GWAS associations were identified for RPR across NAM (nested association mapping) and IBM (intermated B73 × Mo17) families, while numerous weak associations were found in the NCRPIS (North Central Regional Plant Introduction Station) diversity panel for RPR. Mutations, brittle stalk (BK) genes exhibiting a lower proportion of cellulose, had dramatically weakened tissue mechanical strength than that of wild type stalks .
Moreover, whole plant digestibility, which can reflect the feeding value, has been extensively studied in forage maize, and several reports of QTL analyses with low-density markers for stalk digestibility in forage maize were published [24, 25]. Maize mutants and/or genetically engineered plants have highlighted a few genes affecting maize cell wall degradability [26, 27]. Reports have emerged on nucleotide diversity and the extent of linkage disequlibrium (LD) at the gene locus of lignin and cellulose synthesis [28–30].
It was well known that plant breeders are desired to choose lines based on minimizing negative effects of genotype agronomic value, so it was crucial to perform research on the genetic architecture of stalk traits, especially for RPR and IVDMD. In this study, we first used a DH population combined with the corresponding haploid population to identify QTL and observe candidate gene expression about stalk traits. Our objectives were to: (1) explore the genetic architecture of stalk traits; (2) evaluate consistence and difference in phenotype, QTL characters, and gene expression between two different ploidy populations in stalk traits; and (3) preliminary propose and illustrate a ploidy effect mechanism for RPR and IVDMD under one genome dosage situation with the QTL mapping method.
Performance of parental lines, F1 generation and DH and haploid populations derived from each DH line
Performance of parents and derived DH and haploid populations across five environments was presented in Table 1. RPR, water content (WC), acid detergent fiber (ADF), neutral detergent fiber (NDF), and cellulose(Cel) of the male parent Chang7-2 (C7-2) showed significantly higher values than those of the female parent Zheng58 (Z58) in both DH and haploid populations. In contrast, for IVDMD and WSC (water soluble carbohydrate), Z58 had a higher value than the male parent C7-2 in both populations. There was no significant difference in lignin (Lig) content between two parents in the DH and haploid populations. RPR and IVDMD showed a normal distribution in both two ploidy populations (Fig. 1). For all traits investigated in this study, coefficients of variation (CV) in the DH and haploid population ranged from 7.56 to 49.48 % and from 8.28 to 35.28 %, respectively. The genotypic variance (σ G 2 ) was significant at P < 0.01 in both the DH and haploid populations (Table 2). The broad-sense heritability (h B 2 ) of all traits in the DH population were intermediate to high (0.49<h B 2 <0.81) as well as in the haploid population (0.58<h B 2 <0.89). Notably, h B 2 for all traits were higher in the haploid population than in the DH population except for WC, of which h B 2 was slightly lower in haploid population (0.58) than in DH population (0.60).
Inter-population and intra-population phenotypic correlation
The phenotypic correlation coefficients of all stalk traits between the DH and haploid populations ranged from 0.38 to 0.69 (Fig. 2). Coefficients of phenotypic correlation among different traits in DH population showed similar patterns to those in haploid population. In both populations, ADF, NDF and Cel showed high positive correlation among themselves, significantly positively correlated with RPR but negatively correlated with IVDMD, Lig and WSC. RPR negatively correlated with IVDMD but with different correlation coefficients in DH and haploid populations, respectively (Table 3).
Constructing a linkage map and the characteristics of markers
A total of 190 DH lines were used for genotyping with MaizeSNP3K chip, which was carried out on the Illumina Golden-Gate SNP genotyping platform  and 2956 high-quality SNPs were detected. The missing rate for these SNPs ranged from 0 to 20.00 % (average 1.50 %), the heterozygosity ranged from 0 to 14.21 % (average 2.06 %). A total of 4.74 % (9/190) of the DH lines with SNP heterozygosity ≥ 10 % were excluded in further analysis. Minor allele frequency (MAF) for these SNPs ranged from 0 to 0.50 (average 0.42) (Additional file 1: Table S2). Of these high-quality SNPs, 1318 SNPs were polymorphic between the two parental lines, and the marker distribution frequency for the two parents ranged from 30 to 65 % (Additional file 1: Figure S3). After quality control, 1137 SNPs were left and used to construct a linkage map using the Joinmap4.0 instructions . The total length of the linkage map was 1426.83 cM with an average interval of 1.26 cM (Additional file 1: Table S3).
QTL characteristics in the DH and haploid populations
Across five environments, the number and position of QTL detected in the DH and haploid populations was shown in Fig. 3. For each trait evaluated in this study, one or more QTL were identified in one region or even shared the same support intervals with the distance of less than 20 cM between the DH and haploid populations.
In the haploid population, four QTL for RPR were detected on chromosomes 1, 2, 3 and 5, two of which were identified on chromosomes 1 and 5 using the DH population (Table 4 and Fig. 3). The position of QTL identified in the haploid population on chromosome 1 was close to that detected in the DH population, which accounted for 6.60 and 8.00 % of the RPR genetic variation, respectively. The favorable alleles of RPR QTL were contributed by the C7-2 parental line in the DH population. All QTL of RPR detected in the DH and haploid populations together explained 25.90 and 42.90 % of the RPR genetic variation, respectively. The favorable alleles of RPR QTL on chromosomes 1, 2 and 3 were contributed by the RPR-higher parent C7-2, while the RPR-lower parent Z58 donated the alleles on chromosomes 1 and 5.
For IVDMD, three QTL were identified in each population, which explained 8.60–18.50 % of total genetic variation in the DH population and 6.80–18.60 % in haploid population These QTL were detected on chromosomes 1, 2, and 8 in the DH population and on chromosomes 5, 6 and 8 in the haploid population. Two QTL detected on chromosome 8 were tightly linked, which explained the 16.00 and 18.60 % of IVDMD genetic variation, respectively, and both were contributed by the IVDMD-higher parent Z58 in the DH and haploid populations.
In the DH population, QTL of RPR, IVDMD, ADF, NDF and WSC shared the same region ranging from 39.91 cM to 59.43 cM on chromosome 1. The QTL intervals for RPR, ADF, NDF and Cel detected in the haploid population ranging from 46.82 cM to 54.19 cM were completely included in the region described above. On chromosome 2, QTL of IVDMD detected in the DH population and the QTL of RPR detected in the haploid population were located adjacent to each other and shared common regions with the QTL of ADF, NDF and Cel. Two or more QTL located in bin 8.04/8.05 for IVDMD, ADF, NDF, Cel, and WSC clustered in the same chromosome region ranging from 78.11 cM to 94.79 cM in the DH and haploid populations.
Candidate gene identification for RPR and IVDMD in the DH and haploid populations
With a relatively high mapping resolution, some QTL representing the small genomic regions and the linear B73 genome can be used for searching candidate genes related to RPR and IVDMD. Based on the available annotation of the B73 reference sequence Version 5b.60 (http://ftp.maizesequence.org/release-5b/filtered-set/), we applied the MapMan BIN classification  and maizeGDB website (http://www.maizegdb.org/) to search for candidate genes.
In this study, four QTL which equally assigned for RPR and IVDMD in the two populations could account for more than 15.00 % of the genetic variation. The number of genes located in the four QTL of RPR and IVDMD were different from each other based on gene screening using qteller3 (http://qteller.com/qteller3/index.php) (Additional file 1: Table S5–S8). Nineteen genes have previously been demonstrated to be associated with cell wall formation mainly involved with cellulose and lignin synthesis (Table 5) [34–43]. Candidate genes participating same bioprocess were predicated between the DH and haploid populations for RPR and IVDMD, which were consistent with the results proposed by previous studies [20, 21, 44]. Moreover, although some evidence illustrated that some transcription factors, such as NAC, R2R3-MYB, C2H2, C3HC4 transcription factors families and so on, were associated with the cell wall [45–52], there was no clear evidence and further investigation was necessary to confirm the function of other annotated genes encoding similar transcription factors.
Transcriptional expression analyses of key genes involved in lignin and cellulose synthesis for haploid and diploid version of parental lines
To determine whether key genes involved in lignin biosynthesis were ploidy-modulated at a transcriptional level, relative expression levels of five genes, PAL, COMT, ccoAOMT, CCR and CAD, were analyzed (Fig. 4). In parental line Z58, transcript levels of the five genes increased 1.57–5.30folds in haploid plant relative to diploid plants. Particularly, COMT had the largest change fromZ58 in diploids to Z58 in haploids, while with no significant change from C7-2 in diploids to C7-2 in haploids. In the other parental line C7-2, higher expression levels of haploids than diploids were also observed across five genes, however, the changes of expression level (1.61–2.12-fold) from diploids to haploids were lower than what observed in Z58.
We also examined expression of cellulose synthesis genes encoding glycosyltransferase which detected in both populations for RPR and IVDMD (Fig. 4), and CesA11 and CesA12 were consistently co-expressed at all developmental stages of and were predominantly associated with the deposition of the secondary cell wall in maize stems even after the anthesis stage . For C7-2, the two genes, CesA11 and CesA12, were up-regulated 2.08–fold and 3.04–fold, respectively, in the haploid version relative to the diploid version. For Z58, the expression of CesA12 increased 1.48 fold in haploids in comparison to diploids, however, haploids had lower expression level for CesA11.
Performance and heritability of stalk traits in DH and haploid populations
The previous genetic investigations on RPR and IVDMD in diploid populations revealed that the traits were likely polygenic in maize and were affected by several mechanisms and complicated by confounding factors.
In this study, two parent lines presented consistent trends on RPR, as well as other traits, between the DH and haploid populations. This result showed each parent contributed coherent negative or positive allele effects even under different genome dosages and additive effect may played important role for partial phenotypic variation.
Heritability estimation depended on genetic background of the material, population types surveyed, interaction with environments and experimental design . In this study, the heritability of RPR was 0.72 and 0.78 estimated in DH and haploid population respectively, which were all in agreement with a high heritability of RPR reported in previous studies [19–21]. However, Peiffer et al.  reported that the heritability of RPR estimated from 26 RIL families of maize nested association mapping (NAM) population ranged from 0.08 to 0.34 (average 0.21), which may be due to a wide range in the flowering time among NAM families. Likewise, a high heritability was obtained for IVDMD in both DH (0.81) and haploid (0.89) populations, which were different from a moderate heritability, ranging from 0.55 to 0.68, reported in previous studies [55–58]. Moderate or high heritability values were obtained for other stalk traits as well. In conclusion, the high heritability values of stalk traits evaluated in this study could provide solid basis for QTL mapping analysis.
For most of traits evaluated in this study, the heritability estimated from the haploid population was higher than that from the DH population. Our results were consistent with heritability of DH and its haploid populations in maize reported by Geiger et al. , although different traits were studied. The higher heritability estimated from the haploid population can be explained by the relatively smaller σ G × E 2 and σ e 2 in haploid population than in DH population, which can be explained by that haploid lines mainly reacted more sensitively than DH lines to biotic and abiotic stress and therefore effectively uncover susceptibility to diseases and outer constraints, which had been proposed by Chase et al.  and Geiger et al. . In addition, all traits evaluated in this study were measured with high precision and then had a solid genetic basis in two ploidy populations, fundamentally suggesting that the haploid population as well as the DH population could be used in QTL analysis.
Phenotypic correlations and QTL co-localization for the same trait between the DH and haploid populations
Geiger et al.  reported moderate to high correlations between the DH and haploid lines from three material sets (KWS, SWS, and MON) for early vigor, silking, plant height, and stover weight per plant. We also observed a significant (P < 0.0001) moderate to strong positive correlation (r = 0.38-0.69) between the DH and haploid populations for all stalk traits (Fig. 2). This could suggest that moderate to strong correlations can occur independently of material background and trait restrictions. This high correlation between haploids and corresponding DH lines may provide reference information for maize breeders to select desirable lines at haploid stage, which could reduce breeding costs. However, the genetic mechanism on the connection between the DH and haploid populations has not yet been studied and therefore, is still unclear. In this study, through QTL mapping studies conducted in DH and its haploid population, we intended to understand this issue in term of genetic architecture. We first identified common QTL regions between the DH and haploid populations for each stalk trait, which could be considered as the genetic reason for the phenotypic correlation. Other QTL located on different chromosomes or having larger distance (>20 cM) may be partially caused by the change of genome dosage and explained by the different population size. Ming et al.  reported that many QTL for sugar content detected in sugarcane autopolyploids were not consistent with known candidate genes and suggested that other approaches will be necessary to isolate the genetic determinants of high sugar content of vegetative tissues. Until now, QTL detection in haploid population has not been reported.
Phenotypic correlations and QTL co-localizations among different traits
In an attempt to further understand the genetic architecture of RPR and IVDMD in maize, genomic regions for RPR, IVDMD and other stalk component traits were compared and phenotypic correlations between RPR, IVDMD and other stalk component were evaluated. Forty-seven QTL were identified in the DH and haploid populations (Additional file 1: Table S4, Fig. 3 and Table 4). The incidence of QTL clusters in similar genomic regions reflected trait associations .
Two studies have proposed that genes associated with the biosynthesis of cell wall components were considered as candidate genes for RPR [19, 20]. We also observed the positive correlations and QTL co-location of RPR with ADF, NDF and Cel, which were consistent with previous studies. RPR was negatively correlated with WC in a high-oil RIL population . The same correlation trend of RPR with WC and WSC were observed in the DH and haploid populations, except that RPR had no correlation with WC in the haploid population. In addition, Hu et al.  reported that the internode diameter, fresh weight of internode and dry weight of internode were also significantly positively correlated with RPR, and the difference in planting years, densities and maize varieties led to different stalk RPRs .
IVDMD showed the opposite correlation direction as the correlations of RPR with ADF, NDF, Cel and WSC, and had the same correlation direction as the correlations of RPR with WC. Therefore, WC may be one of the improved elements for practical breeding for stalk lodging resistance and forage maize. Several QTL associated with IVDMD and other stalk components were located in the same bins as identified in our studied [57, 58, 63]. Lig was positively correlated with IVDMD and was not correlated with RPR, which was not in agreement with previous studies . This may be due to the no-forage background materials used in this study. No reports were available on QTL both for RPR and IVDMD. Only in the DH population evaluated in this study, we first detected one QTL cluster for RPR and IVDMD at bin 1.10 and bin 1.07, respectively. However, we found more than one QTL of RPR or IVDMD sharing common regions or flanking markers with the QTL of other stalk components, which suggested close linkage or pleiotropy as the explanation for the correlations and some common genes had effects on RPR and IVDMD. The QTL clusters could be deployed for improving RPR and IVDMD in maize through marker-assisted selection.
Compare QTL identified in this study with those identified in previous studies in diploid populations
We have identified additive QTL for RPR on chromosomes 1, 2, 3 and 5. Flint-Garcia et al.  detected one QTL region on chromosome 3 contained overlapping support intervals across four F2:3 maize populations. Also, in less than four populations, other QTL were detected at bins 1.07–1.09, 2.02, 2.06–2.07, 3.04–3.08, and 5.02. Similarly, our mapping study of the DH and haploid populations identified five RPR QTL located near bins 1.07, 1.10, 2.02, 3.09 and 5.01. Hu et al.  investigated RPR in a RIL population derived from a high-oil population and reported that RPR QTL were detected on all chromosomes except for chromosome 5 and the QTL located in bin 3.06 was the most important one and it accounted for 12 % of the phenotypic variation. Li et al.  identified seven RPR-associated QTL in two RIL populations. Among these QTL, the largest-effect QTL accounted for 18.9 % of the phenotypic variation was located at bin 3.06, and other QTL for RPR were observed at bins 2.10, 3.08, 9.03–9.04, 4.06, 6.05, and 6.07, explaining 4.40–13.80 % of the phenotypic variation. In the present study, the QTL location on chromosome 3 were only detected in the haploid populations and accounted for 10.30 % of the genetic variation with highly detected frequency in 1000 runs cross-validation (Additional file 1: Figure S4). Moreover, it is worth noting that RPR-associated QTL, which were observed at bin 2.02 in the haploid populations and explained more than 15.00 % of the contribution to genetic variation, were located in the same region as QTL detected by Flint-Garcia et al. . The QTL detected at bin 5.05 in the DH population could account for the highest percentage of RPR genetic variation (up to 16.90 %) and were not located in the QTL cluster with other traits, and this QTL has not been proposed by previous studies. Since these two newly discovered QTL were also detected with high frequencies in the 1000 cross-validation, this confirmed our conclusion that QTL at bins 2.02 and 5.05 likely carried major candidate genes for RPR (Additional file 1: Figure S4).
Six QTL for IVDMD in total were detected in DH and its haploid population in this study. Two QTL detected in the DH and haploid populations were located in adjacent bins 8.04 and 8.05 with a genetic distance of less than 3 cM. These two QTL also showed high detection frequencies in cross-validation (Additional file 1: Figure S4). Similarly, Wei et al.  reported that a IVDMD QTL located at bin 8.06–8.07 were detected in Pop2 combined analysis, which was adjacent to QTL for IVDMD on chromosome 8 detected in this study. Other QTL for IVDMD identified in the DH and haploid populations were distributed on chromosomes 1, 2, 5 and 6. IVDMD QTL located at bin 1.07 can explain 18.50 % of the genetic variation. Previous reports showed that QTL on chromosome 1 had a great effect on stalk digestibility [57, 63]. The QTL located at bins 5.02–5.03 and 5.03–5.06 were detected in Xuchang and Luoyang Pop2, respectively, by Wei et al. . Wang et al.  suggested that IVDMD QTL explained more than 10 % of the genetic variation in both F3 and F4 generations were mapped on the same genomic position on chromosome 6, which were the same as QTL detected in maize recombinant inbred line progeny of F288 × F271 . One IVDMD QTL detected in our study was also on chromosome 6. These QTL described above were closely linked under high-density SNP markers and deserve further investigation for finding candidate genes underlying IVDMD in a no-forage genetic background.
The role of genome dosage changes on gene expression of lignin and cellulose synthesis in inbred and haploids of two parental lines
Most candidate genes were involved in lignin and cellulose synthesis which affect the stalk cell wall structure. Lignin was a phenolic polymer that imparted mechanical strength of the plant secondary cell wall, and therefore, was considered to confer stalk rot resistance and involve in plant evolution . Particularly, genes participating in lignin synthesis were identified only in haploid population in our study. Therefore, based on the gene function annotations for RPR and IVDMD QTL detected in the DH and haploid populations, we analyzed the key gene expressions of lignin and cellulose synthesis and genome dosage regulation. The expression levels and phenotypes showed several interesting results, suggesting a partial explanation for ploidy effect mechanisms in the haploid condition.
In the one dosage genome, compared to the inbred, CesA11 and CesA112 gene expressions were up-regulated except for the CesA11 and CesA12 gene in Z58, which was consistent with the decreased Cel content in haploid Z58 and higher content in C7-2 haploids (Fig. 4). Unlike other gene expressions in lignin synthesis, the COMT gene showed significantly lower expression levels in C7-2 haploids than that in Z58 haploids. All these results illustrated the existence of genetic variation in morphological responses to ploidy changes, which was consistent with the results which reported by Riddle et al. .
All lignin and synthesis gene expression showed dramatically inversed expression levels, except for CesA11 in Z58, as genome dosage number increased (Fig. 4), which was consistent with higher Lig contents for two parents and increased Cel contents for parent line C7-2 in the haploid condition and can be explained by dosage compensation . These findings implied that the major effect of genome dosage changes were not simply proportional to copy number of regulatory genes and may not be directly related to the phenotype, but rather result from various regulatory components which resulted in compensation for gene expression in the one dosage genome. Previous studies had also reported that many target loci exhibited dosage compensation, such as Adh expression with reduced alleles number . Quantitative traits trended to be regulated by genes exhibiting dosage effects and most were transcription factors [69, 70]. Forty-seven dosage-dependent modifiers operated in macromolecular complexes in a regulatory network . Many other mechanisms, such as siRNA, DNA methylation and so on, may be concerned with expression changes in dosage- sensitive region genes, which can also sharpen the evolution of copy-number varied regions [72–74].
Using DH and its corresponding haploid populations, this analyse revealed important genomic regions associated with eight stalk-related traits. These QTL explained a large extent of phenotypic variance. One or more QTL sharing common region for each stalk-related trait were detected between this two different ploidy populations. The heritabilities in haploid population were higher than DH population in all stalk traits except WC. This study identified candidate genes involving in lignin and cellulose synthesis for RPR and IVDMD, which were the two most important stalk traits. The expression levels of most of these candidate genes were significant higher in haploid parents than that in corresponding diploid parents. Haploid population may be used as one of platforms providing information on the genetic basis for stalk-related traits, and the genetic connection between DH and haploids for traits can facilitate the selection of materials at haploid-stage which can boost the practical maize breeding.
The dosage compensation mechanism and dosage-sensitivity genes may be further examined by analyzing the genetic architecture of certain traits by comparing QTL mapping results between different ploidy populations, which may also have important implications for understanding gene regulatory networks and genome evolution.
Materials and population construction
The parental lines of two populations belong to two distinct maize germplasm groups. The maternal inbred line Zheng58 (Z58) is dent corn of the Reid heterotic group, and the paternal inbred line Chang7-2 (C7-2) is flint corn. The single cross (ZD958) of these two parental lines is one of the most widely planted commercial varieties in China.
The DH population was constructed by production of haploids with an in vivo haploid induction procedure and followed by chromosome doubling (Additional file 1: Figure S1).
In total, the DH population consisted of 190 DH lines. Accordingly, the haploid population was constructed by crossing each DH line, used as source germplasm, with a haploid inducer, used as pollinator, and we finally produced 170 haploid lines. To obtain enough DH and haploid seeds for evaluating phenotypes in field trial, in the experimental station of the China Agricultural University (Yacheng in Hainan Province), we planted two replications of 190 DH lines with 21 plants in each plot that again pollinated with inducer line CAU5 at the anthesis stage . The haploid seeds were manually picked up from the harvested ears based on the R1-nj color markers .
In 2013 and 2014, the two populations, their parental lines (inbred and haploid) and F1 generation were sown in two experimental stations of the China Agricultural University (Shangzhuang in Beijing and Quzhou in Hebei Province), and Shijiazhuang experimental station of Hebei Academy of Agriculture and Forestry Sciences. At each location, the two populations were adjacently planted to reduce the influence of field heterogeneity. We used a randomized complete block design with two blocks for both populations. Within each block, each line was assigned to a single row plot (Additional file 1: Figure S2).
The fourth internode above ground (FIAG) was measured for nine traits described as follows. Three randomly selected plants in each row (i.e., plot) were chosen for trait evaluation. RPR was measured with an electronic penetrometer (AWOS-SL04, Aiwoshi Company, Hebei, China) at the milky stage. The measurement of RPR and water content (WC) of the internode, as well as measurement of IVDMD, acid detergent fiber (ADF), neutral detergent fiber (NDF), lignin (Lig), cellulose (Cel) and water soluble carbohydrate (WSC) followed a standard procedure described by Hu et al. .
Phenotypic data analysis
In order to evaluate phenotypic traits, all DH lines and haploid lines, were planted in 3 locations in 2013 and 2014, which were treated totally as 5 independent macro environments (not including Quzhou 2014 due to large proportion of missing data points,). Then a linear model was used to perform an analysis of variance (ANOVA), genotypic value estimation (BLUE) and variance components estimation for the DH and haploid populations for each trait:
where μ is the grand mean, G is genotypic effect, E is the environment effect, G × E is the genotype-by-environment interaction, R(E) is the effect of block within environment, and ε is the random error. All statistical analyses were performed using SAS version 9.3 . The genotypic effect G was estimated when it was treated as a fixed effect using LSMEANs with PROC GLM. The genotypic value estimates (BLUEs) were used to calculate correlation coefficients by PROC CORR with Pearson’s method. The variance components were estimated using PROC VARCOMP with the method of restricted maximum likelihood (REML). The estimates of genotypic variance (σ G 2 ), genotype-by-environment interaction (σ G × E 2 ) and random error (σ e 2 ) were used to estimate heritability based on the formula :
where l is the number of macro environments and r is the number of blocks within each environment, which equals 5 and 3, respectively in the study.
Genotyping and genetic map construction
Genomic DNA was extracted from young leaves of DH and the parental lines using a CTAB method , and was then purified. The purified DNA was genotyped with the maizeSNP3K chip (3,072SNPs), which is a subset of the Illumina MaizeSNP50 BeadChip . SNP genotyping was performed on the Illumina Golden-Gate SNP genotyping platform at the National Maize Improvement Center of China of the China Agricultural University. Checking the quantity of each SNP was carried out manually as described by Yang et al. . The qualified SNPs were reserved for further screening and creation of linkage maps.
After genotyped, 1228 polymorphic markers between the two parental lines were determined. Subsequently, heterozygosity of each line, the missing rate, minor allele frequency (MAF) and heterozygosity of each SNP were calculated. The DH lines with heterozygosity ≤ 0.1 and the SNPs having polymorphisms between two parents with a missing rate ≤ 0.2 and MAF ≥ 0.05 were selected to construct a genetic linkage map with software package Joinmap4.0. SNPs and lines that shared 100 % similarities were deleted according to the Joinmap4.0 instructions . The linkage groups were composed at a minimum Logarithm of odds (LOD) of 7, and the Kosambi mapping function and regression-mapping algorithm were used for calculating map distances.
Since the segregating populations used in this study were DH and haploid populations, an additive genetic model was chosen for QTL analysis, using the BLUEs across environments as phenotypic data. Composite interval mapping (CIM) with a regression approach  in combination with the use of cofactors was employed . A two-step procedure was utilized for QTL detection as described by Hu et al. . A threshold of LOD =3.0 was used to determine QTL for all traits based on 2000 permutations . 1-LOD support intervals were calculated from the significant peak to a certain position on both sides along the chromosome at which the LOD score had a 1.0 unit decrease . To facilitate comparisons among linked QTL, two QTL were designated as overlapping when they were separated by less than 20 cM . The total proportion of genotypic variation (p G ) explained by the detected QTL was calculated by the formula p G = R adj 2 /h B 2 . Fivefold cross-validation was used to assess the reliability of QTL mapping results with 1000 runs .
Physical positions of QTL were obtained from the physical positions of their flanking markers. All genes between the flanking markers of each QTL were extracted from the filtered gene set of the maize genome sequence (http://ftp.maizesequence.org/release-5b/filtered-set/) by qteller3 (http://qteller.com/qteller3/index.php), assuming a linear relationship between recombination and physical distances within this interval.
RNA extraction and RT-PCR
Stalk rind used in gene expression analysis was collected from plant FIAG at the milky stage. Total RNA was extracted using Ultrapure RNA Reagent (Cat#CW0581, CWbio.Co.Ltd, Beijing, China) according to the manufacturer’s instructions. Total RNA was digested with RNase-Free DNase I to remove any DNA (Cat# CW2090, CWbio.Co.Ltd, Beijing, China). cDNA was synthesized using HiFi-MMLVcDNA (Cat#CW0744, CWbio.Co.Ltd, Beijing, China). Two microliters of diluted cDNA was used for quantitative RT-PCR analysis using the primer pairs listed in Additional file 1: Table S1. An ABI7500 machine was used to conduct RT-PCR analysis . Actin was used as an internal control to calculate the relative expression levels in three biological replications using the analysis method of 2-△△CT .
Availability of supporting data and materials
The data supporting the results of this article can be found in the article and its additional files. The biological sequences used for gene RT-PCR in this study come from National Center for Biotechnology Information (NCBI) under the accession numbers L77912, M73235, AJ242981, NM_001112018, AJ005702, AY372245, AY372246 and DQ492681. All published datasets referred to in the manuscript are cited in the reference list.
Acid detergent fiber
best linear unbiased estimation
- CAD :
Cinnamyl alcohol dehydrogenase gene
- ccoAOMT :
Caffeoyl CoA O-methyltransferase gene
- CCR :
Cinnamoyl CoA reductase1 gene
- CesA11 :
cellulose synthase catalytic subunit 11
- CesA12 :
cellulose synthase catalytic subunit 12
composite interval mapping
- COMT :
coefficients of variation
fourth internode above ground
in vitro dry matter digestion
minor allele frequency
neutral detergent fiber
- PAL :
phenylalanine ammonia-lyase gene
- p G :
proportion of genotypic variation
quantitative trait loci
rind penetrometer resistance
single nucleotide polymorphism
water soluble carbohydrate
Flint-Garcia SA, Darrah LL, McMullen MD, Hibbard BE. Phenotypic versus marker-assisted selection for stalk strength and second-generation European corn borer resistance in maize. Theor Appl Genet. 2003;107:1331–6.
Sibale EM, Darrah LL, Zuber MS. Comparison of two rind penetrometers for measurement of stalk strength in maize. Maydica. 1992;37:111–4.
Abedon BG, Darrah LL, Tracy WF. Developmental changes associated with divergent selection for rind penetrometer resistance in the MoSCSSS maize synthetic. Crop Sci. 1999;39:108–14.
Argiller O, Barrière Y, Lila M, Jeanneteau F, Gélinet K, et al. Genotypic variation in phenolic components of cell-walls in relation to the digestibility of maize stalks. Agronomie. 1996;16:123–30.
Shi C, Uzarowska A, Ouzunova M, Landbeck M, Wenzel G, Lübberstedt T. Identification of candidate genes associated with cell wall digestibility and eQTL (expression quantitative trait loci) analysis in a Flint × Flint maize recombinant inbred line population. BMC Genomics. 2007;8:1471–2164.
Li ZQ, Zhang HM, Wu XP, Sun Y, Liu XH. Quantitative trait locus analysis for ear height in maize based on a recombinant inbred line population. Genet Mol Res. 2014;13:450–6.
Sibov ST, De Souza JRCL, Franco-Garcia AF, Mangolin CA, Benchimol LL, De Souza AP. Molecular mapping in tropical maize (Zea mays L.) using microsatellite markers. 2. Quantitative trait loci (QTL) for grain yield, plant height, ear height and grain moisture. Hereditas. 2003;139:107–15.
Li CH, Li YX, Shi YS, Song YC, Zhang DF, Buckler ES, et al. Genetic control of the leaf angle and leaf orientation value as revealed by ultra-high density maps in three connected maize populations. Plos one. 2015;10(3):e0121624. doi:10.1371/journal.pone.0121624.
Truntzler M, Barrière Y, Sawkins MC, Lespinasse D, Betran J, Charcosset A, et al. Meta-analysis of QTL involved in silage quality of maize and comparison with the position of candidate genes. Theor Appl Genet. 2010;121:1465–82.
Young ND, Debellé F, Oldroyd-Giles ED, Geurts R, Cannon SB, Udvardi MK, et al. The Medicago genome provides insight into the evolution of rhizobial symbioses. Nature. 2012;480:520–4.
Lagunes Espinoza LDC, Julier B. QTL detection for forage quality and stem histology in four connected mapping populations of the model legume Medicago truncatula. Theor Appl Genet. 2013;126:497–509.
Gallais A, Bordes J. The use of Doubled Haploids in recurrent selection and Hybrid development in Maize. Crop Sci. 2007;47:190–201.
Geiger HH, Gordillo GA. Doubled Haploids in Hybrid maize breeding. Maydica. 2009;54:485–99.
Andersen JR, Asp T, Lu YC, Maitz MK, Ouzonova M, Lübberstedt T. Development and mapping of gene-tagged SNP markers in laccases of maize (Zea mays L.). Plant Breeding. 2009;128:423–5.
Dunwell JM. Haploids in flowering plants: origins and exploitation. Plant Biotechnol J. 2010;8:377–424.
Geiger HH, Gordillo GA, Koch S. Genetic Correlations among Haploids, Doubled Haploids, and Testcrosses in Maize. Crop Sci. 2013;53:2313–20.
Ravi M, Marimuthu AMP, Tan EK, Maheshwari S, Henry IM, Marin-Rodriguez B, et al. A haploid genetics toolbox for Arabidopsis thaliana. Nat Commun. 2014;5:5334.
Elling G, Penninger JM. Genome wide functional genetics in haploid cells. FEBS Lett. 2014;588:2415–21.
Flint-Garcia SA, Jampatong C, Darrah LL, McMullen MD. Quantitative Trait Locus Analysis of Stalk Strength in Four Maize Populations. Crop Sci. 2003;43:13–22.
Hu HX, Meng YJ, Wang HW, Liu H, Chen SJ. Identifying quantitative trait loci and determining closely related stalk traits for rind penetrometer resistance in a high-oil maize population. Theor Appl Genet. 2012;124:1439–47.
Li K, Yan JB, Li JS, Yang XH. Genetic architecture of rind penetrometer resistance in two maize recombinant inbred line populations. BMC Plant Biol. 2014;14:152–63.
Peiffer JA, Flint-Garcia SA, De Leon N, McMullen MD, Kaeppler SM, Buckler ES. The Genetic Architecture of Maize Stalk Strength. Plos one. 2013;8:e67066.
Ching A, Dhugga KS, Appenzeller L, Meeley R, Bourett TM, Howard RJ, et al. Brittle stalk 2 encodes a putative glycosylphosphatidylinositolanchored protein that affects mechanical strength of maize tissues by altering the composition and structure of secondary cell walls. Planta. 2006;224:1174–84.
Lübberstedt T, Melchinger AE, Fahr S, Klein D, Dally A, et al. QTL Mapping in Testcrosses of Flint Lines of Maize: III. Comparison across Populations for Forage Traits. Crop Sci. 1998;38:1278–89. https://www.researchgate.net/publication/250117573_QTL_Mapping_in_Testcrosses_of_Flint_Lines_of_Maize_III_Comparison_across_Populations_for_Forage_Traits.
Courtial A, Thomas J, Reymond M, Méchin V, Grima-Pettenati J, Barrière Y. Targeted linkage map densification to improve cell wall related QTL detection and interpretation in maize. Theor Appl Genet. 2013;126:1151–65.
Pichon M, Deswartes C, Gerentes D, Guillaumie S, Lapierre C, Toppan A, et al. Variation in lignin and cell wall digestibility in caffeic acid O-methyltransferase down-regulated maize half-sib progenies in field experiments. Mol Breeding. 2006;18:253–61.
Tamasloukht B, Sarah-Jane Wong Quai Lam M, Martinez Y, Tozo K, Barbier O, Jourda C, et al. Characterization of a cinnamoyl-CoA reductase 1 (CCR1) mutant in maize: effects on lignification, fibre development, and global gene expression. J Exp Bot. 2011;62:3837–48.
Andersen JR, Zein I, Wenzel G, Krützfeldt B, Eder J, Ouzunova M, et al. High levels of linkage disequilibrium and associations with forage quality at a Phenylalanine Ammonia-Lyase locus in European maize (Zea mays L.) inbreds. Theor Appl Genet. 2007;114:307–19.
Brenner EA, Zein I, Chen YS, Andersen JR, Wenzel G, Ouzunova M, et al. Polymorphisms in O-methyltransferase genes are associated with stover cell wall digestibility in European maize (Zea mays L.). BMC Plant Biol. 2010;10:27–39.
Chen YS, Blanco M, Li Q, Frei UK, Lübberstedt T. Extensive genetic diversity and low linkage disequilibrium within the COMT locus in maize exotic populations. Plant Sci. 2014;221–222:69–80.
Fan JB, Gunderson KL, Bibikova M, Yeakley JM, Chen J, Wickham Garcia E, et al.  Illumina Universal Bead Arrays. Methods Enzymol. 2006;410:57–73.
Van Ooijen JW. JoinMap 4. In: Software for the Calculation of Genetic Linkage Maps in Experimental Populations. Wageningen: Kyazma BV; 2006.
Thimm O, Bläsing O, Gibon Y, Nagel A, Meyer S, Krüger P, et al. MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 2004;37:914–39.
Doblin MS, Kurek I, Jacob-Wilk D, Delmer DP. Cellulose Biosynthesis in Plants: from Genes to Rosettes. Plant Cell Physiol. 2002;43:1407–20.
Blanton RL, Fuller D, Iranfar N, Grimson MJ, Loomis WF. The cellulose synthase gene of Dictyostelium. Proc Natl Acad Sci. 2000;97:2391–6.
Grabber JH, Hatfield RD. Methyl esterification divergently affects the degradability of pectic uronosyls in nonlignified and lignified maize cell walls. J Agric Food Chem. 2005;53:1546–9.
Clouse SD, Langford M, McMorris TC. A brassinosteroid-insensitive mutant in Arabidopsis thaliana exhibits multiple defects in growth and development. Plant Physiol. 1996;111:671–8.
Xu SL, Rahman A, Baskin TI, Kieber JJ. Two leucine-rich repeat receptor kinases mediate signaling, linking cell wall biosynthesis and ACC synthase in Arabidopsis. Plant Cell. 2008;20:3065–79.
Meyer K, Shirley AM, Cusumano JC, Bell-Lelong DA, Chapple C. Lignin monomer composition is determined by the expression of a cytochrome P450-dependent monooxygenase in Arabidopsis. Proc Natl Acad Sci. 1998;95:6619–23.
Showalter AM. Structure and function of plant cell wall proteins. Plant Cell. 1993;5:9–23.
Draeger C, Fabrice TN, Gineau E, Mouille G, Kuhn BM, Moller I, et al. Arabidopsis leucine-rich repeat extension (LRX) proteins modify cell wall composition and influence plant growth. BMC Plant Biol. 2015;15:155.
Peleman J, Boerjan W, Engler G, Seurinck J, Botterman J, Alliotte T, et al. Strong cellular preference in the expression of a housekeeping gene of Arabidopsis thafiana encoding S-adenosylmethionine synthetase. Plant Cell. 1989;1:81–93.
Vignols F, Rlgau J, Torres MA, Capellades M, Puigdomènech P. The brown midrib3 (bm3) mutation in maize occurs in the gene encoding caffeic acid O-methyltransferase. Plant Cell. 1995;7:407–16.
Barrière Y, Méchin V, Lefevre B, Maltese S. QTLs for agronomic and cell wall traits in a maize RIL progeny derived from a cross between an old Minnesota13 line and a modern Iodent line. Theor Appl Genet. 2012;125:531–49.
Mitsuda N, Seki M, Shinozaki K, Ohme-Takagi M. The NAC transcription factors NST1 and NST2 of Arabidopsis regulate secondary wall thickenings and are required for anther dehiscence. Plant Cell. 2005;17:2993–3006.
Kubo M, Udagawa M, Nishikubo N, Horiguchi G, Yamaguchi M, Ito J, et al. Transcription switches for protoxylem and metaxylem vessel formation. Genes Dev. 2005;19:1855–60.
Jin H, Cominelli E, Bailey P, Parr A, Mehrtens F, Jones J, et al. Transcriptional repression by AtMYB4 controls production of UV-protecting sunscreens in Arabidopsis. EMBO Journal. 2000;19:6150–61.
Borevitz JO, Xia YJ, Blount J, Dixon RA, Lamb C. Activation Tagging Identifies a Conserved MYB Regulator of Phenylpropanoid Biosynthesis. Plant Cell. 2000;12:2383–94.
Fornalé S, Shi XH, Chai CL, Encina A, Irar S, Capellades M, et al. ZmMYB31 directly represses maize lignin genes and redirects the phenylpropanoid metabolic flux. Plant J. 2010;64:633–44.
Sonbol FM, Fornalé S, Capellades M, Encina A, Touriño S, Torres JL, et al. The maize ZmMYB42 represses the phenylpropanoid pathway and affects the cell wall structure, composition and degradability in Arabidopsis thaliana. Plant Mol Biol. 2009;70:283–96.
Arruda P, Gerhardt IR. Nucleic acid molecules encoding plant proteins in the C3HC4 family and methods for the alteration of plant cellulose and lignin content. US Patent APP. 12/520, 266 US 2010/0077509 A1.
Rengel D, San SH, Servant F, Ladouce N, Paux E, Wincker P, et al. A new genomic resource dedicated to wood formation in Eucalyptus. BMC Plant Biol. 2009;9:1471–2229.
O’Brien BA. Functional analysis of the cellulose synthase gene family in Maize: from bioinformatics to assessments at whole-plant, tissue, cell, and protoplast levels. 3514971. University of Florida; 2011.
Hallauer AR, Miranda JB. Quantitative genetics in maize breeding. Ames: Iowa State University Press; 1981.
Méchin V, Argillier O, Hébert Y, Guingo E, Moreau L, Charcosset A, et al. Genetic analysis and QTL mapping of cell wall digestibility and lignification in silage maize. Crop Sci. 2001;41:690–7.
Papst C, Bohn M, Utz HF, Melchinger AE, Klein D, Eder J. QTL mapping for European corn borer resistance (Ostrinia nubilalis Hb.), agronomic and forage quality traits of testcross progenies in early-maturing European maize (Zea mays L.) germplasm. Theor Appl Genet. 2004;108:1545–54.
Wei MG, Li XH, Li JZ, Fu JF, Wang YZ, Li YL. QTL detection for stover yield and quality traits using two connected populations in high-oil maize. Plant Physiol Biochem. 2009;47:886–94.
Wang HW, Han J, Sun WT, Chen SJ. Genetic analysis and QTL mapping of stalk digestibility and kernel composition in a high-oil maize mutant (Zea mays L.). Plant Breeding. 2010;129:318–26.
Chase SS. Efficient methods of developing and improving inbred lines. The monoploid method of developing inbred lines. In: Report of the 6th Hybrid Corn Industry Research Conference, Chicago, IL. 28–29 Nov. 1951. American Seed Trade Association, Hybrid Corn Division, Chicago, IL. p. 29–34.
Ming R, Liu SC, Moore PH, Irvine JE, Paterson AH. QTL analysis in a complex autopolyploid: genetic control of sugar content in sugarcane. Genome Res. 2001;11:2075–84.
Chen YS, Lübberstedt T. Molecular basis of trait correlations. Cell. 2010;15:454–61.
Gou L, Huang JJ, Sun R, Ding ZS, Dong ZQ, Zhao M. Variation characteristic of stalk penetration strength of maize with different density-tolerance varieties. Trans Chin Soc Agric Eng. 2010;26:156–62.
Barrière Y, Thomas J, Denoue D. QTL mapping for lignin content, lignin monomeric composition, p-hydroxycinnamate content, and cell wall digestibility in the maize recombinant inbred line progeny F838 × F286. Plant Sci. 2008;175:585–95.
Roussel V, Gibelin C, Fontaine AS, Barrière Y. Bin and gene cluster genetic analysis in recombinant inbred lines of early dent forage maize. II - QTL mapping for cell wall constituents and cell wall digestibility from per se value and top cross experiments. Maydica. 2002;47:9–20.
Ruben V, Brecht D, Kris M, John R, Wout B. Lignin biosynthesis and structure. Plant Physiol. 2010;153:895–905.
Riddle NC, Kato A, Birchler JA. Genetic variation for the response to ploidy change in Zea mays L. Theor Appl Genet. 2006;114:101–11.
Casci T. What dosage compensation? Nat Rev Genet. 2010;12:2–2.
Birchler JA. The genetic basis of dosage compensation of Alcohol Dehydrogenase-1 in Maize. Genetics. 1981;97:625–37.
Liu J, Eck JV, Cong B, Tanksley SD. A new class of regulatory genes underlying the cause of pear-shaped tomato fruit. Proc Natl Acad Sci. 2002;99:13302–6.
Cong B, Barrero LS, Tanksley SD. Regulatory change in YABBY-like transcription factor led to evolution of extreme fruit size during tomato domestication. Nat Genet. 2008;40:800–4.
Birchler JA, Bhadra U, Bhadra MP, Auger DL. Dosage-dependent gene regulation in multicellular eukaryotes: implications for dosage compensation, aneuploid syndromes, and quantitative traits. Dev Biol. 2001;234:275–88.
Barbash DA, Siino DF, Tarone AM, Roote J. A rapidly evolving MYB-related protein causes species isolation in Drosophila. Proc Natl Acad Sci. 2003;100:5302–7.
Groszmann M, Greaves IK, Albertyn ZI, Scofield GN, Peacock WJ, Dennis ES. Changes in 24-nt siRNA levels in Arabidopsis hybrids suggest an epigenetic contribution to hybrid vigor. Proc Natl Acad Sci. 2011;108:2617–22.
Gerstein AC, Cleathero LA, Mandegar MA, Otto SP. Haploids adapt faster than diploids across a range of environments. J Evol Biol. 2011;24:531–40.
Xu XW, Li L, Dong X, Jin WW, Melchinger AE, Chen SJ. Gametophytic and zygotic selection leads to segregation distortion through in vivo induction of a maternal haploid in maize. J Exp Bot. 2013;64:1083–96.
Li L, Xu XW, Jin WW, Chen SJ. Morphological and molecular evidences for DNA introgression in haploid induction via a high oil inducer CAUHOI in maize. Planta. 2009;230:367–76.
Institute Inc. SAS user’s guide, version 9.3. SAS institute Cary, NC; 2015.
Utz HF, Melchinger AE, Schön CC. Bias and sampling error of the estimated proportion of genotypic variance explained by quantitative trait loci determined from experimental data in maize using cross validation and validation with independent samples. Genetics. 2000;154:1839–49.
Murray MG, Thompson WF. Rapid isolation of high molecular weight plant DNA. Nucleic Acids Res. 1980;8:4321–5.
Ganal MW, Durstewitz G, Polley A, Be´rard A, Buckler ES, Charcosset A, et al. A large maize (Zea mays L.) SNP genotyping array: development and germplasm genotyping, and genetic mapping to compare with the B73 reference genome. Plos One. 2011;6:e28334.
Yang XH, Xu YB, Shah T, Li HH, Han ZH, Li JS, et al. Comparison of SSRs and SNPs in assessment of genetic relatedness in maize. Genetica. 2011;139:1045–54.
Haley CS, Knott SA. A simple regression method for mapping quantitative trait loci in line crosses using flanking markers. Heredity. 1992;69:315–24.
Zeng ZB. Precision mapping of quantitative trait loci. Genetics. 1994;136:1457–68.
Hu HX, Liu WX, Fu ZY, Homann L, Technow F, Wang HW, et al. QTL mapping of stalk bending strength in a recombinant inbred line maize population. Theor Appl Genet. 2013;126:2257–66.
Doerge RW, Churchill GA. Permutation tests for multiple loci affecting a quantitative character. Genetics. 1996;142:285–94.
Mickelson SM, Stuber CS, Senior L, Kaeppler SM. Quantitative trait loci controlling leaf and tassel traits in a B73 × Mo17 population of maize. Crop Sci. 2002;42:1902–9.
Visscher PM, Thompson R, Haley CS. Confidence intervals in QTL mapping by bootstrapping. Genetics. 1996;143:1013–20.
Würschum T, Liu W, Gowda M, Maurer HP, Fischer S, Schechert A, et al. Comparison of biometrical models for joint linkage association mapping. Heredity. 2012;108:332–40.
Chen ZJ, Wang JL, Tian L, Lee HS, Wang JJ, Chen M, et al. The development of an Arabidopsis model system for genome wide analysis of polyploidy effects. Biol J Linn Soc Lond. 2004;82:689–700.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method. Methods. 2001;25:402–8.
This work was supported by the Modern Maize Industry Technology System (CARS–02–09), National Science and Technology Project (2014ZX08003–002) and the National High-Tech Program of China (2011AA10A103, 2012AA10A305).
The authors declare that they have no competing interests.
CSJ designed the study and assisted in writing the manuscript. LWX conducted the QTL mapping analysis and assisted in writing the manuscript. MYJ carried out the experiments, analyzed data and wrote the manuscript. HHX conducted partial data analysis and proposed suggestions on the manuscript. MYJ, LJH, LJJ and LW collected phenotypic data. All authors read and approved the final manuscript.
Route used to develop DH population and the corresponding haploid populations. Figure S2. Presentation of DH and haploid populations planted in field. Figure S3. Histogram for genotypic percentage in DH population. Figure S4. QTL detected frequency for RPR and IVDMD in DH and haploid populations with five-split cross validation for 1000 times. For RPR, the QTL at bin 2.02 and 5.05 were detected in DH and haploid population with high frequencies in the 1000 cross-validation, respectively. For IVDMD, the QTL at bin 8.04 and 8.05 were detected in DH and haploid population with high frequencies in the 1000 cross-validation, respectively. Table S1. Primer sequences of genes involved in lignin and cellulose synthesis for RT-PCR. Table S2. Summary of SNP characteristics of DH population. Table S3. Summary of the linkage map characteristics DH population. Table S4. QTL detected for other stalk traits in DH and haploid populations. Table S5. Annotation of the 45 predicted genes located within putative RPR QTL interval in DH population. Table S6. Annotation of the 191 predicted genes located within putative RPR QTL interval in haploid population. Table S7. Annotation of the 24 predicted genes located within the putative IVDMD QTL interval in DH population. Table S8. Annotation of the 191 predicted genes located within putative IVDMD QTL interval in haploid population. (DOCX 870 kb)