Skip to main content

Genome-wide association analysis for maize stem Cell Wall-bound Hydroxycinnamates



The structural reinforcement of cell walls by hydroxycinnamates has a significant role in defense against pests and pathogens, but it also interferes with forage digestibility and biofuel production. Elucidation of maize genetic variations that contribute to variation for stem hydroxycinnamate content could simplify breeding for cell wall strengthening by using markers linked to the most favorable genetic variants in marker-assisted selection or genomic selection approaches​.


A genome-wide association study was conducted using a subset of 282 inbred lines from a maize diversity panel to identify single nucleotide polymorphisms (SNPs) associated with stem cell wall hydroxycinnamate content. A total of 5, 8, and 2 SNPs were identified as significantly associated to p-coumarate, ferulate, and total diferulate concentrations, respectively in the maize pith. Attending to particular diferulate isomers, 3, 6, 1 and 2 SNPs were related to 8–O–4 diferulate, 5–5 diferulate, 8–5 diferulate and 8–5 linear diferulate contents, respectively. This study has the advantage of being done with direct biochemical determinations instead of using estimates based on Near-infrared spectroscopy (NIRS) predictions. In addition, novel genomic regions involved in hydroxycinnamate content were found, such as those in bins 1.06 (for FA), 4.01 (for PCA and FA), 5.04 (for FA), 8.05 (for PCA), and 10.03 and 3.06 (for DFAT and some dimers).


The effect of individual SNPs significantly associated with stem hydroxycinnamate content was low, explaining a low percentage of total phenotypic variability (7 to 10%). Nevertheless, we spotlighted new genomic regions associated with the accumulation of cell-wall-bound hydroxycinnamic acids in the maize stem, and genes involved in cell wall modulation in response to biotic and abiotic stresses have been proposed as candidate genes for those quantitative trait loci (QTL). In addition, we cannot rule out that uncharacterized genes linked to significant SNPs could be implicated in dimer formation and arobinoxylan feruloylation because genes involved in those processes have been poorly characterized. Overall, genomic selection considering markers distributed throughout the whole genome seems to be a more appropriate breeding strategy than marker-assisted selection focused in markers linked to QTL.


In terms of agricultural land use and production, maize (Zea mays L.) is one of the most important crops worldwide. In addition to the general uses of maize grain as food, feed, or raw material for generating industrial derivatives, maize stover could be a profitable byproduct for ethanol production, whereas the whole plants can also be used to produce silage for feeding cattle.

Accessibility, extensibility, and digestibility of maize stem tissues may determine important characteristics of the maize stem such as maize resistance to stem borer pests and stem diseases, feedstuff quality and degradability and suitability for ethanol production [1]. These characteristics depend greatly on cell wall functionality and structure, which are controlled by the composition and organization of individual cell wall components [1]. Maize cell walls are mainly composed of cellulose embedded in a hemicellulose matrix and lignin. Among cell wall components, cell wall bound hydroxycinnamates, which are derived from the phenylpropanoid pathway and originate from tyrosine and phenylalanine, play a key role in cell wall structural reinforcement. The main hydroxycinnamates in the maize stem are esters of ferulic acid (FA) and p-coumaric acid (PCA), which both contribute to cell wall stiffening and fortification [2]. PCA acts as a radical transfer agent to promote the formation of sinapyl alcohol and lignin radicals being esterified to the Y-position of the side chains of S lignin units [3, 4]. On the other hand, FA could be found within the cell wall esterified to arabinosyl residues of arabinoxylan chains which are thereafter cross-linked via ether bonds to G units of lignins [5, 6]. In addition, it is important to note that FA forms dimers via peroxidase- or laccase-mediated oxidative coupling, and the resulting diferulates (DFAs) cross-links hemicellulose chains (7).

In pest resistance several studies, focused on the study of differences for cell wall components among contrasting materials for resistance to stem borers, have suggested that cell wall strengthening could be a constitutive resistance mechanism to protect maize plants against the attack [7, 8]. The effects of cell wall-bound hydroxycinnamate content on stem borer resistance have been thoroughly investigated, identifying higher concentrations of hydroxycinnamates in resistant inbred lines [9, 10]. In forage species fiber comprises 300–800 μg/g dry matter, and represents the greatest source of energy for ruminants. Unfortunately, less than 50% of this fiber content is actually digested and used by the animal, mainly because of the recalcitrant role of specific cell wall components. The variability in the digestibility of the maize cell wall, which is often greater than 50%, is due, in part, to variability for lignin content, but also to differences for lignin composition itself and cell wall cross-linkage mediated by ferulates and DFAs [11]. For instance, Jung et al. (2011) observed that dry matter intake, and milk production were greater in cows fed with a diet containing W23sfe, a maize mutant with decreased ferulate content, than in cows fed with forage from the wild-type W23. In the same study in vivo digestibility of cell walls was greater on W23sfe [12].

In ethanol production, the two key parts of the maize plant that can be converted to bioethanol are the grain, which is mainly made up of starch, and the debris (stem and leaves), which predominantly contains lignin and cellulosic components. Ethanol produced from non-grain plant material is defined as cellulosic or lignocellulosic ethanol. Lignocellulose in maize is composed of 33.1% hemicellulose, 39.4% cellulose, and 14.9% lignin [13, 14]. The degradation of carbohydrates to the constituent sugar monomers (saccharification) provides the substrates for the fermentation required for yeast-mediated ethanol production. However, the largest cell wall component, cellulose, is not particularly susceptible to deconstruction and is closely interconnected with hemicelluloses and lignin. Cross-linking of lignin to arabinoxylan by hydroxycinnamates makes cell walls highly recalcitrant to biomass degradation [15, 16].

