High-resolution mapping of epigenetic modifications of the rice genome uncovers interplay between DNA methylation, histone methylation, and gene expression

H3K4me2 functions as a repressive epigenetic mark in plants

BACKGROUND

In animals, H3K4me2 and H3K4me3 are enriched at the transcription start site (TSS) and function as epigenetic marks that regulate gene transcription, but their functions in plants have not been fully characterized.

RESULTS

We used chromatin immunoprecipitation sequencing to analyze the rice genome-wide changes to H3K4me1/H3K4me2/H3K4me3 following the loss of an H3K4-specific methyltransferase, SDG701. The knockdown of SDG701 resulted in a global decrease in H3K4me2/H3K4me3 levels throughout the rice genome. An RNA-sequencing analysis revealed that many genes related to diverse developmental processes were misregulated in the SDG701 knockdown mutant. In rice, H3K4me3 and H3K36me3 are positively correlated with gene transcription; however, surprisingly, the H3K4me2 level was negatively associated with gene transcription levels. Furthermore, the H3K4me3 level at the TSS region decreased significantly in the genes that exhibited down-regulated expression in the SDG701 knockdown mutant. In contrast, the genes with up-regulated expression in the mutant were associated with a considerable decrease in H3K4me2 levels over the gene body region.

CONCLUSION

A comparison of the genome-wide distributions of H3K4me2 in eukaryotes indicated that the H3K4me2 level is not correlated with the gene transcription level in yeast, but is positively and negatively correlated with gene expression in animals and plants, respectively. Our results uncovered H3K4me2 as a novel repressive mark in plants.

OsZIP1 functions as a metal efflux transporter limiting excess zinc, copper and cadmium accumulation in rice

BACKGROUND

Metal homeostasis is critical for plant growth, development and adaptation to environmental stresses and largely governed by a variety of metal transporters. The plant ZIP (Zn-regulated transporter, Iron-regulated transporter-like Protein) family proteins belong to the integral membrane transporters responsible for uptake and allocation of essential and non-essential metals. However, whether the ZIP family members mediate metal efflux and its regulatory mechanism remains unknown.

RESULTS

In this report, we provided evidence that OsZIP1 is a metal-detoxified transporter through preventing excess Zn, Cu and Cd accumulation in rice. OsZIP1 is abundantly expressed in roots throughout the life span and sufficiently induced by excess Zn, Cu and Cd but not by Mn and Fe at transcriptional and translational levels. Expression of OsZIP-GFP fusion in rice protoplasts and tobacco leaves shows that OsZIP1 resides in the endoplasmic reticulum (ER) and plasma membrane (PM). The yeast (Saccharomyces cerevisiae) complementation test shows that expression of OsZIP1 reduced Zn accumulation. Transgenic rice overexpressing OsZIP1 grew better under excess metal stress but accumulated less of the metals in plants. In contrast, both oszip1 mutant and RNA interference (RNAi) lines accumulated more metal in roots and contributed to metal sensitive phenotypes. These results suggest OsZIP1 is able to function as a metal exporter in rice when Zn, Cu and Cd are excess in environment. We further identified the DNA methylation of histone H3K9me2 of OsZIP1 and found that OsZIP1 locus, whose transcribed regions imbed a 242 bp sequence, is demethylated, suggesting that epigenetic modification is likely associated with OsZIP1 function under Cd stress.

CONCLUSION

OsZIP1 is a transporter that is required for detoxification of excess Zn, Cu and Cd in rice.

Mutations in the Rice OsCHR4 Gene, Encoding a CHD3 Family Chromatin Remodeler, Induce Narrow and Rolled Leaves with Increased Cuticular Wax

Leaf blade width, curvature, and cuticular wax are important agronomic traits of rice. Here, we report the rice Oschr4-5 mutant characterized by pleiotropic phenotypes, including narrow and rolled leaves, enhanced cuticular wax deposition and reduced plant height and tiller number. The reduced leaf width is caused by a reduced number of longitudinal veins and increased auxin content. The cuticular wax content was significantly higher in the Oschr4-5 mutant, resulting in reduced water loss rate and enhanced drought tolerance. Molecular characterization reveals that a single-base deletion results in a frame-shift mutation from the second chromodomain of OsCHR4, a CHD3 (chromodomain helicase DNA-binding) family chromatin remodeler, in the Oschr4-5 mutant. Expressions of seven wax biosynthesis genes (GL1-4, WSL4, OsCER7, LACS2, LACS7, ROC4 and BDG) and four auxin biosynthesis genes (YUC2, YUC3, YUC5 and YUC6) was up-regulated in the Oschr4-5 mutant. Chromatin immunoprecipitation assays revealed that the transcriptionally active histone modification H3K4me3 was increased, whereas the repressive H3K27me3 was reduced in the upregulated genes in the Oschr4-5 mutant. Therefore, OsCHR4 regulates leaf morphogenesis and cuticle wax formation by epigenetic modulation of auxin and wax biosynthetic genes expression.

Identification of Histone H3 (HH3) Genes in Gossypium hirsutum Revealed Diverse Expression During Ovule Development and Stress Responses

Histone acts as the core for nucleosomes and is a key protein component of chromatin. Among different histone variants, histone H3 (HH3) variants have been reported to play vital roles in plant development. However, biological information and evolutionary relationships of HH3 genes in cotton remain to be elucidated. The current study identified 34 HH3 genes in Gossypium hirsutum. Phylogenetic analysis classified HH3 genes of 19 plant species into eight distinct clades. Sequence logos analysis among Arabidopsis, rice, and G. hirsutum amino acid residues showed higher conservation in amino acids. Using collinearity analysis, we identified 81 orthologous/paralogous gene pairs among the four genomes (A, D, At, and Dt) of cotton. Further, orthologous/paralogous and the Ka/Ks ratio demonstrated that cotton HH3 genes experienced strong purifying selection pressure with restricted functional divergence resulting from segmental and whole genome duplication. Expression pattern analysis indicated that GhHH3 genes were preferentially expressed in cotton ovule tissues. Additionally, GhHH3 gene expression can be regulated by abiotic stresses (cold, heat, sodium chloride (NaCl), and polyethylene glycol (PEG)) and phytohormonal (brassinolide (BL), gibberellic acid (GA), indole-3-acetic acid (IAA), salicylic acid (SA), and methyl jasmonate (MeJA)) treatments, suggesting that GhHH3 genes might play roles in abiotic and hormone stress resistance. Taken together, this work provides important information to decipher complete molecular and physiological functions of HH3 genes in cotton.

Histone Modification Marks Strongly Regulate CDH1 Promoter in Prostospheres as A Model of Prostate Cancer Stem Like Cells

Objective

Cadherin-1 (CDH1) plays an important role in the metastasis, while expression of this protein is under control of epigenetic changes on its gene promoter. Therefore we evaluated both DNA methylation (DNAmet) and histone modification marks of CDH1 in prostate cancer stem like cells (PCSLCs).

Materials and Methods

In this experimental study, we isolated PCSLCs using cell surface marker and prostaspheroid formation, respectively. The cells isolated from both methods were characterized and then the levels of H3K4me2, H3K27me3, H3K9me2/3 and H3K9ac as well as DNAmet were assessed in CDH1 promoter of the isolated cells.

Results

The CD44⁺ CD49^hi cells were not validated as PCSLCs. However, prostaspheres overexpressed stemness related genes and had higher ability of invasion potential, associated with reduction in CDH1 expression. Epigenetic status analysis showed that CDH1 promoter was hypo-methylated. Histone modifications of H3K9ac and H3K4me3 were significantly reduced, in parallel with an increased level of H3K27me3.

