Genomic distribution of maize facultative heterochromatin marked by trimethylation of H3K27

Establishing an optimized ATAC-seq protocol for the maize

The advent of next-generation sequencing in crop improvement offers unprecedented insights into the chromatin landscape closely linked to gene activity governing key traits in plant development and adaptation. Particularly in maize, its dynamic chromatin structure is found to collaborate with massive transcriptional variations across tissues and developmental stages, implying intricate regulatory mechanisms, which highlights the importance of integrating chromatin information into breeding strategies for precise gene controls. The depiction of maize chromatin architecture using Assay for Transposase Accessible Chromatin with high-throughput sequencing (ATAC-seq) provides great opportunities to investigate cis-regulatory elements, which is crucial for crop improvement. In this context, we developed an easy-to-implement ATAC-seq protocol for maize with fewer nuclei and simple equipment. We demonstrate a streamlined ATAC-seq protocol with four key steps for maize in which nuclei purification can be achieved without cell sorting and using only a standard bench-top centrifuge. Our protocol, coupled with the bioinformatic analysis, including validation by read length periodicity, key metrics, and correlation with transcript abundance, provides a precise and efficient assessment of the maize chromatin landscape. Beyond its application to maize, our testing design holds the potential to be applied to other crops or other tissues, especially for those with limited size and amount, establishing a robust foundation for chromatin structure studies in diverse crop species.

SET DOMAIN GROUP 711-mediated H3K27me3 methylation of cytokinin metabolism genes regulates organ size in rice

Organ size shapes plant architecture during rice (Oryza sativa) growth and development, affecting key factors influencing yield, such as plant height, leaf size, and seed size. Here, we report that the rice Enhancer of Zeste [E(z)] homolog SET DOMAIN GROUP 711 (OsSDG711) regulates organ size in rice. Knockout of OsSDG711 produced shorter plants with smaller leaves, thinner stems, and smaller grains. We demonstrate that OsSDG711 affects organ size by reducing cell length and width and increasing cell number in leaves, stems, and grains. The result of chromatin immunoprecipitation followed by deep sequencing (ChIP-seq) using an antitrimethylation of histone H3 lysine 27 (H3K27me3) antibody showed that the levels of H3K27me3 associated with cytokinin oxidase/dehydrogenase genes (OsCKXs) were lower in the OsSDG711 knockout line Ossdg711. ChIP-qPCR assays indicated that OsSDG711 regulates the expression of OsCKX genes through H3K27me3 histone modification. Importantly, we show that OsSDG711 directly binds to the promoters of these OsCKX genes. Furthermore, we measured significantly lower cytokinin contents in Ossdg711 plants than in wild-type plants. Overall, our results reveal an epigenetic mechanism based on OsSDG711-mediated modulation of H3K27me3 levels to regulate the expression of genes involved in the cytokinin metabolism pathway and control organ development in rice. OsSDG711 may be an untapped epigenetic resource for ideal plant type improvement.

Identification of watermelon H3K4 and H3K27 genes and their expression profiles during watermelon fruit development

BACKGROUND

The ClaH3K4s and ClaH3K27s gene families are subfamilies of the SET family, each with a highly conserved SET structure domain and a PHD structural domain. Both participate in histone protein methylation, which affects the chromosome structure and gene expression, and is essential for fruit growth and development.

METHODS AND RESULTS

In order to demonstrate the structure and expression characteristics of ClaH3K4s and ClaH3K27s in watermelon, members of the watermelon H3K4 and H3K27 gene families were identified, and their chromosomal localization, gene structure, and protein structural domains were analyzed. The phylogeny and covariance of the gene families with other species were subsequently determined, and the expression profiles were obtained by performing RNA-Seq and qRT-PCR. The watermelon genome had five H3K4 genes with 3207-8043 bp nucleotide sequence lengths and four H3K27 genes with a 1107-5499 bp nucleotide sequence. Synteny analysis revealed the close relationship between watermelon and cucumber, with the majority of members displaying a one-to-one covariance. Approximately half of the 'Hua-Jing 13 watermelon' ClaH3K4s and ClaH3K27s genes were expressed more in the late fruit development stages, while the changes were minimal for the remaining half. H3K4-2 expression was observed to be slightly greater on day 21 compared to other periods. Moreover, ClaH3K27-1 and ClaH3K27-2 were hardly expressed throughout the developing period, and ClaH3K27-4 exhibited the highest expression.

CONCLUSION

These results serve as a basis for further functional characterization of the H3K4 and H3K27 genes in the fruit development of watermelon.

SDG26 Is Involved in Trichome Control in Arabidopsis thaliana: Affecting Phytohormones and Adjusting Accumulation of H3K27me3 on Genes Related to Trichome Growth and Development

Plant trichomes formed by specialized epidermal cells play a role in protecting plants from biotic and abiotic stresses and can also influence the economic and ornamental value of plant products. Therefore, further studies on the molecular mechanisms of plant trichome growth and development are important for understanding trichome formation and agricultural production. SET Domain Group 26 (SDG26) is a histone lysine methyltransferase. Currently, the molecular mechanism by which SDG26 regulates the growth and development of Arabidopsis leaf trichomes is still unclear. We found that the mutant of Arabidopsis (sdg26) possessed more trichomes on its rosette leaves compared to the wild type (Col-0), and the trichome density per unit area of sdg26 is significantly higher than that of Col-0. The content of cytokinins and jasmonic acid was higher in sdg26 than in Col-0, while the content of salicylic acid was lower in sdg26 than in Col-0, which is conducive to trichome growth. By measuring the expression levels of trichome-related genes, we found that the expression of genes that positively regulate trichome growth and development were up-regulated, while the negatively regulated genes were down-regulated in sdg26. Through chromatin immunoprecipitation sequencing (ChIP-seq) analysis, we found that SDG26 can directly regulate the expression of genes related to trichome growth and development such as ZFP1, ZFP5, ZFP6, GL3, MYB23, MYC1, TT8, GL1, GIS2, IPT1, IPT3, and IPT5 by increasing the accumulation of H3K27me3 on these genes, which further affects the growth and development of trichomes. This study reveals the mechanism by which SDG26 affects the growth and development of trichomes through histone methylation. The current study provides a theoretical basis for studying the molecular mechanism of histone methylation in regulating leaf trichome growth and development and perhaps guiding the development of new crop varieties.

Epigenetic Changes Occurring in Plant Inbreeding

Inbreeding is the crossing of closely related individuals in nature or a plantation or self-pollinating plants, which produces plants with high homozygosity. This process can reduce genetic diversity in the offspring and decrease heterozygosity, whereas inbred depression (ID) can often reduce viability. Inbred depression is common in plants and animals and has played a significant role in evolution. In the review, we aim to show that inbreeding can, through the action of epigenetic mechanisms, affect gene expression, resulting in changes in the metabolism and phenotype of organisms. This is particularly important in plant breeding because epigenetic profiles can be linked to the deterioration or improvement of agriculturally important characteristics.

Conserved noncoding sequences and de novo Mutator insertion alleles are imprinted in maize

Genomic imprinting is an epigenetic phenomenon in which differential allele expression occurs in a parent-of-origin-dependent manner. Imprinting in plants is tightly linked to transposable elements (TEs), and it has been hypothesized that genomic imprinting may be a consequence of demethylation of TEs. Here, we performed high-throughput sequencing of ribonucleic acids from four maize (Zea mays) endosperms that segregated newly silenced Mutator (Mu) transposons and identified 110 paternally expressed imprinted genes (PEGs) and 139 maternally expressed imprinted genes (MEGs). Additionally, two potentially novel paternally suppressed MEGs are associated with de novo Mu insertions. In addition, we find evidence for parent-of-origin effects on expression of 407 conserved noncoding sequences (CNSs) in maize endosperm. The imprinted CNSs are largely localized within genic regions and near genes, but the imprinting status of the CNSs are largely independent of their associated genes. Both imprinted CNSs and PEGs have been subject to relaxed selection. However, our data suggest that although MEGs were already subject to a higher mutation rate prior to their being imprinted, imprinting may be the cause of the relaxed selection of PEGs. In addition, although DNA methylation is lower in the maternal alleles of both the maternally and paternally expressed CNSs (mat and pat CNSs), the difference between the two alleles in H3K27me3 levels was only observed in pat CNSs. Together, our findings point to the importance of both transposons and CNSs in genomic imprinting in maize.

