Parallel tuning of semi-dwarfism via differential splicing of Brachytic1 in commercial maize and smallholder sorghum.

In the current genomic era, the search and deployment of new semi-dwarf alleles have continued to develop better plant types in all cereals. We characterized an agronomically optimal semi-dwarf mutation in Zea mays L. and a parallel polymorphism in Sorghum bicolor L. We cloned the maize brachytic1 (br1-Mu) allele by a modified PCR-based Sequence Amplified Insertion Flanking Fragment (SAIFF) approach. Histology and RNA-Seq elucidated the mechanism of semi-dwarfism. GWAS linked a sorghum plant height QTL with the Br1 homolog by resequencing a West African sorghum landraces panel. The semi-dwarf br1-Mu allele encodes an MYB transcription factor78 that positively regulates stalk cell elongation by interacting with the polar auxin pathway. Semi-dwarfism is due to differential splicing and low functional Br1 wild-type transcript expression. The sorghum ortholog, SbBr1, co-segregates with the major plant height QTL qHT7.1 and is alternatively spliced. The high frequency of the Sbbr1 allele in African landraces suggests that African smallholder farmers used the semi-dwarf allele to improve plant height in sorghum long before efforts to introduce Green Revolution-style varieties in the 1960s. Surprisingly, variants for differential splicing of Brachytic1 were found in both commercial maize and smallholder sorghum, suggesting parallel tuning of plant architecture across these systems.

Development of an efficient marker-free soybean transformation method using the novel bacterium Ochrobactrum haywardense H1.

We have discovered a novel bacterium, Ochrobactrum haywardense H1 (Oh H1), which is capable of efficient plant transformation. Ochrobactrum is a new host for Agrobacterium-derived vir and T-DNA-mediated transformation. Oh H1 is a unique, non-phytopathogenic species, categorized as a BSL-1 organism. We engineered Oh H1 with repurposed Agrobacterium virulence machinery and demonstrated Oh H1 can transform numerous dicot species and at least one monocot, sorghum. We generated a cysteine auxotrophic Oh H1-8 strain containing a binary vector system. Oh H1-8 produced transgenic soybean plants with an efficiency 1.6 times that of Agrobacterium strain AGL1 and 2.9 times that of LBA4404Thy-. Oh H1-8 successfully transformed several elite Corteva soybean varieties with T0 transformation frequency up to 35%. In addition to higher transformation efficiencies, Oh H1-8 generated high-quality, transgenic events with single-copy, plasmid backbone-free insertion at frequencies higher than AGL1. The SpcN selectable marker gene is excised using a heat shock-inducible excision system resulting in marker-free transgenic events. Approximately, 24.5% of the regenerated plants contained only a single copy of the transgene and contained no vector backbone. There were no statistically significant differences in yield comparing T3 null-segregant lines to wild-type controls. We have demonstrated that Oh H1-8, combined with spectinomycin selection, is an efficient, rapid, marker-free and yield-neutral transformation system for elite soybean.

Characterization of Proteome Variation During Modern Maize Breeding.

The success of modern maize breeding has been demonstrated by remarkable increases in productivity with tremendous modification of agricultural phenotypes over the last century. Although the underlying genetic changes of the maize adaptation from tropical to temperate regions have been extensively studied, our knowledge is limited regarding the accordance of protein and mRNA expression levels accompanying such adaptation. Here we conducted an integrative analysis of proteomic and transcriptomic changes in a maize association panel. The minimum extent of correlation between protein and RNA levels suggests that variation in mRNA expression is often not indicative of protein expression at a population scale. This is corroborated by the observation that mRNA- and protein-based coexpression networks are relatively independent of each other, and many pQTLs arise without the presence of corresponding eQTLs. Importantly, compared with transcriptome, the subtypes categorized by the proteome show a markedly high accuracy to resemble the genomic subpopulation. These findings suggest that proteome evolved under a greater evolutionary constraint than transcriptome during maize adaptation from tropical to temperate regions. Overall, the integrated multi-omics analysis provides a functional context to interpret gene expression variation during modern maize breeding.

Ectopic expression of ARGOS8 reveals a role for ethylene in root-lodging resistance in maize.