Conclusion

Our results suggest that slight decrease of DNAmet of the CpG island in CDH1 promoter does not significantly contribute to the change of CDH1 expression. Therefore, histone modifications are responsible in repressing CDH1 in PCSLCs.

Evolutionary Epigenomics of Retrotransposon-Mediated Methylation Spreading in Rice

Plant genomes contain numerous transposable elements (TEs), and many hypotheses on the evolutionary drivers that restrict TE activity have been postulated. Few models, however, have focused on the evolutionary epigenomic interaction between the plant host and its TE. The host genome recruits epigenetic factors, such as methylation, to silence TEs but methylation can spread beyond the TE sequence and influence the expression of nearby host genes. In this study, we investigated this epigenetic trade-off between TE and proximal host gene silencing by studying the epigenomic regulation of repressing long terminal repeat (LTR) retrotransposons (RTs) in Oryza sativa. Results showed significant evidence of methylation spreading originating from the LTR-RT sequences, and the extent of spreading was dependent on five factors: 1) LTR-RT family, 2) time since the LTR-RT insertion, 3) recombination rate of the LTR-RT region, 4) level of LTR-RT sequence methylation, and 5) chromosomal location. Methylation spreading had negative effects by reducing host gene expression, but only on host genes with LTR-RT inserted in its introns. Our results also suggested high levels of LTR-RT methylation might have a role in suppressing TE-mediated deleterious ectopic recombination. In the end, despite the methylation spreading, no strong epigenetic trade-off was detected and majority of LTR-RT may have only minor epigenetic effects on nearby host genes.

Extensive allele-level remodeling of histone methylation modification in reciprocal F1 hybrids of rice subspecies

Epigenetic mechanisms play a major role in heterosis, partly as a result of the remodeling of epigenetic modifications in F₁ hybrids. Based on chromatin immunoprecipitation-sequencing (ChIP-Seq) analyses, we show that at the allele level extensive histone methylation remodeling occurred for a subset of genomic loci in reciprocal F₁ hybrids of Oryza sativa (rice) cultivars Nipponbare and 93-11, representing the two subspecies japonica and indica. Globally, the allele modification-altered loci in leaf or root of the reciprocal F₁ hybrids involved ˜12-43% or more of the genomic regions carrying either of two typical histone methylation markers, H3K4me3 (>21 000 genomic regions) and H3K27me3 (>11 000 genomic regions). Nevertheless, at the total modification level, the majority (from ˜43 to >90%) of the modification-altered alleles lay within the range of parental additivity in the hybrids because of concerted alteration in opposite directions, consistent with an overall attenuation of allelic differences in the modifications. Importantly, of the genomic regions that did show non-additivity in total modification level by either marker in the two tissues of hybrids, >80% manifested transgressivity, which involved genes enriched in specific functional categories. Extensive allele-level alteration of H3K4me3 alone was positively correlated with genome-wide changes in allele-level gene expression, whereas at the total level, both H3K4me3 and H3K27me3 remodeling, although affecting just a small number of genes, contributes to the overall non-additive gene expression to variable extents, depending on tissue/marker combinations. Our results emphasize the importance of allele-level analysis in hybrids to assess the remodeling of epigenetic modifications and their relation to changes in gene expression.

Tolerance of Iron-Deficient and -Toxic Soil Conditions in Rice

Iron (Fe) deficiency and toxicity are the most widely prevalent soil-related micronutrient disorders in rice (Oryza sativa L.). Progress in rice cultivars with improved tolerance has been hampered by a poor understanding of Fe availability in the soil, the transportation mechanism, and associated genetic factors for the tolerance of Fe toxicity soil (FTS) or Fe deficiency soil (FDS) conditions. In the past, through conventional breeding approaches, rice varieties were developed especially suitable for low- and high-pH soils, which indirectly helped the varieties to tolerate FTS and FDS conditions. Rice-Fe interactions in the external environment of soil, internal homeostasis, and transportation have been studied extensively in the past few decades. However, the molecular and physiological mechanisms of Fe uptake and transport need to be characterized in response to the tolerance of morpho-physiological traits under Fe-toxic and -deficient soil conditions, and these traits need to be well integrated into breeding programs. A deeper understanding of the several factors that influence Fe absorption, uptake, and transport from soil to root and above-ground organs under FDS and FTS is needed to develop tolerant rice cultivars with improved grain yield. Therefore, the objective of this review paper is to congregate the different phenotypic screening methodologies for prospecting tolerant rice varieties and their responsible genetic traits, and Fe homeostasis related to all the known quantitative trait loci (QTLs), genes, and transporters, which could offer enormous information to rice breeders and biotechnologists to develop rice cultivars tolerant of Fe toxicity or deficiency. The mechanism of Fe regulation and transport from soil to grain needs to be understood in a systematic manner along with the cascade of metabolomics steps that are involved in the development of rice varieties tolerant of FTS and FDS. Therefore, the integration of breeding with advanced genome sequencing and omics technologies allows for the fine-tuning of tolerant genotypes on the basis of molecular genetics, and the further identification of novel genes and transporters that are related to Fe regulation from FTS and FDS conditions is incredibly important to achieve further success in this aspect.

An integrated package for bisulfite DNA methylation data analysis with Indel-sensitive mapping

Background

DNA methylation plays crucial roles in most eukaryotic organisms. Bisulfite sequencing (BS-Seq) is a sequencing approach that provides quantitative cytosine methylation levels in genome-wide scope and single-base resolution. However, genomic variations such as insertions and deletions (indels) affect methylation calling, and the alignment of reads near/across indels becomes inaccurate in the presence of polymorphisms. Hence, the simultaneous detection of DNA methylation and indels is important for exploring the mechanisms of functional regulation in organisms.

Results

These problems motivated us to develop the algorithm BatMeth2, which can align BS reads with high accuracy while allowing for variable-length indels with respect to the reference genome. The results from simulated and real bisulfite DNA methylation data demonstrated that our proposed method increases alignment accuracy. Additionally, BatMeth2 can calculate the methylation levels of individual loci, genomic regions or functional regions such as genes/transposable elements. Additional programs were also developed to provide methylation data annotation, visualization, and differentially methylated cytosine/region (DMC/DMR) detection. The whole package provides new tools and will benefit bisulfite data analysis.

Conclusion

BatMeth2 improves DNA methylation calling, particularly for regions close to indels. It is an autorun package and easy to use. In addition, a DNA methylation visualization program and a differential analysis program are provided in BatMeth2. We believe that BatMeth2 will facilitate the study of the mechanisms of DNA methylation in development and disease. BatMeth2 is an open source software program and is available on GitHub (https://github.com/GuoliangLi-HZAU/BatMeth2/).

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2593-4) contains supplementary material, which is available to authorized users.

Dynamic Changes in Genome-Wide Histone3 Lysine27 Trimethylation and Gene Expression of Soybean Roots in Response to Salt Stress

Soybean is an important economic crop for human diet, animal feeds and biodiesel due to high protein and oil content. Its productivity is significantly hampered by salt stress, which impairs plant growth and development by affecting gene expression, in part, through epigenetic modification of chromatin status. However, little is known about epigenetic regulation of stress response in soybean roots. Here, we used RNA-seq and ChIP-seq technologies to study the dynamics of genome-wide transcription and histone methylation patterns in soybean roots under salt stress. Eight thousand seven hundred ninety eight soybean genes changed their expression under salt stress treatment. Whole-genome ChIP-seq study of an epigenetic repressive mark, histone H3 lysine 27 trimethylation (H3K27me3), revealed the changes in H3K27me3 deposition during the response to salt stress. Unexpectedly, we found that most of the inactivation of genes under salt stress is strongly correlated with the de novo establishment of H3K27me3 in various parts of the promoter or coding regions where there is no H3K27me3 in control plants. In addition, the soybean histone modifiers were identified which may contribute to de novo histone methylation and gene silencing under salt stress. Thus, dynamic chromatin regulation, switch between active and inactive modes, occur at target loci in order to respond to salt stress in soybean. Our analysis demonstrates histone methylation modifications are correlated with the activation or inactivation of salt-inducible genes in soybean roots.

