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

Actively Expressed Intergenic Genes Generated by Transposable Element Insertions in Gossypium hirsutum Cotton

The genomes and annotated genes of allotetraploid cotton Gossypium hirsutum have been extensively studied in recent years. However, the expression, regulation, and evolution of intergenic genes (ITGs) have not been completely deciphered. In this study, we identified a novel set of actively expressed ITGs in G. hirsutum cotton, through transcriptome profiling based on deep sequencing data, as well as chromatin immunoprecipitation, followed by sequencing (ChIP-seq) of histone modifications and how the ITGs evolved. Totals of 17,567 and 8249 ITGs were identified in G. hirsutum and Gossypium arboreum, respectively. The expression of ITGs in G. hirsutum was significantly higher than that in G. arboreum. Moreover, longer exons were observed in G. hirsutum ITGs. Notably, 42.3% of the ITGs from G. hirsutum were generated by the long terminal repeat (LTR) insertions, while their proportion in genic genes was 19.9%. The H3K27ac and H3K4me3 modification proportions and intensities of ITGs were equivalent to genic genes. The H3K4me1 modifications were lower in ITGs. Additionally, evolution analyses revealed that the ITGs from G. hirsutum were mainly produced around 6.6 and 1.6 million years ago (Mya), later than the pegged time for genic genes, which is 7.0 Mya. The characterization of ITGs helps to elucidate the evolution of cotton genomes and shed more light on their biological functions in the transcriptional regulation of eukaryotic genes, along with the roles of histone modifications in speciation and diversification.

Whole-genome landscape of histone H3K4me3 modification during sperm cell lineage development in tomato

BACKGROUND

During male gametogenesis of flowering plants, sperm cell lineage (microspores, generative cells, and sperm cells) differentiated from somatic cells and acquired different cell fates. Trimethylation of histone H3 on lysine 4 (H3K4me3) epigenetically contributes to this process, however, it remained unclear how H3K4me3 influences the gene expression in each cell type. Here, we conducted chromatin immunoprecipitation sequencing (ChIP-seq) to obtain a genome-wide landscape of H3K4me3 during sperm cell lineage development in tomato (Solanum lycopersicum).

RESULTS

We show that H3K4me3 peaks were mainly enriched in the promoter regions, and intergenic H3K4me3 peaks expanded as sperm cell lineage differentiated from somatic cells. H3K4me3 was generally positively associated with transcript abundance and served as a better indicator of gene expression in somatic and vegetative cells, compared to sperm cell lineage. H3K4me3 was mutually exclusive with DNA methylation at 3' proximal of the transcription start sites. The microspore maintained the H3K4me3 features of somatic cells, while generative cells and sperm cells shared an almost identical H3K4me3 pattern which differed from that of the vegetative cell. After microspore division, significant loss of H3K4me3 in genes related to brassinosteroid and cytokinin signaling was observed in generative cells and vegetative cells, respectively.

CONCLUSIONS

Our results suggest the asymmetric division of the microspore significantly reshapes the genome-wide distribution of H3K4me3. Selective loss of H3K4me3 in genes related to hormone signaling may contribute to functional differentiation of sperm cell lineage. This work provides new resource data for the epigenetic studies of gametogenesis in plants.

Rice responds to Spodoptera frugiperda infestation via epigenetic regulation of H3K9ac in the jasmonic acid signaling and phenylpropanoid biosynthesis pathways

KEY MESSAGE

This study provided important insights into the complex epigenetic regulatory of H3K9ac-modified genes involved in the jasmonic acid signaling and phenylpropanoid biosynthesis pathways of rice in response to Spodoptera frugiperda infestation. Physiological and molecular mechanisms underlying plant responses to insect herbivores have been well studied, while epigenetic modifications such as histone acetylation and their potential regulation at the genomic level of hidden genes remain largely unknown. Histone 3 lysine 9 acetylation (H3K9ac) is an epigenetic marker widely distributed in plants that can activate gene transcription. In this study, we provided the genome-wide profiles of H3K9ac in rice (Oryza sativa) infested by fall armyworm (Spodoptera frugiperda, FAW) using CUT&Tag-seq and RNA-seq. There were 3269 and 4609 up-regulated genes identified in plants infested by FAW larvae for 3 h and 12 h, respectively, which were mainly enriched in alpha-linolenic acid and phenylpropanoid pathways according to transcriptomic analysis. In addition, CUT&Tag-seq analysis revealed increased H3K9ac in FAW-infested plants, and there were 422 and 543 up-regulated genes enriched with H3K9ac observed at 3 h and 12 h after FAW feeding, respectively. Genes with increased H3K9ac were mainly enriched in the transcription start site (TSS), suggesting that H3K9ac is related to gene transcription. Integrative analysis of both RNA-seq and CUT&Tag-seq data showed that up-expressed genes with H3K9ac enrichment were mainly involved in the jasmonic acid (JA) and phenylpropanoid pathways. Particularly, two spermidine hydroxycinnamoyl transferase genes SHT1 and SHT2 involved in phenolamide biosynthesis were highly modified by H3K9ac in FAW-infested plants. Furthermore, the Ossht1 and Ossht2 transgenic lines exhibited decreased resistance against FAW larvae. Our findings suggest that rice responds to insect herbivory via H3K9ac epigenetic regulation in the JA signaling and phenolamide biosynthesis pathways.

Integrating ATAC-seq and RNA-seq to identify differentially expressed genes with chromatin-accessible changes during photosynthetic establishment in Populus leaves

Leaves are the primary photosynthetic organs, providing essential substances for tree growth. It is important to obtain an anatomical understanding and regulatory network analysis of leaf development. Here, we studied leaf development in Populus Nanlin895 along a development gradient from the newly emerged leaf from the shoot apex to the sixth leaf (L1 to L6) using anatomical observations and RNA-seq analysis. It indicated that mesophyll cells possess obvious vascular, palisade, and spongy tissue with distinct intercellular spaces after L3. Additionally, vacuoles fuse while epidermal cells expand to form pavement cells. RNA-seq analysis indicated that genes highly expressed in L1 and L2 were related to cell division and differentiation, while those highly expressed in L3 were enriched in photosynthesis. Therefore, we selected L1 and L3 to integrate ATAC-seq and RNA-seq and identified 735 differentially expressed genes (DEGs) with changes in chromatin accessibility regions within their promoters, of which 87 were transcription factors (TFs), such as ABI3VP1, AP-EREBP, MYB, NAC, and GRF. Motif enrichment analysis revealed potential regulatory functions for the DEGs through upstream TFs including TCP, bZIP, HD-ZIP, Dof, BBR-BPC, and MYB. Overall, our research provides a potential molecular foundation for regulatory network exploration in leaf development during photosynthesis establishment.

Chromosome-Wide Distribution and Characterization of H3K36me3 and H3K27Ac in the Marine Model Diatom Phaeodactylum tricornutum

