Next Article in Journal
Genomic Predictions for Common Bunt, FHB, Stripe Rust, Leaf Rust, and Leaf Spotting Resistance in Spring Wheat
Next Article in Special Issue
Genetic Architecture of Maize Stalk Diameter and Rind Penetrometer Resistance in a Recombinant Inbred Line Population
Previous Article in Journal
Suicide and Changes in Expression of Neuronal miRNA Predicted by an Algorithm Search through miRNA Databases
Previous Article in Special Issue
Genome-Wide Profiling of Alternative Splicing and Gene Fusion during Rice Black-Streaked Dwarf Virus Stress in Maize (Zea mays L.)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dissection of Maize Drought Tolerance at the Flowering Stage Using Genome-Wide Association Studies

1
Crop Functional Genome Research Center, Biotechnology Research Institute, Chinese Academy of Agricultural Sciences, Beijing 100081, China
2
National Nanfan Research Institute (Sanya), Chinese Academy of Agricultural Sciences, Sanya 572024, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Genes 2022, 13(4), 564; https://doi.org/10.3390/genes13040564
Submission received: 14 March 2022 / Revised: 19 March 2022 / Accepted: 21 March 2022 / Published: 23 March 2022
(This article belongs to the Special Issue Maize Functional Genomics, Genetics and Breeding)

Abstract

:
Drought is one of the most critical environmental factors constraining maize production. When it occurs at the flowering stage, serious yield losses are caused, and often, the damage is irretrievable. In this study, anthesis to silk interval (ASI), plant height (PH), and ear biomass at the silking date (EBM) of 279 inbred lines were studied under both water-stress (WS) and well-water (WW) field conditions, for three consecutive years. Averagely, ASI was extended by 25.96%, EBM was decreased by 17.54%, and the PH was reduced by 12.47% under drought stress. Genome-wide association studies were carried out using phenotypic values under WS, WW, and drought-tolerance index (WS-WW or WS/WW) and applying a mixed linear model that controls both population structure and relative kinship. In total, 71, 159, and 21 SNPs, located in 32, 59, and 12 genes, were significantly (P < 10−5) associated with ASI, EBM, and PH, respectively. Only a few overlapped candidate genes were found to be associated with the same drought-related traits under different environments, for example, ARABIDILLO 1, glycoprotein, Tic22-like, and zinc-finger family protein for ASI; 26S proteasome non-ATPase and pyridoxal phosphate transferase for EBM; 11-ß-hydroxysteroid dehydrogenase, uncharacterised, Leu-rich repeat protein kinase, and SF16 protein for PH. Furthermore, most candidate genes were revealed to be drought-responsive in an association panel. Meanwhile, the favourable alleles/key variations were identified with a haplotype analysis. These candidate genes and their key variations provide insight into the genetic basis of drought tolerance, especially for the female inflorescence, and will facilitate drought-tolerant maize breeding.

1. Introduction

Maize (Zea mays L.) is the most widely grown crop worldwide and has incredible importance for food, feed, and other industrial products [1]. It was estimated that maize production needs to be boosted by 67% to cope with the increased population growth and food demand in 2050 [2]. Thus, it is of great significance to increase maize yield productivity and reduce yield loss caused by biotic and abiotic stresses.
Drought is considered one of the most detrimental curbs of agriculture, which endangers maize production globally because of its erratic nature [3,4]. Maize is an open-pollinating crop that is extremely sensitive to drought stress throughout its life span, especially one week before and three weeks after the flowering period, causing severe yield loss by 15–25% [5]. At this stage, drought inhibits ear development, causing abnormal differentiation in spikelets, and changes in hormone signalling involved in cell division, growth, and primordium development [6], leading to asynchronous development in tassel and ear, a tremendous extension of anthesis and silking interval (ASI), and reduction in silk receptivity [7,8]. Consequently, this time slack between pollen release and silk emergence adversely affects pollination and kernel set and results in the reduction in grain yield [9,10]. Drought during the grain filling stage resulted in starch quality and quantity reduction, leading to a great loss of grain weight and grain quality [11]. Drought stress also reduces plant and ear height and decreases the availability of photosynthate for grain production, ultimately leading to yield reduction [12,13]. Due to global warming and increasing water resources shortage, the water deficit has become more and more prominent. Breeding maize for drought tolerance, especially at the flowering stage, is thus of significant importance for global food production.
Association analysis based on linkage disequilibrium (LD) acquired increasing popularity and accuracy in the genetic architecture of polygenic traits in crops due to broader genetic variation, larger numbers of alleles, and the maximum number of recombinants obtained [14,15]. Genome-wide association study (GWAS), taking advantage of natural variation and historical recombination, have emerged as alternative tools to linkage mapping for identifying superior alleles for complex traits, with reduced time consumption and increased mapping resolution [16,17]. GWAS have been employed to identify numerous SNPs directly associated with drought tolerance, especially at the seedling stage. ZmVPP1 was associated with dehydration tolerance in maize seedlings, encoding a vacuolar-type H+ pyrophosphatase, which improved seedling drought tolerance in maize due to enhanced photosynthetic activity and root development [18]. Significant associations were detected between maize seedling drought tolerance and functional genes promoter variation, such as SNPs located in the promoter of ZmDREB2.7 and miniature inverted-repeat transposable element in ZmNAC111 promoter, which determined gene expression for dehydration tolerance in maize seedlings [19,20]. A 368 maize association panel was used to conduct an association analysis on ZmPP2C-A family genes and found that ZmPP2C-A-10 was closely related to drought stress, through regulating the ABA signalling pathway [21]. This study reveals the correlation between endoplasmic reticulum stress response and drought resistance.
Research on maize drought-tolerant genes for the flowering period is scarce and relatively lagging, which is likely because maize drought tolerance is a set of traits with trade-offs, and plasticity and the internal connections are extremely complicated. However, few genes contributing to drought tolerance in maize were identified through GWAS or linkage mapping using a larger-scale population due to uncontrollable field conditions. NAC transcription factor NUT1 was found specifically expressed in protoxylem at the flowering stage, and functions to manipulate water transport by maintaining protoxylem vessel integrity through activating genes necessary for secondary cell wall reinforcement, thus affecting drought tolerance in NUT1 mutants [22]. A recent study found that ZmEXPA4, identified through transcriptomics analysis, functions during ear growth and silk elongation, alleviating drought-induced ASI elongation without affecting other agronomic traits [23]. More recently, using a high-throughput phenotyping platform, 368 maize inbred lines were continuously non-destructive, tested under normal watering and drought stress conditions at multiple growth stages, resulting in 2318 candidate genes associated with i-traits and drought tolerance in maize [6]. Mutant-based functional validation has shown that ZmcPGM2 (involved in sugar metabolism) and ZmFAB1A (involved in phosphoinositide metabolism) can negatively regulate drought resistance in maize during flowering.
In this study, an association population consisting of 279 inbred lines was evaluated under field WW and WS conditions. We aimed to dissect the drought-tolerant candidate genes at the flowering period through GWAS and to provide potential genetic variations for breeding drought-tolerant maize varieties.

2. Materials and Methods

2.1. Plant Materials and Experimental Design