Effects of ploidy variation on promoter DNA methylation and gene expression in rice (Oryza sativa L.)

BACKGROUND

Polyploidy, or whole-genome duplication (WGD) promotes genetic diversification in plants. However, whether WGD is accompanied by epigenetic regulation especially DNA methylation remains yet elusive. Methylation of different region in genomic DNA play discrete role in gene regulation and developmental processes in plants.

RESULTS

In our study, we used an apomictic rice line (SARII-628) that produces twin seedlings of different ploidy for methylated DNA immunoprecipitation sequencing (MeDIP-seq). We compared the level of methylation and mRNA expression in three different (CG, CHG, and CHH) sequence contexts of promoter region among haploid (1X), diploid (2X), and triploid (3X) seedling. We used MeDIP-Seq analysis of 14 genes to investigate whole genome DNA methylation and found that relative level of DNA methylation across different ploidy was in following order e.g. diploid > triploid > haploid. GO functional classification of differentially methylated genes into 9 comparisons group of promoter, intergenic and intragenic region discovered, these genes were mostly enriched for cellular component, molecular function, and biological process. By the comparison of methylome data, digital gene expression (DGE), mRNA expression profile, and Q-PCR findings LOC_ Os07g31450 and LOC_ Os01g59320 were analyzed for BS-Seq (Bisulphite sequencing).

CONCLUSIONS

We found that (1) The level of the promoter DNA methylation is negatively correlated with gene expression within each ploidy level. (2) Among all ploidy levels, CG sequence context had highest methylation frequency, and demonstrated that the high CG methylation did reduce gene expression change suggesting that DNA methylation exert repressive function and ensure genome stability during WGD. (3) Alteration in ploidy (from diploid to haploid, or diploid to triploid) reveals supreme changes in methylation frequency of CHH sequence context. Our finding will contribute an understanding towards lower stability of CHH sequence context and educate the effect of promoter region methylation during change in ploidy state in rice.

Retrospective and perspective of plant epigenetics in China

Epigenetics refers to the study of heritable changes in gene function that do not involve changes in the DNA sequence. Such effects on cellular and physiological phenotypic traits may result from external or environmental factors or be part of normal developmental program. In eukaryotes, DNA wraps on a histone octamer (two copies of H2A, H2B, H3 and H4) to form nucleosome, the fundamental unit of chromatin. The structure of chromatin is subjected to a dynamic regulation through multiple epigenetic mechanisms, including DNA methylation, histone posttranslational modifications (PTMs), chromatin remodeling and noncoding RNAs. As conserved regulatory mechanisms in gene expression, epigenetic mechanisms participate in almost all the important biological processes ranging from basal development to environmental response. Importantly, all of the major epigenetic mechanisms in mammalians also occur in plants. Plant studies have provided numerous important contributions to the epigenetic research. For example, gene imprinting, a mechanism of parental allele-specific gene expression, was firstly observed in maize; evidence of paramutation, an epigenetic phenomenon that one allele acts in a single locus to induce a heritable change in the other allele, was firstly reported in maize and tomato. Moreover, some unique epigenetic mechanisms have been evolved in plants. For example, the 24-nt siRNA-involved RNA-directed DNA methylation (RdDM) pathway is plant-specific because of the involvements of two plant-specific DNA-dependent RNA polymerases, Pol IV and Pol V. A thorough study of epigenetic mechanisms is of great significance to improve crop agronomic traits and environmental adaptability. In this review, we make a brief summary of important progress achieved in plant epigenetics field in China over the past several decades and give a brief outlook on future research prospects. We focus our review on DNA methylation and histone PTMs, the two most important aspects of epigenetic mechanisms.

Parental experience modifies the Mimulus methylome

BACKGROUND

Transgenerational plasticity occurs when the environmental experience of an organism modifies the growth and development of its progeny. Leaf damage in Mimulus guttatus exhibits transgenerational plasticity mediated through differential expression of hundreds of genes. The epigenetic mechanisms that facilitate this response have yet to be described.

RESULTS

We performed whole genome bisulfite sequencing in the progeny of genetically identical damaged and control plants and developed a pipeline to compare differences in the mean and variance of methylation between treatment groups. We find that parental damage increases the variability of CG and CHG methylation among progeny, but does not alter the overall mean methylation. Instead it has positive effects in some regions and negative in others. We find 3,396 CHH, 203 CG, and 54 CHG Differentially Methylated Regions (DMRs) ranging from tens to thousands of base pairs scattered across the genome. CHG and CHH DMRs tended to overlap with transposable elements. CG DMRs tended to overlap with gene coding regions, many of which were previously found to be differentially expressed.

CONCLUSIONS

Genome-wide increases in methylome variation suggest that parental conditions can increase epigenetic diversity in response to stress. Additionally, the potential association between CG DMRs and differentially expressed genes supports the hypothesis that differential methylation is a mechanistic component of transgenerational plasticity in M. guttatus.

Epigenetic control of defense genes following MeJA-induced priming in rice (O. sativa)

In plant defense, priming is a physiological process by which a plant prepares to more quickly or aggressively respond to a future biotic or abiotic stress. This unique physiological state can be induced by beneficial microbes or by treatment with natural or synthetic compounds. Priming has been described copiously, mainly in dicots, but the study on the molecular mechanisms regulating the phenomenon is still at its infancy. A fascinating possibility is that epigenetic changes may be the main actors on priming establishment. A more accessible chromatin structure, due to DNA methylation and/or histone modifications, could facilitate a quicker and more potent gene response to a subsequent attack. Here, we investigated the impact of methyl jasmonate (MeJA)-induced priming on the response efficacy to mechanical wounding of a monocot (Oryza sativa). In particular, we showed that MeJA primes plants for increased expression of defense-related genes, such as OsBBPI and OsPOX, upon wounding. Finally, we provided evidence that MeJA modulates histone modifications in the promoter region of OsBBPI, as well as changes at genome-wide DNA methylation level. Together these studies corroborate the importance of priming in strengthening plant defense and support the growing evidence that epigenetic regulation plays a pivotal role in priming onset.

Substantial Epigenetic Variation Causing Flower Color Chimerism in the Ornamental Tree Prunus mume Revealed by Single Base Resolution Methylome Detection and Transcriptome Sequencing

Epigenetic changes caused by methylcytosine modification participate in gene regulation and transposable element (TE) repression, resulting in phenotypic variation. Although the effects of DNA methylation and TE repression on flower, fruit, seed coat, and leaf pigmentation have been investigated, little is known about the relationship between methylation and flower color chimerism. In this study, we used a comparative methylomic⁻transcriptomic approach to explore the molecular mechanism responsible for chimeric flowers in Prunus mume "Danban Tiaozhi". High-performance liquid chromatography-electrospray ionization mass spectrometry revealed that the variation in white (WT) and red (RT) petal tissues in this species is directly due to the accumulation of anthocyanins, i.e., cyanidin 3,5-O-diglucoside, cyanidin 3-O-glucoside, and peonidin 3-O-glucoside. We next mapped the first-ever generated methylomes of P. mume, and found that 11.29⁻14.83% of the genomic cytosine sites were methylated. We also determined that gene expression was negatively correlated with methylcytosine level in general, and uncovered significant epigenetic variation between WT and RT. Furthermore, we detected differentially methylated regions (DMRs) and DMR-related genes between WT and RT, and concluded that many of these genes, including differentially expressed genes (DEGs) and transcription factor genes, are critical participants in the anthocyanin regulatory pathway. Importantly, some of the associated DEGs harbored TE insertions that were also modified by methylcytosine. The above evidence suggest that flower color chimerism in P. mume is induced by the DNA methylation of critical genes and TEs.