Histone methylation and acetylation play a crucial role in response to developmental cues and environmental changes. Previously, we employed mass spectrometry to identify histone modifications such as H3K27ac and H3K36me3 in the model diatom Phaeodactylum tricornutum, which have been shown to be important for transcriptional activation in animal and plant species. To further investigate their evolutionary implications, we utilized chromatin immunoprecipitation followed by deep sequencing (ChIP-Seq) and explored their genome-wide distribution in P. tricornutum. Our study aimed to determine their role in transcriptional regulation of genes and transposable elements (TEs) and their co-occurrence with other histone marks. Our results revealed that H3K27ac and H3K36me3 were predominantly localized in promoters and genic regions indicating a high conservation pattern with studies of the same marks in plants and animals. Furthermore, we report the diversity of genes encoding H3 lysine 36 (H3K36) trimethylation-specific methyltransferase in microalgae leveraging diverse sequencing resources including the Marine Microbial Eukaryote Transcriptome Sequencing Project database (MMETSP). Our study expands the repertoire of epigenetic marks in a model microalga and provides valuable insights into the evolutionary context of epigenetic-mediated gene regulation. These findings shed light on the intricate interplay between histone modifications and gene expression in microalgae, contributing to our understanding of the broader epigenetic landscape in eukaryotic organisms.

Epigenetic Landscape Is Largely Shaped by Diversiform Transposons in Aegilops tauschii

Transposons (TEs) account for more than 80% of the wheat genome, the highest among all known crop species. They play an important role in shaping the elaborate genomic landscape, which is the key to the speciation of wheat. In this study, we analyzed the association between TEs, chromatin states, and chromatin accessibility in Aegilops tauschii, the D genome donor of bread wheat. We found that TEs contributed to the complex but orderly epigenetic landscape as chromatin states showed diverse distributions on TEs of different orders or superfamilies. TEs also contributed to the chromatin state and openness of potential regulatory elements, affecting the expression of TE-related genes. Some TE superfamilies, such as hAT-Ac, carry active/open chromatin regions. In addition, the histone mark H3K9ac was found to be associated with the accessibility shaped by TEs. These results suggest the role of diversiform TEs in shaping the epigenetic landscape and in gene expression regulation in Aegilops tauschii. This has positive implications for understanding the transposon roles in Aegilops tauschii or the wheat D genome.

Global Decrease in H3K9 Acetylation in Sorghum Seed Postgermination Stages

Sorghum seed germination is accompanied by increases in nutrient contents and reduced levels of antinutrients and is therefore being applied to food processing. However, the characterization of acetylated histone H3 at lysine residue 9 (H3K9ac) in sorghum postgermination has lagged. In this study, we performed chromatin immunoprecipitation sequencing (ChIP-seq) to identify H3K9ac enrichment and obtained transcriptome in postgermination stages. More than 10,000 hypoacetylated genes gained H3K9ac marks in the postgermination stages. In addition, we observed that the expression of the main histone deacetylase (HDAC) genes was elevated. The application of the HDAC inhibitor trichostatin A (TSA) resulted in seed growth arrest, suggesting that the repression of the H3K9ac modification is critical for postgermination. Additionally, we obtained a comprehensive view of abundant genomic changes in H3K9ac-marked regions and transcription between the mock and TSA treatment groups, which suggested that H3K9ac was required in the late stage of autotrophic seedling establishment. Metabolic profiling, transcriptome analyses, and ChIP-seq revealed that H3K9ac is enriched at genes involved in phenylpropanoid, including lignin and flavonoid, biosynthesis. Our results suggest important roles of H3K9ac in sorghum seed postgermination stages.

Comprehensive Survey of ChIP-Seq Datasets to Identify Candidate Iron Homeostasis Genes Regulated by Chromatin Modifications

Vital biochemical reactions including photosynthesis to respiration require iron, which should be tightly regulated. Although increasing evidence reveals the importance of epigenetic regulation in gene expression and signaling, the role of histone modifications and chromatin remodeling in plant iron homeostasis is not well understood. In this study, we surveyed publicly available ChIP-seq datasets of Arabidopsis wild-type and mutants defective in key enzymes of histone modification and chromatin remodeling and compared the deposition of epigenetic marks on loci of genes involved in iron regulation. Based on the analysis, we compiled a comprehensive list of iron homeostasis genes with differential enrichment of various histone modifications. This report will provide a resource for future studies to investigate epigenetic regulatory mechanisms of iron homeostasis in plants.

Inferring regulatory element landscapes and gene regulatory networks from integrated analysis in eight hulless barley varieties under abiotic stress

BACKGROUND

The cis-regulatory element became increasingly important for resistance breeding. There were many DNA variations identified by resequencing. To investigate the links between the DNA variations and cis-regulatory element was the fundamental work. DNA variations in cis-regulatory elements caused phenotype variations in general.

RESULTS

We used WGBS, ChIP-seq and RNA-seq technology to decipher the regulatory element landscape from eight hulless barley varieties under four kinds of abiotic stresses. We discovered 231,440 lowly methylated regions (LMRs) from the methylome data of eight varieties. The LMRs mainly distributed in the intergenic regions. A total of 97,909 enhancer-gene pairs were identified from the correlation analysis between methylation degree and expression level. A lot of enriched motifs were recognized from the tolerant-specific LMRs. The key transcription factors were screened out and the transcription factor regulatory network was inferred from the enhancer-gene pairs data for drought stress. The NAC transcription factor was predicted to target to TCP, bHLH, bZIP transcription factor genes. We concluded that the H3K27me3 modification regions overlapped with the LMRs more than the H3K4me3. The variation of single nucleotide polymorphism was more abundant in LMRs than the remain regions of the genome.

CONCLUSIONS

Epigenetic regulation is an important mechanism for organisms to adapt to complex environments. Through the study of DNA methylation and histone modification, we found that many changes had taken place in enhancers and transcription factors in the abiotic stress of hulless barley. For example, transcription factors including NAC may play an important role. This enriched the molecular basis of highland barley stress response.

Induction of Somatic Embryogenesis in Plants: Different Players and Focus on WUSCHEL and WUS-RELATED HOMEOBOX (WOX) Transcription Factors

In plants, other cells can express totipotency in addition to the zygote, thus resulting in embryo differentiation; this appears evident in apomictic and epiphyllous plants. According to Haberlandt's theory, all plant cells can regenerate a complete plant if the nucleus and the membrane system are intact. In fact, under in vitro conditions, ectopic embryos and adventitious shoots can develop from many organs of the mature plant body. We are beginning to understand how determination processes are regulated and how cell specialization occurs. However, we still need to unravel the mechanisms whereby a cell interprets its position, decides its fate, and communicates it to others. The induction of somatic embryogenesis might be based on a plant growth regulator signal (auxin) to determine an appropriate cellular environment and other factors, including stress and ectopic expression of embryo or meristem identity transcription factors (TFs). Still, we are far from having a complete view of the regulatory genes, their target genes, and their action hierarchy. As in animals, epigenetic reprogramming also plays an essential role in re-establishing the competence of differentiated cells to undergo somatic embryogenesis. Herein, we describe the functions of WUSCHEL-RELATED HOMEOBOX (WOX) transcription factors in regulating the differentiation-dedifferentiation cell process and in the developmental phase of in vitro regenerated adventitious structures.

Transcriptional Stress Memory and Transgenerational Inheritance of Drought Tolerance in Plants