Ethylene plays a critical role in many diverse processes in plant development. Recent studies have demonstrated that overexpression of the maize ARGOS8 gene reduces the plant's response to ethylene by decreasing ethylene signaling and enhances grain yield in transgenic maize plants. The objective of this study was to determine the effects of ethylene on the development of nodal roots, which are primarily responsible for root-lodging resistance in maize. Exogenous application of the ethylene precursor 1-aminocyclopropane-1-carboxylic acid (ACC) was found to promote the emergence of nodal roots. Transcriptome analysis of nodal tissues revealed that the expression of genes involved in metabolic processes and cell wall biogenesis was upregulated in response to ACC treatment, supporting the notion that ethylene is a positive regulator for the outgrowth of young root primordia. In BSV::ARGOS8 transgenic plants with reduced ethylene sensitivity due to constitutive overexpression of ARGOS8, nodal root emergence was delayed and the promotional effect of ACC on nodal root emergence decreased. Field tests showed that the BSV::ARGOS8 plants had higher root lodging relative to non-transgenic controls. When ARGOS8 expression was controlled by the developmentally regulated promoter FTM1, which conferred ARGOS8 overexpression in adult plants but not in the nodal roots and nodes in juvenile plants, the FTM1::ARGOS8 plants had no significant difference in root lodging compared with the wild type but produced a higher grain yield. These results suggest that ethylene has a role in promoting nodal root emergence and that a delay in nodal root development has a negative effect on root-lodging resistance in maize.

Genome-Wide Analysis of Alternative Splicing during Development and Drought Stress in Maize.

Alternative splicing plays a crucial role in plant development as well as stress responses. Although alternative splicing has been studied during development and in response to stress, the interplay between these two factors remains an open question. To assess the effects of drought stress on developmentally regulated splicing in maize (Zea mays), 94 RNA-seq libraries from ear, tassel, and leaf of the B73 public inbred line were constructed at four developmental stages under both well-watered and drought conditions. This analysis was supplemented with a publicly available series of 53 libraries from developing seed, embryo, and endosperm. More than 48,000 novel isoforms, often with stage- or condition-specific expression, were uncovered, suggesting that developmentally regulated alternative splicing occurs in thousands of genes. Drought induced large developmental splicing changes in leaf and ear but relatively few in tassel. Most developmental stage-specific splicing changes affected by drought were tissue dependent, whereas stage-independent changes frequently overlapped between leaf and ear. A linear relationship was found between gene expression changes in splicing factors and alternative spicing of other genes during development. Collectively, these results demonstrate that alternative splicing is strongly associated with tissue type, developmental stage, and stress condition.

Genome-wide analysis of alternative splicing in Zea mays: landscape and genetic regulation.

Alternative splicing enhances transcriptome diversity in all eukaryotes and plays a role in plant tissue identity and stress adaptation. To catalog new maize (Zea mays) transcripts and identify genomic loci that regulate alternative splicing, we analyzed over 90 RNA-seq libraries from maize inbred lines B73 and Mo17, as well as Syn10 doubled haploid lines (progenies from B73 × Mo17). Transcript discovery was augmented with publicly available data from 14 maize tissues, expanding the maize transcriptome by more than 30,000 and increasing the percentage of intron-containing genes that undergo alternative splicing to 40%. These newly identified transcripts greatly increase the diversity of the maize proteome, sometimes coding for entirely different proteins compared with their most similar annotated isoform. In addition to increasing proteome diversity, many genes encoding novel transcripts gained an additional layer of regulation by microRNAs, often in a tissue-specific manner. We also demonstrate that the majority of genotype-specific alternative splicing can be genetically mapped, with cis-acting quantitative trait loci (QTLs) predominating. A large number of trans-acting QTLs were also apparent, with nearly half located in regions not shown to contain genes associated with splicing. Taken together, these results highlight the currently underappreciated role that alternative splicing plays in tissue identity and genotypic variation in maize.

Genome-wide analysis of leafbladeless1-regulated and phased small RNAs underscores the importance of the TAS3 ta-siRNA pathway to maize development.