Identification and analysis of adenine N6-methylation sites in the rice genome

DNA N⁶-methyladenine (6mA) is a non-canonical DNA modification that is present at low levels in different eukaryotes^1-8, but its prevalence and genomic function in higher plants are unclear. Using mass spectrometry, immunoprecipitation and validation with analysis of single-molecule real-time sequencing, we observed that about 0.2% of all adenines are 6mA methylated in the rice genome. 6mA occurs most frequently at GAGG motifs and is mapped to about 20% of genes and 14% of transposable elements. In promoters, 6mA marks silent genes, but in bodies correlates with gene activity. 6mA overlaps with 5-methylcytosine (5mC) at CG sites in gene bodies and is complementary to 5mC at CHH sites in transposable elements. We show that OsALKBH1 may be potentially involved in 6mA demethylation in rice. The results suggest that 6mA is complementary to 5mC as an epigenomic mark in rice and reinforce a distinct role for 6mA as a gene expression-associated epigenomic mark in eukaryotes.

Seasonal Growth of Zygophyllum dumosum Boiss.: Summer Dormancy Is Associated with Loss of the Permissive Epigenetic Marker Dimethyl H3K4 and Extensive Reduction in Proteins Involved in Basic Cell Functions

Plants thriving in desert environments are suitable for studying mechanisms for plant survival under extreme seasonal climate variation. We studied epigenetic mechanisms underlying seasonal growth cycles in the desert plant Zygophyllum dumosum Boiss., which was previously shown to be deficient in repressive markers of di-methyl and tri-methyl H3K9 and their association with factors regulating basic cell functions. We showed a contingent association between rainfall and seasonal growth and the epigenetic marker of dimethyl H3K4, which disappears upon entry into the dry season and the acquisition of a dormant state. DNA methylation is not affected by a lack of H3K9 di-methyl and tri-methyl. Changes in methylation can occur between the wet and dry season. Proteome analysis of acid soluble fractions revealed an extensive reduction in ribosomal proteins and in proteins involved in chloroplasts and mitochondrial activities during the dry seasons concomitantly with up-regulation of molecular chaperone HSPs. Our results highlight mechanisms underlying Z. dumosum adaptation to seasonal climate variation. Particularly, summer dormancy is associated with a loss of the permissive epigenetic marker dimethyl H3K4, which might facilitate genome compaction concomitantly with a significant reduction in proteins involved in basic cell functions. HSP chaperones might safeguard the integrity of cell components.

Beyond the genetic code in leaf senescence

Leaf senescence is not only genetically programmed but also induced by exogenous stress to ensure completion of the plant life cycle, successful reproduction and environmental adaptability. Genetic reprogramming is a major aspect of leaf senescence, and the senescence signaling that follows is controlled by a complex regulatory network. Recent studies suggest that the activity of transcription factors together with epigenetic mechanisms ensures the robustness of this network, with the latter including chromatin remodeling, DNA modification, and RNA-mediated control of transcription factors and other senescence-associated genes. In this review, we provide an overview of the relevant epigenetic mechanisms and summarize recent findings of epigenetic regulators of plant leaf senescence involved in DNA methylation and histone modification along with the functions of small RNAs in this process.

Analysis of Transposable Elements in Coccidioides Species

Coccidioides immitis and C. posadasii are primary pathogenic fungi that cause disease in immunologically-normal animals and people. The organism is found exclusively in arid regions of the Southwestern United States, Mexico, and South America, but not in other parts of the world. This study is a detailed analysis of the transposable elements (TE) in Coccidioides spp. As is common in most fungi, Class I and Class II transposons were identified and the LTR Gypsy superfamily is the most common. The minority of Coccidioides Gypsy transposons contained regions highly homologous to polyprotein domains. Phylogenetic analysis of the integrase and reverse transcriptase sequences revealed that many, but not all, of the Gypsy reverse transcriptase and integrase domains clustered by species suggesting extensive transposition after speciation of the two Coccidiodies spp. The TEs were clustered and the distribution is enriched for the ends on contigs. Analysis of gene expression data from C. immitis found that protein-coding genes within 1 kB of hAT or Gypsy TEs were poorly expressed. The expression of C. posadasii genes within 1 kB of Gypsy TEs was also significantly lower compared to all genes but the difference in expression was smaller than C. immitis. C. posadasii orthologs of C. immitis Gyspsy-associated genes were also likely to be TE-associated. In both C. immitis and C. posadasii the TEs were preferentially associated with genes annotated with protein kinase gene ontology terms. These observations suggest that TE may play a role in influencing gene expression in Coccidioides spp. Our hope is that these bioinformatic studies of the potential TE influence on expression and evolution of Coccidioides will prompt the development of testable hypotheses to better understand the role of TEs in the biology and gene regulation of Coccidioides spp.

The gymnastics of epigenomics in rice

Epigenomics is represented by the high-throughput investigations of genome-wide epigenetic alterations, which ultimately dictate genomic, transcriptomic, proteomic and metabolomic dynamism. Rice has been accepted as the global staple crop. As a result, this model crop deserves significant importance in the rapidly emerging field of plant epigenomics. A large number of recently available data reveal the immense flexibility and potential of variable epigenomic landscapes. Such epigenomic impacts and variability are determined by a number of epigenetic regulators and several crucial inheritable epialleles, respectively. This article highlights the correlation of the epigenomic landscape with growth, flowering, reproduction, non-coding RNA-mediated post-transcriptional regulation, transposon mobility and even heterosis in rice. We have also discussed the drastic epigenetic alterations which are reported in rice plants grown from seeds exposed to the extraterrestrial environment. Such abiotic conditions impose stress on the plants leading to epigenomic modifications in a genotype-specific manner. Some significant bioinformatic databases and in silico approaches have also been explained in this article. These softwares provide important interfaces for comparative epigenomics. The discussion concludes with a unified goal of developing epigenome editing to promote biological hacking of the rice epigenome. Such a cutting-edge technology if properly standardized, can integrate genomics and epigenomics together with the generation of high-yielding trait in several cultivars of rice.

The Rice Rolled Fine Striped (RFS) CHD3/Mi-2 Chromatin Remodeling Factor Epigenetically Regulates Genes Involved in Oxidative Stress Responses During Leaf Development

In rice (Oryza sativa), moderate leaf rolling increases photosynthetic competence and raises grain yield; therefore, this important agronomic trait has attracted much attention from plant biologists and breeders. However, the relevant molecular mechanism remains unclear. Here, we isolated and characterized Rolled Fine Striped (RFS), a key gene affecting rice leaf rolling, chloroplast development, and reactive oxygen species (ROS) scavenging. The rfs-1 gamma-ray allele and the rfs-2 T-DNA insertion allele of RFS failed to complement each other and their mutants had similar phenotypes, producing extremely incurved leaves due to defective development of vascular cells on the adaxial side. Map-based cloning showed that the rfs-1 mutant harbors a 9-bp deletion in a gene encoding a predicted CHD3/Mi-2 chromatin remodeling factor belonging to the SNF2-ATP-dependent chromatin remodeling family. RFS was expressed in various tissues and accumulated mainly in the vascular cells throughout leaf development. Furthermore, RFS deficiency resulted in a cell death phenotype that was caused by ROS accumulation in developing leaves. We found that expression of five ROS-scavenging genes [encoding catalase C, ascorbate peroxidase 8, a putative copper/zinc superoxide dismutase (SOD), a putative SOD, and peroxiredoxin IIE2] decreased in rfs-2 mutants. Western-blot and chromatin immunoprecipitation (ChIP) assays demonstrated that rfs-2 mutants have reduced H3K4me3 levels in ROS-related genes. Loss-of-function in RFS also led to multiple developmental defects, affecting pollen development, grain filling, and root development. Our results suggest that RFS is required for many aspects of plant development and its function is closely associated with epigenetic regulation of genes that modulate ROS homeostasis.