Histone variants and modifications during abiotic stress response

Plants have developed multiple mechanisms as an adaptive response to abiotic stresses, such as salinity, drought, heat, cold, and oxidative stress. Understanding these regulatory networks is critical for coping with the negative impact of abiotic stress on crop productivity worldwide and, eventually, for the rational design of strategies to improve plant performance. Plant alterations upon stress are driven by changes in transcriptional regulation, which rely on locus-specific changes in chromatin accessibility. This process encompasses post-translational modifications of histone proteins that alter the DNA-histones binding, the exchange of canonical histones by variants that modify chromatin conformation, and DNA methylation, which has an implication in the silencing and activation of hypervariable genes. Here, we review the current understanding of the role of the major epigenetic modifications during the abiotic stress response and discuss the intricate relationship among them.

Will epigenetics be a key player in crop breeding?

If food and feed production are to keep up with world demand in the face of climate change, continued progress in understanding and utilizing both genetic and epigenetic sources of crop variation is necessary. Progress in plant breeding has traditionally been thought to be due to selection for spontaneous DNA sequence mutations that impart desirable phenotypes. These spontaneous mutations can expand phenotypic diversity, from which breeders can select agronomically useful traits. However, it has become clear that phenotypic diversity can be generated even when the genome sequence is unaltered. Epigenetic gene regulation is a mechanism by which genome expression is regulated without altering the DNA sequence. With the development of high throughput DNA sequencers, it has become possible to analyze the epigenetic state of the whole genome, which is termed the epigenome. These techniques enable us to identify spontaneous epigenetic mutations (epimutations) with high throughput and identify the epimutations that lead to increased phenotypic diversity. These epimutations can create new phenotypes and the causative epimutations can be inherited over generations. There is evidence of selected agronomic traits being conditioned by heritable epimutations, and breeders may have historically selected for epiallele-conditioned agronomic traits. These results imply that not only DNA sequence diversity, but the diversity of epigenetic states can contribute to increased phenotypic diversity. However, since the modes of induction and transmission of epialleles and their stability differ from that of genetic alleles, the importance of inheritance as classically defined also differs. For example, there may be a difference between the types of epigenetic inheritance important to crop breeding and crop production. The former may depend more on longer-term inheritance whereas the latter may simply take advantage of shorter-term phenomena. With the advances in our understanding of epigenetics, epigenetics may bring new perspectives for crop improvement, such as the use of epigenetic variation or epigenome editing in breeding. In this review, we will introduce the role of epigenetic variation in plant breeding, largely focusing on DNA methylation, and conclude by asking to what extent new knowledge of epigenetics in crop breeding has led to documented cases of its successful use.

High-quality genome and methylomes illustrate features underlying evolutionary success of oaks

The genus Quercus, which emerged ∼55 million years ago during globally warm temperatures, diversified into ∼450 extant species. We present a high-quality de novo genome assembly of a California endemic oak, Quercus lobata, revealing features consistent with oak evolutionary success. Effective population size remained large throughout history despite declining since early Miocene. Analysis of 39,373 mapped protein-coding genes outlined copious duplications consistent with genetic and phenotypic diversity, both by retention of genes created during the ancient γ whole genome hexaploid duplication event and by tandem duplication within families, including numerous resistance genes and a very large block of duplicated DUF247 genes, which have been found to be associated with self-incompatibility in grasses. An additional surprising finding is that subcontext-specific patterns of DNA methylation associated with transposable elements reveal broadly-distributed heterochromatin in intergenic regions, similar to grasses. Collectively, these features promote genetic and phenotypic variation that would facilitate adaptability to changing environments.

The Importance of Networking: Plant Polycomb Repressive Complex 2 and Its Interactors

Polycomb Repressive Complex 2 (PRC2) is arguably the best-known plant complex of the Polycomb Group (PcG) pathway, formed by a group of proteins that epigenetically represses gene expression. PRC2-mediated deposition of H3K27me3 has amply been studied in Arabidopsis and, more recently, data from other plant model species has also been published, allowing for an increasing knowledge of PRC2 activities and target genes. How PRC2 molecular functions are regulated and how PRC2 is recruited to discrete chromatin regions are questions that have brought more attention in recent years. A mechanism to modulate PRC2-mediated activity is through its interaction with other protein partners or accessory proteins. Current evidence for PRC2 interactors has demonstrated the complexity of its protein network and how far we are from fully understanding the impact of these interactions on the activities of PRC2 core subunits and on the formation of new PRC2 versions. This review presents a list of PRC2 interactors, emphasizing their mechanistic action upon PRC2 functions and their effects on transcriptional regulation.

RLB (RICE LATERAL BRANCH) recruits PRC2-mediated H3K27 tri-methylation on OsCKX4 to regulate lateral branching

Lateral branches such as shoot and panicle are determining factors and target traits for rice (Oryza sativa L.) yield improvement. Cytokinin promotes rice lateral branching; however, the mechanism underlying the fine-tuning of cytokinin homeostasis in rice branching remains largely unknown. Here, we report the map-based cloning of RICE LATERAL BRANCH (RLB) encoding a nuclear-localized, KNOX-type homeobox protein from a rice cytokinin-deficient mutant showing more tillers, sparser panicles, defected floret morphology as well as attenuated shoot regeneration from callus. RLB directly binds to the promoter and represses the transcription of OsCKX4, a cytokinin oxidase gene with high abundance in panicle branch meristem. OsCKX4 over-expression lines phenocopied rlb, which showed upregulated OsCKX4 levels. Meanwhile, RLB physically binds to Polycomb repressive complex 2 (PRC2) components OsEMF2b and co-localized with H3K27me3, a suppressing histone modification mediated by PRC2, in the OsCKX4 promoter. We proposed that RLB recruits PRC2 to the OsCKX4 promoter to epigenetically repress its transcription, which suppresses the catabolism of cytokinin, thereby promoting rice lateral branching. Moreover, antisense inhibition of OsCKX4 under the LOG promoter successfully increased panicle size and spikelet number per plant without affecting other major agronomic traits. This study provides insight into cytokinin homeostasis, lateral branching in plants, and also promising target genes for rice genetic improvement.

Polycomb Repressive Complex 2 in Eukaryotes-An Evolutionary Perspective

Polycomb repressive complex 2 (PRC2) represents a group of evolutionarily conserved multi-subunit complexes that repress gene transcription by introducing trimethylation of lysine 27 on histone 3 (H3K27me3). PRC2 activity is of key importance for cell identity specification and developmental phase transitions in animals and plants. The composition, biochemistry, and developmental function of PRC2 in animal and flowering plant model species are relatively well described. Recent evidence demonstrates the presence of PRC2 complexes in various eukaryotic supergroups, suggesting conservation of the complex and its function. Here, we provide an overview of the current understanding of PRC2-mediated repression in different representatives of eukaryotic supergroups with a focus on the green lineage. By comparison of PRC2 in different eukaryotes, we highlight the possible common and diverged features suggesting evolutionary implications and outline emerging questions and directions for future research of polycomb repression and its evolution.

A novel insight into transcriptional and epigenetic regulation underlying sex expression and flower development in melon (Cucumis melo L.)