Maize leafbladeless1 (lbl1) encodes a key component in the trans-acting short-interfering RNA (ta-siRNA) biogenesis pathway. Correlated with a great diversity in ta-siRNAs and the targets they regulate, the phenotypes conditioned by mutants perturbing this small RNA pathway vary extensively across species. Mutations in lbl1 result in severe developmental defects, giving rise to plants with radial, abaxialized leaves. To investigate the basis for this phenotype, we compared the small RNA content between wild-type and lbl1 seedling apices. We show that LBL1 affects the accumulation of small RNAs in all major classes, and reveal unexpected crosstalk between ta-siRNA biogenesis and other small RNA pathways regulating transposons. Interestingly, in contrast to data from other plant species, we found no evidence for the existence of phased siRNAs generated via the one-hit model. Our analysis identified nine TAS loci, all belonging to the conserved TAS3 family. Information from RNA deep sequencing and PARE analyses identified the tasiR-ARFs as the major functional ta-siRNAs in the maize vegetative apex where they regulate expression of AUXIN RESPONSE FACTOR3 (ARF3) homologs. Plants expressing a tasiR-ARF insensitive arf3a transgene recapitulate the phenotype of lbl1, providing direct evidence that deregulation of ARF3 transcription factors underlies the developmental defects of maize ta-siRNA biogenesis mutants. The phenotypes of Arabidopsis and Medicago ta-siRNA mutants, while strikingly different, likewise result from misexpression of the tasiR-ARF target ARF3. Our data indicate that diversity in TAS pathways and their targets cannot fully account for the phenotypic differences conditioned by ta-siRNA biogenesis mutants across plant species. Instead, we propose that divergence in the gene networks downstream of the ARF3 transcription factors or the spatiotemporal pattern during leaf development in which these proteins act constitute key factors underlying the distinct contributions of the ta-siRNA pathway to development in maize, Arabidopsis, and possibly other plant species as well.

Maize maintains growth in response to decreased nitrate supply through a highly dynamic and developmental stage-specific transcriptional response.

Elucidation of the gene networks underlying the response to N supply and demand will facilitate the improvement of the N uptake efficiency of plants. We undertook a transcriptomic analysis of maize to identify genes responding to both a non-growth-limiting decrease in NO3- provision and to development-based N demand changes at seven representative points across the life cycle. Gene co-expression networks were derived by cluster analysis of the transcript profiles. The majority of NO3--responsive transcription occurred at 11 (D11), 18 (D18) and 29 (D29) days after emergence, with differential expression predominating in the root at D11 and D29 and in the leaf at D18. A cluster of 98 probe sets was identified, the expression pattern of which is similar to that of the high-affinity NO3- transporter (NRT2) genes across the life cycle. The cluster is enriched with genes encoding enzymes and proteins of lipid metabolism and transport, respectively. These are candidate genes for the response of maize to N supply and demand. Only a few patterns of differential gene expression were observed over the entire life cycle; however, the composition of the classes of the genes differentially regulated at individual time points was unique, suggesting tightly controlled regulation of NO3--responsive gene expression.

Transcriptome analysis of near-isogenic lines provides molecular insights into starch biosynthesis in maize kernel.

Starch is the major component in maize kernels, providing a stable carbohydrate source for humans and livestock as well as raw material for the biofuel industry. Increasing maize kernel starch content will help meet industry demands and has the potential to increase overall yields. We developed a pair of maize near-isogenic lines (NILs) with different alleles for a starch quantitative trait locus on chromosome 3 (qHS3), resulting in different kernel starch content. To investigate the candidate genes for qHS3 and elucidate their effects on starch metabolism, RNA-Seq was performed for the developing kernels of the NILs at 14 and 21 d after pollination (DAP). Analysis of genomic and transcriptomic data identified 76 genes with nonsynonymous single nucleotide polymorphisms and 384 differentially expressed genes (DEGs) in the introgressed fragment, including a hexokinase gene, ZmHXK3a, which catalyzes the conversion of glucose to glucose-6-phosphate and may play a key role in starch metabolism. The expression pattern of all DEGs in starch metabolism shows that altered expression of the candidate genes for qHS3 promoted starch synthesis, with positive consequences for kernel starch content. These results expand the current understanding of starch biosynthesis and accumulation in maize kernels and provide potential candidate genes to increase starch content.

GIANT EMBRYO encodes CYP78A13, required for proper size balance between embryo and endosperm in rice.

Among angiosperms there is a high degree of variation in embryo/endosperm size in mature seeds. However, little is known about the molecular mechanism underlying size control between these neighboring tissues. Here we report the rice GIANT EMBRYO (GE) gene that is essential for controlling the size balance. The function of GE in each tissue is distinct, controlling cell size in the embryo and cell death in the endosperm. GE, which encodes CYP78A13, is predominantly expressed in the interfacing tissues of the both embryo and endosperm. GE expression is under negative feedback regulation; endogenous GE expression is upregulated in ge mutants. In contrast to the loss-of-function mutant with large embryo and small endosperm, GE overexpression causes a small embryo and enlarged endosperm. A complementation analysis coupled with heterofertilization showed that complementation of ge mutation in either embryo or endosperm failed to restore the wild-type embryo/endosperm ratio. Thus, embryo and endosperm interact in determining embryo/endosperm size balance. Among genes associated with embryo/endosperm size, REDUCED EMBRYO genes, whose loss-of-function causes a phenotype opposite to ge, are revealed to regulate endosperm size upstream of GE. To fully understand the embryo-endosperm size control, the genetic network of the related genes should be elucidated.

