Recent epigenomic studies point out that epigenetic control of transcription elongation is a widespread regulatory mechanism. Intragenic DNA methylation occurs at higher density 123 and has a larger effect on expression level than promoter methylation 4, inhibiting transcription elongation in filamentous fungi 56, plant protoplasts 7, Arabidopsis 2, and mammalian cells 8. Remarkably, methylation of protein coding regions alone can inhibit gene expression and its inhibitory effects are larger when it occurs nearer the start codon 7. Intriguingly, the methylation map of rice chromosomes exhibits single peaks near start codons 4. It is probable that the methylation peak near the start codon exerts major inhibitory effects on transcription elongation.

It may be that specific organization of nucleosomes surrounding the coding start region is also important in regulating transcription elongation. The maps of H2A.Z-containing nucleosomes in yeast and fly reveal conserved positioning of a nucleosome downstream of the transcription start site (TSS; the +1 nucleosome) 910. Notably, another yeast study provides a hint that the +1 nucleosome peaks very close to the start codon 11. In addition, there seems to be a nucleosomal peak prior to the 3' end of the yeast open reading frame (ORF) as well 12. The +1 nucleosome of fly is located further from the TSS than that of yeast 10. However, fly nucleosome patterns surrounding the start and stop codons have never been examined.

Taken together, anecdotal findings in different species suggest a role for coding region boundaries in maintaining nucleosome deposition and DNA methylation. From a mechanistic perspective, the epigenetic marks can act as 'roadblocks' to RNA polymerase II (Pol II) progression. However, previous studies have focused on transcription initiation or termination, but not elongation. In this study, we attempt to provide insight into elongation inhibition surrounding translation start and end sites, which has never been observed. Given the impact of DNA methylation on nucleosome formation 13, correlations or differences between the two epigenetic marks also deserve systematic investigation.

We first surveyed published nucleosome positions in yeast and fly 910. When aligned at the TSS, there is a significant difference in +1 nucleosome position between the two species 10 (Additional data file 1a). However, aligning by the coding region places the first coding nucleosome in similar positions in the two species - that is, just downstream of the start codon (Additional data file 1b). We also identified a highly conserved nucleosome immediately upstream of the 3' coding end in both species (Additional data file 1c).

Analyzing the H2A.Z map for human T cells 14 also revealed nucleosomal peaks just downstream of start codons and just upstream of stop codons, marking both ends of the coding sequences (Figure 1a). Meanwhile, boundaries at both ends of transcripts are tightly coupled with nucleosome-free regions, potentially allowing access of the initiation and termination complex (Figure 1b). The nucleosome-free region at the TSS was followed by the +1 nucleosome. However, the association of the +1 nucleosome and the TSS appears to be weaker than that of the first coding nucleosome and the start codon. The patterns of nucleosome positioning for some individual genes are shown in Additional data file 2.

We carried out Solexa sequencing of methylated DNA from human T cells and found methylation peaks at the exact same positions (Figure 1c). We also profiled the mouse liver and found the same patterns (Additional data file 3). Together, the nucleosomal peaks are observed in human, fly, and yeast, and the methylation peaks in human, mouse, and plants.

To examine their role in regulating transcription, we first related the level of the epigenetic peaks to expression level. We found that highly expressed genes are depleted of the epigenetic peaks (Figure 2a), consistent with findings of a nucleosomal barrier against high transcription rate 15. However, the overall correlation was not strong. We then estimated elongation efficiency as mRNA production per unit density of elongating Pol II. Upon initiation, Pol II is phosphorylated at Ser5 in its carboxy-terminal domain, switching to an elongation-competent form. Thus, we calculated the ratio of expression level to the density of Ser5-phosphorylated Pol II within the transcript body. Genes with high elongation efficiency will show high expression levels even with a low density of elongating Pol II across the transcribed region, and the opposite for low elongation efficiency. A strong association was found between the level of the epigenetic peaks and elongation efficiency (Figure 2b; Additional data file 4).

Without any interference, elongating Pol II should be distributed evenly across the transcript body except at the initiating and terminating sites. Faced with roadblocks, however, Pol II pauses and a pileup of Pol II forms, which can be observed as a peak of Pol II density. Thus, to demonstrate Pol II pausing at epigenetic marks, three criteria should be met: the presence of a Pol II peak; the presence of an epigenetic peak; and a correspondence between the positions of the two peaks.

