Email updates

Keep up to date with the latest news and content from BMC Plant Biology and BioMed Central.

Open Access Research article

Conservation of ciliary proteins in plants with no cilia

Matthew E Hodges13, Bill Wickstead23, Keith Gull3 and Jane A Langdale1*

Author Affiliations

1 Department of Plant Sciences, University of Oxford, South Parks Rd, Oxford OX1 3RB, UK

2 Centre for Genetics and Genomics, University of Nottingham, Nottingham NG7 2UH, UK

3 Sir William Dunn School of Pathology, University of Oxford, South Parks Rd, Oxford OX1 3RE, UK

For all author emails, please log on.

BMC Plant Biology 2011, 11:185  doi:10.1186/1471-2229-11-185


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2229/11/185


Received:24 August 2011
Accepted:30 December 2011
Published:30 December 2011

© 2011 Hodges et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Eukaryotic cilia are complex, highly conserved microtubule-based organelles with a broad phylogenetic distribution. Cilia were present in the last eukaryotic common ancestor and many proteins involved in cilia function have been conserved through eukaryotic diversification. However, cilia have also been lost multiple times in different lineages, with at least two losses occurring within the land plants. Whereas all non-seed plants produce cilia for motility of male gametes, some gymnosperms and all angiosperms lack cilia. During these evolutionary losses, proteins with ancestral ciliary functions may be lost or co-opted into different functions.

Results

Here we identify a core set of proteins with an inferred ciliary function that are conserved in ciliated eukaryotic species. We interrogate this genomic dataset to identify proteins with a predicted ancestral ciliary role that have been maintained in non-ciliated land plants. In support of our prediction, we demonstrate that several of these proteins have a flagellar localisation in protozoan trypanosomes. The phylogenetic distribution of these genes within the land plants indicates evolutionary scenarios of either sub- or neo-functionalisation and expression data analysis shows that these genes are highly expressed in Arabidopsis thaliana pollen cells.

Conclusions

A large number of proteins possess a phylogenetic ciliary profile indicative of ciliary function. Remarkably, many genes with an ancestral ciliary role are maintained in non-ciliated land plants. These proteins have been co-opted to perform novel functions, most likely before the loss of cilia, some of which appear related to the formation of the male gametes.

Keywords:
Evolution; Cilia; Flagella; Basal Body; Centriole; Land Plants

Background

Centrioles and cilia/flagella are microtubule-based organelles involved in cellular motility and signalling (for review see [1]). The ultrastructural morphology of these ancient organelles is remarkably conserved in extant eukaryotes and their phylogenetic distribution pattern suggests that the last eukaryotic common ancestor (LECA) possessed the ability to produce a centriole plus cilia with both sensory and motility functions [2-4]. Despite this widespread phylogenetic distribution, lineage specific modifications have been shown to occur and numerous instances of independent cilia loss have been reported [5-9].

In recent years, multiple high-throughput proteomic studies and bioinformatic analyses in disparate species have identified a cohort of proteins associated with centriole and ciliary functions [2,3,10-16]. These proteins include both those with ciliary-specific roles (such as intraflagellar transport proteins, outer- and inner-dynein arms and radial spoke proteins) as well as those such as tubulins that are also involved in other cellular functions (for an extensive database see Cildb [17]). To date, however, little is known about what happens to genes involved in ciliary function when an evolutionary transition to cilia loss occurs.

The land plant lineage is a good model for studying the transition to ciliary loss for several reasons. First, this ancient monophyletic group was ancestrally ciliated, but there have been at least two independent loss events within the group, once in gymnosperms and once at the base of the angiosperms [18-20]. Second, sufficient genomic information exists for an in-depth analysis of protein compositional changes during the process of ciliary loss. Third, the land plants are a sister lineage to the Chlorophyta, a group that includes the well-studied ciliary model species Chlamydomonas reinhardtii. The Chlorophytes thus provide a good outgroup for identification and analysis of genes with ciliary function.

Within the ciliated basal land plants, cilia are produced only in specialised sperm cells [4]. Ciliogenesis in these cells occurs de novo, as opposed to the canonical template pathway seen in animal lineages [21]. It can thus be assumed that there are regulatory mechanisms that ensure correct spatial and temporal expression of genes required for ciliary function within this restricted phase of the plant life cycle. On losing the ability to produce cilia, it is unclear whether these regulatory mechanisms are also lost.

Here, we describe an approach that uses phylogenetic profiling to identify a core set of ciliary function proteins. By applying a scoring system to orthologous sets of proteins within 45 eukaryotic species we aimed to identify proteins with an ancestral ciliary function that are maintained in non-ciliated land plants. We then tested this ancestral ciliary function in the flagellated protozoan Trypanosoma brucei, and investigated potential plant roles through analysis of gene expression data and phylogenetic inference. We conclude that at least some of the genes identified in cilia formation in non-seed plants retain a role in male gametophyte development in seed plants.

Results

Identification of ciliary proteins by phylogenetic distribution

