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.

No abstract

SUMMARYRNA editing modifies cytidines (C) to uridines (U) at specific sites in the transcripts of mitochondria and plastids, altering the amino acid specified by the DNA sequence. Here we report the identification of a critical editing factor of mitochondrial nad7 transcript via molecular characterization of a small kernel 1 (smk1) mutant in Zea mays (maize). Mutations in Smk1 arrest both the embryo and endosperm development. Cloning of Smk1 indicates that it encodes an E-subclass pentatricopeptide repeat (PPR) protein that is targeted to mitochondria. Loss of SMK1 function abolishes the C ? U editing at the nad7-836 site, leading to the retention of a proline codon that is edited to encode leucine in the wild type. The smk1 mutant showed dramatically reduced complex-I assembly and NADH dehydrogenase activity, and abnormal biogenesis of the mitochondria. Analysis of the ortholog in Oryza sativa (rice) reveals that rice SMK1 has a conserved function in C ? U editing of the mitochondrial nad7-836 site. T-DNA knock-out mutants showed abnormal embryo and endosperm development, resulting in embryo or seedling lethality. The leucine at NAD7-279 is highly conserved from bacteria to flowering plants, and analysis of genome sequences from many plants revealed a molecular coevolution between the requirement for C ? U editing at this site and the existence of an SMK1 homolog. These results demonstrate that Smk1 encodes a PPR-E protein that is required for nad7-836 editing, and this editing is critical to NAD7 function in complex-I assembly in mitochondria, and hence to embryo and endosperm development in maize and rice.

AUXIN RESPONSE FACTORS (ARFs) are plant-specific transcription factors (TFs) that couple perception of the hormone auxin to gene expression programs essential to all land plants. As with many large TF families, a key question is whether individual members determine developmental specificity by binding distinct target genes. We use DAP-seq to generate genome-wide in vitro TF:DNA interaction maps for fourteen maize ARFs from the evolutionarily conserved A and B clades. Comparative analysis reveal a high degree of binding site overlap for ARFs of the same clade, but largely distinct clade A and B binding. Many sites are however co-occupied by ARFs from both clades, suggesting transcriptional coordination for many genes. Among these, we investigate known QTLs and use machine learning to predict the impact of cis-regulatory variation. Overall, large-scale comparative analysis of ARF binding suggests that auxin response specificity may be determined by factors other than individual ARF binding site selection.

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.

IspH is a key enzyme in the last step of the methyl-D-erythritol-4-phosphate (MEP) pathway. Loss of function of IspH can often result in complete yellow or albino phenotype in many plants. Here, we report the characterization of a recessive mutant of maize, zebra7 (zb7), showing transverse green/yellow striped leaves in young plants. The yellow bands of the mutant have decreased levels of chlorophylls and carotenoids with delayed chloroplast development. Low temperature suppressed mutant phenotype, while alternate light/dark cycle or high temperature enlarged the yellow section. Map-based cloning demonstrated that zb7 encodes the IspH protein with a mis-sense mutation in a conserved region. Transgenic silencing of Zb7 in maize resulted in complete albino plantlets that are aborted in a few weeks, confirming that Zb7 is important in the early stages of maize chloroplast development. Zb7 is constitutively expressed and its expression subject to a 16-h light/8-h dark cycle regulation. Our results suggest that the less effective or unstable IspH in zb7 mutant, together with its diurnal expression, are mechanistically accounted for the zebra phenotype. The increased IspH mRNA in the leaves of zb7 at the late development stage may explain the restoration of mutant phenotype in mature stages.

RNA editing is a posttranscriptional process that is important in mitochondria and plastids of higher plants. All RNA editing-specific trans-factors reported so far belong to PLS-class of pentatricopeptide repeat (PPR) proteins. Here, we report the map-based cloning and molecular characterization of a defective kernel mutant dek39 in maize. Loss of Dek39 function leads to delayed embryogenesis and endosperm development, reduced kernel size, and seedling lethality. Dek39 encodes an E subclass PPR protein that targets to both mitochondria and chloroplasts, and is involved in RNA editing in mitochondrial NADH dehydrogenase3 (nad3) at nad3-247 and nad3-275. C-to-U editing of nad3-275 is not conserved and even lost in Arabidopsis, consistent with the idea that no close DEK39 homologs are present in Arabidopsis. However, the amino acids generated by editing nad3-247 and nad3-275 are highly conserved in many other plant species, and the reductions of editing at these two sites decrease the activity of mitochondria NADH dehydrogenase complex I, indicating that the alteration of amino acid sequence is necessary for Nad3 function. Our results indicate that Dek39 encodes an E sub-class PPR protein that is involved in RNA editing of multiple sites and is necessary for seed development of maize.

Robertson’s Mutator (Mu) system has been used in large scale mutagenesis in maize, exploiting its high mutation frequency, controllability, preferential insertion in genes, and independence of donor location. Eight Mutator elements have been fully characterized (Mu1, Mu2 /Mu1.7, Mu3, Mu4, Mu5, Mu6/7, Mu8, MuDR), and three are defined by TIR (Mu10, Mu11 and Mu12). The genome sequencing revealed a complex family of Mu-like-elements (MULEs) in the B73 genome. In this article, we report the identification of a new Mu element, named Mu13. Mu13 showed typical Mu characteristics by having a ∼220 bp TIR, creating a 9 bp target site duplication upon insertion, yet the internal sequence is completely different from previously identified Mu elements. Mu13 is not present in the B73 genome or a Zea mays subsp. parviglumis accession, but in W22 and several inbreds that found the Robertson’s Mutator line. Analysis of mutants isolated from the UniformMu mutagenic population indicated that the Mu13 element is active in transposition. Two novel insertions were found in expressed genes. To test other unknown Mu elements, we selected six new Mu elements from the B73 genome. Southern analysis indicated that most of these elements were present in the UniformMu lines. From these results, we conclude that Mu13 is a new and active Mu element that significantly contributed to the mutagenesis in the UniformMu population. The Robertson’s Mutator line may harbor other unknown active Mu elements.