tassel-less1 encodes a boron channel protein required for inflorescence development in maize.

tassel-less1 (tls1) is a classical maize (Zea mays) inflorescence mutant. Homozygous mutant plants have no tassels or very small tassels, and ear development is also impaired. Using a positional cloning approach, ZmNIP3;1 (a NOD26-like intrinsic protein) was identified as the candidate gene for tls1. The ZmNIP3;1 gene is completely deleted in the tls1 mutant genome. Two Mutator-insertional TUSC alleles of ZmNIP3;1 exhibited tls1-like phenotypes, and allelism tests confirmed that the tls1 gene encodes ZmNIP3;1. Transgenic plants with an RNA interference (RNAi) construct to down-regulate ZmNIP3;1 also showed tls1-like phenotypes, further demonstrating that TLS1 is ZmNIP3;1. Sequence analysis suggests that ZmNIP3;1 is a boron channel protein. Foliar application of boron could rescue the tls1 phenotypes and restore the normal tassel and ear development. Gene expression analysis indicated that in comparison with that of the wild type or tls1 plants treated with boron, the transition from the vegetative to reproductive phase or the development of the floral meristem is impaired in the shoot apical meristem of the tls1 mutant plants. It is concluded that the tls1 mutant phenotypes are caused by impaired boron transport, and boron is essential for inflorescence development in maize.

Spatial gradients in cell wall composition and transcriptional profiles along elongating maize internodes.

BACKGROUND: The elongating maize internode represents a useful system for following development of cell walls in vegetative cells in the Poaceae family. Elongating internodes can be divided into four developmental zones, namely the basal intercalary meristem, above which are found the elongation, transition and maturation zones. Cells in the basal meristem and elongation zones contain mainly primary walls, while secondary cell wall deposition accelerates in the transition zone and predominates in the maturation zone. RESULTS: The major wall components cellulose, lignin and glucuronoarabinoxylan (GAX) increased without any abrupt changes across the elongation, transition and maturation zones, although GAX appeared to increase more between the elongation and transition zones. Microarray analyses show that transcript abundance of key glycosyl transferase genes known to be involved in wall synthesis or re-modelling did not match the increases in cellulose, GAX and lignin. Rather, transcript levels of many of these genes were low in the meristematic and elongation zones, quickly increased to maximal levels in the transition zone and lower sections of the maturation zone, and generally decreased in the upper maturation zone sections. Genes with transcript profiles showing this pattern included secondary cell wall CesA genes, GT43 genes, some ß-expansins, UDP-Xylose synthase and UDP-Glucose pyrophosphorylase, some xyloglucan endotransglycosylases/hydrolases, genes involved in monolignol biosynthesis, and NAM and MYB transcription factor genes. CONCLUSIONS: The data indicated that the enzymic products of genes involved in cell wall synthesis and modification remain active right along the maturation zone of elongating maize internodes, despite the fact that corresponding transcript levels peak earlier, near or in the transition zone.

Comparative transcriptome profiling of maize coleoptilar nodes during shoot-borne root initiation.

Maize (Zea mays) develops an extensive shoot-borne root system to secure water and nutrient uptake and to provide anchorage in the soil. In this study, early coleoptilar node (first shoot node) development was subjected to a detailed morphological and histological analysis. Subsequently, microarray profiling via hybridization of oligonucleotide microarrays representing transcripts of 31,355 unique maize genes at three early stages of coleoptilar node development was performed. These pairwise comparisons of wild-type versus mutant rootless concerning crown and seminal roots (rtcs) coleoptilar nodes that do not initiate shoot-borne roots revealed 828 unique transcripts that displayed RTCS-dependent expression. A stage-specific functional analysis revealed overrepresentation of "cell wall," "stress," and "development"-related transcripts among the differentially expressed genes. Differential expression of a subset of 15 of 828 genes identified by these microarray experiments was independently confirmed by quantitative real-time-polymerase chain reaction. In silico promoter analyses revealed that 100 differentially expressed genes contained at least one LATERAL ORGAN BOUNDARIES domain (LBD) motif within 1 kb upstream of the ATG start codon. Electrophoretic mobility shift assay experiments demonstrated RTCS binding for four of these promoter sequences, supporting the notion that differentially accumulated genes containing LBD motifs are likely direct downstream targets of RTCS.