To gain a deeper understanding of the key proteins involved in ciliary function across eukaryotes, we analysed the distribution of orthologuos proteins encoded in the genomes of 45 species (Additional file 1). The species were selected to represent a wide evolutionary spread across six major groups of eukaryotes (namely: Archaeplastida, Excavata, Chromalveolata, Holozoa, Fungi and Amoebozoa). For each species, a complete or near-complete genome sequence is publically available. The species analysed include organisms that produce cilia/flagella at some point in their life cycle (which we refer to here as "ciliated species") and those that do not (referred to here as "non-ciliated species"). Proteins which are phylogenetically conserved in ciliated, but not in non-ciliated species, are considered to possess a "ciliary profile" and are predicted to perform a ciliary function.

Additional file 1. Source and version of the genomes used in the analysis.

Format: PDF Size: 264KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

To identify proteins with a ciliary profile, we performed a reciprocal best BLAST (RBB) analysis using the predicted proteome of Chlamydomonas reinhardtii to identify putative orthologues encoded in the genomes of 44 other eukaryotic species. The green alga Chlamydomonas reinhardtii was chosen as the query genome due to its phylogenetic position in the sister lineage to (at the base of) the land plants. The relatively short phylogenetic distance to the land plants, when compared to the spread of the eukaryotic phylogeny, reduces the chance that sequence divergence would prevent accurate identification of orthologues. Furthermore, Chlamydomonas is a well-studied ciliary model species.

The predicted proteome of Chlamydomonas has 15,143 proteins. Of these, 6,982 proteins have more than one putative orthologue. To reduce the noise of false positives in the dataset, only proteins with identifiable hits in at least 5 species within the analysis were retained. This cut-off resulted in a final defined dataset of 4,802 proteins (31.7% of the proteins in the Chlamydomonas predicted proteome) (Additional file 2). Within this dataset, 772 proteins (16.1%) have a putative orthologue in each of the major eukaryotic groups used in this analysis. The most parsimonious explanation for the phylogenetic distribution of these 772 proteins is that these proteins were present in the last eukaryotic common ancestor (LECA), and that they have subsequently been lost in certain lineages through evolutionary divergence.

Additional file 2. Protein information and sequence identification for orthologous sets of each protein identified in the bioinformatic screen. Protein information shown is the annotation from the JGI database for Chlamydomonas reinhardtii. Apime = Apis mellifera, Cioin = Ciona intestinalis, Caeel = Caenorhabditis elegans, Cyame = Cyanidioschyzon merolae, Capsp = Capitella sp., Dicdi = Dictyostelium discoideum, Drome = Drosophila melanogaster, Danre = Danio rerio, Enccu = Encephalitozoon cunculi, Galga = Gallus gallus, Homsa = Homo sapiens, Lotgi = Lottia gigantea, Leima = Leishmania major, Monbr = Monosiga brevicollis, Neucr = Neurospora crassa, Naegr = Naegleria gruberi, Nemve = Nematostella vectensis, Ostta = Ostreococcus tauri, Parte = Paramecium tetraurelia, Sacce = Saccharomyces cerevisiae, Schpo = Schizosaccharomyces pombe, Strpu = Strongylocentrotus purparatus, Triad = Trichoplax adhaerens, Trybr = Trypanosoma brucei, Tetth = Tetrahymena thermophila, Takru = Takifugu rubripes, Ustma = Ustilago maydis.

Format: PDF Size: 4.7MB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

We found that a simple presence or absence analysis of protein phylogenetic distribution in ciliated and non-ciliated species was too stringent to identify proteins with ciliary function because of false negatives. Therefore, we allowed for both false positives and negatives in the dataset by applying a positive score for a hit in a ciliated species, and a negative score for a hit in a non-ciliated species (see Methods). By using this scoring system we found that out of the 4,802 proteins in the dataset, 213 (4.4%, Additional file 3) proteins possess a score indicative of a ciliary profile (examples in Figure 1). As an indication that the method to identify ciliary proteins is accurate, 157 of these 213 (73.7%) proteins have been previously detected within cilia in proteomic analyses. This equates to 33.3% of all proteins (472) found in proteomic analyses of cilia. The remaining 56 proteins with a ciliary profile currently have no known function. We hypothesise that these represent a novel set of proteins with an as yet uncharacterised ciliary function.

Additional file 3. Protein information and sequence identification for ciliary profile proteins.

Format: PDF Size: 2MB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

thumbnailFigure 1. A) Phylogenetic distribution of representative examples of the filtered ciliary profile proteins from the reciprocal best BLAST analysis with different scores. Ciliated species are shown in blue. B) Venn diagram depicting the results of the bioinformatic analysis. A total of 4,802 proteins have five or more identifiable orthologues in 45 queried genomes. 472 of these proteins are detected in ciliary proteomic screens (green). 213 proteins are identified with the ciliary profile (red), of which 157 are detected in proteomic screens (yellow). 21 proteins have been identified as

    c
onserved in
    c