Melon (Cucumis melo L.) is an important cucurbit and has been considered as a model plant for studying sex determination. The four most common sexual morphotypes in melon are monoecious (A-G-M), gynoecious (--ggM-), andromonoecious (A-G-mm), and hermaphrodite (--ggmm). Sex expression in melons is complex, as the genes and associated networks that govern the sex expression are not fully explored. Recently, RNA-seq transcriptomic profiling, ChIP-qPCR analysis integrated with gene ontology annotation and Kyoto Encyclopedia of Genes and Genomes pathways predicted the differentially expressed genes including sex-specific ACS and ACO genes, in regulating the sex-expression, phytohormonal cross-talk, signal transduction, and secondary metabolism in melons. Integration of transcriptional control through genetic interaction in between the ACS7, ACS11, and WIP1 in epistatic or hypostatic manner, along with the recruitment of H3K9ac and H3K27me3, epigenetically, overall determine sex expression. Alignment of protein sequences for establishing phylogenetic evolution, motif comparison, and protein-protein interaction supported the structural conservation while presence of the conserved hydrophilic and charged residues across the diverged evolutionary group predicted the functional conservation of the ACS protein. Presence of the putative cis-binding elements or DNA motifs, and its further comparison with DAP-seq-based cistrome and epicistrome of Arabidopsis, unraveled strong ancestry of melons with Arabidopsis. Motif comparison analysis also characterized putative genes and transcription factors involved in ethylene biosynthesis, signal transduction, and hormonal cross-talk related to sex expression. Overall, we have comprehensively reviewed research findings for a deeper insight into transcriptional and epigenetic regulation of sex expression and flower development in melons.

Epigenome and Epitranscriptome: Potential Resources for Crop Improvement

Crop breeding faces the challenge of increasing food demand, especially under climatic changes. Conventional breeding has relied on genetic diversity by combining alleles to obtain desired traits. In recent years, research on epigenetics and epitranscriptomics has shown that epigenetic and epitranscriptomic diversity provides additional sources for crop breeding and harnessing epigenetic and epitranscriptomic regulation through biotechnologies has great potential for crop improvement. Here, we review epigenome and epitranscriptome variations during plant development and in response to environmental stress as well as the available sources for epiallele formation. We also discuss the possible strategies for applying epialleles and epitranscriptome engineering in crop breeding.

Epigenome plasticity in plants

Plant intra-individual and inter-individual variation can be determined by the epigenome, a set of covalent modifications of DNA and chromatin that can alter genome structure and activity without changes to the genome sequence. The epigenome of plant cells is plastic, that is, it can change in response to internal or external cues, such as during development or due to environmental changes, to create a memory of such events. Ongoing advances in technologies to read and write epigenomic patterns with increasing resolution, scale and precision are enabling the extent of plant epigenome variation to be more extensively characterized and functionally interrogated. In this Review, we discuss epigenome dynamics and variation within plants during development and in response to environmental changes, including stress, as well as between plants. We review known or potential functions of such plasticity and emphasize the importance of investigating the causality of epigenomic changes. Finally, we discuss emerging technologies that may underpin future research into plant epigenome plasticity.

Comprehensive Mechanism of Gene Silencing and Its Role in Plant Growth and Development

Gene silencing is a negative feedback mechanism that regulates gene expression to define cell fate and also regulates metabolism and gene expression throughout the life of an organism. In plants, gene silencing occurs via transcriptional gene silencing (TGS) and post-transcriptional gene silencing (PTGS). TGS obscures transcription via the methylation of 5' untranslated region (5'UTR), whereas PTGS causes the methylation of a coding region to result in transcript degradation. In this review, we summarized the history and molecular mechanisms of gene silencing and underlined its specific role in plant growth and crop production.

There is another: H3K27me3-mediated genomic imprinting

DNA methylation has long been considered the primary epigenetic mediator of genomic imprinting in mammals. Recent epigenetic profiling during early mouse development revealed the presence of domains of trimethylation of lysine 27 on histone H3 (H3K27me3) and chromatin compaction specifically at the maternally derived allele, independent of DNA methylation. Within these domains, genes are exclusively expressed from the paternally derived allele. This novel mechanism of noncanonical imprinting plays a key role in the development of mouse extraembryonic tissues and in the regulation of imprinted X-chromosome inactivation, highlighting the importance of parentally inherited epigenetic histone modifications. Here, we discuss the mechanisms underlying H3K27me3-mediated noncanonical imprinting in perspective of the dynamic chromatin landscape during early mouse development and explore evolutionary origins of noncanonical imprinting.

Genome-wide analysis of long noncoding RNAs, 24-nt siRNAs, DNA methylation and H3K27me3 marks in Brassica rapa

Long noncoding RNAs (lncRNAs) are RNA fragments that generally do not code for a protein but are involved in epigenetic gene regulation. In this study, lncRNAs of Brassica rapa were classified into long intergenic noncoding RNAs, natural antisense RNAs, and intronic noncoding RNAs and their expression analyzed in relation to genome-wide 24-nt small interfering RNAs (siRNAs), DNA methylation, and histone H3 lysine 27 trimethylation marks (H3K27me3). More than 65% of the lncRNAs analyzed consisted of one exon, and more than 55% overlapped with inverted repeat regions (IRRs). Overlap of lncRNAs with IRRs or genomic regions encoding for 24-nt siRNAs resulted in increased DNA methylation levels when both were present. LncRNA did not overlap greatly with H3K27me3 marks, but the expression level of intronic noncoding RNAs that did coincide with H3K27me3 marks was higher than without H3K27me3 marks. The Brassica genus comprises important vegetables and oil seed crops grown across the world. B. rapa is a diploid (AA genome) thought to be one of the ancestral species of both B. juncea (AABB genome) and B. napus (AACC) through genome merging (allotetrapolyploidization). Complex genome restructuring and epigenetic alterations are thought to be involved in these allotetrapolyploidization events. Comparison of lncRNAs between B. rapa and B. nigra, B. oleracea, B. juncea, and B. napus showed the highest conservation with B. oleracea. This study presents a comprehensive analysis of the epigenome structure of B. rapa at multi-epigenetic levels (siRNAs, DNA methylation, H3K27me3, and lncRNAs) and identified a suite of candidate lncRNAs that may be epigenetically regulated in the Brassica genus.

SDG711 Is Involved in Rice Seed Development through Regulation of Starch Metabolism Gene Expression in Coordination with Other Histone Modifications

SDG711 is a histone H3K27me2/3 transmethylase in rice, a homolog of CLF in Arabidopsis, and plays key roles in regulating flowering time and panicle development. In this work, we investigated the role of SDG711 in rice seed development. Overexpression and downregulation of SDG711 lead to a decrease and increase in the expression level of genes related to starch accumulation, resulting in smaller seeds or even seed abortion. ChIP assay showed that SDG711-mediated H3K27me3 changed significantly in genes related to endosperm development, and SDG711 can directly bind to the gene body region of several starch synthesis genes and amylase genes. In addition, H3K4me3 and H3K9ac modifications also cooperate with H3K27me3 to regulate the development of the endosperm. Our results suggest that the crosstalk between SDG711-mediated H3K27me3 and H3K4me3, and H3K9ac are involved in starch accumulation to control normal seed development.

New insights into homoeologous copy number variations in the hexaploid wheat genome