An association panel for GWAS comprising 279 inbred lines with diverse genetic variation (Table S1) was selected from comprehensive geographical regions including tropical, subtropical, and temperate areas of the world, as described by Zhang et al. [24]. This panel was obtained from an available association population for eQTL mapping on maize kernel development [25,26].
The field trials were conducted in Urumqi (Xinjiang province: 43°54′ N, 87°28′ E), from May to September. The flowering date of the association panel was recorded in 2016 under normal growing conditions and then divided into three patches according to growing degree day (GDD) value at −15 days before pollen shedding (−15 D). The daily temperature was recorded, and GDD was calculated for each line using the following formula:
G D D = ( L + H ) 5 50
where L means daily lowest temperature (°F), and H means daily highest temperature (°F).
In order to obtain synchronised flowering, 279 lines were delayed sowing in patches. Each patch was in randomised complete block design, with three independent repeats during the years of 2017, 2018, and 2020. A total of 119 lines belonged to the 1st patch (No. 1–119, Table S1), with GDD ranging from 1612.1 to 1901.9 for the −15 D, and were planted on 8 May 2017. A total of 138 lines were planted 10 days later than those of the first patch (18 May 2017), with GDD values at −15 D ranging from 1324.7 to 1599.3. The remaining 22 lines (No. 258–279, Table S1) with earlier flowering times (GDD = 962–1291.6) were planted on 28 May 2017. Each row was 3.6 m in length, with planting space of 0.24 m and 1.1 m width across rows. Every line was planted side by side, one row for water stress treatment (WS) and the adjacent for a well-watered regime (WW, Figure 1A), with an independent valve to control the irrigation of each row. Drought stress was applied in the WS regime under the Protocol of Irrigation Management for Maize Drought Trials of Syngenta. Briefly, daily moisture level was recorded using the soil moisture sensors WATERMARK MONITOR 900 M (IRROMETER) installed in the field (Figure 1B) and connected to the sensor in every patch at 60 cm and 90 cm depth in multiple locations according to instructions from the trial sponsor in both WS and WW regimes. For the intended drought stress management at the flowering stage, a drought stress trial was conducted by withdrawing irrigation initiated at the −21 D before anthesis estimated according to GDD. The soil moisture was maintained at <40 centibars before −15 D and at 80–120 centibars from −15 D to −7 D. At −7 D, the soil moisture reached 120–150 centibars, which is the trigger point for drought stress, and was kept at 150–200 centibars until 14 days after anthesis when irrigation was resumed. The fertilisers, herbicides, and insecticides were applied as per requirements according to the local recommendation practices. In the WW regime, normal irrigation was provided once a week.

2.2. Phenotyping for Drought-Stress-Related Traits

PH and flowering time, such as anthesis date and silking date, were examined in inbred lines under WS and WW conditions. Days to anthesis (AD) and days to silking (SD) were determined by the number of days from planting to pollen shedding of 50% of plants and 50% of plants having clearly visible silks, respectively. ASI was calculated as the interval between anthesis date and silking date (ASI = SD − AD) and counted in days [9]. PH was measured from the soil level to the lowest tassel branch of each plant and recorded in centimetres (cm). For ear biomass, the ear at its silking date (Figure 1C) of each line was harvested and dried to a consistent weight at 72 °C. The ear biomasses were weighed with an electronic weighing balance and recorded in grams (g). The drought-tolerant index was calculated by dividing the mean values of studied traits in the WS regime by those of the WW regime. To estimate the random errors, each measurement per inbred line comprised 6–8 individual plants with 3 independent repeats.

2.3. Association Analysis for ASI, EBM, and PH