iliated species and in the land
    p
lants (CCP, blue), of which one has been identified in proteomic screens. C) Western blot of whole cell lysates with anti-TY BB2 antibody (above) and loading control anti-PFR2 antibody (below). D) Direct fluorescence of whole cell Trypanosoma brucei showing cellular localisation of TbFBB17 (Tb10.406.0130) N-terminal YPF chimera (left panel) and an overlay of signal from YFP (green), DAPI (magenta) and phase-contrast (right panel). Tb10.406.0130 localises to the flagellum (white arrow). Scale bar = 2 μm.

Ciliary profile proteins localise to the flagellum of trypanosomes

Most of the ciliary profile proteins have been detected in one or more proteomic studies of cilia [10,12,14,22]. To test whether the 56 proteins which have not been previously identified in proteomic studies are likely to play a role in cilia function, we investigated the cellular localisation of one of these proteins in the single-celled excavate Trypanosoma brucei. T. brucei is a notable model organism for cilia studies, with many examples of transgene introduction for localisation studies. Many of the ciliary profile proteins with known function have already been localised in T. brucei (for review see [23]).

We selected a candidate protein, FBB17, which has a putative trypanosome orthologue (herein called TbFBB17) encoded by the gene Tb10.406.0130, for endogenous locus TY-YFP-tagging in T. brucei. TbFBB17 is a putative zinc carboxypeptidase which has not been detected in any cilia proteomic study to date [10,12,14,17,22]. The tagging strategy ensured that the YFP-TbFBB17 chimera is transcribed by the endogenous gene promoter. This strategy has previously been shown to produce protein levels similar to wild-type for other loci [24]. The correct incorporation of the YFP-tag at the 5' end of the targeted Tb10.406.0130 gene was checked by PCR (data not shown) and by immunoblotting of cell lysates with monoclonal antibodies directed against the short epitope TY-tag [25](Figure 1).

To determine the intracellular localisation of TbFBB17, transformed cells were examined by fluorescence microscopy. Figure 1 shows that TbFBB17 is found along the length of the flagellum (white arrow). This observation supports the prediction that the identified ciliary profile proteins do indeed play a role in ciliary function. To determine whether TbFBB17 is stably associated with the cytoskeleton, we assessed the localisation of the tagged protein in cells which had been detergent-extracted to remove non-cytoskeletal components. In such cells the TbFBB17 flagellum signal is removed (data not shown). This explains the lack of identification of TbFBB17 in the trypanosome cilia proteome as the proteomes are derived from cytoskeletal cell fractions [5]. A similar explanation may apply to the other 55 proteins that have been identified in this study but not in ciliary proteomic screens.

Identification of ciliary proteins in non-ciliated land plants and fungi

In land plants cilia have been lost at least twice, once in the gymnosperms and once at the base of the angiosperms [18-20]. Such evolutionary transitions pose the question of what happens to ancestral ciliary proteins when the ability to produce cilia is lost. To investigate this phenomenon, a scoring system was defined that identified proteins with putative orthologues conserved in both ciliated species and also land plants, but not in other non-ciliated species (CCP, for

    c
onserved in
    c
iliated and
    p
lants, see Methods). This analysis identified 21 CCP proteins, out of the dataset of 4,802 proteins (0.44%), that have been retained in land plants which do not possess the ability to produce a cilium (Figure 2, Additional file 4). In agreement with our predictions, one of these proteins (BUG22) has previously been detected in ciliary proteomic studies [12,26]. The remaining proteins have no characterised ciliary function. A similar analysis within the non-ciliated fungal lineages (present in ciliated species and non-ciliated fungi) identifies two proteins (0.04% of the 4,802 proteins in the dataset). Neither of these two proteins are present in the CCP set, and neither have been detected in any proteomic study. The lower number of proteins identified in the fungal analysis may reflect the high level of genomic streamlining in this lineage [27].

thumbnailFigure 2. A) Phylogenetic distribution of the filtered conserved in ciliated species and non-ciliated land plants, but absent in other non-ciliated species (CCP) proteins from the reciprocal best BLAST analysis. Ciliated species are shown in blue and non-ciliated land plants in green. B) Western blot of whole cell lysates with anti-TY antibody (above) and loading control anti-PFR2 antibody (below). C) Direct fluorescence of whole cell Trypanosoma brucei showing cellular localisation of TbCCP1 (Tb11.02.1370), TbBUG22 (Tb10.70.5560) and TbCCP2 (Tb927.8.5380). Left panel showing fluorescence signal from YFP and an overlay in right panel of signal from YFP (green), DAPI (magenta) and phase-contrast. All three possess a ciliary localisation (white arrows point to flagella (CCP1 and BUG22) and basal bodies (CCP2)). D) Direct fluorescence of whole cell Trypanosoma brucei showing cellular localisation of TbCCP2 (Tb927.8.5380). Expression is cell-cycle regulated (white arrows depict basal bodies). Scale bar = 2 μm.

Additional file 4. Protein information and sequence identification for CCP proteins.

