Maize is a major crop and a model plant for studying C4 photosynthesis and leaf development. However, a genomewide regulatory network of leaf development is not yet available. This knowledge is useful for developing C3 crops to perform C4 photosynthesis for enhanced yields. Here, using 22 transcriptomes of developing maize leaves from dry seeds to 192 h post imbibition, we studied gene up-and down-regulation and functional transition during leaf development and inferred sets of strongly coexpressed genes. More significantly, we developed a method to predict transcription factor binding sites (TFBSs) and their cognate transcription factors (TFs) using genomic sequence and transcriptomic data. The method requires not only evolutionary conservation of candidate TFBSs and sets of strongly coexpressed genes but also that the genes in a gene set share the same Gene Ontology term so that they are involved in the same biological function. In addition, we developed another method to predict maize TF-TFBS pairs using known TF-TFBS pairs in Arabidopsis or rice. From these efforts, we predicted 1,340 novel TFBSs and 253 new TF-TFBS pairs in the maize genome, far exceeding the 30 TF-TFBS pairs currently known in maize. In most cases studied by both methods, the two methods gave similar predictions. In vitro tests of 12 predicted TF-TFBS interactions showed that our methods perform well. Our study has significantly expanded our knowledge on the regulatory network involved in maize leaf development. maize transcriptomes | coexpressed genes | cis binding site M aize (Zea mays) is a major crop and a model plant for studying C4 photosynthesis and leaf development. However, the regulatory network that controls maize leaf development is still not well understood. In fact, the number of known maize transcription factor (TF)-binding sites (TFBSs) is far smaller than that for Arabidopsis thaliana (1-3).To understand better the regulation of maize leaf development, several recent studies used next-generation sequencing (NGS) technologies to survey transcriptomic differences among maize leaf cell types and developmental stages. The first large-scale study was by Li et al. (7) investigated the transcriptomes of Kranz (i.e., the foliar leaf blade) and non-Kranz (the husk leaf sheath) maize leaves to identify cohorts of genes associated with procambium initiation and vascular patterning. Recently, Wang et al. (8) conducted comparative transcriptomic and metabolomic analyses of developing leaves in maize and rice and identified putative structural and regulatory components important for C4 and C3 photosynthesis. These studies provided insights into the regulatory mechanisms underlying the development of Kranz leaf anatomy in maize.In the present study, we obtained nine transcriptomes of the second leaf from 84-192 h post imbibition at 12-h (6:00 AM and 6:00 PM) or 24-h (6:00 PM only) intervals. Together with the 13
SignificanceMaize is a major crop and a model plant for studying C4 leaf development. However, its regulatory network of leaf ...