Plants respond to drought stress by producing abscisic acid, a chemical messenger that regulates gene expression and thereby expedites various physiological and cellular processes including the stomatal operation to mitigate stress and promote tolerance. To trigger or suppress gene transcription under drought stress conditions, the surrounding chromatin architecture must be converted between a repressive and active state by epigenetic remodeling, which is achieved by the dynamic interplay among DNA methylation, histone modifications, loop formation, and non-coding RNA generation. Plants can memorize chromatin status under drought conditions to enable them to deal with recurrent stress. Furthermore, drought tolerance acquired during plant growth can be transmitted to the next generation. The epigenetically modified chromatin architectures of memory genes under stressful conditions can be transmitted to newly developed cells by mitotic cell division, and to germline cells of offspring by overcoming the restraints on meiosis. In mammalian cells, the acquired memory state is completely erased and reset during meiosis. The mechanism by which plant cells overcome this resetting during meiosis to transmit memory is unclear. In this article, we review recent findings on the mechanism underlying transcriptional stress memory and the transgenerational inheritance of drought tolerance in plants.

Genome-wide H3K9 acetylation level increases with age-dependent senescence of flag leaf in rice

Flag leaf senescence is an important biological process that drives the remobilization of nutrients to the growing organs of rice. Leaf senescence is controlled by genetic information via gene expression and histone modification, but the precise mechanism is as yet unclear. Here, we analysed genome-wide acetylated lysine residue 9 of histone H3 (H3K9ac) enrichment by chromatin immunoprecipitation-sequencing (ChIP-seq), and examined its association with transcriptomes by RNA-seq during flag leaf aging in rice (Oryza sativa). We found that genome-wide H3K9 acetylation levels increased with age-dependent senescence in rice flag leaf, and there was a positive correlation between the density and breadth of H3K9ac with gene expression and transcript elongation. During flag leaf aging, we observed 1249 up-regulated differentially expressed genes (DEGs) and 996 down-regulated DEGs, showing a strong relationship between temporal changes in gene expression and gain/loss of H3K9ac. We produced a landscape of H3K9 acetylation-modified gene expression targets that include known senescence-associated genes, metabolism-related genes, as well as miRNA biosynthesis-related genes. Our findings reveal a complex regulatory network of metabolism- and senescence-related pathways mediated by H3K9ac, and elucidate patterns of H3K9ac-mediated regulation of gene expression during flag leaf aging in rice.

Characterization of Histone H3 Gene Family Reveals That GmHH3-3 is Associated With Higher Seed Weight in Glycine max

The main function of histone protein is to provide support to the structure of chromosomes. It helps in binding a long thread of DNA into a more condensed shape to fit into the nucleus. From histone variants, histone H3 (HH3) plays a crucial role in plant growth and development. Characterization of histones has not been reported in Glycine max till now. The objective of this study was to characterize the HH3 gene family for molecular breeding of G. max. In this study, 17 HH3 members in G. max were identified by performing local BLASTp using HH3 members from Arabidopsis as a query. Phylogenetic analysis classified HH3 genes in seven clades. Sequence logo analysis among Arabidopsis thaliana, Oryza sativa, and Glycine max showed a higher level of similarity in amino acids. Furthermore, conserveness of G. max HH3 genes was also confirmed by Gene Structure Display. Ten paralogous gene pairs were identified in GmHH3 genes in the Glycine max genome by conducting collinearity analysis. G. max HH3 genes have experienced strong purifying selection pressure, with limited functional divergence originating from the segmental and whole-genome duplication, as evidenced by the Ka/Ks ratio. The KASP marker was developed for GmHH3-3 gene. Genotyping was performed on 46 G. max genotypes. This differentiation was based upon the presence of either GmHH3-3-C or GmHH3-3-T allele in the CDS region. The results showed that G. max accessions containing the GmHH3-3-T allele at respective locus showed higher thousand seed weight than that of those accessions that contain the GmHH3-3-C allele. This research provides the basic information to further decipher the function of HH3 in soybean.

Dynamicity of histone H3K27ac and H3K27me3 modifications regulate the cold-responsive gene expression in Oryza sativa L. ssp. indica

Cultivated in tropical and subtropical regions, Oryza sativa L. ssp. indica is largely affected by cold-stress, especially at the seedling stage. The present model of the stress-responsive regulatory network in plants entails the role of genetic and epigenetic factors in stress-responsive gene expression. Despite extensive transcriptomic studies, the regulation of various epigenetic factors in plants cold-stress response is less explored. The present study addresses the effect of genome-wide changes of H3K27 modifications on gene expression in IR64 rice, during cold-stress. Our results suggest a positive correlation between the changes in H3K27 modifications and stress-responsive gene activation in indica rice. Cold-induced enrichment of H3K27 acetylation promotes nucleosomal rearrangement, thereby facilitating the accessibility of the transcriptional machinery at the stress-responsive loci for transcription activation. Although H3K27ac exhibits uniform distribution throughout the loci of enriched genes; occupancy of H3K27me3 is biased to intergenic regions. Integration of the ChIP-seq data with transcriptome indicated that upregulation of stress-responsive TFs, photosynthesis-TCA-related, water-deficit genes, redox and JA signalling components, was associated with differential changes of H3K27ac and H3K27me3 levels. Furthermore, cold-induced upregulation of histone acetyltransferases and downregulation of DNA methyltransferases was noted through the antagonistic switch of H3K27ac and H3K27me3. Moreover, motif analysis of H3K27ac and H3K27me3 enriched regions are associated with putative stress responsive transcription factors binding sites, GAGA element and histone H3K27demethylase. Collectively our analysis suggests that differential expression of various chromatin and DNA modifiers coupled with increased H3K27ac and depleted H3K27me3 increases DNA accessibility, thereby promoting transcription of the cold-responsive genes in indica rice.

Twenty years of rice genomics research: From sequencing and functional genomics to quantitative genomics

Since the completion of the rice genome sequencing project in 2005, we have entered the era of rice genomics, which is still in its ascendancy. Rice genomics studies can be classified into three stages: structural genomics, functional genomics, and quantitative genomics. Structural genomics refers primarily to genome sequencing for the construction of a complete map of rice genome sequence. This is fundamental for rice genetics and molecular biology research. Functional genomics aims to decode the functions of rice genes. Quantitative genomics is large-scale sequence- and statistics-based research to define the quantitative traits and genetic features of rice populations. Rice genomics has been a transformative influence on rice biological research and contributes significantly to rice breeding, making rice a good model plant for studying crop sciences.

Genomic variants affecting homoeologous gene expression dosage contribute to agronomic trait variation in allopolyploid wheat

Allopolyploidy greatly expands the range of possible regulatory interactions among functionally redundant homoeologous genes. However, connection between the emerging regulatory complexity and expression and phenotypic diversity in polyploid crops remains elusive. Here, we use diverse wheat accessions to map expression quantitative trait loci (eQTL) and evaluate their effects on the population-scale variation in homoeolog expression dosage. The relative contribution of cis- and trans-eQTL to homoeolog expression variation is strongly affected by both selection and demographic events. Though trans-acting effects play major role in expression regulation, the expression dosage of homoeologs is largely influenced by cis-acting variants, which appear to be subjected to selection. The frequency and expression of homoeologous gene alleles showing strong expression dosage bias are predictive of variation in yield-related traits, and have likely been impacted by breeding for increased productivity. Our study highlights the importance of genomic variants affecting homoeolog expression dosage in shaping agronomic phenotypes and points at their potential utility for improving yield in polyploid crops.

Hierarchical chromatin features reveal the toxin production in Bungarus multicinctus

BACKGROUND

