Drought stress is one of the major limiting factors for maize production. With the availability of maize B73 reference genome and whole-genome resequencing of 15 maize inbreds, common variants (CV) and clustering analyses were applied to identify non-synonymous SNPs (nsSNPs) and corresponding candidate genes for drought tolerance. A total of 524 nsSNPs that were associated with 271 candidate genes involved in plant hormone regulation, carbohydrate and sugar metabolism, signaling molecules regulation, redox reaction and acclimation of photosynthesis to environment were detected by CV and cluster analyses. Most of the nsSNPs identified were clustered in bin 1.07 region that harbored six previously reported QTL with relatively high phenotypic variation explained for drought tolerance. Genes Ontology (GO) analysis of candidate genes revealed that there were 35 GO terms related to biotic stimulus and membrane-bounded organelle, showing significant differences between the candidate genes and the reference B73 background. Changes of expression level in these candidate genes for drought tolerance were detected using RNA sequencing for fertilized ovary, basal leaf meristem tissue and roots collected under drought stressed and well-watered conditions. The results indicated that 70% of candidate genes showed significantly expression changes under two water treatments and our strategies for mining candidate genes are feasible and relatively efficient. Our results successfully revealed candidate nsSNPs and associated genes for drought tolerance by comparative sequence analysis of 16 maize inbred lines. Both methods we applied were proved to be efficient for identifying candidate genes for complex traits through the next-generation sequencing technologies (NGS). These selected genes will not only facilitate understanding of genetic basis of drought stress response, but also accelerate genetic improvement through marker-assisted selection in maize.

DNA methylation is important for plant growth, development, and stress response. To understand DNA methylation dynamics in maize roots under water stress (WS), we reanalyzed DNA methylation sequencing data to profile DNA methylation and the gene expression landscape of two inbred lines with different drought sensitivities, as well as two of their derived recombination inbred lines (RILs). Combined with genotyping-by-sequencing, we found that the inheritance pattern of DNA methylation between RILs and parental lines was sequence-dependent. Increased DNA methylation levels were observed under WS and the methylome of drought-tolerant inbred lines were much more stable than that of the drought-sensitive inbred lines. Distinctive differentially methylated genes were found among diverse genetic backgrounds, suggesting that inbred lines with different drought sensitivities may have responded to stress in varying ways. Gene body DNA methylation showed a negative correlation with gene expression but a positive correlation with exon splicing events. Furthermore, a positive correlation of a varying extent was observed between small interfering RNA (siRNA) and DNA methylation, which at different genic regions. The response of siRNAs under WS was consistent with the differential DNA methylation. Taken together, our data can be useful in deciphering the roles of DNA methylation in plant drought-tolerance variations and in emphasizing its function in alternative splicing.

Drought stress is known to significantly limit crop growth and productivity. Lateral organ boundary domain (LBD) transcription factors—particularly class-I members—play essential roles in plant development and biotic stress. However, little information is available on class-II LBD genes related to abiotic stress in maize. Here, we cloned a maize class-II LBD transcription factor, ZmLBD5, and identified its function in drought stress. Transient expression, transactivation, and dimerization assays demonstrated that ZmLBD5 was localized in the nucleus, without transactivation, and could form a homodimer or heterodimer. Promoter analysis demonstrated that multiple drought-stress-related and ABA response cis-acting elements are present in the promoter region of ZmLBD5. Overexpression of ZmLBD5 in Arabidopsis promotes plant growth under normal conditions, and suppresses drought tolerance under drought conditions. Furthermore, the overexpression of ZmLBD5 increased the water loss rate, stomatal number, and stomatal apertures. DAB and NBT staining demonstrated that the reactive oxygen species (ROS) decreased in ZmLBD5-overexpressed Arabidopsis. A physiological index assay also revealed that SOD and POD activities in ZmLBD5-overexpressed Arabidopsis were higher than those in wild-type Arabidopsis. These results revealed the role of ZmLBD5 in drought stress by regulating ROS levels.

