The early maize (Zea mays) seed undergoes several developmental stages after double fertilization to become fully differentiated within a short period of time, but the genetic control of this highly dynamic and complex developmental process remains largely unknown. Here, we report a high temporal-resolution investigation of transcriptomes using 31 samples collected at an interval of 4 or 6 h within the first six days of seed development. These time-course transcriptomes were clearly separated into four distinct groups corresponding to the stages of double fertilization, coenocyte formation, cellularization, and differentiation. A total of 22,790 expressed genes including 1415 transcription factors (TFs) were detected in early stages of maize seed development. In particular, 1093 genes including 110 TFs were specifically expressed in the seed and displayed high temporal specificity by expressing only in particular period of early seed development. There were 160, 22, 112, and 569 seed-specific genes predominantly expressed in the first 16 h after pollination, coenocyte formation, cellularization, and differentiation stage, respectively. In addition, network analysis predicted 31,256 interactions among 1317 TFs and 14,540 genes. The high temporal-resolution transcriptome atlas reported here provides an important resource for future functional study to unravel the genetic control of seed development.

Understanding genetic control of tassel and ear architecture in maize (Zea mays L. ssp. mays) is important due to their relationship with grain yield. High resolution QTL mapping is critical for understanding the underlying molecular basis of phenotypic variation. Advanced populations, such as recombinant inbred lines, have been broadly adopted for QTL mapping; however, construction of large advanced generation crop populations is time-consuming and costly. The rapidly declining cost of genotyping due to recent advances in next-generation sequencing technologies has generated new possibilities for QTL mapping using large early generation populations. A set of 708 F2 progeny derived from inbreds Chang7-2 and 787 were generated and genotyped by whole genome low-coverage genotyping-by-sequencing method (average 0.04×). A genetic map containing 6,533 bin-markers was constructed based on the parental SNPs and a sliding-window method, spanning a total genetic distance of 1,396 cM. The high quality and accuracy of this map was validated by the identification of two well-studied genes, r1, a qualitative trait locus for color of silk (chromosome 10) and ba1 for tassel branch number (chromosome 3). Three traits of tassel and ear architecture were evaluated in this population, a total of 10 QTL were detected using a permutation-based-significance threshold, seven of which overlapped with reported QTL. Three genes (GRMZM2G316366, GRMZM2G492156 and GRMZM5G805008) encoding MADS-box domain proteins and a BTB/POZ domain protein were located in the small intervals of qTBN5 and qTBN7 (~800 Kb and 1.6 Mb in length, respectively) and may be involved in patterning of tassel architecture. The small physical intervals of most QTL indicate high-resolution mapping is obtainable with this method. We constructed an ultra-high-dentisy linkage map for the large early generation population in maize. Our study provides an efficient approach for fast detection of quantitative loci responsible for complex trait variation with high accuracy, thus helping to dissect the underlying molecular basis of phenotypic variation and accelerate improvement of crop breeding in a cost-effective fashion.

Drought is one of the major factors limiting global maize production. Exposure to long-term drought conditions inhibits growth and leads to yield losses. Although several drought-responsive genes have been identified and functionally analyzed, the mechanisms underlying responses to drought and water recovery treatments have not been fully elucidated. To characterize how maize seedling respond to drought stress at the transcriptional level, we analyzed physiological responses and differentially expressed genes (DEGs) in the inbred line B73 under water deficit and recovery conditions. The data for relative leaf water content, leaf size, and photosynthesis-related parameters indicated that drought stress significantly repressed maize seedling growth. Further RNA sequencing analysis revealed that 6107 DEGs were responsive to drought stress and water recovery, with more down-regulated than up-regulated genes. Among the DEGs, the photosynthesis- and hormone-related genes were enriched in responses to drought stress and re-watering. Additionally, transcription factor genes from 37 families were differentially expressed among the three analyzed time-points. Gene ontology enrichment analyses of the DEGs indicated that 50 GO terms, including those related to photosynthesis, carbohydrate metabolism, oxidoreductase activities, nutrient metabolism and other drought-responsive pathways, were over-represented in the drought-treated seedlings. The content of gibberellin in drought treatment seedlings was decreased compared to that of control seedlings, while abscisic acid showed accumulated in the drought treated plants. The deep analysis of DEGs related to cell wall development indicated that these genes were prone to be down-regulated at drought treatment stage. Many genes that are differentially expressed in responses to drought stress and water recovery conditions affect photosynthetic systems and hormone biosynthesis. The identified DEGs, especially those encoding transcription factors, represent potential targets for developing drought-tolerant maize lines.

Doubled haploid technology has been widely applied to multiple plant species and is recognized as one of the most important technologies for improving crop breeding efficiency. Although mutations in MATRILINEAL/Zea mays PHOSPHOLIPASE A1/NOT LIKE DAD (MTL/ZmPLA1/NLD) and Zea mays DOMAIN OF UNKNOWN FUNCTION 679 MEMBRANE PROTEIN (ZmDMP) have been shown to generate haploids in maize, knowledge of the genetic basis of haploid induction (HI) remains incomplete. Therefore, cloning of new genes underlying HI is important for further elucidating its genetic architecture. Here, we found that loss-of-function mutations of Zea mays PHOSPHOLIPASE D3 (ZmPLD3), one of the members from the phospholipase D subfamily, could trigger maternal HI in maize. ZmPLD3 was identified through a reverse genetic strategy based on analysis of pollen-specifically expressed phospholipases, followed by validation through the clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9 (CRISPR–Cas9) system. Mutations of ZmPLD3 resulted in a haploid induction rate (HIR) similar to that of mtl/zmpla1/nld and showed synergistic effects rather than functional redundancy on tripling the HIR (from 1.19% to 4.13%) in the presence of mtl/zmpla1/nld. RNA-seq profiling of mature pollen indicated that a large number of pollen-specific differentially expressed genes were enriched in processes related to gametogenesis development, such as pollen tube development and cell communication, during the double-fertilization process. In addition, ZmPLD3 is highly conserved among cereals, highlighting the potential application of these in vivo haploid-inducer lines for other important crop plant species. Collectively, our discovery identifies a novel gene underlying in vivo maternal HI and provides possibility of breeding haploid inducers with further improved HIR.

Broomcorn millet (Panicum miliaceum L.) has strong tolerance to abiotic stresses, and is probably one of the oldest crops, with its earliest cultivation that dated back to ca. ~10,000 years. We report here its genome assembly through a combination of PacBio sequencing, BioNano, and Hi-C (in vivo) mapping. The 18 super scaffolds cover ~95.6% of the estimated genome (~887.8 Mb). There are 63,671 protein-coding genes annotated in this tetraploid genome. About ~86.2% of the syntenic genes in foxtail millet have two homologous copies in broomcorn millet, indicating rare gene loss after tetraploidization in broomcorn millet. Phylogenetic analysis reveals that broomcorn millet and foxtail millet diverged around ~13.1 Million years ago (Mya), while the lineage specific tetraploidization of broomcorn millet may be happened within ~5.91 million years. The genome is not only beneficial for the genome assisted breeding of broomcorn millet, but also an important resource for other Panicum species.