Cell wall modifications in maize pulvini in response to gravitational stress.

Changes in cell wall polysaccharides, transcript abundance, metabolite profiles, and hormone concentrations were monitored in the upper and lower regions of maize (Zea mays) pulvini in response to gravistimulation, during which maize plants placed in a horizontal position returned to the vertical orientation. Heteroxylan levels increased in the lower regions of the pulvini, together with lignin, but xyloglucans and heteromannan contents decreased. The degree of substitution of heteroxylan with arabinofuranosyl residues decreased in the lower pulvini, which exhibited increased mechanical strength as the plants returned to the vertical position. Few or no changes in noncellulosic wall polysaccharides could be detected on the upper side of the pulvinus, and crystalline cellulose content remained essentially constant in both the upper and lower pulvinus. Microarray analyses showed that spatial and temporal changes in transcript profiles were consistent with the changes in wall composition that were observed in the lower regions of the pulvinus. In addition, the microarray analyses indicated that metabolic pathways leading to the biosynthesis of phytohormones were differentially activated in the upper and lower regions of the pulvinus in response to gravistimulation. Metabolite profiles and measured hormone concentrations were consistent with the microarray data, insofar as auxin, physiologically active gibberellic acid, and metabolites potentially involved in lignin biosynthesis increased in the elongating cells of the lower pulvinus.

Genome-wide expression quantitative trait loci (eQTL) analysis in maize.

BACKGROUND: Expression QTL analyses have shed light on transcriptional regulation in numerous species of plants, animals, and yeasts. These microarray-based analyses identify regulators of gene expression as either cis-acting factors that regulate proximal genes, or trans-acting factors that function through a variety of mechanisms to affect transcript abundance of unlinked genes. RESULTS: A hydroponics-based genetical genomics study in roots of a Zea mays IBM2 Syn10 double haploid population identified tens of thousands of cis-acting and trans-acting eQTL. Cases of false-positive eQTL, which results from the lack of complete genomic sequences from both parental genomes, were described. A candidate gene for a trans-acting regulatory factor was identified through positional cloning. The unexpected regulatory function of a class I glutamine amidotransferase controls the expression of an ABA 8'-hydroxylase pseudogene. CONCLUSIONS: Identification of a candidate gene underlying a trans-eQTL demonstrated the feasibility of eQTL cloning in maize and could help to understand the mechanism of gene expression regulation. Lack of complete genome sequences from both parents could cause the identification of false-positive cis- and trans-acting eQTL.

Digital gene expression signatures for maize development.

Genome-wide expression signatures detect specific perturbations in developmental programs and contribute to functional resolution of key regulatory networks. In maize (Zea mays) inflorescences, mutations in the RAMOSA (RA) genes affect the determinacy of axillary meristems and thus alter branching patterns, an important agronomic trait. In this work, we developed and tested a framework for analysis of tag-based, digital gene expression profiles using Illumina's high-throughput sequencing technology and the newly assembled B73 maize reference genome. We also used a mutation in the RA3 gene to identify putative expression signatures specific to stem cell fate in axillary meristem determinacy. The RA3 gene encodes a trehalose-6-phosphate phosphatase and may act at the interface between developmental and metabolic processes. Deep sequencing of digital gene expression libraries, representing three biological replicate ear samples from wild-type and ra3 plants, generated 27 million 20- to 21-nucleotide reads with frequencies spanning 4 orders of magnitude. Unique sequence tags were anchored to 3'-ends of individual transcripts by DpnII and NlaIII digests, which were multiplexed during sequencing. We mapped 86% of nonredundant signature tags to the maize genome, which associated with 37,117 gene models and unannotated regions of expression. In total, 66% of genes were detected by at least nine reads in immature maize ears. We used comparative genomics to leverage existing information from Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa) in functional analyses of differentially expressed maize genes. Results from this study provide a basis for the analysis of short-read expression data in maize and resolved specific expression signatures that will help define mechanisms of action for the RA3 gene.

A customized gene expression microarray reveals that the brittle stem phenotype fs2 of barley is attributable to a retroelement in the HvCesA4 cellulose synthase gene.