Bungarus multicinctus, from which a classical Chinese medicine is produced, is known as the most venomous land snake in the world, but the chromatin organization and transcription factor activity during venom replenishment progress have not been explored yet. This study aimed to determine the roles of chromatin structure in toxin activity via bioinformatics and experimental validation.

METHODS

Chromosome conformation capture (Hi-C) analysis was used to examine interactions among chromosomes and identify different scales of chromatin during envenomation in B. multicinctus. Correlations between epigenetic modifications and chromatin structure were verified through ChIP-seq analysis. RNA-seq was used to validate the influence of variation in chromatin structure and gene expression levels on venom production and regulation.

RESULTS

Our results suggested that intra-chromosomal interactions are more intense than inter-chromosomal interactions among the control group, 3-day group of venom glands and muscles. Through this, we found that compartmental transition was correlated with chromatin interactions. Interestingly, the up-regulated genes in more compartmental switch regions reflect the function of toxin activity. Topologically associated domain (TAD) boundaries enriched with histone modifications are associated with different distributions of genes and the expression levels. Toxin-coding genes in the same loop are highly expressed, implying that the importance of epigenetic regulation during envenomination. On a smaller scale, the epigenetic markers affect transcriptional regulation by controlling the recruitment/inhibition of transcription initiation complexes.

CONCLUSIONS

Chromatin structure and epigenetic modifications could play a vital status role in the mechanisms of venom regulation in B. multicinctus.

H3.1K27me1 maintains transcriptional silencing and genome stability by preventing GCN5-mediated histone acetylation

Epigenetic mechanisms play diverse roles in the regulation of genome stability in eukaryotes. In Arabidopsis thaliana, genome stability is maintained during DNA replication by the H3.1K27 methyltransferases ARABIDOPSIS TRITHORAX-RELATED PROTEIN 5 (ATXR5) and ATXR6, which catalyze the deposition of K27me1 on replication-dependent H3.1 variants. The loss of H3.1K27me1 in atxr5 atxr6 double mutants leads to heterochromatin defects, including transcriptional de-repression and genomic instability, but the molecular mechanisms involved remain largely unknown. In this study, we identified the transcriptional co-activator and conserved histone acetyltransferase GCN5 as a mediator of transcriptional de-repression and genomic instability in the absence of H3.1K27me1. GCN5 is part of a SAGA-like complex in plants that requires the GCN5-interacting protein ADA2b and the chromatin remodeler CHR6 to mediate the heterochromatic defects in atxr5 atxr6 mutants. Our results also indicate that Arabidopsis GCN5 acetylates multiple lysine residues on H3.1 variants, but H3.1K27 and H3.1K36 play essential functions in inducing genomic instability in the absence of H3.1K27me1. Finally, we show that H3.1K36 acetylation by GCN5 is negatively regulated by H3.1K27me1 in vitro. Overall, this work reveals a key molecular role for H3.1K27me1 in maintaining transcriptional silencing and genome stability in heterochromatin by restricting GCN5-mediated histone acetylation in plants.

Dynamic histone acetylation in floral volatile synthesis and emission in petunia flowers

Biosynthesis of secondary metabolites relies on primary metabolic pathways to provide precursors, energy, and cofactors, thus requiring coordinated regulation of primary and secondary metabolic networks. However, to date, it remains largely unknown how this coordination is achieved. Using Petunia hybrida flowers, which emit high levels of phenylpropanoid/benzenoid volatile organic compounds (VOCs), we uncovered genome-wide dynamic deposition of histone H3 lysine 9 acetylation (H3K9ac) during anthesis as an underlying mechanism to coordinate primary and secondary metabolic networks. The observed epigenome reprogramming is accompanied by transcriptional activation at gene loci involved in primary metabolic pathways that provide precursor phenylalanine, as well as secondary metabolic pathways to produce volatile compounds. We also observed transcriptional repression among genes involved in alternative phenylpropanoid branches that compete for metabolic precursors. We show that GNAT family histone acetyltransferase(s) (HATs) are required for the expression of genes involved in VOC biosynthesis and emission, by using chemical inhibitors of HATs, and by knocking down a specific HAT gene, ELP3, through transient RNAi. Together, our study supports that regulatory mechanisms at chromatin level may play an essential role in activating primary and secondary metabolic pathways to regulate VOC synthesis in petunia flowers.

Refinement of bamboo genome annotations through integrative analyses of transcriptomic and epigenomic data

Bamboo, one of the most crucial nontimber forest resources worldwide, has the capacity for rapid growth. In recent years, the genome of moso bamboo (Phyllostachys edulis) has been decoded, and a large amount of transcriptome data has been published. In this study, we generated the genome-wide profiles of the histone modification H3K4me3 in leaf, stem, and root tissues of bamboo. The trends in the distribution patterns were similar to those in rice. We developed a processing pipeline for predicting novel transcripts to refine the structural annotation of the genome using H3K4me3 ChIP-seq data and 29 RNA-seq datasets. As a result, 12,460 novel transcripts were predicted in the bamboo genome. Compared with the transcripts in the newly released version 2.0 of the bamboo genome, these novel transcripts are tissue-specific and shorter, and most have a single exon. Some representative novel transcripts were validated by semiquantitative RT-PCR and qRT-PCR analyses. Furthermore, we put these novel transcripts back into the ChIP-seq analysis pipeline and discovered that the percentages of H3K4me3 in genic elements were increased. Overall, this work integrated transcriptomic data and epigenomic data to refine the annotation of the genome in order to discover more functional genes and study bamboo growth and development, and the application of this predicted pipeline may help refine the structural annotation of the genome in other species.

Genome-Wide Mapping of Histone H3 Lysine 4 Trimethylation (H3K4me3) and Its Involvement in Fatty Acid Biosynthesis in Sunflower Developing Seeds

Histone modifications are of paramount importance during plant development. Investigating chromatin remodeling in developing oilseeds sheds light on the molecular mechanisms controlling fatty acid metabolism and facilitates the identification of new functional regions in oil crop genomes. The present study characterizes the epigenetic modifications H3K4me3 in relationship with the expression of fatty acid-related genes and transcription factors in developing sunflower seeds. Two master transcriptional regulators identified in this analysis, VIV1 (homologous to Arabidopsis ABI3) and FUS3, cooperate in the regulation of WRINKLED 1, a transcriptional factor regulating glycolysis, and fatty acid synthesis in developing oilseeds.

Asymmetric epigenome maps of subgenomes reveal imbalanced transcription and distinct evolutionary trends in Brassica napus

The complexity of the epigenome landscape and transcriptional regulation is significantly increased during plant polyploidization, which drives genome evolution and contributes to the increased adaptability to diverse environments. However, a comprehensive epigenome map of Brassica napus is still unavailable. In this study, we performed integrative analysis of five histone modifications, RNA polymerase II occupancy, DNA methylation, and transcriptomes in two B. napus lines (2063A and B409), and established global maps of regulatory elements, chromatin states, and their dynamics for the whole genome (including the An and Cn subgenomes) in four tissue types (young leaf, flower bud, silique, and root) of these two lines. Approximately 65.8% of the genome was annotated with different epigenomic signals. Compared with the Cn subgenome, the An subgenome possesses a higher level of active epigenetic marks and lower level of repressive epigenetic marks. Genes from subgenome-unique regions contribute to the major differences between the An and Cn subgenomes. Asymmetric histone modifications between homeologous gene pairs reflect their biased expression patterns. We identified a novel bivalent chromatin state (with H3K4me1 and H3K27me3) in B. napus that is associated with tissue-specific gene expression. Furthermore, we observed that different types of duplicated genes have discrepant patterns of histone modification and DNA methylation levels. Collectively, our findings provide a valuable epigenetic resource for allopolyploid plants.