For association analysis, the 279 lines were genotyped by 776,254 high-density SNP markers, with MAF > 0.05, which was excavated previously [26]. The genome-wide association analysis was conducted with a mixed linear model (MLM) embedded in TASSEL V5.0 [27]. Population structure (Q) and kinship (K) were estimated according to Pang et al. [26]. Briefly, kinship matrices and principal components were estimated based on 236,205 SNPs with MAF > 0.05. The top 3 principal components were selected as population structure, and for kinship estimation, the ‘Normalised_IBS’ method was used. The regression-based coefficient of determination values of all significantly associated SNPs was recorded to determine the variations explained by each SNP locus. Drought tolerance was controlled by several genes with minor effects and non-independence of SNPs in maize genome leading to strong LD. Therefore, a threshold of P < 10−5 was applied in this study to eliminate false negatives. Multiple Manhattan plots were drawn in R (https://www.r-project.org/ (accessed on 8 November 2021)).
Based on the maize B73 reference genome (V4), genes directly colocalised with significantly associated SNPs considered as candidates. Alternatively, genes located within a corresponding LD interval (r2 ≥ 0.2) were also considered [18]. The gene annotation and function were retrieved from the Maize Genomics database (http://www.maizegdb.org/ (accessed on 6 December 2020)), Gramene database (http://gramene.org/ (accessed on 6 December 2020)), NCBI (https://www.ncbi.nlm.nih.gov/ (accessed on 6 December 2020)), and other available sources in the literature.

2.4. Drought Responsive, Linkage Disequilibrium, and Haplotype Analysis of Candidate Genes

In total, 197 inbred lines from a previous association mapping population were planted in pots under well-watered and drought stress treatments [6]. Drought stress conditions were relative to WW conditions, with the soil moisture levels of WS conditions dropping from 50% to 10% and soil moisture levels of WW conditions keeping values around 50% at the flowering stage. Candidate genes—Zm00001d013992, Zm00001d020506, Zm00001d029937, Zm00001d029938, and Zm00001d039319—were commonly identified under two environments. Therefore, to determine whether they were drought-responsive, their expression patterns were analysed according to the RNA profiling data of the abovementioned 197 lines (unpublished data from Dr. Mingqiu Dai, Huazhong Agricultural University). Haplotype analysis normally helps to understand the causal variation in a gene identified by GWAS [20,26], which in turn benefits future functional marker development. To confirm whether the SNPs identified for each gene were meaningful, haplotype analysis was conducted to check the allele effects of the most significant SNPs for the six overlapping candidate genes. A two-pair t-test was used to analyse their allelic effect in corresponding phenotypic performance.

2.5. Statistical Analysis

The phenotypic data were filtered by removing the suspicious value of each replicate by Q-test in excel [24]. The mean value of three replicates in each environment was used for association analysis. Best linear unbiased estimators (BLUEs) were estimated using the genotype and covariate as a fixed factor, and the rest as random factors. Analysis of variance (ANOVA), correlations analysis, and broad-sense heritability (H2) were estimated by SPSS 25.0. The H2 was estimated according to the following equation:
H 2 = σ G 2 σ G 2 + σ G E 2 / n + σ e 2 / n r
where σ2G represents genetic variance, σ2GE shows the genetic and environmental interaction variation, σ2e shows residual error variance, n shows the number of environments, and r shows the number of replicates [28]. LD analysis of these candidate genes was conducted by SNPs within a gene using the ‘LD heatmap’ in the R software (https://www.r-project.org/ (accessed on 8 November 2021)).

3. Results

3.1. Performance of Drought-Tolerant Phenotypes in the Association Panel

The test of normality revealed that the frequency distributions for all of the studied traits were near normal for most of the traits in the association panel (Figure 2; Table 1). The existence of variation might be caused by fluctuating environmental conditions in the field or genotypic differences. The descriptive statistics, heritability analysis, and coefficient of variation (CV) for the phenotypic traits are listed in Table 1. Generally, drought stress significantly decreased the ear biomass and PH (Figure 1B,C) while enlarging the ASI, suggesting that water stress at the flowering stage had diverse effects on drought-related traits.
The mean values of ASI under drought stress were 5.58 days with the CV 42.97%, while under WW condition, ASI was 4.43 days, with a CV value of 48.89% (Table 1). The average ASI across the three years under the WS regime was 1.15 days larger than that of the WW regime. Under drought stress, ASI was increased by 25.96%, indicating that drought stress significantly enlarged silk extrusion at the population level. Drought stress caused a 17.54% reduction in EBM with a minimum of 0.48 g and a maximum of 3.08 g under WS conditions, while they were 0.48 g and 3.81 g under the well-watered regime. The mean value of ear biomass was 1.71 g under the WW condition, with a CV of 27.20%, while it was 1.41 g under the WS condition, with a CV of 29.15% for the three-year environments, respectively (Table 1), suggesting a suppressed ear development by drought stress. Likewise, drought stress caused a 12.47% reduction in PH; the minimum PH was 43.21 cm, and the maximum was 169.72 cm under the WS condition, while they were 51.23 cm and 183.74 cm under the WW condition. The mean value of PH under the WS regime was 111.36 cm with a CV of 19.34%, significantly lower than that in the WW regime, which was 127.22 cm with a CV of 16.92% for the three years, respectively (Table 1). H2 values of ASI, EBM, and PH were higher than 80% under the two water treatments across the three years (Table 1), indicating that these traits are highly heritable in the two water treatment conditions.

3.2. Correlations among Drought-Related Traits

The correlations among ASI, EBM, and PH values under WS and WW conditions are listed in Table 2. Significantly positive correlations (P < 0.05) were found between ASI and PH under WS (0.13) and WW (0.15) conditions. Additionally, significant correlations were found between ear biomass and PH under both WS (0.22) and WW (0.18) conditions. A weak and negative correlation was found between ASI and ear biomass.

3.3. GWAS for Maize Drought Tolerance Genes

A total of 71, 159, and 21 SNPs were significantly associated with ASI, EBM, and PH located in 36, 81, and 16 genes, respectively. These genes were scattered over 10 chromosomes, with R2 ranging from 7.19% to 15.52%. The number of identified SNPs in each chromosome ranged from 9 to 31 on chromosomes 9 and 1, respectively. The information about the SNPs and candidate genes is shown in Tables S2–S4. A total of 17 and 48 SNPs located in 13 and 18 genes were identified for ASI-WS and ASI-WW, with R2 ranging from 7.73% to 14.32% (Figure 3; Table 3 and Table S2). As to the ASI delay, only five SNPs in four genes were identified in 2017 and 2018 (Table S2). In total, 49 and 93 SNPs located in 23 and 43 were associated with EBM-WS and EBM-WW, with R2 ranging from 7.99% to 12.52% (Figure 4; Table 3 and Table S3). In addition, 17 SNPs in 15 genes were associated with the EBM drought tolerance index (WS/WW) across the three years (Table S3). For both conditions, nine SNPs were associated with PH-WS and PH-WW and located in seven and six genes, with R2 ranging from 7.19% to 15.52% across the three-year environments (Figure 5; Table 3 and Table S4). Three genes in Chr. 2, 7, and 9 were identified for the PH drought tolerance index (Table S4).

3.4. Common Genes Identified for Ear Development across Multiple Years or Conditions

Several overlapping genes were identified among years and water treatments (Table 3). Four candidate genes encoding ARABIDILLO 1 protein (Zm00001d029938) and Glycoprotein (Zm00001d029937) on Chr. 1, Tic22-like family protein (Zm00001d039319), and zinc-finger family protein (Zm00001d042997) on Chr. 3, were identified for ASI under both WS and WW conditions in 2018 (Figure 3), suggesting that these genes might be promising candidates, as they can function in maize inflorescence development with or without water deficit. An SNP (S5_27121944), located in pyridoxal phosphate transferase encoding gene (Zm00001d013992), was significantly associated (P < 10−5) with EBM under drought stress condition on Chr. 5 and was consistently detected for two years environments of 2018 and 2020 (Figure 4), suggesting the drought tolerance role of this gene in developing ear. A candidate gene (Zm00001d020506) encoding 26S proteasome non-ATPase regulatory subunit-9 was detected for EBM on Chr. 7 under both WS and WW regimes during the 2017 field trial (Figure 4), suggesting its possible role in maize ear development.

3.5. Candidate Genes Drought Responsive Pattern

Most quantitative traits’ functional genes are responsive at the transcriptional level. Drought tolerance is complex and regulated by many quantitative trait loci (QTLs) with minor effects. In order to determine whether these candidate genes were drought-regulated, we analysed their expression level using the expression data from 197 diverse inbred lines under WS and WW conditions. Unsurprisingly, significant differential expression existed in most candidate genes between WW and WS treatments. Zm00001d013992, Zm00001d029938, and Zm00001d039319 were increased by 45.44%, 17.46%, and 6.01%, while Zm00001d029937 was reduced by 30.26%, compared with their expression under the WW condition (Figure 6). However, no obvious difference was observed for Zm00001d039319.

3.6. Allele Effects of Common Candidate Genes

SNPs within a gene might be critical to its function. Therefore, it was of great interest to find the true association between SNP variation and target traits in a population. Haplotype analysis is a useful strategy to extract more rare causal variants. In our study, inbred lines were grouped into HapA (A) and HapB (G) based on the most significant SNP (S5_27121944) variation in Zm00001d013992 (Table 3, Figure 7A), 226 (2018), and 213 (2020). The rare HapB group had significantly higher EBM (2.10 g) under WS conditions in both 2018 and 2020 (Figure 7B, Table S5) than HapA. Five polymorphisms—namely, chr7.S_116288756, chr7.S_116288791, chr7.S_116288792, chr7.S_116285652, and chr7.S_116285655—were identified in Zm00001d020506 (Figure 7C, Table S6). In total, 101 lines belonged to the HapA (AACCT) with average EBM values of 1.58 g and 1.79 g under WS and WW in 2017, while the other 107 lines carried the GTTTC haplotype with significantly lower EBM values of 1.35 g and 1.55 g under WS and WW (Figure 7D). Three SNPs in Zm00001d029937 were significantly associated with ASI under the WW condition in 2018 (Table 3, Figure 7E), and accordingly, two haplotype AAC (HapA) and GGG (HapB) were obtained. Of the 193 lines harbouring rare haplotype HapA, 15 showed longer ASI (>10 days) than those harbouring HapB, whose ASI values were 5.6 days and 6.6 days under WW and WS conditions, respectively (Figure 7F, Table S7). The same 15 lines belonging to the HapA group (TCGATAATC) presented longer ASI under both WW and WS conditions (Figure 7G,H, Table S8), which was a similar case as that of the adjacent gene Zm00001d029938. PZE-103003226 was the only SNP identified in Zm00001d039319 (Figure 7I). The 165 lines with shorter ASI, grouped into HapA (G), showed weaker sensitivity to drought stress in both WW and WS conditions (Figure 7J, Table S9). For Zm00001d042997, 9 SNPs classified the 237 inbred lines into 3 groups (Figure 7K, Table S10), among which 17 lines belonged to lower frequency of HapC (TGACTTAA) and had a larger ASI under both WS and WW conditions, compared with those belonging to HapA (AACTCCCG) and HapB (AACTCTAA), while no statistical difference was found between HapA and HapB (Figure 7L).

4. Discussion

Maize is highly sensitive to drought stress at the flowering stage. In this study, the silk extrusion was markedly delayed, and plant architecture was greatly affected under drought stress. These results are in line with previous reports that stated significant extension in ASI and a notable reduction in plant and ear height under WS conditions [29,30]. Water deficit at the flowering stage delays or inhibits plant growth and ear development, reducing ear biomass from 1.71 g to 1.41 g (Table 1). It was found that osmotic stress limits the dry matter accumulation by approximately 50% during serious water shortages [31,32]. Thus, considering the higher estimated heritability and vulnerability to drought stress, ear biomass could be an optimised option for improving maize selection under water-scarce conditions [33,34]. Plant height was significantly correlated with yield-related traits under both WW and WS conditions [35]. The average PH decreased less when the plant was exposed to drought, which ensures sufficient ‘source’ availability and exhibits better drought resistance. PH was positively correlated with EBM, suggesting the meaningful role of EBM in drought tolerance. A weak and negative correlation existed between the ASI and ear biomass, indicating that delayed silk extrusion has likely no correlation with ear development. Therefore, the three traits may not be tightly correlated with each other, or the correlation may be disturbed due to the variation of the fluctuating environmental conditions in the field.
Up to now, only few GWAS have been conducted on maize drought tolerance under complex field conditions [30,35]. Researchers tested a number of maize-nested association mapping populations under two contrasting water regimes for seven drought-related traits, including ASI-, PH-, and yield-related traits. Hundreds of promising QTLs and candidate genes were obtained through GWAS and linkage mapping [35]. In addition, many other candidate genes were detected to be associated with drought-tolerance-correlated yield and agronomic traits [30,36,37,38]; however, none of them were found in this study, which may be due to the various growth conditions and drought treatments. Initially, we aimed to find candidate genes contributing to drought tolerance, especially those for the drought tolerance index. However, only a few genes were identified, and no overlapping genes were identified across different environments, which might be due to the inherence of drought tolerance controlled by multiple minor effect genes. Fortunately, few candidate genes identified in this study were colocalised with reported QTLs. The uncharacterised gene GRMZM2G173084, associated with ASI-WW in 2017, overlapped in a QTL for both ASI-WW and ASI-WS, which was detected by joint linkage analysis in a CN-NAM population [35]. Zm00001d003939, encoding 11-ß-hydroxysteroid dehydrogenase 2, is the candidate gene for PH-WS-2017 located in a consistent QTL for PH-WW, which was identified in both CN-NAM and US-NAM populations. In addition, several overlapping genes were identified under water treatments among different years, while only Zm00001d013992 was commonly identified in a two-year environment. This result implied that maize drought tolerance is a complex trait, highly affected by environment and treatment. Hu et al. [39] reported that grain yield in WS (qWS-GY7-1) and an ear setting percentage in WW (qWW-ESP7-1), located in Chr. 7: 132.2–135.6 Mb, and seven other QTLs and one mQTL for drought-related traits also clustered in this region [40,41]. Coincidentally, GRMZM2G173084 encoding an uncharacterised protein was associated with ASI-WW in 2017. Additionally, Hu et al. [39] found that the drought tolerance allele of Chr. 7: 132.2–135.6 Mb improved GY in both WW and WS regions, suggesting that both regions influenced GY performance under water-limited conditions. GRMZM2G173084 colocalised in this interval, which indicated its promising role in drought tolerance. Zm00001d044411 for EBM-WS-18 fell in a QTL hotspot in Chr. 3: 219.8–223.7 Mb for GY, ESP, and ASI [39], and its extended region harbours several QTLs responsible for stay-green, leaf senescence, and chlorophyll content identified under normal growth and water stress conditions in maize [29,40,42,43]. Zm00001d032084 (chaperone protein dnaJ) for ASI-Delay-18 and Zm00001d026286 (40S ribosomal protein S11) for ASI-WS-18 colocalised in QTL clusters in bin 1.07 and 10.05–10.07, respectively, each harbouring at least three maize-flowering-time QTLs under different planting densities [44]. Zm00001d024783 (BHLH transcriptional factor 117) for PH-WW-18 and Zm00001d003939 (11-ß-hydroxysteroid dehydrogenase) for PH-WW-17 were localised in QTL (bin 10.07 and bin 2.05–2.07) for plant height of forage maize [45]. These results indicated that these colocalised genes might play important roles in maize development under different conditions, though further functional validations are needed.
Members of zinc-finger family protein play critical roles in plant growth and developmental processes, including flowering, senescence, and also abiotic stress responses [46,47]. A C2H2 zinc-finger transcription factor determines stomatal closure by modulating H2O2-homeostasis related genes, for example, peroxidases, glutathione S-transferase, and cytochrome P450s [48]. Hence, drought response in rice is regulated. Zm00001d042997 was identified encoding a HIT-type zinc-finger family protein and was associated with ASI-WW and ASI-WS, indicating a potential conserved abiotic-stress-tolerant role of zinc-finger family protein. It was stated that the F-box domain-containing protein ARABIDILLO-1 is conserved in plants, involved in root architecture development, and functions during rice abiotic stress, mainly through regulating root branching and lateral root development [49,50]. It was reported that ARABIDILLO-1-mediated protein degradation, most likely through modulating the GA3 signalling pathway [49]. However, another study revealed that plants with ARABIDILLO-1 knockout and overexpression responded normally to auxin and abscisic acid [51]. Zm00001d029938, encoding ARABIDILLO-1 in maize, was associated with ASI in this study. However, whether it functions in maize drought tolerance was not clearly determined. For EBM, an SNP (S5_27121944), which is annotated as pyridoxal phosphate-dependent transferase (Zm00001d013992), was consistently associated for consecutive two years under drought regime (Table 3) and was upregulated by drought stress. The lead SNP S5_27121944 (A/G) in Zm00001d013992 separated the association panel into two groups, and only around 6% lines (including BY855 and BY4960, etc) carried the favourable haplotype, exhibiting higher ear biomass, 3 g under WS (Figure 7A). Pyridoxal phosphate (PLP) is an active form of pyridoxine (vitamin B6), which functions as a coenzyme in several reactions such as decarboxylation, deamination, and transamination. The PLP dependent enzymes mainly perform in amino acid biosynthesis and the metabolism of its derived metabolites. Therefore, it is interesting to speculate that Zm00001d013992 might be involved in amino acid metabolism and promotes ear development in maize under drought stress conditions. Potential causal SNPs of candidate drought-tolerant genes could be used for drought-tolerant maize improvement through both genome selection and genome editing.

5. Conclusions

The findings of this study provide insights into the genetic basis of drought tolerance at the flowering stage, especially for the female inflorescence’s development. The overlapping genes are proposed as candidate genes for drought tolerance in maize. Moreover, those lines carrying favourable alleles could be used for drought-tolerance marker development, which is of benefit for future marker-assisted or genome-wide selection for drought-tolerant maize breeding. Future investigation is needed to explore the candidate gene function using CRISPR/Cas9 mediated genome editing and the underlying molecular mechanism of maize ear and silk development under water deficit condition.

Supplementary Materials

The following information is available online at https://www.mdpi.com/article/10.3390/genes13040564/s1, Table S1: List of inbred lines used for GWAS analysis, Table S2: Annotation of SNPs associated with ASI under two water regimes, Table S3: Annotation of SNPs associated with ear biomass at silking date (EBM) under two water regimes, Table S4: Annotation of SNPs associated with plant height under two water regimes, Table S5: Haplotype analysis of Zm00001d013992 for EBM in 2018 and 2020, Table S6: Haplotype analysis of Zm00001d020506 for EBM in 2017, Table S7: Haplotype analysis of Zm00001d029937 for ASI in 2018, Table S8: Haplotype analysis of Zm00001d029938 for ASI in 2018, Table S9: Haplotype analysis of Zm00001d039319 for ASI in 2018, Table S10: Haplotype analysis for Zm00001d042997 for ASI in 2018.

Author Contributions

Conceptualisation, J.Z. and S.U.K.; methodology, S.U.K., J.Z. and P.L.; formal analysis, S.U.K., Y.Z. and P.L.; investigation, S.U.K., Z.C., G.Z. and P.L.; resources, J.Z.; data curation, S.U.K., X.Z. and N.Z.; writing—original draft preparation, S.U.K. and Y.Z.; writing—review and editing, J.Z. and P.L.; supervision, J.Z. and P.L.; project administration, J.Z. and P.L.; funding acquisition, J.Z. and P.L. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the National Natural Science Foundation of China (Grant No. 32072067, 31801378), the 2020 Research Program of Sanya Yazhou Bay Science and Technology City (Grant No. SKJC-2020-02-005), and the Joint Research Program of Long Ping High-Tech with BRI-CAAS.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

GWAS results are provided in supplementary data, and the genotypic data of the association panel are available on request from the corresponding author.

Acknowledgments

We would like to thank Mingqiu Dai and Xiaopeng Sun from Huazhong Agricultural University, for kindly offering RNA-seq data of 197 inbred lines.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Shiferaw, B.; Prasanna, B.M.; Hellin, J.; Bänziger, M. Crops that feed the world 6. Past successes and future challenges to the role played by maize in global food security. Food Secur. 2011, 3, 307–327. [Google Scholar] [CrossRef] [Green Version]
  2. Ray, D.K.; Mueller, N.D.; West, P.C.; Foley, J.A. Yield trends are insufficient to double global crop production by 2050. PLoS ONE 2013, 8, e66428. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Daryanto, S.; Wang, L.; Jacinthe, P.A. Global synthesis of drought effects on maize and wheat production. PLoS ONE 2016, 11, e0156362. [Google Scholar] [CrossRef] [PubMed]
  4. Gupta, A.; Rico-Medina, A.; Caño-Delgado, A.I. The physiology of plant responses to drought. Science 2020, 368, 266–269. [Google Scholar] [CrossRef]
  5. Bänziger, M.; Edmeades, G.O.; Beck, D.; Bellon, M. Breeding for Drought and Nitrogen Stress Tolerance in Maize: From Theory to Practice; CIMMYT: Mexico, DF, USA, 2000; pp. 1–68. [Google Scholar]
  6. Wu, X.; Feng, H.; Wu, D.; Yan, S.; Zhang, P.; Wang, W.; Zhang, J.; Ye, J.; Dai, G.; Fan, Y.; et al. Using high-throughput multiple optical phenotyping to decipher the genetic architecture of maize drought tolerance. Genome Biol. 2021, 22, 1–26. [Google Scholar] [CrossRef]
  7. NeSmith, D.S.; Ritchie, J.T. Effects of soil water-deficits during tassel emergence on development and yield component of maize (Zea mays). Field Crops Res. 1992, 28, 251–256. [Google Scholar] [CrossRef]
  8. Saini, H.S.; Westgate, M.E. Reproductive development in grain crops during drought. Adv. Agron. 1999, 68, 59–96. [Google Scholar]
  9. Bolaños, J.; Edmeades, G.O. The importance of the anthesis-silking interval in breeding for drought tolerance in tropical maize. Field Crops Res. 1996, 48, 65–80. [Google Scholar] [CrossRef]
  10. Bruce, W.B.; Edmeades, G.O.; Barker, T.C. Molecular and physiological approaches to maize improvement for drought tolerance. J. Exp. Bot. 2002, 53, 13–25. [Google Scholar] [CrossRef]
  11. Barutçular, C.; Dizlek, H.; El-Sabagh, A.; Sahin, T.; EL-Sabagh, M.; Islam, M.S. Nutritional quality of maize in response to drought stress during grain-filling stages in mediterranean climate condition. J. Exp. Biol. Agric. Sci. 2016, 4, 644–652. [Google Scholar] [CrossRef]
  12. Lopes, M.S.; Araus, J.L.; Van Heerden, P.D.R.; Foyer, C.H. Enhancing drought tolerance in C4 crops. J. Exp. Bot. 2011, 62, 3135–3153. [Google Scholar] [CrossRef] [PubMed]
  13. Sari-Gorla, M.; Krajewski, P.; Di Fonzo, N.; Villa, M.; Frova, C. Genetic analysis of drought tolerance in maize by molecular markers. II. plant height and flowering. Theor. Appl. Genet. 1999, 99, 289–295. [Google Scholar] [CrossRef]
  14. Mackay, T.F.C. Quantitative trait loci in drosophila. Nat. Rev. Genet. 2001, 2, 11–20. [Google Scholar] [CrossRef] [PubMed]
  15. Yu, J.; Buckler, E.S. Genetic association mapping and genome organization of maize. Curr. Opin. Biotechnol. 2006, 17, 155–160. [Google Scholar] [CrossRef]
  16. Rosenberg, N.A.; Huang, L.; Jewett, E.M.; Szpiech, Z.A.; Jankovic, I.; Boehnke, M. Genome-wide association studies in diverse populations. Nat. Rev. Genet. 2010, 11, 356–366. [Google Scholar] [CrossRef] [Green Version]
  17. Yan, J.; Warburton, M.; Crouch, J. Association mapping for enhancing maize (Zea mays L.) genetic improvement. Crop Sci. 2011, 51, 433–449. [Google Scholar] [CrossRef]
  18. Wang, X.; Wang, H.; Liu, S.; Ferjani, A.; Li, J.; Yan, J.; Yang, X.; Qin, F. Genetic variation in ZmVPP1 contributes to drought tolerance in maize seedlings. Nat. Genet. 2016, 48, 1233–1241. [Google Scholar] [CrossRef]
  19. Liu, S.; Wang, X.; Wang, H.; Xin, H.; Yang, X.; Yan, J.; Li, J.; Tran, L.S.; Shinozaki, K.; Yamaguchi-Shinozaki, K.; et al. Genome-wide analysis of ZmDREB genes and their association with natural variation in drought tolerance at seedling stage of Zea mays L. PLoS Genet. 2013, 9, e1003790. [Google Scholar] [CrossRef] [Green Version]
  20. Mao, H.; Wang, H.; Liu, S.; Li, Z.; Yang, X.; Yan, J.; Li, J.; Tran, L.S.; Qin, F. A transposable element in a NAC gene is associated with drought tolerance in maize seedlings. Nat. Commun. 2015, 6, 1–13. [Google Scholar] [CrossRef] [Green Version]
  21. Xiang, Y.; Sun, X.; Gao, S.; Qin, F.; Dai, M. Deletion of an endoplasmic reticulum stress response element in a ZmPP2C-A gene facilitates drought tolerance of maize seedlings. Mol. Plant 2017, 10, 456–469. [Google Scholar] [CrossRef] [Green Version]
  22. Dong, Z.; Xu, Z.; Xu, L.; Galli, M.; Gallavotti, A.; Dooner, H.K.; Chuck, G. Necrotic upper tips1 mimics heat and drought stress and encodes a protoxylem-specific transcription factor in maize. Proc. Natl. Acad. Sci. USA 2020, 117, 20908–20919. [Google Scholar] [CrossRef] [PubMed]
  23. Liu, B.; Zhang, B.; Yang, Z.; Liu, Y.; Yang, S.; Shi, Y.; Jiang, C.; Qin, F. Manipulating ZmEXPA4 expression ameliorates the drought-induced prolonged anthesis and silking interval in maize. Plant Cell 2021, 33, 2058–2071. [Google Scholar] [CrossRef] [PubMed]
  24. Zhang, X.; Liu, H.; Ma, X.; Zhou, G.; Ruan, H.; Cui, H.; Pang, J.; Siffat, U.K.; Zong, N.; Wang, R.; et al. Genome-wide association study and metabolic pathway prediction of barrenness in maize as a response to high planting density. J. Integr. Agric. 2021, in press. [Google Scholar]
  25. Fu, J.; Cheng, Y.; Linghu, J.; Yang, X.; Kang, L.; Zhang, Z.; Zhang, J.; He, C.; Du, X.; Peng, Z. RNA sequencing reveals the complex regulatory network in the maize kernel. Nat. Commun. 2013, 4, 1–12. [Google Scholar] [CrossRef] [Green Version]
  26. Pang, J.; Fu, J.; Zong, N.; Wang, J.; Song, D.; Zhang, X.; He, C.; Fang, T.; Zhang, H.; Fan, Y.; et al. Kernel size-related genes revealed by an integrated eQTL analysis during early maize kernel development. Plant J. 2019, 98, 19–32. [Google Scholar] [CrossRef]
  27. Bradbury, P.J.; Zhang, Z.; Kroon, D.E.; Casstevens, T.M.; Ramdoss, Y.; Buckler, E.S. TASSEL: Software for association mapping of complex traits in diverse samples. Bioinformatics 2007, 23, 2633–2635. [Google Scholar] [CrossRef]
  28. Holland, J.B.; Nyquist, W.E.; Cervantes-Martínez, C.T. Estimating and Interpreting Heritability for Plant Breeding: An Update; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2010. [Google Scholar]
  29. Messmer, R.; Fracheboud, Y.; Bänziger, M.; Vargas, M.; Stamp, P.; Ribaut, J.M. Drought stress and tropical maize: QTL-by-environment interactions and stability of QTLs across environments for yield components and secondary traits. Theor. Appl. Genet. 2009, 119, 913–930. [Google Scholar] [CrossRef] [Green Version]
  30. Xue, Y.; Warburton, M.L.; Sawkins, M.; Zhang, X.; Setter, T.; Xu, Y.; Grudloyma, P.; Gethi, J.; Ribaut, J.M.; Li, W. Genome-wide association analysis for nine agronomic traits in maize under well-watered and water-stressed conditions. Theor. Appl. Genet. 2013, 126, 2587–2596. [Google Scholar] [CrossRef]
  31. Hejnák, V.; Tatar, Ö.; Atasoy, G.D.; Martinková, J.; Çelen, A.E.; Hnilička, F.; Skalický, M. Growth and Photosynthesis of upland and pima cotton: Response to drought and heat stress. Plant Soil Environ. 2015, 62, 507–514. [Google Scholar] [CrossRef] [Green Version]
  32. Zafar, S.A.; Patil, S.B.; Uzair, M.; Fang, J.; Zhao, J.; Guo, T.; Yuan, S.; Uzair, M.; Luo, Q.; Shi, J.; et al. DEGENERATED PANICLE AND PARTIAL STERILITY 1 (DPS 1) encodes a cystathionine β–synthase domain containing protein required for anther cuticle and panicle development in rice. New Phytol. 2020, 225, 356–375. [Google Scholar] [CrossRef] [Green Version]
  33. Edmeades, G.O.; Bolanos, J.; Elings, A.; Ribaut, J.M.; Bänziger, M.; Westgate, M.E. The role and regulation of the anthesis-silking interval in maize. Physiol. Modeling Kernel Set Maize 2000, 29, 43–73. [Google Scholar]
  34. Setter, T.L. Analysis of constituents for phenotyping drought tolerance in crop improvement. Front. Physiol. 2012, 3, 180. [Google Scholar] [CrossRef] [Green Version]
  35. Li, C.; Sun, B.; Li, Y.; Liu, C.; Wu, X.; Zhang, D.; Shi, Y.; Song, Y.; Buckler, E.S.; Zhang, Z.; et al. Numerous genetic loci identified for drought tolerance in the maize nested association mapping populations. BMC Genom. 2016, 17, 1–11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Farfan, I.D.; De La Fuente, G.N.; Murray, S.C.; Isakeit, T.; Huang, P.C.; Warburton, M.; Williams, P.; Windham, G.L.; Kolomiets, M. Genome wide association study for drought, aflatoxin resistance, and important agronomic traits of maize hybrids in the sub-tropics. PLoS ONE 2015, 10, e0117737. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Hao, Z.; Li, X.; Xie, C.; Weng, J.; Li, M.; Zhang, D. Identification of functional genetic variations underlying drought tolerance in maize using SNP markers. J. Integr. Plant Biol. 2011, 53, 641–652. [Google Scholar] [CrossRef] [PubMed]
  38. Lu, Y.; Zhang, S.; Shah, T.; Xie, C.; Hao, Z.; Li, X.; Farkhari, M.; Ribaut, J.M.; Cao, M.; Rong, T. Joint linkage-linkage disequilibrium mapping is a powerful approach to detecting quantitative trait loci underlying drought tolerance in maize. Proc. Natl. Acad. Sci. USA 2010, 107, 19585–19590. [Google Scholar] [CrossRef] [Green Version]
  39. Hu, X.; Wang, G.; Du, X.; Zhang, H.; Xu, Z.; Wang, J.; Chen, G.; Wang, B.; Li, X.; Chen, X.; et al. QTL analysis across multiple environments reveals promising chromosome regions associated with yield-related traits in maize under drought conditions. Crop J. 2021, 9, 759–766. [Google Scholar] [CrossRef]
  40. Almeida, G.D.; Nair, S.; Borém, A.; Cairns, J.; Trachsel, S.; Ribaut, J.M.; Bänziger, M.; Prasanna, B.M.; Crossa, J.; Babu, R. Molecular mapping across three populations reveals a QTL hotspot region on chromosome 3 for secondary traits associated with drought tolerance in tropical maize. Mol. Breed. 2014, 34, 701–715. [Google Scholar] [CrossRef] [Green Version]
  41. Xiao, Y.N.; Li, X.H.; George, M.L.; Li, M.S.; Zhang, S.H.; Zheng, Y.L. Quantitative trait locus analysis of drought tolerance and yield in maize in China. Plant Mol. Biol. Rep. 2005, 23, 155–165. [Google Scholar] [CrossRef]
  42. Cai, H.; Chu, Q.; Gu, R.; Yuan, L.; Liu, J.; Zhang, X.; Chen, F.; Mi, G.; Zhang, F. Identification of QTLs for plant height, ear height and grain yield in maize (Zea mays L.) in response to nitrogen and phosphorus supply. Plant Breed. 2012, 131, 502–510. [Google Scholar] [CrossRef]
  43. Messmer, R.; Fracheboud, Y.; Bänziger, M.; Stamp, P.; Ribaut, J.M. Drought stress and tropical maize: QTLs for leaf greenness, plant senescence, and root capacitance. Field Crops Res. 2011, 124, 93–103. [Google Scholar] [CrossRef]
  44. Wang, L.; Zhou, Z.; Li, R.; Weng, J.; Zhang, Q.; Li, X.; Wang, B.; Zhang, W.; Song, W.; Li, X. Mapping QTL for flowering time-related traits under three plant densities in maize. Crop J. 2021, 9, 372–379. [Google Scholar] [CrossRef]
  45. Leng, P.; Ouzunova, M.; Landbeck, M.; Wenzel, G.; Eder, J.; Darnhofer, B.; Lübberstedt, T. Quantitive trait loci mapping of forage agronomic traits in six mapping populations derived from European elite maize germplasm. Plant Breed. 2018, 137, 370–378. [Google Scholar] [CrossRef]
  46. de Lorenzo, L.; Merchan, F.; Blanchet, S.; Megias, M.; Frugier, F.; Crespi, M.; Sousa, C. Differential expression of the TFIIIA regulatory pathway in response to salt stress between Medicago truncatula genotypes. Plant Physiol. 2007, 145, 1521–1532. [Google Scholar] [CrossRef] [Green Version]
  47. Yan, Z.; Jia, J.; Yan, X.; Shi, H.; Han, Y. Arabidopsis KHZ1 and KHZ2, two novel non-tandem CCCH zinc-finger and K-homolog domain proteins, have redundant roles in the regulation of flowering and senescence. Plant Mol. Biol. 2017, 95, 549–565. [Google Scholar] [CrossRef]
  48. Huang, X.; Chao, D.; Gao, J.; Zhu, M.; Shi, M.; Lin, H. A previously unknown zinc finger protein DST regulates drought and salt tolerance in rice via stomatal aperture control. Gene. Dev. 2009, 23, 1805–1817. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Mu, C.; Chen, N.; Li, X. F-box Protein arabidillo-1 promotes lateral root development by depressing the functioning of GA3 in Arabidopsis. J. Plant Biol. 2010, 53, 374–380. [Google Scholar] [CrossRef]
  50. Sharma, M.; Singh, A.; Shankar, A.; Pandey, A.; Baranwal, V.; Kapoor, S.; Tyagi, A.K.; Pandey, G.K. Comprehensive expression analysis of rice armadillo gene family during abiotic stress and development. DNA Res. 2014, 21, 267–283. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Nibau, C.; Gibbs, D.J.; Bunting, K.A.; Moody, L.A.; Smiles, E.J.; Tubby, J.A.; Bradshaw, S.J.; Coates, J.C. ARABIDILLO proteins have a novel and conserved domain structure important for the regulation of their stability. Plant Mol. Biol. 2011, 75, 77–92. [Google Scholar] [CrossRef]
Figure 1. Soil moisture monitor and plant performance under different water treatment: (A) WATERMARK MONITOR installed in the field and connected to the sensor; (B,C) plant height and ear development of drought-sensitive lines under WW and WS condition.
Figure 1. Soil moisture monitor and plant performance under different water treatment: (A) WATERMARK MONITOR installed in the field and connected to the sensor; (B,C) plant height and ear development of drought-sensitive lines under WW and WS condition.
Genes 13 00564 g001
Figure 2. Frequency distributions of ASI, EBM, and PH under two water treatments.
Figure 2. Frequency distributions of ASI, EBM, and PH under two water treatments.
Genes 13 00564 g002
Figure 3. Manhattan plots of ASI using MLM model in three years.
Figure 3. Manhattan plots of ASI using MLM model in three years.
Genes 13 00564 g003
Figure 4. Manhattan plots of ear biomass using MLM model in three years.
Figure 4. Manhattan plots of ear biomass using MLM model in three years.
Genes 13 00564 g004
Figure 5. Manhattan plots of plant height using MLM model in three years.
Figure 5. Manhattan plots of plant height using MLM model in three years.
Genes 13 00564 g005
Figure 6. Drought responsive pattern of candidate genes under WS and WW conditions.
Figure 6. Drought responsive pattern of candidate genes under WS and WW conditions.
Genes 13 00564 g006
Figure 7. LD patterns and the allele effects of the most significant SNPs for candidate genes Zm00001d013992 (A), Zm00001d020506 (C), Zm00001d029937 (E), Zm00001d029938 (G), Zm00001d039319 (I), and Zm00001d042997 (K). Results of haplotype analysis for Zm00001d013992 (B), Zm00001d020506 (D), Zm00001d029937 (F), Zm00001d029938 (H), Zm00001d039319 (J), and Zm00001d042997 (L). Blue: HapA; pink: HapB; green: HapC.
Figure 7. LD patterns and the allele effects of the most significant SNPs for candidate genes Zm00001d013992 (A), Zm00001d020506 (C), Zm00001d029937 (E), Zm00001d029938 (G), Zm00001d039319 (I), and Zm00001d042997 (K). Results of haplotype analysis for Zm00001d013992 (B), Zm00001d020506 (D), Zm00001d029937 (F), Zm00001d029938 (H), Zm00001d039319 (J), and Zm00001d042997 (L). Blue: HapA; pink: HapB; green: HapC.
Genes 13 00564 g007
Table 1. Descriptive statistics and heritability estimate for the traits of association panel.
Table 1. Descriptive statistics and heritability estimate for the traits of association panel.
TraitsRangeMean ± SDCV (%)(H2) %p Value
ASI-WS1.67–21.425.58 ± 2.3942.9789.312.6 × 10−3
ASI-WW0.56–18.154.43 ± 2.1748.8987.45
EBM-WS0.48–3.081.42 ± 0.4129.1586.783.2 × 10−5
EBM-WW0.48–3.811.71 ± 0.4727.2086.64
PH-WS43.21–169.72111.36 ± 21.5419.3492.522.7 × 10−4
PH-WW51.23–183.74127.22 ± 21.5316.9294.60
ASI, anthesis-silking interval, days; EBM, ear biomass, g; PH, plant height, cm; SD, standard deviation; CV, coefficient of variance; H2, Broad-sense heritability.
Table 2. Correlation analysis based on BLUP values across three-year environments under drought (WS, above diagonal) and well-watered (WW, under diagonal) regimes.
Table 2. Correlation analysis based on BLUP values across three-year environments under drought (WS, above diagonal) and well-watered (WW, under diagonal) regimes.
Trait ASIEBMPH
ASI −0.090.13 *
EBM−0.04 0.22 **
PH0.15 *0.18 **
ASI, anthesis-silking interval; EBM, ear biomass; PH, plant height; *, ** significant at P < 0.05, 0.01, respectively.
Table 3. Annotation of SNPs associated with ASI, EBM, and PH under multiple environments.
Table 3. Annotation of SNPs associated with ASI, EBM, and PH under multiple environments.
Traits MarkerChr.Positionp ValueR2Gene IDAnnotation
ASI-WS-18S1_935135641935135645.38 × 10−60.08207Zm00001d029938Protein ARABIDILLO 1
S1_932776411932776416.06 × 10−60.08113Zm00001d029937Glycoprotein
PZE-103003226324499131.03 × 10−70.14322Zm00001d039319Tic22-like family protein
chr3.S_18326319231833192921.01 × 10−50.07963Zm00001d042997HIT-type zinc-finger family protein
ASI-WW-18S1_932776411932776412.20 × 10−70.1079Zm00001d029937Glycoprotein
S1_932777751932777753.28 × 10−70.10511
S1_932781501932781507.29 × 10−70.09852
S1_935135641935135641.01 × 10−60.09549Zm00001d029938Protein ARABIDILLO 1
S1_935070461935070462.48 × 10−60.08831
S1_935058551935058553.76 × 10−60.08489
S1_935098921935098923.76 × 10−60.08489
S1_935106461935106463.76 × 10−60.08489
S1_935111551935111553.76 × 10−60.08489
S1_935100581935100588.64 × 10−60.07831
S1_935115211935115218.64 × 10−60.07831
S1_935130961935130968.64 × 10−60.07831
PZE-103003226324499131.64 × 10−60.10835Zm00001d039319Tic22-like family protein
chr3.S_18326319231833192921.66 × 10−60.09449Zm00001d042997HIT-type zinc-finger family protein
S3_18331545731833154571.91 × 10−50.09027
S3_18331565831833156581.91 × 10−60.09027
S3_18331691631833169161.91 × 10−60.09027
S3_18331864231833186421.91 × 10−60.09027
S3_18331540031833154005.78 × 10−60.08148
S3_18331173331833117337.14 × 10−60.07982
S3_18331177731833117777.14 × 10−60.07982
EBM-WS-17chr7.S_11628875671163167095.92 × 10−60.1034Zm00001d02050626S proteasome non-ATPase regulatory subunit 9
chr7.S_11628879171163167445.92 × 10−60.1034
chr7.S_11628879271163167455.92 × 10−60.1034
chr7.S_11628565271163136051.01 × 10−50.09798
chr7.S_11628565571163136081.01 × 10−50.09798
EBM-WW-17S7_11631557671163155761.17 × 10−60.11793Zm00001d02050626S proteasome non-ATPase regulatory subunit 9
S7_11631642571163164251.17 × 10−60.11793
S7_11631655971163165591.17 × 10−60.11793
chr7.S_11628875671163167091.29 × 10−60.12102
chr7.S_11628879171163167441.29 × 10−60.12102
chr7.S_11628879271163167451.29 × 10−60.12102
chr7.S_11628565271163136053.42 × 10−60.11084
chr7.S_11628565571163136083.42 × 10−60.11084
S7_11631442371163144231.81 × 10−60.11403
S7_11631666771163166672.11 × 10−60.11193
EBM-WS-18S5_271219445271219443.25 × 10−60.08944Zm00001d013992Pyridoxal phosphate-dependent transferase family protein
EBM-WS-20S5_271219445271219449.15 × 10−60.09491Zm00001d013992Pyridoxal phosphate-dependent transferase family protein
PH-WS-17chr2.S_686916182693219212.98 × 10−70.14098Zm00001d00393911-ß-hydroxysteroid dehydrogenase
chr2.S_686916212693219242.98 × 10−70.14098
S2_21802677022180267701.11 × 10−60.11601Zm00001d007189Uncharacterised
S2_22644987022264498702.08 × 10−60.10972GRMZM2G070937Leu-rich repeat protein kinase family protein
PH-WW-17chr2.S_686916182693219217.15 × 10−80.15528Zm00001d00393911-ß-hydroxysteroid dehydrogenase
chr2.S_686916212693219247.15 × 10−80.15528
S2_21802677022180267701.95 × 10−60.11069Zm00001d007189Uncharacterised
S2_22644987022264498709.07 × 10−60.09541GRMZM2G070937Leu-rich repeat protein kinase family protein
PH-WS-18S8_16392701181639270114.47 × 10−60.07836Zm00001d012167Silk fibroin (SF16) protein
PH-WW-18S8_16392701181639270118.67 × 10−60.07275Zm00001d012167Silk fibroin (SF16) protein
S8_16392701281639270129.66 × 10−60.07196
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Khan, S.U.; Zheng, Y.; Chachar, Z.; Zhang, X.; Zhou, G.; Zong, N.; Leng, P.; Zhao, J. Dissection of Maize Drought Tolerance at the Flowering Stage Using Genome-Wide Association Studies. Genes 2022, 13, 564. https://doi.org/10.3390/genes13040564

AMA Style

Khan SU, Zheng Y, Chachar Z, Zhang X, Zhou G, Zong N, Leng P, Zhao J. Dissection of Maize Drought Tolerance at the Flowering Stage Using Genome-Wide Association Studies. Genes. 2022; 13(4):564. https://doi.org/10.3390/genes13040564

Chicago/Turabian Style

Khan, Siffat Ullah, Yanxiao Zheng, Zaid Chachar, Xuhuan Zhang, Guyi Zhou, Na Zong, Pengfei Leng, and Jun Zhao. 2022. "Dissection of Maize Drought Tolerance at the Flowering Stage Using Genome-Wide Association Studies" Genes 13, no. 4: 564. https://doi.org/10.3390/genes13040564

APA Style

Khan, S. U., Zheng, Y., Chachar, Z., Zhang, X., Zhou, G., Zong, N., Leng, P., & Zhao, J. (2022). Dissection of Maize Drought Tolerance at the Flowering Stage Using Genome-Wide Association Studies. Genes, 13(4), 564. https://doi.org/10.3390/genes13040564

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop