Maize (Zea mays) is an excellent cereal model for research on seed development because of its relatively large size for both embryo and endosperm. Despite the importance of seed in agriculture, the genome-wide transcriptome pattern throughout seed development has not been well characterized. Using high-throughput RNA sequencing, we developed a spatiotemporal transcriptome atlas of B73 maize seed development based on 53 samples from fertilization to maturity for embryo, endosperm, and whole seed tissues. A total of 26,105 genes were found to be involved in programming seed development, including 1,614 transcription factors. Global comparisons of gene expression highlighted the fundamental transcriptomic reprogramming and the phases of development. Coexpression analysis provided further insight into the dynamic reprogramming of the transcriptome by revealing functional transitions during maturation. Combined with the published nonseed high-throughput RNA sequencing data, we identified 91 transcription factors and 1,167 other seed-specific genes, which should help elucidate key mechanisms and regulatory networks that underlie seed development. In addition, correlation of gene expression with the pattern of DNA methylation revealed that hypomethylation of the gene body region should be an important factor for the expressional activation of seed-specific genes, especially for extremely highly expressed genes such as zeins. This study provides a valuable resource for understanding the genetic control of seed development of monocotyledon plants.

These authors contributed equally to this work. SUMMARYPlants can respond to environmental changes with various mechanisms occurred at transcriptional and translational levels. Thus far, there have been relatively extensive understandings of stress responses of plants on transcriptional level, while little information is known about that on translational level. To uncover the landscape of translation in plants in response to drought stress, we performed the recently developed ribosome profiling assay with maize seedlings growing under normal and drought conditions. Comparative analysis of the ribosome profiling data and the RNA-seq data showed that the fold changes of gene expression at transcriptional level were moderately correlated with that of translational level globally (R 2 = 0.69). However, less than half of the responsive genes were shared by transcription and translation under drought condition, suggesting that drought stress can introduce transcriptional and translational responses independently. We found that the translational efficiencies of 931 genes were changed significantly in response to drought stress. Further analysis revealed that the translational efficiencies of genes were highly influenced by their sequence features including GC content, length of coding sequences and normalized minimal free energy. In addition, we detected potential translation of 3063 upstream open reading frames (uORFs) on 2558 genes and these uORFs may affect the translational efficiency of downstream main open reading frames (ORFs). Our study indicates that plant can respond to drought stress with highly dynamic translational mechanism, that acting synergistically with that of transcription.

Genetic imprinting is a specific epigenetic phenomenon in which a subset of genes is expressed depending on their parentof-origin. Two types of chromatin modifications, DNA methylation and histone modification, are generally believed to be involved in the regulation of imprinting. However, the genome-wide correlation between allele-specific chromatin modifications and imprinted gene expression in maize remains elusive. Here we report genome-wide high resolution allelespecific maps of DNA methylation and histone H3 lysine 27 trimethylation (H3K27me3) in maize endosperm. For DNA methylation, thousands of parent-of-origin dependent differentially methylated regions (pDMRs) were identified. All pDMRs were uniformly paternally hypermethylated and maternally hypomethylated. We also identified 1131 allele-specific H3K27me3 peaks that are preferentially present in the maternal alleles. Maternally expressed imprinted genes (MEGs) and paternally expressed imprinted genes (PEGs) had different patterns of allele-specific DNA methylation and H3K27me3. Allele-specific expression of MEGs was not directly related to allele-specific H3K27me3, and only a subset of MEGs was associated with maternal-specific DNA demethylation, which was primarily located in the upstream and 59 portion of gene body regions. In contrast, allele-specific expression of a majority of PEGs was related to maternal-specific H3K27me3, with a subgroup of PEGs also associated with maternal-specific DNA demethylation. Both pDMRs and maternal H3K27me3 peaks associated with PEGs are enriched in gene body regions. Our results indicate highly complex patterns of regulation on genetic imprinting in maize endosperm.

BackgroundUnderstanding 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.ResultsA 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.ConclusionsWe 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.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-15-433) contains supplementary material, which is available to authorized users.

Maize (Zea mays) root system architecture (RSA) mediates the key functions of plant anchorage and acquisition of nutrients and water. In this study, a set of 204 recombinant inbred lines (RILs) was derived from the widely adapted Chinese hybrid ZD958(Zheng58 × Chang7‐2), genotyped by sequencing (GBS) and evaluated as seedlings for 24 RSA related traits divided into primary, seminal and total root classes. Significant differences between the means of the parental phenotypes were detected for 18 traits, and extensive transgressive segregation in the RIL population was observed for all traits. Moderate to strong relationships among the traits were discovered. A total of 62 quantitative trait loci (QTL) were identified that individually explained from 1.6% to 11.6% (total root dry weight/total seedling shoot dry weight) of the phenotypic variation. Eighteen, 24 and 20 QTL were identified for primary, seminal and total root classes of traits, respectively. We found hotspots of 5, 3, 4 and 12 QTL in maize chromosome bins 2.06, 3.02‐03, 9.02‐04, and 9.05‐06, respectively, implicating the presence of root gene clusters or pleiotropic effects. These results characterized the phenotypic variation and genetic architecture of seedling RSA in a population derived from a successful maize hybrid.