Format: PDF Size: 584KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

CCP proteins localise to the trypanosome flagellum

To test the hypothesis that CCP proteins possess an ancestral ciliary function, localisation of three proteins from this subset was examined in T. brucei cells. These three proteins, named TbCCP1 (the putative p60 subunit of katanin), TbBUG22 (the trypanosome orthologue of BUG22), and TbCCP2 (a proteins of hypothetical function) are encoded by genes Tb11.02.1370, Tb10.70.5560 and Tb927.8.5380, respectively. Orthologues of CCP1 and CCP2 have not been detected in any previous ciliary proteomic study, whereas BUG22 has been identified in a proteomic analysis of Chlamydomonas cilia [12]. The localisation of each protein was investigated using endogenous locus N-terminal TY-YFP tagging as for TbFBB17. Correct insertion of the tag was checked by PCR (data not shown) and immunoblotting of cell lysates with monoclonal antibodies directed against the TY tag (Figure 2).

Fluorescence microscopy of transformed cells demonstrated that all three proteins have a localisation pattern consistent with a role in ciliary function. Both TbCCP1 and TbBUG22 localise to the flagellum of T. brucei throughout the cell cycle (Figure 2, white arrows). Interestingly, TbCCP2 was found to localise to a position consistent with the basal body of the flagellum in a cell-cycle regulated fashion (Figure 2, white arrows). Trypanosomes undergo morphological changes which can be used to position cells within the cell cycle (see cartoons in Figure 2). Early in the cell cycle, the pro-basal body matures to form a second basal body from which a new flagellum will extend. Shortly thereafter, the kinetoplast (an organelle containing all the mitochondrial DNA) divides, coupled to the separation of the two basal bodies and extension of a new flagellum from the posterior basal body. Once the new flagellum is fully extended, nuclear mitosis occurs and finally cytokinesis. Examination of TbCCP2 fluorescence in the light of this positional information reveals that TbCCP2 is not detected just before kinetoplast division, nor during nuclear division, but it is present on or near the basal bodies throughout the rest of the cell cycle (Figure 2). In conclusion, these data demonstrate that the identified CCP proteins do indeed perform a ciliary function in a species which can produce flagella.

Evolutionary history of CCP proteins

To investigate possible functions of the CCP proteins in non-ciliated land plants we performed phylogenetic analysis of the identified protein sets. Understanding the evolutionary relationships of the CCP proteins provides insight into possible mechanisms of neo- or subfunctionalisation. Bayesian phylogenetic inference of the identified homologous proteins of BUG22 (Figure 3) suggest that the gene duplicated at the base of the land plants because Physcomitrella patens a basal land plant with the ability to produce cilia, possesses four copies of the gene that resolve in 2 clades. Subfunctionalisation may then have occurred with one copy retaining the ancestral ciliary role and the other being co-opted for alternative roles. Subsequently the ancestral copy was lost within the non-ciliated land plants (represented by Arabidopsis thaliana, Oryza sativa and Populus trichocarpa). In contrast, the Bayesian phylogenies of the other CCP orthologues are consistent with the hypothesis of neofunctionalisation/co-option of an ancestral ciliary protein upon the evolutionary transition to loss of ability to produce cilia because the genes resolve into a single land plant clade (Figure 3).

thumbnailFigure 3. A) Phylogenetic Bayesian inferences for four of the CCP proteins, demonstrating different evolutionary scenarios. Land plants in green. For each protein, Pfam domain structure of human orthologue is shown. B) Heatmap illustration of microarray expression data for the CCP proteins (left) and controls (right) in different developmental stages scored by fold change from median.

CCP proteins are highly expressed in male floral organs and mature pollen

While the CCP proteins have a ciliary role in ciliated species, the question of what their function is within the non-ciliated land plants remains. To this end, we analysed microarray data mined from a public database of the non-ciliated flowering plant Arabidopsis thaliana [28]. The expression pattern of each of the genes at different developmental stages is shown in Figure 3, (expression data in Additional file 5). Clustering of the CCP expression data by fold change shows that the genes are highly expressed in mature pollen, with an average fold increase of 2.8. In contrast, the expression pattern of a number of control genes shows a more ubiquitous expression pattern. We performed a Monte Carlo estimation analysis to calculate the sum fold change of pollen expression data for 22,640 random sets of 21 genes within the Arabidopsis genome to determine how likely such a fold increase in expression is. This showed that pollen expression levels of the CCP proteins was very highly significant (p = 0.0006).

Additional file 5. mRNA Expression data for Arabidopsis thaliana mined from and used in analysis of CCP protein functions.

Format: PDF Size: 176KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Discussion

Over 300 proteins have previously been associated with ciliary function (for example see [17]) with most being identified through proteomic analysis of isolated cilia (e.g. [10,12,14,15]). However, such studies fail to detect proteins which are not stably associated with the protein preparation used. Here, we have identified a novel set of 56 proteins possessing a ciliary profile, through phylogenetic distribution analysis and experimental validation in T. brucei (Figure 1). In addition we identified 21 proteins with a predicted ancestral ciliary role that have been maintained in non-ciliated land plants (Figures 2, 3).