Genome-wide shifts in histone modifications at early stage of rice infection with Meloidogyne graminicola

Epigenetic processes play a crucial role in the regulation of plant stress responses, but their role in plant-pathogen interactions remains poorly understood. Although histone-modifying enzymes have been observed to be deregulated in galls induced by root-knot nematodes (RKN, Meloidogyne graminicola) in rice, their influence on plant defence and their genome-wide impact has not been comprehensively investigated. First, the role of histone modifications in plant-nematode interactions was confirmed by pharmacological inhibition of histone-modifying enzymes, which all significantly affected rice susceptibility to RKN. For a more specific view, three histone marks, H3K9ac, H3K9me2, and H3K27me3, were subsequently studied by chromatin-immunoprecipitation-sequencing on RKN-induced galls at 3 days postinoculation. While levels of H3K9ac and H3K27me3 were strongly enriched, H3K9me2 was generally depleted in galls versus control root tips. Differential histone peaks were generally associated with plant defence-related genes. Transcriptome analysis using RNA-Seq and RT-qPCR-based validation revealed that genes marked with H3K9ac or H3K9me2 showed the expected activation or repression gene expression pattern, but this was not the case for H3K27me3 marks. Our results indicate that histone modifications respond dynamically to RKN infection, and that posttranslational modifications mainly at H3K9 specifically target plant defence-related genes.

Characterization of the GRAS gene family reveals their contribution to the high adaptability of wheat

GRAS transcription factors play important roles in many processes of plant development as well as abiotic and biotic stress responses. However, little is known about this gene family in bread wheat (Triticum aestivum), one of the most important crops worldwide. The completion of a quality draft genome allows genome-wide detection and evolutionary analysis of the GRAS gene family in wheat. In this study, 188 TaGRAS genes were detected and divided into 12 subfamilies based on phylogenetic analyses: DELLA, DLT, HAM, LISCL, SCL3, SCL4/7, SCR, SHR, PAT1, Os19, Os4 and LAS. Tandem and segmental duplications are the main contributors to the expansion of TaGRAS, which may contribute to the adaptation of wheat to various environmental conditions. A high rate of homoeolog retention during hexaploidization was detected, suggesting the nonredundancy and biological importance of TaGRAS homoeologs. Systematic analyses of TaGRAS indicated the conserved expression pattern and function of the same subfamily during evolution. In addition, we detected five genes belonging to the LISCL subfamily induced by both biotic and abiotic stresses and they may be potential targets for further research through gene editing. Using degradome and ChIP-seq data, we identified the targets of miR171 and histone modifications and further analyzed the contribution of epigenetic modification to the subfunctionalization of TaGRAS. This study laid a foundation for further functional elucidation of TaGRAS genes.

Genome-Wide Identification and Characterization of Main Histone Modifications in Sorghum Decipher Regulatory Mechanisms Involved by mRNA and Long Noncoding RNA Genes

Post-translational modifications of histones play an important chromatic role of a transcript activity in eukaryotes. Even though mRNA and long noncoding RNA (lncRNA) genes share similar biogenetic processes, these transcript classes may differ in many ways. However, knowledge about the crosstalk between histone methylations and the two types of sorghum genes is still ambiguous. In the present study, we reveal the genome-wide distribution of six histone modifications, namely, di- and trimethylation of H3K4 (H3K4me2 and H3K4me3), H3K27 (H3K27me2 and H3K27me3), and H3K36 (H3K36me2 and H3K36me3) in sorghum and analyze their functional relationships. Unlike other histone methylation, the codecoration of H3K4me3 and H3K36me3 is negatively associated with the production of lincRNAs in the context of active expression of mRNA genes. Our data demonstrated that H3K4me3 may act as a complementary component to H3K36me3 in the transcriptional regulatory process. Moreover, we observe that both H3K4me3 and H3K36me3 are involved in the negative-going regulation of plant lincRNA and mRNA genes. Our data provide a genome-wide landscape of histone methylation in sorghum, decrypt its reciprocity, and shed light on its transcriptional regulation roles in mRNA and lncRNA genes.

An optimised chromatin immunoprecipitation (ChIP) method for starchy leaves of Nicotiana benthamiana to study histone modifications of an allotetraploid plant

All flowering plants have evolved through multiple rounds of polyploidy throughout the evolutionary process. Intergenomic interactions between subgenomes in polyploid plants are predicted to induce chromatin modifications such as histone modifications to regulate expression of gene homoeologs. Nicotiana benthamiana is an ancient allotetraploid plant with ecotypes collected from climatically diverse regions of Australia. Studying the chromatin landscape of this unique collection will likely shed light on the importance of chromatin modifications in gene regulation in polyploids as well its implications in adaptation of plants in environmentally diverse conditions. Generally, chromatin immunoprecipitation and high throughput DNA sequencing (ChIP-seq) is used to study chromatin modifications. However, due to the starchy nature of mature N. benthamiana leaves, previously published protocols were unsuitable. The higher amounts of starch in leaves that co-precipitated with nuclei hindered downstream processing of DNA. Here we present an optimised ChIP protocol for N. benthamiana leaves to facilitate comparison of chromatin modifications in two closely related ecotypes. Several steps of ChIP were optimised including tissue harvesting, nuclei isolation, nuclei storage, DNA shearing and DNA recovery. Commonly available antibodies targeting histone 3 lysine 4 trimethylation (H3K4me3) and histone 3 lysine 9 dimethylation (H3K9me2) histone modifications were used and success of ChIP was confirmed by PCR and next generation sequencing. Collectively, our optimised method is the first comprehensive ChIP method for mature starchy leaves of N. benthamiana to enable studies of chromatin landscape at the genome-wide scale.

Chromatin regulation in plant hormone and plant stress responses

The gene expression is tightly regulated temporally and spatially to ensure the plant and animal proper development, function, growth, and survival under different environmental conditions. Chromatin regulation plays a central role in the gene expression by providing transcription factors and the transcription machinery with dynamic access to an otherwise tightly packaged genome. In this review, we will summarize recent progress in understanding the roles of chromatin regulation in the gene expression, and their contribution to the plant hormone and stress responses. We highlight the most recent publications within this topic and underline the roles of chromatin regulation in gene expression.

Whole-genome landscape of H3K4me3, H3K36me3 and H3K9ac and their association with gene expression during Paulownia witches' broom disease infection and recovery processes

