The Role of Chromatin Modifications in the Evolution of Giant Plant Genomes.

Angiosperm genome sizes (GS) range ~2400-fold and comprise genes and their regulatory regions, repeats, semi-degraded repeats, and 'dark matter'. The latter represents repeats so degraded that they can no longer be recognised as repetitive. In exploring whether the histone modifications associated with chromatin packaging of these contrasting genomic components are conserved across the diversity of GS in angiosperms, we compared immunocytochemistry data for two species whose GS differ ~286-fold. We compared published data for Arabidopsis thaliana with a small genome (GS = 157 Mbp/1C) with newly generated data from Fritillaria imperialis, which has a giant genome (GS = 45,000 Mbp/1C). We compared the distributions of the following histone marks: H3K4me1, H3K4me2, H3K9me1, H3K9me2, H3K9me3, H3K27me1, H3K27me2, and H3K27me3. Assuming these histone marks are associated with the same genomic features across all species, irrespective of GS, our comparative analysis enables us to suggest that while H3K4me1 and H3K4me2 methylation identifies genic DNA, H3K9me3 and H3K27me3 marks are associated with 'dark matter', H3K9me1 and H3K27me1 mark highly homogeneous repeats, and H3K9me2 and H3K27me2 mark semi-degraded repeats. The results have implications for our understanding of epigenetic profiles, chromatin packaging and the divergence of genomes, and highlight contrasting organizations of the chromatin within the nucleus depending on GS itself

Arabidopsis SWR1-associated protein methyl-CpG-binding domain 9 is required for histone H2A.Z deposition.

Deposition of the histone variant H2A.Z by the SWI2/SNF2-Related 1 chromatin remodeling complex (SWR1-C) is important for gene regulation in eukaryotes, but the composition of the Arabidopsis SWR1-C has not been thoroughly characterized. Here, we aim to identify interacting partners of a conserved Arabidopsis SWR1 subunit ACTIN-RELATED PROTEIN 6 (ARP6). We isolate nine predicted components and identify additional interactors implicated in histone acetylation and chromatin biology. One of the interacting partners, methyl-CpG-binding domain 9 (MBD9), also strongly interacts with the Imitation SWItch (ISWI) chromatin remodeling complex. MBD9 is required for deposition of H2A.Z at a distinct subset of ARP6-dependent loci. MBD9 is preferentially bound to nucleosome-depleted regions at the 5' ends of genes containing high levels of activating histone marks. These data suggest that MBD9 is a SWR1-C interacting protein required for H2A.Z deposition at a subset of actively transcribing genes

Identification of the expressome by machine learning on omics data.

Accurate annotation of plant genomes remains complex due to the presence of many pseudogenes arising from whole-genome duplication-generated redundancy or the capture and movement of gene fragments by transposable elements. Machine learning on genome-wide epigenetic marks, informed by transcriptomic and proteomic training data, could be used to improve annotations through classification of all putative protein-coding genes as either constitutively silent or able to be expressed. Expressed genes were subclassified as able to express both mRNAs and proteins or only RNAs, and CG gene body methylation was associated only with the former subclass. More than 60,000 protein-coding genes have been annotated in the reference genome of maize inbred B73. About two-thirds of these genes are transcribed and are designated the filtered gene set (FGS). Classification of genes by our trained random forest algorithm was accurate and relied only on histone modifications or DNA methylation patterns within the gene body; promoter methylation was unimportant. Other inbred lines are known to transcribe significantly different sets of genes, indicating that the FGS is specific to B73. We accurately classified the sets of transcribed genes in additional inbred lines, arising from inbred-specific DNA methylation patterns. This approach highlights the potential of using chromatin information to improve annotations of functional genes

DNA methylation patterns of Brachypodium distachyon chromosomes and their alteration by 5-azacytidine treatment

Sequential immunolocalisation of 5-methylcytosine (5-MeC) and fluorescence in situ hybridisation with chromosome-specific BAC clones were performed on Brachypodium distachyon mitotic metaphase chromosomes to determine specific DNA methylation patterns of each chromosome in the complement. In the majority of cells examined, chromosomes Bd4 and Bd5, which bear the loci of 5S and 35S ribosomal DNA, respectively, had characteristic 5-MeC patterns. In contrast, the distribution of 5-MeC along the metacentric chromosome pairs Bd1, Bd2 and Bd3 was more variable. There were numerous differences in distribution of methylated sites between homologous chromosomes as well as between chromosome arms. Some chromosome sites, such as pericentromeric regions, were highly methylated in all chromosomes. Additionally, the influence of a hypomethylating agent, 5-azacytidine, on B. distachyon chromosome methylation patterns was confirmed. It was found that some chromosome pairs underwent demethylation more easily than others, but there was no apparent regularity in demethylation of particular chromosome segments

ChIP-Hub provides an integrative platform for exploring plant regulome