Bread wheat is an allohexaploid species originating from two successive and recent rounds of hybridization between three diploid species that were very similar in terms of chromosome number, genome size, TE content, gene content and synteny. As a result, it has long been considered that most of the genes were in three pairs of homoeologous copies. However, these so-called triads represent only one half of wheat genes, while the remaining half belong to homoeologous groups with various number of copies across subgenomes. In this study, we examined and compared the distribution, conservation, function, expression and epigenetic profiles of triads with homoeologous groups having undergone a deletion (dyads) or a duplication (tetrads) in one subgenome. We show that dyads and tetrads are mostly located in distal regions and have lower expression level and breadth than triads. Moreover, they are enriched in functions related to adaptation and more associated with the repressive H3K27me3 modification. Altogether, these results suggest that triads mainly correspond to housekeeping genes and are part of the core genome, while dyads and tetrads belong to the Triticeae dispensable genome. In addition, by comparing the different categories of dyads and tetrads, we hypothesize that, unlike most of the allopolyploid species, subgenome dominance and biased fractionation are absent in hexaploid wheat. Differences observed between the three subgenomes are more likely related to two successive and ongoing waves of post-polyploid diploidization, that had impacted A and B more significantly than D, as a result of the evolutionary history of hexaploid wheat.

Polycomb Repressive Complex 2-mediated histone modification H3K27me3 is associated with embryogenic potential in Norway spruce

Epigenetic reprogramming during germ cell formation is essential to gain pluripotency and thus embryogenic potential. The histone modification H3K27me3, which is catalysed by the Polycomb repressive complex 2 (PRC2), regulates important developmental processes in both plants and animals, and defects in PRC2 components cause pleiotropic developmental abnormalities. Nevertheless, the role of H3K27me3 in determining embryogenic potential in gymnosperms is still elusive. To address this, we generated H3K27me3 profiles of Norway spruce (Picea abies) embryonic callus and non-embryogenic callus using CUT&RUN, which is a powerful method for chromatin profiling. Here, we show that H3K27me3 mainly accumulated in genic regions in the Norway spruce genome, similarly to what is observed in other plant species. Interestingly, H3K27me3 levels in embryonic callus were much lower than those in the other examined tissues, but markedly increased upon embryo induction. These results show that H3K27me3 levels are associated with the embryogenic potential of a given tissue, and that the early phase of somatic embryogenesis is accompanied by changes in H3K27me3 levels. Thus, our study provides novel insights into the role of this epigenetic mark in spruce embryogenesis and reinforces the importance of PRC2 as a key regulator of cell fate determination across different plant species.

Identification of polycomb repressive complex 1 and 2 core components in hexaploid bread wheat

BACKGROUND

Polycomb repressive complexes 1 and 2 play important roles in epigenetic gene regulation by posttranslationally modifying specific histone residues. Polycomb repressive complex 2 is responsible for the trimethylation of lysine 27 on histone H3; Polycomb repressive complex 1 catalyzes the monoubiquitination of histone H2A at lysine 119. Both complexes have been thoroughly studied in Arabidopsis, but the evolution of polycomb group gene families in monocots, particularly those with complex allopolyploid origins, is unknown.

RESULTS

Here, we present the in silico identification of the Polycomb repressive complex 1 and 2 (PRC2, PRC1) subunits in allohexaploid bread wheat, the reconstruction of their evolutionary history and a transcriptional analysis over a series of 33 developmental stages. We identified four main subunits of PRC2 [E(z), Su(z), FIE and MSI] and three main subunits of PRC1 (Pc, Psc and Sce) and determined their chromosomal locations. We found that most of the genes coding for subunit proteins are present as paralogs in bread wheat. Using bread wheat RNA-seq data from different tissues and developmental stages throughout plant ontogenesis revealed variable transcriptional activity for individual paralogs. Phylogenetic analysis showed a high level of protein conservation among temperate cereals.

CONCLUSIONS

The identification and chromosomal location of the Polycomb repressive complex 1 and 2 core components in bread wheat may enable a deeper understanding of developmental processes, including vernalization, in commonly grown winter wheat.

Differential chromatin accessibility landscape reveals structural and functional features of the allopolyploid wheat chromosomes

Background

Our understanding of how the complexity of the wheat genome influences the distribution of chromatin states along the homoeologous chromosomes is limited. Using a differential nuclease sensitivity assay, we investigate the chromatin states of the coding and repetitive regions of the allopolyploid wheat genome.

Results

Although open chromatin is found to be significantly enriched around genes, the majority of MNase-sensitive regions are located within transposable elements (TEs). Chromatin of the smaller D genome is more accessible than that of the larger A and B genomes. Chromatin states of different TEs vary among families and are influenced by the TEs’ chromosomal position and proximity to genes. While the chromatin accessibility of genes is influenced by proximity to TEs, and not by their position on the chromosomes, we observe a negative chromatin accessibility gradient along the telomere-centromere axis in the intergenic regions, positively correlated with the distance between genes. Both gene expression levels and homoeologous gene expression bias are correlated with chromatin accessibility in promoter regions. The differential nuclease sensitivity assay accurately predicts previously detected centromere locations. SNPs located within more accessible chromatin explain a higher proportion of genetic variance for a number of agronomic traits than SNPs located within more closed chromatin.

Conclusions

Chromatin states in the wheat genome are shaped by the interplay of repetitive and gene-encoding regions that are predictive of the functional and structural organization of chromosomes, providing a powerful framework for detecting genomic features involved in gene regulation and prioritizing genomic variation to explain phenotypes.

Systematic Analysis of Differential H3K27me3 and H3K4me3 Deposition in Callus and Seedling Reveals the Epigenetic Regulatory Mechanisms Involved in Callus Formation in Rice

Plant growth and development occurs through meristematic cell activity, and cell fate transition is accompanied by epigenetic modifications. Callus with cell pluripotency exhibits the ability to undergo continued cell division, and is ideal for studying plant meristematic differentiation. By comparing the differential epigenetic modifications between callus and seedling, the changes in chromatin state and effects of various epigenetic modifications on the growth and development of plants can be revealed, and the key genes related to plant growth and development can be identified, providing novel insights into the regulation of plant growth and development. In this study, we performed ChIP assays using various antibodies in rice seed-induced callus and seedlings grown for about 15 days to examine the differential deposition of H3K27me3 and H3K4me3. Furthermore, data for DNase I-hypersensitive sites in the corresponding tissues were downloaded from National Center for Biotechnology Information. We analyzed 4,562 callus H3K27me3-decreased genes especially those encoding transcription factors in callus, and found that most of the transcription factors, including AP2-ERREBP, NAC, and HB gene families, were related to growth and development. Genes related to meristemization, such as OsWOX9, OsWOX11, OsPLT4, OsPLT5, and OsSHR, were also included. In contrast, H3K4me3 positively regulated callus characteristics through its higher deposition in the callus than in the seedling. We further performed transcriptomic analysis on 45 sets of Affymetrix GeneChip arrays and identified 1,565 genes preferentially expressed in the callus. Callus development and root development in rice were found to share a common regulatory mechanism. We found that these genes, which are associated with meristems, require the removal of H3K27me3 and the deposition of H3K4me3, and DNase I-hypersensitive sites to maintain a relatively active state in the callus than in the seedling. The present study provides novel data about the epigenetic mechanisms involved in callus formation and additional resources for the study of cell division and differentiation in plants.

Maize E2F transcription factors. Expression, association to promoters of S-phase genes and interaction with the RBR1 protein in chromatin during seed germination

For seed germination, it is necessary to restart the cell cycle, a process regulated at multiple levels including transcriptional control, that is executed by the E2F family of transcription factors. We identified 12 genes of the E2F family in maize that are expressed differentially during the first 28 h post imbibition (HAI). E2Fa/b1;1 and E2Fc proteins were characterized as an activator and a putative repressor respectively, both forming heterodimers with DPb2 that bind differentially to consensus E2F response elements in promoters of E2F target genes. Transcripts of target genes for these transcription factors accumulate during germination; in dry seeds E2Fc protein is enriched in the target promoters and is replaced by E2Fa/b1;1 as germination advances. RBR1 is found in the same promoters in non-imbibed and 28 HAI seeds, when DNA replication has concluded, and transcription of the E2F targets should stop. During germination promoters of these target genes seem to be decorated with histone marks related to relaxed chromatin structure. Therefore, E2Fs appear to occupy their target genes in a context of open chromatin, with RBR1 fine tuning the progression between the phases.