Overall, the mechanical resistance due to higher hydroxycinnamate content and lignification makes maize tissues more recalcitrant to damage by insects, less digestible for ruminants, and less suitable for biofuel production [17,18,19]. The study of the maize functional genetic variability for each hydroxycinnamate component could be crucial to identify relevant genetic variants that may be incorporated into selection programs to breed maize varieties for multiple uses.

Some maize genomic regions mediating hydroxycinnamate accumulation have been detected throughout the genome [17, 20, 21]. Most of these previous studies involved biparental populations and were useful for detecting large regions of interest by analyzing those populations and relatively few markers. Genome-Wide Association Studies (GWAS) enable high-resolution mapping of QTL to narrow genomic regions where the searching for genes contributing to trait variability is feasible. This technique has been extensively used to identify maize genetic variants significantly associated with yield and agronomic traits [22], resistance to biotic and abiotic stresses [23], specific cell wall components such as lignin, cellulose, hemicellulose, detergent fibers, and in vitro digestibility of dry matter [24]. This study constitutes the first high resolution association mapping analysis made to detect QTLs associated to hydroxycinnamate content and has been performed using a diversity panel that contents most genetic variability of maize adapted to temperate areas. Additionally, materials were phenotyped based on wet chemistry determinations, which are more precise than estimations based on Near-Infrared Spectroscopy or Fourier-Transform Infrared Spectroscopy predictions. The objective of this study was to identify genomic regions that could make relevant contributions to the genetic variability for cell wall-bound hydroxycinnamate content in a diverse inbred panel and to propose candidate genes for those regions.


A subset of 270 inbred lines from this diversity panel was evaluated for cell wall-bound hydroxycinnamates. We quantified by liquid chromatography the two main cell wall hydroxcinnamic acids, p-coumaric and ferulic acid, and four ferulic acid dimers (DFA 5–5, DFA 8–O–4, DFA 8–5b, and DFA 8–5 l). As DFA 8–5 l and 8–5b are different isoforms of DFA 8–5, we calculated DFA 8–5 as the sum of both isoforms. Besides, the sum of total dimers (DFAT), total monomers (MONOTOT) and total hydroxycinnamates (FENTOT) were also calculated.

Means, analysis of variance, and Heritabilities

Significant differences were detected among inbred lines for all cell wall-bound hydroxycinnamates (see Additional file 1). Significant genotype-by-environment interactions were observed for individual dimers and DFAT (Data not shown). High heritability values (> 0.70) were estimated for every trait (Table 1).

Table 1 Heritability estimates for the traits under study

Correlation analysis

All correlation (phenotypic and genotypic) coefficients between traits were positive and significant except for the correlation coefficient between DFA 8–5 l and PCA and between DFA 8–5 l and FA. Additionally, PCA was strongly correlated with the total monomer content (MONOTOT) (rg = rp = 0.98) and the total cell wall-bound phenolic acid content (FENTOT) (rg = rp = 0.98). Both correlation coefficients between FA and dimers, as well as the correlation coefficients among the diverse dimers were high and significant (r > 0.70) (Table 2). The strong correlations suggested that MONOTOT and FENTOT are mainly determined by the PCA concentration. Thus, MONOTOT and FENTOT are not further addressed.

Table 2 Genotypic1 (above diagonal) and phenotypic2 (below diagonal) correlation coefficient estimates for each pair of traits

Association analysis

A marker was considered to be significantly associated with a trait at RMIP (Resample Model Inclusion Probability) values more than 0.50. We considered a +/− 150 kbp region around the significant SNP. Two SNPs were assigned to the same QTL when their confidence intervals overlapped. Consequently, 27 SNPs, which corresponded to 22 QTLs, were identified as significantly associated with cell wall-bound hydroxycinnamates, 5 SNPs were associated with PCA, 8 SNPs with FA, 3 SNPs with DFA 8–O–4, 6 SNPs with DFA 5–5, 1 SNP with DFA 8–5, 2 SNPs with DFA 8–5 linear and 2 SNPs with DFAT (Table 3). Among these associations, novel genomic regions involved in hydroxycinnamate content were found in bins 1.06 (for FA), 4.01 (for PCA and FA), 5.04 (for FA), 8.05 (for PCA), and 10.03 and 3.06 (for DFAT and some dimers). Minor and major frequency alleles contributed almost equally to increased PCA levels, while minor frequency alleles generally increased cell wall-bound FA and DFA concentrations. Differences between homozygous for minor and major frequency alleles at each significant SNP varied from 569 to 851 μg/g DW for PCA concentration, from 177 to 321 μg/g DW for FA concentration and from 55 to 56 μg/g DW for DFAT concentration. Differences for dimer concentrations were as follows: 15 to 17 μg/g DW for DFA 8–O–4, 6 to 15 μg/g DW for DFA 5–5, 12 μg/g DW for DFA 8–5 and 5 to 8 μg/g DW for DFA 8–5 l (Table 3). The percentages of variances explained by each significant SNP ranged from 7 to 10%. The significant SNPs found in the current study were distributed in bins 1.05, 1.06, 1.07, 1.11, 3.04, 3.06, 4.01, 5.04, 5.04, 8.05 and 10.03.

Table 3 SNP identification (SNP ID), additive effect and allelic variants for the SNP, proportion of total variance explained by the SNPs significantly associated with cell wall traits (PCA, FA, DFAT, DFA 8–O–4, DFA 8–5 l, DFA 8–5, DFA 5–5), and significance values for the association between the SNP and the phenotype (P-value and RMIP)

Candidate gene selection