Flowering plants have retained CCP proteins despite loss of ability to produce cilia

The phylogenetic distribution of cilia in land plants shows that cilia loss occurred at least twice: once in gymnosperms at the divergence of Ginkgo and the Cycads from the Conifers and Gnetals and once at the base of the flowering plants [18,19]. All other land plants produce cilia only within specialised sperm cells ([29-32], reviewed in [20]). During the evolutionary transition to loss of cilia, many ciliary function proteins were also lost (Additional file 3). However, we have identified 21 proteins with identifiable orthologues that are

    c
onserved in
    c
iliated species and in the land
    p
lants (CCP proteins). At least a subset of these proteins possess a ciliary function in the flagellated Trypanosoma brucei, but the function of orthologues within the non-ciliated land plants has presumably diverged. Interestingly, a similar analysis for those proteins conserved during the evolutionary loss of cilia in fungal lineages identifies only two proteins, neither of which is present in the CCP set. Such low (21 and 2) numbers suggest that while loss of proteins is widespread during evolutionary transitions, it is relatively rare for a protein to be co-opted to a different function without undergoing significant sequence divergence.

The flagella localisation of all three tested CCP proteins in trypanosomes (Figure 2) supports their hypothesised ancestral ciliary role and validates the scoring system that we used for CCP identification. While there is limited data on the putative roles of the 21 CCP proteins in non-ciliated plants, it is of note that where information is available they are generally associated with cytoskeletal/microtubule roles. For example, CCP1 is a putative katanin, a class of known AAA+ microtubule-severing proteins [33] and when mutated the CCP gene tonneau2/fass1 markedly affects cell shape and arrangement in Arabidopsis thaliana, due to cytoskeletal abnormalities [34]. Another CCP protein, FUSED, has been shown to be involved in the hedgehog signalling pathway in Drosophila, and hedgehog receptors localise to the cilium in many species [35]. Together, these observations suggest that CCP proteins play regulatory roles that were ancestrally associated with cilia, but are now maintained for cytoskeletal functions in non-ciliated plants. Examination of the CCP proteins for Pfam domains [36] which include ubiquitanases, AAA+ domain containing proteins, protein kinases, tubulin-tryosine ligases and phosphatidic acid phosphatases, supports the suggestion that ancestral ciliary function proteins have been co-opted in land plants to perform cytoskeletal roles. Picket-Heaps [37] proposed an autogenous origin of cilia from a microtubule-organising centre ancestrally involved in the regulation of cytoplasmic and mitotic microtubules [4]. If this was the case, it is equally possible that upon loss of cilia evolutionary pressures drove the co-option of cytoskeletal functions in the opposite direction.

CCP protein function can also be inferred from the observation that CCP genes are highly expressed within mature pollen [28] (Figure 3, Additional file 5). This expression pattern is curious given the ancestral role of these proteins in ciliary formation in the sperm cells of basal land plants and the fact that sperm cells in pollen are not ciliated. It is conceivable that genes involved in the cell-specific formation of cilia in basal land plants may have acquired additional (new) functions in gametes before the loss of cilia, and are thus maintained in pollen on loss of cilia to perform similar cell-specific (possibly cytoskeletal) roles in a non-ciliary context. Such specificity would explain why only a small number of ciliary function genes were maintained or co-opted on the evolutionary transition to loss of cilia in seed plants, whereas the majority were lost or diverged.

Conclusions

The work presented here demonstrates the existence of a set of ancestral ciliary proteins within non-ciliated land plants, and suggests that the process of cilia loss can be accompanied by the co-option of ciliary function proteins into novel contexts.

Methods

Bioinformatic definition of protein sets

Every predicted protein encoded in the Chlamydomonas genome was used as a query in a reciprocal best BLAST (RBB--see below) search with an e-value threshold of < 10-5 against each of the predicted proteomes of each of 45 eukaryotic organisms. During RBB, a BLAST search is performed using a query sequence against a target genome. If the top hit from the target genome identifies the query sequence in the reciprocal BLAST, the sequence identified is a RBB hit. The species used in the analysis have a broad phylogenetic distribution and include 29 organisms which produce centriole/cilia/flagella at some point in their life-cycle (termed "ciliated") and 16 organisms which do not (termed "non-ciliated"). Importantly, each of these species possess a complete or near-complete genome sequence which is publically available. If the top hit from the RBB was the initial query sequence from Chlamydomonas, the identified target sequence was inferred to be an orthologue and was thus included in the dataset. Non-reciprocating sequences were discarded. The ciliary profile proteins were determined as those proteins whose evolutionary distribution of orthologues within the 45 species scored above 10 for a metric of β - 10α, (where α was the number of non-ciliated species possessing an orthologue and β was the number of ciliated species possessing an orthologue). The conserved in ciliated species and plants (CCP) protein set was defined as proteins with plant orthologues where β > 18 and μ < 2, (where μ was the number of non-ciliated species, excluding plant species, possessing an orthologue).