Changes and Associations of Genomic Transcription and Histone Methylation with Salt Stress in Castor Bean

Soil salinity is a major source of abiotic plant stress, adversely affecting plant growth, development and productivity. Although the physiological and molecular mechanisms that underlie plant responses to salt stress are becoming increasingly understood, epigenetic modifications, such as histone methylations and their potential regulation of the transcription of masked genes at the genome level in response to salt stress, remain largely unclear. Castor bean, an important nonedible oil crop, has evolved the capacity to grow under salt stress. Here, based on high-throughput RNA-seq and ChIP-seq data, we systematically investigated changes in genomic transcription and histone methylation using typical histone H3 lysine 4 trimethylation (H3K4me3) and histone H3 tri-methylated lysine 27 (H3K27me3) markers in castor bean leaves subjected to salt stress. The results showed that gain or loss of histone methylation was closely associated with activated or repressed gene expression, though variations in both transcriptome and histone methylation modifications were relatively narrow in response to salt stress. Diverse salt responsive genes and switched histone methylation sites were identified in this study. In particular, we found for the first time that the transcription of the key salt-response regulator RADIALIS-LIKE SANT (RSM1), a MYB-related transcription factor involved in ABA(abscisic acid)-mediated salt stress signaling, was potentially regulated by bivalent H3K4me3-H3K27me3 modifications. Combining phenotypic variations with transcriptional and epigenetic changes, we provide a comprehensive profile for understanding histone modification, genomic transcription and their associations in response to salt stress in plants.

Targeted reprogramming of H3K27me3 resets epigenetic memory in plant paternal chromatin

Epigenetic marks are reprogrammed in the gametes to reset genomic potential in the next generation. In mammals, paternal chromatin is extensively reprogrammed through the global erasure of DNA methylation and the exchange of histones with protamines1,2. Precisely how the paternal epigenome is reprogrammed in flowering plants has remained unclear since DNA is not demethylated and histones are retained in sperm3,4. Here, we describe a multi-layered mechanism by which H3K27me3 is globally lost from histone-based sperm chromatin in Arabidopsis. This mechanism involves the silencing of H3K27me3 writers, activity of H3K27me3 erasers and deposition of a sperm-specific histone, H3.10 (ref. 5), which we show is immune to lysine 27 methylation. The loss of H3K27me3 facilitates the transcription of genes essential for spermatogenesis and pre-configures sperm with a chromatin state that forecasts gene expression in the next generation. Thus, plants have evolved a specific mechanism to simultaneously differentiate male gametes and reprogram the paternal epigenome.

Epigenomic landscape and epigenetic regulation in maize

Epigenetic regulation has been implicated in the control of multiple agronomic traits in maize. Here, we review current advances in our understanding of epigenetic regulation, which has great potential for improving agronomic traits and the environmental adaptability of crops. Epigenetic regulation plays vital role in the control of complex agronomic traits. Epigenetic variation could contribute to phenotypic diversity and can be used to improve the quality and productivity of crops. Maize (Zea mays L.), one of the most widely cultivated crops for human food, animal feed, and ethanol biofuel, is a model plant for genetic studies. Recent advances in high-throughput sequencing technology have made possible the study of epigenetic regulation in maize on a genome-wide scale. In this review, we discuss recent epigenetic studies in maize many achieved by Chinese research groups. These studies have explored the roles of DNA methylation, posttranslational modifications of histones, chromatin remodeling, and noncoding RNAs in the regulation of gene expression in plant development and environment response. We also provide our future prospects for manipulating epigenetic regulation to improve crops.

Epigenetic changes and their relationship to somaclonal variation: a need to monitor the micropropagation of plantation crops

Chromatin modulation plays important roles in gene expression regulation and genome activities. In plants, epigenetic changes, including variations in histone modification and DNA methylation, are linked to alterations in gene expression. Despite the significance and potential of in vitro cell and tissue culture systems in fundamental research and marketable applications, these systems threaten the genetic and epigenetic networks of intact plant organs and tissues. Cell and tissue culture applications can lead to DNA variations, methylation alterations, transposon activation, and finally, somaclonal variations. In this review, we discuss the status of the current understanding of epigenomic changes that occur under in vitro conditions in plantation crops, including coconut, oil palm, rubber, cotton, coffee and tea. It is hoped that comprehensive knowledge of the molecular basis of these epigenomic variations will help researchers develop strategies to enhance the totipotent and embryogenic capabilities of tissue culture systems for plantation crops.

The histone modification H3 lysine 27 tri-methylation has conserved gene regulatory roles in the triplicated genome of Brassica rapa L

Brassica rapa L. is an important vegetable and oilseed crop. We investigated the distribution of the histone mark tri-methylation of H3K27 (H3K27me3) in B. rapa and its role in the control of gene expression at two stages of development (2-day cotyledons and 14-day leaves) and among paralogs in the triplicated genome. H3K27me3 has a similar distribution in two inbred lines, while there was variation of H3K27me3 sites between tissues. Sites that are specific to 2-day cotyledons have increased transcriptional activity, and low levels of H3K27me3 in the gene body region. In 14-day leaves, levels of H3K27me3 were associated with decreased gene expression. In the triplicated genome, H3K27me3 is associated with paralogs that have tissue-specific expression. Even though B. rapa and Arabidopsis thaliana are not closely related within the Brassicaceae, there is conservation of H3K27me3-marked sites in the two species. Both B. rapa and A. thaliana require vernalization for floral initiation with FLC being the major controlling locus. In all four BrFLC paralogs, low-temperature treatment increases H3K27me3 at the proximal nucleation site reducing BrFLC expression. Following return to normal temperature growth conditions, H3K27me3 spreads along all four BrFLC paralogs providing stable repression of the gene.

Genomic imprinting in plants-revisiting existing models

In this review, Batista and Köhler revisit the current models explaining imprinting regulation in plants, and discuss novel regulatory mechanisms that could function independently of parental DNA methylation asymmetries in the establishment of imprinting.

Genomic imprinting is an epigenetic phenomenon leading to parentally biased gene expression. Throughout the years, extensive efforts have been made to characterize the epigenetic marks underlying imprinting in animals and plants. As a result, DNA methylation asymmetries between parental genomes emerged as the primary factor controlling the imprinting status of many genes. Nevertheless, the data accumulated so far suggest that this process cannot solely explain the imprinting of all genes. In this review, we revisit the current models explaining imprinting regulation in plants, and discuss novel regulatory mechanisms that could function independently of parental DNA methylation asymmetries in the establishment of imprinting.

Epigenetics Regulates Reproductive Development in Plants

Seed, resulting from reproductive development, is the main nutrient source for human beings, and reproduction has been intensively studied through genetic, molecular, and epigenetic approaches. However, how different epigenetic pathways crosstalk and integrate to regulate seed development remains unknown. Here, we review the recent progress of epigenetic changes that affect chromatin structure, such as DNA methylation, polycomb group proteins, histone modifications, and small RNA pathways in regulating plant reproduction. In gametogenesis of flowering plants, epigenetics is dynamic between the companion cell and gametes. Cytosine DNA methylation occurs in CG, CHG, CHH contexts (H = A, C, or T) of genes and transposable elements, and undergoes dynamic changes during reproduction. Cytosine methylation in the CHH context increases significantly during embryogenesis, reaches the highest levels in mature embryos, and decreases as the seed germinates. Polycomb group proteins are important transcriptional regulators during seed development. Histone modifications and small RNA pathways add another layer of complexity in regulating seed development. In summary, multiple epigenetic pathways are pivotal in regulating seed development. It remains to be elucidated how these epigenetic pathways interplay to affect dynamic chromatin structure and control reproduction.