The genes containing or physically close to SNPs significantly (RMIP > 0.5) associated with traits were identified and characterized according to the maize B73 reference genome assembly (version 4). Analyses of +/− 150kbp regions surrounding significant SNPs resulted in the identification of 111 genes listed in Table 4. Sixty-three of these genes have not been so far characterized.

Table 4 Complete list of candidate genes for hydroxycinnamatesa identified in a maize diversity panel. Genes in bold are the ones mentioned in the manuscript


Heritabilities and correlations

The high heritabilities observed in the current study are consistent with the results of previous studies [17, 21]. Our data indicate that additive effects are more important than additive × environment interaction effects in the inheritance of hydroxycinnamates in the pith of maize stems, and, consequently, high responses to selection, using any of these traits as selection criteria, would be expected [27]. Based on the evidenced correlation between greater strengthening of the cell wall and the increased hydroxycinnamate content, these traits could be used as indirect selection criteria for the improvement of more complex traits such as pest resistance or forage digestibility [28].

The correlations between traits follow as well the trends reported in the literature [29, 30]. FA, particular dimers and DFAT show co-variation, meanwhile co-variations of PCA and any other hydroxycinnamate compounds were not relevant: Therefore, most genetic variability detected would not be consequence of variations for genes of the common metabolic pathway of these hydroxycinnamates.

QTL co-localization

Significant SNPs were usually located in bins where QTL for hydroxycinnamates (PCA, FA, total diferulates and dimers) were found in previous studies (detailed below). However, novel genomic regions involved in hydroxycinnamate content were also found, such as those in bins 1.06 (for FA), 4.01 (for PCA and FA), 5.04 (for FA), 8.05 (for PCA), and 10.03 and 3.06 (for DFAT and some dimers).

SNPs significantly associated to FA content were found within the supporting intervals of QTLs detected for the same trait by Barriére et al. [29, 31], whereas markers associated with PCA content were included within QTLs for PCA detected by Santiago et al. [21], Barrière et al. [17, 29, 31], and Courtial et al. [32]. Similarly, significant SNPs for DFAT and individual dimers co-localized with QTLs previously published for those traits [31, 33, 34]. It is important to note that previously detected QTLs were mostly identified in bi-parental populations with lower level of resolution. Regarding the correlation coefficients among different traits, we would expect extensive co-localizations between QTLs for FA and individual and total dimers because the genotypic correlation coefficients among those traits were high. However, we only found one genomic region at bin 1.07 where QTLs for FA, DFA 8–O–4, DFA 8–5 l and DFA 5–5 co-localized. This could be explained by the high percentage the unexplained phenotypic variability, factors such as low frequency of minor alleles, small additive effects of genes and/or low density of markers could be responsible.

In the same regions, QTLs for resistance to pests [23, 35,36,37], animal degradability [24, 29, 38] and fuel production [29, 31] have already been reported (See Adittional File 2), supporting that the variation for cell wall bound hydroxycinnamates could have a significant role on these three aspects. Marker-assisted breeding programs to indirectly improve maize resistance to insect attack could focus on genomic regions where QTLs for hydroxycinnamates across different populations are found, as well as QTLs for increased resistance to stem borers; especially, if associated unfavorable effects for biofuel production and/or animal digestibility are not found as in bins 1.06 and 3.06. However, the lack of significant co-localizations does not mean that linked and/or pleiotropic genes do not exist in that region because detection power in the present study would not be enough to uncover all genomic regions involved in any trait. Therefore, we propose to implement genomic selection for increased cell-wall hydroxycinnamate content in a population with a high genetic diversity in order to better establish indirect effects of selection for increased cell-wall strength on resistance to stem borers, animal digestibility and/or biofuel production. In addition, effects on other cell-wall components such as lignin content or composition that also have impact on cell-wall strengthening would be desirable.

Candidate genes

Key genes implicated in monolignols and p-hydroxycinnamic acid biosynthesis are known, and none of them were located at distances shorter than +/− 0.15 Mbp from SNPs significantly associated to hydroxycinnamates. Even though, little is known about the genes implicated in dimerization by peroxidase-mediated oxidative coupling of FA residues to form ferulic acid oligomers [39, 40]. Similarly, genes involved in feruloylation of arabonoxylans (AX) are poorly described; there are some hypothesis that propose CoA-acyl transferases in the PF02458 family as good candidates for feruloylation as the modification of CoA-acyl transferase genes changes the level of esterified ferulates and diferulates in the cell wall, although their roles are inconclusive [41, 42]. Recent studies describe that suppression of a single BAHD gene in Setaria viridis causes large, stable decreases in cell wall feruloylation and increases biomass digestibility [43]. However, we cannot rule out that some of the uncharacterized candidate genes and hypothetical proteins noted in the current study could be genes involved in the dimerization, specifically those located in bins 3.06 and 10.03, where significant SNPs for dimers have been found but not for FA, while uncharacterized candidate genes in bins 1.07 and 1.11 that contain significant SNPs for FA and dimer contents could be implicated in feruloylation of arabinoxylans.

Among the genes found in the confidence interval of QTL significantly associated with PCA, we highlight three of them: two possible 4-coumarate ligases (Zm00001d031090, Zm00001d034486) and a cell wall kinase (Zm00001d039926). The 4-coumarate ligase is a key enzyme involved in the general phenylpropanoid metabolism, that catalyzes the conversion of 4-coumaric acid into 4-coumaroyl-CoA leading to the formation of precursors that serve as structural components for the biosynthesis of cell wall associated phenolics [44].