Phylogenetic inference

For phylogenetic analysis, the sequences of each protein set were aligned using MAFFT6.24 [38] adopting the L-INS-I strategy, and trimmed by eye to well aligned blocks. A Bayesian phylogeny was inferred from the protein alignment using the metropolis-coupled Markov chain Monte Carlo (MCMCMC) method as implemented in the program MrBayes3.1.2, with 2 independent runs from random start trees using the WAG substitution matrix with a gamma-distributed variation in substitution rate approximated to 4 discrete categories [39].

Endogenous-locus tagging

Endogenous-locus tagging was performed using procyclic-form Trypanosoma brucei Lister 427 cells [40]. Cells were maintained as axenic cultures in SDM79 supplemented with 10% v/v foetal calf serum (FCS) at 28°C, as described by Brun and Schӧnenberger [41]. Sequences were integrated at endogenous gene loci following the strategy outlined in [42]. Briefly, genomic DNA for PCR was isolated from strain Lister 427 T. brucei procyclic cells as in [43]. PCR amplicons containing ~200 bp from the 5'end of the target cDNA sequence (CDS) with a linearisation site added to the 5' end of the 3' primer, and ~200 from the 3'end of the 5' intergenic region for the target CDS with the same linearization site incorporated into the 5' primer for this fragment, were cloned together into the XbaI-BamHI sites downstream of TY-YFP in the pEnT6P vector, such that the N-terminal end of the CDS was in frame with YFP. The vector was then linearised at the site between the 5'end of the CDS and intergenic region fragments prior to transfection into procyclic-form Lister 427 T. brucei cells as described previously [44]. Post-electroporation, cells were allowed to recover in normal growth media for 14 h, after which time a stable transformant population was selected by the addition of 1 μg ml-1 puromycin.

Fluorescence microscopy

For analysis of YFP-tagged proteins in fixed whole cells, procyclic trypanosomes were harvested from mid-log phase cultures by centrifugation, washed once in phosphate buffered saline (PBS) and allowed to adhere to ethanol-washed plain glass slides for 5 min in PBS at ~2 × 107 cells ml-1. Cells were then fixed for 5 min in 2% w/v formaldehyde in PBS, followed by permeabilisation in -20°C methanol for 10 min and rehydrated in PBS. Samples were mounted in 1% w/v 1,4-Diazabicyclo[2.2.2]octane, 90% glycerol, 50 mM sodium phosphate pH 8.0 containing 0.25 μg ml-1 4',6-diamidino-2-phenylindole.

Western blot analysis

For whole-cell protein samples, cells were harvested, washed in PBS buffer, pelleted and immediately resuspended in boiling Laemmli buffer (2% sodium dodecyl sulphate (SDS), 10% glycerol, 400 mM β-mercaptoethanol, 50 mM Tris-HCl pH 7.2). SDS-polyacrylamide gel electrophoresis and immersion transfer to nitrocellulose membrane were performed using standard techniques described elsewhere [43]. For immunoblotting, the membranes were blocked with 5% (w/v) skimmed milk powder in TBS-T (20 mM Tris-HCl, 500 mM NaCl, 0.05% Tween 20), labelled with anti-TY BB2 monoclonal antibody [25] in TBS-T/1% skimmed milk powder, and developed with horseradish peroxidase-conjugated anti-mouse immunoglobulins (Sigma). As a loading control, membranes were re-probed with anti-PFR2 (L8C4) monoclonal antibody and developed as above.

Abbreviations

LECA: Last Eukaryotic Common Ancestor; BLAST: Basic Local Alignment Search Tool; RBB: Reciprocal Best BLAST; CCP: Conserved in Ciliated and Plants.

Authors' contributions

MEH carried out the bioinformatic analysis and molecular experiments. MEH, BW, JAL and KG conceived the study and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We would like to thank Nicole Scheumman (Dunn School of Pathology, University of Oxford) for expert guidance on the transformation of T. brucei. We are grateful to Robert Scotland and Jim Fouracre (Dept. Plant Sciences, University of Oxford) for useful discussion. We would like to thank Steven Kelly for advice and assistance with informatics and computational resources. This work was funded by grants from the Wellcome Trust and BBSRC (BB/day/1) to KG and the Gatsby Charitable foundation to JAL. MEH is supported by a graduate studentship from BBSRC.

