Evolution of lysine acetylation in the RNA polymerase II C-terminal domain
- Corinne N Simonti†1,
- Katherine S Pollard†2, 3,
- Sebastian Schröder†2,
- Daniel He2,
- Benoit G Bruneau2,
- Melanie Ott2 and
- John A Capra1, 4Email author
© Simonti et al.; licensee BioMed Central. 2015
Received: 24 November 2014
Accepted: 24 February 2015
Published: 10 March 2015
RPB1, the largest subunit of RNA polymerase II, contains a highly modifiable C-terminal domain (CTD) that consists of variations of a consensus heptad repeat sequence (Y1S2P3T4S5P6S7). The consensus CTD repeat motif and tandem organization represent the ancestral state of eukaryotic RPB1, but across eukaryotes CTDs show considerable diversity in repeat organization and sequence content. These differences may reflect lineage-specific CTD functions mediated by protein interactions. Mammalian CTDs contain eight non-consensus repeats with a lysine in the seventh position (K7). Posttranslational acetylation of these sites was recently shown to be required for proper polymerase pausing and regulation of two growth factor-regulated genes.
To investigate the origins and function of RPB1 CTD acetylation (acRPB1), we computationally reconstructed the evolution of the CTD repeat sequence across eukaryotes and analyzed the evolution and function of genes dysregulated when acRPB1 is disrupted. Modeling the evolutionary dynamics of CTD repeat count and sequence content across diverse eukaryotes revealed an expansion of the CTD in the ancestors of Metazoa. The new CTD repeats introduced the potential for acRPB1 due to the appearance of distal repeats with lysine at position seven. This was followed by a further increase in the number of lysine-containing repeats in developmentally complex clades like Deuterostomia. Mouse genes enriched for acRPB1 occupancy at their promoters and genes with significant expression changes when acRPB1 is disrupted are enriched for several functions, such as growth factor response, gene regulation, cellular adhesion, and vascular development. Genes occupied and regulated by acRPB1 show significant enrichment for evolutionary origins in the early history of eukaryotes through early vertebrates.
Our combined functional and evolutionary analyses show that RPB1 CTD acetylation was possible in the early history of animals, and that the K7 content of the CTD expanded in specific developmentally complex metazoan lineages. The functional analysis of genes regulated by acRPB1 highlight functions involved in the origin of and diversification of complex Metazoa. This suggests that acRPB1 may have played a role in the success of animals.
KeywordsRNA pol II C-terminal domain Heptad repeats RNA pol II acetylation Posttranslational modification Metazoa Phylogenetics
Eukaryotic RNA polymerases evolved from a single ancestral enzyme into three structurally related RNA polymerase enzymes (I–III) with specialized functions in eukaryotes. RNA polymerase II generates all protein-coding mRNAs as well as a large number of non-coding microRNAs (miRNAs), small nuclear RNAs (snRNAs) and small nucleolar RNAs (snoRNAs). The enzyme is composed of 12 subunits (RPB1–12), five of which are shared among the three eukaryotic polymerase complexes. The largest subunit, called RPB1, is unique to RNA polymerase II and is involved in its catalytic activity.
The C-terminal domain (CTD) of RPB1 is essential for the proper regulation of RNA polymerase II . The CTD consists of a largely unstructured, repetitive stretch of tandem heptad amino acid repeats with a consensus sequence of tyrosine-serine-proline-threonine-serine-proline-serine (Y1S2P3T4S5P6S7). The CTD is connected to the core polymerase enzyme via a flexible linker close to the RNA exit site  and contains a protein-protein interaction surface for cofactors involved in the regulation of transcription initiation, elongation and RNA processing—highly specialized functions that determine the speed and reliability of the polymerase enzyme traversing a gene during transcription [3,4].
In mammals, the CTD heptad repeats undergo a sequence of characteristic posttranslational modifications during transcription  and these modifications determine which cofactors can bind . Phosphorylation of the S5 residue occurs when the polymerase binds to gene promoters, and transcription is initiated . The phospho-S2 modification is a hallmark of elongating polymerase complexes and is maintained until transcription termination occurs . Phosphorylation of S7 residues in linker-proximal heptad repeats is important for the transcription of snRNAs and facilitates the interaction of the RNA polymerase II complex with the Integrator complex at these genes [9,10]. Methylation of a single arginine residue plays a role in the production of snoRNAs as well as snRNAs and serves as a binding site for the tudor domain–containing protein TDRD3 . Phosphorylation of T4 residues was shown to play a novel function in 3′-end processing of histone pre-mRNAs . Phosphorylation of Y1 was coupled to antitermination of transcription , and was recently associated with promoters, enhancers, and degradation of RNA polymerase II [14,15]. However, in spite of our increasing knowledge of these marks, the full number and identity of repeats that undergo modification during transcription are unknown.
The CTD repeat motif and tandem repeat orientation are observed in most eukaryotic model organisms, and this configuration evolved early in the history of eukaryotes . However, there is considerable diversity in the number and sequence of the repeats across eukaryotic clades, and the complexity of the CTD is roughly correlated with developmental complexity in animal, plant, and fungal multicellular lineages . It has been proposed that this diversity in CTD sequences reflects functional constraint due to lineage-specific CTD–protein interactions . In mammals, non-consensus repeats are required for the stability of the CTD . However, deletion of most of the non-consensus repeats does not affect the housekeeping functions of mammalian cells, but it does impair activator-induced gene transcription [18,19].
A novel CTD modification of non-consensus repeats necessary for the proper regulation of polymerase pausing was recently discovered . Lysine residues at the seventh position (K7) in non-consensus distal RPB1 CTD repeats are acetylated (acRPB1) in mouse and human, but not in yeast (Saccharomyces cerevisiae), where the repeats lack lysines. This modification, performed by the acetyltransferase (KAT) enzyme p300/KAT3B, which is absent in yeast, distinctly marks promoter-proximally paused polymerases. Disruption of this mechanism interfered with the expression of two growth-factor-induced genes regulated by polymerase pausing, but did not influence expression or polymerase occupancy at two non-paused genes .
Knowledge of the evolutionary origins of a gene or pathway can inform analysis of its functions [21,22]. Thus, we undertook an evolutionary and functional analysis of the RPB1 CTD sequence and the genes influenced by acRPB1. We found that the presence of multiple K7 residues, and thus the potential for RPB1 CTD acetylation, arose with animal multicellularity during an expansion in the overall number of CTD repeats in Metazoa. Our phylogenetic analysis further showed that p300/KAT3B, the acetyltransferase that modifies the RPB1 CTD, was present at the appearance K7-containing repeats. We then performed a genome-wide survey of acRPB1 occupancy and its influence on gene regulation in mice. Genes with enrichment for acRPB1 at their promoters and genes dysregulated when acRPB1 was disrupted were specifically enriched for functions in growth-factor signaling, cell adhesion, vascular development, and cell-cell interaction. In addition, the two sets of acRPB1 sensitive genes were enriched for evolutionary origins in early eukaryotes through the ancestor of Euteleostomi (bony vertebrates). Together, our functional and evolutionary results support a model in which K7-containing CTDs were selected for in the early history of animal multicellularity. Given the association between polymerase pausing and acRPB1, the potential to acetylate these residues may have enabled tighter control of gene expression as animals grew in complexity and diversified. Indeed, acRPB1 now influences the regulation of growth factor-target genes and genes involved in lineage-specific processes, such as cell adhesion and vasculature development, in mammals.
The RPB1 CTD experienced an expansion in repeat number and lysine content in the early history of animals
Ancestral state reconstruction using symmetric Wagner parsimony on the species phylogeny revealed an expansion of CTD repeats in the common ancestor of all Metazoa (Figure 2), as expected from previous studies [16,17]. We estimate that the ancestor of all Metazoa had 44 repeats, while the last common ancestor of Metazoa and their closest relatives, the choanoflagellates, had only 31 repeats. K7-containing repeats also first appeared consistently in the ancestor of Metazoa, and this was followed by an increase in the number of K7-containing repeats (from 3 to 7) in the last common ancestor of Deuterosomia. All the Cnidaria and Ecdysozoa examined have K7-containing repeats, but with the exception of the deer tick (Ixodes scapularis) and hydra (Hydra magnipapillata), the number is comparatively small (see Discussion). The plants, algae, and other eukaryotes we examined have half as many or fewer repeats as the Metazoa, and only one has more than two lysine-containing repeats: Plasmodium falciparum, a human pathogen. (See the Discussion and  for more on the evolution of the Plasmodium falciparum CTD.) Thus, we conclude that CTD repeat length markedly increased with the origin of animals, and the distal repeats gained lysine residues and expanded further in different animal lineages. These results are in agreement with an extensive study of CTD repeat number and content that appeared while this manuscript was in preparation .
The acetyltransferase p300 was present for the metazoan increase in CTD repeat count and lysine content
The p300/KAT3B acetyltransferase was shown to efficiently target and acetylate the mammalian RPB1 CTD, while another KAT enzyme, PCAF/KAT2B, did not . Recent analysis of the Capsaspora owczarzaki genome revealed that p300/KAT3B and its paralog, CBP/KAT3A, were present before the divergence of Capsaspora from choanoflagellates and metazoans . To further refine this estimate, we searched for homologs in the additional early branching holozoan species considered in our CTD analysis and other eukaryotic genomes. All tested holozoans contain a gene with significant homology (all have BLAST E-values < 7E-155) to human/mouse p300 and CBP; thus, the ancestor of this gene was present for the metazoan expansion of the CTD and introduction of K7-containing repeats. Homologs of this family are present in plants , but we did not find clear homologs in the fungi. This suggests that the KAT3 family may be a more ancient eukaryotic family that has been lost in some lineages and maintained in others.
We also determined the phylogenetic age of 17 KAT enzymes with expressed in both mouse and humans  using ProteinHistorian , a web server that estimates protein ages, and found a wide range of predicted evolutionary ages (Additional file 3).
Acetylation of RPB1 regulates many mammalian genes
Mutating the RPB1 CTD to prevent acetylation was shown to disrupt the expression of two immediate-early genes (c-Fos and Egr2) in the epidermal growth factor (EGF) and ERK/MAPK pathways with paused polymerases, but it did not influence the expression of several housekeeping genes . To explore the functional relevance of acetylation of the RPB1 CTD genome-wide, we analyzed two sets of genes influenced, directly and indirectly, by acRPB1. First, we analyzed the total and acetylated RNA Pol II promoter occupancy in mouse embryonic stem cell (ESC) ChIP-Seq data from Schroeder et al. . We found 10,078 genes with acetylated RNA Pol II occupancy at their promoters, and 707 genes with significant enrichment for acetylated RNA Pol II at their promoters over total RNA Pol II levels (Methods; Additional file 4). We will refer to these gene sets as “acRPB1 occupied” and “acRPB1 enriched,” respectively. Second, we expressed a mutated murine HA-tagged RPB1 in which all K7 residues were substituted by arginines (8KR) in mouse NIH/3 T3 fibroblasts (Methods; ). This mutation prevents acetylation while preserving the positive charge at these positions. Gene expression profiling using microarrays identified 1787 RNAs that were significantly (False Discovery Rate (FDR) < 0.05) up- or down-regulated in 8KR cells compared to wild type cells. We refer to these as “acRPB1 dysregulated genes”. In the following, we focus on the acRPB1 enriched genes, which reflect the direct influence of acRPB1, and the acRPB1 dysregulated genes, which summarize the broader downstream effects of RPB1 acetylation.
There was significant overlap between the acRPB1 enriched and dysregulated genes (83 genes; p < 0.0001, chi-squared test). The relatively modest magnitude of this overlap likely reflects the different origins of the gene sets. The ChIP-Seq was performed on ESCs, while the microarrays are from fibroblasts. AcRPB1 is likely involved in many context-dependent responses in different cells, and a gene with enrichment for acRPB1 at its promoter in one context is not necessarily dysregulated in all contexts if the ability to acetylate RPB1 is lost. It also suggests that some genes dysregulated in the disruption of acRPB1 may not be directly regulated by acRPB1. However, it is also possible that some differences in the gene sets result from technical artifacts due to the different experimental techniques used to define them.
Genes influenced by acetylation of RPB1 are enriched for functions in growth-factor response, cell adhesion, regulation of gene expression, and vasculature development
Functional annotations and pathways enriched among acRPB1 enriched genes
Number of genes
Protein amino acid phosphorylation (GO BP)
Neg. regulation of gene expression (GO BP)
Neg. regulation of macromol. biosynthetic process (GO BP)
Functional annotations and pathways enriched among acRPB1 dysregulated genes
Number of genes
Cell adhesion (GO BP)
Vasculature development (GO BP)
Blood vessel development (GO BP)
Focal adhesion (KEGG)
ECM-receptor interaction (KEGG)
Blood vessel morphogenesis (GO BP)
The disruption of acRPB1 was previously shown to abolish expression of two epidermal growth-factor-induced genes, c-Fos and Egr2, in mouse fibroblasts . In order to further explore the connection between EGF-induced genes and RPB1 CTD acetylation, we analyzed the expression and promoter acRPB1 occupancy of several curated sets of growth factor responsive genes. First, we identified 49 mouse homologs of human genes whose transcription is directly induced by EGF signaling (Additional file 6) . Overall, EGF-induced genes overlap significantly with acRPB1 enriched genes (5 of 48 with promoter occupancy data, p = 0.001, Yates’ chi-squared test), as well as the dysregulated genes (10 of 49; p = 0.001). The association between acRPB1 and EGF became even stronger when we considered promoter acRPB1 occupancy (regardless of whether it was significantly higher than total RPB1 occupancy): 73% (35 of 48) of EGF genes had acRPB1 occupancy (p = 1E-8). Similarly, the number of EGF genes nominally dysregulated (p < 0.05) was also significant (19 of 49; p = 0.003).
Given the association between acRPB1 and EGF-induced genes, we investigated whether acRPB1 regulated genes induced by other growth factors by analyzing 100 genes induced by platelet-derived growth factor (PDGF; Additional file 6) . There was a similar enrichment for PDGF-induced genes among acRPB1 occupied (65 of 100; p = 0) and dysregulated (23 of 100; p = 8E-8) genes, but not among acRPB1 enriched genes (2 of 100; p = 0.906). Thus, genes induced by both growth factors are regulated by acetylation of the RPB1 CTD, though perhaps mostly indirectly in the case of PDGF. Supporting this conclusion, EGF-induced genes were more strongly enriched among acRPB1-enriched genes than PDGF-induced genes (p = 0.008, Fisher’s exact test), but both growth factors were similarly represented among dysregulated genes (p = 0.381). Fifteen genes were present in both the EGF- and PDGF-induced sets; the above results were similar when these genes were not considered in the analysis.
Growth factor-induced genes can be split into two categories: immediate-early genes (IEGs) and delayed primary response genes (DPRGs), which are transcribed after the IEGs, but before secondary response genes . These two classes are relevant, because their transcription occurs before de novo protein synthesis, and thus they reflect the direct effects of growth factor signaling versus those further downstream. Regardless of the growth factor, both IEGs and DPRGs were enriched among acRPB1-occupied and -dysregulated genes, but the effect was stronger for IEGs than DPRGs (Fisher’s exact test p = 0.047) supporting the hypothesis that acRPB1 targets many IEGs.
Many genes regulated by acetylation of RPB1 originated early in the evolution of eukaryotes and animals
AcRPB1 likely first occurred early in the development of animals
Our integrated evolutionary and functional analyses of RPB1, its modifiers, regulators, and targets suggest that the potential for acRPB1 was present in early animals. While we cannot definitively establish the timing of the first acRPB1 event, several lines of evidence suggest that acRPB1 began to occur soon after K7-containing heptad repeats appeared during the increase in RPB1 CTD repeat length in the metazoan ancestor. First, elevated repeat count and K7-containing repeats are present in nearly all Metazoa. Second, K7-containing repeats are rare outside of Metazoa. Third, the ancestor of the p300 KAT enzyme, which acetylates CTD K7 residues, was present when the CTD expanded. Fourth, the genes bound and regulated by acRPB1 are enriched for evolutionary origins before and during the appearance and diversification of animals and depleted for origins more recent than the last common ancestor of all bony vertebrates (Euteleostomi). Finally, the functions enriched among genes regulated by acRPB1 suggest involvement in traits essential to the origin and development of complex animal multicellularity.
AcRPB1 may have served as a basis for establishing the regulation of some functions essential to complex multicellular animals
Transitions to multicellular lifestyles require organisms to perform a range of new functions that facilitate cell adhesion, cell-cell signalling, and context-dependent gene expression [31,32]. We found enrichment for functions in cell adhesion, extracellular matrix signaling, gene regulation, and phosphorylation among genes influenced by acRPB1 (Tables 1 and 2). This suggests a connection between new regulatory mechanisms, such as polymerase pausing, enabled by acRPB1 and the development of animal multicellularity. In addition, the enrichment for vascular development among the genes sensitive to acRPB1 supports the involvement of acRPB1 in the later development and regulation of animal body plan complexity. Indeed, intercellular regulation of proliferation, differentiation, and migration by paracrine factors such as EGF is a unique feature of higher eukaryotes and lies at the root of body complexity in animals . Our analyses show that proper regulation of genes induced by two growth factors, EGF and PDGF, is dependent on the presence of K7 residues in the CTD. While our study explicitly provides evidence for the importance of K7 residues in EGF- and PDGF-induced transcription, it also suggests that other signal transduction pathways in higher eukaryotes may rely on K7 residues for proper gene activation.
What is the role of acRPB1 in different animal lineages?
The eukaryotic RPB1 CTD experiences posttranslational phosphorylation and methylation modifications that are essential for proper regulation of its activity. The initial expansion in CTD repeat number during the origin of Metazoa likely increased the potential for these RPB1 modifications. Our work suggests that acRPB1 appeared soon after this initial expansion. The integration of K7 residues into the RPB1 heptad repeat appears to predate the divergence of the Cnidaria, Ecdysozoa, and Deuterostomia; however, these clades exhibit differences in the number and type of repeats. Deuterostomes exhibit dramatic and consistent high K7 count (Figure 2). Nonetheless, the presence of at least one K7-containing repeat is conserved in all species analyzed in these clades, and some species (H. vulgaris and I. scapularis) have high K7 counts. Given the role of acRPB1 in promoter-proximal polymerase pausing  and the importance of pausing in non-deuterostome species, such as Drosophila melanogaster , it is possible that acRPB1 occurs commonly in these species in spite of their lower K7 content. Indeed, a single K7-containing repeat in a modified mouse CTD was sufficient to produce detectable acetylation in vitro . Most pausing research has been performed in human, mouse, or fly, and the full phylogenetic range of polymerase pausing is not known. However, it appears that pausing is rare in Caenorhabditis elegans , a species with only one K7 repeat, but systematic studies of pausing and acRPB1 in a more diverse array of species are necessary to resolve this question.
AcRPB1 may be associated with developmental complexity in different eukaryotic lineages
Recent work has established the plasticity of the RPB1 CTD across eukaryotes and suggested that lineage-specific CTD modifications are associated with increased complexity in different multicellular lineages . While our work does not prove the involvement of acRPB1 in the development of animal multicellularity, it suggests that this novel regulatory mechanism may have played a role in the establishment of the complex gene regulatory programs necessary for the proper cooperation of different cells within complex multicellular animals, in addition to known expansions and rearrangements of gene families .
The presence of a non-consensus K7 repeat in many plant species  suggests the possibility that acRPB1 could have evolved in parallel in another developmentally complex multicellular lineage. However, it is not known whether plant K7 repeats are acetylated.
Why are many K7-containing repeats seen in some species of malaria?
K7-containing repeats are rare outside of animals, yet in agreement with previous studies, we see many K7 repeats in P. falciparum, a parasite that causes malaria in humans (Figure 2). In general, Plasmodium species have very short repeat regions with high variability between and within species [16,23]. Primate-infecting Plasmodium species have many K7 repeats, yet this is not true of those that infect other mammals, like rodents, and birds. Analysis of the CTDs of many Plasmodium species suggests that expansion of K7 repeats in the CTD occurred twice in parallel in different lineages of primate parasites; however, the forces driving these independent primate-specific expansions are not known. The establishment of the acetylation of the CTD  and its importance in regulating many ancient genes suggests that CTD acetylation may be involved. However, since acRPB1 occurs in mice, the presence of K7 repeats and acRPB1 in the host does not necessarily result in high K7 content in the Plasmodium CTD. Horizontal transfer of epigenetic regulators, like the Set2 and Set8 methyltransferase domains, occurred in the ancestor of Apicomplexans and has been linked to their transition to parasitism . The Set2 domain interacts with the CTD, but it is absent in rodent infecting Plasmodium species. This suggests that there may be differences in the epigenetic methylation and acetylation of the CTD in Plasmodium species infecting different mammals, but the nature and effects of these modifications remain uncertain.
In this work, we trace the evolutionary origin of K7-containing RPB1 CTD repeats to the early history of animals. Our integrated evolutionary and functional analyses suggest that the potential for acRPB1 was present at this time and suggest, due to the conserved increased CTD repeat count and K7 content, that acRPB1 is a common regulatory mechanism in many animals. Most studies of the genetic changes involved in transitions to animal multicellularity have implicated expansions or rearrangements of gene families . While our work does not prove the involvement of acRPB1 in this transition, it suggests that this novel regulatory mechanism may have provided a foundation on which gene regulatory programs involved in the proper function and cooperation of different cells within complex multicellular animals could be built.
We identified species with sequenced RPB1 genes using BLAST searches and previous studies of RPB1 . We confirmed that all sequenced vertebrates have the same pattern of repeats as human and selected human, mouse, and zebrafish (Danio rerio) as representative vertebrate species, based on their phylogenetic placement and sequence data quality. Removal of other vertebrate sequences avoids problems with oversampling, as well as problems associated with poor quality or missing data in some species. In addition to the three vertebrates, we included RPB1 amino acid sequences for 35 other eukaryotes, including worms, insects, plants, algae, amoebas, fungi, and several recently sequenced, early branching animals. For each species’ RPB1 sequence, we manually counted the number of CTD heptad repeats and the number of CTD repeats with lysine residues (Additional file 1). To be classified as an acceptable heptad, three of the seven amino acids had to match the consensus sequence or common alterations seen in multiple species.
We modified a species tree downloaded from the NCBI Taxonomy Database to reflect recent research on the evolutionary relationships between the clades considered and the approximate timing of divergence events [36,39]. To model losses and gains of CTD repeats and lysine content within each ancestral clade, we used Wagner parsimony (with an equal weight for gains and losses) as implemented in the Count program for analysis of numerical observations on a phylogeny .
Estimates of the phylogenetic age of human and mouse genes were made with ProteinHistorian using asymmetric Wagner parsimony on the PPOD-PANTHER protein-family database . ProteinHistorian was also used to identify significant differences in the distribution of gene ages between gene sets of interest, such as those occupied by different polymerase forms or in different functional classes.
Sequence logo generation
Sequence logos for clades were created using the online tool LogOddsLogo . In generation of the summary logos for clades, all heptad repeats from the specified species were input into the generator.
Identification of acRPB1 occupied and enriched genes
To determine promoter RPB1 occupancy genome-wide, we used ChIP-Seq data collected with antibodies for total RPB1 and acetylated RPB1 in mouse embryonic stem cells from . Following their definitions, we considered a gene to be “RPB1 occupied” if it had total RPB1 signal at its promoter (2 kilobases upstream of the transcription start site) greater than twice the input signal at the promoter (total/input > 2). Similarly, we defined “acRPB1 occupied” genes as those with acRPB1 occupancy twice that of the promoter input signal (ac/input > 2). Since the amount of acRPB1 at a promoter strongly correlated with the amount of RPB1 , we defined a set of “acRPB1 enriched” genes that were RPB1 occupied and had acRPB1 promoter occupancy at least twice the RPB1 occupancy (total/input > 2 and ac/total > 2). The genes in each set are listed in Additional file 4.
Identification of genes dysregulated with the disruption acRPB1
We used a HA-tagged mouse RPB1 construct in which all K7 residues were substituted with arginines (8KR) . This mutation resembles unacetylated lysines by conserving the positive charge at these positions, but preventing acetylation. To examine the potential functions of K7 acetylation in regulating gene expression, we stably expressed wildtype or 8KR HA-RPB1 in murine NIH/3 T3 fibroblasts and cultured these cells in media containing α-amanitin. Both were expressed at equivalent levels, but acetylation was present only in wildtype, and not mutant HA-RPB1. We then performed gene expression profiling using the Affymetrix Mouse Gene 1.0 ST microarray with three biological replicates. Array values were normalized and log2 scaled. For statistical analyses, we removed all array probe sets in which no experimental groups had an average log2 intensity greater than 3.0. This is a standard cutoff, below which expression is indistinguishable from background noise. Linear models were fit for each gene with the Bioconductor “limma” package in R [42,43]. Moderated t-statistics, fold-change and the associated p-values were calculated for each gene. To account for the fact that thousands of genes were tested, we controlled the false discovery rate (FDR) using the Benjamini-Hochberg method . Genes with p-values corresponding to expected FDRs of 0.05 or less were considered dysregulated compared to wildtype. These genes are listed in Additional file 4. The raw array data are available in the Gene Expression Omnibus (GSE66088).
Functional annotation of genes
We used the online functional annotation tool DAVID [27,28] to calculate KEGG pathway and gene ontology functional annotation enrichment for each set of genes of interest. EGF-induced and PDGF-induced genes were taken from  and . The human gene identifiers from these studies were mapped to their mouse homologs using the HGND database.
Availability of supporting data
The data sets supporting the results of this article are available as additional files. The microarray data we collected are available from the Gene Expression Omnibus (accession number GSE66088), and the previously collected ChIP-Seq data  are available from the Sequence Read Archive (accession number SRX338012).
We thank John Stiller for providing RPB1 sequences. We thank Sean Thomas and Alex Williams for help with data processing and analysis. We thank Antonis Rokas for helpful discussions. This work was supported by the National Institutes of Health (5T32EY021453 to CNS, R01GM082901 to KSP, and R01AI083139 to MO), the Boehringer Ingelheim Fund (SS), institutional funds from the Gladstone Institutes (KSP), and institutional funds from Vanderbilt University (JAC).
- Stiller JW, Cook MS. Functional unit of the RNA polymerase II C-terminal domain lies within heptapeptide pairs. Eukaryot Cell. 2004;3(3):735–40.View ArticlePubMed CentralPubMedGoogle Scholar
- Cramer P, Bushnell DA, Kornberg RD. Structural Basis of Transcription: RNA Polymerase II at 2.8 Ångstrom Resolution. Science. 2001;292(5523):1863–76.View ArticlePubMedGoogle Scholar
- Eick D, Geyer M. The RNA polymerase II carboxy-terminal domain (CTD) code. Chem Rev. 2013;113(11):8456–90.View ArticlePubMedGoogle Scholar
- Corden JL. RNA polymerase II C-terminal domain: Tethering transcription to transcript and template. Chem Rev. 2013;113(11):8423–55.View ArticlePubMed CentralPubMedGoogle Scholar
- Egloff S, Dienstbier M, Murphy S. Updating the RNA polymerase CTD code: adding gene-specific layers. Trends Genet: TIG. 2012;28(7):333–41.View ArticlePubMedGoogle Scholar
- Phatnani HP, Greenleaf AL. Phosphorylation and functions of the RNA polymerase II CTD. Genes Dev. 2006;20(21):2922–36.View ArticlePubMedGoogle Scholar
- Komarnitsky P, Cho E-J, Buratowski S. Different phosphorylated forms of RNA polymerase II and associated mRNA processing factors during transcription. Genes Dev. 2000;14:2452–60.View ArticlePubMed CentralPubMedGoogle Scholar
- Fujita T, Ryser S, Tortola S, Piuz I, Schlegel W. Gene-specific recruitment of positive and negative elongation factors during stimulated transcription of the MKP-1 gene in neuroendocrine cells. Nucleic Acids Res. 2007;35(3):1007–17.View ArticlePubMed CentralPubMedGoogle Scholar
- Egloff S, O’Reilly D, Chapman RD, Taylor A, Tanzhaus K, Pitts L, et al. Serine-7 of the RNA Polymerase II CTD Is Specifically Required for snRNA Gene Expression. Science. 2007;318(5857):1777–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Baillat D, Hakimi M-A, Näär AM, Shilatifard A, Cooch N, Shiekhattar R. Integrator, a multiprotein mediator of small nuclear RNA processing, associates with the C-terminal repeat of RNA polymerase II. Cell. 2005;123(2):265–76.View ArticlePubMedGoogle Scholar
- Sims RJ, Rojas LA, Beck D, Bonasio R, Schüller R, Drury WJ, et al. The C-terminal domain of rna polymerase II Is modified by site-specific methylation. Science. 2011;332(6025):99–103.View ArticlePubMed CentralPubMedGoogle Scholar
- Hsin J-P, Sheth A, Manley JL. RNAP II CTD phosphorylated on threonine-4 is required for histone mRNA 3′ end processing. Science. 2011;334(6056):683–6.View ArticlePubMed CentralPubMedGoogle Scholar
- Mayer A, Heidemann M, Lidschreiber M, Schreieck A, Sun M, Hintermair C, et al. CTD tyrosine phosphorylation impairs termination factor recruitment to RNA polymerase II. Science. 2012;336(6089):1723–5.View ArticlePubMedGoogle Scholar
- Hsin JP, Li W, Hoque M, Tian B, Manley JL. RNAP II CTD tyrosine 1 performs diverse functions in vertebrate cells. ELife. 2014;3:e02112.View ArticlePubMed CentralPubMedGoogle Scholar
- Descostes N, Heidemann M, Spinelli L, Schüller R, Maqbool MA, Fenouil R, et al. Tyrosine phosphorylation of RNA polymerase II CTD is associated with antisense promoter transcription and active enhancers in mammalian cells. Elife. 2014;3:e02105.View ArticlePubMed CentralPubMedGoogle Scholar
- Yang C, Stiller JW. Evolutionary diversity and taxon-specific modifications of the RNA polymerase II C-terminal domain. Proc Natl Acad Sci U S A. 2014;111(16):5920–5.View ArticlePubMed CentralPubMedGoogle Scholar
- Liu P, Kenney JM, Stiller JW, Greenleaf AL. Genetic organization, length conservation, and evolution of RNA polymerase II carboxyl-terminal domain. Mol Biol Evol. 2010;27(11):2628–41.View ArticlePubMed CentralPubMedGoogle Scholar
- Chapman RD, Conrad M, Eick D. Role of the mammalian RNA polymerase II C-terminal domain (CTD) nonconsensus repeats in CTD stability and cell proliferation. Mol Cell Biol. 2005;25(17):7665–74.View ArticlePubMed CentralPubMedGoogle Scholar
- Gerber H-P, Hagmann M, Seipel K, Georgiev O, West MAL, Litingtung Y, et al. RNA polymerase II C-terminal domain required for enhancer-driven transcription. Nature. 1995;374(6523):660–2.View ArticlePubMedGoogle Scholar
- Schroder S, Herker E, Itzen F, He D, Thomas S, Gilchrist DA, et al. Acetylation of RNA polymerase II regulates growth-factor-induced gene transcription in mammalian cells. Mol Cell. 2013;52(3):314–24.View ArticlePubMed CentralPubMedGoogle Scholar
- Capra JA, Stolzer M, Durand D, Pollard KS. How old is my gene? Trends Genet: TIG. 2013;29(11):659–68.View ArticlePubMedGoogle Scholar
- Capra JA, Williams AG, Pollard KS. ProteinHistorian: Tools for the comparative analysis of eukaryote protein origin. PLoS Comput Biol. 2012;8(6):e1002567.View ArticlePubMed CentralPubMedGoogle Scholar
- Kishore SP, Perkins SL, Templeton TJ, Deitsch KW. An unusual recent expansion of the C-terminal domain of RNA polymerase II in primate malaria parasites features a motif otherwise found only in mammalian polymerases. J Mol Evol. 2009;68(6):706–14.View ArticlePubMed CentralPubMedGoogle Scholar
- Sebé-Pedrós A, de Mendoza A, Lang BF, Degnan BM, Ruiz-Trillo I. Unexpected repertoire of metazoan transcription factors in the unicellular holozoan Capsaspora owczarzaki. Mol Biol Evol. 2011;28(3):1241–54.View ArticlePubMed CentralPubMedGoogle Scholar
- Bordoli L, Netsch M, Luthi U, Lutz W, Eckner R. Plant orthologs of p300/CBP: conservation of a core domain in metazoan p300/CBP acetyltransferase-related proteins. Nucleic Acids Res. 2001;29(3):589–97.View ArticlePubMed CentralPubMedGoogle Scholar
- Allis CD, Berger SL, Cote J, Dent S, Jenuwien T, Kouzarides T, et al. New nomenclature for chromatin-modifying enzymes. Cell. 2007;131(4):633–6.View ArticlePubMedGoogle Scholar
- da Huang W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37(1):1–13.View ArticlePubMed CentralPubMedGoogle Scholar
- da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.View ArticlePubMedGoogle Scholar
- Amit I, Citri A, Shay T, Lu Y, Katz M, Zhang F, et al. A module of negative feedback regulators defines growth factor signaling. Nat Genet. 2007;39(4):503–12.View ArticlePubMedGoogle Scholar
- Tullai JW, Schaffer ME, Mullenbrock S, Sholder G, Kasif S, Cooper GM. Immediate-early and delayed primary response genes are distinct in function and genomic architecture. J Biol Chem. 2007;282(33):23981–95.View ArticlePubMed CentralPubMedGoogle Scholar
- Rokas A. The molecular origins of multicellular transitions. Curr Opin Genet Dev. 2008;18(6):472–8.View ArticlePubMedGoogle Scholar
- Rokas A. The origins of multicellularity and the early history of the genetic toolkit for animal development. Annu Rev Genet. 2008;42:235–51.View ArticlePubMedGoogle Scholar
- Shilo B-Z. Regulating the dynamics of EGF receptor signaling in space and time. Development. 2005;132(18):4017–27.View ArticlePubMedGoogle Scholar
- Adelman K, Lis JT. Promoter-proximal pausing of RNA polymerase II: emerging roles in metazoans. Nat Rev Genet. 2012;13(10):720–31.View ArticlePubMed CentralPubMedGoogle Scholar
- Kruesi WS, Core LJ, Waters CT, Lis JT, Meyer BJ. Condensin controls recruitment of RNA polymerase II to achieve nematode X-chromosome dosage compensation. ELife. 2013;2:e00808.View ArticlePubMed CentralPubMedGoogle Scholar
- Richter DJ, King N. The genomic and cellular foundations of animal origins. Annu Rev Genet. 2013;47:509–37.View ArticlePubMedGoogle Scholar
- Kishore SP, Stiller JW, Deitsch KW. Horizontal gene transfer of epigenetic machinery and evolution of parasitism in the malaria parasite Plasmodium falciparum and other apicomplexans. BMC Evol Biol. 2013;13:37.View ArticlePubMed CentralPubMedGoogle Scholar
- Stiller JW, Hall BD. Evolution of the RNA polymerase II C-terminal domain. Proc Natl Acad Sci U S A. 2002;99(9):6091–6.View ArticlePubMed CentralPubMedGoogle Scholar
- Hedges SB, Dudley J, Kumar S. TimeTree: a public knowledge-base of divergence times among organisms. Bioinformatics. 2006;22(23):2971–2.View ArticlePubMedGoogle Scholar
- Csuros M. Count: evolutionary analysis of phylogenetic profiles with parsimony and likelihood. Bioinformatics. 2010;26(15):1910–2.View ArticlePubMedGoogle Scholar
- Yu Y-K, Capra JA, Stojmirović A, Landsman D, Altschul SF. Log-odds sequence logos. Bioinformatics. 2015;31(3):324–31.View ArticlePubMedGoogle Scholar
- Gentleman R, Carey V, Bates D, Bolstad B, Dettling M, Dudoit S, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5(10):R80.View ArticlePubMed CentralPubMedGoogle Scholar
- Smyth GK. Linear Models and Empirical Bayes Methods for Assessing Differential Expression in Microarray Experiments. Stat Appl Genet Mol Biol. 2004;3(1):1–25. ISSN (Online) 1544-6115, doi:10.2202/1544-6115.1027.Google Scholar
- Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B. 1995;57(1):289–300.Google Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.