Cell wall kinases (WAKs) are members of the RLCKs (Receptor-like cytoplasmic kinase) super family which are one of the main candidates for sensing and controlling cell wall integrity [45, 46]. Some WAKs display carbohydrate binding activity interacting with the cell wall probably via pectin binding, which appears to influence a WAK-dependent signaling pathway regulating cell expansion and also involved in signaling during stress response and pathogen infection resistance. A WAK gene is proposed as the best candidate gene for the QTL for PCA at bin 3.04 although the specific function of this gene needs to be tested.

Genes annotated as probably involved in glycerolipid metabolism within the supporting intervals of four QTLs for FA were also found: a glycerol-3-phosphate acyltransferase 1 (Zm00001d034732), an Epoxide hydrolase (Zm00001d034740), a Monogalactosyldiacylglycerol synthase 2 (Zm00007a00030071), a glycerophosphodiester phosphodiesterase (Zm00001d048970) and a digalactosyldiacylglycerol synthase 1 (Zm00001d016731). These enzymes could be also be involved in biosynthesis of steryl glycosides and acylated steryl glycosides [47]. Increased acylated steryl glycosides has been associated to matured fiber cells compared to elongating fiber cells [48]. Total FA content was highest or increased during cell elongation and was lower or decreased thereafter [48, 49], likely reflecting the assembly of the expanding cell membranes during elongation and the shift to membrane maintenance (and increase in secondary cell wall content) in maturing fibers [48]. Moreover, for the QTL for FA, DFA 8–O–4, DFA 8–5 l and DFA 5–5 in bin 1.07, we propose an UDP-glucose 4-epimerase 4 (UGE4) gene (Zm00001d032346) as probable candidate gene. In Arabidopsis, UGE4 participates in cell carbohydrate biosynthesis and growth by catalyzing the interconversion between UDP-glucose and UDP-galactose. This enzyme is involved in channeling UDP-D-galactose into cell wall polymers which is required for the galactosylation of xyloglucan (XyG) and type II arabinogalactan (AGII) [50]. Finally, we propose genes encoding for a polygalacturonase (Zm00001d034727) for the QTL for FA at 1.11 and two genes encoding pectinesterase inhibitors (Zm00001d016738 and Zm00001d016741) for QTL in bin 5.04. These enzymes are involved in physiological processes such as modification of different cell wall polymers and cell wall modulation processes during development and stress responses [51, 52]. In relation to DFAT, we propose five genes as the most probable candidates out of the ones found in the confidence intervals of QTL. We highlight an aldose-1-epimerase gene as probable candidate gene for QTLs for DFA 8–5 l and DFA 8–5 at bin 1.11 because accumulation of the aldose-1-epimerase has been related to cell wall modification, rearrangement and stiffening [53]. For the QTL at bin 3.04, we propose a gene enconding a Gibbellerin 2-oxidase (Zm00001d040737) that irreversibly catalyze the deactivation of bioactive gibbellerin because Wuddineh et al. [54] demonstrated that lignification and biomass recalcitrance could be optimized by targeting gibberellin biosynthesis. We also suggest a candidate gene that encodes an endoplasmic reticulum transporter vesicle protein (Zm00001d043052) for the QTL in bin 3.06. The dynamic structure of the cell wall is controlled by polysaccharide modification, and the pectic and non-cellulosic polysaccharide constituents of plant cell walls are made within endoplasmic reticulum (ER)-Golgi apparatus and exported to the cell surface [55]. Finally, a polygacturonase gene is proposed as candidate gene for the QTL at bin 10.3 because its implication in cell wall loosening during growth and development by the modification and/or breakdown of different polymers [56].


Overall, new and known genomic regions associated with the accumulation of cell wall-bound hydroxycinnamic acids in the maize pith were revealed, and could have impact on hydroxycinnamate content across different genetic backgrounds. Genes involved in cell wall modulation were proposed as candidate genes for cell wall-bound hydroxycinnamate accumulation. However, we cannot rule out that some others uncharacterized genes linked to significant SNPs could be implicated in dimers formation and arobinoxylan feruloylation.

The effects of individual SNPs significantly associated with hydroxycinnamate content were low, and each SNP explained a low percentage of total genetic variability. Therefore, the most appropriate marker-based breeding strategy to increase hydroxycinnamate content would be the genomic selection, because markers evenly distributed throughout the entire genome would better estimate a breeding value than few markers linked to QTL with a modest contribution to the genetic variability.


Plant materials and experimental design

Most diversity available in the worldwide public breeding sector is represented in the diversity panel, comprising 302 inbred lines, evaluated [57, 58]. A subset of 270 inbred lines from this diversity panel was evaluated for cell wall-bound hydroxycinnamates (mostly due to seeds germinations problems). Field trials using 18 × 15 α-lattice experimental designs with two replicates were performed in Pontevedra, Spain, in 2010 and 2011. Plot characteristics and agronomic practices were as described by Samayoa et al. [23]. The plant material was obtained from the North Central Regional Plant Introduction Station from Iowa State University following the seeds import law applicable for Spain. The list of inbred lines used in this study along with their accession identifications at the North Central Regional Plant Introduction Station (NCRPIS) is available as additional file in Samayoa et al. [59].

Biochemical determinations

The second internodes below the main ear from five plants per plot were collected 30 days after silking. Samples were frozen at − 20 °C. The pith was manually detached, lyophilized, and ground in a Wiley mill with a 0.75 mm screen. Ground pith samples were maintained at 5 °C until biochemical analyses. Hydroxycinnamates quantification was done following a recently optimized protocol by Santiago et al. [60]. The identities of FA dimers were confirmed by a comparison with the authentic 5–5 standard or published retention times and UV spectra [61]. The total diferulate content (DFAT) was calculated as the sum of the following three identified and quantified DFA isomers: DFA 8–O–4, DFA 5–5-, and DFA 8–5. The DFA 8–5 concentration was calculated as the sum of 8–5-cyclic (or benzofuran)-DFA and 8–5-noncyclic (or open). Additionally, MONOTOT represents the sum of all monomers (PCA and FA), while FENTOT refers to the sum of all hydroxycinnamate (monomers and dimers).

