Comparative genomics of cyclin-dependent kinases suggest co-evolution of the RNAP II C-terminal domain and CTD-directed CDKs

Abstract

Background: Cyclin-dependent kinases (CDKs) are a large family of proteins that function in a variety of key regulatory pathways in eukaryotic cells, including control over the cell cycle and gene transcription. Among the most important and broadly studied of these roles is reversible phosphorylation of the C-terminal domain (CTD) of RNA polymerase II, part of a complex array of CTD/protein interactions that coordinate the RNAP II transcription cycle. The RNAP CTD is strongly conserved in some groups of eukaryotes, but highly degenerate or absent in others; the reasons for these differences in stabilizing selection on CTD structure are not clear. Given the importance of reversible phosphorylation for CTD-based transcription, the distribution and evolutionary history of CDKs may be a key to understanding differences in constraints on CTD structure; however, the origins and evolutionary relationships of CTD kinases have not been investigated thoroughly. Moreover, although the functions of most CDKs are reasonably well studied in mammals and yeasts, very little is known from most other eukaryotes.

Results: Here we identify 123 CDK family members from animals, plants, yeasts, and four protists from which genome sequences have been completed, and 10 additional CDKs from incomplete genome sequences of organisms with known CTD sequences. Comparative genomic and phylogenetic analyses suggest that cell-cycle CDKs are present in all organisms sampled in this study. In contrast, no clear orthologs of transcription-related CDKs are identified in the most putatively ancestral eukaryotes, Trypanosoma or Giardia. Kinases involved in CTD phosphorylation, CDK7, CDK8 and CDK9, all are recovered as well-supported and distinct orthologous families, but their relationships to each other and other CDKs are not well-resolved. Significantly, clear orthologs of CDK7 and CDK8 are restricted to only those organisms belonging to groups in which the RNAP II CTD is strongly conserved.

Conclusions: The apparent origins of CDK7 and CDK8, or at least their conservation as clearly recognizable orthologous families, correlate with strong stabilizing selection on RNAP II CTD structure. This suggests co-evolution of the CTD and these CTD-directed CDKs. This observation is consistent with the hypothesis that CDK7 and CDK8 originated at about the same time that the CTD was canalized as the staging platform RNAP II transcription. Alternatively, extensive CTD phosphorylation may occur in only a subset of eukaryotes and, when present, this interaction results in greater stabilizing selection on both CTD and CDK sequences. Overall, our results suggest that transcription-related kinases originated after cell-cycle related CDKs, and became more evolutionarily and functionally diverse as transcriptional complexity increased.

Figure 1

Unrooted 50% majority consensus tree from 4,000 ML trees sampled from the Bayesian posterior probability distribution. Support values are shown above the internode from Bayesian inference/distance bootstrap respectively. Only values above 50% are reported and values under 50% are indicated by (-). 100% values are indicated by (+). CDK names in blue are from organisms that fall into the "CTD-clade" in RPB1 phylogenetic analyses (see Fig. 2); and those in red are from groups in which the CTD is not strongly conserved. Inferred groups of CTD-directed CDKs 7, 8 and 9 are shown in bold. A large group of unidentified CDKs from Arabidopsis and Oryza, which appear to represent a plant-specific amplification of CTK9, were excluded from this analysis to determine whether identified plant CDK9s show a specific phylogenetic affinity to either the BUR1 or CTK1 subgroup. All identified plant sequences are included in an expanded analysis shown in additional file 1.

Figure 2

Hypothesis of RNA polymerase II evolution inferred from phylogenetic analyses of RPB1 sequences conserved regions A-H. The tree displayed, after Stiller and Cook [60] had the highest likelihood of all trees sampled from the posterior probability distribution in 10⁶ generations of Bayesian inference. Organisms with genomes included in this study are in larger/bold font, and whether each of the three primary CTD kinases (CDKs7,8,9) are present in this genome, as inferred from global phylogenetic analyses and distributions of CDK kinases (see Fig. 1), are indicated next to the name. Names in blue indicate the presence of repeated heptads at the RPB1 C-terminus, which includes several from protistan organisms that do not conform to the consensus sequence or known structural requirements of the canonical CTD [60]; names in red have no tandem-heptapeptide structure whatsoever. The node supporting a "CTD-clade," in which the consensus sequence and repetitive structure of the CTD are invariably conserved, occurred in 98% of the 8000 trees sampled from the Bayesian posterior probability distribution. See references 31 and 60 for a more complete phylogenetic treatment of the origin and conservation of the CTD.