Plant growth and development are closely related to phosphate (Pi) and auxin. However, data regarding auxin response factors (ARFs) and their response to phosphate in maize are limited. Here, we isolated ZmARF4 in maize and dissected its biological function response to Pi stress. Overexpression of ZmARF4 in Arabidopsis confers tolerance of Pi deficiency with better root morphology than wild-type. Overexpressed ZmARF4 can partially restore the absence of lateral roots in mutant arf7 arf19. The ZmARF4 overexpression promoted Pi remobilization and up-regulated AtRNS1, under Pi limitation while it down-regulated the expression of the anthocyanin biosynthesis genes AtDFR and AtANS. A continuous detection revealed higher activity of promoter in the Pi-tolerant maize P178 line than in the sensitive 9782 line under low-Pi conditions. Meanwhile, GUS activity was specifically detected in new leaves and the stele of roots in transgenic offspring. ZmARF4 was localized to the nucleus and cytoplasm of the mesophyll protoplast and interacted with ZmILL4 and ZmChc5, which mediate lateral root initiation and defense response, respectively. ZmARF4 overexpression also conferred salinity and osmotic stress tolerance in Arabidopsis. Overall, our findings suggest that ZmARF4, a pleiotropic gene, modulates multiple stress signaling pathways, and thus, could be a candidate gene for engineering plants with multiple stress adaptation.

Maize (Zea mays) is an important multi-functional crop. The growth and yield of maize are severely affected by drought stress. Previous studies have shown that microRNAs (miRNAs) in maize play important roles in response to abiotic stress; however, their roles in response to drought stress in maize roots is unclear. In our study, we found 375 miRNAs in the roots of 16 inbred lines. Of the 16 lines, zma-MIR168, zma-MIR156, and zma-MIR166 were highly expressed, whereas zma-MIR399, zma-MIR2218, and zma-MIR2275 exhibited low expression levels. The expression patterns of miRNA in parental lines and their derived RILs are different. Over 50% of miRNAs exhibited a lower expression in recombinant inbred lines than in parents. The expression of 50 miRNAs was significantly altered under water stress (WS) in at least three inbred lines, and the expression of miRNAs in drought-tolerant lines changed markedly. To better understand the reasons for miRNA response to drought, the degree of histone modifications for miRNA genes was estimated. The methylation level of H3K4 and H3K9 in miRNA precursor regions changed more noticeably after WS, but no such phenomenon was seen for DNA methylation and m6A modification. After the prediction of miRNA targets using psRNATarget and psRobot, we used correlation analysis and qRT-PCR to further investigate the relationship between miRNAs and target genes. We found that 87 miRNA–target pairs were significantly negatively correlated. In addition, a weighted gene co-expression network analysis using miRNAs, as well as their predicted targets, was conducted to reveal that miR159, miR394, and miR319 may be related to maize root growth. The results demonstrated that miRNAs might play essential roles in the response to drought stress.

Maize was originally domesticated in a tropical environment but is now widely cultivated at temperate latitudes. Temperate and tropical maize populations have diverged both genotypically and phenotypically. Tropical maize lines grown in temperate environments usually exhibit delayed flowering, pollination, and seed set, which reduces their grain yield relative to temperate adapted maize lines. One potential mechanism by which temperate maize may have adapted to a new environment is novel transposable element insertions, which can influence gene regulation. Recent advances in sequencing technology have made it possible to study variation in transposon content and insertion location in large sets of maize lines. In total, 274,408 non-redundant TEs (NRTEs) were identified using resequencing data generated from 83 maize inbred lines. The locations of DNA TEs and copia-superfamily retrotransposons showed significant positive correlations with gene density and genetic recombination rates, whereas gypsy-superfamily retrotransposons showed a negative correlation with these two parameters. Compared to tropical maize, temperate maize had fewer unique NRTEs but higher insertion frequency, lower background recombination rates, and higher linkage disequilibrium, with more NRTEs close to flowering and stress-related genes in the genome. Association mapping demonstrated that the presence/absence of 48 NRTEs was associated with flowering time and that expression of neighboring genes differed between haplotypes where a NRTE was present or absent. This study suggests that NRTEs may have played an important role in creating the variation in gene regulation that enabled the rapid adaptation of maize to diverse environments.