Skip to main content
  • Research article
  • Open access
  • Published:

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

Abstract

Background

Corn borers are the primary maize pest; their feeding on the pith results in stem damage and yield losses. In this study, we performed a genome-wide association study (GWAS) to identify SNPs associated with resistance to Mediterranean corn borer in a maize diversity panel using a set of more than 240,000 SNPs.

Results

Twenty five SNPs were significantly associated with three resistance traits: 10 were significantly associated with tunnel length, 4 with stem damage, and 11 with kernel resistance. Allelic variation at each significant SNP was associated with from 6 to 9% of the phenotypic variance. A set of genes containing or physically close to these SNPs are proposed as candidate genes for borer resistance, supported by their involvement in plant defense-related mechanisms in previously published evidence. The linkage disequilibrium decayed (r 2 < 0.10) rapidly within short distance, suggesting high resolution of GWAS associations.

Conclusions

Most of the candidate genes found in this study are part of signaling pathways, others act as regulator of expression under biotic stress condition, and a few genes are encoding enzymes with antibiotic effect against insects such as the cystatin1 gene and the defensin proteins. These findings contribute to the understanding the complex relationship between plant-insect interactions.

Background

Corn borers are the primary maize pest in many environments [1,2]. Corn borers feeding on the pith of the stem results in yield losses because stem damage interferes with assimilate movement to developing kernels. They can also attack the ears, promoting secondary fungal infection, leading to contamination of grain with mycotoxins that may affect human and animal health [3,4].

There are different species of borers that attack maize in different parts of the world. The most economically important species are classified into two families: Crambidae and Noctuidae. Within the Crambidae family, the species with economic importance are: Ostrinia nubilalis Hubner in North America, Europe and North Africa; Ostrinia furnacalis Guenée in Asia; Diatraea saccharalis Fabricius from USA to Argentina; Chilo partellus Swinhoe in Southern USA, Central America and the Caribbean; Diatraea lineolata Walker in Central America, the Caribbean region and South America; and Diatraea grandiosella Dyar in North and Central America. Within the Noctuidae family, the main maize borers are: Sesamia nonagrioides Lefebvre in the Mediterranean region, Busseola fusca Fuller in sub-Saharan Africa, and Sesamia calamistis Hampson in West Africa. This study focuses on the noctuid Mediterranean corn borer (MCB) Sesamia nonagroides Lefebvre, the most important insect pest of maize in the Mediterranean region that includes Southern Europe [2,5,6].

The use of transgenic corn which produces Bacillus thutingiensis (Bt) toxins is a good method for controlling these pests, but transgenic crops are not authorized in several European countries under any agricultural system [7] and are not allowed for organic production [8]. In addition, recent studies have reported a reduction of efficacy of Bt transgenes caused by evolved resistance of some important pests [9-11]. The stacking of several resistant genes has been proposed as one means to delay insect adaptation [12]. Natural sources of resistance to stem borers in maize could reveal promising genes for use in either breeding or transgenic approaches to resistance.

In Spain, there are three MCB generations per year and the second and subsequent generations are able of making significant damage on the stem and, secondarily, on the ear. Following artificial infestation, the level of maize host resistance to stem borer is measured by the tunnel length made by the larvae in the stem as well as by a visual scale for kernel damage. These traits have a complex genetic architecture because resistance depends on the plant-insect relationship, which is influenced by environmental conditions and the developmental stage of the host plant [13,14]. The line mean heritability estimates for tunnel length under corn borer infestation ranged among studies from moderate to relatively high (h 2 = 0.50 to 0.78) [15-18], depending on the genetic background.

At present, commercial materials with high levels of native resistance to these insects are not available even though breeding for increasing maize resistance to corn borers has been conducted during the last three decades in different regions around the world. Klenke et al. [19] reduced tunnel length by 4 and 6 cm after four cycles of recurrent selection for resistance to the first and second generations of ECB, respectively. Bosque-Pérez et al. [20] described successful results in the development of materials with resistance to S. calamistis and African sugarcane borer (Eldana saccharina Walker). Sandoya et al. [21] achieved a reduction of 1.8 cm per cycle of recurrent selection for tunnel length by MCB. In addition, a negative relationship between resistance to stem borers and yield has been found [22] when selection for improved yield under infestation was practiced [23]. In summary, classical breeding experiments have demonstrated some successful improvement in corn borer resistance, but natural levels of resistance in elite cultivars remain insufficient to manage the pest. Detection of stem borer resistance QTL could enhance breeding for this trait via marker-assisted breeding or genomic selection.

Genetic effects for resistance against borer attack fit an additive – dominant model, although additive effects appear to be the most important in determining resistance to tunnel length and kernel damage [14,24-28]. Therefore, the study of genetic factors involved in maize resistance to borers can be performed using highly inbred lines.