Histone methylation and acetylation participate in the modulation of gene expression. Here, chromatin immunoprecipitation sequencing (ChIP-Seq) was used to determine genome-wide patterns of three histone modifications, H3K4me3, H3K36me3, and H3K9ac (associated with actively expressed genes) and their associations with gene expression in Paulownia fortunei following phytoplasma infection and recovery from Paulownia witches' broom (PaWB) disease after methyl methane sulfonate treatment. The three histone marks were preferentially deposited in genic regions, especially downstream of transcription start sites, and were highly concurrent with gene expression. Genes with all three histone marks exhibited the highest expression levels. Based on the comparison scheme, we detected 365, 2244, and 752 PaWB-associated genes with H3K4me3, H3K36me3, and H3K9ac marks, separately. KEGG pathway analysis showed that these genes were involved in plant-pathogen interaction, plant hormone signal transduction, and starch and sucrose metabolism. A small proportion of differentially modified genes showed changes in expression in response to phytoplasma infection, including genes involved in calcium ion signal transduction, abscisic acid signal transduction, and ethylene biosynthesis. This comprehensive analysis of genome-wide histone modifications and gene expression in Paulownia following phytoplasma infection provides new insights into the epigenetic responses to phytoplasma infection and will be useful for further studies on epigenetic regulation mechanisms in plants under biotic stress.

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

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

The transcription factor GATA10 regulates fertility conversion of a two-line hybrid tms5 mutant rice via the modulation of UbL40 expression

The thermosensitive genic male sterile 5 (tms5) mutation causes thermosensitive genic male sterility in rice (Oryza sativa) through loss of RNase Z^S1 function, which influences ubiquitin fusion ribosomal protein L40 (Ub_L40 ) messenger RNA levels during male development. Here, we used ATAC-seq, combined with analysis of H3K9ac and H3K4me2, to identify changes in accessible chromatin during fertility conversion of the two-line hybrid rice Wuxiang S (WXS) derived from a mutant tms5 allele. Furthermore, RNA-seq and bioinformatic analyses identified specific transcription factors (TFs) in differentially accessible chromatin regions. Among these TFs, only GATA10 targeted Ub_L40 . Osgata10 knockout mutations, which resulted in low expression of Ub_L40 and a tendency toward male fertility, confirmed that GATA10 regulated fertility conversion via the modulation of Ub_L40 . Meanwhile, GATA10 acted as a mediator for interactions with ERF65, which revealed that transcriptional regulation is a complex process involving multiple complexes of TFs, namely TF modules. It appears that the ERF141/MADS7/MADS50/MYB modules affect metabolic processes that control anther and pollen development, especially cell wall formation. Our analysis revealed that these modules directly or indirectly affect metabolic pathway-related genes to coordinate plant growth with proper anther development, and furthermore, that GATA10 regulates fertility conversion via the modulation of Ub_L40 expression.

Integrative analysis of reference epigenomes in 20 rice varieties

Epigenomic modifications are instrumental for transcriptional regulation, but comprehensive reference epigenomes remain unexplored in rice. Here, we develop an enhanced chromatin immunoprecipitation (eChIP) approach for plants, and generate genome-wide profiling of five histone modifications and RNA polymerase II occupancy with it. By integrating chromatin accessibility, DNA methylation, and transcriptome datasets, we construct comprehensive epigenome landscapes across various tissues in 20 representative rice varieties. Approximately 81.8% of rice genomes are annotated with different epigenomic properties. Refinement of promoter regions using open chromatin and H3K4me3-marked regions provides insight into transcriptional regulation. We identify extensive enhancer-like promoters with potential enhancer function on transcriptional regulation through chromatin interactions. Active and repressive histone modifications and the predicted enhancers vary largely across tissues, whereas inactive chromatin states are relatively stable. Together, these datasets constitute a valuable resource for functional element annotation in rice and indicate the central role of epigenomic information in understanding transcriptional regulation.

Understanding epigenomics based on the rice model

The purpose of this paper provides a comprehensive overview of the recent researches on rice epigenomics, including DNA methylation, histone modifications, noncoding RNAs, and three-dimensional genomics. The challenges and perspectives for future research in rice are discussed. Rice as a model plant for epigenomic studies has much progressed current understanding of epigenetics in plants. Recent results on rice epigenome profiling and three-dimensional chromatin structure studies reveal specific features and implication in gene regulation during rice plant development and adaptation to environmental changes. Results on rice chromatin regulator functions shed light on mechanisms of establishment, recognition, and resetting of epigenomic information in plants. Cloning of several rice epialleles associated with important agronomic traits highlights importance of epigenomic variation in rice plant growth, fitness, and yield. In this review, we summarize and analyze recent advances in rice epigenomics and discuss challenges and directions for future research in the field.

Drought-inducible changes in the histone modification H3K9ac are associated with drought-responsive gene expression in Brachypodium distachyon

H3K9ac, an epigenetic marker, is widely distributed in plant genomes. H3K9ac enhances gene expression, which is highly conserved in eukaryotes. However, genome-wide studies of H3K9ac in monocot species are limited, and the changes in H3K9ac under drought stress for individual genes are still not clear. We analysed changes in the H3K9ac level of Brachypodium distachyon under 20% PEG-6000-simulated drought stress conditions. We also performed chromatin immunoprecipitation, followed by next generation sequencing (ChIP-seq) on H3K9ac to reveal changes in H3K9ac for individual genes at the genome-wide level. Our study showed that H3K9ac was mainly enriched in gene exon regions. Drought increased or decreased the H3K9ac level at specific genomic loci. We identified 40 genes associated with increased H3K9ac levels and 36 genes associated with decreased H3K9ac levels under drought stress. Further, RT-qPCR analyses showed that H3K9ac was positively associated with gene expression of those drought-responsive genes. We conclude that H3K9ac enhances the expression level of a large number of drought-responsive genes under drought stress in B. distachyon. The data presented here will help to reveal the correlation of some specific drought-responsive genes and their enriched H3K9ac levels in the model plant B. distachyon.

Homeostasis of histone acetylation is critical for auxin signaling and root morphogenesis

The auxin signaling and root morphogenesis are harmoniously controlled by two counteracted teams including (1) auxin/indole-3-acetic acid (AUX/IAA)-histone deacetylase (HDA) and (2) auxin response factor (ARF)-histone acetyltransferase (HAT). The involvement of histone acetylation in the regulation of transcription was firstly reported a few decades ago. In planta, auxin is the first hormone group that was discovered and it is also the most studied phytohormone. Current studies have elucidated the functions of histone acetylation in the modulation of auxin signaling as well as in the regulation of root morphogenesis under both normal and stress conditions. Based on the recent outcomes, this review is to provide a hierarchical view about the functions of histone acetylation in auxin signaling and root morphogenesis. In this report, we suggest that the auxin signaling must be controlled harmoniously by two counteracted teams including (1) auxin/indole-3-acetic acid (AUX/IAA)-histone deacetylase (HDA) and (2) auxin response factor (ARF)-histone acetyltransferase (HAT). Moreover, the balance in auxin signaling is very critical to contribute to normal root morphogenesis.

UNBRANCHED3 Expression and Inflorescence Development is Mediated by UNBRANCHED2 and the Distal Enhancer, KRN4, in Maize

Enhancers are cis-acting DNA segments with the ability to increase target gene expression. They show high sensitivity to DNase and contain specific DNA elements in an open chromatin state that allows the binding of transcription factors (TFs). While numerous enhancers are annotated in the maize genome, few have been characterized genetically. KERNEL ROW NUMBER4 (KRN4), an intergenic quantitative trait locus for kernel row number, is assumed to be a cis-regulatory element of UNBRANCHED3 (UB3), a key inflorescence gene. However, the mechanism by which KRN4 controls UB3 expression remains unclear. Here, we found that KRN4 exhibits an open chromatin state, harboring sequences that showed high enhancer activity toward the 35S and UB3 promoters. KRN4 is bound by UB2-centered transcription complexes and interacts with the UB3 promoter by three duplex interactions to affect UB3 expression. Sequence variation at KRN4 enhances ub2 and ub3 mutant ear fasciation. Therefore, we suggest that KRN4 functions as a distal enhancer of the UB3 promoter via chromatin interactions and recruitment of UB2-centered transcription complexes for the fine-tuning of UB3 expression in meristems of ear inflorescences. These results provide evidence that an intergenic region helps to finely tune gene expression, providing a new perspective on the genetic control of quantitative traits.