Genome-wide analysis of the H3K27me3 epigenome and transcriptome in Brassica rapa

BACKGROUND

Genome-wide maps of histone modifications have been obtained for several plant species. However, most studies focus on model systems and do not enforce FAIR data management principles. Here we study the H3K27me3 epigenome and associated transcriptome of Brassica rapa, an important vegetable cultivated worldwide.

FINDINGS

We performed H3K27me3 chromatin immunoprecipitation followed by high-throughput sequencing and transcriptomic analysis by 3'-end RNA sequencing from B. rapa leaves and inflorescences. To analyze these data we developed a Reproducible Epigenomic Analysis pipeline using Galaxy and Jupyter, packaged into Docker images to facilitate transparency and reuse. We found that H3K27me3 covers roughly one-third of all B. rapa protein-coding genes and its presence correlates with low transcript levels. The comparative analysis between leaves and inflorescences suggested that the expression of various floral regulatory genes during development depends on H3K27me3. To demonstrate the importance of H3K27me3 for B. rapa development, we characterized a mutant line deficient in the H3K27 methyltransferase activity. We found that braA.clf mutant plants presented pleiotropic alterations, e.g., curly leaves due to increased expression and reduced H3K27me3 levels at AGAMOUS-like loci.

CONCLUSIONS

We characterized the epigenetic mark H3K27me3 at genome-wide levels and provide genetic evidence for its relevance in B. rapa development. Our work reveals the epigenomic landscape of H3K27me3 in B. rapa and provides novel genomics datasets and bioinformatics analytical resources. We anticipate that this work will lead the way to further epigenomic studies in the complex genome of Brassica crops.

Updating and interaction of polycomb repressive complex 2 components in maize (Zea mays)

The information on core components in maize polycomb repressive complex 2 (PRC2) are updated at a genome-wide scale, and the protein-protein interaction networks of PRC2 components are further provided in maize. The evolutionarily conserved polycomb group (PcG) proteins form multi-subunits polycomb repressive complexes (PRCs) that repress gene expression via chromatin condensation. In Arabidopsis, three distinct PRC2s have been identified, each determining a specific developmental program with partly functional redundancy. However, the core components and biological functions of PRC2 in cereals remain obscure. Here, we updated the information on maize PRC2 components at a genome-wide scale. Maize PRC2 subunits are highly duplicated, with five MSI1, three E(z), two ESC and two Su(z)12 homologs. ZmFIE1 is preferentially expressed in the endosperm, whereas the remaining are broadly expressed in many tissues. ZmCLF/MEZ1 and ZmFIE1 are maternally expressed imprinted genes, in contrast to the paternal-dominantly expression of ZmFIE2 in the endosperm. In maize, E(z) members likely provide a scaffold for assembling PRC2 complexes, whereas Su(z)12 and p55/MSI1-like proteins together reinforce the complex; ESC members probably determine its specificity: FIE1-PRC2 regulates endosperm cell development, whereas FIE2-PRC2 controls other cell types. The duplicated Brassicaceae-specific MEA and FIS2 also directly interact with maize PRC2 members. Together, this study establishes a roadmap for protein-protein interactions of maize PRC2 components, providing new insights into their functions in the growth and development of cereals.

WOX11 recruits a histone H3K27me3 demethylase to promote gene expression during shoot development in rice

WUSCHEL-related homeobox (WOX) genes are key regulators of meristem activity and plant development, the chromatin mechanism of which to reprogram gene expression remains unclear. Histone H3K27me3 is a chromatin mark of developmentally repressed genes. How the repressive mark is removed from specific genes during plant development is largely unknown. Here, we show that WOX11 interacts with the H3K27me3 demethylase JMJ705 to activate gene expression during shoot development in rice. Genetic analysis indicates that WOX11 and JMJ705 cooperatively control shoot growth and commonly regulate the expression of a set of genes involved in meristem identity, chloroplast biogenesis, and energy metabolism in the shoot apex. Loss of WOX11 led to increased H3K27me3 and overexpression of JMJ705 decreased the methylation levels at a subset of common targets. JMJ705 is associated with most of the WOX11-binding sites found in the tested common targets in vivo, regardless of presence or absence of the JMJ705-binding motif. Furthermore, wox11 mutation reduced JMJ705-binding to many targets genome-wide. The results suggest that recruitment of JMJ705 to specific developmental pathway genes is promoted by DNA-binding transcription factors and that WOX11 functions to stimulate shoot growth through epigenetic reprogramming of genes involved in meristem development and energy-generating pathways.

Genome-wide Analysis of Transcriptional Variability in a Large Maize-Teosinte Population

Gene expression regulation plays an important role in controlling plant phenotypes and adaptation. Here, we report a comprehensive assessment of gene expression variation through the transcriptome analyses of a large maize-teosinte experimental population. Genome-wide mapping identified 25 660 expression quantitative trait loci (eQTL) for 17 311 genes, capturing an unprecedented range of expression variation. We found that local eQTL were more frequently mapped to adjacent genes, displaying a mode of expression piggybacking, which consequently created co-regulated gene clusters. Genes within the co-regulated gene clusters tend to have relevant functions and shared chromatin modifications. Distant eQTL formed 125 significant distant eQTL hotspots with their targets significantly enriched in specific functional categories. By integrating different sources of information, we identified putative trans- regulators for a variety of metabolic pathways. We demonstrated that the bHLH transcription factor R1 and hexokinase HEX9 might act as crucial regulators for flavonoid biosynthesis and glycolysis, respectively. Moreover, we showed that domestication or improvement has significantly affected global gene expression, with many genes targeted by selection. Of particular interest, the Bx genes for benzoxazinoid biosynthesis may have undergone coordinated cis-regulatory divergence between maize and teosinte, and a transposon insertion that inactivates Bx12 was under strong selection as maize spread into temperate environments with a distinct herbivore community.

Functional Divergence between Subgenomes and Gene Pairs after Whole Genome Duplications

Gene loss following whole genome duplication (WGD) is often biased, with one subgenome retaining more ancestral genes and the other sustaining more gene deletions. While bias toward the greater expression of gene copies on one subgenome can explain bias in gene loss, this raises the question to what drives differences in gene expression levels between subgenomes. Differences in chromatin modifications and epigenetic markers between subgenomes in several model species are now being identified, providing an explanation for bias in gene expression between subgenomes. WGDs can be classified into duplications with higher, biased gene loss and bias in gene expression between subgenomes versus those with lower, unbiased rates of gene loss and an absence of detectable bias between subgenomes; however, the originally proposed link between these two classes and whether WGD results from an allo- or autopolyploid event is inconsistent with recent data from the allopolyploid Capsella bursa-pastoris. The gene balance hypothesis can explain bias in the functional categories of genes retained following WGD, the difference in gene loss rates between unbiased and biased WGDs, and how plant genomes have avoided being overrun with genes encoding dose-sensitive subunits of multiprotein complexes. Comparisons of gene expression patterns between retained transcription factor pairs in maize suggest the high degree of retention for WGD-derived pairs of transcription factors may instead be explained by the older duplication-degeneration-complementation model.

Patterns and Consequences of Subgenome Differentiation Provide Insights into the Nature of Paleopolyploidy in Plants