Statistical analysis

Inbred lines were previously genotyped with a unique set of SNP markers derived from the Illumina maize 50 k array and a genopying-by-sequencing (GBS) strategy (90, 91). The two resulting datasets were combined. The SNPs with more than 20% missing data and a minor allele frequency less than 5% were omitted. Heterozygous genotypes were considered missing data. After filtering, 246,477 SNPs distributed across the maize genome were retained [23, 62, 63].

Best linear unbiased estimator

Each trial was analyzed separately according to the mixed model procedure (PROC MIXED) of the SAS program (version 9.4) [64] and the best linear unbiased estimators for each inbred line were calculated based on the combined data for the 2-year analysis. Inbred lines were considered as the fixed effect, while years, replicates and blocks within replicates were random effects.

Heritabilities and correlations

Heritabilities (ĥ2) were estimated for each trait on a family-mean basis as previously described [65]. The genetic (rg) and phenotypic (rp) correlation coefficients between traits were calculated using REML estimates according to a published SAS mixed model procedure [66].

Association analysis

A genome-wide association analysis was completed with Tassel 5 [67] based on a mixed linear model using a genotype-phenotype matrix and a kinship matrix [26] as a covariate. Among the mixed linear model options, we used the optimum compression level and P3D to estimate the variance component. An established subsampling method [23] was applied to identify SNPs with the most robust associations (RMIP).

The genetic kinship matrix used for this GWAS was previously published [26], and was estimated using a subset of 5000 SNPs (with no missing genotypes) that were distributed almost evenly across the whole maize genome.

Candidate gene selection

Following previous studies performed with the same genetic material [23] and the LD decay observed in the regions surrounding significant SNPs, we considered a +/− 150 kbp confidence interval region around each significant SNP. LD between SNPs separated by more than 150 Kbp is inappreciable (r2 < 0.1) in at least 90% of cases [62]. In case confidence intervals of two SNPs overlapped, they were assigned to a single QTL. The two described genes that delimit the 300 kbp region around the SNP in the reference genome assembly version 2 were positioned in version 4 of the reference genome, and all genes contained in the region delimited by those genes were then identified and characterized based on the maize B73 reference genome assembly (version 4) available on the MaizeGDB browser [25].

Availability of data and materials