Low-Phosphate Chromatin Dynamics Predict a Cell Wall Remodeling Network in Rice Shoots

Phosphorus (P) is an essential plant macronutrient vital to fundamental metabolic processes. Plant-available P is low in most soils, making it a frequent limiter of growth. Declining P reserves for fertilizer production exacerbates this agricultural challenge. Plants modulate complex responses to fluctuating P levels via global transcriptional regulatory networks. Although chromatin structure plays a substantial role in controlling gene expression, the chromatin dynamics involved in regulating P homeostasis have not been determined. Here we define distinct chromatin states across the rice (Oryza sativa) genome by integrating multiple chromatin marks, including the H2A.Z histone variant, H3K4me3 modification, and nucleosome positioning. In response to P starvation, 40% of all protein-coding genes exhibit a transition from one chromatin state to another at their transcription start site. Several of these transitions are enriched in subsets of genes differentially expressed under P deficiency. The most prominent subset supports the presence of a coordinated signaling network that targets cell wall structure and is regulated in part via a decrease of H3K4me3 at transcription start sites. The P starvation-induced chromatin dynamics and correlated genes identified here will aid in enhancing P use efficiency in crop plants, benefitting global agriculture.

An inferred fitness consequence map of the rice genome

The extent to which sequence variation impacts plant fitness is poorly understood. High-resolution maps detailing the constraint acting on the genome, especially in regulatory sites, would be beneficial as functional annotation of noncoding sequences remains sparse. Here, we present a fitness consequence (fitCons) map for rice (Oryza sativa). We inferred fitCons scores (ρ) for 246 inferred genome classes derived from nine functional genomic and epigenomic datasets, including chromatin accessibility, messenger RNA/small RNA transcription, DNA methylation, histone modifications and engaged RNA polymerase activity. These were integrated with genome-wide polymorphism and divergence data from 1,477 rice accessions and 11 reference genome sequences in the Oryzeae. We found ρ to be multimodal, with ~9% of the rice genome falling into classes where more than half of the bases would probably have a fitness consequence if mutated. Around 2% of the rice genome showed evidence of weak negative selection, frequently at candidate regulatory sites, including a novel set of 1,000 potentially active enhancer elements. This fitCons map provides perspective on the evolutionary forces associated with genome diversity, aids in genome annotation and can guide crop breeding programs.

Epigenetic signatures of stress adaptation and flowering regulation in response to extended drought and recovery in Zea mays

During their lifespan, plants respond to a multitude of stressful factors. Dynamic changes in chromatin and concomitant transcriptional variations control stress response and adaptation, with epigenetic memory mechanisms integrating environmental conditions and appropriate developmental programs over the time. Here we analyzed transcriptome and genome-wide histone modifications of maize plants subjected to a mild and prolonged drought stress just before the flowering transition. Stress was followed by a complete recovery period to evaluate drought memory mechanisms. Three categories of stress-memory genes were identified: i) "transcriptional memory" genes, with stable transcriptional changes persisting after the recovery; ii) "epigenetic memory candidate" genes in which stress-induced chromatin changes persist longer than the stimulus, in absence of transcriptional changes; iii) "delayed memory" genes, not immediately affected by the stress, but perceiving and storing stress signal for a delayed response. This last memory mechanism is described for the first time in drought response. In addition, applied drought stress altered floral patterning, possibly by affecting expression and chromatin of flowering regulatory genes. Altogether, we provided a genome-wide map of the coordination between genes and chromatin marks utilized by plants to adapt to a stressful environment, describing how this serves as a backbone for setting stress memory.

Genome-wide Profiling of Histone Lysine Butyrylation Reveals its Role in the Positive Regulation of Gene Transcription in Rice

BACKGROUND

Histone modifications play important roles in growth and development of rice (Oryza sativa L.). Lysine butyrylation (Kbu) with a four-carbon chain is a newly-discovered histone acylation modification in rice.

MAIN BODY

In this study, we performed chromatin immunoprecipitation sequencing (ChIP-seq) analyses, the result showed that major enrichment of histone Kbu located in genebody regions of rice genome, especially in exons. The enrichment level of Kbu histone modification is positively correlated with gene expression. Furthermore, we compared Kbu with DNase-seq and other histone modifications in rice. We found that 60.06% Kub enriched region co-located with DHSs in intergenic regions. The similar profiles were detected among Kbu and several acetylation modifications such as H3K4ac, H3K9ac, and H3K23ac, indicating that Kbu modification is an active signal of transcription. Genes with both histone Kbu and one other acetylation also had significantly increased expression compared with genes with only one acetylation. Gene Ontology (GO) enrichment analysis revealed that these genes with histone Kbu can regulate multiple metabolic process in different rice varieties.

CONCLUSION

Our study showed that the lysine butyrylation modificaiton may promote gene expression as histone acetylation and will provide resources for futher studies on histone Kbu and other epigenetic modifications in plants.

Identification and functional evaluation of accessible chromatin associated with wood formation in Eucalyptus grandis

Accessible chromatin changes dynamically during development and harbours functional regulatory regions which are poorly understood in the context of wood development. We explored the importance of accessible chromatin in Eucalyptus grandis in immature xylem generally, and MYB transcription factor-mediated transcriptional programmes specifically. We identified biologically reproducible DNase I Hypersensitive Sites (DHSs) and assessed their functional significance in immature xylem through their associations with gene expression, epigenomic data and DNA sequence conservation. We identified in vitro DNA binding sites for six secondary cell wall-associated Eucalyptus MYB (EgrMYB) transcription factors using DAP-seq, reconstructed protein-DNA networks of predicted targets based on binding sites within or outside DHSs and assessed biological enrichment of these networks with published datasets. 25 319 identified immature xylem DHSs were associated with increased transcription and significantly enriched for various epigenetic signatures (H3K4me3, H3K27me3, RNA pol II), conserved noncoding sequences and depleted single nucleotide variants. Predicted networks built from EgrMYB binding sites located in accessible chromatin were significantly enriched for systems biology datasets relevant to wood formation, whereas those occurring in inaccessible chromatin were not. Our study demonstrates that DHSs in E. grandis immature xylem, most of which are intergenic, are of functional significance to gene regulation in this tissue.

Genome-wide identification of histone methylation (H3K9me2) and acetylation (H4K12ac) marks in two ecotypes of switchgrass (Panicum virgatum L.)

BACKGROUND

Histone modifications play a significant role in the regulation of transcription and various biological processes, such as development and regeneration. Though a few genomic (including DNA methylation patterns) and transcriptomic studies are currently available in switchgrass, the genome-wide distribution of histone modifications has not yet been studied to help elucidate gene regulation and its application to switchgrass improvement.

RESULTS

