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.

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.

Superior field performance of waxy corn engineered using CRISPR-Cas9.

We created waxy corn hybrids by CRISPR-Cas9 editing of a waxy allele in 12 elite inbred maize lines, a process that was more than a year faster than conventional trait introgression using backcrossing and marker-assisted selection. Field trials at 25 locations showed that CRISPR-waxy hybrids were agronomically superior to introgressed hybrids, producing on average 5.5 bushels per acre higher yield.

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.

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.

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.

Southern-by-Sequencing: A Robust Screening Approach for Molecular Characterization of Genetically Modified Crops.

Molecular characterization of events is an integral part of the advancement process during genetically modified (GM) crop product development. Assessment of these events is traditionally accomplished by polymerase chain reaction (PCR) and Southern blot analyses. Southern blot analysis can be time-consuming and comparatively expensive and does not provide sequence-level detail. We have developed a sequence-based application, Southern-by-Sequencing (SbS), utilizing sequence capture coupled with next-generation sequencing (NGS) technology to replace Southern blot analysis for event selection in a high-throughput molecular characterization environment. SbS is accomplished by hybridizing indexed and pooled whole-genome DNA libraries from GM plants to biotinylated probes designed to target the sequence of transformation plasmids used to generate events within the pool. This sequence capture process enriches the sequence data obtained for targeted regions of interest (transformation plasmid DNA). Taking advantage of the DNA adjacent to the targeted bases (referred to as next-to-target sequence) that accompanies the targeted transformation plasmid sequence, the data analysis detects plasmid-to-genome and plasmid-to-plasmid junctions introduced during insertion into the plant genome. Analysis of these junction sequences provides sequence-level information as to the following: the number of insertion loci including detection of unlinked, independently segregating, small DNA fragments; copy number; rearrangements, truncations, or deletions of the intended insertion DNA; and the presence of transformation plasmid backbone sequences. This molecular evidence from SbS analysis is used to characterize and select GM plants meeting optimal molecular characterization criteria. SbS technology has proven to be a robust event screening tool for use in a high-throughput molecular characterization environment.

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.

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.

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.

Analytical method evaluation and discovery of variation within maize varieties in the context of food safety: transcript profiling and metabolomics.

Profiling techniques such as microarrays, proteomics, and metabolomics are used widely to assess the overall effects of genetic background, environmental stimuli, growth stage, or transgene expression in plants. To assess the potential regulatory use of these techniques in agricultural biotechnology, we carried out microarray and metabolomic studies of 3 different tissues from 11 conventional maize varieties. We measured technical variations for both microarrays and metabolomics, compared results from individual plants and corresponding pooled samples, and documented variations detected among different varieties with individual plants or pooled samples. Both microarray and metabolomic technologies are reproducible and can be used to detect plant-to-plant and variety-to-variety differences. A pooling strategy lowered sample variations for both microarray and metabolomics while capturing variety-to-variety variation. However, unknown genomic sequences differing between maize varieties might hinder the application of microarrays. High-throughput metabolomics could be useful as a tool for the characterization of transgenic crops. However, researchers will have to take into consideration the impact on the detection and quantitation of a wide range of metabolites on experimental design as well as validation and interpretation of results.

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.

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.

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.

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.

Maize global transcriptomics reveals pervasive leaf diurnal rhythms but rhythms in developing ears are largely limited to the core oscillator.

BACKGROUND: Plant diurnal rhythms are vital environmental adaptations to coordinate internal physiological responses to alternating day-night cycles. A comprehensive view of diurnal biology has been lacking for maize (Zea mays), a major world crop. METHODOLOGY: A photosynthetic tissue, the leaf, and a non-photosynthetic tissue, the developing ear, were sampled under natural field conditions. Genome-wide transcript profiling was conducted on a high-density 105 K Agilent microarray to investigate diurnal rhythms. CONCLUSIONS: In both leaves and ears, the core oscillators were intact and diurnally cycling. Maize core oscillator genes are found to be largely conserved with their Arabidopsis counterparts. Diurnal gene regulation occurs in leaves, with some 23% of expressed transcripts exhibiting a diurnal cycling pattern. These transcripts can be assigned to over 1700 gene ontology functional terms, underscoring the pervasive impact of diurnal rhythms on plant biology. Considering the peak expression time for each diurnally regulated gene, and its corresponding functional assignment, most gene functions display temporal enrichment in the day, often with distinct patterns, such as dawn or midday preferred, indicating that there is a staged procession of biological events undulating with the diurnal cycle. Notably, many gene functions display a bimodal enrichment flanking the midday photosynthetic maximum, with an initial peak in mid-morning followed by another peak during the afternoon/evening. In contrast to leaves, in developing ears as few as 47 gene transcripts are diurnally regulated, and this set of transcripts includes primarily the core oscillators. In developing ears, which are largely shielded from light, the core oscillator therefore is intact with little outward effect on transcription.