Figures
Abstract
Background
NIMA-related kinases (Neks) have been studied in diverse eukaryotes, including the fungus Aspergillus and the ciliate Tetrahymena. In the former, a single Nek plays an essential role in cell cycle regulation; in the latter, which has more than 30 Neks in its genome, multiple Neks regulate ciliary length. Mammalian genomes encode an intermediate number of Neks, several of which are reported to play roles in cell cycle regulation and/or localize to centrosomes. Previously, we reported that organisms with cilia typically have more Neks than organisms without cilia, but were unable to establish the evolutionary history of the gene family.
Methodology/Principle Findings
We have performed a large-scale analysis of the Nek family using Bayesian techniques, including tests of alternate topologies. We find that the Nek family had already expanded in the last common ancestor of eukaryotes, a ciliated cell which likely expressed at least five Neks. We suggest that Neks played an important role in the common ancestor in regulating cilia, centrioles, and centrosomes with respect to mitotic entry, and that this role continues today in organisms with cilia. Organisms that lack cilia generally show a reduction in the number of Nek clades represented, sometimes associated with lineage specific expansion of a single clade, as has occurred in the plants.
Citation: Parker JDK, Bradley BA, Mooers AO, Quarmby LM (2007) Phylogenetic Analysis of the Neks Reveals Early Diversification of Ciliary-Cell Cycle Kinases. PLoS ONE 2(10): e1076. https://doi.org/10.1371/journal.pone.0001076
Academic Editor: Rosemary Redfield, University of British Columbia, Canada
Received: September 4, 2007; Accepted: October 4, 2007; Published: October 24, 2007
Copyright: © 2007 Parker et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was supported by a Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery Grant to LMQ. BAB was supported by fellowships from NSERC and the Michael Smith Foundation for Health Research. Neither of these agencies had any role in any aspect of this study or the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
The NIMA-related family of serine/threonine kinases (Neks) are widespread among eukaryotes. They are defined by similarity in their N-terminal catalytic domains to the founding member, Never In Mitosis A from the fungus Aspergillus [1], [2]. NIMA, the sole Nek in the Aspergillus genome, plays multiple roles in cell cycle progression and localizes to the fungal equivalent of the centrosome, the spindle pole body [1], [3]. There are 11 NIMA-related kinase (Nek) genes in the human genome. Those that have been studied appear to have cell cycle-related functions and some localize to centrosomes [4]–[11].
Centrosomes are typically located near the nucleus in the centre of the cell, where they serve as organizers of the microtubular cytoskeleton during both interphase and mitosis [12]. In addition to their well-established role as microtubule organizers, centrosomes serve as important scaffolds for the integration of multiple signaling pathways that coordinate cell cycle progression [13]. Of the centrosomal Neks, only for Nek2 has a function been clearly defined. Nek2 phosphorylates proteins which physically link centrioles, the microtubule-based structures at the core of centrosomes [14]–[16].
In mammalian cells, the elder centriole of the centrosomal pair often directly nucleates a cilium. Cilia have functions in both motility and sensory signaling [17]. The recent demonstration that flies without centrioles develop normally but then die due to multiple sensory defects [18] suggests that the most important function of centrioles may be to nucleate sensory cilia, at least in some lineages. An indication of the ciliary function of some Neks derives from their association with polycystic kidney diseases, which are caused by defective ciliary signaling [19]: The causative mutations in two mouse models of polycystic kidney diseases are in the mNek1 and mNek8 genes [20], [21]. We have shown that Nek8 is ciliary [11] and mutations that affect ciliary localization are associated with a rare form of a cystic kidney disease that affects children, Nephronophthisis type 9 [22]. Although the functions of Neks in mammalian cilia are unknown, it has been established that Neks in Chlamydomonas and Tetrahymena regulate ciliary length and/or disassembly [23]–[26]. We have previously noted a correlation between the number of Neks a eukaryotic organism expresses, and the presence of ciliated cells which re-enter the cell cycle in that organism [2]. Together, these observations have led us to propose that Neks evolved with centrosomes and serve to coordinate ciliary and cell cycle functions [2], [27].
The breadth of organisms in which Neks are found and the lack of sequence conservation outside the kinase domains has meant no large-scale phylogenetic analysis has been possible via traditional methods. Many authors refer to NIMA as “ancestral” to the mammalian Neks, yet no eukaryotic phylogeny would ever place Aspergillus as an ancestor to mammals. The assumption by these authors is that the last common ancestor of Fungi and Metazoa had one Nek, similar to NIMA, and this gene underwent an expansion in Metazoa and especially in mammals. This conflicts with previously published phylogenies [26], [28] which suggest that there are orthologies between mammalian Neks and Neks in both Chlamydomonas and Tetrahymena. However, these earlier Nek phylogenies had few taxa and poor support at many nodes. We set out to establish a more complete Nek family tree, by sampling distantly-related phyla, by using Bayesian inference methods, and by testing specific hypotheses in a Bayesian framework. To our knowledge, no similar analyses have been published for a eukaryotic kinase family. Our results suggest that the last common ancestor of the eukaryotes, which must have been a ciliated cell able to divide [29], had at least five Neks. This finding and new cell biological data further strengthen the three-way bond between Neks, cilia/centrosomes, and the cell cycle.
Materials and Methods
Species from a broad range of eukaryotic lineages, with publicly-available sequenced genomes, were selected for analysis (Figure 1). We generated a dataset of Nek kinase domains from these species via the following methods: with the NIMA kinase domain as a query sequence, we used BLASTP analysis [30] to identify Neks from the predicted proteomes, using their publicly-available genomes. Top scoring BLASTP hits from the predicted proteomes were then reciprocally used as queries of the non-redundant NCBI GenBank database. Only those queries that produced Nek proteins as top scoring BLASTP hits were retained in the dataset. Only kinase domains, as predicted by the SMART HMMER search tool [31], [32] were included in our phylogenetic analysis, leading us to exclude the predicted human Nek10 (GenBank id NP_001026911), amongst other weak hits from various other organisms. Multiple sequence alignments were produced using ClustalX software [33] and manually inspected. Any kinase domain sequences containing large insertions or truncated kinase domains, likely the result of incorrect protein prediction, were excluded from the dataset at this stage. Our master list of kinase domains is available as supplemental data (Table S1).
Relationships among taxa based on [29], [55]. The “Ciliated cells divide” column indicates presence (Y) or absence (N) of cells in the species which both bear cilia and re-enter the cell cycle: gametes and other terminally differentiated ciliated cells are excluded. The “Neks” column indicates both how many genes encoding Nek kinase domains were identified in the indicated species. The number in parentheses is a minimal estimate of how many Nek subfamilies are represented in the genome, as determined from our analysis (see text). N.D. indicates that we did not determine the number of subfamilies for the particular species.
In the course of our BLAST searching, we found some prokaryotic kinases which were reciprocal best hits with Neks. However, these sequences do not represent bona fide Neks, as they are outgroup to a clade containing multiple eukaryotic kinase families (not shown).
We used the Metropolis-coupled Markov chain Monte Carlo program MrBayes v. 3.1.2 [34], [35] to infer phylogenies. We allowed the program to determine the most probable model of protein evolution (aamodel), and specified an inverse gamma distribution of variability in substitution rate. Each analysis was composed of two simultaneous runs of either four or eight MCMC chains run for at least 2×106 generations, sampled every 100 generations, and run until the average standard deviation of split frequencies between the two runs remained below the 0.01 convergence threshold for at least 1×105 generations. Ninety percent of each run was burned in before the parameters were analyzed. Human Cdk2 was used to root all trees in either Treeview [36] or HyperTree [37]. HsCdk2 was selected as the outgroup because on a tree including Cdk2, Mlk3, AuroraA, Plk1 and the mammalian Neks, Cdk2 was the kinase most closely related to the Neks (data not shown). See Supplementary Table 1 for detailed tree statistics.
To test hypotheses on the unikonts plus Chlamydomonas dataset we first analyzed the unconstrained dataset and obtained a harmonic mean of the log likelihood values of the MCMC samples (Table 1). We then re-ran the analysis, constraining the topology to reflect various hypotheses (for example, constraining nodes or pairs of nodes, Table 1), and recorded the harmonic mean of log likelihood values. Twice the difference in mean harmonic likelihood (known as the Bayes Factor) can be used to determine relative support for alternative hypotheses; values for the Bayes Factors of >10 indicate strong support for the less constrained tree over the a priori imposed constraint (i.e. the hypothesis). [38], [39]. Analyses were otherwise run as described above. Constraining the amino acid model did not affect either the topology or likelihood of the resulting tree (Table 1).
Results
We use the terminology of Cavalier-Smith [29] to refer to the major eukaryotic lineages: unikonts include Metazoa, Fungi, and Amoebozoa, while bikonts include all other organisms considered. We found Neks from a broad sampling of eukaryotic lineages (Figure 1), but we identified no Neks in the genomes of the diatom Thalassiosira or the green alga Ostreococcus.
Large numbers of Neks were found in the genomes of the kinetoplastids Leishmania and Trypanosoma, as well in the genome of Giardia (Figure 1). Inclusion of sequences from either the kinetoplastids or Giardia in our analyses lead to failure of convergence, so we excluded these species from further analysis. We note that these species are parasitic, and recent analyses consider them to be highly derived, rather than basal, eukaryotes [40]. This, in combination with their large number of Nek sequences, accounts for the difficulty in resolving Nek trees including these organisms. However, our analysis indicates that most of the kinetoplastid Neks form a lineage-specific expansion in a clade with human Nek8 ( Figure S1). Our SMART analysis showed that with rare exceptions (all of which are in Tetrahymena) all sequences recognized as having a Nek-like kinase domain had that domain at the N-terminus of the predicted protein; when C-terminal domains were predicted by SMART, these were always protein-protein interaction domains. These results are consistent with the view that Neks have rapidly-evolving C-terminal regulatory/interaction domains, allowing differential localization, regulation, or function in various lineages.
A dataset consisting of the unikont Neks plus the Chlamydomonas Neks was the largest dataset we could analyze in a reasonable time (∼1 week on a desktop computer), while inclusion of the bikont Neks to produce the tree in Figure 3 required much more analysis time (∼3 months) on a computer cluster using the equivalent of eighteen desktop processors. Due to these computational limitations, we performed stringent analyses by Bayesian testing of alternate topologies on the former dataset, represented in Figure 2 and Table 1 (see below).
Numbers at nodes are interpreted as a posteriori support values. Bold digits at some nodes indicate nodes used for hypothesis testing, and the Bayes Factor (BF) where indicated. Right brackets include sequences joined by well-supported internal nodes, which we propose to refer to as Nek subfamilies, named after the human member with the shortest branch length.
Neks from all species shown in Figure 1 having Neks in their genome were included in this dataset except for Leishmania, Trypanosoma, and Giardia. This represents an overview of the tree which is examined in detail in Figures 4–7.
The topology of the tree in Figure 2 is consistent with previous studies (e.g. [41]), but the inclusion of several additional species in our tree allows us to differentiate more clearly orthology and paralogy. For example, analysis of the unikont Neks (Figure 2) suggests that hNek2 is indeed the NIMA orthologue, based on the observations that the nodes connecting these taxa are strongly-supported (posterior probability typically ∼1), and contain sequences from all the unikont taxa included, including distantly-related taxa such as humans and yeast. Consequently, we refer to this group of sequences as the Nek2 clade. Similarly, we note that the human Neks are scattered among six a posteriori-defined clades, each of which contains representatives from most species tested and most of which have well-supported nodes (posterior probability >0.85, Figure 2). Importantly, the above statement is valid despite the inclusion of Chlamydomonas; that is, most Chlamydomonas Neks fall into the same well-supported clades as the Neks of the unikont species included in this analysis.
Four of the ten Chlamydomonas Neks are not members of any clades in Figure 2. In order to test whether this was a bona fide result, we tested various alternative hypotheses for the topology. We constrained Cnk1p and Cnk2p (node 5 in Figure 2) to be sister to the Nek2 clade (node 1 in Figure 2), and calculated the resulting Bayes Factor (Table 1, “node 5+node 1”). The Bayes Factor of 42 strongly rejects a composite Cnk2/Nek2 clade in favour of the topology seen in Figure 2. As controls, we repeated this analysis to test alternate topologies for the Chlamydomonas sequences which fell within well-supported clades in Figure 2. Reassuringly, the topology reported in Figure 2 was not rejected for any of these hypotheses, including Cnk10p membership in the same clade as vertebrate Nek1 (Table 1). We repeated this experiment for less well-supported nodes, e.g., the Nek4 and Nek11 clades. Cnk5p and Cnk8p are sister to unikont Nek4's with a posterior probability of only 0.54, but could not be rejected as members of the Nek4 clade by Bayes Factor Analysis (BF = 7.04); a similar result was obtained indicating that Cnk9p is a tentative member of the Nek11 clade.
In only three instances did Bayes Factor analysis reveal a lack of support for the tree reported in Figure 2: the Dictyostelium DDB0229345 sequence was rejected from the Fa2p/Cnk4p clade, the C.elegans sequence D1044.8 was rejected from the Nek2 clade, and the Cnk2p/Cnk1p and Fa2p/Cnk4p clades cannot be distinguished as independent clades (Table 1). In summary, the analyses presented in Figure 2 and Table 1 indicate at least five bona fide Nek clades with members from both the mammal Homo and the green alga Chlamydomonas, plus two additional clades, the Chlamydomonas-specific Fa2p/Cnk2p clade, and the unikont-specific Nek2 clade.
Addition of Nek sequences from other bikonts, as well as the (unikont) chytrid fungus Batrachochytrium to the data set used to generate Figure 2 did not significantly affect the topology of the tree (Figure 3). Several of the newly-added sequences arise in the tree as polytomies, a result consistent with the difficulties inherent in constructing a multigene phylogeny across a broad sampling of eukaryotes. However, with the exception of the Nek1/3/5 clade, but including the Nek2 and Cnk2p/Fa2p clades, our previously-defined clades now include sequences from distantly related species such as Arabidopsis, Naegleria, and Tetrahymena (Figures 4–7). The previously unikont-only Nek2 clade has members from Tetrahymena, while two apicocomplexan sequences (Cp cgd1 1490 and Pf Nek1) branch near the well-supported Nek2 clade (Figure 4). The Nek6/7 clade is well-supported, and like the Nek2 clade, now contains sequences from Tetrahymena and Batrachochytrium (Figure 5). Satisfyingly, Cnk2p and Fa2p now form a clade with several Tetrahymena Neks (consistent with Wloga et al. [26]), as well as a Naegleria Nek (Figure 5). Similarly, the Nek8/9 clade remains intact (Figure 6), and now contains Naegleria “Nek2”, which has a similar C-terminal domain composition to its orthologue in Chlamydomonas, Cnk3p (not shown).
For ease of visualization, large nodes have been collapsed (triangles) and are shown expanded in subsequent figures (Figures 5–7) as indicated. In these figures sequences from bikont species are cyan, while sequences from unikont species are red. Right brackets include sequences joined by well-supported internal nodes, which we propose to refer to as Nek subfamilies, named after the human member with the shortest branch length where appropriate.
For details see legend to Figure 4. The Nek4 clade is marked with an asterisk due to its poor support (see text).
The Nek1/3/5 clade did not resolve with sequences from both unikonts and bikonts (Figure 6). However, a tree from a smaller dataset (Figure S2) shows that Chlamydomonas Cnk10p, Tetrahymena Nrk29 and Cryptosporidium cgd7 3760 are members of a well-supported Nek1/3/5 clade (posterior probability of 0.94). This suggests that some of the clades arising from the polytomy could be orthologous to members of the unikont Nek1/3/5 clade. Above, we provided additional evidence that Cnk10p is a Nek1/3/5 family member (Figure 2, Table 1, Figure S1, Figure S2). Unfortunately, tests of hypotheses on this larger dataset are not feasible with our current methods and available computational power.
Consistent with our previous results (Figure 2 and Table 1), the Nek4 and Nek11 clades are closely related to each other and thus difficult to resolve (Figure 7). We note that human Nek11 (a protein localized to the nucleolus [9]) and the higher plant Neks are all in the Nek11 clade, consistent with the hypothesis that a non-centrosomal Nek family expanded in a lineage which does not have centrosomes; however, these Neks still play roles in cell cycle regulation [42], [43]. In contrast, we find Tetrahymena Nrk30p, which localizes to cilia [26], to be in the Nek4 clade (Figure 7). We therefore hypothesize that in contrast to the Nek11 clade, members of the Nek4 group may retain ciliary function.
Neks were found in nearly every genome sampled, and the numbers of Neks we identified roughly correlates with whether the organism has non-terminally differentiated ciliated cells (Figure 1). We included both the raw number of Neks found in the organisms' genomes, and the minimal number of Nek subfamilies represented, to further clarify the Nek/cilia relationship. We define our Nek subfamilies as moderately well-supported clades (posterior probability >0.85) with members from both unikont and bikont species. The Nek2, Nek6/7, Nek8/9, and Nek11 clades meet this definition (Figures 4–7). As discussed above, the Nek1/3/5 clade is well-supported in Figure S2. Also, the Nek 11 clade in Figure 7 meets our definition, but the Nek4 clade does not; however, a composite Nek4/11 clade is well-supported in our smaller dataset (Figure S2). Thus, we define five Nek subfamilies, named after their human members: Nek1/3/5, Nek2, Nek4/11, Nek6/7, and Nek8/9. Our analysis thus allows identification of lineage-specific expansions, for example, the plant Arabidopsis has seven Neks in its genome, but all seven are part of an expansion in the same Nek subfamily (Figures 3). Conversely, Batrachochytrium has only four Neks, but these all belong to separate subfamilies. Thus a correlation is evident amongst the organisms sampled between the presence of ciliated cells which divide and the diversity of Neks (Figure 1).
Discussion
We draw two major conclusions from our phylogenetic analysis. First, tree reconciliation suggests that the last common ancestor (LCA) of unikonts and bikonts most likely had at least five Neks and variations in the numbers of Neks in different lineages was due to expansion or loss of multiple Neks; this is a novel finding and changes our view of this kinase family and the relationships of its members. Second, this work confirms and extends our previous suggestion that Neks are more abundant in the genomes of organisms which must co-ordinate ciliary assembly and disassembly with respect to the cell cycle.
Currently, there is a paucity of functional data relevant to the cell biology of the Nek family. Consequently, there is insufficient data available to discern clear roles for particular Nek subfamilies from comparison of the Nek repertoires. However, the availability of the relationships between the various Neks will make future research on any particular Nek more relevant to the field at large. While most Neks that have been studied are associated with cilia/centrosomes/cell cycle, we do not propose that all Neks subserve this triumvirate. Indeed, expansion of the subfamilies was inevitably associated with both loss of function and acquisition of new roles [44], but analysis of which Neks have ancestral functions and which have derived functions requires the availability of phylogeny such as we present here. For example, Nek3, which is not centrosomal and has functions seemingly unrelated to microtubular-based structures [45], [46], is now revealed to be a relatively recent innovation (Figure 2), which likely underwent neofunctionalization.
Although we cannot discern cellular roles for individual Nek subfamilies, it is clear that Neks are not required for the formation of centrioles or cilia. Diatoms, for example, make both of these organelles de novo when differentiating into gametes [47] but have no Neks. Nor are the Neks required for mitosis, as exemplified by the many species whose genomes do not encode Neks. Importantly, we have yet to find a species that does not encode any Neks but does have ciliated cells that divide. Conversely, the Giardia genome has >70 Neks; Giardia flagella undergo a complex maturation process which takes four cell cycles [48], which may represent an extreme case of requirement for Neks during division in a flagellated cell. We predict that as functional studies progress, it will be revealed that many Neks participate in the important cellular role of coordinating the activities of centrioles as basal bodies during interphase and as spindle pole foci during mitosis.
Given our phylogenetic results, ciliary and centrosomal Nek functions are likely to be ancestral and thus may be retained in phylogenetically distant organisms. The wide distribution of the Neks from both unikonts and bikonts suggests strongly that the Neks had already expanded in the eukaryotic LCA. We suggest deviations of the topologies of the Nek subfamilies from our analysis from the species tree may represent different evolutionary constraints in the various lineages.
Our conclusion that the LCA, which must have been a non-terminally differentiated ciliated cell, had at least five Neks, should be interpreted in two ways: first, we are discerning the nature of the eukaryotic LCA, which is distinct from the first eukaryote; second, the presence of several Neks in the LCA conforms well with our idea that this kinase family is ancestrally correlated with the presence of ciliated cells which divide. The cell biology of the eukaryotic LCA can be partially inferred by applying the principal of parsimony to diverse extant eukaryotes: it is thus likely that the LCA was ciliated with a 9+2 axoneme and centrioles with triplet microtubules, and that the LCA could undergo both mitosis and meiosis [49]. The feature of extant eukaryotes most predictive for the number of retained Nek subfamilies appears to be the presence of ciliated cells which can re-enter the cell cycle (Figure 1 and [2]). Our data is consistent with a growing body of data which suggests that the eukaryotic last common ancestor had a sophisticated cellular organization and possessed many proteins required for centriole/cilium assembly [50], [51]. It is also consistent with the notion that early eukaryotes underwent a period of rapid evolution, especially of cytoskeletal elements [52]. We suggest that the Neks may have expanded during this time.
Although the relationship between cilia and the cell cycle is a long-standing problem in cell biology [53] which is receiving renewed interest [27], [54], our understanding of the complex regulatory relationship that must exist remains rudimentary. Functional studies of the Nek family may provide new points of entry into this important but undoubtedly complex signaling network.
Supporting Information
Figure S1.
Phylogram of the largest and most diverse dataset containing Leishmania and Trypanosoma sequences that was able to converge. This tree is rooted on HsCdk2. Please note that the majority of Leishmania and Trypanosoma Neks are members of a well-supported clade that includes HsNek8 (posterior probability = 0.90).
https://doi.org/10.1371/journal.pone.0001076.s002
(1.05 MB TIF)
Figure S2.
Phylogram of Nek kinase domains from species listed in Figure 1 with the exception of Batrachochytrium, Plasmodium, Naegleria, Arabidopsis, and Oryza. This tree is rooted on HsCdk2.
https://doi.org/10.1371/journal.pone.0001076.s003
(1.16 MB TIF)
Acknowledgments
We thank Martin Seigert of the SFU High Performance Computing Facility for assistance.
Author Contributions
Conceived and designed the experiments: JP BB. Performed the experiments: JP BB. Analyzed the data: JP BB. Wrote the paper: JP. Other: Provided expertise for the design and analysis of the experiments: AM. Conceived the project: LQ. Contributed to data analysis: LQ. Contributed to writing of the paper: LQ.
References
- 1. O'Connell MJ, Krien MJ, Hunter T (2003) Never say never. The NIMA-related protein kinases in mitotic control. Trends Cell Biol 13: 221–228.
- 2. Quarmby LM, Mahjoub MR (2005) Caught Nek-ing: cilia and centrioles. J Cell Sci 118: 5161–5169.
- 3. Osmani AH, McGuire SL, Osmani SA (1991) Parallel activation of the NIMA and p34cdc2 cell cycle-regulated protein kinases is required to initiate mitosis in A. nidulans. Cell 67: 283–291.
- 4. Belham C, Roig J, Caldwell JA, Aoyama Y, Kemp BE, et al. (2003) A mitotic cascade of NIMA family kinases. Nercc1/Nek9 activates the Nek6 and Nek7 kinases. J Biol Chem 278: 34897–34909.
- 5. Feige E, Shalom O, Tsuriel S, Yissachar N, Motro B (2006) Nek1 shares structural and functional similarities with NIMA kinase. Biochim Biophys Acta 1763: 272–281.
- 6. Fry AM (2002) The Nek2 protein kinase: a novel regulator of centrosome structure. Oncogene 21: 6184–6194.
- 7. Tan BC, Lee SC (2004) Nek9, a novel FACT-associated protein, modulates interphase progression. J Biol Chem 279: 9321–9330.
- 8. Yissachar N, Salem H, Tennenbaum T, Motro B (2006) Nek7 kinase is enriched at the centrosome, and is required for proper spindle assembly and mitotic progression. FEBS Lett 580: 6489–6495.
- 9. Noguchi K, Fukazawa H, Murakami Y, Uehara Y (2004) Nucleolar Nek11 is a novel target of Nek2A in G1/S-arrested cells. J Biol Chem 279: 32716–32727.
- 10. Yin MJ, Shao L, Voehringer D, Smeal T, Jallal B (2003) The serine/threonine kinase Nek6 is required for cell cycle progression through mitosis. J Biol Chem 278: 52454–52460.
- 11. Mahjoub MR, Trapp ML, Quarmby LM (2005) NIMA-related kinases defective in murine models of polycystic kidney diseases localize to primary cilia and centrosomes. J Am Soc Nephrol 16: 3485–3489.
- 12. Doxsey S, McCollum D, Theurkauf W (2005) Centrosomes in cellular regulation. Annu Rev Cell Dev Biol 21: 411–434.
- 13. Doxsey S, Zimmerman W, Mikule K (2005) Centrosome control of the cell cycle. Trends Cell Biol 15: 303–311.
- 14. Bahe S, Stierhof YD, Wilkinson CJ, Leiss F, Nigg EA (2005) Rootletin forms centriole-associated filaments and functions in centrosome cohesion. J Cell Biol 171: 27–33.
- 15. Mayor T, Hacker U, Stierhof YD, Nigg EA (2002) The mechanism regulating the dissociation of the centrosomal protein C-Nap1 from mitotic spindle poles. J Cell Sci 115: 3275–3284.
- 16. Fry AM, Mayor T, Meraldi P, Stierhof YD, Tanaka K, et al. (1998) C-Nap1, a novel centrosomal coiled-coil protein and candidate substrate of the cell cycle-regulated protein kinase Nek2. J Cell Biol 141: 1563–1574.
- 17. Marshall WF, Nonaka S (2006) Cilia: tuning in to the cell's antenna. Curr Biol 16: R604–614.
- 18. Basto R, Lau J, Vinogradova T, Gardiol A, Woods CG, et al. (2006) Flies without centrioles. Cell 125: 1375–1386.
- 19. Pazour GJ (2004) Intraflagellar transport and cilia-dependent renal disease: the ciliary hypothesis of polycystic kidney disease. J Am Soc Nephrol 15: 2528–2536.
- 20. Upadhya P, Birkenmeier EH, Birkenmeier CS, Barker JE (2000) Mutations in a NIMA-related kinase gene, Nek1, cause pleiotropic effects including a progressive polycystic kidney disease in mice. Proc Natl Acad Sci U S A 97: 217–221.
- 21. Liu S, Lu W, Obara T, Kuida S, Lehoczky J, et al. (2002) A defect in a novel Nek-family kinase causes cystic kidney disease in the mouse and in zebrafish. Development 129: 5839–5846.
- 22. Otto EA, Trapp ML, Schultheiss UT, Quarmby LM, Hildebrandt F (2007) Mutations in NIMA-related kinase NEK8 causes nephronophthisis in humans and affects ciliary and centrosomal localization. J. Am. Soc. Nephrology. In press.
- 23. Bradley BA, Quarmby LM (2005) A NIMA-related kinase, Cnk2p, regulates both flagellar length and cell size in Chlamydomonas. J Cell Sci 118: 3317–3326.
- 24. Mahjoub MR, Montpetit B, Zhao L, Finst RJ, Goh B, et al. (2002) The FA2 gene of Chlamydomonas encodes a NIMA family kinase with roles in cell cycle progression and microtubule severing during deflagellation. J Cell Sci 115: 1759–1768.
- 25. Mahjoub MR, Qasim Rasi M, Quarmby LM (2004) A NIMA-related kinase, Fa2p, localizes to a novel site in the proximal cilia of Chlamydomonas and mouse kidney cells. Mol Biol Cell 15: 5172–5186.
- 26. Wloga D, Camba A, Rogowski K, Manning G, Jerka-Dziadosz M, et al. (2006) Members of the NIMA-related kinase family promote disassembly of cilia by multiple mechanisms. Mol Biol Cell 17: 2799–2810.
- 27. Quarmby LM, Parker JD (2005) Cilia and the cell cycle? J Cell Biol 169: 707–710.
- 28. Bradley BA, Wagner JJ, Quarmby LM (2004) Identification and sequence analysis of six new members of the NIMA-related kinase family in Chlamydomonas. J Eukaryot Microbiol 51: 66–72.
- 29. Richards TA, Cavalier-Smith T (2005) Myosin domain evolution and the primary divergence of eukaryotes. Nature 436: 1113–1118.
- 30. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ (1990) Basic local alignment search tool. J Mol Biol 215: 403–410.
- 31. Letunic I, Copley RR, Pils B, Pinkert S, Schultz J, et al. (2006) SMART 5: domains in the context of genomes and networks. Nucleic Acids Res 34: D257–260.
- 32. Schultz J, Copley RR, Doerks T, Ponting CP, Bork P (2000) SMART: a web-based tool for the study of genetically mobile domains. Nucleic Acids Res 28: 231–234.
- 33. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG (1997) The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 25: 4876–4882.
- 34. Huelsenbeck JP, Ronquist F (2001) MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics 17: 754–755.
- 35. Ronquist F, Huelsenbeck JP (2003) MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572–1574.
- 36. Page RD (1996) TreeView: an application to display phylogenetic trees on personal computers. Comput Appl Biosci 12: 357–358.
- 37. Bingham J, Sudarsanam S (2000) Visualizing large hierarchical clusters in hyperbolic space. Bioinformatics 16: 660–661.
- 38. Kass R, Raftery A (1995) Bayes Factors. Journal of the American Statistical Association 90: 773–795.
- 39. Nylander JA, Ronquist F, Huelsenbeck JP, Nieves-Aldrey JL (2004) Bayesian phylogenetic analysis of combined data. Syst Biol 53: 47–67.
- 40. Simpson AG, Stevens JR, Lukes J (2006) The evolution and diversity of kinetoplastid flagellates. Trends Parasitol 22: 168–174.
- 41. Rellos P, Ivins FJ, Baxter JE, Pike A, Nott TJ, et al. (2007) Structure and regulation of the human Nek2 centrosomal kinase. J Biol Chem 282: 6833–6842.
- 42. Cloutier M, Vigneault F, Lachance D, Seguin A (2005) Characterization of a poplar NIMA-related kinase PNek1 and its potential role in meristematic activity. FEBS Lett 579: 4659–4665.
- 43. Zhang H, Scofield G, Fobert P, Doonan JH (1996) A nimA-like protein kinase transcript is highly expressed in meristems of Antirrhinum majus. J Microsc 181: 186–194.
- 44. Lynch M, Katju V (2004) The altered evolutionary trajectories of gene duplicates. Trends Genet 20: 544–549.
- 45. Miller SL, DeMaria JE, Freier DO, Riegel AM, Clevenger CV (2005) Novel association of Vav2 and Nek3 modulates signaling through the human prolactin receptor. Mol Endocrinol 19: 939–949.
- 46. Miller SL, Antico G, Raghunath PN, Tomaszewski JE, Clevenger CV (2007) Nek3 kinase regulates prolactin-mediated cytoskeletal reorganization and motility of breast cancer cells. Oncogene.
- 47. Chepurnov VA, Mann DG, Sabbe K, Vyverman W (2004) Experimental studies on sexual reproduction in diatoms. Int Rev Cytol 237: 91–154.
- 48. Nohynkova E, Tumova P, Kulda J (2006) Cell division of Giardia intestinalis: flagellar developmental cycle involves transformation and exchange of flagella between mastigonts of a diplomonad cell. Eukaryot Cell 5: 753–761.
- 49. Cavalier-Smith T (2002) The phagotrophic origin of eukaryotes and phylogenetic classification of Protozoa. Int J Syst Evol Microbiol 52: 297–354.
- 50. Inglis PN, Boroevich KA, Leroux MR (2006) Piecing together a ciliome. Trends Genet 22: 491–500.
- 51. Jekely G, Arendt D (2006) Evolution of intraflagellar transport from coated vesicles and autogenous origin of the eukaryotic cilium. Bioessays 28: 191–198.
- 52. Erickson HP (2007) Evolution of the cytoskeleton. Bioessays 29: 668–677.
- 53. Tucker RW, Pardee AB, Fujiwara K (1979) Centriole ciliation is related to quiescence and DNA synthesis in 3T3 cells. Cell 17: 527–535.
- 54. Qin H, Wang Z, Diener D, Rosenbaum J (2007) Intraflagellar transport protein 27 is a small G protein involved in cell-cycle control. Curr Biol 17: 193–202.
- 55. Keeling PJ, Burger G, Durnford DG, Lang BF, Lee RW, et al. (2005) The tree of eukaryotes. Trends Ecol Evol 20: 670–676.