Elongating Pol II appears to pile up immediately upstream of the nucleosomal peaks at both ends of protein coding units (Figure 2c), satisfying the three criteria. However, there seem to be confounding effects from Pol II enriched at nearby transcription initiation and termination sites. For example, the Pol II tail downstream of the stop codon (arrow above the right panel of Figure 2c) might reflect Pol II awaiting to be released from the transcription termination site. We thus selected genes with a long (> 5 kb) 5' untranslated region (UTR) and examined Pol II density around the TSS and start codon separately. Both unphosphorylated and elongating Pol II were enriched at the TSS, but only elongating Pol II showed high downstream density (Additional data file 5) with a pileup upstream of the start codon (left panel of Figure 2d). Another Pol II peak upstream of the start codon (arrow above the left panel of Figure 2d) seems to reflect Pol II at the TSS. A pileup of Pol II was also found before a long (> 5 kb) 3' UTR (right panel of Figure 2d), indicative of Pol II blockage that occurs independently of the transcription termination site.

Pol II pausing was not observed with low nucleosome occupancy (Figure 2e), indicating that elongating Pol II is indeed impeded by the boundary nucleosome. To observe the specific effect of the boundary nucleosome, we computed relative occupancy at the coding end compared to the surrounding region. Higher or lower nucleosome occupancy near the coding end directly led to higher or lower Pol II density in the immediate upstream region (Additional data file 6).

We roughly estimated the percentage of genes that are affected by Pol II pausing by comparing the average Pol II density around boundaries and that across surrounding regions. We found that 54% of genes exhibit higher Pol II density near the start codon than in the flanking region and 41% of genes have a Pol II peak near the stop codon.

Nucleosome positioning is governed by DNA sequences 1116. Methylation level is dependent on the CpG content of the target sequence 17. Given the distinctive patterns of nucleosome positioning and methylation maintained in specific regions among different species, there should be strong constraints on the underlying DNA sequences. Being under strong natural selection, protein coding sequences could be better candidates than UTRs for conserved epigenetic targets downstream of transcription initiation and upstream of termination. Coding region boundaries might be subject to considerable negative selection that purifies sequence changes that are detrimental to nucleosome deposition or DNA methylation.

We examined two sequence characteristics deemed to be involved in epigenetic programming: DNA-bending propensity and CpG density. DNA-bending propensity, the ability of nucleotide sequences to wrap around a histone complex, is an important determinant of nucleosome formation 1819. DNase I digestion experiments indicate that bending parameters for the start codon and three stop codons are among the 8 highest out of those for the 32 trinucleotides 20. Therefore, they can significantly contribute to the high bendability of coding boundaries (Additional data file 7). The boundary sequences with higher bendability tend to be more enriched for nucleosomes (upper panel in Figure 3a). Unexpectedly, DNA methylation level was also proportional to bending propensity.

CpG density should be a determinant of methylation level. The boundary sequences with intermediate CpG density were densely methylated (lower panel in Figure 3a). In contrast, nucleosome occupancy was dominant among genes with lower CpG density. While the two marks commonly have affinity for base compositions with high bending propensity, DNA methylation at CpG sites might affect structural DNA bending and nucleosome formation. The proportions of genes marked by both nucleosomes and methylation or by just one of these are shown in Figure 3b. More genes are specifically marked by nucleosomes than methylation, possibly because many boundary regions have relatively low CpG density (gray curve in lower panels of Figure 3a).

A group of genes had highest CpG density at the 5' end (gray curve at CpG density > 0.8 in bottom left panel of Figure 3a). These genes showed a markedly reduced level of DNA methylation, reflecting the fact that CpG islands are typically unmethylated. Indeed, 97.2% of these genes contained a CpG island within their promoter (-1,000 bp to 500 bp from the TSS) and 92.8% had a short (< 500 bp) UTR (P < 10^-100), an indication that promoter CpG islands are overlapping or located very close to the start codon. These genes exhibited high expression compared to the rest of the genes (P < 10^-80), even higher than the genes with a promoter CpG island (P = 1.4 × 10^-10) (Additional data file 8), indicating additional effects of elongation control.

Next, we explored the intragenic distribution of the marks. Although significantly higher than its flanking region, the 5' peak is generally lower than the 3' peak (Additional data file 9a). Meanwhile, k-means clustering shows that most genes have higher peaks at both ends compared to the central region (Additional data file 9b). We then examined these patterns according to the size of the coding region (Figure 4a). We found that genes with a short coding sequence (< 1 kb) have nucleosomes and methylation in their inner region over a large portion of the gene body. In particular, their 3' ends lack both marks, in sharp contrast to most other genes, in which the marks are shifted toward both ends with a bias toward the 3' end. Unlike at the 5' end, both marks commonly peaked at the 3' end, especially in many genes of intermediate size (Figure 4b).