Bioinformatics resources for deciphering the biogenesis and action pathways of plant small RNAs

The next-generation sequencing (NGS) technology has revolutionized our previous understanding of the plant genomes, relying on its innate advantages, such as high throughput and deep sequencing depth. In addition to the protein-coding gene loci, massive transcription signals have been detected within intergenic or intragenic regions. Most of these signals belong to non-coding ones, considering their weak protein-coding potential. Generally, these transcripts could be divided into long non-coding RNAs and small non-coding RNAs (sRNAs) based on their sequence length. In addition to the well-known microRNAs (miRNAs), many plant endogenous sRNAs were collectively referred to as small interfering RNAs. However, an increasing number of unclassified sRNA species are being discovered by NGS. The high heterogeneity of these novel sRNAs greatly hampered the mechanistic studies, especially on the clear description of their biogenesis and action pathways. Fortunately, public databases, bioinformatics softwares and NGS datasets are increasingly available for plant sRNA research. Here, by summarizing these valuable resources, we proposed a general workflow to decipher the RDR (RNA-dependent RNA polymerase)- and DCL (Dicer-like)-dependent biogenesis pathways, and the Argonaute-mediated action modes (such as target cleavages and chromatin modifications) for specific sRNA species in plants. Taken together, we hope that by summarizing a list of the public resources, this work will facilitate the plant biologists to perform classification and functional characterization of the interesting sRNA species.

Dynamic changes in methylome and transcriptome patterns in response to methyltransferase inhibitor 5-azacytidine treatment in citrus

DNA methylation is known to play an important role in various developmental processes in plants. However, there is a general lack of understanding about the possible functions of DNA methylation in fruit trees. Using callus as a model, methylome, transcriptome and metabolite changes were assessed after treatment with the DNA methyltransferase inhibitor 5-azacytidine (5azaC). Genome-wide methylome analysis revealed the demethylation of a diverse of genes, including many genes encoding transcription factors (TFs), genes involved in biological processes, and the up-regulation of a wide range of transposable elements (TEs). Combined with the RNA-seq data, we observed no obvious genome-wide correlation between the changes in methylation status and expression levels. Furthermore, 5azaC treatment induced carotenoid degradation along with strong activation of carotenoid cleavage dioxygenases 1 (CpCCD1). Functional complementation analysis in bacterial system showed that CpCCD1 exhibited strong catalytic activities toward zeaxanthin, β-carotene and lycopene. In summary, 5azaC treatments induced carotenoid degradation by CpCCD1 activation and led to a genome-wide demethylation effect.

Small RNA Activity in Archeological Barley Shows Novel Germination Inhibition in Response to Environment

The recovery of ancient RNA from archeological material could enable the direct study of microevolutionary processes. Small RNAs are a rich source of information because their small size is compatible with biomolecular preservation, and their roles in gene regulation make them likely foci of evolutionary change. We present here the small RNA fraction from a sample of archeological barley generated using high-throughput sequencing that has previously been associated with localized adaptation to drought. Its microRNA profile is broadly similar to 19 globally distributed modern barley samples with the exception of three microRNAs (miRNA159, miRNA319, and miR396), all of which are known to have variable expression under stress conditions. We also found retrotransposon activity to be significantly reduced in the archeological barley compared with the controls, where one would expect the opposite under stress conditions. We suggest that the archeological barley's conflicting stress signals could be the result of long-term adaptation to its local environment.

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.

SET DOMAIN GROUP701 encodes a H3K4-methytransferase and regulates multiple key processes of rice plant development

Chromatin-based epigenetic information plays an important role in developmental gene regulation, in response to environment, and in natural variation of gene expression levels. Histone H3 lysine 4 di/trimethylation (H3K4me2/3) is abundant in euchromatin and is generally associated with transcriptional activation. Strikingly, however, enzymes catalyzing H3K4me2/3 remain poorly characterized in crops so far. Here, we investigated the function of the rice SET DOMAIN GROUP 701 (SDG701) gene by molecular and biochemical characterization of the gene product, and by studying effects of its loss or gain of function on plant growth and development. We demonstrated that SDG701 encodes a methytransferase specifically catalyzing H3K4 methylation. Overexpression and knockdown experiments showed that SDG701 is crucial for proper sporophytic plant development as well as for gametophytic transmission that directly impacts rice grain production. In-depth analysis of plant flowering time revealed that SDG701 promotes rice flowering under either long-day or short-day photoperiods. Consistently, the SDG701 protein was found to bind chromatin to promote H3K4me3 and to enhance expression of the rice Hd3a and RFT1 florigens. Collectively, our results establish SDG701 as a major rice H3K4-specific methyltransferase and provide important insights into function of H3K4me3 deposition in transcription activation of florigens in promoting plant flowering.

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.

Epigenomic and functional analyses reveal roles of epialleles in the loss of photoperiod sensitivity during domestication of allotetraploid cottons

Background

Polyploidy is a pervasive evolutionary feature of all flowering plants and some animals, leading to genetic and epigenetic changes that affect gene expression and morphology. DNA methylation changes can produce meiotically stable epialleles, which are transmissible through selection and breeding. However, the relationship between DNA methylation and polyploid plant domestication remains elusive.

Results

We report comprehensive epigenomic and functional analyses, including ~12 million differentially methylated cytosines in domesticated allotetraploid cottons and their tetraploid and diploid relatives. Methylated genes evolve faster than unmethylated genes; DNA methylation changes between homoeologous loci are associated with homoeolog-expression bias in the allotetraploids. Significantly, methylation changes induced in the interspecific hybrids are largely maintained in the allotetraploids. Among 519 differentially methylated genes identified between wild and cultivated cottons, some contribute to domestication traits, including flowering time and seed dormancy. CONSTANS (CO) and CO-LIKE (COL) genes regulate photoperiodicity in Arabidopsis. COL2 is an epiallele in allotetraploid cottons. COL2A is hypermethylated and silenced, while COL2D is repressed in wild cottons but highly expressed due to methylation loss in all domesticated cottons tested. Inhibiting DNA methylation activates COL2 expression, and repressing COL2 in cultivated cotton delays flowering.

Conclusions

We uncover epigenomic signatures of domestication traits during cotton evolution. Demethylation of COL2 increases its expression, inducing photoperiodic flowering, which could have contributed to the suitability of cotton for cultivation worldwide. These resources should facilitate epigenetic engineering, breeding, and improvement of polyploid crops.

Electronic supplementary material

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

Maintaining the factory: the roles of the unfolded protein response in cellular homeostasis in plants