Several studies performed with segregating biparental populations have reported genomic regions containing minor and major quantitative traits loci (QTL) for resistance to stem and leaf attack by European corn borer (ECB; [15-17,29-31]. Fewer studies have reported QTLs for resistance to other borer species such as the Mediterranean corn borer (MCB) [18,32] or some tropical borers [33-35]. In maize, QTLs for resistance to tunneling by corn borers have been detected on all chromosomes, with the most commonly detected regions occurring on chromosomes 1, 2, 3, 5, and 9 [1]. Identification of causal genes underlying the QTL for resistance in these genomic regions could help breeders to transfer the allelic variants that confer resistance from the lines that carry them to elite breeding lines that lack these resistance alleles. The identification of the causal genetic variants or markers in high LD with causal variants in diverse materials would minimize the risk of dragging other genes with negative effect on the agronomic value during the transfer process.

Although conventional QTL mapping based on linkage maps of biparental populations has been an efficient approach to detect regions related to the resistance to corn borers, higher resolution is needed to detect the genes involved in the defense mechanism of the plant. Genome wide association study (GWAS) based on linked disequilibrium (LD) in diverse genetic samples is a relatively new approach which offers higher resolution mapping that under optimal conditions can pinpoint causal genes underlying quantitative trait variation. Exploiting advances in genotyping and sequencing technology, this approach has been successful in detecting genes associated with diseases in humans [36-38], animals [39-42] and different quantitative traits in plants [43-46]. In contrast to the conventional QTL mapping approach based on linkage in a biparental population, GWAS is based on LD among extant lines from different populations, such that a large number of markers covering the whole genome are required [45,47]. In diverse maize samples, LD is low, therefore, many more markers are needed than in autogamous species (with higher LD) to adequately explore the genetic architecture of complex traits [48]. The low LD offers the benefit of better resolution to delineate potential causal genes within small LD blocks.

Many single-nucleotide polymorphism (SNP) markers have been identified and scored on a maize diversity panel (composed of 302 inbred line) that represents the diversity available in public breeding sector around the world [49,50]. The population has been successfully used by the maize community to perform GWAS in economically important quantitative traits such as kernel composition [51], hypersensitive response [52] and Fusarium ear rot resistance [53].

To date, no GWAS for insect resistance in maize has been reported, although a few GWA studies that deal with plant defense mechanisms against insect attack have been reported in other plant species [54,55]. In this study, we performed GWAS to identify SNPs associated with resistance to MCB. GWAS was done in a subset (270 inbreds) of the maize diversity panel using the maize 50 k SNP genotyping array [56] plus a more recently developed set of 425 k SNPs found through genotyping by sequencing [57,58].

Results

Means, analysis of variance and heritabilities

Differences among inbred lines were highly significant (P < 0.01) for all resistance (TL, SD, and KR) and agronomic traits (PH, DTA, and DTS); while significant (P < 0.05) genotype by environment (G × E) interactions were observed for all traits as shown in Additional file 1: Table S1 and S2. The inbred means for TL ranged from 5.2 to 49.2 cm with an overall mean of 20.9 cm, and for SD from 4.1 to 22.6% with an overall mean of 11.9%, respectively. The inbred mean for KR ranged from 5.4 to 9 with an overall mean of 7.8 in the subjective scale. Higher values for TL were observed in 2012 (overall mean = 43.3 cm) compared with 2010 (overall mean = 17.5 cm) and 2011 (overall mean = 14.7 cm, Additional file 1: Table S3 and S4). Intermediate values for heritability on a line mean-basis were estimated for TL (h 2 = 0.60) and KR (h 2 = 0.52) across the three years, whereas a low heritability value was obtained for SD (h 2 = 0.25, Additional file 1: Table S1).

Correlation analysis

A significant and high (r > 0.50) phenotypic correlation was observed between TL and SD. Genetic correlation coefficients were also significant and high between TL and SD and between TL and PH (Table 1). KR showed significant phenotypic correlations with other resistance traits but the correlation coefficients were not higher than 0.5. Genetic correlations between KR and the agronomic traits (DTA, DTS and PH) were significant and high (Table 1).

Table 1 Genotypic 1 (above diagonal) and phenotypic 2 (below diagonal) correlation coefficient estimates for each pair of traits

Association analysis of maize resistance to MCB and agronomic traits

The compressed mixed linear model computed for each trait in Tassel reduced the pairwise kinship matrix by clustering the 267 lines into 32 groups for TL, 201 groups for SD, 166 groups for KR, 192 for PH, and 240 for DTA, and DTS (Table 2). The proportion of the total phenotypic variation explained by background genetic effects was 39, 13, and 48% for TL, SD, and KR, respectively. By comparison, the background polygenic effects modeled by the K matrix accounted for 60, 94 and 90% of total variation for PH, DTA, and DTS, respectively.

Table 2 Summary of the compressed mixed linear model analysis for three traits related to resistance to MCB attack and three agronomic traits in an inbred association panel evaluated in three years

Ten SNPs were identified as significantly associated with length of tunnels made by MCB (Figure 1, Table 3). Based on the additive effects, the major allele reduce TL for all significant SNPs except for the SNP located on chromosome 10 (Table 3). Four SNPs were significantly associated with SD made by MCB (Figure 1). The minor allele for those SNPs increases the SD from 1.1 to 2%. The total variance explained (R 2) by each SNP associated with TL and SD ranged from 7 to 9%. Eleven SNPs were significantly associated with KR (Figure 1) and the total variance (R 2) explained by each of them ranged from 6 to 8%. The minor allele for those SNP reduced from 0.15 to 0.40 points the ratio that accounts for kernel resistance.

Figure 1
figure 1

GWAS results for the three resistance traits to MCB attack in a maize association panel. Each graph (TL, SD, and KR) represent the P-values of the 246,477 SNPs tested for each resistance trait. Each row indicates the SNP significantly associated (RMIP ≥ 0.30) to each resistance trait analyzed.

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 resistance traits (TL, SD, and KR), and significance values for the association between the SNP and the phenotype ( P -value and RMIP)

Nineteen SNPs were significantly associated with PH and 48 and 43 SNPs were significantly associated with DTA and DTS, respectively (Additional file 1: Table S5, S6, and S7). But none of those SNPs coincide or are close to those detected for resistance traits.

Candidate genes selection

The filtered predicted gene set from the annotated B73 reference maize genome [60] was used to characterize the gene containing or nearby the SNP declared significant. Seven candidate genes containing or adjacent to the SNPs significantly associated with TL, four candidate genes containing or adjacent to the significant SNPs associated to SD, and ten candidate genes containing or adjacent to the significant SNPs associated to KR were proposed (Table 4).

Table 4 Candidate genes for each significantly SNP associated with TL, SD, and KR and its respective encoding product

In general, the LD (r 2) between significantly associated SNPs and SNPs around them decays (r 2 ≤ 0.10) rapidly; but, in some cases, the LD spans as much as 0.5 Mb at an r 2 value greater than 0.20 (Figure 2). For each trait, LD estimates between significant SNPs were always below 0.1 except for SNPs significant for TL located on chromosome 7 as shown in Additional file 1: Table S8.

Figure 2
figure 2

Linkage disequilibrium heat chart showing LD measure ( r 2 ) between the SNP significantly associated with traits related to resistance to MCB attack and the closest 60 SNPs. Each bar represent a region (ranging from ~0.15 to ~1 Mbp) containing each significant SNP associated (black square) to resistance traits. The LD (value of r 2) with the 30 upstream SNPs were shown at the right side of the black square, and the LD with the 30 downstream SNPs were shown at the left side of the black square; on each bar, the extreme distances (in kbp) covered by the upstream and downstream SNPs are indicated.

Discussion

Means, heritabilities, and correlations

Substantial differences among yearly means for TL made by MCB were observed across environments, with the highest values of TL observed in 2012. Artificial infestation is used to ensure the contact of the insect with the plant, but the severity of the attack is conditioned by environmental factors that in turn influence natural infestation. The prevalence of the pest can vary greatly between and within locations [5,6]. Data on the monthly samplings at different locations in Pontevedra (unpublished data) indicate a higher rate of natural infestation in 2012 compared with 2010 and 2011 in a plot adjacent to the GWAS trial. In addition, the experiment was harvested earlier in 2011 compared with 2010 and 2012, limiting the time that larvae could damage the plants. Therefore, differences on infestation levels and harvest time could be major causes of the observed differences in TL among years.

Significant genotype × environment interaction was observed for all traits, although no G × E interaction for TL was found in previous studies under infestation with MCB [18,61]. However, G × E interaction for resistance traits was not significant, except for KR, when we made the analysis discarding data from 2012. Therefore, the high values for TL in 2012 could be obscuring the genotype effect. In previous studies, typically no G × E interaction for KR was found [61-63], with some exceptions [18].

The heritabilities for resistance traits ranged from low to moderate while PH and flowering time were higher as expected. The heritability for TL estimated herein is within the range (h 2 = 0.40-0.77) of those obtained in other works with biparental crosses under infestation with MCB or ECB [16,18,64].

A significant genetic relationship was observed between TL and PH (r g  = 0.51) in a panel of diverse origin. Nevertheless, we did not find a single SNP or linked SNPs associated with both traits, suggesting that it is primarily due to the polygenic background. This relationship has also been important in some previous experiments with MCB and ECB [16,32], but negligible in others [17,30].

A high genetic correlation coefficients between KR and flowering time (r g  = 0.75 for DTA and 0.71 for DTS) and between KR and PH were found, which suggest that late and taller genotypes will be the healthiest; but these results have to be taken with caution because infestation was made simultaneously for all genotypes placing the MCB eggs close to the ground and the higher distance between the eggs and the ear of the taller and later plants could have impeded the larval arrival to the ear. In addition, when stem tissue is more lignified at time of infestation (as would be the case for earlier flowering genotypes), the preference of MCB larvae for the stem tissue compared to ear tissue could be less evident.

Association analysis

There was a minimal variation in model fit of the compressed MLM among different traits because similar compression level values were observed for all traits (c = 1.1 to 1.6), except for TL, which had the highest value for compression level (c = 8.3). However, Zhang et al. [65] shown that this method controls the false positive rate well when the compression levels ranged from 1.5 to 10.

SNPs located on chromosomes 3 and 7 for TL co-localized with previously reported QTLs for TL by corn borers in genome bins 3.02 and 7.03 [30]. The proportion of the phenotypic variance explained (R 2 = 7 - 9%) by each SNP was comparable to that explained by QTLs for TL made by ECB and MCB (R 2 = 3.5 – 15.7%) in biparental crosses [16-18,64]. No QTLs for SD and KR have been previously reported in biparental crosses in the same regions where we located the significant SNPs, except one QTL for KR made by MCB at the bin 5.04 [66]. Therefore, association mapping uncovers additional genomic regions involved in maize resistance to corn borers that were not detected using biparental populations.

Candidate genes

We used the maize B73 genome v2 (RefGen_v2) available from the Maize GDB [67] (http://www.maizegdb.org/) to identify genes that either include or are close to the significantly associated SNPs. A region of approximately 0.2 Mb around the SNP was checked for annotated genes putatively involved in plant response to wounding and/or damage by microbes (insects or pathogens) based on bibliographic records.

Genes associated with TL made by MCB

The candidate gene adjacent to the significant SNP associated to TL on chromosome 2 encodes a Tetratricopeptide repeat (TPR) protein containing. The TPR is one of many repeat motifs that form structural domain mediating protein-protein interactions in several cellular process including translocation and degradations of proteins [68,69]. A recent study has proposed that the presence of those protein-protein interaction motifs could be acting as a modulator of the gene function and protein expression during the stress response caused by invading pathogens [70]. The hex1 gene containing the significant SNP on chromosome 3 encodes hexokinase, a sugar sensor with numerous physiological functions within the cell including response to oxidative-stress and pathogen resistance [71-73]. A gene located on chromosome 4 that encodes a Double Clp-N motif-containing P-loop nucleoside triphosphate hydrolase superfamily protein contains a significant SNP. This gene shows a weak similarity to the AtHSP101 gen in Arabidopsis, that codifies for a heat shock protein required for acclimation to high temperatures, and probably could be involved in response to other stresses [74]. Furthermore the ofp44 (OVATE-transcription factor 22) gene is relatively close to this SNP. It is known that some transcription factors from the Ovate family protein interact with other transcription factor families such as NAC domain protein1, MYB transcription factors, and KNOX homeodomain protein to regulate the synthesis of the three major components of secondary cell wall (lignin, cellulose and hemicellulose) in A. thaliana [75-78] and plants with a fortified cell wall are more resistance to corn borer attack [79].

The candidate gene for the other significant SNP on chromosome 4 encodes a Phospholipase A2 (PLA2) which plays a very important role in signal transduction in plants since it is the precursor of oxylipins and jasmonated acid, two hormones which regulates defense genes against herbivores [80-82].

A 1 Mbp region on chromosome 7 contains three SNPs that were significantly associated with TL and were in significant LD with at least one of the other associated SNPs in the region (r 2 > 0.2; Additional file 1: Table S8). The SNPs at 154,739,818 and 154,741,622 bp are both located within a gene that putatively encodes for a Histidine kinase, hybrid-type, ethylene sensor; five other genes encoding Serine/Threonine kinase receptor and receptor-like ser/thre kinases family proteins (RLK) were also physically nearby. It is well known that both kinases and RLK proteins play a central role in signaling during pathogen recognition and the subsequent activation of plant defense mechanisms [83-85]. They are also involved in wound-mediated defense response [86], and maintenance of plant cell wall integrity [87]. Polymorphisms at the proposed kinase genes could also be responsible for the QTLs at bin 7.03 detected for TL in a biparental population [30]. The candidate gene encoding the maize Calcium/calmodulin-dependent protein kinase (CDPK) close to the three significant SNPs on chromosome 10, could be playing an important role in the activation of defense against the attack of MCB since it has been known that the CDPK is induced by mechanical wounding by herbivore attack inducing accumulation of jasmonic acid in maize and other species [88-91].

Genes associated with SD made by MCB

A gene encoding an Actin polymerizing factor 4 (APF4) contains the second SNP at chromosome 1 significantly associated to SD. One of the functions of APF4 protein is remodeling the actin of cytoskeleton under different stimulus, including wounding and pathogen attacks [92]. Some studies in Arabidopsis indicated that the APF4 mediated defense signal and it is also relates with actin dynamic of cytoskeleton during the innate immune signaling [93-95]. The candidate gene for the SNP significantly associated to SD on chromosome 2 encodes a Leucine-rich receptor-like protein kinase (LRR-RLK). LRR-RLK protein family plays an important role in cell-cell signaling and other signals involving peptide in ligands. They are involved in systemic activation of protease inhibitors in response to wounding by insect feeding [96,97]. On the other hand, the significant SNP located on chromosome 5 is within a gene encoding a protein with unknown function, and it is interesting that the SNP is close to the gene nactf30 which encodes a NAC domain protein transcription factor. As already mentioned, this gene and other transcription factor family members regulate the synthesis of secondary cell wall [98,99]. In addition, other authors have associated the NAC domain proteins with response to stresses made by herbivore attack [100].

Genes associated with KR made by MCB

The SNP on chromosome 3 significantly associated to KR was located nearby the cystatin1 gene, whose product is the corn kernel cysteine proteinase inhibitor//cysteine proteinase inhibitor I (psei1), an anti-metabolic protein synthetized and stored in the maize kernel [101]. The expression of some proteinase inhibitor genes are induced in response to mechanical wounding and insect damage [102], and it has been shown that cysteine proteinase inhibitor interferes with the digestive process of insects [103,104]. LD is low in this region, but it is interesting that the SNP significantly associated with KR was in LD with a SNP positioned in the exonic region of the gene (Figure 3). Therefore, if the association found between the SNP at position 187,742,562 and KR is due to the linkage between that SNP and a certain undetected polymorphism in the cysteine1 gene, the effect of this polymorphism would be expected to be especially large because the linkage between the significant SNP and SNPs in the cystatin1 gene is low, although significant.

Figure 3
figure 3

A region of approximately 280 kpb in chromosome 3 where a SNP significantly associated with borer kernel resistance (KR) has been found. The red points represent the P-values of the SNPs in the region mentioned above. The orange strip indicates the location of the SNP significantly associated with KR and the dark blue strips indicate genomic locations in LD (r 2 > 0.2) with the SNP significantly associated with KR. The solid black line and the green rectangle indicate the intronic and exonic region of the cystatin1 gene, respectively. The SNP significantly associated with KR herein is in LD with the exonic region of the cystatin1 gene, which encodes a maize cysteine proteinase inhibitor involved in plant resistance to insects.

Another significant SNP associated to KR was located close to a gene encoding a basic Helix-Loop-Helix (bHLH) transcription factor family member. Recent studies in other plant species have demonstrated that this family gene protein has an important role in regulating jasmonic acid response [105-109]. The third significant SNP associated with KR located on chromosome 3 was within a gene encoding a starch branching enzyme interacting protein, whose deficiency leads to decreased digestibility of maize kernel [110]. No previous evidence about its involvement in resistance to insect was found, however. A gene that encodes a Germin 1–2 protein is 106 kb downstream from this associated SNP; these proteins are implicated in the response to abiotic and biotic stresses including the response to mechanical wounding and insect damage [111-114]. Another significant SNP located on chromosome 3 at position 222,733,400 is within a gene that encodes a Glycine rich protein (GRP), a structural protein which could be playing an important role in the fortification of plant cell wall [115-117], also this protein is wound-inducible [118,119]. A set of plant Thionin family protein precursor genes were found nearby the SNP on chromosome 3 at position 222,733,400. This finding is particularly interesting since these proteins belong to the Defensin family protein which has antibacterial, antifungal, and insecticide activities [120-124].

The two significant SNPs associated to KR on chromosome 4 are adjacent to a gene that encodes for an Ubiquitin thiolesterase. These proteins form complex systems of selective protein degradation [125,126] and mediate the biosynthesis of plant hormone signaling such as salicylic, jasmonic, and abscisic acid and auxins [127,128].

The candidate gene containing the significant SNP associated to KR on chromosome 5 encodes a Phytosulfokine receptor (PSK), a recognition hormone whose level of expression has been increased by pathogens elicitors [129,130].

The candidate gene, adjacent to the significant SNP located on chromosome 6, encodes a Polyamine oxidase (propa-1,3-aimine-forming) (PAO). PAO plays an important role in stress tolerance by generating H2O2 which is a key component in signal transduction pathways leading to stress responses. In Zea mays, the function of PAO in wound-healing is likely due to increased lignin and suberine deposition as consequence of H2O2 release [131]. In addition, it has been described to play a role in cell wall stiffening and mediating abiotic and biotic stresses [132]. The candidate adjacent gene to the significant SNP associated to KR on chromosome 7 encodes a Catalase//L-ascorbate peroxidase, they are two major hydrogen peroxide-detoxifying enzymes whose activity is very important in the reduction of the oxidative stress caused by H2O2 [133,134], and the importance of these detoxifying enzymes in resistance to insect attack has been recently reported [135].

A candidate gene for the significant SNP on chromosome 7 at position 19,347,596 is adjacent to the SNP that encodes a Kinase associated protein phosphatase. This gene family has been proposed to regulate the response to different type of stresses including pathogens and herbivore attack [136]. It is interesting to note that another close gene to this SNP is the ipt4 (isopentenyl transferase4) gene, which is involved in the regulation of cytokinin biosynthesis pathway [137]. The expression of an itp gen (fused with a wound-inducible promoter) in transgenic plants of Nicotiana plumbaginifolia decreased leaf consumption by Manduca sexta (lepidopterous) and reduced survival of Myzus persicae (aphid) [138]. Although there are several reports about the importance of cytokinins in the modulation of plant defense against pathogens and insect attack, their role is not clear [139].

Unlike markers linked to QTLs for resistance to corn borers, the SNPs associated with resistance in the present study could be incorporated in whole-genome predictor models in order to improve genomic selection [140]. Marker–assisted selection has proved an useful tool for improving resistance to the European corn borer [141] of segregating bi-parental populations related to the mapping populations used for QTL detection, but could be inappropriate in non-structured populations. As additive effects are the most important genetic effects for resistance traits, crossing inbreds with improved resistance will render hybrids more resistant to attack by MCB larvae.

Conclusion

We conducted a genome wide association study for resistance traits to MCB with more than 245 kb SNP distributed through the whole genome. We found a set of significant SNPs associated to the three resistance traits to MCB attack. Of which 10 SNPs were significant associated to TL, 4 SNPs were associated to SD, and 11 SNPs were associated to KR. In general, each of these SNPs explain a considerable proportion of the phenotypic variance (R 2 = 6-9%). No co-localized SNPs were found for resistance and agronomic traits that could underlie the genetic correlations found between these traits.

Twenty one candidate genes were proposed for the three resistant traits, they were either containing or adjacent (within a window of ±130 kbp) to each significantly associated SNP.

Most of the candidate genes proposed herein are part of the signaling pathway, others act as regulator of expression under biotic stress condition, and a few genes are encoding enzymes with antibiotic effect against insects such as the cystatin1 gen and the defensin proteins.

The identification of these polymorphisms associated to resistance traits to MCB attack can be useful to understand the molecular mechanisms that affect resistance and susceptibility of host plants to insect attack, in order to contribute to advance in the understanding of plant-insect interactions. Nevertheless further studies are necessary to validate the candidate genes identified herein.

Methods

Plant material and phenotypic data

The maize diversity panel (composed of 302 inbred lines) represents much of the diversity available in public breeding sector around the world. A subset of the maize diversity panel (henceforth we will refer to this population as “association panel”) was evaluated for resistance to MCB attack at Pontevedra (42°24’ N, 8°38’ W, and 20 m above sea level), Spain, through three years (2010, 2011, and 2012). A subset of 270 inbred lines was assayed in an 18 × 15 α-lattice design with two replications in 2010 and 2011. In the third year a subset of 255 inbred lines was assayed (because we did not have enough seed for the remaining15 lines) in a 17 × 15 α-lattice design with two replications.

The trials were hand-planted and each experimental plot consisted of one row spaced 0.8 m apart from the other row with 29 two-kernel hills spaced 0.18 m apart. Plots were overplanted and thinned, obtaining a final density of ~70,000 plant ha−1. The evaluations were performed under artificial infestation with eggs of MCB. The eggs for inoculation were obtained at the Misión Biológica de Galicia by rearing the insect as described by Eizaguirre and Albajes [142]. Five plants of each plot were infested with ~ 40 MCB eggs placed between the stem and the sheath of a basal leaf.

Data collected were: tunnel length (TL), the mean length of stem tunnels made by borers on the five infested plants, which were longitudinally split at the time of harvest; stem damage (SD) as the percentage of the stem damaged by MCB larvae on the five infested plants; kernel resistance to borer attack (KR) recorded at harvest as the damage on the main ear of the five infested plants according to a subjective visual resistance scale of 1 to 9 in which 1 indicates completely damaged and 9 indicates no damage; days to anthesis (DTA) and days to silking (DTS) as the days from planting to the date on which 50% of plants were shedding pollen or showing silks, respectively; and plant height (PH) on five representative plants as the distance from the ground to the top of the plant.

Genotypic data

We used a set of unique SNP markers derived from a Illumina maize 50 k array [56] and a genotyping-by-sequencing (GBS) strategy [48]. The two data sets were combined and filtered to exclude SNPs with more than 20% missing genotype data and minor allele frequency (MAF) less than 5%. Heterozygous genotypes were considered as missing data in the analysis. After filtering, a total of 246,477 SNPs (Additional file 2) distributed across the maize genome were used in this study.

A genetic kinship matrix (K) previously published by Olukolu et al. [52] was used for GWAS. The kinship matrix was estimated using a subset of 5000 SNPs without any missing genotypes and distributed approximately uniformly across the entire genome (at least 60 kbp between any two markers).

Statistical analyses

Best linear unbiased estimator (BLUE)

Each trial was analyzed separately with the SAS mixed model procedure (PROC MIXED) in SAS software version 9.3 [143] considering inbred lines as a fixed effect and replications and block within replication as random effects. Then, trials were combined using a mixed linear model across years and considering inbred lines as the only fixed effect. As large predicted values for stem damage and tunnel length were associated with larger residuals, a natural logarithmic transformation of TL and SD scores was used for obtaining BLUEs. The logarithmic transformation eliminated the relationship between residual variance and predicted values. Line BLUEs were back-transformed and then used to perform GWAS.

Heritabilities

Heritabilities (ĥ 2) for each year were estimated for each trait on a family-mean basis as described previously by Holland et al. [144]. The model for these analyses was similar to the model mentioned above with the exception that inbred lines were treated as random effects. The genetic (r g ) and phenotypic (r p ) correlations between traits were computed using REML estimation in SAS mixed model procedure following Holland [145].

Association analysis

Genome-wide association analysis based on mixed linear model (MLM) was performed in Tassel 4.1.26 [146]. The MLM used by Tassel was

$$ \mathrm{y}=\mathrm{X}\beta +\mathrm{Z}\mathrm{u}+\mathrm{e} $$

where y is the vector of phenotypes (BLUEs), β is a vector of fixed effects, including the SNP marker tested, u is a vector of random additive effects (inbred lines), X and Z represents matrices, and e is a vector of random residuals. The variance of random line effects was modeled as \( \mathrm{V}\mathrm{a}\mathrm{r}\left(\mathrm{u}\right)=\mathbf{K}{\sigma}_a^2 \), where K is the n × n matrix of pairwise kinship coefficient and \( {\sigma}_a^2 \) is the estimated additive genetic variance [147].

Restricted maximum likelihood estimates of variance components were obtained by using the optimum compression level (compressed MLM) and population parameters previously determined options (P3D) in Tassel [65]. The optimum compression level option reduces the computation demand by clustering the (n) total individuals into (s) groups based on their realized genomic relationships, allowing the original K matrix to be replaced by a smaller relationship matrix. The P3D option uses iteration to estimate population parameters such as genetic and residual variance only once in a model with no fixed marker effects, then uses those estimates without iteration in subsequent association tests for each marker. The combination of these two methods reduces computational time and improves model fit [65].

Threshold for GWAS

To identify SNPs with the most robust associations with traits, a subsampling or subagging procedure was employed in GWAS analysis [148,149]. Each of 100 subsampled datasets generated using the R software [150] comprised a random sample of 80% of inbred lines from the diversity population. Only SNP markers determined as significant at P < 1 × 10−4 and subsequently detected in ≥ 30 subsamples, i.e. resample model inclusion probability (RMIP) threshold of 0.30, were considered as significantly associated to the trait under study.

Linkage disequilibrium and candidate gene selection

We examined the linkage disequilibrium (LD) measure (r 2) with each SNP significantly associated with resistance traits in a region of 60 SNPs (ranging from 0.15 to 1 Mbp). For each trait, the linkage disequilibrium between significant SNPs was also calculated. The genes containing or adjacent to SNPs significantly associated with traits were identified and characterized by the use of MaizeGDB genome browser [67]. We examined a region around each significant SNP in order to identify candidate genes of interest. For most SNPs, more than two genes that could be involved in plant defense mechanism (based on previously published evidence) were preselected; the closest gene to the SNPs significantly associated with resistance was selected as candidate gene and the remaining genes were presented as other interesting genes that could be putatively involved in plant defense against insect herbivores.

Abbreviations

APF4:

Actin polymerizing factor 4

bHLH:

Basic Helix-Loop-Helix

BLUE:

Best linear unbiased estimator

Bt:

Bacillus thuringiensis

CDPK:

Calcium/calmodulin-dependent protein kinase

DTA:

Days to anthesis

DTS:

Days to silking

ECB:

European corn borer

GBS:

Genotyping by sequencing

GRP:

Glycine rich protein

GWAS:

Genome-wide association study

KR:

Kernel resistance

LD:

Linkage disequilibrium

LRR-RLK:

Leucine-rich receptor-like protein kinase

MAF:

Minor allele frequency

MCM:

Mediterranean corn borer

PAO:

Polyamine oxidase

PH:

Plant height

PLA2:

Phospholipase A2

PSK:

Phytosulfokine receptor

QTL:

Quantitative trait locus/loci

REML:

Restricted maximum likelihood

RLK:

Receptor like kinase

RMIP:

Resample model inclusion probability

SD:

Stem damage

SNP:

Single nucleotide polymorphism

TL:

Tunnel length

TPR:

Tetratricopeptide repeat

References

  1. Meihls LN, Kaur H, Jander G. Natural variation in maize defense against insect herbivores. Cold Spring Harb Sym. 2012;77:269–83.

    CAS  Google Scholar 

  2. Malvar RA, Cartea ME, Revilla P, Ordás A, Alvarez A, Mansilla JP. Sources of resistance to pink stem borer and European corn borer in maize. Maydica. 1993;38:313–9.

    Google Scholar 

  3. Avantaggiato G, Quaranta F, Desiderio E, Visconti A. Fumonisin contamination of maize hybrids visibly damaged by Sesamia. J Sci Food Agric. 2002;83:13–8.

    Google Scholar 

  4. Visconti A, Marasas WFO, Miller JD, Riley R. Mycotoxins of growing interest: Fumonisins. In: Book Mycotoxins of growing interest: Fumonisins. 1999. p. 1–22.

    Google Scholar 

  5. Cordero A, Malvar RA, Butrón A, Revilla P, Velasco P, Ordás A. Population dynamics and life-cycle of corn borers in South Atlantic European Coast. Maydica. 1998;43:5–12.

    Google Scholar 

  6. Velasco P, Revilla P, Monetti L, Butrón A, Ordás A, Malvar RA. Corn borers (Lepidoptera: Noctuidae; Crambidae) northwestern Spain: population dinamics and distribution. Maydica. 2007;52:195–203.

    Google Scholar 

  7. Meissle M, Romeis J, Bigler F. Bt maize and integrated pest management - a European perspective. Pest Manag Sci. 2011;67:1049–58.

    CAS  PubMed  Google Scholar 

  8. Speiser B, Tamm L. Regulation of plant protection in organic farming. In: Ehlers R-U, editor. Regulation of Biological Control Agents. Kiel, Germany: Springer; 2011. p. 113–25.

    Google Scholar 

  9. Campagne P, Kruger M, Pasquet R, Le Ru B, Van den Berg J. Dominant inheritance of field-evolved resistance to Bt corn in Busseola fusca. PLoS One. 2013;8:e69675.

    PubMed Central  CAS  PubMed  Google Scholar 

  10. Tabashnik BE, Brévault T, Carrière Y. Insect resistance to Bt crops: lessons from the first billion acres. Nat Biotechnol. 2013;31:510–21.

    CAS  PubMed  Google Scholar 

  11. González-Cabrera J, García M, Hernández-Crespo P, Farinós GP, Ortego F, Castañera P. Resistance to Bt maize in Mythimna unipuncta (Lepidoptera: Noctuidae) is mediated by alteration in Cry1Ab protein activation. Insect Biochem Mol Biol. 2013;43:635–43.

    PubMed  Google Scholar 

  12. Gould F. Sustainability of transgenic insecticidal cultivars: integrating pest genetics and ecology. Annu Rev Entomol. 1998;43:701–26.

    CAS  PubMed  Google Scholar 

  13. Butrón A, Malvar R, Velasco P, Vales M, Ordás A. Combining abilities for maize stem antibiosis, yield loss, and yield under infestation and non infestation with pink stem borer. Crop Sci. 1999;39:691–6.

    Google Scholar 

  14. Cartea ME, Malvar RA, Butrón A, Vales MI, Ordás A. Inheritance of antibiosis to Sesamia nonagrioides (Lepidoptera: Noctuidae) in maize. J Econ Entomol. 1999;92:994–8.

    Google Scholar 

  15. 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.

    CAS  PubMed  Google Scholar 

  16. Schön CC, Lee M, Melchinger AE, Guthrie WD, Woodman WL. Mapping and characterization of quantitative trait loci affecting resistancie against second-generation European corn borer in maize with the aid of RFLPs. Heredity. 1993;70:648–59.

    Google Scholar 

  17. Cardinal AJ, Lee M, Sharopova N, Woodman WL, Long MJ. Genetic mapping and analysis of quantitative trait loci for resistance to stalk tunneling by European corn borer in maize. Crop Sci. 2001;41:835–45.

    CAS  Google Scholar 

  18. Ordás B, Malvar RA, Santiago R, Sandoya G, Romay MC, Butrón A. Mapping of QTL for resistance to the Mediterranean corn borer attack using the intermated B73 × Mo17 (IBM) population in maize. Theor Appl Genet. 2009;119:1451–9.

    PubMed  Google Scholar 

  19. Klenke J, Russel W, Guthrie W. Recurrent selection for resistance to European corn borer in a corn synthetic and correlated effects on agronomic traits. Crop Sci. 1986;26:864–8.

    Google Scholar 

  20. Bosque-Pérez N, Kling J, Odubiyi S. Recent advances in the development of sources of resistance to pink stalk borer and African sugarcane borer. In: Insect resistant maize: Recent advances and utilization Proc of an International Symposium. 1997. p. 234–40.

    Google Scholar 

  21. Sandoya G, Butrón A, Alvarez A, Ordás A, Malvar RA. Direct response of a maize synthetic to recurrent selection for resistance to stem borers. Crop Sci. 2008;48:113–8.

    Google Scholar 

  22. Butrón A, Romay MC, Peña-Asin J, Alvarez A, Malvar RA. Genetic relationship between maize resistance to corn borer attack and yield. Crop Sci. 2012;52:1176–80.

    Google Scholar 

  23. Ordás B, Butrón A, Alvarez A, Revilla P, Malvar R. Comparison of two methods of reciprocal recurrent selection in maize (Zea mays L.). Theor Appl Genet. 2012;124:1183–91.

    PubMed  Google Scholar 

  24. Malvar RA, Butrón A, Revilla P, Ordás A. Resistance to the pink stem borer, Sesamia nonagrioides, in maize. Recent Res Devel Plant Sci. 2004;2:1–32.

    Google Scholar 

  25. Cartea ME, Malvar RA, Vales I, Butrón A, Ordás A. Inheritance of resistance to ear damage caused by Sesamia nonagrioides (Lepidoptera: Noctuidae) in maize. J Econ Entomol. 2001;94:277–83.

    CAS  PubMed  Google Scholar 

  26. Hudon M, Chiang MS. Evaluation of resistance of maize germplasm to the univoltine European corn borer Ostrinia nubilalis (Hübner) and relationship with maize maturity in Quebec. Maydica. 1991;36:69–74.

    Google Scholar 

  27. Guthrie W, Russell W. Breeding methodologies and genetic basis of resistance in maize to the European corn borer. In: International Symposium on Methodologies for Developing Host Plant Resistance to Maize Insects Mexico, DF (Mexico) 9–14 Mar 1987. 1989.

    Google Scholar 

  28. Krakowsky MD, Brinkman MJ, Woodman WL, Lee M. Genetic components of resistance to stalk tunneling by the European corn borer in maize. Crop Sci. 2002;42:1309–15.

    Google Scholar 

  29. Krakowsky MD, Lee M, Holland JB. Genotypic correlation and multivariate QTL analyses for cell wall components and resistance to stalk tunneling by the European corn borer in maize. Crop Sci. 2007;47:485–8.

    Google Scholar 

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

    CAS  Google Scholar 

  31. Orsini E, Krchov L, Uphaus J, Melchinger A. Mapping of QTL for resistance to first and second generation of European corn borer using an integrated SNP and SSR linkage map. Euphytica. 2012;183:197–206.

    CAS  Google Scholar 

  32. Ordás B, Malvar RA, Santiago R, Butrón A. QTL mapping for Mediterranean corn borer resistance in European flint germplasm using recombinant inbred lines. BMC Genomics. 2010;11:1–10.

    Google Scholar 

  33. Bohn M, Khairallah MM, González-de-León D, Hoisington DA, Utz HF, Deutsch JA, et al. QTL mapping in tropical maize: I. Genomic regions affecting leaf feeding resistance to sugarcane borer and other traits. Crop Sci. 1996;36:1352–61.

    Google Scholar 

  34. Groh S, González-de-León D, Khairallah MM, Jiang C, Bergvinson D, Bohn M, et al. QTL mapping in tropical maize III. Genomic regions for resistance to Diatraea spp and associated traits in two RIL populations. Crop Sci. 1998;1998:1062–72.

    Google Scholar 

  35. Bohn M, Khairallah M, Jiang C, González-de-León D, Hoisington D, Utz H, et al. QTL mapping in tropical maize: II. Comparison of genomic regions for resistance to Diatraea spp. Crop Sci. 1997;37:1892–902.

    CAS  Google Scholar 

  36. Weiss LA, Arking DE, Daly MJ, Chakravarti A, Brune CW, West K, et al. A genome-wide linkage and association scan reveals novel loci for autism. Nature. 2009;461:802–8.

    PubMed Central  CAS  PubMed  Google Scholar 

  37. Duerr RH, Taylor KD, Brant SR, Rioux JD, Silverberg MS, Daly MJ, et al. A genome-wide association study identifies IL23R as an inflammatory bowel disease gene. Science. 2006;314:1461–3.

    CAS  PubMed  Google Scholar 

  38. Sladek R, Rocheleau G, Rung J, Dina C, Shen L, Serre D, et al. A genome-wide association study identifies novel risk loci for type 2 diabetes. Nature. 2007;445:881–5.

    CAS  PubMed  Google Scholar 

  39. Barendse W, Reverter A, Bunch RJ, Harrison BE, Barris W, Thomas MB. A validated whole-genome association study of efficient food conversion in cattle. Genetics. 2007;176:1893–905.

    PubMed Central  CAS  PubMed  Google Scholar 

  40. Kijas JW, Townley D, Dalrymple BP, Heaton MP, Maddox JF, McGrath A, et al. A genome wide survey of SNP variation reveals the genetic structure of sheep breeds. PLoS One. 2009;4:e4668.

    PubMed Central  PubMed  Google Scholar 

  41. Bolormaa S, Hayes B, Savin K, Hawken R, Barendse W, Arthur P, et al. Genome-wide association studies for feedlot and growth traits in cattle. J Anim Sci. 2011;89:1684–97.

    CAS  PubMed  Google Scholar 

  42. Fan B, Onteru SK, Du Z-Q, Garrick DJ, Stalder KJ, Rothschild MF. Genome-wide association study identifies loci for body composition and structural soundness traits in pigs. PLoS One. 2011;6:e14726.

    PubMed Central  CAS  PubMed  Google Scholar 

  43. Tian F, Bradbury PJ, Brown PJ, Hung H, Sun Q, Flint-Garcia S, et al. Genome-wide association study of leaf architecture in the maize nested association mapping population. Nat Genet. 2011;43:159–62.

    CAS  PubMed  Google Scholar 

  44. Neumann K, Kobiljski B, Denčić S, Varshney R, Börner A. Genome-wide association mapping: a case study in bread wheat (Triticum aestivum L.). Mol Breed. 2011;27:37–58.

    Google Scholar 

  45. Zhu C, Gore M, Buckler ES, Yu J. Status and prospects of association mapping in plants. Plant Genome. 2008;1:5–20.

    CAS  Google Scholar 

  46. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  47. Oraguzie NC, Rikkerink EH, Gardiner SE, Silva H. Association mapping in plants. New York: Springer-Verlag GmbH; 2007.

    Google Scholar 

  48. Romay MC, Millard MJ, Glaubitz JC, Peiffer J, Swarts K, Casstevens TM, et al. Comprehensive genotyping of the USA national maize inbred seed bank. Genome Biol. 2013;14:1–18.

    Google Scholar 

  49. Flint-Garcia SA, Thuillet A-C, 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:1054–64.

    CAS  PubMed  Google Scholar 

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

    PubMed Central  CAS  PubMed  Google Scholar 

  51. Cook JP, McMullen MD, Holland JB, Tian F, Bradbury P, Ross-Ibarra J, et al. Genetic architecture of maize kernel composition in the nested association mapping and inbred association panels. Plant Physiol. 2012;158:824–34.

    PubMed Central  CAS  PubMed  Google Scholar 

  52. Olukolu BA, Negeri A, Dhawan R, Venkata B, 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:609–20.

    PubMed Central  CAS  PubMed  Google Scholar 

  53. Zila CT, Samayoa LF, Santiago R, Butrón A, Holland JB. A genome-wide association study reveals genes associated with fusarium ear rot resistance in a maize core diversity panel. G3-Genes Genomes Genet. 2013;3:2095–104.

    Google Scholar 

  54. Atwell S, Huang YS, Vilhjálmsson BJ, Willems G, Horton M, Li Y, et al. Genome-wide association study of 107 phenotypes in Arabidopsis thaliana inbred lines. Nature. 2010;465:627–31.

    PubMed Central  CAS  PubMed  Google Scholar 

  55. Chan EK, Rowe HC, Kliebenstein DJ. Understanding the evolution of defense metabolites in Arabidopsis thaliana using genome-wide association mapping. Genetics. 2010;185:991–1007.

    PubMed Central  CAS  PubMed  Google Scholar 

  56. 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:e28334.

    PubMed Central  CAS  PubMed  Google Scholar 

  57. Elshire RJ, Glaubitz JC, Sun Q, Poland JA, Kawamoto K, Buckler ES, et al. A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS One. 2011;6:e19379.

    PubMed Central  CAS  PubMed  Google Scholar 

  58. Glaubitz JC, Casstevens TM, Lu F, Harriman J, Elshire RJ, Sun Q, et al. TASSEL-GBS: A High Capacity Genotyping by Sequencing Analysis Pipeline. PLoS One. 2014;9:e90346.

    PubMed Central  PubMed  Google Scholar 

  59. Steel RG, Torrie JH. Principles and procedures of statistics. New York, Toronto, London: McGraw-Hill Book Company; 1960.

    Google Scholar 

  60. Schnable PS, Ware D, Fulton RS, Stein JC, Wei F, Pasternak S, et al. The B73 maize genome: complexity, diversity, and dynamics. Science. 2009;326:1112–5.

    CAS  PubMed  Google Scholar 

  61. Butrón A, Malvar RA, Velasco P, Cartea ME, Ordás A. Combining abilities and reciprocal effects for maize ear resistance to pink stem borer. Maydica. 1998;43:117–22.

    Google Scholar 

  62. Velasco P, Revilla P, Butrón A, Ordás B, Ordás A, Malvar RA. Ear damage of sweet corn inbreds and their hybrids under multiple corn borer infestation. Crop Sci. 2002;42:724–9.

    Google Scholar 

  63. Velasco P, Malvar RA, Butrón A, Revilla P, Ordás A. Ear feeding resistance of sweet corn inbreds to pink stem borer. J Am Soc Hortic Sci. 1999;124:268–72.

    Google Scholar 

  64. Bohn M, Schulz B, Kreps R, Klein D, Melchinger AE. QTL mapping for resistance against the European corn borer (Ostrinia nubilalis H.) in early maturing European dent germplasm. Theor Appl Genet. 2000;101:907–17.

    CAS  Google Scholar 

  65. Zhang Z, Ersoz E, Lai C-Q, Todhunter RJ, Tiwari HK, Gore MA, et al. Mixed linear model approach adapted for genome-wide association studies. Nat Genet. 2010;42:355–60.

    PubMed Central  CAS  PubMed  Google Scholar 

  66. Samayoa LF, Butrón A, Malvar RA. QTL mapping for maize resistance and yield under infestation with Sesamia nonagrioides. Mol Breeding. 2014;34:1331–44.

    Google Scholar 

  67. Harper LC, Schaeffer ML, Thistle J, Gardiner J, Andorf C, Campbell D, et al. The MaizeGDB Genome Browser tutorial: one example of database outreach to biologists via video. Database. 2011;2011:1–7.

    Google Scholar 

  68. Allan RK, Ratajczak T. Versatile TPR domains accommodate different modes of target protein recognition and function. Cell Stress Chaperones. 2011;16:353–67.

    PubMed Central  CAS  PubMed  Google Scholar 

  69. Blatch GL, Lässle M. The tetratricopeptide repeat: a structural motif mediating protein-protein interactions. Bioessays. 1999;21:932–9.

    CAS  PubMed  Google Scholar 

  70. Pogorelko GV, Mokryakova M, Fursova OV, Abdeeva I, Piruzian ES, Bruskin SA. Characterization of three Arabidopsis thaliana immunophilin genes involved in the plant defense response against P. syringae. Gene. 2014;538:12–22.

    CAS  PubMed  Google Scholar 

  71. Hellmann HA, Smeekens S. Sugar sensing and signaling in plants. Front Plant Sci. 2014;5:113.

    PubMed Central  PubMed  Google Scholar 

  72. Granot D, David-Schwartz R, Kelly G. Hexose kinases and their role in sugar-sensing and plant development. Front Plant Sci. 2013;4:1–17.

    Google Scholar 

  73. Morkunas I, Ratajczak L. The role of sugar signaling in plant defense responses against fungal pathogens. Acta Physiol Plant. 2014;8:409–14.

    Google Scholar 

  74. Lee U, Rioflorido I, Hong SW, Larkindale J, Waters ER, Vierling E. The Arabidopsis ClpB/Hsp100 family of proteins: chaperones for stress and chloroplast development. Plant J. 2007;49:115–27.

    CAS  PubMed  Google Scholar 

  75. Zhong R, Richardson EA, Ye Z-H. The MYB46 transcription factor is a direct target of SND1 and regulates secondary wall biosynthesis in Arabidopsis. Plant Cell. 2007;19:2776–92.

    PubMed Central  CAS  PubMed  Google Scholar 

  76. Zhong R, Lee C, Zhou J, McCarthy RL, Ye Z-H. A battery of transcription factors involved in the regulation of secondary cell wall biosynthesis in Arabidopsis. Plant Cell. 2008;20:2763–82.

    PubMed Central  CAS  PubMed  Google Scholar 

  77. Liu Y. Investigation of a KNAT7-BLH-OFP transcription factor complex involved in regulation of secondary cell wall biosynthesis in Arabidopsis thaliana. In: MSc Thesis. Vancouver, Canada: University of British Columbia, Botany; 2010.

    Google Scholar 

  78. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  79. Barros-Rios J, Malvar RA, Jung H-JG, Santiago R. Cell wall composition as a maize defense mechanism against corn borers. Phytochemistry. 2011;72:365–71.

    CAS  PubMed  Google Scholar 

  80. Ryan CA. The systemin signaling pathway: differential activation of plant defensive genes. BBA-Protein Struct M. 2000;1477:112–21.

    CAS  Google Scholar 

  81. Meijer HJ, Munnik T. Phospholipid-based signaling in plants. Annu Rev Plant Biol. 2003;54:265–306.

    CAS  PubMed  Google Scholar 

  82. Zhu-Salman K, Bi J-L, Liu T-X. Molecular strategies of plant defense and insect counter-defense. Insect Sci. 2005;12:3–15.

    Google Scholar 

  83. Afzal AJ, Wood AJ, Lightfoot DA. Plant receptor-like serine threonine kinases: Roles in signaling and plant defense. Mol Plant-Microbe Interact. 2008;21:507–17.

    CAS  PubMed  Google Scholar 

  84. Nürnberger T, Kemmerling B. Receptor protein kinases–pattern recognition receptors in plant immunity. Trends Plant Sci. 2006;11:519–22.

    PubMed  Google Scholar 

  85. Duan Y, Ge C, Liu S, Wang J, Zhou M. A two-component histidine kinase Shk1 controls stress response, sclerotial formation and fungicide resistance in Sclerotinia sclerotiorum. Mol Plant Pathol. 2013;14:708–18.

    CAS  PubMed  Google Scholar 

  86. G-i A, Kost C, Boland W. Herbivore-induced, indirect plant defences. BBA-Mol Cell Biol L. 2005;1734:91–111.

    Google Scholar 

  87. Engelsdorf T, Hamann T. An update on receptor-like kinase involvement in the maintenance of plant cell wall integrity. Ann Bot. 2014;114:1339–47.

    PubMed  Google Scholar 

  88. Szczegielniak J, Borkiewicz L, Szurmak B, Lewandowska-Gnatowska E, Statkiewicz M, Klimecka M, et al. Maize calcium-dependent protein kinase (ZmCPK11): local and systemic response to wounding, regulation by touch and components of jasmonate signaling. Physiol Plant. 2012;146:1–14.

    CAS  PubMed  Google Scholar 

  89. Yang D-H, Hettenhausen C, Baldwin IT, Wu J. Silencing Nicotiana attenuata calcium-dependent protein kinases, CDPK4 and CDPK5, strongly up-regulates wound-and herbivory-induced jasmonic acid accumulations. Plant Physiol. 2012;159:1591–607.

    PubMed Central  CAS  PubMed  Google Scholar 

  90. Szczegielniak J, Klimecka M, Liwosz A, Ciesielski A, Kaczanowski S, Dobrowolska G, et al. A wound-responsive and phospholipid-regulated maize calcium-dependent protein kinase. Plant Physiol. 2005;139:1970–83.

    PubMed Central  CAS  PubMed  Google Scholar 

  91. Ma F, Lu R, Liu H, Shi B, Zhang J, Tan M, et al. Nitric oxide-activated calcium/calmodulin-dependent protein kinase regulates the abscisic acid-induced antioxidant defence in maize. J Exp Bot. 2012;63:4835–47.

    PubMed Central  CAS  PubMed  Google Scholar 

  92. Staiger CJ, Gibbon BC, Kovar DR, Zonia LE. Profilin and actin-depolymerizing factor: modulators of actin organization in plants. Trends Plant Sci. 1997;2:275–81.

    Google Scholar 

  93. Tian M, Chaudhry F, Ruzicka DR, Meagher RB, Staiger CJ, Day B. Arabidopsis actin-depolymerizing factor AtADF4 mediates defense signal transduction triggered by the Pseudomonas syringae effector AvrPphB. Plant Physiol. 2009;150:815–24.

    PubMed Central  CAS  PubMed  Google Scholar 

  94. Henty-Ridilla JL, Li J, Day B, Staiger CJ. Actin depolymerizing factor4 regulates actin dynamics during innate immune signaling in Arabidopsis. Plant Cell. 2014;26:340–52.

    PubMed Central  CAS  PubMed  Google Scholar 

  95. Porter K, Shimono M, Tian M, Day B. Arabidopsis Actin-Depolymerizing Factor-4 links pathogen perception, defense activation and transcription to cytoskeletal dynamics. PLoS Pathog. 2012;8:e1003006.

    PubMed Central  CAS  PubMed  Google Scholar 

  96. Wu Y, Zhou JM. Receptor-Like Kinases in Plant Innate Immunity. J Integr Plant Biol. 2013;55:1271–86.

    CAS  PubMed  Google Scholar 

  97. Jones DA, Jones J. The role of leucine-rich repeat proteins in plant defences. Adv Bot Res. 1997;24:90–168.

    Google Scholar 

  98. McCann MC, Carpita NC. Designing the deconstruction of plant cell walls. Curr Opin Plant Biol. 2008;11:314–20.

    CAS  PubMed  Google Scholar 

  99. Zhong R, Demura T, Ye Z-H. SND1, a NAC domain transcription factor, is a key regulator of secondary wall synthesis in fibers of Arabidopsis. Plant Cell. 2006;18:3158–70.

    PubMed Central  CAS  PubMed  Google Scholar 

  100. Nogueira FT, Schlögl PS, Camargo SR, Fernandez JH. SsNAC23, a member of the NAC domain protein family, is associated with cold, herbivory and water stress in sugarcane. Plant Sci. 2005;169:93–106.

    CAS  Google Scholar 

  101. Abe M, Abe K, Kuroda M, Arai S. Corn kernel cysteine proteinase inhibitor as a novel cystatin superfamily member of plant origin. Eur J Biochem. 1992;209:933–7.

    CAS  PubMed  Google Scholar 

  102. McMullen MD, Frey M, Degenhardt J. Genetics and biochemistry of insect resistance in maize. In: Bennetzen JL, Hake SC, editors. Handbook of maize: Its biology. NY, USA: Springer; 2009. p. 271–89.

    Google Scholar 

  103. Ussuf K, Laxmi N, Mitra R. Proteinase inhibitors: plant-derived genes of insecticidal protein for developing insect-resistant transgenic plants. Curr Sci India. 2001;80:847–53.

    CAS  Google Scholar 

  104. Fabrick J, Behnke C, Czapla T, Bala K, Rao A, Kramer K, et al. Effects of a potato cysteine proteinase inhibitor on midgut proteolytic enzyme activity and growth of the southern corn rootworm, Diabrotica undecimpunctata howardi (Coleoptera: Chrysomelidae). Insect Biochem Mol Biol. 2002;32:405–15.

    CAS  PubMed  Google Scholar 

  105. Zhang H, Hedhili S, Montiel G, Zhang Y, Chatel G, Pré M, et al. The basic helix-loop-helix transcription factor CrMYC2 controls the jasmonate-responsive expression of the ORCA genes that regulate alkaloid biosynthesis in Catharanthus roseus. Plant J. 2011;67:61–71.

    CAS  PubMed  Google Scholar 

  106. Niu Y, Figueroa P. Characterization of JAZ-interacting bHLH transcription factors that regulate jasmonate responses in Arabidopsis. J Exp Bot. 2011;62:2143–54.

    PubMed Central  CAS  PubMed  Google Scholar 

  107. Sasaki-Sekimoto Y, Jikumaru Y, Obayashi T, Saito H, Masuda S, Kamiya Y, et al. Basic Helix-Loop-Helix transcription factors JASMONATE-ASSOCIATED MYC2-LIKE1 (JAM1), JAM2, and JAM3 are negative regulators of jasmonate responses in Arabidopsis. Plant Physiol. 2013;163:291–304.

    PubMed Central  CAS  PubMed  Google Scholar 

  108. Qi T, Huang H, Wu D, Yan J, Qi Y, Song S, et al. Arabidopsis DELLA and JAZ Proteins Bind the WD-Repeat/bHLH/MYB Complex to Modulate Gibberellin and Jasmonate Signaling Synergy. Plant Cell. 2014;26:1118–33.

    CAS  PubMed  Google Scholar 

  109. Song S, Qi T, Fan M, Zhang X, Gao H, Huang H, et al. The bHLH subgroup IIId factors negatively regulate jasmonate-mediated plant defense and development. PLoS Genet. 2013;9:e1003653.

    PubMed Central  CAS  PubMed  Google Scholar 

  110. Xia H, Yandeau-Nelson M, Thompson DB, Guiltinan MJ. Deficiency of maize starch-branching enzyme i results in altered starch fine structure, decreased digestibility and reduced coleoptile growth during germination. BMC Plant Biol. 2011;11:1–13.

    Google Scholar 

  111. Dunwell JM, Gibbings JG, Mahmood T, Saqlan Naqvi S. Germin and germin-like proteins: evolution, structure, and function. Crit Rev Plant Sci. 2008;27:342–75.

    CAS  Google Scholar 

  112. Davidson RM, Reeves PA, Manosalva PM, Leach JE. Germins: A diverse protein family important for crop improvement. Plant Sci. 2009;177:499–510.

    CAS  Google Scholar 

  113. Breen J, Bellgard M. Germin-like proteins (GLPs) in cereal genomes: gene clustering and dynamic roles in plant defence. Funct Integr Genomics. 2010;10:463–76.

    CAS  PubMed  Google Scholar 

  114. Armstrong JS, Abdel-Mageed H, Fokar M, Allen R, Adamczyk Jr JJ. Dietary effects of cotton tissue expressing germin like protein on beet armyworm (Lepidoptera: Noctuidae) growth, survival and pupation. Fla Entomol. 2013;96:693–700.

    CAS  Google Scholar 

  115. Ryser U, Schorderet M, Guyot R, Keller B. A new structural element containing glycine-rich proteins and rhamnogalacturonan I in the protoxylem of seed plants. J Cell Sci. 2004;117:1179–90.

    CAS  PubMed  Google Scholar 

  116. Ringli C, Keller B, Ryser U. Glycine-rich proteins as structural components of plant cell walls. Cell Mol Life Sci. 2001;58:1430–41.

    CAS  PubMed  Google Scholar 

  117. Mangeon A, Junqueira RM, Sachetto-Martins G. Functional diversity of the plant glycine-rich proteins superfamily. Plant Signal Behav. 2010;5:99–104.

    PubMed Central  CAS  PubMed  Google Scholar 

  118. Showalter AM. Structure and function of plant cell wall proteins. Plant Cell. 1993;5:9.

    PubMed Central  CAS  PubMed  Google Scholar 

  119. Sturm A. A wound-inducible glycine-rich protein from Daucus carota with homology to single-stranded nucleic acid-binding proteins. Plant Physiol. 1992;99:1689–92.

    PubMed Central  CAS  PubMed  Google Scholar 

  120. Bohlmann H, Broekaert W. The role of thionins in plant protection. Crit Rev Plant Sci. 1994;13:1–16.

    CAS  Google Scholar 

  121. De Coninck B, Cammue B, Thevissen K. Modes of antifungal action and in planta functions of plant defensins and defensin-like peptides. Fungal Biol Rev. 2013;26:109–20.

    Google Scholar 

  122. Shiau Y-S, Horng S-B, Chen C-S, Huang P-T, Lin C, Hsueh Y-C, et al. Structural analysis of the unique insecticidal activity of novel mungbean defensin VrD1 reveals possibility of homoplasy evolution between plant defensins and scorpion neurotoxins. J Mol Recognit. 2006;19:441–50.

    CAS  PubMed  Google Scholar 

  123. Bohlmann H, Apel K. Thionins. Annu Rev Plant Biol. 1991;42:227–40.

    CAS  Google Scholar 

  124. Kessler A, Baldwin IT. Plant responses to insect herbivory: The emerging molecular analysis. Annu Rev Plant Biol. 2002;53:299–328.

    CAS  PubMed  Google Scholar 

  125. Hershko A, Ciechanover A. The ubiquitin system for protein degradation. Annu Rev Biochem. 1992;61:761–807.

    CAS  PubMed  Google Scholar 

  126. Zeng L-R, Vega-Sánchez ME, Zhu T, Wang G-L. Ubiquitination-mediated protein degradation and modification: an emerging theme in plant-microbe interactions. Cell Res. 2006;16:413–26.

    CAS  PubMed  Google Scholar 

  127. Kelley DR, Estelle M. Ubiquitin-mediated control of plant hormone signaling. Plant Physiol. 2012;160:47–55.

    PubMed Central  CAS  PubMed  Google Scholar 

  128. Santner A, Estelle M. The ubiquitin-proteasome system regulates plant hormone signaling. Plant J. 2010;61:1029–40.

    PubMed Central  CAS  PubMed  Google Scholar 

  129. Igarashi D, Tsuda K, Katagiri F. The peptide growth factor, phytosulfokine, attenuates pattern-triggered immunity. Plant J. 2012;71:194–204.

    CAS  PubMed  Google Scholar 

  130. Rodriguez PA. Analysis of AtPSKR1, an LRR receptor protein kinase, and other PSK-signalling components in plant defence responses. In: PhD Thesis. Tübingen, Germany: University of Tübingen, Mathematics and Natural Sciences; 2011.

    Google Scholar 

  131. Angelini R, Tisi A, Rea G, Chen MM, Botta M, Federico R, et al. Involvement of polyamine oxidase in wound healing. Plant Physiol. 2008;146:162–77.

    PubMed Central  CAS  PubMed  Google Scholar 

  132. Wimalasekera R, Tebartz F, Scherer GF. Polyamines, polyamine oxidases and nitric oxide in development, abiotic and biotic stresses. Plant Sci. 2011;181:593–603.

    CAS  PubMed  Google Scholar 

  133. Rizhsky L, Hallak-Herr E, Van Breusegem F, Rachmilevitch S, Barr JE, Rodermel S, et al. Double antisense plants lacking ascorbate peroxidase and catalase are less sensitive to oxidative stress than single antisense plants lacking ascorbate peroxidase or catalase. Plant J. 2002;32:329–42.

    CAS  PubMed  Google Scholar 

  134. Willekens H, Chamnongpol S, Davey M, Schraudner M, Langebartels C, Van Montagu M, et al. Catalase is a sink for H2O2 and is indispensable for stress defence in C3 plants. EMBO J. 1997;16:4806–16.

    PubMed Central  CAS  PubMed  Google Scholar 

  135. Kaur R, Gupta AK, Taggar GK. Role of catalase, H2O2 and phenolics in resistance of pigeonpea towards Helicoverpa armigera (Hubner). Acta Physiol Plant. 2014;36:1513–27.

    CAS  Google Scholar 

  136. Bajsa J, Pan Z, Duke SO. Serine/threonine protein phosphatases: Multi-purpose enzymes in control of defense mechanisms. Plant Signal Behav. 2011;6:1921.

    PubMed Central  CAS  PubMed  Google Scholar 

  137. Vyroubalová Š, Václavíková K, Turečková V, Novák O, Šmehilová M, Hluska T, et al. Characterization of new maize genes putatively involved in cytokinin metabolism and their expression during osmotic stress in relation to cytokinin levels. Plant Physiol. 2009;151:433–47.

    PubMed Central  PubMed  Google Scholar 

  138. Smigocki A, Neal Jr J, McCanna I, Douglass L. Cytokinin-mediated insect resistance in Nicotiana plants transformed with the ipt gene. Plant Mol Biol. 1993;23:325–35.

    CAS  PubMed  Google Scholar 

  139. Giron D, Frago E, Glevarec G, Pieterse CM, Dicke M. Cytokinins as key regulators in plant-microbe-insect interactions: connecting plant growth and defence. Funct Ecol. 2013;27:599–609.

    Google Scholar 

  140. De Los CG, Gianola D, Allison DB. Predicting genetic predisposition in humans: the promise of whole-genome markers. Nat Rev Genet. 2010;11:880–6.

    Google Scholar 

  141. 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.

    CAS  PubMed  Google Scholar 

  142. Eizaguirre M, Albajes R. Diapause induction in the stem corn borer Sesamia nonagrioides (Lepidoptera, Noctuidae). Entomol Gen. 1992;17:277–83.

    Google Scholar 

  143. SAS Institute Inc. SAS 9.3 Guide to software updates. Cary, NC, USA, SAS Institute Inc. 2011; 314 p

  144. Holland JB, Nyquist WE, Cervantes-Martínez CT. Estimated an interpreting heritability for plant breeding: An update. In: Janick J, editor. Plant Breeding Reviews, vol. Volume 22. Hoboken, New Jersey, USA: Jonh Wiley & Sons; 2003. p. 9–112.

    Google Scholar 

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

    Google Scholar 

  146. 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:2633–5.

    CAS  PubMed  Google Scholar 

  147. Yu J, Pressoir G, Briggs WH, Vroh I, Yamasaki I, Doebley JF, et al. A unified mixed-model method for association mapping that accounts for multiple levels of relatedness. Nat Genet. 2006;38:203–8.

    CAS  PubMed  Google Scholar 

  148. Valdar W, Holmes CC, Mott R, Flint J. Mapping in structured populations by resample model averaging. Genetics. 2009;182:1263–77.

    PubMed Central  PubMed  Google Scholar 

  149. Panagiotou OA, Ioannidis JP. What should the genome-wide significance threshold be? Empirical replication of borderline genetic associations. Int J Epidemiol. 2012;41:273–86.

    PubMed  Google Scholar 

  150. R Core Team. R: A languange and environment for statistical computing. In: Book R: A languange and environment for statistical computing (Editor ed.^eds.). Vienna, Austria: City: R Foundation for Statistical Computing; 2013.

    Google Scholar 

Download references

Acknowledgements

This work was supported by the National Plan for Research and Development of Spain (projects AGL2012-33415). L.F. Samayoa acknowledges a contract JAE-Predoc from the Spanish Council for Scientific Research (CSIC). We also acknowledges to the Super Computation Center of Galicia (CESGA) for allowed to use their serves to do the preliminary analysis of this study. We thank Mike McMullen for inspiring this study and reviewing the manuscript.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Luis Fernando Samayoa.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

LFS carried out field experiments, performed statistical and bioinformatics analysis of the data and made the draft of the manuscript. AB and RAM conceived the study, assisted LFS in field experiment design, and data collection and analysis. BAO and JBH prepared some arrays of data and assisted LFS in statistical and bioinformatics analysis. All authors read and approved the final manuscript.

Additional files

Additional file 1:

Table S1. Mean squares and heritability estimates for MCB pest resistance traits evaluated in an association panel in three years. Table S2 Mean squares and heritability estimates for agronomic traits evaluated in an association panel in three years. Table S3 Phenotypic data of three resistance traits to MCB attack and three agronomic traits in a maize diversity panel. Table S4 Annual and average means for traits related to resistance to MCB and three agronomic traits. Table S5 Summary of the SNPs significantly associated to plant height. SNP identification (SNP ID), additive effect and allelic variants for the SNP, proportion of total variance explained by the SNPs significantly associated with plant height (PH) and significance values for the association between the SNP and the phenotype (P-value and RMIP). Table S6 Summary of the SNPs significantly associated to days to anthesis. SNP identification (SNP ID), additive effect and allelic variants for the SNP, proportion of total variance explained by the SNPs significantly associated with days to anthesis (DTA) and significance values for the association between the SNP and the phenotype (P-value and RMIP). Table S7 Summary of the SNPs significantly associated to days to silking. SNP identification (SNP ID), additive effect and allelic variants for the SNP, proportion of total variance explained by the SNPs significantly associated with days to silking (DTS) and significance values for the association between the SNP and the phenotype (P-value and RMIP). Table S8 Linked disequilibrium (r 2) between the ten SNPs significantly associated to TL made by MCB.

Additional file 2:

Genotypic data of a maize diversity panel with more than 245 kb SNP (GBS + Illumina). The data were filtered to exclude SNPs with more than 20% missing genotype data and minor allele frequency (MAF) less than 5%. Also the heterozygous are excluded. The array has 290 columns and 246,477 rows.

Rights and permissions

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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Samayoa, L.F., Malvar, R.A., Olukolu, B.A. et al. 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 15, 35 (2015). https://doi.org/10.1186/s12870-014-0403-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-014-0403-3

Keywords