Notably, nucleosome composition within the coding region of yeast or fly genes is not sharply shifted to the coding boundaries - the 3' peak is especially not very prominent (upper panel in Figure 4c) - a similar pattern to that seen for human genes of similar size (1 to 2 kb). Intragenic DNA methylation in plant genomes is also concentrated in the central region of protein coding sequences (lower panel in Figure 4c). Arabidopsis genes share a similar pattern with human genes of similar size (1 to 3 kb). Rice genes are longer than Arabidopsis genes (Figure 4d) and have detectable, if not complete, peak patterns, which can explain the observed 5'-end peak 4.

Evolutionary processes seem to have maintained the boundaries of protein coding sequences as targets for nucleosome binding or DNA methylation. Incorporating novel protein domains or introns in inner sequences might have concentrated the epigenetic marks in the 5'- and 3'-end segments. It could be a good strategy to diversify DNA sequences in the middle of coding regions to encode various protein functions while constraining the boundary sequences, including the start and stop codons, to create an epigenetically favorable environment.

What is the physiological role of the conserved 5'- and 3'-end peaks? A few studies have hinted at their role in RNA splicing. First, a crosstalk between nucleosomes and RNA splicing has been suggested 2122. The chromatin remodeling complex SWI/SNF was shown to create roadblocks of nucleosomes that slow Pol II elongation, which facilitates the inclusion of the exon during RNA splicing. Another study showed that H3K46me3, which impedes Pol II elongation, was specifically enriched in constitutive exons compared to alternative exons, thus suggesting it has a role in controlling RNA splicing 2324.

These findings suggest a novel role for nucleosome deposition and DNA methylation at coding boundaries in controlling RNA splicing. The first and last coding exons should be constitutively included in a mature transcript since skipping these exons causes translation failure. Slowing elongation by these epigenetic roadblocks might facilitate the inclusion of these indispensable boundary exons. DNA sequences surrounding the start and stop codons are highly conserved to ensure a favorable epigenetic environment, explaining why the epigenetic peaks are found in the specific loci of the first and last exons. This mechanism seems to be more common among higher organisms with longer genes.

Previous epigenetic studies have focused on promoter regions in an attempt to associate epigenetic patterns with transcription initiation activity. However, recent genome-wide epigenetic studies challenge the traditional viewpoint and have shed new light on the epigenetic control of transcription elongation. Notably, epigenetic inhibition of Pol II elongation has been proposed to facilitate the inclusion of constitutive exons during RNA splicing. The evolutionarily conserved epigenetic patterns we identified here seem to ensure the inclusion of the first and last coding exons as they are indispensable for accurate translation. Further mechanistic studies are required to support this hypothesis.

GEO: Gene Expression Omnibus; GST: glutathione S-transferase; MIRA: methylated CpG island recovery assay; ORF: open reading frame; Pol II: RNA polymerase II; TSS: transcription start site; UTR: untranslated region.

The authors declare that they have no competing interests.

JKC conceived of the study, carried out the analysis, and wrote the manuscript. JBB carried out the purification and sequencing of methylated DNA. JL processed the raw sequencing data and participated in data analysis. TYK participated in sample preparation. YJK participated in study design and coordination, and finalized the manuscript.

The following additional data are available with the online version of this paper: a figure showing nucleosome patterns surrounding the TSS, start codon, and stop codon in yeast and fly (Additional data file 1); a figure showing illustrative genes with nucleosomal peaks at coding boundaries (Additional data file 2); a figure showing DNA methylation level surrounding the transcript and coding region boundaries in the mouse liver (Additional data file 3); a figure showing nucleosome occupancy according to differential Pol II elongation efficiency (Additional data file 4); a figure comparing densities of Ser5-phosphorylated and unphosphorylated Pol II (Additional data file 5); a figure showing Pol II density with higher and lower nucleosome occupancy (Additional data file 6); a figure showing DNA bending propensity at the start and stop codons (Additional data file 7); a figure demonstrating the length of the 5' UTR and gene expression level for genes with high CpG density around the start codon (Additional data file 8); a figure showing the overall patterns of nucleosome occupancy and DNA methylation level inside the protein coding region (Additional data file 9).