This study provides a comprehensive epigenomic analyses of two contrasting switchgrass ecotypes, lowland (AP13) and upland (VS16), by employing chromatin immunoprecipitation sequencing (ChIP-Seq) with two histone marks (suppressive- H3K9me2 and active- H4K12ac). In this study, most of the histone binding was in non-genic regions, and the highest enrichment was seen between 0 and 2 kb regions from the transcriptional start site (TSS). Considering the economic importance and potential of switchgrass as a bioenergy crop, we focused on genes, transcription factors (TFs), and pathways that were associated with C4-photosynthesis, biomass, biofuel production, biotic stresses, and abiotic stresses. Using quantitative real-time PCR (qPCR) the relative expression of five genes selected from the phenylpropanoid-monolignol pathway showed preferential binding of acetylation marks in AP13 rather than in VS16.

CONCLUSIONS

The genome-wide histone modifications reported here can be utilized in understanding the regulation of genes important in the phenylpropanoid-monolignol biosynthesis pathway, which in turn, may help understand the recalcitrance associated with conversion of biomass to biofuel, a major roadblock in utilizing lignocellulosic feedstocks.

PCSD: a plant chromatin state database

Genome-wide maps of chromatin states have become a powerful representation of genome annotation and regulatory activity. We collected public and in-house plant epigenomic data sets and applied a Hidden Markov Model to define chromatin states, which included 290 553 (36 chromatin states), 831 235 (38 chromatin states) and 3 936 844 (26 chromatin states) segments across the whole genome of Arabidopsis thaliana, Oryza sativa and Zea mays, respectively. We constructed a Plant Chromatin State Database (PCSD, http://systemsbiology.cau.edu.cn/chromstates) to integrate detailed information about chromatin states, including the features and distribution of states, segments in states and related genes with segments. The self-organization mapping (SOM) results for these different chromatin signatures and UCSC Genome Browser for visualization were also integrated into the PCSD database. We further provided differential SOM maps between two epigenetic marks for chromatin state comparison and custom tools for new data analysis. The segments and related genes in SOM maps can be searched and used for motif and GO analysis, respectively. In addition, multi-species integration can be used to discover conserved features at the epigenomic level. In summary, our PCSD database integrated the identified chromatin states with epigenetic features and may be beneficial for communities to discover causal functions hidden in plant chromatin.

Transcriptional and epigenetic adaptation of maize chromosomes in Oat-Maize addition lines

By putting heterologous genomic regulatory systems into contact, chromosome addition lines derived from interspecific or intergeneric crosses allow the investigation of transcriptional regulation in new genomic environments. Here, we report the transcriptional and epigenetic adaptation of stably inherited alien maize chromosomes in two oat–maize addition (OMA) lines. We found that the majority of maize genes displayed maize-specific transcription in the oat genomic environment. Nevertheless, a quarter of the expressed genes encoded by the two maize chromosomes were differentially expressed genes (DEGs). Notably, highly conserved orthologs were more severely differentially expressed in OMAs than less conserved orthologs. Additionally, syntenic genes and highly abundant genes were over-represented among DEGs. Gene suppression was more common than activation among the DEGs; however, the genes in the former maize pericentromere, which expanded to become the new centromere in OMAs, were activated. Histone modifications (H3K4me3, H3K9ac and H3K27me3) were consistent with these transcriptome results. We expect that cis regulation is responsible for unchanged expression in OMA versus maize; and trans regulation is the predominant mechanism behind DEGs. The genome interaction identified here reveals the important consequences of interspecific/intergeneric crosses and potential mechanisms of plant evolution when genomic environments interact.

The unique epigenetic features of Pack-MULEs and their impact on chromosomal base composition and expression spectrum

Acquisition and rearrangement of host genes by transposable elements (TEs) is an important mechanism to increase gene diversity as exemplified by the ∼3000 Pack-Mutator-like TEs in the rice genome which have acquired gene sequences (Pack-MULEs), yet remain enigmatic. To identify signatures of functioning Pack-MULEs and Pack-MULE evolution, we generated transcriptome, translatome, and epigenome datasets and compared Pack-MULEs to genes and other TE families. Approximately 40% of Pack-MULEs were transcribed with 9% having translation evidence, clearly distinguishing them from other TEs. Pack-MULEs exhibited a unique expression profile associated with specificity in reproductive tissues that may be associated with seed traits. Expressed Pack-MULEs resemble regular protein-coding genes as exhibited by a low level of DNA methylation, association with active histone marks and DNase I hypersensitive sites, and an absence of repressive histone marks, suggesting that a substantial fraction of Pack-MULEs are potentially functional in vivo. Interestingly, the expression capacity of Pack-MULEs is independent of the local genomic environment, and the insertion and expression of Pack-MULEs may have altered the local chromosomal expression pattern as well as counteracted the impact of recombination on chromosomal base composition, which has profound consequences on the evolution of chromosome structure.

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.

Brachypodium histone deacetylase BdHD1 positively regulates ABA and drought stress responses

Despite recent evidence that HDACs are involved in the environmental stress responses of plants, their roles in the abiotic stress responses of monocot plants remain largely unexplored. We investigated a HDAC gene, Bradi3g08060 (BdHD1), in Brachypodium distachyon. The Brachypodium BdHD1-overexpression plants displayed a hypersensitive phenotype to ABA and exhibited better survival under drought conditions. On the other hand, the RNA-interference plants were insensitive to ABA and showed low survival under drought stress. At the genome-wide level, overexpression of BdHD1 led to lower H3K9 acetylation at the transcriptional start sites of 230 genes than in the wild type plants under the drought treatment. We validated our ChIP-Seq data on 10 selected transcription factor genes from the 230 drought-specific genes. These genes exhibited much lower expression in the BdHD1-overexpression plants compared to the wild type plants under drought stress. We further identified an ABA-inducible transcription factor gene BdWRKY24 that was repressed in BdHD1-OE plants, but highly expressed in RNA-interference plants under drought stress. These results indicate that BdHD1 plays a positive role in ABA sensitivity and drought stress tolerance and they provide a link between the role of BdHD1 and the drought stress response at a genome-wide level in Brachypodium.

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.

Dynamic control of enhancer activity drives stage-specific gene expression during flower morphogenesis

Enhancers are critical for developmental stage-specific gene expression, but their dynamic regulation in plants remains poorly understood. Here we compare genome-wide localization of H3K27ac, chromatin accessibility and transcriptomic changes during flower development in Arabidopsis. H3K27ac prevalently marks promoter-proximal regions, suggesting that H3K27ac is not a hallmark for enhancers in Arabidopsis. We provide computational and experimental evidence to confirm that distal DNase І hypersensitive sites are predictive of enhancers. The predicted enhancers are highly stage-specific across flower development, significantly associated with SNPs for flowering-related phenotypes, and conserved across crucifer species. Through the integration of genome-wide transcription factor (TF) binding datasets, we find that floral master regulators and stage-specific TFs are largely enriched at developmentally dynamic enhancers. Finally, we show that enhancer clusters and intronic enhancers significantly associate with stage-specific gene regulation by floral master TFs. Our study provides insights into the functional flexibility of enhancers during plant development, as well as hints to annotate plant enhancers.

Enhancer elements can control spatial and temporal patterns of gene expression. Here the authors profile chromatin accessibility, histone modifications and gene expression during Arabidopsis flower development providing evidence for sets of distal enhancers acting in a highly stage-specific manner.