Much like a factory, the endoplasmic reticulum (ER) assembles simple cellular building blocks into complex molecular machines known as proteins. In order to protect the delicate protein folding process and ensure the proper cellular delivery of protein products under environmental stresses, eukaryotes have evolved a set of signaling mechanisms known as the unfolded protein response (UPR) to increase the folding capacity of the ER. This process is particularly important in plants, because their sessile nature commands adaptation for survival rather than escape from stress. As such, plants make special use of the UPR, and evidence indicates that the master regulators and downstream effectors of the UPR have distinct roles in mediating cellular processes that affect organism growth and development as well as stress responses. In this review we outline recent developments in this field that support a strong relevance of the UPR to many areas of plant life.

Chromatin modification contributes to the expression divergence of three TaGS2 homoeologs in hexaploid wheat

Plastic glutamine synthetase (GS2) is responsible for ammonium assimilation. The reason that TaGS2 homoeologs in hexaploid wheat experience different selection pressures in the breeding process remains unclear. TaGS2 were minimally expressed in roots but predominantly expressed in leaves, and TaGS2-B had higher expression than TaGS2-A and TaGS2-D. ChIP assays revealed that the activation of TaGS2-B expression in leaves was correlated with increased H3K4 trimethylation. The transcriptional silencing of TaGS2 in roots was correlated with greater cytosine methylation and less H3K4 trimethylation. Micrococcal nuclease and DNase I accessibility experiments indicated that the promoter region was more resistant to digestion in roots than leaves, which indicated that the closed nucleosome conformation of the promoter region was important to the transcription initiation for the spatial-temporal expression of TaGS2. In contrast, the transcribed regions possess different nuclease accessibilities of three TaGS2 homoeologs in the same tissue, suggesting that nucleosome conformation of the transcribed region was part of the fine adjustment of TaGS2 homoeologs. This study provides evidence that histone modification, DNA methylation and nuclease accessibility coordinated the control of the transcription of TaGS2 homoeologs. Our results provided important evidence that TaGS2-B experienced the strongest selection pressures during the breeding process.

Dynamic and Antagonistic Allele-Specific Epigenetic Modifications Controlling the Expression of Imprinted Genes in Maize Endosperm

Genomic imprinting is often associated with allele-specific epigenetic modifications. Although many reports suggested potential roles of DNA methylation and H3K27me3 in regulating genomic imprinting, the contributions of allele-specific active histone modifications to imprinting remain still unclear in plants. Here, we report the identification of 337 high-stringency allele-specific H3K4me3 and H3K36me3 peaks in maize endosperm. Paternally preferred H3K4me3 and H3K36me3 peaks mostly co-localized with paternally expressed genes (PEGs), while endosperm-specific maternally expressed genes (endo-MEGs) were associated with maternally preferred H3K4me3 and H3K36me3 peaks. A unique signature for PEGs was observed, where the active H3K4me4 and H3K36me3 as well as repressive H3K27me3 appeared together. At the gene body of con-PEGs (constitutively expressed PEG), H3K27me3 and H3K36me3 were specifically deposited on hypomethylated maternal alleles and hypermethylated paternal alleles, respectively. Around the transcription start sites of endo-MEGs, DNA methylation and H3K4me3 specifically marked paternal and maternal alleles, respectively. In addition, 35 maternally expressed non-coding RNAs exhibited the same allele-specific epigenetic features as endo-MEGs, indicating similar mechanisms for the regulation of imprinted genes and non-coding RNAs. Taken together, our results uncover the complex patterns of mutually exclusive epigenetic modifications deposited at different alleles of imprinted genes that are required for genomic imprinting in maize endosperm.

EpiHRMAssay, in tube and in silico combined approach for the scanning and epityping of heterogeneous DNA methylation

