Genetic architecture of rind penetrometer resistance in two maize recombinant inbred line populations
© Li et al.; licensee BioMed Central Ltd. 2014
Received: 24 February 2014
Accepted: 29 May 2014
Published: 3 June 2014
Maize (Zea Mays L.) is one of the most important cereal crops worldwide and provides food for billions of people. Stalk lodging can greatly undermine the standability of maize plants and therefore decrease crop yields. Rind penetrometer resistance is an effective and reliable method for evaluating maize stalk strength, which is highly correlated with stalk lodging resistance. In this study, two recombinant inbred line populations were constructed from crosses between the H127R and Chang7-2 lines, and between the B73 and By804 lines. We genotyped these two populations and their parents using 3,072 single nucleotide polymorphism markers and performed phenotypic assessment of rind penetrometer resistance in multiple environments to dissect the genetic architecture of rind penetrometer resistance in maize.
Based on two linkage maps of 1,397.1 and 1,600.4 cM with average interval of 1.7 and 2.1 cM between adjacent makers, respectively, seven quantitative trait loci (QTL) for rind penetrometer resistance were detected in the two recombinant inbred line populations. These QTL were distributed in seven genomic regions, and each accounted for 4.4–18.9% of the rind penetrometer resistance variation. The QTL with the largest effect on rind penetrometer resistance, qRPR3-1, was located on chromosome 3 with the flanking markers PZE-103123325 and SYN23245. This locus was further narrowed down to a 3.1-Mb interval by haplotype analysis using high-density markers in the target region. Within this interval, four genes associated with the biosynthesis of cell wall components were considered as potential candidate genes for the rind penetrometer resistance effect.
The inheritance of rind penetrometer resistance is rather complex. A few large-effect quantitative trait loci, together with a several minor-effect QTL, contributed to the phenotypic variation in rind penetrometer resistance in the two recombinant inbred line populations that were examined. A potential approach for improving stalk strength and crop yields in commercial maize lines may be to introgress favorable alleles of the locus that was found to have the largest effect on rind penetrometer resistance (qRPR3-1).
KeywordsMaize Rind penetrometer resistance QTL SNP
Plant lodging is a complicated phenomenon that is affected by several factors, including genetics, environment and field management. Lodging is a considerable challenge for main crops during the growth as it often causes severe reduction in yields. In maize (Zea Mays L.), stalk lodging, breakage that occurs at or below the ear, can lead to loss of ears at harvest [1, 2]. It is estimated that yield losses caused by stalk lodging range from 5 to 20% worldwide [1, 3]. Additionally, stalk lodging poses an obstacle to mechanized harvesting, and consequently increases labor costs. Thus, improving stalk-lodging resistance has become a key target for maize breeding programs.
Developing an effective and accurate way to evaluate stalk-lodging resistance is a critical issue in improving maize stalk strength. Numerous quantitative methods have been developed to predict stalk lodging resistance potential, which mainly include chemical methods based on analysis of stalk chemical composition and anatomical structures, and mechanical methods based on measurements of stalk breaking, bending, penetration and crushing [4–9]. Among the mechanical methods, rind thickness and crushing strength have been useful in increasing lodging resistance in maize as they have shown a strong relationship with stalk lodging [7, 9]. However, these two methods are not ideal because they require the destruction of maize stalk. More recently, an efficient and non-destructive measure, rind penetrometer resistance (RPR), was developed to assess stalk strength . Increased RPR shows a high correlation with stalk-lodging resistance [11–15] and this method has been widely applied in estimating stalk lodging resistance potential in maize [2, 16–18] and in breeding maize hybrids that are highly resistant to stalk lodging [1, 11, 12, 19, 20].
Despite these advances in measuring stalk lodging, little was known about the genetic basis of stalk lodging and RPR variation until quantitative trait loci (QTL) mapping was applied to RPR. The first of these studies identified 35 individual QTL and 11 pairs of epistatic interactions associated with RPR in four F2:3 populations derived from B73, Mo47 and four inbred maize lines selected for stalk strength diversity . The majority of these QTL explained <15% of the phenotypic variation in RPR. An additional nine individual QTL and one more pair of epistatic interactions were detected in a recombinant inbred line (RIL) population crossed with a high-oil inbred line and B73, which account for another 1.15–12.43% of the phenotypic variation . Recently, 18 QTL and 141 significant associations for RPR were identified using a nested association mapping panel containing 4,536 lines and 174 intermated B73 × Mo17 RILs. Only 10 QTL were shared between two populations or two studies, reflecting the complex nature of stalk lodging [2, 17, 18].
The availability of the complete maize genome sequence  and haplotype maps [22, 23] have facilitated the development of single nucleotide polymorphism (SNP) genotyping technologies for maize . SNPs have become widely used markers in investigations of genetic variation and in linkage and association analyses in maize [25–32]. Compared with simple sequence repeat markers, SNPs are more accurate, less time-consuming and lower cost to identify, and better suited to high throughput genotyping platforms [33–35]. Currently, two GoldenGate assays containing 1,536 SNPs each and one Infinium BeadChip containing 56,110 SNPs have been developed for maize [34, 36, 37]. These SNP assays have been successfully used to examine the population structure and estimate genetic diversity of maize populations [34, 36, 37], and to identify variants associated with maize agronomic and quality traits [25–32].
In this study, we developed two RIL populations, H127R × Chang7-2 (referred to as POP-HRC) and B73 × By804 (referred to as POP-BYB), from four inbred lines with varying stalk strengths, and genotyped them using a GoldenGate maize SNP assay containing 3,072 SNPs to increase QTL resolution. Our objectives were to (1) identify QTL associated with RPR of maize stalks; (2) dissect the main-effect QTL with detailed haplotype in the target region; and (3) mine candidate genes associated with maize stalk strength.
Phenotypic variation in RPR
Descriptive statistics and broad-sense heritability for RPR in two RIL populations
Means ± SD (N/mm2)
37.64 ± 6.18
21.67 ± 2.63
23.25 ± 2.21
21.09 ± 2.82
Population mean ± SD (N/mm2)
27.19 ± 2.12
22.95 ± 1.03
F value Environment
Environment × Genotype
Heritability (H2) a (%)
Summary of SNPs and genetic linkage maps
Summary of SNP characteristics in the two RIL populations
Mean ± SD
Mean ± SD
High-quality SNP number
SNP missing rate (%)
0.83 ± 1.44
0.90 ± 1.60
MAF of SNPs
0.21 ± 0.22
0.23 ± 0.22
SNP heterozygosity (%)
2.38 ± 2.86
0.92 ± 1.10
SNP missing rate in each line (%)
0.83 ± 0.53
0.90 ± 2.46
SNP heterozygosity in each line (%)
2.40 ± 1.48
0.92 ± 1.12
Summary of the linkage map characteristics of the two RIL populations
Number of markers
Average interval (cM)
Minimum interval (cM)
Maximum interval (cM)
Number of markers
Average interval (cM)
Minimum interval (cM)
Maximum interval (cM)
RPR-associated QTL in the two RIL populations
Genetic interval (cM)
In POP-HRC, four of the RPR-associated QTL, located on chromosomes 2, 3 and 9 (Table 4, Additional file 1), together explained 50.4% of the phenotypic variation. The QTL on chromosome 3 flanked by the PZE-103123325 and SYN23245 markers, qRPR3-1, had the largest effect and accounted for 18.9% of the phenotypic variation in RPR. The H127R allele at this locus was correlated with a 1.05 N/mm2 increase in RPR. The second largest-effect QTL for RPR, qRPR9, which explained 8.1% of the phenotypic variation, was located between PZE-109058177 and PZE-109076761 on chromosome 9. The remaining two QTL, qRPR2 on chromosome 2 and qPPR3-2 on chromosome 3, explained 4.4% and 6.7% of the phenotypic variation, respectively. The alleles that were associated with increased RPR at these two loci also came from H127R.
In POP-BYB, the remaining three RPR-associated QTL, qRPR4, qRPR6-1 and qRPR6-2, accounted for 31.7% of the phenotypic variation. The B73 alleles at qRPR4 and qRPR6-2 were correlated with similar 0.39 N/mm2 increases in RPR. The By804 allele at qRPR6-1, had an additive effect of 0.27 N/mm2 for increased RPR.
To further confirm the seven RPR-associated QTL identified using the best linear unbiased prediction (BLUP) values, we also mapped RPR-associated QTL in the RIL populations grown in different environments and replications grown in the same environments (Additional file 2). The association with RPR was stable for qRPR3-1 in all environments/replications, whereas the remaining six QTL differed significantly in at least two environments/replications or showed obvious LOD peaks in different environments/replications. In addition to the original seven QTL, 13 QTL were identified that associated with RPR in one or two environments/replications.
Beyond individual QTL, one pair of epistatic QTL between qRPR3-1 and qRPR3-2 was detected in POP-HRC. The type of epistasis between qRPR3-1 and qRPR3-2 was additive interacted by additive. This pair of epistatic QTL explains 2.5% of the phenotypic variation with a positive effect on RPR coming from the parental digenic combination. None epistatic QTL were identified in POP-BYB.
Fine mapping of qRPR3-1in POP-HRC
Candidate genes in the target QTL region
Genetic characterization of RPR in maize
Precise phenotypic measures are crucial for genotype-phenotype association analysis . Previous studies have shown that RPR is highly associated with stalk-lodging resistance in maize [8, 11, 12, 20]. For example, divergent selection for stalk crushing strength in synthetic maize populations has resulted in increased RPR [6, 40]. In addition, RPR of the internodes below the uppermost ear has been found to be highly correlated with the RPR of internodes between the last ear and ground on maize plants . Thus, measuring RPR of the internode below the uppermost ear is suggested to be one of the best ways to evaluate stalk-lodging resistance in maize in the current status.
Although RPR is a complex quantitative trait that can be affected by environment, most of the phenotypic variation appears to be due to genetic factors. The broad-sense heritability of RPR in maize, estimated in two previous studies as well as our study, reached over 90% in some segregating populations [2, 17]. The high broad-sense heritability reflects the accuracy and feasibility of the method used to quantify RPR in these studies. Whereas, broad-sense heritability values of RPR in maize estimated from nested association map families are far lower than the values we estimated, ranging from 8 to 34% (averaged 21%) across 26 RIL populations . The reduced heritability values may be attributable to the different populations surveyed, differences in the growing environments, or to the relatively low number of replications examined for each line in their study . Further characterization of RPR in more bi-parent segregating populations is needed to reconcile these differences in heritability values.
The complex nature of RPR in maize
The present study identified seven RPR-associated QTL were identified in two RIL populations. Among these QTL, only the largest-effect QTL, qRPR3-1, was also identified in two previous studies by Flint-Garcia et al.  and Hu et al. ; qRPR3-2 and qRPR6-2 were also detected in the Flint-Garcia et al. study and were found to explain 6.7% and 13.8% of the phenotypic variation, respectively . Our study revealed that a few large-effect QTL, together with some minor-effect QTL, provide most of the genetic basis of RPR, consistent with previous studies [2, 17, 18]. Together with this study, a total of 69 RPR-associated QTL have been identified in 33 segregating populations. The phenotypic variation explained by the largest-effect QTL in each population ranged from 5.6 to 20.2%. Among these QTL, only ~10 were common in at least two populations. The low repeatability across populations may be due to the complex nature of RPR in maize, and the fact that most individual loci have small effects, which results in relatively small differences in RPR between parent strains [2, 17, 18].
In addition to single-effect QTL for RPR, Flint-Garcia et al.  and Hu et al.  detected 11 pairs of epistatic QTL in three F2:3 populations and one pair of epistatic QTL in one RIL population. The majority of these pairs of epistatic QTL explained <10% of phenotypic variations in each population. These findings, together with one pair of epistatic QTL identified in our study, suggest that epistasis also contributes to the genetic basis of RPR, but the effect of epistasis is relatively weak and not consistent in different bi-parent populations.
Pleiotropic loci for stalk components
RPR is a physical measure of maize stalk strength, and significantly correlates with levels of some cell wall components [17, 42]. Specifically, selection for stalk strength in a maize synthetic line using RPR found that divergent selection for increased RPR leads to increases in crude fiber, cellulose and lignin content in the plants . Similarly, RPR significantly correlates with the levels of cell wall components, including acid detergent fiber, acid detergent lignin, crude fiber, cellulose, and neutral detergent fiber in internodes of a Ce03005/B73 RIL population . The high correlation indicates that RPR-associated QTL may overlap, at least to some degree with QTL related to cell wall components of maize stalks.
Because RPR results from the combination of different cell wall component levels, it is likely that many of the same polymorphisms underlie the QTL associated with these two traits. Previous studies have identified 55 genomic regions related to cell wall components [43–45], including three loci in chromosome bins 3.06, 3.08 and 6.07 that are located near RPR-related genomic regions detected in our study. The largest-effect QTL for RPR in chromosome bin 3.06, qRPR3-1, co-localized with QTL clusters for acid detergent lignin and hemicelluloses detected in a F288 × F271 RIL × F286 top-cross population . Another RPR-associated QTL in chromosome bin 3.08 in POP-HRC, qRPR3-2, is close to a QTL for Klason lignin in an RIL population derived from a cross between an old Minnesota13 line and a modern Iodent line . This QTL also co-localizes with the hemicellulose-associated QTL in the F288 × F271 RIL × F286 top-cross population . The POP-BYB QTL in chromosome bin 6.07, qRPR6-2, co-localized with a QTL that controls acid detergent fiber in De811 × B73 RIL populations . Taken together, this information will help us to mine the candidate genes underlying QTL for RPR based on pleiotropy.
Co-localization of RPR-related QTL and candidate genes
Generally, the final goal of primary QTL mapping is to clone genes of interest. Association analysis combined with traditional linkage analysis can speed up the process of cloning genes . Recently, Xu et al.  reported that haplotype analysis using high-density markers within the target QTL interval in segregation populations is an effective way to narrow down the primary QTL with large effects. Applying this strategy, we narrowed qRPR3-1 to a 3.1-Mb window predicted to contain 86 genes. Although this still leaves many possible genes within the loci that could be responsible for the RPR association of qRPR3-1, knowledge of connection between RPR and cell wall components can guide the choice of candidate genes for further study and potential cloning. Consequently, the four candidate genes known to be involved in cell wall component metabolism are considered the most likely to be the RPR-related qRPR3-1 gene. GRMZM2G132169 encodes a laccase that has been shown to play a role in lignin polymerization in Arabidopsis thaliana [47, 48]. GRMZM2G359234 encodes a UDP-glucuronic acid decarboxylase that produces UDP-xylose, a substrate for xylan biosynthesis , and decreased xylan content leads to reduced stem mechanical strength . GRMZM2G440016 encodes pectin methylesterase, which catalyses the de-esterification of pectin, and is significantly associated with stem strength in Arabidopsis thaliana . GRMZM2G126077 encodes the precursor of pectate lyase, and is involved in the pectin biodegradation pathway. Although some evidence from Arabidopsis supports the idea that GRMZM2G359234 and GRMZM2G440016 are the top candidates to be the RPR-associated genes, further investigation is necessary to confirm this connection, such as candidate-gene association mapping, fine mapping, and functional validation. Besides these four candidate genes, the other genes of unknown function or genes with function not linked to cell wall metabolism may also be the true variant for qRPR3-1.
Application of RPR-related QTL to the improvement of maize stalk strength
Stalk strength is an important factor in breeding maize varieties to maintain grain yield. Phenotypic selection based on RPR has been successful in improving stalk strength in several maize synthetic populations [11, 12, 19, 20, 42, 52, 53]. Marker-assisted selection (MAS) is an alternative way to improve target traits [54–57], including disease resistance in maize . Flint-Garcia et al. compared the efficiency of MAS and phenotypic selection for RPR, and found that MAS was more efficient . Therefore, understanding the genetic architecture of stalk strength will enhance efforts to optimize stalk strength, and ultimately mitigate stalk lodging. The large-effect QTL for RPR, qRPR3-1, is a potential QTL for improving stalk strength via MAS. Additionally, qRPR3-1 was narrowed to a relatively small QTL interval, which reduces the probability of linkage drag with deleterious alleles. The genetic effects of qRPR3-1 associated with maize stalk cell wall components will provide additional phenotypic markers to guide the introgression of favorable alleles at the qRPR3-1 locus.
In this study, we mapped RPR-associated QTL in two RIL populations using medium SNP-density based linkage maps. We found that four QTL in POP-HRC and three in POP-BYB explained 50.4% and 31.7% of RPR variation, respectively. Only one of the seven QTL accounted for >15% of the RPR variations. These findings indicate that a few large-effect QTL and additional minor-effect QTL contribute to the phenotypic variation in RPR in the two RIL populations, reflecting the complex nature of stalk strength. The largest-effect QTL in chromosome bin 3.06 in POP-HRC, qRPR3-1, was narrowed to a 3.1-Mb interval by haplotype analysis using high-density markers in the target QTL interval. Within this interval, four genes associated with the biosynthesis of cell wall component were considered the most likely candidate genes for the qRPR3-1 locus. This information will be valuable for introgressing favourable alleles of qRPR3-1 into elite inbred lines to enhance stalk strength, and in turn mitigate stalk lodging.
One maize F9 RIL population, consisting of 200 lines, was derived from a cross between the B73 and By804 lines. B73 is an elite inbred line derived from the Iowa Stiff Stalk Synthetic maize population. By804 is an inbred line developed from a Beijing high-oil population. Due to the high heterozygosity of 12 RILs in this population (>10%), only the remaining 188 lines were selected for subsequent analysis. Another F6 RIL population, containing 215 lines, was constructed by crossing the inbred lines H127R and Chang7-2. H127R is a parental line of the elite hybrid Zhongnongda 4, and Chang7-2 is the male parent line of the hybrid Zhengdan985. H127R is more resistant to stalk lodging than Chang7-2. For simplicity, we refer to the B73 × By804 RIL population as POP-BYB, and the H127R × Chang7-2 RIL population as POP-HRC.
Field experiments and phenotyping
All 415 RILs, together with the four parent lines, were planted in a randomized complete block design from 2011 to 2013. For POP-HRC, two replications were planted in each of three environments, including Beijing in 2012 and in 2013 and Henan in 2013. For POP-BYB, one replication was planted in each of six environments, including Hainan in 2011 and in 2012, and Beijing, Henan, Chongqing and Yunnan in 2012. Each line was grown in a single 2.5 m row, rows were 0.67 m apart, and planting density was 45,000 plants/ha. The RPR of six randomly selected plants in each row was evaluated in the middle of the flat side of the internodes below the primary ear with an electronic penetrometer (AWOS-SL04, Aiwoshi Company, Hebei, China) at two weeks after flowering at the average level of each population, which roughly corresponded with the milk stage.
Phenotypic data analysis
The variance components of RPR were estimated using PROC GLM in SAS 9.2 (SAS Institute). The model for variance analysis for POP-HRC was: y ijk = μ + e l + r k(l) + f i + (fe) il + ϵ lik , where μ is the grand mean of RPR, f i is the genetic effect of the “i”th line, e l is environmental effect of the “l”th environment, (fe) il is the interaction effect between genetic and environmental effects, r k(l) is effect of replications within environments, and ϵ lik is the residual error. For POP-BYB, the interaction effect between environment and genotype was treated as residual error due to the fact that there were no replications within each environment. These variance components were used to calculate broad-sense heritability based on the population means . The broad-sense heritability in POP-HRC was estimated as , where is the genetic variance, is the interaction of genotype with environment, is the residual error, e and r represent the number of environments and replications in each environment. In POP-BYB, the broad-sense heritability was estimated as , where is the genetic variance, is the residual error, e stands for the number of environments. Confidence interval of h2 were calculated according the method described by Knapp et al. .
A mixed linear model was fitted to each RIL to obtain the BLUP for RPR: y i = μ + f i + e i + ϵ i , where y i is the phenotypic value of individual i, μ is the grand mean for all environments, f i is the genetic effect, e i is effect of different environments, and ϵ i is the random error. The grand mean was fitted as a fixed effect, and genotype and environment were considered random effects. The MIXED procedure in SAS9.2 (SAS Institute) was used to obtain the BLUP value.
Genotyping and genetic map construction
Genomic DNA was extracted from leaf tissue of the RILs and parent lines using the modified CTAB method  and used for genotyping with the MaizeSNP3K subset (3,072 SNPs) of the Illumina MaizeSNP50 BeadChip . SNP genotyping was performed on the Illumina GoldenGate SNP genotyping platform  at the National Maize Improvement Center of China, China Agricultural University. The quality of each SNP was checked manually as described by Yan et al. , and SNPs with poor quality were excluded for further analysis.
In each RIL population, the missing rate, MAF and heterozygosity for each SNP and the missing rate and heterozygosity for each line were calculated using PLINK packages . The SNPs with missing rates ≤20% and MAFs ≥0.05 were used to construct the genetic linkage map with JoinMap 4.0 , using the Kosambi mapping function for calculating map distances. Linkage groups were formed at a minimum LOD of 6, and a regression-mapping algorithm was used to calculate map distances.
Windows QTL Cartographer 2.5  was used for QTL detection with the RPR BLUP values across the different populations, environments and replications. The whole genome scan was performed using composite interval mapping with a 0.5 cM scanning interval between markers, and the window size was set at 10 cM. Model 6 of the Zmapqtl module was selected for detecting QTL and estimating their effects. Forward–backward stepwise regression with five controlling markers was used to control for background from flanking makers. After 1,000 permutations, the threshold LOD value was determined at a significance level of P < 0.05. The confidence interval of QTL position was determined with one-LOD support interval method . To estimate the interactions of significant QTL and their total phenotypic variation, multiple interval mapping (MIM) in Windows QTL Cartographer 2.5 was performed with Bayesian Information Criteria (BIC-M0) as criteria of MIM model .
Annotation of candidate genes
Based on the information available in the MaizeSequence database , the function of each gene within the largest-effect QTL interval was inferred from orthologues in Arabidopsis or rice. Additional protein prediction information was obtained from the InterPro module in the European Bioinformatics Institute database (http://www.ebi.ac.uk/interpro/) .
Rind penetrometer resistance
Quantitative trait loci
Recombinant inbred line
Single nucleotide polymorphism
Minor allele frequency
Logarithm of odds
Best linear unbiased prediction.
This study was supported by the National High-Tech Research and Development Program of China (2012AA101104) and the National Natural Science Foundation of China (31171561).
- 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-1336. 10.1007/s00122-003-1387-9.View ArticlePubMedGoogle Scholar
- 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. 10.2135/cropsci2003.0013.View ArticleGoogle Scholar
- Tuberosa R, Salvi S: QTL for agronomic traits in maize production. Handbook of Maize: Its Biology. Edited by: Bennetzen JL, Hake SC. New York: Springer, 2009:501-541.View ArticleGoogle Scholar
- Berzonsky W, Hawk J, Pizzolato T: Anatomical characteristics of three inbred lines and two maize synthetics recurrently selected for high and low stalk crushing strength. Crop Sci. 1986, 26: 482-488. 10.2135/cropsci1986.0011183X002600030009x.View ArticleGoogle Scholar
- Martin M, Russell W: Correlated responses of yield and other agronomic traits to recurrent selection for stalk quality in a maize synthetic. Crop Sci. 1984, 24: 746-750. 10.2135/cropsci1984.0011183X002400040028x.View ArticleGoogle Scholar
- Sibale E, Darrah L, Zuber M, Sleper D: Correlated responses to selection for stalk crushing strength in two maize populations. Maydica. 1992, 37: 115-122.Google Scholar
- Thompson DL: Stalk strength of corn as measured by crushing strength and rind thickness. Crop Sci. 1963, 3: 323-329. 10.2135/cropsci1963.0011183X000300040013x.View ArticleGoogle Scholar
- Twumasi-Afriyie S, Hunter R: Evaluation of quantitative methods for determining stalk quality in short-season corn genotypes. Can J Plant Sci. 1982, 62: 55-60. 10.4141/cjps82-008.View ArticleGoogle Scholar
- Zuber M, Grogan C: A new technique for measuring stalk strength in corn. Crop Sci. 1961, 1: 378-380.View ArticleGoogle Scholar
- Thompson D: Recurrent selection for lodging susceptibility and resistance in corn. Crop Sci. 1972, 12: 631-634. 10.2135/cropsci1972.0011183X001200050023x.View ArticleGoogle Scholar
- Chesang-Chumo JJ: Direct and correlated responses to divergent selection for rind penetrometer resistance in MoSCSS maize synthetic. PhD thesis. University of Missouri-Columbia 1993.Google Scholar
- Dudley J: Selection for rind puncture resistance in two maize populations. Crop Sci. 1994, 34: 1458-1460. 10.2135/cropsci1994.0011183X003400060007x.View ArticleGoogle Scholar
- Jampatong S, Darrah L, Krause G, Barry B: Effect of one-and two-eared selection on stalk strength and other characters in maize. Crop Sci. 2000, 40: 605-611. 10.2135/cropsci2000.403605x.View ArticleGoogle Scholar
- Kang MS, Din AK, Zhang Y, Magari R: Combining ability for rind puncture resistance in maize. Crop Sci. 1999, 39: 368-371. 10.2135/cropsci1999.0011183X0039000200011x.View ArticleGoogle Scholar
- McDevitt CM: Responses to three cycles of selection with a clamp for plant standability in two maize populations. PhD thesis. University of Missouri-Columbia 1999.Google Scholar
- Sibale E, Darrah L, Zuber M: Comparison of two rind penetrometers for measurement of stalk strength in maize. Maydica. 1992, 37: 111-114.Google Scholar
- Hu H, Meng Y, Wang H, Liu H, Chen S: 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-1447. 10.1007/s00122-012-1799-5.View ArticlePubMedGoogle Scholar
- 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-10.1371/journal.pone.0067066.PubMed CentralView ArticlePubMedGoogle Scholar
- 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-114. 10.2135/cropsci1999.0011183X003900010017x.View ArticleGoogle Scholar
- Albrecht B, Dudley J: Divergent selection for stalk quality and grain yield in an adapted× exotic maize population cross. Crop Sci. 1987, 27: 487-494. 10.2135/cropsci1987.0011183X002700030012x.View ArticleGoogle Scholar
- Schnable PS, Ware D, Fulton RS, Stein JC, Wei F, Pasternak S, Liang C, Zhang J, Fulton L, Graves TA: The B73 maize genome: complexity, diversity, and dynamics. Science. 2009, 326: 1112-1115. 10.1126/science.1178534.View ArticlePubMedGoogle Scholar
- Chia J-M, Song C, Bradbury PJ, Costich D, de Leon N, Doebley J, Elshire RJ, Gaut B, Geller L, Glaubitz JC: Maize HapMap2 identifies extant variation from a genome in flux. Nat Genet. 2012, 44: 803-807. 10.1038/ng.2313.View ArticlePubMedGoogle Scholar
- Gore MA, Chia J-M, Elshire RJ, Sun Q, Ersoz ES, Hurwitz BL, Peiffer JA, McMullen MD, Grills GS, Ross-Ibarra J: A first-generation haplotype map of maize. Science. 2009, 326: 1115-1117. 10.1126/science.1177837.View ArticlePubMedGoogle Scholar
- Rafalski A: Applications of single nucleotide polymorphisms in crop genetics. Curr Opin Plant Biol. 2002, 5: 94-100. 10.1016/S1369-5266(02)00240-6.View ArticlePubMedGoogle Scholar
- Buckler ES, Holland JB, Bradbury PJ, Acharya CB, Brown PJ, Browne C, Ersoz E, Flint-Garcia S, Garcia A, Glaubitz JC: The genetic architecture of maize flowering time. Science. 2009, 325: 714-718. 10.1126/science.1174276.View ArticlePubMedGoogle Scholar
- Poland JA, Bradbury PJ, Buckler ES, Nelson RJ: Genome-wide nested association mapping of quantitative resistance to northern leaf blight in maize. Proc Natl Acad Sci U S A. 2011, 108: 6893-6898. 10.1073/pnas.1010894108.PubMed CentralView ArticlePubMedGoogle Scholar
- Tian F, Bradbury PJ, Brown PJ, Hung H, Sun Q, Flint-Garcia S, Rocheford TR, McMullen MD, Holland JB, Buckler ES: Genome-wide association study of leaf architecture in the maize nested association mapping population. Nat Genet. 2011, 43: 159-162. 10.1038/ng.746.View ArticlePubMedGoogle Scholar
- Cook JP, McMullen MD, Holland JB, Tian F, Bradbury P, Ross-Ibarra J, Buckler ES, Flint-Garcia SA: Genetic architecture of maize kernel composition in the nested association mapping and inbred association panels. Plant Physiol. 2012, 158: 824-834. 10.1104/pp.111.185033.PubMed CentralView ArticlePubMedGoogle Scholar
- Li H, Peng Z, Yang X, Wang W, Fu J, Wang J, Han Y, Chai Y, Guo T, Yang N: Genome-wide association study dissects the genetic architecture of oil biosynthesis in maize kernels. Nat Genet. 2013, 45: 43-50.View ArticlePubMedGoogle Scholar
- Li Q, Yang X, Xu S, Cai Y, Zhang D, Han Y, Li L, Zhang Z, Gao S, Li J: Genome-wide association studies identified three independent polymorphisms associated with α-tocopherol content in maize kernels. PLoS One. 2012, 7: e36807-10.1371/journal.pone.0036807.PubMed CentralView ArticlePubMedGoogle Scholar
- Pan Q, Ali F, Yang X, Li J, Yan J: Exploring the genetic characteristics of two recombinant inbred line populations via high-density SNP markers in maize. PLoS One. 2012, 7: e52777-10.1371/journal.pone.0052777.PubMed CentralView ArticlePubMedGoogle Scholar
- Shutu X, Dalong Z, Ye C, Yi Z, Shah T, Ali F, Qing L, Zhigang L, Weidong W, Jiansheng L: Dissecting tocopherols content in maize (Zea mays L.), using two segregating populations and high-density single nucleotide polymorphism markers. BMC Plant Biol. 2012, 12: 201-10.1186/1471-2229-12-201.PubMed CentralView ArticlePubMedGoogle Scholar
- Jones E, Sullivan H, Bhattramakki D, Smith J: A comparison of simple sequence repeat and single nucleotide polymorphism marker technologies for the genotypic analysis of maize (Zea mays L.). Theor Appl Genet. 2007, 115: 361-371. 10.1007/s00122-007-0570-9.View ArticlePubMedGoogle Scholar
- Yan J, Yang X, Shah T, Sánchez-Villeda H, Li J, Warburton M, Zhou Y, Crouch JH, Xu Y: High-throughput SNP genotyping with the GoldenGate assay in maize. Mol Breed. 2010, 25: 441-451. 10.1007/s11032-009-9343-2.View ArticleGoogle Scholar
- Yang X, Xu Y, Shah T, Li H, Han Z, Li J, Yan J: Comparison of SSRs and SNPs in assessment of genetic relatedness in maize. Genetica. 2011, 139: 1045-1054. 10.1007/s10709-011-9606-9.View ArticlePubMedGoogle Scholar
- Yan J, Shah T, Warburton ML, Buckler ES, McMullen MD, Crouch J: Genetic characterization and linkage disequilibrium estimation of a global maize collection using SNP markers. PLoS One. 2009, 4: e8451-10.1371/journal.pone.0008451.PubMed CentralView ArticlePubMedGoogle Scholar
- Ganal MW, Durstewitz G, Polley A, Berard A, Buckler ES, Charcosset A, Clarke JD, Graner E-M, Hansen M, Joets J: 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-10.1371/journal.pone.0028334.PubMed CentralView ArticlePubMedGoogle Scholar
- MaizeSequence. http://ensembl.gramene.org/Zea_mays/Info/Index,
- Cobb JN, DeClerck G, Greenberg A, Clark R, McCouch S: Next-generation phenotyping: requirements and strategies for enhancing our understanding of genotype–phenotype relationships and its relevance to crop improvement. Theor Appl Genet. 2013, 126: 867-887. 10.1007/s00122-013-2066-0.PubMed CentralView ArticlePubMedGoogle Scholar
- Colbert T, Darrah L, Zuber M: Effect of recurrent selection for stalk crushing strength on agronomic characteristics and soluble stalk solids in maize. Crop Sci. 1984, 24: 473-478. 10.2135/cropsci1984.0011183X002400030010x.View ArticleGoogle Scholar
- Ling G, Jianjun H, Rui S, Zaisong D, Zhiqiang D, Ming Z: Variation characteristic of stalk penetration strength of maize with different density-tolerance varieties. Trans Chin Soc Agr Eng. 2010, 11: 156-162.Google Scholar
- Martin SA, Darrah LL, Hibbard BE: Divergent selection for rind penetrometer resistance and its effects on European corn borer damage and stalk traits in corn. Crop Sci. 2004, 44: 711-717. 10.2135/cropsci2004.0711.View ArticleGoogle Scholar
- 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-549. 10.1007/s00122-012-1851-5.View ArticlePubMedGoogle Scholar
- Krakowsky M, Lee M, Coors J: Quantitative trait loci for cell-wall components in recombinant inbred lines of maize (Zea mays L.) I: stalk tissue. Theor Appl Genet. 2005, 111: 337-346. 10.1007/s00122-005-2026-4.View ArticlePubMedGoogle Scholar
- Roussel V, Gibelin C, Fontaine A, Barriere Y: 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.Google Scholar
- Li L, Li H, Li Q, Yang X, Zheng D, Warburton M, Chai Y, Zhang P, Guo Y, Yan J: An 11-bp Insertion in Zea mays fatb reduces the palmitic acid content of fatty acids in maize grain. PLoS One. 2011, 6: e24699-10.1371/journal.pone.0024699.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhao Q, Nakashima J, Chen F, Yin Y, Fu C, Yun J, Shao H, Wang X, Wang Z-Y, Dixon RA: LACCASE is necessary and nonredundant with PEROXIDASE for lignin polymerization during vascular development in Arabidopsis. Plant Cell. 2013, 25: 3976-3987. 10.1105/tpc.113.117770.PubMed CentralView ArticlePubMedGoogle Scholar
- Berthet S, Demont-Caulet N, Pollet B, Bidzinski P, Cézard L, Le Bris P, Borrega N, Hervé J, Blondet E, Balzergue S: Disruption of LACCASE4 and 17 results in tissue-specific alterations to lignification of Arabidopsis thaliana stems. Plant Cell. 2011, 23: 1124-1137. 10.1105/tpc.110.082792.PubMed CentralView ArticlePubMedGoogle Scholar
- Suzuki K, Suzuki Y, Kitamura S: Cloning and expression of a UDP‒glucuronic acid decarboxylase gene in rice. J Exp Bot. 2003, 54: 1997-1999. 10.1093/jxb/erg203.View ArticlePubMedGoogle Scholar
- Lee C, Teng Q, Zhong R, Ye Z-H: Arabidopsis GUX proteins are glucuronyltransferases responsible for the addition of glucuronic acid side chains onto xylan. Plant Cell Physiol. 2012, 53: 1204-1216. 10.1093/pcp/pcs064.View ArticlePubMedGoogle Scholar
- Hongo S, Sato K, Yokoyama R, Nishitani K: Demethylesterification of the primary wall by PECTIN METHYLESTERASE35 provides mechanical support to the Arabidopsis stem. Plant Cell. 2012, 24: 2624-2634. 10.1105/tpc.112.099325.PubMed CentralView ArticlePubMedGoogle Scholar
- Alsirt A: Evaluation of frequency distribution changes resulting from bi-directional selection for rind penetrometer resistance in maize. PhD thesis. 1993, University of Missouri-ColumbiaGoogle Scholar
- Masole H: Evaluation of high and low divergent rind penetrometer resistance selection at three plant densities in maize. PhD thesis. 1993, University of Missouri-ColumbiaGoogle Scholar
- Lande R, Thompson R: Efficiency of marker-assisted selection in the improvement of quantitative traits. Genetics. 1990, 124: 743-756.PubMed CentralPubMedGoogle Scholar
- Mohan M, Nair S, Bhagwat A, Krishna T, Yano M, Bhatia C, Sasaki T: Genome mapping, molecular markers and marker-assisted selection in crop plants. Mol Breed. 1997, 3: 87-103. 10.1023/A:1009651919792.View ArticleGoogle Scholar
- Collard BC, Mackill DJ: Marker-assisted selection: an approach for precision plant breeding in the twenty-first century. Philos Trans R Soc Lond Ser B Biol Sci. 2008, 363: 557-572. 10.1098/rstb.2007.2170.View ArticleGoogle Scholar
- Huang N, Angeles E, Domingo J, Magpantay G, Singh S, Zhang G, Kumaravadivel N, Bennett J, Khush G: Pyramiding of bacterial blight resistance genes in rice: marker-assisted selection using RFLP and PCR. Theor Appl Genet. 1997, 95: 313-320. 10.1007/s001220050565.View ArticleGoogle Scholar
- Zhao X, Tan G, Xing Y, Wei L, Chao Q, Zuo W, Lübberstedt T, Xu M: Marker-assisted introgression of qHSR1 to improve maize resistance to head smut. Mol Breed. 2012, 30: 1077-1088. 10.1007/s11032-011-9694-3.View ArticleGoogle Scholar
- Holland JB, Nyquist WE, Cervantes-Martínez CT: Estimating and interpreting heritability for plant breeding: An update. Plant Breed Rev. 2003, 22: 9-112.Google Scholar
- Knapp S, Stroup W, Ross W: Exact confidence intervals for heritability on a progeny mean basis. Crop Sci. 1985, 25: 192-194. 10.2135/cropsci1985.0011183X002500010046x.View ArticleGoogle Scholar
- Murray M, Thompson WF: Rapid isolation of high molecular weight plant DNA. Nucleic Acids Res. 1980, 8: 4321-4326. 10.1093/nar/8.19.4321.PubMed CentralView ArticlePubMedGoogle Scholar
- Fan JB, Gunderson KL, Bibikova M, Yeakley JM, Chen J, Wickham Garcia E, Lebruska LL, Laurent M, Shen R, Barker D:  Illumina Universal Bead Arrays. Methods Enzymol. 2006, 410: 57-73.View ArticlePubMedGoogle Scholar
- Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, De Bakker PI, Daly MJ: PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007, 81: 559-575. 10.1086/519795.PubMed CentralView ArticlePubMedGoogle Scholar
- Van Ooijen J: JoinMap 4. Software for the Calculation of Genetic Linkage Maps in Experimental Populations Kyazma BV, Wageningen, Netherlands. 2006Google Scholar
- Basten CJ, Weir BS, Zeng Z-B: QTL cartographer: a reference manual and tutorial for QTL mapping. Department of Statistics. 1997, Raleigh, NC: North Carolina State UniversityGoogle Scholar
- Lander ES, Botstein D: Mapping mendelian factors underlying quantitative traits using RFLP linkage maps. Genetics. 1989, 121: 185-199.PubMed CentralPubMedGoogle Scholar
- Kao C-H, Zeng Z-B, Teasdale RD: Multiple interval mapping for quantitative trait loci. Genetics. 1999, 152: 1203-1216.PubMed CentralPubMedGoogle Scholar
- Quevillon E, Silventoinen V, Pillai S, Harte N, Mulder N, Apweiler R, Lopez R: InterProScan: protein domains identifier. Nucleic Acids Res. 2005, 33: W116-W120. 10.1093/nar/gki442.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.