Polyploidy is an important feature of plant genomes, but the nature of many polyploidization events remains to be elucidated. Here, we demonstrate that the evolutionary fates of the subgenomes in maize (Zea mays) and soybean (Glycine max) have followed different trajectories. One subgenome has been subject to relaxed selection, lower levels of gene expression, higher rates of transposable element accumulation, more small interfering RNAs and DNA methylation around genes, and higher rates of gene loss in maize, whereas none of these features were observed in soybean. Nevertheless, individual gene pairs exhibit differentiation with respect to these features in both species. In addition, we observed a higher number of chromosomal rearrangements and higher frequency of retention of duplicated genes in soybean than in maize. Furthermore, soybean "singletons" were found to be more frequently tandemly duplicated than "duplicates" in soybean, which may, to some extent, counteract the genome imbalance caused by gene loss. We propose that unlike in maize, in which two subgenomes were distinct prior to the allotetraploidization event and thus experienced global differences in selective constraints, in soybean, the two subgenomes were far less distinct prior to polyploidization, such that individual gene pairs, rather than subgenomes, experienced stochastic differences over longer periods of time, resulting in retention of the majority of duplicates.

Comparative epigenomics in the Brassicaceae reveals two evolutionarily conserved modes of PRC2-mediated gene regulation

Background

Polycomb Repressive Complexes 2 (PRC2) are multi-protein chromatin modifiers that are evolutionarily conserved among eukaryotes and play key roles in the regulation of gene expression, notably through the trimethylation of lysine 27 of histone H3 (H3K27me3). Although PRC2-mediated gene regulation has been studied in many organisms, few studies have explored in depth the evolutionary conservation of PRC2 targets.

Results

Here, we compare the H3K27me3 epigenomic profiles for the two closely related species Arabidopsis thaliana and Arabidopsis lyrata and the more distant species Arabis alpina, three Brassicaceae that diverged from each other within the past 24 million years.

Using a robust set of gene orthologs present in the three species, we identify two classes of evolutionarily conserved PRC2 targets, which are characterized by either developmentally plastic or developmentally constrained H3K27me3 marking across species. Constrained H3K27me3 marking is associated with higher conservation of promoter sequence information content and higher nucleosome occupancy compared to plastic H3K27me3 marking. Moreover, gene orthologs with constrained H3K27me3 marking exhibit a higher degree of tissue specificity and tend to be involved in developmental functions, whereas gene orthologs with plastic H3K27me3 marking preferentially encode proteins associated with metabolism and stress responses. In addition, gene orthologs with constrained H3K27me3 marking are the predominant contributors to higher-order chromosome organization.

Conclusions

Our findings indicate that developmentally plastic and constrained H3K27me3 marking define two evolutionarily conserved modes of PRC2-mediated gene regulation that are associated with distinct selective pressures operating at multiple scales, from DNA sequence to gene function and chromosome architecture.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1333-9) contains supplementary material, which is available to authorized users.

Genomic Analysis of the DNA Replication Timing Program during Mitotic S Phase in Maize (Zea mays) Root Tips

All plants and animals must replicate their DNA, using a regulated process to ensure that their genomes are completely and accurately replicated. DNA replication timing programs have been extensively studied in yeast and animal systems, but much less is known about the replication programs of plants. We report a novel adaptation of the "Repli-seq" assay for use in intact root tips of maize (Zea mays) that includes several different cell lineages and present whole-genome replication timing profiles from cells in early, mid, and late S phase of the mitotic cell cycle. Maize root tips have a complex replication timing program, including regions of distinct early, mid, and late S replication that each constitute between 20 and 24% of the genome, as well as other loci corresponding to ∼32% of the genome that exhibit replication activity in two different time windows. Analyses of genomic, transcriptional, and chromatin features of the euchromatic portion of the maize genome provide evidence for a gradient of early replicating, open chromatin that transitions gradually to less open and less transcriptionally active chromatin replicating in mid S phase. Our genomic level analysis also demonstrated that the centromere core replicates in mid S, before heavily compacted classical heterochromatin, including pericentromeres and knobs, which replicate during late S phase.

Natural antisense transcripts are significantly involved in regulation of drought stress in maize

Natural antisense transcripts (NATs) are a prominent and complex class of regulatory RNAs. Using strand-specific RNA sequencing, we identified 1769 sense and antisense transcript pairs (NAT pairs) in two maize inbreds with different sensitivity to drought, as well as in two derivative recombination inbred lines (RILs). A significantly higher proportion of NATs relative to non-NATs are specifically expressed under water stress (WS). Surprisingly, expression of sense and antisense transcripts produced by NAT pairs is significantly correlated, particularly under WS. We found an unexpected large proportion of NATs with protein coding potential, as estimated by ribosome release scores. Small RNAs significantly accumulate within NAT pairs, with 21 nt smRNA particularly enriched in overlapping regions of these pairs of genes. The abundance of these smRNAs is significantly altered in the leafbladeless1 mutant, suggesting that these genes may be regulated by the tasiRNA pathway. Further, NATs are significantly hypomethylated and include fewer transposable element sequences relative to non-NAT genes. NAT gene regions also exhibit higher levels of H3K36me3, H3K9ac, and H3K4me3, but lower levels of H3K27me3, indicating that NAT gene pairs generally exhibit an open chromatin configuration. Finally, NAT pairs in 368 diverse maize inbreds and 19 segregating populations were specifically enriched for polymorphisms associated with drought tolerance. Taken together, the data highlight the potential impact of that small RNAs and histone modifications have in regulation of NAT expression, and the significance of NATs in response to WS.

Genome-wide mapping of transcriptional enhancer candidates using DNA and chromatin features in maize

Background

While most cells in multicellular organisms carry the same genetic information, in each cell type only a subset of genes is being transcribed. Such differentiation in gene expression depends, for a large part, on the activation and repression of regulatory sequences, including transcriptional enhancers. Transcriptional enhancers can be located tens of kilobases from their target genes, but display characteristic chromatin and DNA features, allowing their identification by genome-wide profiling. Here we show that integration of chromatin characteristics can be applied to predict distal enhancer candidates in Zea mays, thereby providing a basis for a better understanding of gene regulation in this important crop plant.

Result

To predict transcriptional enhancers in the crop plant maize (Zea mays L. ssp. mays), we integrated available genome-wide DNA methylation data with newly generated maps for chromatin accessibility and histone 3 lysine 9 acetylation (H3K9ac) enrichment in young seedling and husk tissue. Approximately 1500 intergenic regions, displaying low DNA methylation, high chromatin accessibility and H3K9ac enrichment, were classified as enhancer candidates. Based on their chromatin profiles, candidate sequences can be classified into four subcategories. Tissue-specificity of enhancer candidates is defined based on the tissues in which they are identified and putative target genes are assigned based on tissue-specific expression patterns of flanking genes.

Conclusions

Our method identifies three previously identified distal enhancers in maize, validating the new set of enhancer candidates and enlarging the toolbox for the functional characterization of gene regulation in the highly repetitive maize genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1273-4) contains supplementary material, which is available to authorized users.

Spatio-temporal analysis of coding and long noncoding transcripts during maize endosperm development

The maize endosperm consists of three major compartmentalized cell types: the starchy endosperm (SE), the basal endosperm transfer cell layer (BETL), and the aleurone cell layer (AL). Differential genetic programs are activated in each cell type to construct functionally and structurally distinct cells. To compare gene expression patterns involved in maize endosperm cell differentiation, we isolated transcripts from cryo-dissected endosperm specimens enriched with BETL, AL, or SE at 8, 12, and 16 days after pollination (DAP). We performed transcriptome profiling of coding and long noncoding transcripts in the three cell types during differentiation and identified clusters of the transcripts exhibiting spatio-temporal specificities. Our analysis uncovered that the BETL at 12 DAP undergoes the most dynamic transcriptional regulation for both coding and long noncoding transcripts. In addition, our transcriptome analysis revealed spatio-temporal regulatory networks of transcription factors, imprinted genes, and loci marked with histone H3 trimethylated at lysine 27. Our study suggests that various regulatory mechanisms contribute to the genetic networks specific to the functions and structures of the cell types of the endosperm.

Integrated analysis and transcript abundance modelling of H3K4me3 and H3K27me3 in developing secondary xylem