Reliable and cost-effective assays with adequate sensitivity are required to detect the DNA methylation profile in plants for scientific and industrial purposes. The proposed novel assay, named EpiHRMAssay, allows to quantify the overall methylation status at target loci and to enable high-throughput analyses. It combines in tube High Resolution Melting Analysis on bisulphite-treated templates with the in silico prediction of the melting profile of virtual epialleles using uMELT^SM software. The predicted melting temperatures (T_m-s) of a set of epialleles characterized by different numbers of methylated cytosines (#^mC) or different ^mC configurations were obtained and used to build calibration models, enabling the quantification of methylation in unknown samples using only the in tube observed melting temperature (T_m-o). EpiHRMAssay was validated by analysing the promoter region of CMT3, DDM1, and ROS1 genes involved in the regulation of methylation/demethylation processes and chromatin remodelling within a population of peach plants. Results demonstrate that EpiHRMAssay is a sensitive and reliable tool for locus-specific large-scale research and diagnostic contexts of the regulative regions of genes, in a broad range of organisms, including mammals. EpiHRMAssay also provides complementary information for the assessment of heterogeneous methylation and can address an array of biological questions on epigenetic regulation for diversity studies and for large-scale functional genomics.

The DNA methylation level against the background of the genome size and t-heterochromatin content in some species of the genus Secale L

Methylation of cytosine in DNA is one of the most important epigenetic modifications in eukaryotes and plays a crucial role in the regulation of gene activity and the maintenance of genomic integrity. DNA methylation and other epigenetic mechanisms affect the development, differentiation or the response of plants to biotic and abiotic stress. This study compared the level of methylation of cytosines on a global (ELISA) and genomic scale (MSAP) between the species of the genus Secale. We analyzed whether the interspecific variation of cytosine methylation was associated with the size of the genome (C-value) and the content of telomeric heterochromatin. MSAP analysis showed that S. sylvestre was the most distinct species among the studied rye taxa; however, the results clearly indicated that these differences were not statistically significant. The total methylation level of the studied loci was very similar in all taxa and ranged from 60% in S. strictum ssp. africanum to 66% in S. cereale ssp. segetale, which confirmed the lack of significant differences in the sequence methylation pattern between the pairs of rye taxa. The level of global cytosine methylation in the DNA was not significantly associated with the content of t-heterochromatin and did not overlap with the existing taxonomic rye relationships. The highest content of 5-methylcytosine was found in S. cereale ssp. segetale (83%), while very low in S. strictum ssp. strictum (53%), which was significantly different from the methylation state of all taxa, except for S. sylvestre. The other studied taxa of rye had a similar level of methylated cytosine ranging from 66.42% (S. vavilovii) to 74.41% in (S. cereale ssp. afghanicum). The results obtained in this study are evidence that the percentage of methylated cytosine cannot be inferred solely based on the genome size or t-heterochromatin. This is a significantly more complex issue.

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.

Analyses of Methylomes Derived from Meso-American Common Bean (Phaseolus vulgaris L.) Using MeDIP-Seq and Whole Genome Sodium Bisulfite-Sequencing

Common bean (Phaseolus vulgaris L.) is economically important for its high protein, fiber, and micronutrient contents, with a relatively small genome size of ∼587 Mb. Common bean is genetically diverse with two major gene pools, Meso-American and Andean. The phenotypic variability within common bean is partly attributed to the genetic diversity and epigenetic changes that are largely influenced by environmental factors. It is well established that an important epigenetic regulator of gene expression is DNA methylation. Here, we present results generated from two high-throughput sequencing technologies, methylated DNA immunoprecipitation-sequencing (MeDIP-seq) and whole genome bisulfite-sequencing (BS-Seq). Our analyses revealed that this Meso-American common bean displays similar methylation patterns as other previously published plant methylomes, with CG ∼50%, CHG ∼30%, and CHH ∼2.7% methylation, however, these differ from the common bean reference methylome of Andean origin. We identified higher CG methylation levels in both promoter and genic regions than CHG and CHH contexts. Moreover, we found relatively higher CG methylation levels in genes than in promoters. Conversely, the CHG and CHH methylation levels were highest in promoters than in genes. This is the first genome-wide DNA methylation profiling study in a Meso-American common bean cultivar ("Sierra") using NGS approaches. Our long-term goal is to generate genome-wide epigenomic maps in common bean focusing on chromatin accessibility, histone modifications, and DNA methylation.

Global Methylation Patterns and Their Relationship with Gene Expression and Small RNA in Rice Lines with Different Ploidy

Whole genome duplication (WGD) is a major force in angiosperm evolution. Whether WGD is accompanied by the evolution of epigenetic regulators remains to be explored. Here we investigate whole genome methylation, gene expression, and miRNA regulation among monoploid, diploid, and triploid rice plants isolated from a twin-seedling population. The DNA methylation patterns in the three different ploidy plants were highly similar, with DNA methylation primarily enriched in the promoters. We examined the methylation of single genes and detected around 25,500 methylated genes, of which 22,751 were methylated in all three lines. Significantly divergent DNA methylation patterns between each pair of three lines were only detected in 64 genes, though more genes were found to exhibit differential expression. Analysis of DNA methylation and expression patterns showed that higher DNA methylation levels upstream of the transcription start sites are correlated with higher levels of expression of related genes; whereas higher DNA methylation levels in gene body regions are correlated with lower levels of expression. We also carried out high-throughput sequencing of small RNA libraries and identified 36 new miRNAs. These miRNAs have different expression levels depending on the ploidy.

Whole-genome DNA methylation profiling with nucleotide resolution

In many eukaryotic organisms, methylation at the fifth carbon of cytosine (5mC) is a stable epigenetic mark crucial for many biological processes, including cell differentiation, X-chromosome inactivation, transposon silencing, and genomic imprinting. DNA methylation can be stably inherited to the subsequent generation. It can also change dynamically in response to developmental cues or environmental stimuli, and is an important regulator for developmental switch and cell fate determination. Consequently, many human diseases are associated with aberrant DNA methylation. Gene-specific methylation analysis by sequencing of bisulfite-treated genomic DNA has been instrumental in understanding how DNA methylation affects gene transcription. In recent years, techniques have been developed for genome-wide 5mC detection, and complete methylome at single base resolution has been reported for several organisms, providing unprecedented details on the dynamic nature of DNA methylation during development. With the advance in high-throughput sequencing and the availability of genome sequences, mapping the methylome for species with complex genomes has become increasingly feasible.

Crop epigenetics and the molecular hardware of genotype × environment interactions

Crop plants encounter thermal environments which fluctuate on a diurnal and seasonal basis. Future climate resilient cultivars will need to respond to thermal profiles reflecting more variable conditions, and harness plasticity that involves regulation of epigenetic processes and complex genomic regulatory networks. Compartmentalization within plant cells insulates the genomic central processing unit within the interphase nucleus. This review addresses the properties of the chromatin hardware in which the genome is embedded, focusing on the biophysical and thermodynamic properties of DNA, histones and nucleosomes. It explores the consequences of thermal and ionic variation on the biophysical behavior of epigenetic marks such as DNA cytosine methylation (5mC), and histone variants such as H2A.Z, and how these contribute to maintenance of chromatin integrity in the nucleus, while enabling specific subsets of genes to be regulated. Information is drawn from theoretical molecular in vitro studies as well as model and crop plants and incorporates recent insights into the role epigenetic processes play in mediating between environmental signals and genomic regulation. A preliminary speculative framework is outlined, based on the evidence of what appears to be a cohesive set of interactions at molecular, biophysical and electrostatic level between the various components contributing to chromatin conformation and dynamics. It proposes that within plant nuclei, general and localized ionic homeostasis plays an important role in maintaining chromatin conformation, whilst maintaining complex genomic regulation that involves specific patterns of epigenetic marks. More generally, reversible changes in DNA methylation appear to be consistent with the ability of nuclear chromatin to manage variation in external ionic and temperature environment. Whilst tentative, this framework provides scope to develop experimental approaches to understand in greater detail the internal environment of plant nuclei. It is hoped that this will generate a deeper understanding of the molecular mechanisms underlying genotype × environment interactions that may be beneficial for long-term improvement of crop performance in less predictable climates.

Chromatin state analysis of the barley epigenome reveals a higher-order structure defined by H3K27me1 and H3K27me3 abundance

Combinations of histones carrying different covalent modifications are a major component of epigenetic variation. We have mapped nine modified histones in the barley seedling epigenome by chromatin immunoprecipitation next-generation sequencing (ChIP-seq). The chromosomal distributions of the modifications group them into four different classes, and members of a given class also tend to coincide at the local DNA level, suggesting that global distribution patterns reflect local epigenetic environments. We used this peak sharing to define 10 chromatin states representing local epigenetic environments in the barley genome. Five states map mainly to genes and five to intergenic regions. Two genic states involving H3K36me3 are preferentially associated with constitutive gene expression, while an H3K27me3-containing genic state is associated with differentially expressed genes. The 10 states display striking distribution patterns that divide barley chromosomes into three distinct global environments. First, telomere-proximal regions contain high densities of H3K27me3 covering both genes and intergenic DNA, together with very low levels of the repressive H3K27me1 modification. Flanking these are gene-rich interior regions that are rich in active chromatin states and have greatly decreased levels of H3K27me3 and increasing amounts of H3K27me1 and H3K9me2. Lastly, H3K27me3-depleted pericentromeric regions contain gene islands with active chromatin states separated by extensive retrotransposon-rich regions that are associated with abundant H3K27me1 and H3K9me2 modifications. We propose an epigenomic framework for barley whereby intergenic H3K27me3 specifies facultative heterochromatin in the telomere-proximal regions and H3K27me1 is diagnostic for constitutive heterochromatin elsewhere in the barley genome.

Global Analysis Reveals the Crucial Roles of DNA Methylation during Rice Seed Development

Seed development is an important process of reproductive development and consists of embryo and endosperm development; both comprise several key processes. To determine and investigate the functions of the dynamic DNA methylome during seed development, we profiled the DNA methylation genome wide in a series of developmental stages of rice (Oryza sativa) embryo and endosperm by methylcytosine immunoprecipitation followed by Illumina sequencing. The results showed that embryo is hypermethylated predominantly around non-transposable element (TE) genes, short DNA-TEs, and short interspersed TEs compared with endosperm, and non-TE genes have the most diverse methylation status across seed development. In addition, lowly expressed genes are significantly enriched in hypermethylated genes, but not vice versa, confirming the crucial role of DNA methylation in suppressing gene transcription. Further analysis revealed the significantly decreased methylation at early developing stages (from 2 to 3 d after pollination), indicating a predominant role of demethylation during early endosperm development and that genes with a consistent negative correlation between DNA methylation change and expression change may be potentially directly regulated by DNA methylation. Interestingly, comparative analysis of the DNA methylation profiles revealed that both rice indica and japonica subspecies showed robust fluctuant profiles of DNA methylation levels in embryo and endosperm across seed development, with the highest methylation level at 6 d after pollination (2 d after pollination of endosperm in japonica as well), indicating that a complex and finely controlled methylation pattern is closely associated with seed development regulation. The systemic characterization of the dynamic DNA methylome in developing rice seeds will help us understand the effects and mechanism of epigenetic regulation in seed development.

Lineage-specific chromatin signatures reveal a regulator of lipid metabolism in microalgae

Alga-derived lipids represent an attractive potential source of biofuels. However, lipid accumulation in algae is a stress response tightly coupled to growth arrest, thereby imposing a major limitation on productivity. To identify transcriptional regulators of lipid accumulation, we performed an integrative chromatin signature and transcriptomic analysis to decipher the regulation of lipid biosynthesis in the alga Chlamydomonas reinhardtii. Genome-wide histone modification profiling revealed remarkable differences in functional chromatin states between the algae and higher eukaryotes and uncovered regulatory components at the core of lipid accumulation pathways. We identified the transcription factor, PSR1, as a pivotal switch that triggers cytosolic lipid accumulation. Dissection of the PSR1-induced lipid profiles corroborates its role in coordinating multiple lipid-inducing stress responses. The comprehensive maps of functional chromatin signatures in a major clade of eukaryotic life and the discovery of a transcriptional regulator of algal lipid metabolism will facilitate targeted engineering strategies to mediate high lipid production in microalgae.

DNA methylation and gene expression in Mimulus guttatus

Background

The presence of methyl groups on cytosine nucleotides across an organism’s genome (methylation) is a major regulator of genome stability, crossing over, and gene regulation. The capacity for DNA methylation to be altered by environmental conditions, and potentially passed between generations, makes it a prime candidate for transgenerational epigenetic inheritance. Here we conduct the first analysis of the Mimulus guttatus methylome, with a focus on the relationship between DNA methylation and gene expression.

Results

We present a whole genome methylome for the inbred line Iron Mountain 62 (IM62). DNA methylation varies across chromosomes, genomic regions, and genes. We develop a model that predicts gene expression based on DNA methylation (R² = 0.2). Post hoc analysis of this model confirms prior relationships, and identifies novel relationships between methylation and gene expression. Additionally, we find that DNA methylation is significantly depleted near gene transcriptional start sites, which may explain the recently discovered elevated rate of recombination in these same regions.

Conclusions

The establishment here of a reference methylome will be a useful resource for the continued advancement of M. guttatus as a model system. Using a model-based approach, we demonstrate that methylation patterns are an important predictor of variation in gene expression. This model provides a novel approach for differential methylation analysis that generates distinct and testable hypotheses regarding gene expression.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1668-0) contains supplementary material, which is available to authorized users.