Plant genomes encode a complex and evolutionary diverse regulatory grammar that forms the basis for most life on earth. A wealth of regulome and epigenome data have been generated in various plant species, but no common, standardized resource is available so far for biologists. Here, we present ChIP-Hub, an integrative web-based platform in the ENCODE standards that bundles >10,000 publicly available datasets reanalyzed from >40 plant species, allowing visualization and meta-analysis. We manually curate the datasets through assessing ~540 original publications and comprehensively evaluate their data quality. As a proof of concept, we extensively survey the co-association of different regulators and construct a hierarchical regulatory network under a broad developmental context. Furthermore, we show how our annotation allows to investigate the dynamic activity of tissue-specific regulatory elements (promoters and enhancers) and their underlying sequence grammar. Finally, we analyze the function and conservation of tissue-specific promoters, enhancers and chromatin states using comparative genomics approaches. Taken together, the ChIP-Hub platform and the analysis results provide rich resources for deep exploration of plant ENCODE. ChIP-Hub is available at https://biobigdata.nju.edu.cn/ChIPHub/.Peer Reviewe

Locus-Specific Ribosomal RNA Gene Silencing in Nucleolar Dominance

The silencing of one parental set of rRNA genes in a genetic hybrid is an epigenetic phenomenon known as nucleolar dominance. We showed previously that silencing is restricted to the nucleolus organizer regions (NORs), the loci where rRNA genes are tandemly arrayed, and does not spread to or from neighboring protein-coding genes. One hypothesis is that nucleolar dominance is the net result of hundreds of silencing events acting one rRNA gene at a time. A prediction of this hypothesis is that rRNA gene silencing should occur independent of chromosomal location. An alternative hypothesis is that the regulatory unit in nucleolar dominance is the NOR, rather than each individual rRNA gene, in which case NOR localization may be essential for rRNA gene silencing. To test these alternative hypotheses, we examined the fates of rRNA transgenes integrated at ectopic locations. The transgenes were accurately transcribed in all independent transgenic Arabidopsis thaliana lines tested, indicating that NOR localization is not required for rRNA gene expression. Upon crossing the transgenic A. thaliana lines as ovule parents with A. lyrata to form F1 hybrids, a new system for the study of nucleolar dominance, the endogenous rRNA genes located within the A. thaliana NORs are silenced. However, rRNA transgenes escaped silencing in multiple independent hybrids. Collectively, our data suggest that rRNA gene activation can occur in a gene-autonomous fashion, independent of chromosomal location, whereas rRNA gene silencing in nucleolar dominance is locus-dependent

Identification and computational analysis of differential H3K27me3 targets between <em>Arabidopsis thaliana</em> accessions

Histone H3 lysine 27 trimethylation (H3K27me3) and lysine 9 dimethylation (H3K9me2) are two independent repressive chromatin modifications in Arabidopsis thaliana. H3K27me3 is established and maintained by Polycomb repressive complexes whereas H3K9me2 is catalyzed by histone methyltransferases SUVH(4-6). H3K27me3 mostly targets at protein coding genes in euchromatin which are reversible in repression. H3K9me2 mainly targets at transposons and repetitive sequences which should be constitutively silenced. Both marks can spread to flanking regions after initialization and they have been shown to be mutually exclusive in distribution in the Arabidopsis genome. In this study, the extent of natural variation of H3K27me3 in the two accessions of Arabidopsis thaliana, Landsberg erecta (Ler) and Columbia (Col), and their hybrids was analyzed using chromatin immunoprecipitation followed by microarray or sequencing analysis (ChIP-chip and ChIP-seq). A computational workflow was implemented that includes remapping of probes to the Col and Ler genome assemblies in order to exclude differential signals due to genome polymorphisms. The majority of genes that are H3K27me3 targets in Col are also targets in Ler and the F1 of reciprocal crosses. A small number of Ler-specific H3K27me3 targets were detected and well validated with independent ChIP-PCR whereas the Col-specific targets have not been confirmed so far. Ler-specific H3K27me3 targets showed an allele-specific H3K27me3 in both hybrids, consistent with a cis-regulatory mechanism for establishing H3K27me3. Five Ler-specific H3K27me3 targets were marked by H3K4me3 in Col. Consistent with the activation role of H3K4me3 during transcription, the differential H3K27me3 of the these five genes accords to the expression variation between the two accessions. For the majority of Ler-specific H3K27me3 targets, no expression could be detected in Col, Ler or 17 other Arabidopsis accessions. Instead of H3K27me3, the antagonistic mark H3K9me2 and other heterochromatic features were observed at these loci in Col. More frequently than expected, transposable elements were found neighboring these loci in Col, and in many cases these transposable elements are missing in the Ler genome assembly. We propose a model where a transposon insertion specific to Col results in recruitment of H3K9me2, which spreads to neighboring genes already in a repressed state through H3K27me3, resulting in Ler-specific H3K27me3 as the ancestral state