Despite the considerable contribution of xylem development (xylogenesis) to plant biomass accumulation, its epigenetic regulation is poorly understood. Furthermore, the relative contributions of histone modifications to transcriptional regulation is not well studied in plants. We investigated the biological relevance of H3K4me3 and H3K27me3 in secondary xylem development using ChIP-seq and their association with transcript levels among other histone modifications in woody and herbaceous models. In developing secondary xylem of the woody model Eucalyptus grandis, H3K4me3 and H3K27me3 genomic spans were distinctly associated with xylogenesis-related processes, with (late) lignification pathways enriched for putative bivalent domains, but not early secondary cell wall polysaccharide deposition. H3K27me3-occupied genes, of which 753 (~31%) are novel targets, were enriched for transcriptional regulation and flower development and had significant preferential expression in roots. Linear regression models of the ChIP-seq profiles predicted ~50% of transcript abundance measured with strand-specific RNA-seq, confirmed in a parallel analysis in Arabidopsis where integration of seven additional histone modifications each contributed smaller proportions of unique information to the predictive models. This study uncovers the biological importance of histone modification antagonism and genomic span in xylogenesis and quantifies for the first time the relative correlations of histone modifications with transcript abundance in plants.

Concerted Divergence after Gene Duplication in Polycomb Repressive Complexes

Duplicated genes are a major contributor to genome evolution and phenotypic novelty. There are multiple possible evolutionary fates of duplicated genes. Here, we provide an example of concerted divergence of simultaneously duplicated genes whose products function in the same complex. We studied POLYCOMB REPRESSIVE COMPLEX2 (PRC2) in Brassicaceae. The VERNALIZATION (VRN)-PRC2 complex contains VRN2 and SWINGER (SWN), and both genes were duplicated during a whole-genome duplication to generate FERTILIZATION INDEPENDENT SEED2 (FIS2) and MEDEA (MEA), which function in the Brassicaceae-specific FIS-PRC2 complex that regulates seed development. We examined the expression of FIS2, MEA, and their paralogs, compared their cytosine and histone methylation patterns, and analyzed the sequence evolution of the genes. We found that FIS2 and MEA have reproductive-specific expression patterns that are correlated and derived from the broadly expressed VRN2 and SWN in outgroup species. In vegetative tissues of Arabidopsis (Arabidopsis thaliana), repressive methylation marks are enriched in FIS2 and MEA, whereas active marks are associated with their paralogs. We detected comparable accelerated amino acid substitution rates in FIS2 and MEA but not in their paralogs. We also show divergence patterns of the PRC2-associated VERNALIZATION5/VIN3-LIKE2 that are similar to FIS2 and MEA These lines of evidence indicate that FIS2 and MEA have diverged in concert, resulting in functional divergence of the PRC2 complexes in Brassicaceae. This type of concerted divergence is a previously unreported fate of duplicated genes. In addition, the Brassicaceae-specific FIS-PRC2 complex modified the regulatory pathways in female gametophyte and seed development.

Epigenetic processes in flowering plant reproduction

Seeds provide up to 70% of the energy intake of the human population, emphasizing the relevance of understanding the genetic and epigenetic mechanisms controlling seed formation. In flowering plants, seeds are the product of a double fertilization event, leading to the formation of the embryo and the endosperm surrounded by maternal tissues. Analogous to mammals, plants undergo extensive epigenetic reprogramming during both gamete formation and early seed development, a process that is supposed to be required to enforce silencing of transposable elements and thus to maintain genome stability. Global changes of DNA methylation, histone modifications, and small RNAs are closely associated with epigenome programming during plant reproduction. Here, we review current knowledge on chromatin changes occurring during sporogenesis and gametogenesis, as well as early seed development in major flowering plant models.

Cooperation between the H3K27me3 Chromatin Mark and Non-CG Methylation in Epigenetic Regulation

H3K27me3 is a repressive chromatin mark of genes and is catalyzed by homologs of Enhancer of zeste [E(z)], a component of Polycomb-repressive complex 2 (PRC2), while DNA methylation that occurs in CG and non-CG (CHG and CHH, where H is A, C, or T) contexts is a hallmark of transposon silencing in plants. However, the relationship between H3K27me3 and DNA methylation in gene repression remains unclear. In addition, the mechanism of PRC2 recruitment to specific genes is not known in plants. Here, we show that SDG711, a rice (Oryza sativa) E(z) homolog, is required to maintain H3K27me3 of many developmental genes after shoot meristem to leaf transition and that many H3K27me3-marked developmental genes are also methylated at non-CG sites in the body regions. SDG711-binding and SDG711-mediated ectopic H3K27me3 also target genes methylated at non-CG sites. Conversely, mutation of OsDRM2, a major rice CHH methyltransferase, resulted in loss of SDG711-binding and H3K27me3 from many genes and their de-repression. Furthermore, we show that SDG711 physically interacts with OsDRM2 and a putative CHG methylation-binding protein. These results together suggest that the repression of many developmental genes may involve both DRM2-mediated non-CG methylation and PRC2-mediated H3K27me3 and that the two marks are not generally mutually exclusive but may cooperate in repression of developmentally regulated genes in rice.

LHP1 Regulates H3K27me3 Spreading and Shapes the Three-Dimensional Conformation of the Arabidopsis Genome

Precise expression patterns of genes in time and space are essential for proper development of multicellular organisms. Dynamic chromatin conformation and spatial organization of the genome constitute a major step in this regulation to modulate developmental outputs. Polycomb repressive complexes (PRCs) mediate stable or flexible gene repression in response to internal and environmental cues. In Arabidopsis thaliana, LHP1 co-localizes with H3K27me3 epigenetic marks throughout the genome and interacts with PRC1 and PRC2 members as well as with a long noncoding RNA. Here, we show that LHP1 is responsible for the spreading of H3K27me3 towards the 3' end of the gene body. We also identified a subset of LHP1-activated genes and demonstrated that LHP1 shapes local chromatin topology in order to control transcriptional co-regulation. Our work reveals a general role of LHP1 from local to higher conformation levels of chromatin configuration to determine its accessibility to define gene expression patterns.

Parental epigenetic asymmetry of PRC2-mediated histone modifications in the Arabidopsis endosperm

Parental genomes in the endosperm are marked by differential DNA methylation and are therefore epigenetically distinct. This epigenetic asymmetry is established in the gametes and maintained after fertilization by unknown mechanisms. In this manuscript, we have addressed the key question whether parentally inherited differential DNA methylation affects de novo targeting of chromatin modifiers in the early endosperm. Our data reveal that polycomb-mediated H3 lysine 27 trimethylation (H3K27me3) is preferentially localized to regions that are targeted by the DNA glycosylase DEMETER (DME), mechanistically linking DNA hypomethylation to imprinted gene expression. Our data furthermore suggest an absence of de novo DNA methylation in the early endosperm, providing an explanation how DME-mediated hypomethylation of the maternal genome is maintained after fertilization. Lastly, we show that paternal-specific H3K27me3-marked regions are located at pericentromeric regions, suggesting that H3K27me3 and DNA methylation are not necessarily exclusive marks at pericentromeric regions in the endosperm.

Evolution and function of genomic imprinting in plants

Genomic imprinting, an inherently epigenetic phenomenon defined by parent of origin-dependent gene expression, is observed in mammals and flowering plants. Genome-scale surveys of imprinted expression and the underlying differential epigenetic marks have led to the discovery of hundreds of imprinted plant genes and confirmed DNA and histone methylation as key regulators of plant imprinting. However, the biological roles of the vast majority of imprinted plant genes are unknown, and the evolutionary forces shaping plant imprinting remain rather opaque. Here, we review the mechanisms of plant genomic imprinting and discuss theories of imprinting evolution and biological significance in light of recent findings.