The barley (Hordeum vulgare) brittle stem mutants, fs2, designated X054 and M245, have reduced levels of crystalline cellulose compared with their parental lines Ohichi and Shiroseto. A custom-designed microarray, based on long oligonucleotide technology and including genes involved in cell wall metabolism, revealed that transcript levels of very few genes were altered in the elongation zone of stem internodes, but these included a marked decrease in mRNA for the HvCesA4 cellulose synthase gene of both mutants. In contrast, the abundance of several hundred transcripts changed in the upper, maturation zones of stem internodes, which presumably reflected pleiotropic responses to a weakened cell wall that resulted from the primary genetic lesion. Sequencing of the HvCesA4 genes revealed the presence of a 964-bp solo long terminal repeat of a Copia-like retroelement in the first intron of the HvCesA4 genes of both mutant lines. The retroelement appears to interfere with transcription of the HvCesA4 gene or with processing of the mRNA, and this is likely to account for the lower crystalline cellulose content and lower stem strength of the mutants. The HvCesA4 gene maps to a position on chromosome 1H of barley that coincides with the previously reported position of fs2.

Gene network reconstruction from transcriptional dynamics under kinetic model uncertainty: a case for the second derivative.

MOTIVATION: Measurements of gene expression over time enable the reconstruction of transcriptional networks. However, Bayesian networks and many other current reconstruction methods rely on assumptions that conflict with the differential equations that describe transcriptional kinetics. Practical approximations of kinetic models would enable inferring causal relationships between genes from expression data of microarray, tag-based and conventional platforms, but conclusions are sensitive to the assumptions made. RESULTS: The representation of a sufficiently large portion of genome enables computation of an upper bound on how much confidence one may place in influences between genes on the basis of expression data. Information about which genes encode transcription factors is not necessary but may be incorporated if available. The methodology is generalized to cover cases in which expression measurements are missing for many of the genes that might control the transcription of the genes of interest. The assumption that the gene expression level is roughly proportional to the rate of translation led to better empirical performance than did either the assumption that the gene expression level is roughly proportional to the protein level or the Bayesian model average of both assumptions. AVAILABILITY: http://www.oisb.ca points to R code implementing the methods (R Development Core Team 2004). SUPPLEMENTARY INFORMATION: http://www.davidbickel.com.

Allelic genome structural variations in maize detected by array comparative genome hybridization.

DNA polymorphisms such as insertion/deletions and duplications affecting genome segments larger than 1 kb are known as copy-number variations (CNVs) or structural variations (SVs). They have been recently studied in animals and humans by using array-comparative genome hybridization (aCGH), and have been associated with several human diseases. Their presence and phenotypic effects in plants have not been investigated on a genomic scale, although individual structural variations affecting traits have been described. We used aCGH to investigate the presence of CNVs in maize by comparing the genome of 13 maize inbred lines to B73. Analysis of hybridization signal ratios of 60,472 60-mer oligonucleotide probes between inbreds in relation to their location in the reference genome (B73) allowed us to identify clusters of probes that deviated from the ratio expected for equal copy-numbers. We found CNVs distributed along the maize genome in all chromosome arms. They occur with appreciable frequency in different germplasm subgroups, suggesting ancient origin. Validation of several CNV regions showed both insertion/deletions and copy-number differences. The nature of CNVs detected suggests CNVs might have a considerable impact on plant phenotypes, including disease response and heterosis.

Mapping regulatory variants controlling gene expression in drought response and tolerance in maize.

BACKGROUND: Gene expression is a key determinant of cellular response. Natural variation in gene expression bridges genetic variation to phenotypic alteration. Identification of the regulatory variants controlling the gene expression in response to drought, a major environmental threat of crop production worldwide, is of great value for drought-tolerant gene identification. RESULTS: A total of 627 RNA-seq analyses are performed for 224 maize accessions which represent a wide genetic diversity under three water regimes; 73,573 eQTLs are detected for about 30,000 expressing genes with high-density genome-wide single nucleotide polymorphisms, reflecting a comprehensive and dynamic genetic architecture of gene expression in response to drought. The regulatory variants controlling the gene expression constitutively or drought-dynamically are unraveled. Focusing on dynamic regulatory variants resolved to genes encoding transcription factors, a drought-responsive network reflecting a hierarchy of transcription factors and their target genes is built. Moreover, 97 genes are prioritized to associate with drought tolerance due to their expression variations through the Mendelian randomization analysis. One of the candidate genes, Abscisic acid 8'-hydroxylase, is verified to play a negative role in plant drought tolerance. CONCLUSIONS: This study unravels the effects of genetic variants on gene expression dynamics in drought response which allows us to better understand the role of distal and proximal genetic effects on gene expression and phenotypic plasticity. The prioritized drought-associated genes may serve as direct targets for functional investigation or allelic mining.