References

  1. Mitchell DR: The evolution of eukaryotic cilia and flagella as motile and sensory organelles.

    Eukaryot Membr Cytoskeleton: Orig Evol 2007, 607:130-140. Publisher Full Text OpenURL

  2. Pereira-Leal JB, Carvalho-Santos Z, Machado P, Branco P, Tavares-Cadete F, Rodrigues-Martins A, Bettencourt-Dias M: Stepwise evolution of the centriole-assembly pathway.

    J Cell Sci 2010, 123:1414-1426. PubMed Abstract | Publisher Full Text OpenURL

  3. Hodges ME, Scheumann N, Wickstead B, Langdale JA, Gull K: Reconstructing the evolutionary history of the centriole from protein components.

    J Cell Sci 2010, 123:1407-1413. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Cavalier-Smith T: Basal body and flagellar development during the vegetative cell cycle and the sexual cycle of Chlamydomonas reinhardii.

    J Cell Sci 1974, 16:529-556. PubMed Abstract | Publisher Full Text OpenURL

  5. Broadhead R, Dawe HR, Farr H, Griffiths S, Hart SR, Portman N, Shaw MK, Ginger ML, Gaskell SJ, McKean PG, Gull K: Flagellar motility is required for the viability of the bloodstream trypanosome.

    Nature 2006, 440:224-227. PubMed Abstract | Publisher Full Text OpenURL

  6. Inglis PN, Ou G, Leroux MR, Scholey JM: The sensory cilia of Caenorhabditis elegans.

    WormBook 2007, 1-22. OpenURL

  7. Satir P, Christensen ST: Overview of structure and function of mammalian cilia.

    Annu Rev Physiol 2007, 69:377-400. PubMed Abstract | Publisher Full Text OpenURL

  8. Wickstead B, Gull K: Dyneins across eukaryotes: a comparative genomic analysis.

    Traffic 2007, 8:1708-1721. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Cavalier-Smith T: The origin of eukaryotic and archaebacterial cells.

    Ann N Y Acad Sci 1987, 503:17-54. PubMed Abstract | Publisher Full Text OpenURL

  10. Pazour GJ, Agrin N, Leszyk J, Witman GB: Proteomic analysis of a eukaryotic cilium.

    J Cell Biol 2005, 170:103-113. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  11. Keller LC, Romijn EP, Zamora IP, Yates JR, Marshall WF: Dissecting the centriole proteome.

    FASEB J 2006, 20:A493-A493. OpenURL

  12. Keller LC, Romijn EP, Zamora I, Yates JR, Marshall WF: Proteomic analysis of isolated Chlamydomonas centrioles reveals orthologs of ciliary-disease genes.

    Curr Biol 2005, 15:1090-1098. PubMed Abstract | Publisher Full Text OpenURL

  13. Kilburn CL, Pearson CG, Romijn EP, Meehl JB, Giddings TH, Culver BP, Yates JR, Winey M: New Tetrahymena basal body protein components identify basal body domain structure (vol 178, pg 905, 2007).

    J Cell Biol 2007, 179:167-167. Publisher Full Text OpenURL

  14. Inglis PN, Boroevich KA, Leroux MR: Piecing together a ciliome.

    Trends Genet 2006, 22:491-500. PubMed Abstract | Publisher Full Text OpenURL

  15. Merchant SS, Prochnik SE, Vallon O, Harris EH, Karpowicz SJ, Witman GB, Terry A, Salamov A, Fritz-Laylin LK, Marechal-Drouard L, et al.: The Chlamydomonas genome reveals the evolution of key animal and plant functions.

    Science 2007, 318:245-251. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Woodland HR, Fry AM: Pix proteins and the evolution of centrioles.

    PloS One 2008, 3:e3778. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Arnaiz O, Malinowska A, Klotz C, Sperling L, Dadlez M, Koll F, Cohen J: Cildb: a knowledgebase for centrosomes and cilia.

    Database (Oxford) 2009, 2009:bap022. Publisher Full Text OpenURL

  18. Finet C, Timme RE, Delwiche CF, Marleta F: Multigene phylogeny of the green lineage reveals the origin and diversification of land plants.

    Curr Biol 2010, 20:2217-2222. PubMed Abstract | Publisher Full Text OpenURL

  19. Bremer K, Humphries CJ, Mishler BD, Churchill SP: On cladistic relationships in green plants.

    Taxon 1987, 36:339-349. Publisher Full Text OpenURL

  20. Renzaglia KS, Garbary DJ, Motile gametes of land plants: Diversity, development, and evolution.

    Crit Rev Plant Sci 2001, 20:107-213. OpenURL

  21. Beisson J, Wright M: Basal body/centriole assembly and continuity.

    Curr Opin Cell Biol 2003, 15:96-104. PubMed Abstract | Publisher Full Text OpenURL

  22. Andersen JS, Wilkinson CJ, Mayor T, Mortensen P, Nigg EA, Mann M: Proteomic characterization of the human centrosome by protein correlation profiling.

    Nature 2003, 426:570-574. PubMed Abstract | Publisher Full Text OpenURL

  23. Bastin P, Pullen TJ, Moreira-Leite FF, Gull K: Inside and outside of the trypanosome flagellum: a multifunctional organelle.

    Microbes Infect 2000, 2:1865-1874. PubMed Abstract | Publisher Full Text OpenURL

  24. Wickstead B, Carrington JT, Gluenz E, Gull K: The expanded kinesin-13 Repertoire of trypanosomes contains only one mitotic kinesin indicating multiple extra-nuclear roles.

    PloS One 2010, 5:e15020. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Bastin P, Bagherzadeh Z, Matthews KR, Gull K: A novel epitope tag system to study protein targeting and organelle biogenesis in Trypanosoma brucei.

    Mol Biochem Parasitol 1996, 77:235-239. PubMed Abstract | Publisher Full Text OpenURL

  26. Laligne C, Klotz C, de Loubresse NG, Lemullois M, Hori M, Laurent FX, Papon JF, Louis B, Cohen J, Koll F: Bug22p, a conserved centrosomal/ciliary protein also present in higher plants, is required for an effective ciliary stroke in Paramecium.

    Eukaryot Cell 2010, 9:645-655. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Soanes DM, Alam I, Cornell M, Wong HM, Hedeler C, Paton NW, Rattray M, Hubbard SJ, Oliver SG, Talbot NJ: Comparative genome analysis of filamentous fungi reveals gene family expansions associated with fungal pathogenesis.

    PloS One 2008, 3:15. OpenURL

  28. Schmid M, Davison TS, Henz SR, Pape UJ, Demar M, Vingron M, Scholkopf B, Weigel D, Lohmann JU: A gene expression map of Arabidopsis thaliana development.

    Nat Genet 2005, 37:501-506. PubMed Abstract | Publisher Full Text OpenURL

  29. Graham LE: Coleochaete and the origin of land plants.

    Am J Bot 1984, 71:603-608. Publisher Full Text OpenURL

  30. Graham LE: The occurrence, evolution, and phylogenetic significance of parenchyma in Coleochaete breb(Chlorophyta).

    Am J Bot 1982, 69:447-454. Publisher Full Text OpenURL

  31. Graham LE, Mcbride GE: Occurrence and phylogenetic significance of a multilayered structure in Coleochaete spermatozoids.

    Am J Bot 1979, 66:887-894. Publisher Full Text OpenURL

  32. Graham LE: Ultrastructure of spermatozoids of Coleochaete scutata.

    J Phycol 1978, 14:30-30. OpenURL

  33. Mcnally FJ, Vale RD: Identification of katanin, an ATPase that severs and disassembles stable microtubules.

    Cell 1993, 75:419-429. PubMed Abstract | Publisher Full Text OpenURL

  34. Camilleri C, Azimzadeh J, Pastuglia M, Bellini C, Grandjean O, Bouchez D: The Arabidopsis TONNEAU2 gene encodes a putative novel protein phosphatase 2A regulatory subunit essential for the control of the cortical cytoskeleton.

    Plant Cell 2002, 14:833-845. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. Oh SA, Johnson A, Smertenko A, Rahman D, Park SK, Hussey PJ, Twell D: A divergent cellular role for the FUSED kinase family in the plant-specific cytokinetic phragmoplast.

    Curr Biol 2005, 15:2107-2111. PubMed Abstract | Publisher Full Text OpenURL

  36. Finn RD, Mistry J, Tate J, Coggill P, Heger A, Pollington JE, Gavin OL, Gunasekaran P, Ceric G, Forslund K, et al.: The Pfam protein families database.

    Nucleic Acids Res 2010, 38:D211-D222. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  37. Pickett-Heaps J: Autonomy of the centriole--fact or fallacy?

    Cytobios 1971, 3:205-214. OpenURL

  38. Katoh K, Toh H: Recent developments in the MAFFT multiple sequence alignment program.

    Brief Bioinform 2008, 9:286-298. PubMed Abstract | Publisher Full Text OpenURL

  39. Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models.

    Bioinformatics 2003, 19:1572-1574. PubMed Abstract | Publisher Full Text OpenURL

  40. Cross GAM: Identification, purification and properties of clone-specific glycoprotein antigens constituting surface coat of Trypanosoma brucei.

    Parasitology 1975, 71:393-417. PubMed Abstract | Publisher Full Text OpenURL

  41. Brun R, Schonenberger M: Cultivation and invitro cloning of procyclic culture forms of Trypanosoma brucei in a semi-defined medium.

    Acta Trop 1979, 36:289-292. PubMed Abstract OpenURL

  42. Kelly S, Reed J, Kramer S, Ellis L, Webb H, Sunter J, Salje J, Marinsek N, Gull K, Wickstead B, Carrington M, Functional genomics in Trypanosoma brucei: A collection of vectors for the expression of tagged proteins from endogenous and ectopic gene loci.

    Mol Biochem Parasitol 2007, 154:103-109. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  43. Chong L: Molecular cloning--A laboratory manual.

    Science 3rd edition. 2001, 292:446-446. Publisher Full Text OpenURL

  44. Wirtz E, Leal S, Ochatt C, Cross GAM: A tightly regulated inducible expression system for conditional gene knock-outs and dominant-negative genetics in Trypanosoma brucei.

    Mol Biochem Parasitol 1999, 99:89-101. PubMed Abstract | Publisher Full Text OpenURL