The roles of cross-talk epigenetic patterns in Arabidopsis thaliana

The epigenetic mechanisms, including histone modifications, DNA cytosine methylation, histone variants and noncoding RNAs (ncRNAs), play a key role in determining transcriptional outcomes. Recently, many studies have demonstrated that the different epigenetic mechanisms interplay with each other rather than work independently. In this article, we outline a framework for how different epigenetic mechanisms work with each other in Arabidopsis thalianaWe build a network of cross-talk between chromatin marks based on six classes of cross-talk interactions. The first pattern details coordinated modifications that act together to enhance or repress gene expression. The second pattern details bivalent modifications that act antagonistically toward gene expression. The third pattern is for unilateral promotion of one modification by the existence of another modification. The fourth pattern is for unilateral inhibition of one modification by another modification. The fifth pattern is for mutual inhibitory patterns. The sixth pattern is for epigenetic modifications that appear independent.We also explore the mutual regulation between chromatin marks and ncRNAs in various ways. These regulations can be divided into six parts: how ncRNA affects the binding of chromatin mark, such as miR2Epi, siR2Epi and lncR2Epi; how chromatin mark regulates ncRNA, such as Epi2miR, Epi2siR and Epi2lncR.A comprehensive network of cross-talk between different epigenetic mechanisms will help in fully understanding the functional roles and biological impacts of epigenetic regulation.

Conserved Epigenetic Mechanisms Could Play a Key Role in Regulation of Photosynthesis and Development-Related Genes during Needle Development of Pinus radiata

Needle maturation is a complex process that involves cell growth, differentiation and tissue remodelling towards the acquisition of full physiological competence. Leaf induction mechanisms are well known; however, those underlying the acquisition of physiological competence are still poorly understood, especially in conifers. We studied the specific epigenetic regulation of genes defining organ function (PrRBCS and PrRBCA) and competence and stress response (PrCSDP2 and PrSHMT4) during three stages of needle development and one de-differentiated control. Gene-specific changes in DNA methylation and histone were analysed by bisulfite sequencing and chromatin immunoprecipitation (ChIP). The expression of PrRBCA and PrRBCS increased during needle maturation and was associated with the progressive loss of H3K9me3, H3K27me3 and the increase in AcH4. The maturation-related silencing of PrSHMT4 was correlated with increased H3K9me3 levels, and the repression of PrCSDP2, to the interplay between AcH4, H3K27me3, H3K9me3 and specific DNA methylation. The employ of HAT and HDAC inhibitors led to a further determination of the role of histone acetylation in the regulation of our target genes. The integration of these results with high-throughput analyses in Arabidopsis thaliana and Populus trichocarpa suggests that the specific epigenetic mechanisms that regulate photosynthetic genes are conserved between the analysed species.

Genome-wide mapping of histone H3 lysine 4 trimethylation in Eucalyptus grandis developing xylem

BACKGROUND

Histone modifications play an integral role in plant development, but have been poorly studied in woody plants. Investigating chromatin organization in wood-forming tissue and its role in regulating gene expression allows us to understand the mechanisms underlying cellular differentiation during xylogenesis (wood formation) and identify novel functional regions in plant genomes. However, woody tissue poses unique challenges for using high-throughput chromatin immunoprecipitation (ChIP) techniques for studying genome-wide histone modifications in vivo. We investigated the role of the modified histone H3K4me3 (trimethylated lysine 4 of histone H3) in gene expression during the early stages of wood formation using ChIP-seq in Eucalyptus grandis, a woody biomass model.

RESULTS

Plant chromatin fixation and isolation protocols were optimized for developing xylem tissue collected from field-grown E. grandis trees. A "nano-ChIP-seq" procedure was employed for ChIP DNA amplification. Over 9 million H3K4me3 ChIP-seq and 18 million control paired-end reads were mapped to the E. grandis reference genome for peak-calling using Model-based Analysis of ChIP-Seq. The 12,177 significant H3K4me3 peaks identified covered ~1.5% of the genome and overlapped some 9,623 protein-coding genes and 38 noncoding RNAs. H3K4me3 library coverage, peaking ~600 - 700 bp downstream of the transcription start site, was highly correlated with gene expression levels measured with RNA-seq. Overall, H3K4me3-enriched genes tended to be less tissue-specific than unenriched genes and were overrepresented for general cellular metabolism and development gene ontology terms. Relative expression of H3K4me3-enriched genes in developing secondary xylem was higher than unenriched genes, however, and highly expressed secondary cell wall-related genes were enriched for H3K4me3 as validated using ChIP-qPCR.

CONCLUSIONS

In this first genome-wide analysis of a modified histone in a woody tissue, we optimized a ChIP-seq procedure suitable for field-collected samples. In developing E. grandis xylem, H3K4me3 enrichment is an indicator of active transcription, consistent with its known role in sustaining pre-initiation complex formation in yeast. The H3K4me3 ChIP-seq data from this study paves the way to understanding the chromatin landscape and epigenomic architecture of xylogenesis in plants, and complements RNA-seq evidence of gene expression for the future improvement of the E. grandis genome annotation.

MethylC-seq library preparation for base-resolution whole-genome bisulfite sequencing

Current high-throughput DNA sequencing technologies enable acquisition of billions of data points through which myriad biological processes can be interrogated, including genetic variation, chromatin structure, gene expression patterns, small RNAs and protein-DNA interactions. Here we describe the MethylC-sequencing (MethylC-seq) library preparation method, a 2-d protocol that enables the genome-wide identification of cytosine DNA methylation states at single-base resolution. The technique involves fragmentation of genomic DNA followed by adapter ligation, bisulfite conversion and limited amplification using adapter-specific PCR primers in preparation for sequencing. To date, this protocol has been successfully applied to genomic DNA isolated from primary cell culture, sorted cells and fresh tissue from over a thousand plant and animal samples.