- Research article
- Open Access
Functional mechanisms of drought tolerance in subtropical maize (Zea mays L.) identified using genome-wide association mapping
BMC Genomics volume 15, Article number: 1182 (2014)
Earlier studies were focused on the genetics of temperate and tropical maize under drought. We identified genetic loci and their association with functional mechanisms in 240 accessions of subtropical maize using a high-density marker set under water stress.
Out of 61 significant SNPs (11 were false-discovery-rate-corrected associations), identified across agronomic traits, models, and locations by subjecting the accessions to water stress at flowering stage, 48% were associated with drought-tolerant genes. Maize gene models revealed that SNPs mapped for agronomic traits were in fact associated with number of functional traits as follows: stomatal closure, 28; flowering, 15; root development, 5; detoxification, 4; and reduced water potential, 2. Interactions of these SNPS through the functional traits could lead to drought tolerance. The SNPs associated with ABA-dependent signalling pathways played a major role in the plant’s response to stress by regulating a series of functions including flowering, root development, auxin metabolism, guard cell functions, and scavenging reactive oxygen species (ROS). ABA signalling genes regulate flowering through epigenetic changes in stress-responsive genes. ROS generated by ABA signalling are reduced by the interplay between ethylene, ABA, and detoxification signalling transductions. Integration of ABA-signalling genes with auxin-inducible genes regulates root development which in turn, maintains the water balance by regulating electrochemical gradient in plant.
Several genes are directly or indirectly involved in the functioning of agronomic traits related to water stress. Genes involved in these crucial biological functions interacted significantly in order to maintain the primary as well as exclusive functions related to coping with water stress. SNPs associated with drought-tolerant genes involved in strategic biological functions will be useful to understand the mechanisms of drought tolerance in subtropical maize.
Drought at the flowering stage produces infertile pollen in maize (Zea mays L.) , resulting in substantially lower yield. A shorter anthesis-to-silking interval (ASI) was observed under drought in a lowland tropical maize population , and the interval was considered to be an indirect selection criterion for grain yield . Grain yield, a complex trait dependent on several environmental factors , was increased through genetic modifications despite stress .
Several researchers have identified drought-related genes in different tissues that regulate molecular and physiological responses under stress [6–8]. In the ears and silks of maize, some genomic regions that control the levels of abscisic acid (ABA) and sugar showed the signal transduction involved in stress-related pathways and regulated kernel size and productivity of the plant under drought . Different tissues showed the expression of an NAC-transcription-factor-encoding gene , which was observed to be strongly associated with the relative ear position, 100-kernel weight, and flowering time across the experimental maize population . Some drought response has also been recorded recently in reproductive and leaf meristem tissues in maize .
Furthermore, a high-resolution association map created using a high-density marker set has the potential to unravel stress-associated genetic variability in a genome . Considerable efforts have been made to understand the association of SNPs with different phenotypic traits in maize [10–12], and strong SNP associations have been identified for flowering time , kernel shape , 100-kernel weight , and kernel quality . Moreover, target genes for crop improvement have been successfully identified using genome-wide association (GWA) mapping in maize [6, 10, 12].
In this study, we assembled a GWA mapping panel using elite subtropical maize genotypes and phenotyped it under water stress (WS) at flowering time at multiple locations to identify the SNPs associated with key agronomic traits. The association of genetic loci with the agronomic traits and their recurrence level across different locations and models were examined, to identify the SNPs significantly associated with the candidate genes that interact to confer tolerance to WS.
An association mapping panel of 240 subtropical maize lines was subjected to WS and the performance of the stressed maize lines was compared with that of the lines grown under a well-watered (WW) environment. The panel was phenotyped for anthesis-to-silking interval (ASI), ear length (EL), ear girth (EG), kernels per row (KR), the number of kernel rows (KRN), 100-kernel weight (HKW), and grain yield (GY) at three locations: Indian Agriculture Research Institute (IARI), New Delhi; Acharya N G Ranga Agricultural University (ANGRAU), Hyderabad; and Agricultural Research Station (ARS), Karimnagar for two years. The overall performance of the hybrids under WS was significantly impaired [Additional file 1: Table S1]. The phenotypic correlation coefficient showed that GY and the traits contributing to it were positively and significantly correlated with each other, and ASI was negatively but significantly correlated with other agronomic traits [Additional file 2: Table S2].
Genome-wide association mapping
Genome-wide association models
We detected GWA signals from GenABEL and Genome Association and Prediction Integrated Tool (GAPIT) models. Model fitness varied over the data sets (traits + models + locations) and an example from GenABEL is given below. For ASI–Hyderabad, a higher model (Q10 + K + Admixture) from GenABEL showed greater fitness than the lower hierarchy models (K and Q10 + K) did. The genomic inflation factor (λ) was 1.04 in the latter model but was reduced to 1.02 when admixture was added as a covariate (Figure 1a); Q10 + K and Q10 + K + Admixture from GenABEL had equal effect on correcting cryptic relationship for ASI–Karimnagar (Figure 1b); and all the models proved equally good for GY–Hyderabad (Figure 1c).
Comparison of GWA signals: WW and WS
Of 61 SNPs strongly (P < 5.8 × 10−4) associated with WS, 11 were false-discovery-rate (FDR)-corrected associations and the rest (50) were consistent across the data sets (traits + models + locations). The maximum number (16) of association signals were identified on chromosome 5 followed by those on chromosome 3 (15), whereas chromosomes 2 and 8 shared the least number of significant associations (2) across traits measured under WS.
Genome-wide analysis was also assessed for the WW data sets, which mapped 70 significant associations (P < 5.8 × 10−4) including 9 FDR-corrected associations. The maximum number of GWA signals for WW was contributed by the SNPs mapped on chromosome 10 (15) followed by those on chromosomes 1 (12) and 7 (10), whereas the minimum number (2) of GWA signals was shared by the SNPs mapped on chromosome 8. The recurrence of these signals across all data sets (traits + models + locations) varied with those observed in WS. Also the stringency of these results was observed through the number of significantly associated SNPs for EG trait being 17 in WS and 24 in WW. But, for traits ASI, GY, and KR, the number of GWA signals in WW were equivalent to that in WS. SNPs significantly associated in both conditions were mapped near both drought tolerant genes and regulatory genes having their importance in both stressed and normal environment–PZE-107021672 and PZE-107021673 for GY; PZE-101135368 for ASI and KR; PZE-107110985 for ASI, GY, and KR. However rest of the SNPs were also significantly associated in WW and WS but were specific to environmental locations.
GWA signals: WS
On chromosome 1, a significant signal (PZE-101135368) associated with WS-ASI was located 12 kb from the nuclear factor-YA transcription factor (NF-YA) [Additional file 3: Table S3]. This SNP also showed an FDR-corrected P value of 5.4 × 10−9 for WS-ASI and was repeated 75% of the times across all data sets (traits + models + locations) under WS (Figure 2). For ASI–New Delhi, PZE-101135368 was significantly associated (at 5% FDR) with a stronger P value of 5.4 × 10−9 for WS [Additional file 3: Table S3] than that for the non-significant signal observed under WW (P = 1.4 × 10−4) [Additional file 4: Figure S1]. In addition, this locus was co-localized with a quantitative trait locus (QTL) earlier mapped for ASI by Almeida et al. Furthermore, significant associations for WS also comprised DnaJ-49-like chaperone protein and a domain of unknown function (231) that were co-localized with two other QTLs for ASI on chromosome 1 under WS .
Two FDR-corrected associations for WS-ASI, one on chromosome 7 (PZE-107110985) and the other on chromosome 10 (PZE-110016300), were found near a gene encoding a basic helix-loop-helix (bHLH) transcription factor [Additional file 3: Table S3]. The first locus identified at 92 kb from bHLH was also significantly associated with WS-GY, -EG, -HKW, -EL, -KRN, and -KR at specific locations [Additional file 3: Table S3]. The second locus, mapped at 118 kb from bHLH, also showed strong associations with WS-HKW, -KRN, and -KR [Additional file 3: Table S3] and significant associations for ASI across all locations and mean data (Figures 2 and 3a), while the rest of the traits were confined to specific locations (Figure 2).
SNP PZE-103149505 mapped on chromosome 3, indicated strong GWA signals for WS-ASI across multiple locations (Figure 2). These signals showed a stronger and significant association under WS than that for WW, in which the P values were less significant and less repeatable. This SNP was mapped 30 kb from a stress-related gene encoding NADP-malic enzyme [Additional file 3: Table S3].
SNPs mapped on chromosome 5 showed five FDR-corrected association signals for WS-GY [Additional file 3: Table S3]: three SNPs (PZE-105073248, PZE-105073274, and PZE-105073275) were consistent across all locations, whereas two (PZE-105073295 and PZE-105076764) were specific to the locations and models (Figure 3b). Of these five SNPs, three (PZE-105073248, PZE-105073275, and PZE-105073274) were clustered in a 131 kb region, which enclosed genes encoding MYB-related (MYB) transcription factor and squamous promoter binding protein (SBP) transcription factor [Additional file 3: Table S3]. This cluster was identified 65% of the times across all WS data sets, including GY (Figure 2), and showed strong associations with WS-GY that were more recurrent across locations and models than those in the WW data sets. In addition to this, two FDR-corrected SNPs (PZE-107021673 and PZE-107110985) on chromosome 7 were also strongly associated with WS-GY. The first SNP was mapped nearer to the glycerol-3-phosphate dehydrogenase gene and an ethylene responsive factor (ERF) transcription factor, and the second was close to the bHLH encoding gene [Additional file 3: Table S3]. Another SNP on chromosome 9, 16 kb from the MYB, was also associated with WS-GY [Additional file 3: Table S3]. All SNPs on chromosomes 7 and 9 also had a significant association with WS-EG [Additional file 3: Table S3].
Two SNPs on chromosome 5 (PZE-105025225 and PZE-105076764) mapped close to a gene encoding a C2H2-type zinc finger (C2H2) transcription factor were independently associated for WS-ASI and -GY when measured at Hyderabad. The first SNP (PZE-105025225), 31 kb from C2H2, showed a strong association (P = 2.4 × 10−5) for WS-ASI. The second SNP (PZE-105076764), 7 kb from C2H2, was associated with a strong signal (P = 1.2 × 10−5) for WS-GY [Additional file 3: Table S3], which was non-significant under WW. Also, SNPs mapped on chromosomes 2 (PZE-102185808) and 7 (PZE-107110985) were associated with both WS-ASI and -GY (Figure 3a, b) where the association of the first locus with ASI was more significant in WS than that in non-significant signals observed under WW at Hyderabad. The first locus was mapped near the genes encoding MYB, a bZIP transcription factor, and a zinc finger-homeodomain (ZF-HD) transcription factor, and the second was near bHLH [Additional file 3: Table S3].
For WS-HKW, four strong GWA signals were observed on different chromosomes close to the stress-related genes. The first two loci (PZE-101135368 on chromosome 1 and PZE-105025225 on chromosome 5) were observed 12 kb and 20 kb from NF-YA respectively [Additional file 3: Table S3]. The third locus (PZE-106063888 on chromosome 6), also showing a significant association for WS-HKW (P = 6.4 × 10−5), was mapped close to one gene, namely vacuolar proton ATPase pump (V-type H+pump), and two transcription factors, namely NAC and calmodulin-binding transcription activator (CAMTA) [Additional file 3: Table S3]. The fourth locus (PZE-103035609 mapped on chromosome 3) was close to a universal stress protein [Additional file 3: Table S3].
Metabolic mechanisms common to WS and WW conditions
Strong GWA signals were detected near the genes responsible for drought tolerance as well as for plant metabolism across both WS and WW environments. This commonality highlights the importance of genes that are responsible for maintaining yield and were associated with yield-related traits under both the environments. Two SNPs (PZE-107021672 and PZE-107021673) on chromosome 7 were strongly associated with GY for both WS (Figure 2) and WW data sets [Additional file 4: Figure S1, Additional file 5: Figure S2, Additional file 6: Figure S3, and Additional file 7: Figure S4]. These SNPs showed a drought-tolerant ERF transcription factor in their vicinity, which promotes ABA-dependent stomatal closure under WS  (Figure 4). Apart from this drought-tolerant gene, these SNPs also included two more genes encoding MADS-box and glycerol-3-phosphate dehydrogenase (G3PDH), which regulate generalized functions in plants. During root development, the MADS-box gene regulates auxin transport , cell proliferation , the transition to flowering and flower development . G3PDH is part of the mitochondrial glycerol-3-phosphate shuttle system that is induced at different stages of plant development when the demand for redox adjustment is high . G3PDH is also one of the enzymes involved in glycerol metabolism, where the role of glycerol in regulating root development through multiple pathways in Arabidopsis has been studied recently .
A strong GWA signal (PZE-101135368) for ASI and KRN was identified for both WS [Additional file 3: Table S3] and WW across all data sets although the P value varied with the data set [Additional file 4: Figure S1, Additional file 5: Figure S2, Additional file 6: Figure S3, and Additional file 7: Figure S4]. This association pointed to a drought-tolerant NF-YA transcription factor 12 kb from the SNP. Under WS, NF-YA promotes ABA-dependent stomatal closure, which contributes to drought tolerance  (Figure 4). NF-YA is also important to several other functions including male gametogenesis, embryogenesis, and seed development under WW conditions . The consistency of these GWA signals shows that NF-YA is important to high yields irrespective of whether if water is scarce or abundant.
Another strong GWA signal (PZE-107110985) was that for ASI, GY, and KR across both environments and at all locations. BHLH, a drought-tolerant gene, was mapped close to this SNP, and has been known to interact with two other two drought-tolerant genes, bZIP and MYB, which are responsible for tissue-specific flavonoids production and light responsiveness under normal conditions : bZIP regulates auxin transport  and MYB promotes light signalling transductions including photosynthesis in Arabidopsis thaliana. The importance of these genes under WS was noticed from the specific GWA signals for WS that were mapped close to three drought-tolerant gene families: PZE-103061107 near MYB, PZE-104061181 and SYNGENTA14972 near bHLH, and PZE-104109358 near bZIP. Apart from drought tolerance, these genes can also interact strongly to regulate other generalized functions, as can be inferred from the GWA signals mapped across both the environments.
Mechanisms of drought tolerance
Strong GWA signals were recorded near the interacting drought-tolerant genes under WS [Additional file 3: Table S3] (Figure 4). These mechanisms, namely stomatal closure, detoxification signalling, root development, reduced water potential, and flowering, are influenced at different physiological and functional levels and enhance drought tolerance (Figure 4).
A genomic region encoding two neighbouring SNPs 253 bp apart (PZE-105073274 and PZE-105073275) in high linkage disequilibrium (LD) (r2 = 1)  was seen on chromosome 5, near two CAMTA regulating factors, the MYB and SBP (Figure 5). These SNPs showed significant P values for several traits including GY, EG, EL, and KR under WS [Additional file 3: Table S3]. The extended regions with strong LD exhibited a steep selection and genetic drift  of those traits in the population. CAMTA was also identified near a locus mapped on chromosome 6 that was strongly associated with WS-ASI, -EG, -EL, -KRN, -KR, and -HKW [Additional file 3: Table S3]. CAMTA (at 115 Mb on chromosome 6), MYB, and SBP (at 77.9 Mb on chromosome 5) contributed to stomatal closure, which promotes drought tolerance, presumably through strong inter-chromosomal LD. Other CAMTA-regulating factors such as WRKY and bZIP were mapped near an SNP (PZE-104109358) on chromosome 4 that was strongly associated with WS-ASI, -KRN, and -HKW [Additional file 3: Table S3].
An SNP cluster on chromosome 5 with a significant association for WS-GY and -HKW included the genes encoding C2H2 and MYB [Additional file 3: Table S3]. Although the function of C2H2 in maize is not clear, a C2H2-type zinc finger was shown to play a defensive role in oxidative stress generated under different forms of environmental stress in Arabidopsis. Moreover, a gene of the MYB family is known to integrate ABA and auxin signalling under WS .
In the present study, two SNPs were identified, one on chromosome 3 and one on chromosome 4, with strong associations for WS-ASI and WS-KRN [Additional file 3: Table S3] (Figure 2). The SNP on chromosome 3 encompassed the gene encoding 9-cis-epoxycarotenoid dioxygenase (NCED) whereas the one on chromosome 4 was localized at a point close to a molybdenum cofactor sulfurase encoding gene [Additional file 3: Table S3]. Because of possible inter-chromosomal LD, this region may contribute to ABA-dependent stomatal closure. This defensive mechanism can be explained by an integrative ABA synthesis, accumulation , and signalling in guard cells . We also identified two SNPs, one on chromosome 3 and another on chromosome 10, close to an NADP-malic enzyme encoding gene [Additional file 3: Table S3], with strong P values for WS-ASI and -GY, respectively (Figure 3a, b). This gene promotes stomatal closure, which induces tolerance during WS, by removing malate from guard cells .
Detoxification signalling is a major response to environmental forms of stress, and WS is known to stimulate oxidative stress in chloroplasts, mitochondria, and peroxisomes , thereby leading to the production of reactive oxygen species (ROS) that are sensitive to ABA-signalling pathways  and induce programmed cell death in plants . Therefore, removal of free oxygen radicals may represent the primary defence mechanism in drought tolerance. Glutathione-S-transferase (GST) and glutathione peroxidase (GP), which are responsible for ROS homeostasis, change turgor pressure through an ABA-mediated pathway [35, 36] in order to sustain the plant under drought. These genes were observed to be strongly associated with WS-GY, -HKW, and -KRN [Additional file 3: Table S3], thereby highlighting the role of ROS homeostasis under WS (Figure 4).
Phospholipid hydrolysis is also known to play a role in detoxification signalling under WS [37, 38]. In the present study, a genomic region (chromosome 3) encoding phospholipase D (PLD) displayed a strong and WS-specific association signal for GY [Additional file 3: Table S3]. PLD is a signalling messenger that facilitates stomatal closure, cell viability, and enhanced root growth during stress .
Under WS, auxin signalling results in a deeper root system, which increases water use efficiency [39, 40]. We noticed two strong association signals in the genomic regions that included the genes for auxin transport (NAC) and auxin biosynthesis (flavin monooxygenase) on chromosomes 6 and 10 respectively [Additional file 3: Table S3]. These regions showed high P values for WS-ASI, -KR, and -HKW, suggesting that maintenance of the root system was in fact associated with flowering and seed development stages.
Reduced water potential
Vacuolar proton pumps reduce water potential under WS: their expression in the root system enhances water-absorbing efficiency . The strong signal of the vacuolar proton pyrophosphatase pump (V-type PPase H+pump) for WS-KRN and that of the V-type ATPase H+pump for WS-ASI, -EG, -HKW, -EL, -KRN, and -KR [Additional file 3: Table S3] establish a proton electrochemical gradient across the vacuolar membrane [41, 42], thereby enhancing ion homeostasis, which is responsible for maintaining osmotic balance and in turn, coping with drought (Figure 4).
Drought tolerance also involves changes in flowering through epigenetic mechanisms and the ABA-dependent pathway – an observation borne out by the present study (Figure 4) –and DNA cytosine methylation and histone modifications (ABA levels are regulated by histone) are innate to epigenomes [43, 44]. Our study identified 15 significant SNPs [Additional file 3: Table S3] close to genes encoding DNA-cytosine-5-methyltransferase (DCM) and MYB transcription factor that are involved in epigenetic mechanisms and ABA signalling pathways under WS, thereby pointing to their role in drought tolerance.
Functional relationship of stress-responsive genes
Interactions of several functional mechanisms associated with drought tolerance are discussed below. Epigenetic changes such as DNA methylation respond differently at tissue and developmental levels. In a study of the association of methylation-specific sites with drought tolerance, methylation was reduced in response to water stress at the tillering stage whereas drought-induced DNA methylation and demethylation were higher at the tillering stage than at the booting and heading stages . Epigenetic changes in MYB, an ABA-signalling transcription factor, have been found to be specific to tissues including pericarp and cob in maize . These epigenetic changes include site-specific methylation patterns, which have been studied in ABA-signalling genes including the MYB gene family, where stomatal development genes such as MYB44 are demethylated and expressed in response to WS . MYB proteins also promote stamen development and ensure adaptation to WS . Therefore the epigenetic changes in ABA-signalling genes can possibly lead to changes in flowering under WS (Figure 4).
ROS are produced in an ABA-signalling pathway in which two ABA-signalling transcription factors C2H2 and MAPK are negatively regulated by ERF, which is repressed under WS . ROS are scavenged by a detoxification signalling process , in which PLD is reported to influence ABA responses and opening and closing of stomata through bifurcating pathways . When PLD activity was studied under water stress, three ABA-signalling genes, namely MYB, NAC (PZE-108058161; chromosome 8), and bHLH, showed differential regulation in both wild-type and antisense-PLD Arabidopsis. Therefore, networking among ABA and PLD signals is one of the important defensive mechanisms against WS. This suggests that the ROS generated by ABA signalling are reduced by the interplay between ethylene signalling, ABA signalling, and detoxification signalling under WS.
Guard cell signal transductions regulate CO2 influx for photosynthesis as well as water loss through stomatal closure under WS [31, 32]: ABA signalling promotes the closing of stomata, which reduces CO2 intake in guard cells and thereby decreases the rate of photosynthesis. To maintain the rate of photosynthesis under WS, plants adopt an alternative pathway for delivering CO2. This alternative pathway includes malate degradation, catalyzed by a respiratory enzyme, NADP-malic enzyme, which releases CO2 as a secondary source for CO2 influx under WS .
ABA is a hormone released in response to various forms of abiotic stress and regulates tolerance through a network of interconnected genes. The mechanisms of drought tolerance are stimulated by both ABA-dependent and ABA-independent pathways: C2H2, for instance, is part of an ABA-dependent pathway as well as that of an ABA-independent pathway. C2H2 zinc finger proteins with ERF-associated amphiphilic repression (EAR) motifs have been reported to be responsive to ABA  and to WS .These findings suggest that ABA and ethylene interact under drought and enhance tolerance by promoting stomatal closure. Under WS, C2H2 acts as a transcriptional regulator through an ABA-independent pathway : ABA-independent C2H2 triggers the genes that are related to H2O2 homeostasis and thus decreases the amount of H2O2 in guard cells . However, stomatal closure as part of the drought response is promoted by ABA and H2O2 accumulation in guard cells .
The electrochemical gradient generated by the proton pump (V-type PPase H+pump mapped near PZE-110027802; chromosome 10) promotes the secondary active transport of sugar molecules (hexokinase 3 mapped near PZE-108106212; chromosome 8) into the vacuole. This helps the plant in maintaining its internal water balance . In addition, the proton pump has also been reported in roots , which suggest that this pump regulates water balance in roots and thus contributes to drought tolerance. It is well understood that the water balance in plants is interrelated to root development, where the auxin-inducible genes that promote root development are regulated by the C2H2. This suggests that the integration of root development and ABA signalling molecule C2H2 allows the maintenance of an electrochemical gradient in the plant.
Genome-wide association analysis, using multiple locations and two models, identified SNPs from agronomic traits linked with the genes directly or indirectly associated with drought tolerance in subtropical maize lines. These genes uncovered physiological responses and molecular mechanisms related to drought tolerance. Genes governing several functional traits were identified, including stomatal closure, reduced water potential, root development, signalling pathways, and flowering. These genes interact extensively to help the plant cope with drought. SNPs and their functional association with several drought-responsive genes will be useful in elucidating the mechanism of drought tolerance in subtropical maize.
Our experiment was based on a panel of 240 elite inbred lines of subtropical maize  from several Indian breeding programmes and also included genotypes from CIMMYT, Mexico adapted to subtropical climates. The lines were separated into three groups, namely early, medium, and late, based on the number of days to flowering.
Field experiments were laid out in an alpha-lattice design at three locations: IARI, New Delhi (28°N 77°E; 229 Amsl); Maize Research Centre, Acharya N G Ranga Agricultural University, Hyderabad (17°N 78°E; 536 Amsl), and ARS, Acharya N G Ranga Agricultural University, Karimnagar (18° N 79°E; 264 Amsl) during the post-rainy seasons of 2010/11 and 2011/12. The design comprised 16 incomplete blocks, each made up of 15 plots with three replicates. Sowing was staggered to ensure that all the lines – whether early, medium, or late – are exposed to stress at the same stage of their growth. WS was induced at the flowering stage and continued until grain filling by withholding irrigation, which was resumed in all the maturity groups at the end of the stress treatment. The following observations were recorded for all the 240 lines: anthesis-to-silking interval (ASI, in days), ear girth (EG, in centimetres), ear length (EL, in centimetres), the number of kernels per row (KR), the number of kernel rows (KRN), 100-kernel weight (HKW, in grams), and grain yield (GY, kilograms per plot).
All the 240 genotypes were genotyped with Infinium Maize SNP50 BeadChip (Illumina, San Diego, California, USA) containing 56 110 SNPs published earlier by Nepolean et al. A set of 29 619 high-quality SNPs  was used for association analysis.
Mixed model analysis using restricted maximum likelihood approach was performed to estimate the best linear unbiased predictors (BLUPs) of genotypes at each location. Plot-level data from each location were analyzed using the following linear mixed-effects model:
where y ijk is the observed value of the kth genotype in the jth incomplete block within the ith replication, μ is the grand mean, r i the main effect of the ith replication, r(b) ij the nested effect of the jth incomplete block in ith replication, g k is the main effect of the kth genotype, and ϵ ijk is the error term of each y ijk with N(0,σ2). At each location, the nested blocks and genotype effects were treated as random and replicated as fixed effects. The residual diagnostic plots indicated that the data satisfied the assumption of normality. To understand the effect of genotypes across locations in greater detail, the data across the three locations were analyzed using the model
where l m is the main effect of mth location and (gl) mk is the interaction effect of kth genotype in mth location. Residual variances of individual locations were modelled in combined analysis using a mixed model procedure. BLUPs for the genotypes were also estimated across combined locations. Data were analyzed using a proc mixed procedure of the software package SAS ver. 9.3 for Windows (SAS Institute Inc.) .
Genome-wide association analysis
Two models, Q10 + K + Admixture from GenABEL and Q10 + K + Admixture from GAPIT were used to identify significant SNPs from BLUPs. Here, admixture values were obtained from earlier data for the same panel . In GenABEL, a mixed linear model (MLM) was employed with polygenic (maximum likelihood method)  and “mmscore” functions . In the present study, the K estimate represents the identity by state (IBS) shared by individuals computing to IBS or kinship matrix generated from GenABEL . The population structure was corrected for the first ten components (Q10) [12, 62] by implementing the “cmdscale” function in GenABEL . P values with 1 degree of freedom were corrected at 5% FDR using the “qvaluebh95” function  in GenABEL. In GAPIT, GWAS was investigated by the compressed mixed linear model (CMLM) approach with “EMMAx” (Efficient Mixed Model Association-eXpedited) , and P values were adjusted at 5% FDR to determine significant associations . To compute kinship matrix, the VanRaden algorithm  was implemented. Admixture and principal component analysis (PCA) were used to correct the population structure.
Best linear unbiased predictors
Calmodulin binding transcription activator
Ethylene responsive factor
Mixed linear model
Quantitative trait locus
Reactive oxygen species
Squamous promoter binding protein
Single nucleotide polymorphism
Schoper JB, Lambert RJ, Vasilas BL: Maize pollen viability and ear receptivity under water and high temperature stress. Crop Sci. 1986, 26: 1029-1033. 10.2135/cropsci1986.0011183X002600050038x.
Edmeades GO, Bolaños J, Hernandez M, Bello S: Causes for silk delay in lowland tropical maize. Crop Sci. 1993, 33: 1029-1035. 10.2135/cropsci1993.0011183X003300050031x.
Magorokosho C, Pixley KV, Tongoona P: Selection for drought tolerance in two tropical maize populations. African Crop Sci J. 2003, 11: 151-161.
Russell WA: Genetic improvements of maize yields. Adv Agron. 1991, 46: 245-298.
Bolaños J, Edmeades GO: The importance of the anthesis-silking interval in breeding for drought tolerance in tropical maize. Field Crops Res. 1996, 48: 65-80. 10.1016/0378-4290(96)00036-6.
Setter TL, Yan J, Warburton M, Ribaut JM, Xu Y, Sawkins M, Buckler ES, Zhang Z, Gore MA: Genetic association mapping identifies single nucleotide polymorphisms in genes that affect abscisic acid levels in maize floral tissues during drought. J Exp Bot. 2011, 62: 701-716. 10.1093/jxb/erq308.
Kakumanu A, Ambavaram MM, Klumas C, Krishnan A, Batlang U, Myers E, Grene R, Pereira A: Effects of drought on gene expression in maize reproductive and leaf meristem tissue revealed by RNA-Seq. Plant Physiol. 2012, 160: 846-867. 10.1104/pp.112.200444.
Pandey N, Ranjan A, Pant P, Tripathi RK, Ateek F, Pandey HP, Patre UV, Sawant SV: CAMTA 1 regulates drought responses in Arabidopsis thaliana. BMC Genomics. 2013, 14: 216-10.1186/1471-2164-14-216.
Lu M, Ying S, Zhang D-F, Shi Y-S, Song Y-C, Wang T-Y, Li Y: A maize stress-responsive NAC transcription factor, ZmSNAC1, confers enhanced tolerance to dehydration in transgenic Arabidopsis. Plant Cell Rep. 2012, 31: 1701-1711. 10.1007/s00299-012-1284-2.
Xue Y, Warburton ML, Sawkins M, Zhang X, Setter T, Xu Y, Grudloyma P, Gethi J, Ribaut JM, Li W, Zhang X, Zheng Y, Yan J: Genome-wide association analysis for nine agronomic traits in maize under well-watered and water-stressed conditions. Theor Appl Genet. 2013, 10.1007/s00122-013-2158-x
Tian F, Bradbury PJ, Brown PJ, Hung H, Sun Q, Flint-Garcia S, Rocheford TR, McMullen MD, Holland JB, Buckler ES: Genome-wide association study of leaf architecture in the maize nested association mapping population. Nat Genet. 2011, 43: 159-162. 10.1038/ng.746.
Riedelsheimer C, Lisec J, Czedik-Eysenberg A, Sulpice R, Flis A, Grieder C, Altmann T, Stitt M, Willmitzer L, Melchinger AE: Genome-wide association mapping of leaf metabolic profiles for dissecting complex traits in maize. Proc Natl Acad Sci U S A. 2012, 109: 8872-8877. 10.1073/pnas.1120813109.
Salvi S, Sponza G, Morgante M, Tomes D, Niu X, Fengler KA, Meeley R, Ananiev EV, Svitashev S, Bruggemann E, Li B, Hainey CF, Radovic S, Zaina G, Rafalski JA, Tingey SV, Miao G-H, Phillips RL, Tuberosa R: Conserved noncoding genomic sequences associated with a flowering-time quantitative trait locus in maize. Proc Natl Acad Sci U S A. 2007, 104: 11376-11381. 10.1073/pnas.0704145104.
Li Q, Li L, Yang XH, Warburton ML, Bai GH, Dai JR, Li JS, Yan JB: Relationship, evolutionary fate and function of two maize orthologous genes of rice GW2 associated with kernel size and weight. BMC Plant Biol. 2010, 10: 143-10.1186/1471-2229-10-143.
Manicacci D, Camus-Kulandaivelu L, Fourmann M, Arar C, Barrault S, Rousselet A, Feminias N, Consoli L, Francès L, Méchin V, Murigneux A, Prioul JL, Charcosset A, Damerval C: Epistatic interactions between opaque2 transcriptional activator and its target gene CyPPDK1 control kernel trait variation in maize. Plant Physiol. 2009, 150: 506-520. 10.1104/pp.108.131888.
Almeida GD, Makumbi D, Magorokosho C, Nair S, Borem A, Ribaut JM, Banziger M, Prasanna BM, Crossa J, Babu R: QTL mapping in three tropical maize populations reveals a set of constitutive and adaptive genomic regions for drought tolerance. Theor Appl Genet. 2013, 126: 583-600. 10.1007/s00122-012-2003-7.
Cheng MC, Liao PM, Kuo WW, Lin TP: The Arabidopsis ETHYLENE-RESPONSE-FACTOR1 regulates abiotic-stress-responsive gene expression by binding to different cis-acting elements in response to different stress signals. Plant Physiol. 2013, 162: 1566-1582. 10.1104/pp.113.221911.
Garay-Arroyo A, Ortiz-Moreno E, de la Paz SM, Murphy AS, García-Ponce B, Marsch-Martínez N, de Folter S, Corvera-Poiré A, Jaimes-Miranda F, Pacheco-Escobedo MA, Dubrovsky JG, Pelaz S, Álvarez-Buylla ER: The MADS transcription factor XAL2/AGL14 modulates auxin transport during Arabidopsis root development by regulating PIN expression. EMBO J. 2013, 32: 2884-2895. 10.1038/emboj.2013.216.
Tapia-Lopez R, García-Ponce B, Dubrovsky JG, Garay-Arroyo A, Pérez-Ruiz RV, Kim SH, Acevedo F, Pélaz S, Alvarez-Buylla ER: An AGAMOUS-related MADS-box gene, XAL1 (AGL12), regulates root meristem cell proliferation and flowering transition in Arabidopsis. Plant Physiol. 2008, 146: 1182-1192. 10.1104/pp.107.108647.
Shen W, Wei Y, Dauk M, Tan Y, Taylor DC, Selvaraj G, Zou J: Involvement of a glycerol-3-phosphate dehydrogenase in modulating the NADH/NAD+ ratio provides evidence of a mitochondrial glycerol-3-phophate shuttle in Arabidopsis. Plant Cell. 2006, 18: 422-441. 10.1105/tpc.105.039750.
Hu J, Zhang Y, Wang J, Zhou Y: Glycerol affects root development through regulation of multiple pathways in Arabidopsis. PLoS One. 2014, 9: e86269-10.1371/journal.pone.0086269.
Li WX, Oono Y, Zhu J, He XJ, Wu JM, Iida K, Lu XY, Cui X, Jin H, Zhu JK: The Arabidopsis NFYA5 transcription factor is regulated transcriptionally and post transcriptionally to promote drought resistance. The Plant Cell. 2008, 20: 2238-2251. 10.1105/tpc.108.059444.
Mu J, Tan H, Hong S, Liang Y, Zuo J: Arabidopsis transcription factor genes NF-YA1, 5, 6, and 9 play redundant roles in male gametogenesis, embryogenesis, and seed development. Mol Plant. 2012, 6: 188-201.
Hartmann U, Sagasser M, Mehrtens F, Stracke R, Weisshaar B: Differential combinatorial interactions of cis-acting elements recognized by R2R3-MYB, BZIP, and BHLH factors control light-responsive and tissue-specific activation of phenylpropanoid biosynthesis genes. Plant Mol Biol. 2005, 57: 155-171. 10.1007/s11103-004-6910-0.
Brown DE, Rashotte AM, Murphy AS, Normanly J, Tague BW, Peer WA, Taiz L, Muday GK: Flavonoids act as negative regulators of auxin transport in vivo in Arabidopsis. Plant Physiol. 2001, 126: 524-535. 10.1104/pp.126.2.524.
Nepolean T, Hossain F, Shiriga K, Mittal S, Arora K, Rathore A, Mohan S, Shah T, Sharma R, Namratha PM, Amitha Mithra SV, Mohaptara T, Gupta HS: Unraveling the genetic architecture of subtropical maize (Zea mays L.) lines and their utility in breeding programs. BMC Genomics. 2013, 14: 877-10.1186/1471-2164-14-877.
Barton NH: Genetic hitchhiking. Philos Trans R Soc Lond B Biol Sci. 2000, 355: 1553-1562. 10.1098/rstb.2000.0716.
Davletova S, Schlauch K, Coutu J, Mittler R: The zinc-finger protein zat12 plays a central role in reactive oxygen and abiotic stress signaling in Arabidopsis. Plant Physiol. 2005, 139: 847-856. 10.1104/pp.105.068254.
Seo PJ, Xiang F, Qiao M, Park JY, Lee YN, Kim SG, Lee YH, Park WJ, Park CM: The MYB96 transcription factor mediates abscisic acid signaling during drought stress response in Arabidopsis. Plant Physiol. 2009, 151: 275-289. 10.1104/pp.109.144220.
Lu Y, Li Y, Zhang J, Xiao Y, Yue Y, Duan L, Zhang M, Li Z: Overexpression of Arabidopsis molybdenum cofactor sulfurase gene confers drought tolerance in maize (Zea mays L.). PLoS One. 2013, 8: e52126-10.1371/journal.pone.0052126.
Schroeder JI, Kwak JM, Allen GJ: Guard cell abscisic acid signalling and engineering drought hardiness in plants. Nature. 2001, 410: 327-330. 10.1038/35066500.
Laporte MM, Shen B, Tarczynski MC: Engineering for drought avoidance: expression of maize NADP-malic enzyme in tobacco results in altered stomatal function. J Exp Bot. 2002, 53: 699-705. 10.1093/jexbot/53.369.699.
Cruz de Carvalho MH: Drought stress and reactive oxygen species. Plant Signal Behav. 2008, 3: 156-165. 10.4161/psb.3.3.5536.
Van Breusegem F, Dat JF: Reactive oxygen species in plant cell death. Plant Physiol. 2006, 141: 384-390. 10.1104/pp.106.078295.
Miao Y, Lv D, Wang P, Wang XC, Chen J, Miao C, Song CP: An Arabidopsis glutathione peroxidase functions as both a redox transducer and a scavenger in abscisic acid and drought stress responses. The Plant Cell. 2006, 18: 2749-2766. 10.1105/tpc.106.044230.
Chen JH, Jiang HW, Hsieh EJ, Chen HY, Chien CT, Hsieh HL, Lin TP: Drought and salt stress tolerance of an Arabidopsis glutathione S-transferase U17 knockout mutant are attributed to the combined effect of glutathione and abscisic acid. Plant Physiol. 2012, 158: 340-351. 10.1104/pp.111.181875.
Zhu JK: Salt and drought signal transduction in plants. Annu Rev Plant Biol. 2002, 53: 247-273. 10.1146/annurev.arplant.53.091401.143329.
Hong Y, Zhang W, Wang X: Phospholipase D and phosphatidic acid signalling in plant response to drought and salinity. Plant Cell Environ. 2010, 33: 627-635. 10.1111/j.1365-3040.2009.02087.x.
Hund A, Ruta N, Liedgens M: Rooting depth and water use efficiency of tropical maize inbred lines, differing in drought tolerance. Plant Soil. 2009, 318: 311-325. 10.1007/s11104-008-9843-6.
Overvoorde P, Fukaki H, Beeckman T: Auxin control of root Development. Cold Spring Harb Perspect Biol. 2010, 2: a001537-
Pasapula V, Shen G, Kuppu S, Paez-Valencia J, Mendoza M, Hou P, Chen J, Qiu X, Zhu L, Zhang X, Auld D, Blumwald E, Zhang H, Gaxiola R, Payton P: Expression of an Arabidopsis vacuolar H + -pyrophosphatase gene (AVP1) in cotton improves drought- and salt tolerance and increases fibre yield in the field conditions. Plant Biotechnol J. 2011, 9: 88-99. 10.1111/j.1467-7652.2010.00535.x.
Li B, Wei A, Song C, Li N, Zhang J: Heterologous expression of the TsVP gene improves the drought resistance of maize. Plant Biotechnol J. 2008, 6: 146-159. 10.1111/j.1467-7652.2007.00301.x.
Ratcliffe OJ, Riechmann JL: Arabidopsis transcription factors and the regulation of flowering time: a genomic perspective. Curr Issues Mol Biol. 2002, 4: 77-91.
Chinnusamy V, Gong Z, Zhu JK: Abscisic acid-mediated epigenetic processes in plant development and stress responses. J Integr Plant Biol. 2008, 50: 1187-1195. 10.1111/j.1744-7909.2008.00727.x.
Wang WS, Pan YJ, Zhao XQ, Dwivedi D, Zhu LH, Ali J, Fu BY, Li ZK: Drought-induced site-specific DNA methylation and its association with drought tolerance in rice (Oryza sativa L.). J Exp Bot. 2010, 62: 1951-1960.
Cocciolone SM, Chopra S, Flint-Garcia SA, McMullen MD, Peterson T: Tissue-specific patterns of a maize Myb transcription factor are epigenetically regulated. The Plant J. 2001, 27: 467-478. 10.1046/j.1365-313X.2001.01124.x.
Colaneri AC, Jones AM: Genome-wide quantitative identification of DNA differentially methylated sites in Arabidopsis seedlings growing at different water potential. PLoS One. 2013, 8: e59878-10.1371/journal.pone.0059878.
Su Z, Ma X, Guo H, Sukiran NL, Guo B, Assmann SM, Ma H: Flower development under drought stress: morphological and transcriptomic analyses reveal acute responses and long-term acclimation in Arabidopsis. Plant Cell. 2013, 25: 3785-3807. 10.1105/tpc.113.115428.
Sewelam N, Kazan K, Thomas-Hall SR, Kidd BN, Manners JM, Schenk PM: Ethylene response factor 6 is a regulator of reactive oxygen species signaling in Arabidopsis. PLoS One. 2013, 8: e70289-10.1371/journal.pone.0070289.
Joo J, Choi HJ, Lee YH, Kim YK, Song SI: A transcriptional repressor of the ERF family confers drought tolerance to rice and regulates genes preferentially located on chromosome 11. Planta. 2013, 238: 155-170. 10.1007/s00425-013-1880-6.
Mishra G, Zhang W, Deng F, Zhao J, Wang X: A bifurcating pathway directs abscisic acid effects on stomatal closure and opening in Arabidopsis. Science. 2006, 312: 264-266. 10.1126/science.1123769.
Mane S, Vasquez-Robinet C, Sioson A, Heath L, Grene R: Early PLDα-mediated events in response to progressive drought stress in Arabidopsis: a transcriptome analysis. J Exp Bot. 2007, 58: 241-252.
Kam J, Gresshoff P, Shorter R, Xue GP: The Q-type C2H2 zinc finger subfamily of transcription factors in Triticum aestivum is predominantly expressed in roots and enriched with members containing an EAR repressor motif and responsive to drought stress. Plant Mol Biol. 2008, 67: 305-322. 10.1007/s11103-008-9319-3.
Sakamoto H, Maruyama K, Sakuma Y, Meshi T, Iwabuchi M, Shinozaki K, Yamaguchi-Shinozaki K: Arabidopsis Cys2/His2-type zinc-finger proteins function as transcription repressors under drought, cold, and high-salinity stress conditions. Plant Physiol. 2004, 136: 2734-2746. 10.1104/pp.104.046599.
Huang J, Wang JF, Wang QH, Zhang HS: Identification of a rice zinc finger protein whose expression is transiently induced by drought, cold but not by salinity and abscisic acid. DNA Seq. 2005, 16: 130-136.
Huang XY, Chao DY, Gao JP, Zhu MZ, Shi M, Lin HX: A previously unknown zinc finger protein, DST, regulates drought and salt tolerance in rice via stomatal aperture control. Genes Dev. 2009, 23: 1805-1817. 10.1101/gad.1812409.
Leigh RA: The Plant Vacuole. Volume 25. Edited by: Leigh RA, Sanders D. 1997, San Diego: Academic Press, 171-194.
Kodaira KS, Qin F, Tran LS, Maruyama K, Kidokoro S, Fujita Y, Shinozaki K, Yamaguchi-Shinozaki K: Arabidopsis Cys2/His2 zinc-finger proteins AZF1 and AZF2 negatively regulate abscisic acid-repressive and auxin-inducible genes under abiotic stress conditions. Plant Physiol. 2011, 157: 742-756. 10.1104/pp.111.182683.
Base SAS® 9.3 Procedures Guide. 2011, Cary: SAS Institute Inc.
Aulchenko YS, Ripke S, Isaacs A, van Duijn CM: GenABEL: an R library for genome-wide association analysis. Bioinformatics. 2007, 23: 1294-1296. 10.1093/bioinformatics/btm108.
Chen WM, Abecasis GR: Family-based association tests for genomewide association scans. Am J Hum Genet. 2007, 81: 913-926. 10.1086/521580.
Yu J, Pressoir G, Briggs WH, Vroh Bi I, Yamasaki M, Doebley JF, McMullen MD, Gaut BS, Nielsen DM, Holland JB, Kresovich S, Buckler ES: A unified mixed-model method for association mapping that accounts for multiple levels of relatedness. Nat Genet. 2006, 38: 203-208. 10.1038/ng1702.
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc B. 1995, 57: 289-300.
Lipka AE, Tian F, Wang Q, Peiffer J, Li M, Bradbury PJ, Gore M, Buckler ES, Zhang Z: GAPIT: genome association and prediction integrated tool. Bioinformatics. 2012, 28: 2397-2399. 10.1093/bioinformatics/bts444.
VanRaden PM: Efficient methods to compute genomic predictions. J Dairy Sci. 2008, 91: 4414-4423. 10.3168/jds.2007-0980.
We thank the breeders from different breeding centres in India and from CIMMYT for providing the seeds, and the National Agricultural Innovation Project (NAIP, Component IV) and the ICAR Network Project on Transgenics in Crop Plants (Maize Functional Genomics Component) for funding the study.
The authors declare that they have no competing interest.
NT and HSG conceived and designed the experiments; KS, SMt, SMn and PMN generated genotyping data; FH, SD, TSR, and SK generated phenotyping data; RS, KA, KS, SMt, SMn, AR, and TS analyzed the data; TM helped with SNP data generation; and NT and KA drafted the paper. HSG coordinated the research. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Table S1: Agronomic performance of genotypes in well-watered and water-stressed conditions. (XLSX 10 KB)
Additional file 2: Table S2: Phenotypic correlation coefficients for agronomic traits in well-watered and water-stressed conditions. Significant phenotypic correlations (p < 0.01) are marked with an asterisk. (XLSX 10 KB)
Additional file 3: Table S3: GWA signals observed for seven agronomic traits under water stress. Significant associations (P < 5.8 × 10−4 as the cut-off) for all agronomic traits were observed at all the three locations (New Delhi, Hyderabad, and Karimnagar) and for the mean of all phenotypic data and four models (M1, M2, M3, and M4) (Refer Methods for details). All FDR-corrected SNP associations are in bold. The drought-tolerant gene annotations (refer Figure 4) are marked with the respective references. (XLSX 48 KB)
Additional file 4: Figure S1: Manhattan plots showing SNPs significantly associated (log10P < −4.0) with seven agronomic traits under well-watered conditions measured across M2 (see Methods for details). (TIFF 643 KB)
Additional file 5: Figure S2: Manhattan plots showing SNPs significantly associated (log10P < −4.0) to seven agronomic traits under well-watered conditions measured across M1 (see Methods for details). (TIFF 634 KB)
Additional file 6: Figure S3: Manhattan plots showing SNPs significantly associated (log10P < −4.0) to seven agronomic traits under well-watered conditions measured across M3 (see Methods for details). (TIFF 635 KB)
Additional file 7: Figure S4: Manhattan plots showing SNPs significantly associated (log10P < −4.0) to seven agronomic traits under well-watered conditions measured across M4 (see Methods for details). (TIFF 638 KB)
Additional file 8: Figure S5: Manhattan plots showing SNPs significantly associated (log10P < −4.0) to seven agronomic traits under water stress measured across M1 (see Methods for details). (TIFF 634 KB)
Additional file 9: Figure S6: Manhattan plots showing SNPs significantly associated (log10P < −4.0) to seven agronomic traits under water stress measured across M3 (see Methods for details). (TIFF 633 KB)
Additional file 10: Figure S7: Manhattan plots showing SNPs significantly associated (log10P < −4.0) to seven agronomic traits under water stress measured across M4 (see Methods for details). (TIFF 626 KB)
Additional file 11: Figure S8: Manhattan plots showing SNPs significantly associated (log10P < −4.0) with seven agronomic traits under water stress measured across M2 (see Methods for details). (TIFF 628 KB)
About this article
Cite this article
Thirunavukkarasu, N., Hossain, F., Arora, K. et al. Functional mechanisms of drought tolerance in subtropical maize (Zea mays L.) identified using genome-wide association mapping. BMC Genomics 15, 1182 (2014). https://doi.org/10.1186/1471-2164-15-1182
- Genome-wide SNPs
- Association mapping
- Functional mechanisms
- Candidate SNPs
- Water stress
- Drought tolerance