The data sets used and/or analysed during the current study will be available upon reasonable request to the corresponding author. Vegetal materials are distributed to the scientific community by the NCRPIS upon request (https://.


DFA 5–5:

5–5 diferulic acid

DFA 8–5:

sum (DFA8–5 l + DFA8–5b)

DFA 8–5 l:

8–5 linear diferulic acid

DFA 8–5b:

8–5 benzofuran diferulic acid

DFA 8–O–4:

8–O–4 diferulic acid


total diferulates (DFA 5–5 + DFA 8–O–4 + DFA 8–5b + DFA 8–5 l)


ferulic acid


total hydroxycinnamates (PCA+ FA+ DFAT)


Genome-Wide Association


Linkage disequilibrium


Mediterranean Corn Borer


total monomers (PCA + FA)


p-coumaric acid


Single Nucleotide Polimorfism


  1. Fry SC. Cross-linking of matrix polymers in the growing cell walls of angiosperms. Annu Rev Plant Physiol. 1986;37(1):165–86.

    Article  CAS  Google Scholar 

  2. Grabber JH, Ralph J, Lapierre C, Barrière Y. Genetic and molecular basis of grass cell-wall degradability. I. Lignin-cell wall matrix interactions. Comptes Rendus - Biol. 2004;327(5):455–65.

    Article  CAS  Google Scholar 

  3. Wolf DP, Coors JG, Albrecht KA, Undersander DJ, Carter PR. Forage quality of maize genotypes selected for extreme fiber concentrations. Crop Sci. 1993;33(6):1353–9.

    Article  Google Scholar 

  4. Grabber JH, Quideau S, Ralph J. P-coumaroylated syringyl units in maize lignin: implications for β-ether cleavage by thioacidolysis. Phytochemistry. 1996;43(6):1189–94.

    Article  CAS  Google Scholar 

  5. Bunzel M, Ralph J, Steinhart H. Phenolic compounds as cross-links of plant derived polysaccharides. Czech J Food Sci. 2004;22(6):39–42.

    Google Scholar 

  6. Hatfield RD, Marita JM. Enzymatic processes involved in the incorporation of hydroxycinnamates into grass cell walls. Phytochem Rev. 2010;9(1):35–45.

    Article  CAS  Google Scholar 

  7. Malvar RA, Butrón A, Ordás B, Santiago R. Causes of natural resistance to stem borers in maize. E.N. Burton, P.V. Williams (Eds.) Crop Prot Res Adv Nova Science Publishers, Inc. 2008;57–100.

  8. Meihls LN, Higdon ML, Ellersieck MR, Tabashnik BE, Hibbard BE. Greenhouse-Selected Resistance to Cry3Bb1-Producing Corn in Three Western Corn Rootworm Populations. PLoS One. 2012;7(12):e51055.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Santiago R, Butrón A, Reid LM, Arnason JT, Sandoya G, Souto XC, et al. Diferulate content of maize sheaths is associated with resistance to the Mediterranean corn borer Sesamia nonagrioides (Lepidoptera: Noctuidae). J Agric Food Chem. 2006;54(24):9140–4.

    Article  CAS  PubMed  Google Scholar 

  10. Barros-Rios J, Santiago R, Jung HJG, Malvar RA. Covalent cross-linking of cell-wall polysaccharides through esterified diferulates as a maize resistance mechanism against corn borers. J Agric Food Chem. 2015;63(8):2206–14.

    Article  CAS  PubMed  Google Scholar 

  11. Barrière Y, Laperche A, Barrot L, Aurel G, Briand M, Jouanin L. QTL analysis of lignification and cell wall digestibility in the Bay-0 x Shahdara RIL progeny of Arabidopsis thaliana as a model system for forage plant. Plant Sci. 2005;168(5):1235–45.

    Article  CAS  Google Scholar 

  12. Jung HG, Mertens DR, Phillips RL. Effect of reduced ferulate-mediated lignin/arabinoxylan cross-linking in corn silage on feed intake, digestibility, and milk production. J Dairy Sci. 2011;94(10):5124–37.

    Article  CAS  PubMed  Google Scholar 

  13. Chundawat SPS, Venkatesh B, Dale BE. Effect of particle size based separation of milled corn Stover on AFEX pretreatment and enzymatic digestibility. Biotechnol Bioeng. 2007;96(2):219–31.

    Article  CAS  PubMed  Google Scholar 

  14. Pauly M, Keegstra K. Cell-wall carbohydrates and their modification as a resource for biofuels. Plant J. 2008;54(4):559–68.

    Article  CAS  PubMed  Google Scholar 

  15. Torres AF, Visser RGF, Trindade LM. Bioethanol from maize cell walls: Genes, molecular tools, and breeding prospects. GCB Bioenergy. 2015;7:591–607.

    Article  CAS  Google Scholar 

  16. Courtial A, Soler M, Chateigner-Boutin A-L, Reymond M, Mechin V, Wang H, et al. Breeding grasses for capacity to biofuel production or silage feeding value: an updated list of genes involved in maize secondary cell wall biosynthesis and assembly. Maydica. 2013;58(1):67–102.

    Google Scholar 

  17. 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(4):585–95.

    Article  CAS  Google Scholar 

  18. Ralph J. Hydroxycinnamates in lignification. Phytochem Rev. 2010;9(1):65–83.

    Article  CAS  Google Scholar 

  19. Barros-Rios J, Santiago R, Malvar RA, Jung HJG. Chemical composition and cell wall polysaccharide degradability of pith and rind tissues from mature maize internodes. Anim Feed Sci Technol. 2012;172(3–4):226–36.

    Article  CAS  Google Scholar 

  20. Riboulet C, Fabre F, Dénoue D, Martinantä JP, Lefèvre B, Barrière Y. QTL mapping and candidate gene research for lignin content and cell wall digestibility in a top-cross of a flint maize recombinant inbred line progeny harvested at silage stage. Maydica. 2008;53(1):1–9.

    Google Scholar 

  21. Santiago R, Malvar RA, Barros-Rios J, Samayoa LF, Butron A. Hydroxycinnamate synthesis and association with Mediterranean corn borer resistance. J Agric Food Chem. 2016;64(3):539–51.

    Article  CAS  PubMed  Google Scholar 

  22. Suwarno WB, Pixley KV, Palacios-Rojas N, Kaeppler SM, Babu R. Genome-wide association analysis reveals new targets for carotenoid biofortification in maize. Theor Appl Genet. 2015;128(5):851–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Samayoa LF, Malvar RA, Olukolu BA, Holland JB, Butrón A. Genome-wide association study reveals a set of genes associated with resistance to the Mediterranean corn borer (Sesamia nonagrioides L.) in a maize diversity panel. BMC Plant Biol. 2015;15(1):35.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  24. Wang H, Li K, Hu X, Wu Y, Huang C. Genome-wide association analysis of forage quality in maize mature stalk. BMC Plant Biol. 2016;16:227.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. Cho KT, Portwood J, Harper LC, Gardiner JM, Lawrence-Dill CJ, Friedberg I, Andorf CM. (2019) MaizeDIG: Maize Database of Images and Genomes. Frontiers in Plant Science. 2019.

  26. Olukolu BA, Negeri A, Dhawan R, Venkata BP, Sharma P, Garg A, et al. A connected set of genes associated with programmed cell death implicated in controlling the hypersensitive response in maize. Genetics. 2013;193(2):609–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Santiago R, Sandoya G, Butrón A, Barros J, Malvar RA. Changes in phenolic concentrations during recurrent selection for resistance to the Mediterranean corn borer (Sesamia nonagrioides Lef.). J Agric Food Chem. 2008;56(17):8017–22.

    Article  CAS  PubMed  Google Scholar 

  28. Barros-Rios J, Malvar RA, Jung HJG, Bunzel M, Santiago R. Divergent selection for ester-linked diferulates in maize pith stalk tissues. Effects on cell wall composition and degradability. Phytochemistry. 2012;83:43–50.

    Article  CAS  PubMed  Google Scholar 

  29. Barrière Y, Méchin V, Riboulet C, Guillaumie S, Thomas J, Bosio M, et al. Genetic and genomic approaches for improving biofuel production from maize. Euphytica. 2009;170(1):183–202.

    Article  CAS  Google Scholar 

  30. Piston F, Uauy C, Fu L, Langston J, Labavitch J, Dubcovsky J. Down-regulation of four putative arabinoxylan feruloyl transferase genes from family PF02458 reduces ester-linked ferulate content in rice cell walls. Planta. 2010;231(3):677–91.

    Article  CAS  PubMed  Google Scholar 

  31. Barriere Y, Courtial A, Chateigner-Boutin A-L, Denoue D, Grima-Pettenati J. Breeding maize for silage and biofuel production, an illustration of a step forward with the genome sequence. Plant Sci. 2016;242:310–29.

    Article  CAS  PubMed  Google Scholar 

  32. Courtial A, Méchin V, Reymond M, Grima-Pettenati J, Barrière Y. Colocalizations between several QTLs for Cell Wall degradability and composition in the F288 × F271 early maize RIL progeny raise the question of the nature of the possible underlying determinants and breeding targets for biofuel capacity. Bioenergy Res. 2014;7(1):142–56.

    Article  Google Scholar 

  33. García-Lara S, Burt AJ, Arnason JT, Bergvinson DJ. QTL mapping of tropical maize grain components associated with maize weevil resistance. Crop Sci. 2010;50(3):815–25.

    Article  CAS  Google Scholar 

  34. 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(3):531–49.

    Article  PubMed  CAS  Google Scholar 

  35. Ordas B, Malvar RA, Santiago R, Sandoya G, Romay MC, Butron A. Mapping of QTL for resistance to the Mediterranean corn borer attack using the intermated B73 × Mo17 (IBM) population of maize. Theor Appl Genet. 2009;119(8):1451–9.

    Article  CAS  PubMed  Google Scholar 

  36. Krakowsky MD, Lee M, Woodman-Clikeman WL, Long MJ, Sharopova N. QTL mapping of resistance to stalk tunneling by the European corn borer in RILs of maize population B73 X De811. Crop Sci. 2004;44(1):274–82.

    CAS  Google Scholar 

  37. Khairallah MM, Jiang C, González-De-Leon D, Hoisington DA, Bohn M, Melchinger AE, et al. Molecular mapping of QTL for southwestern corn borer resistance, plant height and flowering in tropical maize. Plant Breed. 1998;117(4):309–18.

    Article  Google Scholar 

  38. Meng Y, Li J, Liu J, Hu H, Li W, Liu W, et al. Ploidy effect and genetic architecture exploration of stalk traits using DH and its corresponding haploid populations in maize. BMC Plant Biol. 2016;16(1):1–15.

    Article  CAS  Google Scholar 

  39. Saha BC. Hemicellulose bioconversion. J Ind Microbiol Biotechnol. 2003;30:279–91.

    Article  CAS  PubMed  Google Scholar 

  40. Mosier N, Wyman C, Dale B, Elander R, Lee YY, Holtzapple M, et al. Features of promising technologies for pretreatment of lignocellulosic biomass. Bioresour Technol. 2005;96(6):673–86.

    Article  CAS  PubMed  Google Scholar 

  41. Buanafina MMO, Langdon T, Hauck B, Dalton S, Timms-Taravella E, Morris P. Targeting expression of a fungal ferulic acid esterase to the apoplast, endoplasmic reticulum or golgi can disrupt feruloylation of the growing cell wall and increase the biodegradability of tall fescue (Festuca arundinacea). Plant Biotechnol J. 2010;8(3):316–31.

    Article  CAS  PubMed  Google Scholar 

  42. Bartley LE, Peck ML, Kim S-R, Ebert B, Manisseri C, Chiniquy DM, et al. Overexpression of a BAHD Acyltransferase, OsAt10, alters Rice Cell Wall Hydroxycinnamic acid content and Saccharification. Plant Physiol. 2013;161(4):1615–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. de Souza WR, Martins PK, Freeman J, Pellny TK, Michaelson LV, Sampaio BL, et al. Suppression of a single BAHD gene in Setaria viridis causes large, stable decreases in cell wall feruloylation and increases biomass digestibility. New Phytol. 2018;218(1):81–93.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  44. Stuible HP, Büttner D, Ehlting J, Hahlbrock K, Kombrink E. Mutational analysis of 4-coumarate:CoA ligase identifies functionally important amino acids and verifies its close relationship to other adenylate-forming enzymes. FEBS Lett. 2000;467(1):117–22.

    Article  CAS  PubMed  Google Scholar 

  45. He Z-H. A cluster of five cell wall-associated receptor kinase genes, Wak1–5, are expressed in specific organs of Arabidopsis. Plant Mol Biol. 1999;39:1189–96.

  46. Kohorn BD, Kohorn SL. The cell wall-associated kinases, WAKs, as pectin receptors. Front Plant Sci. 2012;3(May):1–5.

    Google Scholar 

  47. Tevini M, Lichtenthaler HK. Lipids and lipid polymers in higher plants. Springer. 2012.

  48. Wanjie SW, Welti R, Moreau RA, Chapman KD. Identification and quantification of glycerolipids in cotton fibers: reconciliation with metabolic pathway predictions from DNA databases. Lipids. 2005;40(8):773–85.

    Article  CAS  PubMed  Google Scholar 

  49. MacAdam JW, Grabber JH. Relationship of growth cessation with the formation of diferulate cross-links and p-coumaroylated lignins in tall fescue leaf blades. Planta. 2002;215(5):785–93.

    Article  CAS  PubMed  Google Scholar 

  50. Andème-Onzighi C, Sivaguru M, Judy-March J, Baskin TI, Driouich A. The reb1-1 mutation of Arabidopsis alters the morphology of trichoblasts, the expression of arabinogalactan-proteins and the organization of cortical microtubules. Planta. 2002;215(6):949–58.

    Article  PubMed  CAS  Google Scholar 

  51. Minic Z. Physiological roles of plant glycoside hydrolases. Planta. 2008;227(4):723–40.

    Article  CAS  PubMed  Google Scholar 

  52. Ferrari S, Sella L, Janni M, De Lorenzo G, Favaron F, D’Ovidio R. Transgenic expression of polygalacturonase-inhibiting proteins in Arabidopsis and wheat increases resistance to the flower pathogen Fusarium graminearum. Plant Biol. 2012;14(SUPPL. 1):31–8.

    Article  CAS  PubMed  Google Scholar 

  53. Holden HM, Rayment I, Thoden JB. Structure and function of enzymes of the Leloir pathway for Galactose metabolism. J Biol Chem. 2003;278(45):43885–8.

    Article  CAS  PubMed  Google Scholar 

  54. Wuddineh WA, Mazarei M, Zhang J, Poovaiah CR, Mann DGJ, Ziebell A, et al. Identification and overexpression of gibberellin 2-oxidase (GA2ox) in switchgrass (Panicum virgatum L.) for improved plant architecture and reduced biomass recalcitrance. Plant Biotechnol J. 2015;13(5):636–47.

    Article  CAS  PubMed  Google Scholar 

  55. Nebenführ A, Staehelin LA. Mobile factories: Golgi dynamics in plant cells. Trends Plant Sci. 2001;6(4):160–7.

    Article  PubMed  Google Scholar 

  56. Bilska-Kos A, Solecka D, Dziewulska A, Ochodzki P, Jończyk M, Bilski H, et al. Low temperature caused modifications in the arrangement of cell wall pectins due to changes of osmotic potential of cells of maize leaves (Zea mays L.). Protoplasma. 2017;254(2):713–24.

    Article  CAS  PubMed  Google Scholar 

  57. Liu K, Goodman M, Muse S, Smith JS, Buckler E, Doebley J. Genetic structure and diversity among maize inbred lines as inferred from DNA microsatellites. Genetics. 2003;165(4):2117–28.

    CAS  PubMed  PubMed Central  Google Scholar 

  58. Flint-Garcia SA, Thuillet AC, Yu J, Pressoir G, Romero SM, Mitchell SE, et al. Maize association population: a high-resolution platform for quantitative trait locus dissection. Plant J. 2005;44(6):1054–64.

    Article  CAS  PubMed  Google Scholar 

  59. Samayoa LF, Cao A, Santiago R, Malvar RA, Butrón A. Genome-wide association analysis for fumonisin content in maize kernels. BMC Plant Biol. 2019;19(1):1–11.

    Article  Google Scholar 

  60. Santiago R, López-Malvar A, Souto C, Barros-Ríos J. Methods for determining Cell Wall-bound Phenolics in maize stem tissues. J Agric Food Chem. 2018;66(5):1279–84.

    Article  CAS  PubMed  Google Scholar 

  61. Waldron KW, Parr AJ, Ng A, Ralph J. Cell wall esterified phenolic dimers: identification and quantification by reverse phase high performance liquid chromatography and diode array detection. Phytochem Anal. 1996;7(6):305–12.

    Article  CAS  Google Scholar 

  62. Romay MC, Millard MJ, Glaubitz JC, Peiffer JA, Swarts KL, Casstevens TM, et al. Comprehensive genotyping of the USA national maize inbred seed bank. Genome Biol. 2013;14(6):R55.

  63. Ganal MW, Durstewitz G, Polley A, Bé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(12):e28334.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. SAS/STAT; SAS Institute Inc.: Cary, NC, 2007.

  65. Holland JB, Nyquist WE, Cervantes-Martínez CT. Estimated and interpreting heritability for plant breeding: an update. Plant Breed Rev. 2003;22:9–122.

    Google Scholar 

  66. Holland JB. Estimating genotypic correlations and their standard errors using multivariate restricted maximum likelihood estimation with SAS Proc MIXED. Crop Sci. 2006;46(2):642–54.

    Article  Google Scholar 

  67. Bradbury PJ, Zhang Z, Kroon DE, Casstevens TM, Ramdoss Y, Buckler ES. TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics. 2007;23(19):2633–5.

    Article  CAS  PubMed  Google Scholar 

Download references


We thank Edanz Group (www,edanzediting,com/ac) for editing a draft of this manuscript. We thank Ana Carballeda for technical assistance in laboratory experimentations.


This research has been developed in the frame of the ‘Agri-Food Research and Transfer Centre of the Water Campus (CITACA) at the University of Vigo (Spain), which is economically supported by the Galician Government. It was funded by the “Plan Estatal de Ciencia y Tecnología de España” (project RTI2018–096776-B-C21, and RTI2018–096776-B-C22co-financed with European Union funds under the FEDER program) and by “Xunta de Galicia” (project ED431F 2016/014). A. López-Malvar’s scholarship for the PhD fulfillment has been granted by University of Vigo and R. Santiago’s postdoctoral contract “Ramón y Cajal” has been financed by the Ministry of Economy and Competiveness (Spain), Vigo University, and the European Social Fund. The funding body played no role in study design, data analysis, and manuscript preparation.

Author information

Authors and Affiliations



AB, RAM and RS conceived the study, participated in its design, carried out particular analyses, and assisted in data analysis; AL and LFS performed data analysis and wrote the manuscript; and DJF assisted in laboratory analysis. All authors read and approved the final manuscript.

Corresponding author

Correspondence to A. López-Malvar.

Ethics declarations

Ethics approval and consent to participate

Not applicable in this study.

Consent for publication

Not applicable in this study.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1:

Means of 270 association-panel inbreds for contents of cell-wall bound hydroxycinnamates evaluated in two years.

Additional file 2:

Colocalizations between genomic regions associated with cell-wall bound hydroxycinnamates in the current study and published QTLs for cell-wall traits, resistance to pests, forage digestibility and bioethanol production.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

López-Malvar, A., Butrón, A., Samayoa, L.F. et al. Genome-wide association analysis for maize stem Cell Wall-bound Hydroxycinnamates. BMC Plant Biol 19, 519 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: