Genome-wide analysis of histone modifications: H3K4me2, H3K4me3, H3K9ac, and H3K27ac in Oryza sativa L. Japonica

Chromatin State-Based Analysis of Epigenetic H3K4me3 Marks of Arabidopsis in Response to Dark Stress

Light is essential to plant growth and development. Extended darkness causes dramatic gene expression changes, leading to leaf senescence, hypocotyl growth, petiole elongation, reduced leaf area, and early flowering, etc. However, the underlying mechanism of response to darkness at epigenetic levels remains largely unknown. In this study, we conducted ChIP-seq to generate global epigenomic profiles of H3K4me3 under 3-day extended darkness and normal light conditions in Arabidopsis. We applied chromatin state analysis together with self-organization mapping (SOM) to study the combination of epigenetic regulation under dark stress. The SOM map clusters the segments on the genome according to multiple diverse epigenomic datasets, which breaks the limit of dispersed distribution of epigenetic marks on the genome. Through SOM analysis, we also found that the signals of H3K4me3 were mainly increased after darkness. Analysis of H3K4me3-changed genes together with differentially expressed genes indicated that the genes showing dark-increased H3K4me3 were most involved in senescence and autophagy, and cross-talk existed between dark-induced and natural senescence. In summary, we studied the regulation of the epigenetic H3K4me3 marks of Arabidopsis in response to dark stress using chromatin state and SOM analyses. Our study revealed the regulatory mechanisms of the epigenome in response to dark stress, and SOM analysis based on chromatin states used in our study will also be helpful for other studies on dynamic changes of multiple epigenetic marks.

Genome-wide analysis of three histone marks and gene expression in Paulownia fortunei with phytoplasma infection

BACKGROUND

Paulownia withes'-broom (PaWB) disease caused by phytoplasma is a serious infectious disease for Paulownia. However, the underlying molecular pathogenesis is not fully understood. Recent studies have demonstrated that histone modifications could play a role in plant defense responses to pathogens. But there is still no available genome-wide histone modification data in non-model ligneous species infected with phytoplasma.

RESULTS

Here, we provided the first genome-wide profiles of three histone marks (H3K4me3, H3K36me3 and H3K9ac) in Paulownia fortunei under phytoplasma stress by using chromatin immunoprecipitation sequencing (ChIP-Seq). We found that H3K4me3, H3K36me3 and H3K9ac were mainly enriched in the genic regions in P. fortunei with (PFI) and without (PF) phytoplasma infection. ChIP-Seq analysis revealed 1738, 986, and 2577 genes were differentially modified by H3K4me3, H3K36me3 and H3K9ac marks in PFI under phytoplasma infection, respectively. The functional analysis of these genes suggested that most of them were mainly involved in metabolic pathways, biosynthesis of secondary metabolites, phenylpropanoid biosynthesis, plant-pathogen interaction and plant hormone signal transduction. In addition, the combinational analysis of ChIP-Seq and RNA-Seq showed that differential histone methylation and acetylation only affected a small subset of phytoplasma-responsive genes.

CONCLUSIONS

Taken together, this is the first report of integrated analysis of histone modifications and gene expression involved in Paulownia-phytoplasma interaction. Our results will provide the valuable resources for the mechanism studies of gene regulation in non-model plants upon pathogens attack.

Global Proteome Analysis Links Lysine Acetylation to Diverse Functions in Oryza Sativa

Lysine acetylation (Kac) is an important protein post-translational modification in both eukaryotes and prokaryotes. Herein, we report the results of a global proteome analysis of Kac and its diverse functions in rice (Oryza sativa). We identified 1353 Kac sites in 866 proteins in rice seedlings. A total of 11 Kac motifs are conserved, and 45% of the identified proteins are localized to the chloroplast. Among all acetylated proteins, 38 Kac sites are combined in core histones. Bioinformatics analysis revealed that Kac occurs on a diverse range of proteins involved in a wide variety of biological processes, especially photosynthesis. Protein-protein interaction networks of the identified proteins provided further evidence that Kac contributes to a wide range of regulatory functions. Furthermore, we demonstrated that the acetylation level of histone H3 (lysine 27 and 36) is increased in response to cold stress. In summary, our approach comprehensively profiles the regulatory roles of Kac in the growth and development of rice.

BABA-Primed Histone Modifications in Potato for Intergenerational Resistance to Phytophthora infestans

In this paper we analyzed β-aminobutyric acid (BABA)-primed epigenetic adjustment of potato cv. "Sarpo Mira" to Phytophthora infestans. The first stress-free generation of the potato genotype obtained from BABA-primed parent plants via tubers and seeds showed pronounced resistance to the pathogen, which was tuned with the transcriptional memory of SA-responsive genes. During the early priming phase before the triggering stress, we found robust bistable deposition of histone marks (H3K4me2 and H3K27me3) on the NPR1 (Non-expressor of PR genes) and the SNI1 gene (Suppressor of NPR1, Inducible), in which transcription antagonized silencing. Switchable chromatin states of these adverse systemic acquired resistance (SAR) regulators probably reprogrammed responsiveness of the PR1 and PR2 genes and contributed to stress imprinting. The elevated levels of heritable H3K4me2 tag in the absence of transcription on SA-dependent genes in BABA-primed (F0) and its vegetative and generative progeny (F1) before pathogen challenge provided evidence for the epigenetic mark for intergenerational memory in potato. Moreover, our study revealed that histone acetylation was not critical for maintaining BABA-primed defense information until the plants were triggered with the virulent pathogen when rapid and boosted PRs gene expression probably required histone acetyltransferase (HAT) activity both in F0 and F1 progeny.

Overexpressing lncRNA LAIR increases grain yield and regulates neighbouring gene cluster expression in rice

Long noncoding RNAs (lncRNAs) are essential regulators of gene expression in eukaryotes. Despite increasing knowledge on the function of lncRNAs, little is known about their effects on crop yield. Here, we identify a lncRNA transcribed from the antisense strand of neighbouring gene LRK (leucine-rich repeat receptor kinase) cluster named LAIR (LRK Antisense Intergenic RNA). LAIR overexpression increases rice grain yield and upregulates the expression of several LRK genes. Additionally, chromatin immunoprecipitation assay results indicate H3K4me3 and H4K16ac are significantly enriched at the activated LRK1 genomic region. LAIR binds histone modification proteins OsMOF and OsWDR5 in rice cells, which are enriched in LRK1 gene region. Moreover, LAIR is demonstrated to bind 5' and 3' untranslated regions of LRK1 gene. Overall, this study reveals the role of lncRNA LAIR in regulating rice grain yield and lncRNAs may be useful targets for crop breeding.

Global Involvement of Lysine Crotonylation in Protein Modification and Transcription Regulation in Rice

Lysine crotonylation (Kcr) is a newly discovered posttranslational modification (PTM) existing in mammals. A global crotonylome analysis was undertaken in rice (Oryza sativa L. japonica) using high accuracy nano-LC-MS/MS in combination with crotonylated peptide enrichment. A total of 1,265 lysine crotonylation sites were identified on 690 proteins in rice seedlings. Subcellular localization analysis revealed that 51% of the crotonylated proteins identified were localized in chloroplasts. The photosynthesis-associated proteins were also mostly enriched in total crotonylated proteins. In addition, a genomic localization analysis of histone Kcr by ChIP-seq was performed to assess the relevance between histone Kcr and the genome. Of the 10,923 identified peak regions, the majority (86.7%) of the enriched peaks were located in gene body, especially exons. Furthermore, the degree of histone Kcr modification was positively correlated with gene expression in genic regions. Compared with other published histone modification data, the Kcr was co-located with the active histone modifications. Interestingly, histone Kcr-facilitated expression of genes with existing active histone modifications. In addition, 77% of histone Kcr modifications overlapped with DNase hypersensitive sites (DHSs) in intergenic regions of the rice genome and might mark other cis-regulatory DNA elements that are different from IPA1, a transcription activator in rice seedlings. Overall, our results provide a comprehensive understanding of the biological functions of the crotonylome and new active histone modification in transcriptional regulation in plants.

Gene expression and promoter characterization of heat-shock protein 90B gene (HSP90B) in the model unicellular green alga Chlamydomonas reinhardtii

Molecular chaperones or heat shock proteins are a large protein family with important functions in every cellular organism. Among all types of the heat shock proteins, information on the ER-localized HSP90 protein (HSP90B) and its encoding gene is relatively scarce in the literature, especially in photosynthetic organisms. In this study, expression profiles as well as promoter sequence of the HSP90B gene were investigated in the model green alga Chlamydomonas reinhardtii. We have found that HSP90B is strongly induced by heat and ER stresses, while other short-term exposure to abiotic stresses, such as salinity, dark-to-light transition or light stress does not appear to affect the expression. Promoter truncation analysis as well as chromatin immunoprecipitation using the antibodies recognizing histone H3 and acetylated histone H3, revealed a putative core constitutive promoter sequence between -1 to -253 bp from the transcription start site. Our results also suggested that the nucleotides upstream of the core promoter may contain repressive elements such as putative repressor binding site(s).

Rice Functional Genomics Research: Past Decade and Future

Rice (Oryza sativa) is a major staple food crop for more than 3.5 billion people worldwide. Understanding the regulatory mechanisms of complex agronomic traits in rice is critical for global food security. Rice is also a model plant for genomics research of monocotyledons. Thanks to the rapid development of functional genomic technologies, over 2000 genes controlling important agronomic traits have been cloned, and their molecular biological mechanisms have also been partially characterized. Here, we briefly review the advances in rice functional genomics research during the past 10 years, including a summary of functional genomics platforms, genes and molecular regulatory networks that regulate important agronomic traits, and newly developed tools for gene identification. These achievements made in functional genomics research will greatly facilitate the development of green super rice. We also discuss future challenges and prospects of rice functional genomics research.

ChIP-ping the branches of the tree: functional genomics and the evolution of eukaryotic gene regulation

Advances in the methods for detecting protein-DNA interactions have played a key role in determining the directions of research into the mechanisms of transcriptional regulation. The most recent major technological transformation happened a decade ago, with the move from using tiling arrays [chromatin immunoprecipitation (ChIP)-on-Chip] to high-throughput sequencing (ChIP-seq) as a readout for ChIP assays. In addition to the numerous other ways in which it is superior to arrays, by eliminating the need to design and manufacture them, sequencing also opened the door to carrying out comparative analyses of genome-wide transcription factor occupancy across species and studying chromatin biology in previously less accessible model and nonmodel organisms, thus allowing us to understand the evolution and diversity of regulatory mechanisms in unprecedented detail. Here, we review the biological insights obtained from such studies in recent years and discuss anticipated future developments in the field.

Relationship between histone modifications and transcription factor binding is protein family specific

The very small fraction of putative binding sites (BSs) that are occupied by transcription factors (TFs) in vivo can be highly variable across different cell types. This observation has been partly attributed to changes in chromatin accessibility and histone modification (HM) patterns surrounding BSs. Previous studies focusing on BSs within DNA regulatory regions found correlations between HM patterns and TF binding specificities. However, a mechanistic understanding of TF-DNA binding specificity determinants is still not available. The ability to predict in vivo TF binding on a genome-wide scale requires the identification of features that determine TF binding based on evolutionary relationships of DNA binding proteins. To reveal protein family-dependent mechanisms of TF binding, we conducted comprehensive comparisons of HM patterns surrounding BSs and non-BSs with exactly matched core motifs for TFs in three cell lines: 33 TFs in GM12878, 37 TFs in K562, and 18 TFs in H1-hESC. These TFs displayed protein family-specific preferences for HM patterns surrounding BSs, with high agreement among cell lines. Moreover, compared to models based on DNA sequence and shape at flanking regions of BSs, HM-augmented quantitative machine-learning methods resulted in increased performance in a TF family-specific manner. Analysis of the relative importance of features in these models indicated that TFs, displaying larger HM pattern differences between BSs and non-BSs, bound DNA in an HM-specific manner on a protein family-specific basis. We propose that TF family-specific HM preferences reveal distinct mechanisms that assist in guiding TFs to their cognate BSs by altering chromatin structure and accessibility.

Comparative analysis of Histone modifications and DNA methylation at OsBZ8 locus under salinity stress in IR64 and Nonabokra rice varieties

Rice being an important cereal crop is highly sensitive to salinity stress causing growth retardation and loss in productivity. However, certain rice genotypes like Nonabokra and Pokkali show a high level of tolerance towards salinity stress compared to IR64 variety. This differential response of tolerant varieties towards salinity stress may be a cumulative effect of genetic and epigenetic factors. In this study, we have compared the salinity-induced changes in chromatin modifications at the OsBZ8 locus in salt-tolerant Nonabokra and salt-sensitive IR64 rice varieties. Expression analysis indicates that the OsBZ8 gene is highly induced in Nonabokra plants even in the absence of salt stress, whereas in IR64, the expression significantly increases only during salt stress. Sequence analysis and nucleosomal arrangement within the region -2000 to +1000 of OsBZ8 gene show no difference between the two rice varieties. However, there was a considerable difference in histone modifications and DNA methylation at the locus between these varieties. In Nonabokra, the upstream region was hyperacetylated at H3K9 and H3K27, and this acetylation did not change during salt stress. However, in IR64, histone acetylation was observed only during salt stress. Moreover, the upstream region of OsBZ8 gene has highly dynamic nucleosome arrangement in Nonabokra, compared to IR64. Furthermore, loss of DNA methylation was observed at OsBZ8 locus in Nonabokra control plants along with low H3K27me3 and high H3K4me3. Control IR64 plants show high DNA methylation and enriched H3K27me3. Collectively these results indicate a significant difference in chromatin modifications between the rice varieties that regulates differential expression of OsBZ8 gene during salt stress.

Genome-wide comparative analysis of H3K4me3 profiles between diploid and allotetraploid cotton to refine genome annotation

Polyploidy is a common evolutionary occurrence in plants. Recently, published genomes of allotetraploid G. hirsutum and its donors G. arboreum and G. raimondii make cotton an accessible polyploid model. This study used chromatin immunoprecipitation with high-throughput sequencing (ChIP-Seq) to investigate the genome-wide distribution of H3K4me3 in G. arboreum and G. hirsutum, and explore the conservation and variation of genome structures between diploid and allotetraploid cotton. Our results showed that H3K4me3 modifications were associated with active transcription in both cottons. The H3K4me3 histone markers appeared mainly in genic regions and were enriched around the transcription start sites (TSSs) of genes. We integrated the ChIP-seq data of H3K4me3 with RNA-seq and ESTs data to refine the genic structure annotation. There were 6,773 and 12,773 new transcripts discovered in G. arboreum and G. hirsutum, respectively. Furthermore, co-expression networks were linked with histone modification and modularized in an attempt to explain differential H3K4me3 enrichment correlated with changes in gene transcription during cotton development and evolution. Taken together, we have combined epigenomic and transcriptomic datasets to systematically discover functional genes and compare them between G. arboreum and G. hirsutum, which may be beneficial for studying diploid and allotetraploid plants with large genomes and complicated evolution.

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.

ATX3, ATX4, and ATX5 Encode Putative H3K4 Methyltransferases and Are Critical for Plant Development

Methylation of Lys residues in the tail of the H3 histone is a key regulator of chromatin state and gene expression, conferred by a large family of enzymes containing an evolutionarily conserved SET domain. One of the main types of SET domain proteins are those controlling H3K4 di- and trimethylation. The genome of Arabidopsis (Arabidopsis thaliana) encodes 12 such proteins, including five ARABIDOPSIS TRITHORAX (ATX) proteins and seven ATX-Related proteins. Here, we examined three until-now-unexplored ATX proteins, ATX3, ATX4, and ATX5. We found that they exhibit similar domain structures and expression patterns and are redundantly required for vegetative and reproductive development. Concurrent disruption of the ATX3, ATX4, and ATX5 genes caused marked reduction in H3K4me2 and H3K4me3 levels genome-wide and resulted in thousands of genes expressed ectopically. Furthermore, atx3/atx4/atx5 triple mutants resulted in exaggerated phenotypes when combined with the atx2 mutant but not with atx1 Together, we conclude that ATX3, ATX4, and ATX5 are redundantly required for H3K4 di- and trimethylation at thousands of sites located across the genome, and genomic features associated with targeted regions are different from the ATXR3/SDG2-controlled sites in Arabidopsis.

The quest for epigenetic regulation underlying unisexual flower development in Cucumis melo

BACKGROUND

Melon (Cucumis melo) is an important vegetable crop from the Cucurbitaceae family and a reference model specie for sex determination, fruit ripening and vascular fluxes studies. Nevertheless, the nature and role of its epigenome in gene expression regulation and more specifically in sex determination remains largely unknown.

RESULTS

We have investigated genome wide H3K27me3 and H3K9ac histone modifications and gene expression dynamics, in five melon organs. H3K9ac and H3K27me3 were mainly distributed along gene-rich regions and constrained to gene bodies. H3K9ac was preferentially located at the TSS, whereas H3K27me3 distributed uniformly from TSS to TES. As observed in other species, H3K9ac and H3K27me3 correlated with high and low gene expression levels, respectively. Comparative analyses of unisexual flowers pointed out sex-specific epigenetic states of TFs involved in ethylene response and flower development. Chip-qPCR analysis of laser dissected carpel and stamina primordia, revealed sex-specific histone modification of MADS-box genes. Using sex transition mutants, we demonstrated that the female promoting gene, CmACS11, represses the expression of the male promoting gene CmWIP1 via deposition of H3K27me3.

CONCLUSIONS

Our findings reveal the organ-specific landscapes of H3K9ac and H3K27me3 in melon. Our results also provide evidence that the sex determination genes recruit histone modifiers to orchestrate unisexual flower development in monoecious species.

[Reciprocity between active transcription and histone methylation]

In the nucleus of eukaryotic cells, the chromatin states dictated by the different combinations of histone post-translational modifications, such as the methylation of lysine residues, are an integral part of the multitude of epigenomes involved in the fine tuning of all genome functions, and in particular transcription. Over the last decade, an increasing number of factors have been identified as regulators involved in the establishment, reading or erasure of histone methylations. Their characterization in model organisms such as Arabidopsis has thus unraveled their fundamental roles in the control and regulation of essential developmental processes such as the floral transition, cell differentiation, gametogenesis, and/or the response/adaptation of plants to environmental stresses. In this review, we will focus on the methylation of histones functioning as a mark of activate transcription and we will try to highlight, based on recent findings, the more or less direct links between this mark and gene expression. Thus, we will discuss the different mechanisms allowing the dynamics and the integration of the chromatin states resulting from the different histone methylations in connection with the transcriptional machinery of the RNA polymerase II.

Epigenetic silencing of JMJD5 promotes the proliferation of hepatocellular carcinoma cells by down-regulating the transcription of CDKN1A 686

Proteins that contain jumonji C (JmjC) domains have recently been identified as major contributors to various malignant human cancers through epigenetic remodeling. However, the roles of these family members in the pathogenesis of hepatocellular carcinoma (HCC) are obscure. By mining public databases, we found that the HCC patients with lower JmjC domain-containing protein 5 (JMJD5) expression exhibited shorter survival time. We then confirmed that JMJD5 expression was indeed decreased in HCC specimens, which was caused by the altered epigenetic histone modifications, the decreased H3K9ac, H3K27ac and H3K4me2/3 together with the increased trimethylation of H3K27 and H3K9 on the JMJD5 promoter. Functional experiments revealed that JMJD5 knockdown promoted HCC cell proliferation and in vivo tumorigenicity by accelerating the G1/S transition of the cell cycle; in contrast, ectopic JMJD5 expression had the opposite effects. At molecular mechanism, we found that, in HCC cell lines including TP53-null Hep3B, JMJD5 knockdown led to the down-regulation of CDKN1A and ectopic expression of JMJD5 not only increased but also rescued CDKN1A transcription. Moreover, CDKN1A knockdown could abrogate the effect of JMJD5 knockdown or overexpression on cell proliferation, suggesting that JMJD5 inhibits HCC cell proliferation mainly by activating CDKN1A expression. We further revealed that JMJD5 directly enhances CDKN1A transcription by binding to CDKN1A's promoter independent of H3K36me2 demethylase activity. In short, we first prove that JMJD5 is a tumor suppressor gene in HCC pathogenesis, and the epigenetic silencing of JMJD5 promotes HCC cell proliferation by directly down-regulating CDKN1A transcription.

Epigenetic regulation and epigenomic landscape in rice

Epigenetic regulation has been implicated in the control of complex agronomic traits in rice (Oryza sativa), a staple food crop and model monocot plant. Recent advances in high-throughput sequencing and the moderately complex genome of rice have made it possible to study epigenetic regulation in rice on a genome-wide scale. This review discusses recent advances in our understanding of epigenetic regulation in rice, with an emphasis on the roles of key epigenetic regulators, the epigenomic landscape, epigenetic variation, transposon repression, and plant development.

High-resolution mapping of H4K16 and H3K23 acetylation reveals conserved and unique distribution patterns in Arabidopsis and rice

Histone acetylation and deacetylation are key epigenetic gene regulatory mechanisms that play critical roles in eukaryotes. Acetylation of histone 4 lysine 16 (H4K16ac) is implicated in many cellular processes. However, its biological function and relationship with transcription are largely unexplored in plants. We generated first genome-wide high-resolution maps of H4K16ac in Arabidopsis thaliana and Oryza sativa. We showed that H4K16ac is preferentially enriched around the transcription start sites and positively correlates with gene expression levels. Co-existence of H4K16ac and H3K23ac is correlated with high gene expression levels, suggesting a potentially combinatorial effect of H4K16ac and H3K23ac histone 3 lysine 23 acetylation on gene expression. Our data further revealed that while genes enriched with both H4K16ac and H3K23ac are ubiquitously expressed, genes enriched with only H4K16ac or H3K23ac showed significantly distinct expression patterns in association with particular developmental stages. Unexpectedly, and unlike in Arabidopsis, there are significant levels of both H4K16ac and H3K23ac in the lowly expressed genes in rice. Furthermore, we found that H4K16ac-enriched genes are associated with different biological processes in Arabidopsis and rice, suggesting a potentially species-specific role of H4K16ac in plants. Together, our genome-wide profiling reveals the conserved and unique distribution patterns of H4K16ac and H3K23ac in Arabidopsis and rice and provides a foundation for further understanding their function in plants.

Knockdown of WHIRLY1 Affects Drought Stress-Induced Leaf Senescence and Histone Modifications of the Senescence-Associated Gene HvS40

The plastid-nucleus located protein WHIRLY1 has been described as an upstream regulator of leaf senescence, binding to the promoter of senescence-associated genes like HvS40. To investigate the impact of WHIRLY1 on drought stress-induced, premature senescence, transgenic barley plants with an RNAi-mediated knockdown of the HvWHIRLY1 gene were grown under normal and drought stress conditions. The course of leaf senescence in these lines was monitored by physiological parameters and studies on the expression of senescence- and drought stress-related genes. Drought treatment accelerated leaf senescence in WT plants, whereas WHIRLY 1 knockdown lines (RNAi-W1) showed a stay-green phenotype. Expression of both senescence-associated and drought stress-responsive genes, was delayed in the transgenic plants. Notably, expression of transcription factors of the WRKY and NAC families, which are known to function in senescence- and stress-related signaling pathways, was affected in plants with impaired accumulation of WHIRLY1, indicating that WHIRLY1 acts as an upstream regulator of drought stress-induced senescence. To reveal the epigenetic indexing of HvS40 at the onset of drought-induced senescence in WT and RNAi-W1 lines, stress-responsive loading with histone modifications of promoter and coding sequences of HvS40 was analyzed by chromatin immunoprecipitation and quantified by qRT-PCR. In the wildtype, the euchromatic mark H3K9ac of the HvS40 gene was low under control conditions and was established in response to drought treatment, indicating the action of epigenetic mechanisms in response to drought stress. However, drought stress caused no significant increase in H3K9ac in plants impaired in accumulation of WHIRLY1. The results show that WHIRLY1 knockdown sets in motion a delay in senescence that involves all aspects of gene expression, including changes in chromatin structure.

Identification of Regulatory DNA Elements Using Genome-wide Mapping of DNase I Hypersensitive Sites during Tomato Fruit Development

Development and ripening of tomato fruit are precisely controlled by transcriptional regulation, which depends on the orchestrated accessibility of regulatory proteins to promoters and other cis-regulatory DNA elements. This accessibility and its effect on gene expression play a major role in defining the developmental process. To understand the regulatory mechanism and functional elements modulating morphological and anatomical changes during fruit development, we generated genome-wide high-resolution maps of DNase I hypersensitive sites (DHSs) from the fruit tissues of the tomato cultivar "Moneymaker" at 20 days post anthesis as well as break stage. By exploring variation of DHSs across fruit development stages, we pinpointed the most likely hypersensitive sites related to development-specific genes. By detecting binding motifs on DHSs of these development-specific genes or genes in the ascorbic acid biosynthetic pathway, we revealed the common regulatory elements contributing to coordinating gene transcription of plant ripening and specialized metabolic pathways. Our results contribute to a better understanding of the regulatory dynamics of genes involved in tomato fruit development and ripening.

Towards genome-wide prediction and characterization of enhancers in plants

Enhancers are important cis-regulatory DNA elements that regulate transcription programs by recruiting transcription factors and directing them to the promoters of target genes in a cell-type/tissue-specific manner. The expression of a gene can be regulated by one or multiple enhancers at different developmental stages and/or in different tissues. Enhancers are difficult to identify because of their unpredictable positions relative to their cognate promoters. Remarkably, only a handful of enhancers have been identified in plant species largely due to the lack of general approaches for enhancer identification. Extensive genomic and epigenomic research in mammalian species has revealed that the genomic locations of enhancers can be predicted based on the binding sites of transcriptional co-factors and several distinct features associated with open chromatin. Here we review the methodologies used in enhancer prediction in mammalian species. We also review the recent applications of these methodologies in Arabidopsis thaliana and discuss the future directions of enhancer identification in plants. This article is part of a Special Issue entitled: Plant Gene Regulatory Mechanisms and Networks, edited by Dr. Erich Grotewold and Dr. Nathan Springer.

The analysis of Gene Regulatory Networks in plant evo-devo

We provide an overview of methods and workflows that can be used to investigate the topologies of Gene Regulatory Networks (GRNs) in the context of plant evolutionary-developmental (evo-devo) biology. Many of the species that occupy key positions in plant phylogeny are poorly adapted as laboratory models and so we focus here on techniques that can be efficiently applied to both model and non-model species of interest to plant evo-devo. We outline methods that can be used to describe gene expression patterns and also to elucidate the transcriptional, post-transcriptional, and epigenetic regulatory mechanisms underlying these patterns, in any plant species with a sequenced genome. We furthermore describe how the technique of Protein Resurrection can be used to confirm inferences on ancestral GRNs and also to provide otherwise-inaccessible points of reference in evolutionary histories by exploiting paralogues generated in gene and whole genome duplication events. Finally, we argue for the better integration of molecular data with information from paleobotanical, paleoecological, and paleogeographical studies to provide the fullest possible picture of the processes that have shaped the evolution of plant development.

Molecular aspects of zygotic embryogenesis in sunflower (Helianthus annuus L.): correlation of positive histone marks with HaWUS expression and putative link HaWUS/HaL1L

The link HaWUS/ HaL1L , the opposite transcriptional behavior, and the decrease/increase in positive histone marks bond to both genes suggest an inhibitory effect of WUS on HaL1L in sunflower zygotic embryos. In Arabidopsis, a group of transcription factors implicated in the earliest events of embryogenesis is the WUSCHEL-RELATED HOMEOBOX (WOX) protein family including WUSCHEL (WUS) and other 14 WOX protein, some of which contain a conserved WUS-box domain in addition to the homeodomain. WUS transcripts appear very early in embryogenesis, at the 16-cell embryo stage, but gradually become restricted to the center of the developing shoot apical meristem (SAM) primordium and continues to be expressed in cells of the niche/organizing center of SAM and floral meristems to maintain stem cell population. Moreover, WUS has decisive roles in the embryonic program presumably promoting the vegetative-to-embryonic transition and/or maintaining the identity of the embryonic stem cells. However, data on the direct interaction between WUS and key genes for seed development (as LEC1 and L1L) are not collected. The novelty of this report consists in the characterization of Helianthus annuus WUS (HaWUS) gene and in its analysis regarding the pattern of the methylated lysine 4 (K4) of the Histone H3 and of the acetylated histone H3 during the zygotic embryo development. Also, a parallel investigation was performed for HaL1L gene since two copies of the WUS-binding site (WUSATA), previously identified on HaL1L nucleotide sequence, were able to be bound by the HaWUS recombinant protein suggesting a not described effect of HaWUS on HaL1L transcription.

AAL exacerbates pro-inflammatory response in macrophages by regulating Mincle/Syk/Card9 signaling along with the Nlrp3 inflammasome assembly

Previously, we demonstrated that Agrocybe aegerita lectin (AAL), a galectin isolated from edible mushroom Agrocybe aegerita, exerts potent anti-tumor activity, while the mechanisms by which AAL suppresses tumor growth are yet to be elucidated. Here, we conducted studies with focus for its impact on the cecal ligation and puncture (CLP)-induced innate immune response. Administration of AAL significantly exacerbated the severity of CLP-induced septic shock as manifested the increased lethality. AAL promoted inflammatory cytokine production by preferentially regulating macrophage activation and recruitment. Mechanistic studies revealed that AAL likely targets macrophages through receptor Mincle to activate Syk/Card9 signaling, which then couples to the Nlrp3 inflammasome assembly. It was further noted that AAL markedly promotes H3K4 di- and trimethylation, by which it enhances Hmgb1 expression. Specifically, AAL induced macrophages secretion of copious amount of Hmgb1 as manifested the Hmgb1 cytoplasmic translocation along with the detection of extracellular Hmgb1. AAL also stimulated a significant increase for nuclear Hmgb1, which then formed a complex with RelA, and thereby enhancing NF-κB transcriptional activity. Together, our data suggest that AAL may possess important pharmaceutical properties in the regulation of innate immune response.

Alterations of histone modifications at the senescence-associated gene HvS40 in barley during senescence

The barley gene HvS40, encoding a putative regulator of leaf senescence, is strongly induced during leaf senescence. As shown by chromatin immunoprecipitation, euchromatic histone modification H3K9ac is added at promoter close to ATG and coding sequence of HvS40 after onset of senescence. In parallel, level of heterochromatic H3K9me2 decreases at this gene. Bisulfite sequencing revealed no DNA-methylation in this region, but a heavily methylated DNA-island, starting 664 bp upstream from translational start site in both, mature and senescent leaves. A decrease in DNA methylation in senescing leaves could be shown at one specific CpG motif at the end of this methylation island. In addition, global changes in chromatin structure during senescence were analyzed via immunocytology, revealing senescence-associated changes in spatial distribution of heterochromatic H3K9me2 patterns in the nuclei. Our results prove a senescence-specific mechanism, altering histone modification marks at distinct sequences of the senescence-associated gene HvS40 and altering distribution of heterochromatic areas in the nuclei.

Photosynthetic Genes and Genes Associated with the C4 Trait in Maize Are Characterized by a Unique Class of Highly Regulated Histone Acetylation Peaks on Upstream Promoters

Histone modifications contribute to gene regulation in eukaryotes. We analyzed genome-wide histone H3 Lysine (Lys) 4 trimethylation and histone H3 Lys 9 acetylation (two modifications typically associated with active genes) in meristematic cells at the base and expanded cells in the blade of the maize (Zea mays) leaf. These data were compared with transcript levels of associated genes. For individual genes, regulations (fold changes) of histone modifications and transcript levels were much better correlated than absolute intensities. When focusing on regulated histone modification sites, we identified highly regulated secondary H3 Lys 9 acetylation peaks on upstream promoters (regulated secondary upstream peaks [R-SUPs]) on 10% of all genes. R-SUPs were more often found on genes that were up-regulated toward the blade than on down-regulated genes and specifically, photosynthetic genes. Among those genes, we identified six genes encoding enzymes of the C4 cycle and a significant enrichment of genes associated with the C4 trait derived from transcriptomic studies. On the DNA level, R-SUPs are frequently associated with ethylene-responsive elements. Based on these data, we suggest coevolution of epigenetic promoter elements during the establishment of C4 photosynthesis.

Genome-wide histone acetylation correlates with active transcription in maize

Gene expression is regulated at many different levels during the life cycle of all plant species. Recent investigations have taken advantage of next-generation sequencing to study the relevance of DNA methylation and sRNAs in controlling tissue-specific gene expression in maize at the genome-wide level. Here, we profiled H3K27ac in maize, which has one of the largest sequenced plant genomes due to the amplification of retrotransposons. Because transcribed genes represent only a small proportion of its genome, gene-specific epigenetic modifications are concentrated in a relatively small percentage of the genome. Indeed, H3K27ac marks are mostly in gene-rich, in contrast to gene-poor regions. A large proportion of those marks are located in transcribed regions of genes, including 111 out of 458 known genetic loci. Moreover, increased transcription correlates with the presence of H3K27ac modification in gene bodies. Using maize as an example, we suggest that H3K27ac marks actively transcribed genes in plants.

Important biological information uncovered in previously unaligned reads from chromatin immunoprecipitation experiments (ChIP-Seq)

Establishing the architecture of gene regulatory networks (GRNs) relies on chromatin immunoprecipitation followed by massively parallel sequencing (ChIP-Seq) methods that provide genome-wide transcription factor binding sites (TFBSs). ChIP-Seq furnishes millions of short reads that, after alignment, describe the genome-wide binding sites of a particular TF. However, in all organisms investigated an average of 40% of reads fail to align to the corresponding genome, with some datasets having as much as 80% of reads failing to align. We describe here the provenance of previously unaligned reads in ChIP-Seq experiments from animals and plants. We show that a substantial portion corresponds to sequences of bacterial and metazoan origin, irrespective of the ChIP-Seq chromatin source. Unforeseen was the finding that 30%–40% of unaligned reads were actually alignable. To validate these observations, we investigated the characteristics of the previously unaligned reads corresponding to TAL1, a human TF involved in lineage specification of hemopoietic cells. We show that, while unmapped ChIP-Seq read datasets contain foreign DNA sequences, additional TFBSs can be identified from the previously unaligned ChIP-Seq reads. Our results indicate that the re-evaluation of previously unaligned reads from ChIP-Seq experiments will significantly contribute to TF target identification and determination of emerging properties of GRNs.

Eukaryotic tRNA paradox

tRNAs are widely believed to segregate into two classes, I and II. Computational analysis of eukaryotic tRNA entries in Genomic tRNA Database, however, leads to new, albeit paradoxical, presence of more than a thousand class-I tRNAs with uncharacteristic long variable arms (V-arms), like in class-II. Out of 62,202 tRNAs from 69 eukaryotes, as many as 1431 class-I tRNAs have these novel extended V-arms, and we refer to them as paradoxical tRNAs (pxtRNAs). A great majority of these 1431 pxtRNA genes are located in intergenic regions, about 18% embedded in introns of genes or ESTs, and just one in 3'UTR. A check on the conservations of 2D and 3D base pairs for each position of these pxtRNAs reveals a few variations, but they seem to have almost all the known features (already known identity and conserved elements of tRNA). Analyses of the A-Box and B-Box of these pxtRNA genes in eukaryotes display salient deviations from the previously annotated conserved features of the standard promoters, whereas the transcription termination signals are just canonical and non-canonical runs of thymidine, similar to the ones in standard tRNA genes. There is just one such pxtRNA(ProAGG) gene in the entire human genome, and the availability of data allows epigenetic analysis of this human pxtRNA(ProAGG) in three different cell lines, H1 hESC, K562, and NHEK, to assess the level of its expression. Histone acetylation and methylation of this lone pxtRNA(ProAGG) gene in human differ from that of the nine standard human tRNA(ProAGG) genes. The V-arm nucleotide sequences and their secondary structures in pxtRNA differ from that of class-II tRNA. Considering these differences, hypotheses of alternative splicing, non-canonical intron and gene transfer are examined to partially improve the Cove scores of these pxtRNAs and to critically question their antecedence and novelty.

Cytomixis doesn't induce obvious changes in chromatin modifications and programmed cell death in tobacco male meiocytes

Cytomixis is a poorly studied process of nuclear migration between plant cells. It is so far unknown what drives cytomixis and what is the functional state of the chromatin migrating between cells. Using immunostaining, we have analyzed the distribution of posttranslational histone modifications (methylation, acetylation, and phosphorylation) that reflect the functional state of chromatin in the tobacco microsporocytes involved in cytomixis. We demonstrate that the chromatin in the cytomictic cells does not differ from the chromatin in intact microsporocytes according to all 14 analyzed histone modification types. We have also for the first time demonstrated that the migrating chromatin contains normal structures of the synaptonemal complex (SC) and lacks any signs of apoptosis. As has been shown, the chromatin migrating between cells in cytomixis is neither selectively heterochromatized nor degraded both before its migration to another cell and after it enters a recipient cell as micronuclei. We also showed that cytomictic chromatin contains marks typical for transcriptionally active chromatin as well as heterochromatin. Moreover, marks typical for chromosome condensation, SC formation and key proteins required for the formation of bivalents were also detected at migrated chromatin.

Identification and expression analysis of rice histone genes

Histones, acting as the core of nucleosome, are the chief protein component of chromatin. They play an important role in gene regulation by covalent modification at several sites and histone variants replacement. Five major families of histones exist: H1, H2A, H2B, H3 and H4. The protein sequences within each family appear to be highly conserved. In this paper, we identified 60 histone proteins in rice (Oryza sativa) including 14 H2A, 15 H2B, 16 H3, 11 H4 and 4 H1. Sequence analysis indicates that histone protein sequences in plant are more variable than in animal. Interestingly, we found a rice-specific H4 variant which showed several amino acid substitutions with canonical protein and was expressed in different tissues in a low level. Expression analysis indicates that a subset of histone genes were expressed in a similar pattern and many of them responded to stress conditions. Specifically, we found that two H2A.Z genes were down-regulated by stress in leaves but not in roots suggesting that they might be involved in stress response.

PNRD: a plant non-coding RNA database

The first ncRNA found was an alanine tRNA in baker's yeast, and the first detected microRNAs (miRNAs) promoted ncRNA research to a whole new level. Research on ncRNAs in animals has focused on the medical field, while in plant scientists are more concerned with improving agronomic traits. In 2010, we constructed a plant miRNA database named PMRD to meet the demand for miRNA research in plants. To provide a way to do fundamental research on plant ncRNAs and take full advantage of tremendous public resources, we designed an updated platform called plant ncRNA database (PNRD) based on its predecessor PMRD, which is accessible at http://structuralbiology.cau.edu.cn/PNRD. We collected a total of 25739 entries of 11 different types of ncRNAs from 150 plant species. Targets of miRNAs were extended to 178138 pairs in 46 species, while the number of miRNA expression profiles reached 35. Improvements in PNRD are not only the larger amounts of data, but also better service, such as a more user-friendly interface, more multifunctional and browsing options and more background data for users to download. We also integrated currently prevalent technologies and toolkits to strengthen the capability of the database and provide a one-stop service for scientific users.

Nucleosome-specific, time-dependent changes in histone modifications during activation of the early growth response 1 (Egr1) gene

Histone post-translational modifications and nucleosome remodeling are coordinate events involved in eukaryotic transcriptional regulation. There are relatively few data on the time course with which these events occur in individual nucleosomes. As a contribution to fill this gap, we first describe the nature and time course of structural changes in the nucleosomes -2, -1, and +1 of the murine Egr1 gene upon induction. To initiate the transient activation of the gene, we used the stimulation of MLP29 cells with phorbol esters and the in vivo activation after partial hepatectomy. In both models, nucleosomes -1 and +1 are partially evicted, whereas nucleosomes +1 and -2 slide downstream during transcription. The sliding of the latter nucleosome allows the EGR1 protein to bind its site, resulting in the repression of the gene. To decide whether EGR1 is involved in the sliding of nucleosome -2, Egr1 was knocked down. In the absence of detectable EGR1, the nucleosome still slides and remains downstream longer than in control cells, suggesting that the product of the gene may be rather involved in the returning of the nucleosome to the basal position. Moreover, the presence of eight epigenetic histone marks has been determined at a mononucleosomal level in that chromatin region. H3S10phK14ac, H3K4me3, H3K9me3, and H3K27me3 are characteristic of nucleosome +1, and H3K9ac and H4K16ac are mainly found in nucleosome -1, and H3K27ac predominates in nucleosomes -2 and -1. The temporal changes in these marks suggest distinct functions for some of them, although changes in H3K4me3 may result from histone turnover.

Elucidating the molecular bases of epigenetic inheritance in non-model invertebrates: the case of the root-knot nematode Meloidogyne incognita

Root-knot nematodes of the genus Meloidogyne are biotrophic plant parasites that exhibit different life cycles and reproduction modes, ranging from classical amphimixis to obligatory mitotic parthenogenesis (apomixis), depending on the species. Meloidogyne incognita, an apomictic species, exhibits a worldwide distribution and a wide host range affecting more than 3000 plant species. Furthermore, evidences suggest that apomixis does not prevent M. incognita from adapting to its environment in contrast to what is expected from mitotic parthenogenesis that should theoretically produce clonal progenies. This raises questions about mechanisms of genome plasticity leading to genetic variation and adaptive evolution in apomictic animals. We reasoned that epigenetic mechanisms might in part be responsible for the generation of phenotypic variants that provide potential for rapid adaptation. We established therefore a pipeline to investigate the principal carriers of epigenetic information, DNA methylation and post-translational histone modifications. Even if M. incognita possesses the epigenetic machinery i.e., chromatin modifying enzymes, 5-methyl-cytosine and 5-hydroxy-methyl-cytosine content is absent or very weak. In contrast, we demonstrated that the canonical histone modifications are present and chromatin shows typical nucleosome structure. This work is the first characterization of carriers of epigenetic information in M. incognita and constitutes a preamble to further investigate if M. incognita development and its adaptation to plant hosts are under epigenetic control. Our pipeline should allow performing similar types of studies in any non-model organism.

Epigenetic regulation of rice flowering and reproduction

Current understanding of the epigenetic regulator roles in plant growth and development has largely derived from studies in the dicotyledonous model plant Arabidopsis thaliana. Rice (Oryza sativa) is one of the most important food crops in the world and has more recently becoming a monocotyledonous model plant in functional genomics research. During the past few years, an increasing number of studies have reported the impact of DNA methylation, non-coding RNAs and histone modifications on transcription regulation, flowering time control, and reproduction in rice. Here, we review these studies to provide an updated complete view about chromatin modifiers characterized in rice and in particular on their roles in epigenetic regulation of flowering time, reproduction, and seed development.