<title>Introduction</title> <p>Post-translational modification of core histones considerably extends the information potential of the genetic code ##REF##10638745##[1]##,##REF##10944586##[2]##. Methylation of specific residues within the amino-terminal tail of nucleosomal histone H3, in particular, has been tied to activation or repression of transcription within the associated gene(s). For example, where studied in budding yeast and human, nucleosomes containing H3 trimethylated at lysine-4 (H3K4me3) are globally enriched near the transcriptional start sites (TSS) and 5′ regions of active genes, with the degree of enrichment correlating with gene activity ##REF##12353038##[3]##–##REF##17030614##[6]##. In yeast, this pattern is thought to be an outcome of cotranscriptional recruitment of the histone methyltransferase SET1 during the early elongation phase ##REF##12667454##[7]##,##REF##12667453##[8]##. SET1 and homologous methyltransferases such as Trithorax (Trx) in fruit fly and mixed lineage leukemia 1 (MLL1) in human target nucleosomal H3K4 for methylation as components of larger protein complexes ##REF##11742990##[9]##–##REF##14730313##[11]##. However, methylated H3K4 likely serves an instructive and promotive role in transcription as well: methylated H3K4 is required for efficient chromatin remodeling at promoters ##REF##14636589##[12]##,##REF##16728976##[13]##, and potentially enhances interaction with the SET1-related complexes ##REF##15960974##[14]##. Thus, H3K4me3 may define a mechanism that reinforces the active state of transcription. Di- and trimethylated H3K36 (H3K36me2/me3) is prevalent within transcribed regions in yeast and human, especially near the 3′ ends ##REF##12773564##[15]##,##REF##15760899##[16]##, reflecting cotranscriptional activity of the H3K36-specific SET2 methyltransferase during elongation ##REF##12381723##[17]##. Although localized within active genes, H3K36 methylation probably has an overall negative influence on transcription that is mediated at least in part through recruitment of histone deacetylase activity and consequent maintenance of low acetylation levels ##REF##16286007##[18]##–##REF##16286008##[20]##. Repressing histone acetylation within transcribed regions is expected to promote internucleosomal interactions and/or chromatin assembly in the wake of PolII, thus minimizing inappropriate intragenic transcriptional initiation at cryptic sites.</p> <p>Methylation of H3K27 is an elaboration seen only in higher eukaryotes, where it has been linked with transcriptional repression. Dimethylated H3K27 (H3K27me2) is abundant within heterochromatin, whereas in human and fruit fly, trimethylated H3K27 (H3K27me3) is found in frequent islands scattered throughout euchromatin, with extended domains surrounding Polycomb-Group (PcG) protein binding sites including the Hox loci ##REF##12351676##[21]##–##REF##16628213##[25]##. In plants, H3K27me3 marks weakly expressed and/or developmentally silenced genes, including known targets of plant PcG proteins ##REF##16878125##[26]##–##REF##17439305##[29]##. H3K27 methylation may repress transcription through several mechanisms, including recruitment of PRC1 in metazoans ##REF##12351676##[21]##, and in plants, direct binding to LHP1, the homolog of Heterochromatin Protein 1 ##REF##17542647##[28]##. The conserved PcG protein Enhancer of zeste [E(z)] and associated proteins, designated Polycomb Repressive Complex 2 (PRC2) mediate methylation of H3K27, thus connecting this modification to the maintenance of gene silencing ##REF##17173055##[30]##.</p> <p>The PolII-associating factor 1 complex (Paf1C), minimally composed of Paf1, Ctr9, Cdc73, Rtf1, and Leo1 has an important role in establishing patterns of methylated H3K4 and H3K36 by promoting ubiquitination of histone H2B ##REF##12876293##[31]##,##REF##12876294##[32]## and linking elongating PolII with SET1 and SET2 ##REF##12667454##[7]##,##REF##12667453##[8]##,##REF##12773564##[15]##. Paf1C also has transcription-related roles potentially independent of its function in histone modification, related to elongation ##REF##11927560##[33]##, suppression of intragenic initiation ##REF##17948059##[34]##, poly(A) site selection ##REF##16246724##[35]##, mRNA polyadenylation ##REF##15149594##[36]##, and 3′ end formation on nonpolyadenylated PolII-generated transcripts ##REF##16246725##[37]##.</p> <p>Components of Paf1C are also conserved in higher eukaryotes. The product of the human HRPT2 gene, parafibromin, shows moderate homology with Cdc73, and interacts with human homologs of Paf1, Ctr9, and Leo1 as well as elongating (Ser-2/Ser-5 phosphorylated) PolII <italic>in vivo</italic>
##REF##15632063##[38]##–##REF##16024656##[40]##. The human Paf1C complex (hPAF) also contains hSki8, a protein that physically associates with the exosome, required for 3′–5′ mRNA degradation ##REF##16024656##[40]##. Similar to yeast Paf1C, hPAF was localized to transcriptionally active genes, and disruption of hPAF led to global reduction in H3K4me3 levels ##REF##16024656##[40]##. Both parafibromin and the human Paf1 homolog are known to be disrupted associated with cancers, although potential mechanisms have not been well described ##REF##12434154##[41]##,##REF##16491129##[42]##. In fruit fly, homologs of Paf1, Rtf1 (dRtf1) and Cdc73 (<italic>hyrax</italic>) colocalize with transcribing PolII ##REF##16354696##[43]##,##REF##16630820##[44]##, and at least dRtf1 is required to maintain global H3K4me3 in chromatin ##REF##16882721##[45]##.</p> <p>In <italic>Arabidopsis thaliana</italic>, the <italic>VERNALIZATION INDEPENDENCE (VIP)</italic> genes; <italic>VIP2</italic> (now called <italic>ELF7</italic>), <italic>VIP4</italic>, <italic>VIP5</italic>, and <italic>VIP6/ELF8</italic> encode proteins closely related to Paf1, Leo1, Rtf1, and Ctr9, respectively ##REF##12207655##[46]##–##REF##15472079##[48]##, whereas VIP3 shows homology with hSki8. VIP3 physically interacts with VIP4 and VIP6 <italic>in vivo</italic> suggesting that these proteins comprise a complex analogous to Paf1C ##REF##15472079##[48]##. The <italic>VIP</italic> genes are required for proper expression of a common subset of genes including <italic>FLC</italic>, a MADS-box gene that represses the transition from vegetative growth to flowering. In <italic>vip</italic> mutants, <italic>FLC</italic> is ectopically silenced, allowing flowering soon after germination. <italic>FLC</italic> has emerged as a plant model for understanding the relationship between histone modifications and gene activity ##REF##17709278##[49]##. Activity of <italic>FLC</italic> early in development also requires the SET-domain proteins SDG8/EFS and ATX1, and is associated with methylation of H3K4 and H3K36 within <italic>FLC</italic> chromatin ##REF##16258034##[50]##–##REF##18375658##[52]##. Silencing of <italic>FLC</italic> in response to growth in cold temperatures (vernalization) is associated with loss of H3 acetylation and H3K4me3, and concomitant accumulation of H3K27me2/me3, within the <italic>FLC</italic> promoter and transcribed region ##REF##16957776##[27]##,##REF##14712277##[53]##,##REF##16682972##[54]##. The VRN2 protein, related to the PRC2 component Su(z)12, participates in K27 methylation at <italic>FLC</italic> and is required to maintain <italic>FLC</italic> silence in vernalized plants ##REF##14712277##[53]##,##REF##14712276##[55]##,##REF##16983073##[56]##.</p> <p>Unlike specific Paf1C components in yeast, null <italic>vip</italic> mutants (<italic>vip3/4/5/6</italic>) did not exhibit discernible reduction in the amount of H3K4me2/3 or H3K36me2 ##REF##15472079##[48]## when assayed at a whole-organism and whole-genome level, indicating that at least the bulk of these modifications is not dependent on Paf1C. However, H3K4me3 levels were reduced within <italic>FLC</italic> chromatin ##REF##17709278##[49]##. This loss of H3K4me3 could result indirectly from transcriptional inactivity, or could reveal a locus-specific role for plant Paf1C in mediating H3K4me3 deposition. To investigate potential mechanisms of histone modification and gene regulation modulated by Paf1C-related proteins in plants, and the relationship between Paf1C activity and H3K27 methylation/PcG-associated gene silencing, we mapped histone H3 modifications (trimethylation of H3K4/27, dimethylation of H3K36) and H3 occupancy in the entire Arabidopsis genome from wild-type and <italic>vip3</italic> mutant plants using ChIP combined with high-density tiling microarrays, and linked this information with Paf1C-dependent gene expression.</p>
<title>Materials and Methods</title> <title>Chromatin Immunoprecipitation</title> <p>Antibodies specific for histone H3 or H3 modifications were as follows: anti-H3-CT (Upstate, Lake Placid, NY; catalog no. 07-690), anti-H3 K4me3 antibody (Abcam, Cambridge, MA; catalog no. ab8580), anti-H3 K36me2 (Upstate 07-369), and anti-H3 K27me3 (Upstate 07-449). The specificity of these anti-H3 K4me3, anti-H3 K36me2, and anti-H3 K27me3 antibodies has been previously documented ##REF##17567990##[78]##–##REF##14690609##[80]##. ChIP was carried out essentially as previously described ##REF##15315638##[81]##, using two grams of aerial tissues from 14-d-old soil-grown plants. Immunoprecipitated DNA fragments were ∼500–600 bp and thus were expected to span two to three nucleosomal units.</p> <p>Chromatin immunoprecipitation followed by microarray analysis employed the Affymetrix GeneChip Arabidopsis Tiling 1.0R Array (Affymetrix, Santa Clara, CA), as described in the Affymetrix Chromatin Immunoprecipitation Assay Protocol. The Arabidopsis Tiling 1.0R Array represents ∼97% of the Arabidopsis genome with probes spaced every 35 bp. Signal intensities [perfect match (PM)-mismatch (MM)] from two independent biological replicates were quantile-normalized after log<sub>2</sub>-transformation using the TileMap package (<ext-link ext-link-type="uri" xlink:href="http://biogibbs.stanford.edu/jihk/TileMap/index.htm">http://biogibbs.stanford.edu/jihk/TileMap/index.htm</ext-link>) ##REF##16046496##[82]##. Subsequently, for each experiment, signals from immunoprecipitated (IP) or control (input) DNA were linearly scaled to the same mean. We computed a log ratio of the average IP to input value for each probe for further analyses. MvA plots and the correlation values (-R) of two replicates showed high reproducibility (R≥0.979) (##SUPPL##10##Figure S11## and ##SUPPL##13##Table S3##). To verify enrichments detected by microarray analysis, we carried out standard ChIP followed by semi-quantitative PCR for selected genes (##SUPPL##10##Figure S11## and data not shown).</p> <title>Derivation and Analyses of Gene-Level Modification Patterns</title> <p>For analysis of genic H3 occupancy or H3 modification profiles, we included only those genes spaced 350 bp or greater from an adjacent gene at the 5′ end, and 150 bp or greater from an adjacent gene at the 3′ end. This subset contained 14,485 protein-coding genes and 2,989 transposon-related/pseudogenes, from a total of 31,762 annotated nuclear genes. Gene annotations were taken from release 7 of The Arabidopsis Information Resource (TAIR) genome (<ext-link ext-link-type="ftp" xlink:href="ftp://ftp.arabidopsis.org/home/tair/Genes">ftp://ftp.arabidopsis.org/home/tair/Genes</ext-link>).</p> <p>Genic profiles were derived by analyzing probe signals for 100-bp windows within the proximal promoter (−350 to −50 relative to the TSS), TSS region (−49 bp to 0 bp to 5% of transcribed region), transcribed region (intervals of 10% of transcribed region from 5% to 95%), 3′ end region (from 95% to 100% of transcribed region to +50 bp relative to the 3′ end), and 3′ flanking region (51 bp to 150 bp relative to the 3′ end). To assess significance of differences in enrichment for H3 (##FIG##0##Figure 1##) or H3 modifications (##FIG##1##Figure 2##) within genes between wild-type and <italic>vip3</italic>, we treated positional signals within gene subsets as populations and computed P values using both Student's t-test and Wilcoxon rank sum test.</p> <p>To identify genomic regions substantially enriched for specific H3 modifications relative to H3 content as described in ##FIG##4##Figure 5##, probe-level t-statistics were computed for each probe based on ChIP vs. input. Neighboring probe signals were integrated by applying a hidden Markov model (HMM) to the probe level statistics with a maximal gap of 1000 bp, a minimal run of 200 bp, and posterior probability cutoff of 0.5. All procedures were performed using the TileMap package ##REF##16046496##[82]## and custom Perl scripts. As a result, combined enriched regions for H3K4me3 spanned a total of 18.4 Mbp; for H3K36me2, 19 Mbp; for H3K27me3, 18.2 Mbp. The identified genomic regions were filtered for annotated genes.</p> <title>Estimation of Gene Expression Levels and Tissue-Specificity in Wild-Type Plants</title> <p>AtGenExpress data sets 490, 491 and 492, corresponding to 21-, 22-, and 23-d-old whole plants, respectively (<ext-link ext-link-type="uri" xlink:href="http://www.arabidopsis.org/portals/expression/microarray/ATGenExpress.jsp">http://www.arabidopsis.org/portals/expression/microarray/ATGenExpress.jsp</ext-link>) were utilized to estimate gene expression levels in wild-type plants. Analysis using data set replicate 475, corresponding to 7-d-old seedlings, gave essentially identical results (data not shown). AtGenExpress data sets 469–547, corresponding to 79 samples representing various tissues and developmental stages, were used to estimate tissue-specificity of expression ##REF##15806101##[83]##. We computed Shannon entropy of genes as described ##REF##16949657##[84]##.</p> <title>Analysis of Expression and Chromatin Signatures of Genes Misregulated in <italic>vip3</italic> Mutants</title> <p>Expression of ∼22,600 genes was analyzed in wild-type and <italic>vip3-1</italic>
##REF##12750345##[57]## mutants using the Affymetrix ATH1 GeneChip. Data from CEL files were adjusted for background and normalized using the using the Bioconductor GCRMA package (<ext-link ext-link-type="uri" xlink:href="http://www.bioconductor.org">http://www.bioconductor.org</ext-link>). Statistically significant changes (p&lt;0.001 and 2-fold change) in gene expression between wild type and <italic>vip3</italic> were detected using the Bioconductor LIMMA package ##UREF##0##[85]##. Of 218 upregulated genes and 241 downregulated genes in <italic>vip3</italic> relative to wild type, 139 (upregulated) and 159 (downregulated) were also included in the gene set evaluated for chromatin modifications. To assess statistical significance of differences in chromatin signatures for Paf1C-regulated genes between wild-type and <italic>vip3</italic> mutants, as shown in ##SUPPL##6##Figure S7## we computed the 95<sup>th</sup> percentile confidence intervals for differences in mean positional signals within 1,000 randomly resampled gene sets, each containing 139 (for upregulated) or 159 (for downregulated) genes. To assess the statistical significance of differences in genic chromatin signatures for Paf1C-dependent genes relative to typical genic chromatin signatures as shown in ##SUPPL##7##Figure S8##, we computed the 95th percentile confidence intervals for the mean positional signals within these gene sets. To assess significance of the lower entropy values observed in genes misregulated in <italic>vip3</italic> mutants, we generated random gene combinations with mean wild-type expression level values similar to those of the upregulated or downregulated gene sets (7.2–7.4 and 7.8–8.0, respectively), and used these populations to compute P values using both Student's t-test and Wilcoxon rank sum test.</p> <title>Accession Numbers</title> <p>Raw data from these experiments has been deposited in the NCBI Gene Expression Omnibus (GEO), accession number GSE7907 (genomic tiling arrays) and GSE10928 (expression arrays).</p>
<title>Results</title> <title>Mapping of Wild-Type and Paf1C-Dependent H3 Occupancy and H3 Modifications in the Arabidopsis Genome</title> <p>To investigate the influence of plant Paf1C activity on histone modifications in plants, we mapped H3 occupancy and distribution of specific histone H3 methylations at high resolution throughout the genome of both wild-type Arabidopsis plants and mutants homozygous for a null allele of the Paf1C-related gene <italic>VIP3</italic>. We targeted H3K4me3 and H3K36me2 because previous analyses showed that both modifications are associated with transcriptional activity of the flowering regulatory gene <italic>FLC</italic>
##REF##15520273##[47]##,##REF##16299497##[51]##, which is silenced in <italic>vip</italic> mutants ##REF##12750345##[57]##, and because of the observation that methylation of H3K4 and H3K36 are associated with Paf1C activity in budding yeast ##REF##12667454##[7]##,##REF##12667453##[8]##,##REF##12773564##[15]##. We also analyzed H3K27me3, as this modification was reported to be associated with <italic>FLC</italic> silencing in vernalized plants ##REF##16957776##[27]##.</p> <p>We first estimated chromatin occupancy by H3 using an antibody (H3-CT) specific for the carboxyl-terminal domain. Consistent with previous reports of H3 occupancy in Arabidopsis and other species ##REF##16122420##[58]##–##REF##17439305##[60]##, when mean positional signals were calculated for a set of 17,771 annotated genes from a variety of classes (see <xref ref-type="sec" rid="s4">Materials and Methods</xref>), a pattern of mean H3 signal was evident, characterized by enrichment within transcribed regions concomitant with depletion within both promoter and 3′ regions, relative to the genomic median (##FIG##0##Figure 1A##). Protein-coding genes generally showed lower signals than transposon-related genes or pseudogenes. Also consistent with previous reports, we found a clear association between estimated transcriptional activity (see <xref ref-type="sec" rid="s4">Materials and Methods</xref>) and H3 signal depletion within the proximal promoter/TSS (##FIG##0##Figure 1B##).</p> <p>When mean positional signals were calculated for the 17,771 gene set for <italic>vip3</italic> mutant plants relative to those for wild-type plants, H3 signals were significantly lower (P value&lt;0.0006; Student's t-test and Wilcoxon rank sum test) across the transcribed region (##FIG##0##Figure 1C##). The ∼3,000 transposon-related and pseudogenes analyzed within this set showed a slight, but insignificant, increase in H3 signals across the transcribed region (##FIG##0##Figure 1D##). We sorted protein coding genes into ten-percentile bins based on expression level (see <xref ref-type="sec" rid="s4">Materials and Methods</xref>), and analyzed positional signals for <italic>vip3</italic> mutants relative to wild-type plants within each bin. Interestingly, H3 signals were lower in <italic>vip3</italic> plants within transcribed regions for the five highest expression bins, with the two top expression bins showing a significant (P&lt;0.0001) relative loss of H3 (##FIG##0##Figure 1E##). When Paf1C-dependence of H3 occupancy was considered as a function of gene length, we noted a subtle relationship between gene length and degree of depletion in <italic>vip3</italic> plants relative to wild-type, with the shortest genes analyzed showing little or no relative H3 loss, and genes with transcribed regions &gt;1 kb in length showing significant (P&lt;0.0001) H3 depletion throughout the transcribed region (##FIG##0##Figure 1E##). Although these effects were slight when analyzed on a genome level, our results suggests plant Paf1C may have transcription-dependent activity in maintaining H3 and/or nucleosomal density, especially in long genes.</p> <p>For subsequent analysis of H3K4me3, H3K36me2, or H3K27me3, we expressed signal values relative to total H3 signal. For H3K4me3 and H3K36me2, signals showed a similar pattern at the chromosomal level, being generally above the genomic median in gene-rich chromosome arms and below the genomic median in heterochromatic regions (##SUPPL##0##Figure S1## and data not shown). Considering only protein-coding genes, H3K4me3 increased in the promoter region, peaked near the TSS, and decreased throughout most of the body of the genes (##FIG##1##Figure 2A##). In contrast, H3K36me2 showed an increase throughout the transcribed regions, peaking near the end of the transcribed regions. Transposon-related/pseudogenes showed depletion of both modifications throughout the extent of the transcribed regions (##FIG##1##Figure 2A##). Further analysis of genes within expression percentile bins revealed a striking relationship between the enrichment for H3K4me3 or H3K36me2 and transcriptional activity. For example, those genes included in Bin 10, which exhibited the top 10% of expression values, also exhibited the highest peak of H3K4me3 modification near the TSS, whereas those genes with the lowest expression levels showed no obvious peak (Bins 1–3) or only a subtle peak (Bin 4) (##FIG##1##Figure 2B##). For H3K36me2, genes included in Bins 8 and 9 showed the strongest 3′ peaks, Bin 10 genes showed slightly lower peaks, and genes in Bins 1–4 exhibited low levels of signal throughout the transcriptional unit with no apparent 3′ peak (##FIG##1##Figure 2B##). Peaks of these two H3 modifications also varied with gene length. When H3K4me3 and H3K36me2 were plotted for genes representing each of five length-assigned bins, the longest genes (Bins 4 and 5, containing genes with transcribed regions &gt;3 kb) displayed the strongest H3K4me3 and H3K36me2 signals at the 5′ and 3′ ends, respectively, whereas Bin 2 (1–2 kb) showed only weak peak signals (##FIG##1##Figure 2C##). For genes &gt;2 kb in length, the peak of H3K4me3 occurred in a ∼1-kb region at the 5′ end of genes, independent of gene length, whereas the H3K36me2 peak occupied the 3′ ∼one-half of genes, being much broader in longer genes (##FIG##1##Figure 2C##). These genic patterns and relationship with transcription are similar to those found previously for these modifications in yeast, humans, and where studied, in plants ##REF##16122420##[4]##,##REF##15680324##[5]##,##REF##12667454##[7]##,##REF##12773564##[15]##,##REF##15175259##[61]##,##REF##18263775##[62]## and are consistent with an evolutionarily conserved role for H3K4me3 marking transcriptional engagement, and H3K36me2 as a mark of transcriptional elongation ##REF##12667454##[7]##,##REF##14636589##[12]##,##REF##12773564##[15]##.</p> <p>When signal values were averaged over extended euchromatic or heterochromatic regions, the abundance of H3K4me3 and H3K36me2 was not perturbed in <italic>vip3</italic> plants relative to wild-type (##SUPPL##1##Figure S2## and data not shown). This is consistent with our previous finding that disruption of plant Paf1C did not affect total cellular levels of H3K4me3 or H3K36me2 ##REF##15472079##[48]##. However, the effects of loss of Paf1C on distribution of both H3K4me3 and H3K36me2 within genes were substantial: H3K4me3 signals from <italic>vip3</italic> were significantly lower than in wild-type within the peak of this modification near the TSS, but were elevated in the 3′ ∼one-half of the transcribed region (P&lt;2E-16), whereas H3K36me2 signals from <italic>vip3</italic> were significantly lower within the 3′ half, and elevated within the 5′ half of the transcribed region (P&lt;2E-16) (##FIG##1##Figure 2E##). The most highly expressed genes showed the greatest relative loss of H3K4me3 in 5′ regions, gain in H3K4me3 in 3′ regions, gain of H3K36me2 in 5′ regions, and loss of H3K36me2 in 3′ regions in <italic>vip3</italic> plants (P&lt;0.001 for bins 7–10) (##FIG##1##Figure 2G##). Transposon-related genes and pseudogenes showed a significant loss of H3K4me3 across the transcribed region (P&lt;1E-16) (##FIG##1##Figure 2F##).</p> <p>We also found a substantial length-associated relative increase in H3K4me3 in genes ≥2 kb in length, extending from ∼1 kb downstream of the TSS to the 3′ end (P&lt;2E-06) (##FIG##1##Figure 2H##). H3K36me2 showed more moderate changes within domains proportional to gene length: an increase within the 5′ ∼one half and decrease within the 3′ ∼one-half of the transcribed region (P&lt;0.005 for genes ≥2 kb in length) (##FIG##1##Figure 2H##). To analyze this effect independently of any potential relationship between gene length and expression level, we considered differences in relative H3K4me3 or H3K36me2 enrichment within restricted subsets of genes showing similar expression levels. Genes comprising the top 10% expression level bin showed a striking relationship between 3′ H3K4me3 enrichment/H3K36me2 depletion and length (##SUPPL##2##Figure S3##). Weakly expressed genes (Expression Level Bin 2) did not show this relationship. We interpret this data as showing that plant Paf1C is required globally not only to maintain a ∼1-kb 5′ peak of H3K4me3 and 3′ enriched region of H3K36me2, respectively, but also for exclusion of these modifications from the remainder of the transcribed region. Because the extent of 5′ H3K4me3 and 3′ H3K36me2 enrichment is generally related to expression level ##REF##16122420##[4]##,##REF##12667454##[7]##,##REF##12773564##[15]##, see above, these data imply that the role of Paf1C to maintain appropriate patterns of these modifications is linked to transcriptional activity.</p> <p>Our analyses of H3K27me3 distribution in wild-type plants support recent reports detailing the global pattern of this modification in Arabidopsis ##REF##17542647##[28]##,##REF##17439305##[29]##,##REF##15457214##[63]##. We found generally strong signals along euchromatic chromosome arms and relatively weak signals within heterochromatic regions (##SUPPL##0##Figure S1## and data not shown) with domains of enrichment occupying ∼8,000 of ∼32,000 annotated Arabidopsis genes (see below), and encompassing the transcribed regions of genes known to be subject to PcG repression, including <italic>FLC</italic> (##SUPPL##3##Figure S4##; see below). At the genic level, H3K27me3 was relatively enriched near the TSS and 3′ ends, with weak signals seen across the transcribed region (##FIG##2##Figure 3A##). Like H3K4me3 and H3K36me2, levels of H3K27me3 were much higher across protein-coding genes than transposon-related and pseudogenes. Similar to the previously reported findings of Zhang et al. ##REF##17439305##[29]##, and in contrast to our results obtained for H3K4me3 and H3K36me2, H3K27me3 was highest in those genes exhibiting the lowest expression levels (##FIG##2##Figure 3B##). Analysis of H3K27me3 within gene length bins also revealed a negative relationship between levels of this modification and gene length, with the greatest degree of depletion in the longest genes, and relative enrichment in the shortest genes (##FIG##2##Figure 3C##). Our results also support the finding of Zhang et al. ##REF##17439305##[29]## that genes with low Shannon expression entropy values (tending toward very specific expression patterns) tended to be highly enriched in H3K27me3, whereas genes with high entropy values (widespread expression) generally showed very low H3K27me3 signals (##FIG##2##Figures 3D## and ##SUPPL##4##S5##). In contrast, our results do not reveal a strong relationship between entropy and modification by H3K4me3 or H3K36me2 (##FIG##1##Figures 2D## and ##SUPPL##4##S5##).</p> <p>At the chromosomal level, H3K27me3 signals were not noticeably disrupted in <italic>vip3</italic> plants (##SUPPL##1##Figure S2## and data not shown). At the genic level, H3K27me3 signals were essentially unchanged across the transcriptional unit (##FIG##2##Figure 3E## and data not shown). Highly expressed genes also showed relative gains in H3K27me3 signal within transcribed regions in the <italic>vip3</italic> mutant, but the magnitude of these gains was small compared to the effect seen with H3K4me3 and H3K36me2 (##FIG##2##Figure 3G##). We saw no obvious gene length-related dependence on Paf1C for H3K27me3 (##FIG##2##Figure 3H##).</p> <p>In addition, we examined the relationship between Paf1C-dependence of H3 occupancy or distribution of H3 modifications and tissue-specificity of expression. H3 occupancy was decreased in <italic>vip3</italic> plants within genes showing the highest entropy values (tending toward ubiquitous expression patterns) (##SUPPL##5##Figure S6##). This was not unexpected, because genes with highest expression values, which show the greatest depletion of H3, also tend to be ubiquitously expressed (see below). Similarly, entropy-associated profiles of changes in H3 modifications in <italic>vip3</italic> plants were similar to those seen when analyzed for expression level (##SUPPL##5##Figures S6##, ##FIG##1##2B## and ##FIG##2##3B##). Enrichment for H3K4me3 and H3K36me2, and to a lesser extent H3K27me3, is dependent on expression in subsets of genes with similar entropy values (##FIG##1##Figures 2D## and ##FIG##2##3D##). We did not observe a convincing relationship between Paf1C-dependence of H3 occupancy or modification enrichment and entropy when genes within similar expression windows were considered (data not shown). Thus the apparent entropy-associated changes in these modifications dependent on Paf1C may be driven largely by levels of gene expression.</p> <title>Mapping of Paf1C-Dependent H3 Modifications within Paf1C-Targeted Genes</title> <p>We then considered Paf1C-dependent changes in H3 modifications within those genes whose normal expression depends on Paf1C. Gene expression profiling utilizing microarrays representing most canonical Arabidopsis genes identified a small subset of genes, including the previously identified Paf1C target <italic>FLC</italic>, that were regulated by <italic>VIP3</italic> (data not shown). We observed a statistically significant, average loss of H3K4me3 and H3K36me2 across most of the extent of downregulated genes (##FIG##3##Figures 4A## and ##SUPPL##6##S7##). Because enrichment for H3K4me3 and H3K36me2 is correlated with degree of gene expression, this result is consistent with the expected changes in these modifications associated with decreased gene activity. This reveals at least an indirect role for Paf1C in mediating these modifications. For H3K27me3, we observed a slight decrease across upregulated genes, and increase across downregulated genes, but this was not found to be statistically significant (##FIG##3##Figures 4A## and ##SUPPL##6##S7##).</p> <p>Expression of <italic>FLC</italic> is promoted through a mechanism involving plant Paf1C, but is also subject to developmental silencing through a PcG-like mechanism that includes the Su(z)12-like protein VRN2 and accumulation of H3K27me2/3 within the <italic>FLC</italic> gene ##REF##16957776##[27]##,##REF##14712277##[53]##,##REF##16983073##[56]##. To explore the relationship between Paf1C and activating or silencing modifications at <italic>FLC</italic>, we analyzed the chromatin profile of <italic>FLC</italic> in wild-type plants, and the differences in chromatin profiles between wild-type and <italic>vip3</italic> mutant plants. In wild-type plants, H3K4me3 showed a pronounced peak near the TSS and beginning of the first intron (##FIG##3##Figure 4B##). H3K36me2 showed relatively low levels throughout the <italic>FLC</italic> gene, increasing slightly through the transcribed region and peaking near the 3′ end. Substantial H3K27me3 was seen throughout most of the transcribed region (##FIG##3##Figure 4B##). In <italic>vip3</italic> mutants, similar to the effect of loss of <italic>VIP3</italic> on the average signal in protein-coding genes, the 5′ peak of H3K4me3 was reduced and levels of H3K4me3 increased in more 3′ regions. H3K36me2 decreased further throughout most of the transcribed region, including the 3′ end (##FIG##3##Figure 4B##). These observations are consistent with those of Xu et al. ##REF##18070919##[64]## who found decreases in both H3K4m3 and H3K36me2 at the 5′ end of <italic>FLC</italic> in plants dysfunctional for the Paf1C-related factor VIP4. In striking contrast, H3K27me3 increased substantially within the proximal promoter, TSS, and the 3′ ∼one-half of the gene (##FIG##3##Figure 4B##). Thus, at <italic>FLC</italic>, loss of expression is associated with chromatin changes both typical (loss of H3K4me3 and H3K36me2) and atypical (substantial gain of H3K27me3) for Paf1C-regulated genes.</p> <title>Chromatin and Expression Characteristics of Paf1C-Targeted Genes</title> <p>Common distinguishing features of genes dependent on Paf1C for appropriate activity have not been identified. To explore the involvement of chromatin structure in predisposing genes to regulation by Paf1C, we examined the wild-type pattern of H3 modifications among genes that were misexpressed in Paf1C mutants. We found that these genes showed several unique chromatin signatures when compared with the entire set of protein-coding genes. Strikingly, genes misregulated in <italic>vip3</italic> mutants showed much greater enrichment for H3K27me3 across most of the transcribed region (upregulated genes) or the entire transcriptional unit (downregulated genes), relative to average levels for the entire gene set (##FIG##3##Figures 4A## and ##SUPPL##7##S8##). These genes were also typified by significantly greater H3K4me3 enrichment throughout much of the transcribed region, with a peak of enrichment 3′ to that seen for the typical gene. Wild-type levels of H3K36me2 were higher throughout the transcribed region, especially near the 5′ end of the transcribed region where levels in the typical gene are lowest. Levels of these three modifications were not dramatically different from the typical gene in the promoter and 3′ regions (##FIG##3##Figures 4A## and ##SUPPL##7##S8##). We also analyzed the chromatin signatures of the subsets of gene that we previously found to be misregulated in loss-of-function mutants for the <italic>VIP5</italic> and <italic>VIP6/ELF8</italic> genes, encoding plant homologs of the Paf1C components Rtf1 and Ctr9, respectively ##REF##15520273##[47]##,##REF##15472079##[48]##. As expected from the observation that the subsets of misregulated genes in <italic>vip5</italic> or <italic>vip6</italic> were largely overlapping with that of <italic>vip3</italic> (data not shown), the chromatin signatures of <italic>VIP5</italic>- or <italic>VIP6</italic>-regulated genes were similarly characterized by a conspicuous enrichment for H3K27me3 across the transcribed region (##SUPPL##8##Figure S9##). These genes also showed enhanced H3K4me3 in the 3′ transcribed region, and elevated H3K36me2 in the 5′ transcribed region; this was most apparent for genes downregulated in the mutants (##SUPPL##8##Figure S9##).</p> <p>The finding that Paf1C targeted genes were typically distinguished by combinatorial enrichment for H3K27me3, H3K4me3, and H3K36me2 was intriguing, because we found that H3K4me3, and to a lesser extent H3K36me2, co-occupies only a small subset of genes with H3K27me3 (##FIG##4##Figure 5##). Indeed, when assigned to groups defined by substantial enrichment for each H3 modification within the transcriptional unit, genes whose expression are positively or negatively influenced by <italic>VIP3</italic> were most significantly overrepresented within a group of genes strongly enriched for both H3K4me3 and H3K27me3 [P value &lt;1E-10 or &lt;1E-05, respectively; Fisher's exact test) (##FIG##4##Figure 5## and ##SUPPL##11##Table S1##). Additionally, when genes were clustered based on distinctions in genic profiles of the modifications as well as enrichment levels, Paf1C-regulated genes were overrepresented in a group of genes showing strong enrichment for H3K4me3 and H3K27me3 and moderate enrichment for H3K36me2, and with H3K4me3 and H3K36me2 distributed broadly across the transcriptional unit rather than in discrete 5′/3′ peaks (##SUPPL##9##Figure S10##, ##SUPPL##12##Table S2## and data not shown).</p> <p>High levels of H3K27me3 mark transcriptionally quiescent and developmentally silenced genes, including known targets of plant PcG-like machinery ##REF##17542647##[28]##,##REF##17439305##[29]##; above. To determine if plant Paf1C has a special role in the regulation of such genes, we compared expression level and entropy of Paf1C-regulated genes with that of the entire gene set. We found that, in wild-type plants, those genes strongly upregulated or downregulated in the Paf1C mutants tended to show low wild-type entropy values, even relative to genes expressed to similar levels (##FIG##5##Figure 6## and ##SUPPL##8##Figure S9##). Scatter plots of expression level and entropy for Paf1C regulated genes, in the context of the entire gene set, clearly showed that within a specific expression level, Paf1C regulated genes tended to show lower entropy values; this was especially significant for genes downregulated in the mutant (P&lt;0.01 and P&lt;0.05 for downregulated or upregulated genes, respectively; Wilcoxon signed-rank test and Student's t-test; see <xref ref-type="sec" rid="s4">Materials and Methods</xref>) (##FIG##5##Figures 6B## and ##SUPPL##8##S9##). This is not a trivial result of a potential tissue-specific expression pattern of these genes, as they are expressed relatively ubiquitously ##REF##15472079##[48]##; Expression Entropy Bins 8–10 (data not shown). Taken together, these observations suggest that plant Paf1C has an important role in maintaining appropriate expression of developmentally regulated genes.</p>
<title>Discussion</title> <title>Global and Locus-Specific Roles for Plant Paf1C in Chromatin Modifications</title> <p>In accordance with our earlier report that loss of the plant Paf1C subunits VIP3, VIP4, VIP5, or VIP6 did not affect total cellular levels of H3K4me3 or H3K36me2, here we show that the overall abundance of these modifications within chromatin is not obviously altered in a <italic>vip3</italic> mutant. Instead, loss of plant Paf1C led to redistribution of these modifications within genes: a 3′ shift in the distribution of H3K4me3 and 5′ shift in the distribution of H3K36me2. In yeast, the spatially restricted pattern of genic methyl-K4 and methyl-K36 is thought to depend on the transition from Ser5-phosphorylation to Ser2-phosphorylation within the heptapeptide repeat of the PolII CTD and recruitment of Set1 and Set2, a mechanism that requires Paf1C ##REF##15489290##[65]##. If an analogous mechanism linking PolII with these chromatin modifications exists in plants, then a plausible explanation for the general spreading of both H3K4me3 and H3K36me2 within genes in <italic>vip3</italic> mutants is an irregular transition from the Ser-5 to Ser-2 phosphorylated form of PolII. For example, plant Paf1C could be required for interaction of CTD phosphatases with PolII, leading to accumulation of hyperphosphorylated PolII when dysfunctional.</p> <p>We postulate several scenarios by which this disruption of H3K4me3 and H3K36me2 patterning could directly affect gene activity. Reduction in H3K4me3 near the promoter and TSS may disrupt recruitment of chromatin remodeling machinery needed for efficient initiation, as has been demonstrated for the NURF complex at a Hox promoter ##REF##16728976##[13]## or lead to defective pre-mRNA processing ##REF##18042460##[66]##. Enhanced levels of H3K4me3 within transcribed regions may promote transcriptional initiation at cryptic sites; the aberrant RNAs thus formed could trigger gene repression through RNAi-related pathways. Loss of H3K36 methylation within transcribed regions may similarly promote aberrant initiation ##REF##16286007##[18]##–##REF##16286008##[20]##, and RNAi-related silencing, or may alter elongation ##REF##12773564##[15]##,##REF##12511561##[67]##. Enhanced H3K36 methylation within the 5′/TSS region may disrupt initiation ##REF##11839797##[68]##,##REF##16977311##[69]##.</p> <p>The apparent global role for Paf1C in maintaining methyl K4/K36 patterns in plants is consistent with observations that budding yeast Paf1C components are abundant and ubiquitously associated with promoters and open reading frames ##REF##15149594##[36]##,##REF##11983171##[70]##. Interestingly, however, only a small subset of genes are misregulated in yeast or plant Paf1C mutants ##REF##16246724##[35]##,##REF##15472079##[48]##,##REF##9032243##[71]##,##REF##12455700##[72]##; this study. What are the distinguishing features of genes whose expression is dependent on plant Paf1C activity? Here, we showed that genes either positively or negatively regulated by plant Paf1C were generally enriched in H3K4me3, H3K27me3, and H3K36me2. We observed some distinction in the degree and pattern of enhanced H3K4me3 and H3K36me2 enrichment between those genes misregulated in <italic>vip3</italic>, and those genes misregulated in <italic>vip5</italic> or <italic>vip6</italic>. This may be attributed to the fact that our gene expression information for <italic>vip5</italic> and <italic>vip6</italic> was archival and derived from plants of a slightly different developmental stage. The observation that both up- and down-regulated genes show similar chromatin profiles may be explained by the fact that some of these genes may be targeted by Paf1C only indirectly. If plant Paf1C primarily targets developmentally regulated genes, then genes positively or negatively regulated by these genes would also be expected to show developmental regulation and accordant chromatin signatures.</p> <p>The apparent co-occurrence of H3K4me3/H3K36me2 and H3K27me3 domains seen in these genes may result from the net observation of distinct chromatins marked predominantly by H3K4me3/H3K36me2 (in cells where the gene is mostly active) or by H3K27me3 (in cells where the gene is repressed). Another possibility, not mutually exclusive, is that these modifications could be juxtaposed within contiguous chromatin, as seen for H3K4me3 and H3K27me3 in numerous developmentally important genes in mammalian embryonic stem (ES) cells ##REF##16630819##[73]##. In ES cells, this so-called H3K4me3/H3K27me3 bivalent domain has been considered as a mechanism to facilitate switching from a repressed to active state, and can resolve to a H3K4me3-dominated or H3K27me3-dominated signature in differentiated cells where the locus is active or repressed, respectively ##REF##16630819##[73]##,##REF##17603471##[74]##.</p> <title>Chromatin Dynamics at the <italic>FLC</italic> Gene</title> <p>The MADS-box gene <italic>FLC</italic> is promoted through a mechanism involving Paf1C during early plant development, and is targeted for repression by a PRC2-like mechanism in response to cold. In wild-type plants, the <italic>FLC</italic> locus showed H3 modification profiles typical for Paf1C-regulated genes: high levels of H3K4me3 at the 5′ end, relatively low enrichment for H3K36me2 at the 3′ end, and a domain of H3K27me3 enrichment throughout much of the transcribed region. However, unlike other genes whose expression is promoted by plant Paf1C, <italic>FLC</italic> exhibited a substantial further accumulation of H3K27me3 when silenced in mutant plants. This suggests a role for Paf1C in antagonizing PcG repression at this locus. Plant Paf1C may also function to antagonize silencing of the several additional known PcG targets, including genes with homeotic functions in flower development, and this could explain the misregulation of these genes and floral abnormalities seen in mutants for various Paf1C-related genes ##REF##12750345##[57]##.</p> <p>How might such antagonism be mediated? One of several possibilities is that a role for Paf1C in linking transcription with H3K4/K36 methylation may be elaborated in higher eukaryotes through transcription-associated histone replacement, in which canonical H3 assembled into nascent chromatin is exchanged for ‘variant’ histone H3 (H3.3, also called H3.2 in plants) ##REF##12086617##[75]## known to be enriched for methyl-K4 and/or methyl-K36 ##REF##15598823##[76]##,##REF##14732680##[77]##. Random distribution of H3.3 nucleosomes during replication of active loci would result in a relatively high proportion of methyl-K4-modified H3.3 in nascent chromatin. This content may be further increased during pioneering rounds of transcription through Set1-like H3K4 methyltransferase activity, effectively resetting chromatin to the active state. In contrast, nascent chromatin at silenced loci is expected to be enriched for nucleosomes containing canonical H3, known to be preferentially modified by methyl-H3K27 ##REF##15598823##[76]##. H3K27me3 occupancy may be actively reinforced, or passively sustained by modification of canonical H3 in nascent chromatin upon successive replication events. Disrupting Paf1C, and thus the linkage between transcription and H3K4 methylation, would negatively influence resetting of chromatin to the active state and shift the balance of modification to H3K27me3. For some genes, such as <italic>FLC</italic>, even a small disruption of such a balance may then have qualitative effects on chromatin structure and expression.</p>
<p>Conceived and designed the experiments: SO SP SvN. Performed the experiments: SO SP. Analyzed the data: SO SP SvN. Wrote the paper: SO SP SvN.</p> <p>In budding yeast, intragenic histone modification is linked with transcriptional elongation through the conserved regulator Paf1C. To investigate Paf1C-related function in higher eukaryotes, we analyzed the effects of loss of Paf1C on histone H3 density and patterns of H3 methylated at K4, K27, and K36 in Arabidopsis genes, and integrated this with existing gene expression data. Loss of Paf1C did not change global abundance of H3K4me3 or H3K36me2 within chromatin, but instead led to a 3′ shift in the distribution of H3K4me3 and a 5′ shift in the distribution of H3K36me2 within genes. We found that genes regulated by plant Paf1C showed strong enrichment for both H3K4me3 and H3K27me3 and also showed a high degree of tissue-specific expression. At the Paf1C- and PcG-regulated gene <italic>FLC</italic>, transcriptional silencing and loss of H3K4me3 and H3K36me2 were accompanied by expansion of H3K27me3 into the promoter and transcriptional start regions and further enrichment of H3K27me3 within the transcribed region. These results highlight both genic and global functions for plant Paf1C in histone modification and gene expression, and link transcriptional activity with cellular memory.</p> <title>Author Summary</title> <p>In eukaryotes, DNA is packaged with histones and other proteins into a dynamic fabric called chromatin. Specific modifications of histones—including methylation of key lysine residues—provide genetic information that acts synergistically with the DNA code. In yeast, the conserved transcriptional regulator Paf1C is required for marking histone H3 within active genes by methylation of Lysine-4, a modification thought to promote gene activity. In higher eukaryotes, this mechanism is elaborated through Polycomb-Group (PcG), which maintains transcriptional repression through cell divisions and involves methylation of Lysine-27 of H3. In this study, we mapped these and other key H3 modifications throughout the genome of the plant <italic>Arabidopsis thaliana</italic> and evaluated the effects of loss of Paf1C on these modifications and gene expression. We found that Paf1C acts globally to maintain histone modification within genes, but is required for appropriate expression of only a handful of genes. These typically showed a high degree of developmental regulation in both Lysine-4 and Lysine-27 methylation. At the flowering regulator <italic>FLC</italic>, targeted by both Paf1C and PcG, loss of activating (Lysine-4) methylation was accompanied by further accumulation of repressive (Lysine-27) methylation. These results provide a link between transcriptional activity and cellular memory.</p>
<title>Supporting Information</title>
<p>We thank members of the Gene Expression in Development and Disease group (Michigan State University) for helpful discussion.</p>
<fig id="pgen-1000077-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000077.g001</object-id><label>Figure 1</label><caption><title>H3 Signal Profiles within Genes and Effects of Paf1C Disruption.</title><p>(A, B) H3 signal profiles associated with gene type and expression level. (A) Mean genic positional signal for H3 was calculated independently for ∼14,500 likely protein-coding genes and ∼3,000 transposon-related or likely pseudogenes from our ∼18,000 gene set, and is depicted across the promoter regions (shown in bp from −300 to 0 relative to the presumed transcriptional start site), transcribed regions (shown proportionally from 0 to 100% of total length), and 3′ regions (shown in bp from 0 to +100 relative to the presumed 3′ end). (B) Protein-coding genes were sorted into ten-percentile bins according to their expression level, as estimated from publicly available microarray data (see <xref ref-type="sec" rid="s4">Materials and Methods</xref>). The averaged positional signals for H3 within each bin were plotted. Data is shown for absolute positions across the 5′ end including the presumed transcriptional start site (TSS). (C–E) Genic patterns of H3 signals associated with Paf1C activity. (C) In the left column, signals for all 17,771 genes evaluated are depicted for wild-type plants (WT) or <italic>vip3</italic> mutants across the transcriptional unit as described above for ##FIG##0##Figure 1A##. In the center and right columns, data is shown for absolute positions across the 5′ end including the presumed transcriptional start site (TSS) (center), or across the 3′ end (right) of a subset of 6,180 genes with transcribed regions &gt;2 kb in length. (D) Mean positional H3 signals for transposon-related genes and pseudogenes. (E) Genic patterns of Paf1C-dependent H3 signals with respect to expression level and length. Genic positional signals for H3 were averaged separately for protein coding genes within ten-percentile bins according to expression level (left panel), or twenty-percentile bins according to length of transcribed region (right panel) for <italic>vip3</italic> plants relative to wild-type.</p></caption></fig> <fig id="pgen-1000077-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000077.g002</object-id><label>Figure 2</label><caption><title>Enrichment of H3K4me3 and H3K36me2 within Genes and Effects of Paf1C Disruption.</title><p>(A) Mean genic positional enrichment for H3K4me3 or H3K36me2 was calculated independently for protein-coding genes or transposon-related/pseudogenes as described above for ##FIG##0##Figure 1A##. (B) H3K4me3 or H3K36me2 enrichment is depicted for genes within ten-percentile expression level bins as described for ##FIG##0##Figure 1B##. (C) H3K4me3 or H3K36me2 enrichment is depicted for genes within twenty-percentile bins according to length of transcribed region. Enrichment is shown across the transcriptional unit (left panels) or the TSS/5′end (for H3K4me3) or 3′ end (for H3K36me2) (right panels) (D) Protein-coding genes were assigned to ten bins according to tissue-specificity of expression, as estimated by Shannon entropy (see <xref ref-type="sec" rid="s4">Materials and Methods</xref>). Genes in Bin 10 (high entropy) show the most ubiquitous expression across various plant parts, whereas genes in Bin 1 (low entropy) show the most specific expression domains. Mean positional signals were calculated for genes within specific ten-percentile expression (Exp) and entropy (Ent) bins, as indicated. Lines were smoothed using a three-point sliding window. (E) Enrichment for all 17,771 genes evaluated are depicted across the transcriptional unit (left panel), or across the 5′ end/TSS (center) or 3′ end (right) of genes with transcribed regions &gt;2 kb in length, for wild-type plants (WT) or <italic>vip3</italic> mutants. (F) Enrichment within transposon-related genes and pseudogenes. (G) Paf1C-dependent H3K4me3 or H3K36me2 enrichment with respect to expression level, as determined for ##FIG##0##Figure 1E##. (H) Paf1C-dependent enrichment for H3K4me3 or H3K36me2 with respect to gene length, as determined for ##FIG##0##Figure 1E##.</p></caption></fig> <fig id="pgen-1000077-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000077.g003</object-id><label>Figure 3</label><caption><title>Enrichment of H3K27me3 within Genes and Effects of Paf1C Disruption.</title><p>(A) Mean genic positional enrichment for H3K27me3 for protein-coding genes or transposon-related/pseudogenes. (B) H3K27me3 enrichment for genes within ten-percentile expression level bins as described for ##FIG##0##Figure 1B##. (C) Mean genic positional enrichment for genes within twenty-percentile bins according to length of transcribed region. (D) Enrichment for H3K27me3 for genes within specific expression level and expression entropy bins as described for ##FIG##1##Figure 2D##. (E) Mean enrichment across the transcriptional unit (left panel) for the ∼18,000-gene set, or across the 5′ end/TSS (center) or 3′ end (right) of genes with transcribed regions &gt;2 kb in length, for wild-type plants (WT) or <italic>vip3</italic> mutants. (F) Enrichment for transposon-related genes and pseudogenes. (G) Paf1C-dependent H3K27me3 enrichment with respect to expression level, as determined for ##FIG##0##Figure 1E##. (H) Paf1C-dependent enrichment for H3K27me3 with respect to gene length, as determined for ##FIG##0##Figure 1E##.</p></caption></fig> <fig id="pgen-1000077-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000077.g004</object-id><label>Figure 4</label><caption><title>H3 Methylation Profiles and Paf1C-Dependent H3 Methylation Profiles in Paf1C-Targeted Genes.</title><p>(A) Genic positional signals for H3 lysine methylations as indicated were averaged separately for genes upregulated in <italic>vip3</italic> mutants (top row of panels, red) or downregulated in <italic>vip3</italic> mutants (lower panels, green) for both wild-type plants (solid lines) and <italic>vip3</italic> mutants (dashed lines). Averaged signals for all genes in the protein-coding gene set presented in ##FIG##1##Figures 2## and ##FIG##2##3## are shown in black. (B) Signals for H3 lysine methylations are shown within a ∼14-kb region encompassing the plant Paf1C-dependent gene <italic>FLC</italic> from wild-type (WT) plants (top panel). Lower panels show the relative difference in signals between <italic>vip3</italic> and wild-type chromatin. Horizontal colored bars in these panels indicate regions where significant (2.5-fold change in <italic>vip3</italic>/WT; P&lt;10<sup>−3</sup> in either <italic>vip3</italic> or WT) differences in signals were observed. A depiction of the <italic>FLC</italic> gene within this region is shown at bottom. Data depicted in this figure were corrected for total H3.</p></caption></fig> <fig id="pgen-1000077-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000077.g005</object-id><label>Figure 5</label><caption><title>Relationships Among H3 Lysine Methylation Domains.</title><p>Regions of the genome containing substantial enrichment for H3K4me3, H3K36me2 or H3K27me3 were identified using the TileMap package ##REF##16046496##[82]## and linked with genome annotation to identify substantially enriched genes. A Venn diagram indicating the number of annotated genes containing substantial enrichment for single or combinatorial modifications is shown. The shaded area represents the subset of genes enriched in both H3K4me3 and H3K27me3, which shows the most significant overrepresentation for Paf1C-regulated genes.</p></caption></fig> <fig id="pgen-1000077-g006" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000077.g006</object-id><label>Figure 6</label><caption><title>Expression Level, Size and Entropy of Paf1C-Regulated Genes.</title><p>(A) Box plots showing the distribution of expression level (left panel), expression entropy (middle panel), and gene size (right panel) for ten-percentile subsets of genes according to misregulation in <italic>vip3</italic> mutants. Distribution of genes strongly downregulated in <italic>vip3</italic> relative to wild-type is shown in column 1 of each panel; distribution for genes most strongly upregulated in <italic>vip3</italic> is shown in column 10 of each panel. Colored boxes indicate the 25<sup>th</sup>, 50<sup>th</sup>, and 75<sup>th</sup> percentiles (bottom, center line, and top of box, respectively). (B) Scatter plot relating gene expression levels with entropy for Arabidopsis genes. Genes strongly upregulated in <italic>vip3</italic> mutants are depicted as red circles, whereas strongly downregulated genes are shown as blue triangles. Locally weighted regression (Lowess) fit lines were superimposed onto the scatterplot (gray, all genes; red, upregulated; blue, downregulated).</p></caption></fig>
<supplementary-material content-type="local-data" id="pgen.1000077.s001"><label>Figure S1</label><caption><p>Representative Views of H3 and H3 Lysine Methylation Signal Profiles across Arabidopsis Chromosome IV. Histone H3 and H3 lysine methylation in the Arabidopsis genome were quantified by ChIP combined with microarray analysis (ChIP-on-chip) using antibodies directed against the carboxyl terminus of H3, H3K4me3, H3K36me2, or H3K27me3, and the Affymetrix GeneChip Arabidopsis Tiling 1.0R Array. The entire sequenced region of Chromosome IV is shown above, with numbers indicating the approximate distance (Mbp) from the end of the nonsequenced telomeric rDNA repeats. The centromere is depicted as a blue oval; HK: heterochromatic knob. Raw array data were quantile-normalized and analyzed using Affymetrix Tiling Analysis Software (TAS), and visualized using the Affymetrix Integrated Genome Browser (IGB). For each profile, signals for immunoprecipitate for H3 or each H3 modification are shown relative to the corresponding input signal, and with respect to the genomic median (horizontal colored line). Profiles in the region of the gypsy-class retrotransposon <italic>At4g06591</italic> and an active protein-coding gene, <italic>At4g27760</italic>, are depicted below to illustrate typical genic patterns.</p><p>(3.44 MB EPS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000077.s002"><label>Figure S2</label><caption><p>Representative Views of Changes in H3 and H3 Lysine Methylation Signal Profiles across Arabidopsis Chromosome IV Related to Disruption of Paf1C. Data is shown for <italic>vip3</italic> mutants with respect to wild-type plants. Profiles in the region of the gypsy-class retrotransposon, <italic>At4g06566</italic>, and the active protein-coding gene <italic>At4g190202</italic>, <italic>CHROMOMETHYLASE 2</italic>, are depicted below to illustrate typical genic pattern changes.</p><p>(3.19 MB EPS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000077.s003"><label>Figure S3</label><caption><p>Relationship Between Paf1C-Dependent, Positional H3 Occupancy, and H3 Methylation Levels and Gene Length for Strongly and Weakly Expressed Genes. Data as analyzed in ##FIG##0##Figures 1E##, ##FIG##1##2H##, and ##FIG##2##3H## were filtered to include only those genes in expression bin 10 (blue lines) or expression bin 2 (black lines).</p><p>(0.31 MB EPS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000077.s004"><label>Figure S4</label><caption><p>Positional Profiles for H3 Methylation and DNA Methylation within PcG-Regulated Genes in Arabidopsis. ChIP-on-chip was carried out using antibodies directed against H3, H3K4me3, H3K36me2, and H3K27me3. Signals were corrected for H3 and plotted across genomic regions containing <italic>AG</italic>, <italic>AGL19</italic>, <italic>PHE</italic>, and <italic>STM</italic>. A heterochromatic region containing the retrotransposon Ta3 is shown as well. Publicly available information for cytosine methylation is plotted in each frame in parallel with ChIP-on-chip results.</p><p>(2.07 MB EPS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000077.s005"><label>Figure S5</label><caption><p>H3 Modifications Related to Gene Expression Level and Entropy. (A) Protein-coding genes were assigned to ten bins according to entropy of expression (see <xref ref-type="sec" rid="s4">Materials and Methods</xref>). The mean positional signals for each modification were calculated independently for each bin. (B) Scores representing the level of each modification within protein-coding genes were calculated (see <xref ref-type="sec" rid="s4">Materials and Methods</xref>), and are illustrated by a color gradient of green (low) to black to red (high) on a plot relating expression level to entropy. Expression level and entropy are indicated on the X and Y axes, respectively. (C) Mean positional signals were calculated independently for genes expressed strongly (top 20% of expression values) in both root and aerial tissues (High, High), strongly only in aerial tissues (High, Low), strongly only in root tissues (Low, High), or expressed to a low level or silenced (lowest 30% of expression values) in both root and aerial tissues (Low, Low).</p><p>(6.09 MB EPS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000077.s006"><label>Figure S6</label><caption><p>Genic Patterns of Paf1C-Dependent H3 Occupancy and Methylations with Respect to Tissue-Specificity. Genic positional signals for H3 and H3 lysine methylations as indicated were averaged separately for protein coding genes within ten-percentile bins according to Shannon entropy for vip3 plants relative to wild-type (WT). Data are depicted across the promoter region (shown in bp from −300 to 0 relative to the presumed transcriptional start site), transcribed region (shown proportionally from 0 to 100% of total length), and 3′ region (shown in bp from 0 to +100 relative to the presumed 3′ end).</p><p>(0.32 MB EPS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000077.s007"><label>Figure S7</label><caption><p>Significance Analysis of Differences in Chromatin Signatures for Paf1C-Regulated Genes Between Wild-Type and <italic>vip3</italic> Mutants. The 95<sup>th</sup> percentile confidence intervals were calculated for differences in mean positional signals within 1,000 randomly resampled gene sets for upregulated (left) and downregulated (right) genes.</p><p>(0.35 MB EPS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000077.s008"><label>Figure S8</label><caption><p>Significance Analysis of Differences in Genic Chromatin Signatures for Paf1C-Dependent Genes Relative to Typical Genic Chromatin Signatures. The 95<sup>th</sup> percentile confidence intervals were determined for mean positional signals within 1,000 randomly resampled gene sets for upregulated (left) and downregulated (right) genes.</p><p>(0.36 MB EPS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000077.s009"><label>Figure S9</label><caption><p>H3 Methylation Profiles, Expression Level, Size, and Entropy for Genes Misregulated in the <italic>vip5</italic> or <italic>vip6</italic> Mutants. (A) Wild-type genic positional signals for H3 lysine methylations as indicated were averaged separately for those genes upregulated in <italic>vip5</italic> or <italic>vip6</italic> mutants (top row of panels) or downregulated in <italic>vip5</italic> or <italic>vip6</italic> mutants (lower panels). The 95<sup>th</sup> percentile confidence interval of signals for all genes is depicted with dotted lines. (B) Box plots show the distribution of expression level (left panel), gene size (center panel) and expression entropy (right panel) for ten-percentile subsets of genes according to misregulation in <italic>vip5</italic> or <italic>vip6</italic> mutants. Distribution of genes strongly downregulated in <italic>vip5</italic> or <italic>vip6</italic> relative to wild-type is shown in column 1 of each panel; distribution for genes most strongly upregulated in <italic>vip5</italic> is shown in column 10 of each panel. Colored boxes indicate the 25<sup>th</sup>, 50<sup>th</sup>, and 75<sup>th</sup> percentiles (bottom, center line, and top of box, respectively). (C) Scatter plot relating gene expression levels with entropy for Arabidopsis genes. Genes strongly upregulated in <italic>vip5</italic> or <italic>vip6</italic> mutants are depicted as red circles, whereas strongly downregulated genes are shown as blue triangles. Lowess fit lines were superimposed onto the scatterplot (gray, all genes; red, upregulated; blue, downregulated).</p><p>(1.71 MB EPS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000077.s010"><label>Figure S10</label><caption><p>Clustering and Analyses of Genes According to Chromatin Modification Profile Class. (A) Cluster analyses were performed for the ∼18,000-gene set based on genic positional signals for H3K4me3, H3K36me2, H3K27me3, H3, and DNA methylation. Data was plotted across promoter regions (columns 1–3 in each modification panel), TSS (column 4), transcribed regions (columns 5–14) and 3′ end (column 15). (B) Averaged positional profiles for H3, H3 modifications and DNA methylation ##REF##16949657##[84]## are shown separately for each of the four groups. Positions are relative to a representative transcriptional unit shown at bottom. (C) Box plots showing the percentile level of expression (top), expression entropy (middle) and gene length (bottom) for each group. Boxes indicate the 25<sup>th</sup>, 50<sup>th</sup>, and 75<sup>th</sup> percentiles (bottom, center line, and top of box, respectively).</p><p>(1.94 MB EPS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000077.s011"><label>Figure S11</label><caption><p>Reproducibility and Confirmation of ChIP-on-Chip Data. (A) An M versus A (MvA) plot representing signal intensities from the two biological replicates of each genotype is shown. The X-axis is defined as the average of the log base 2 of the intensities from the two replicates, and the Y-axis is the difference of the log base 2 of the intensities from the two replicates. The color bar at right indicates the number of probes on the plots. (B) ChIP analysis of H3 and H3 modifications within the FLC locus is shown. ChIP was carried out using antibodies recognizing H3K4me3 (upper left), H3K36me2 (upper right), H3K27me3 (lower left), and the H3 carboxyl terminus (lower right) within a promoter segment (red), 5′ region (yellow), or 3′ region (green). A 5′ and 3′ region of the ACTIN7 gene was used for an internal control. Band intensities from gel images were quantified and normalized based on those for ACTIN7 (lower band in each gel image). ChIP analysis was performed twice using biologically independent samples and yielded essentially identical results.</p><p>(3.96 MB EPS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000077.s012"><label>Table S1</label><caption><p>Representation of <italic>VIP3</italic>-Dependent Genes within Chromatin Enrichment Groups Depicted in ##FIG##4##Figure 5##.</p><p>(0.03 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000077.s013"><label>Table S2</label><caption><p>Representation of <italic>VIP3</italic>-Dependent Genes within Chromatin Enrichment/Profile Groups Depicted in ##SUPPL##9##Figure S10##.</p><p>(0.03 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000077.s014"><label>Table S3</label><caption><p>Reproducibility of Tiling Microarray Results (Correlation Coefficient, R).</p><p>(0.03 MB DOC)</p></caption></supplementary-material>
<fn-group><fn fn-type="COI-statement"><p>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p>This work was supported by funding from the National Science Foundation (Grant MCB 0445867) to SVN.</p></fn></fn-group>
<graphic xlink:href="pgen.1000077.g001"/> <graphic xlink:href="pgen.1000077.g002"/> <graphic xlink:href="pgen.1000077.g003"/> <graphic xlink:href="pgen.1000077.g004"/> <graphic xlink:href="pgen.1000077.g005"/> <graphic xlink:href="pgen.1000077.g006"/>
<media xlink:href="pgen.1000077.s001.eps"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000077.s002.eps"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000077.s003.eps"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000077.s004.eps"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000077.s005.eps"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000077.s006.eps"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000077.s007.eps"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000077.s008.eps"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000077.s009.eps"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000077.s010.eps"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000077.s011.eps"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000077.s012.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000077.s013.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000077.s014.doc"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["85"], "element-citation": ["\n"], "surname": ["Smyth"], "given-names": ["GK"], "year": ["2004"], "article-title": ["Linear models and empirical Bayes methods for assessing differential expression in microarray experiments."], "comment": ["Statistical Applications in Genetics and Molecular Biology 3, No. 1, Article 3"]}]
PLoS Genet. 2008 Aug 22; 4(8):e1000077