- Open Access
Basal body positioning and anchoring in the multiciliated cell Paramecium tetraurelia: roles of OFD1 and VFL3
Ciliavolume 6, Article number: 6 (2017)
The development of a ciliary axoneme requires the correct docking of the basal body at cytoplasmic vesicles or plasma membrane. In the multiciliated cell Paramecium, three conserved proteins, FOR20, Centrin 2, and Centrin 3 participate in this process, FOR20 and Centrin 2 being involved in the assembly of the transition zone. We investigated the function of two other evolutionary conserved proteins, OFD1 and VFL3, likely involved in this process.
In Paramecium tetraurelia, a single gene encodes OFD1, while four genes encode four isoforms of VFL3, grouped into two families, VFL3-A and VFL3-B. Depletion of OFD1 and the sole VFL3-A family impairs basal body docking. Loss of OFD1 yields a defective assembly of the basal body distal part. Like FOR20, OFD1 is recruited early during basal body assembly and localizes at the transition zone between axoneme and membrane at the level of the microtubule doublets. While the recruitment of OFD1 and Centrin 2 proceed independently, the localizations of OFD1 and FOR20 at the basal body are interdependent. In contrast, in VFL3-A depleted cells, the unanchored basal bodies harbor a fully organized distal part but display an abnormal distribution of their associated rootlets which mark their rotational asymmetry. VFL3-A, which is required for the recruitment of Centrin 3, is transiently present near the basal bodies at an early step of their duplication. VFL3-A localizes at the junction between the striated rootlet and the basal body.
Our results demonstrate the conserved role of OFD1 in the anchoring mechanisms of motile cilia and establish its relations with FOR20 and Centrin 2. They support the hypothesis of its association with microtubule doublets. They suggest that the primary defect of VFL3 depletion is a loss of the rotational asymmetry of the basal body which specifies the sites of assembly of the appendages which guide the movement of basal bodies toward the cell surface. The localization of VFL3 outside of the basal body suggests that extrinsic factors could control this asymmetry.
Generally absent in land plants or higher fungi, one or several cilia involved in motility, feeding, sensation, and sexual process, localize at the cell surface of numerous protists. This diversity is also observed in metazoa in which the number of cilia protruding from the surface depends on the cellular type: either a solitary cilium, also called primary cilium or several, often hundred cilia, can develop at the surface of specialized differentiated cells. Cilia are anchored at the cell surface by the basal body. In the case of the primary cilium, it is a derivative of the mother centriole, one of the two centrioles which organize the centrosome and contains two specific structures, the distal and sub-distal appendages.
Depending on the cellular type, the development of the primary cilium is initiated in the cytoplasm after docking of the mother centriole, via its distal appendages to Golgi derived vesicles or, directly in the external cellular space, after docking at the plasma membrane [1, 2]. This second pathway is also observed during the formation of motile cilia in multiciliated cells which occurs after multiplication of the centrioles/basal bodies deep in the cytoplasm and their migration towards the apical membrane. In multiciliated cells including epithelial cells of the human oviduct , of the epidermal Xenopus embryo  and of the mouse trachea , observations of basal bodies whose distal parts were surrounded by large or small vesicles suggest that in some cases, a basal body/vesicle interaction could also precede the docking of the basal body at the cell surface. The key role of the distal appendages in the anchoring process is supported by several results demonstrating that the depletion of proteins involved in their assembly impairs the docking .
In the multiciliated organism Paramecium, basal body anchoring at the cell surface does not involve vesicular intermediates . The ciliogenesis is initiated by the development, from an already anchored basal body, of a new basal body which directly docks at the surface. Ciliated and non-ciliated basal bodies are observed indicating that axoneme extension is not necessarily associated with the basal body docking event . Three typical plates organize the transition zone of the ciliated basal body, which also displays transition fibers, a ciliary necklace and Y links . These three plates appear more closely apposed at the distal end of non-ciliated basal body and form the pro-transition zone . These structures cap the tip of the basal body before its docking at the cell surface .
Three appendages protrude asymmetrically from the Paramecium basal bodies, one striated rootlet and two microtubular ribbons. They are thought to act as a scaffold for defining the site of new basal body assembly and to maintain the organization of the basal body row at the cell surface [11, 12]. A transient appendage of the mother basal body, the anterior left filament (ALF) also guides the movement of the daughter basal body toward the cell surface .
We have previously shown that Centrin 2, Centrin 3, and FOR20 are involved in the positioning and anchoring of basal bodies at the surface [8, 14]. Centrin 2 and FOR20 localize, respectively, in the basal body lumen and at the distal part of anchored basal body. The assembly of the structural elements of the pro-transition zone is impaired upon Centrin 2 and FOR20 depletion. Thus, in Paramecium as in metazoan cells, there is a correlation between defects in the structure of the basal body tip and in the anchoring process. In contrast, the depletion of Centrin 3 which is mainly located anteriorly to the proximal part of the centriolar cylinder  does not lead to defects in this assembly. The fully assembled new basal bodies are disoriented with respect to the mother, suggesting that the initial defect resides in the movement of the new basal body towards its docking site. This is supported by the fact that inactivation of Centrin 3 inhibits the formation of the ALF which may contribute to tilt up the basal body toward the cell surface .
In order to go further in the anchoring process in Paramecium, we undertook a functional analysis of homologs of conserved protein likely involved in this mechanism such as OFD1 and VFL3. OFD1, a protein whose mutations induce ciliopathies such as oro-facial-digital syndrome type 1  is well known to be involved in the basal body anchoring process during the formation of the primary cilium. Located at the distal end of the mother centriole in mammalian cells , it is essential for the formation of the distal appendages which anchor the ciliary membrane vesicles before axoneme extension. Like FOR20, OFD1 possesses a conserved TOF/LisH motif  required for centrosomal localization . OFD1 was recently shown to form a ternary complex with FOR20 and OFIP at the centriolar satellites of mammalian cells . However the function of OFD1 in the assembly of the motile cilium remains relatively unstudied. VFL3 was thoroughly studied in Chlamydomonas. Its mutation which affects basal body positioning at the cell surface is correlated with the lack of appendages . Its ortholog in metazoa, CCDC61, is associated with the centrosome in human cells [21, 22]. We present here the analysis of OFD1 and VFL3 in Paramecium.
Strains and culture conditions
Stock d4-2 of Paramecium tetraurelia, the wild-type reference strain, was used in RNAi experiments. The mutant nd7-1  which carries a recessive monogenic mutation preventing trichocyst discharge, a dispensable function under laboratory conditions, was used for transformations. Cells were grown at 27 °C in a wheat grass powder infusion, (Bio Herbe de Blé L’arbre de vie, Luçay Le Male, France) bacterized with Klebsiella pneumoniae and supplemented with 0.8 µg/ml β-sitosterol according to standard procedures .
Constitutive expression of the Myc-tagged VFL3-1 and -VFL3-3 genes: the tag was added at the 5′ end of each gene by PCR amplification, using 5′-end specific primers in which the Myc coding sequence was added. After restriction digests, the fragments were cloned into the SpeI restriction site of the pPXV vector, the recombinant genes being under the control of the Paramecium calmodulin constitutive regulators.
Expression of the GFP-VFL3-3 and GFP-OFD1 under the control of their own regulators: the putative promoter of each gene was first amplified by PCR and cloned at the 3′ end of the GFP-coding fragment inserted into the pZZ-GFP vector (a derivative of the pPXV-GFP vector kindly provided by Jean Cohen) using the BamH1 and SpeI restriction sites. Then, each gene, followed by its putative 3′ UTR, was amplified by PCR and cloned at the 3′ end of the GFP sequence, between the Xho1 and XmaI sites.
Constitutive expression of GFP-OFD1: the OFD1 gene was amplified by PCR from genomic DNA and cloned downstream the GFP sequence into the KpnI site of the pPXV-GFP plasmid which contains the constitutive regulators of the calmodulin gene. Three glycine codons were added between the GFP and the OFD1 sequences.
Gene silencing: all RNAi plasmids are derivatives of the vector L4440  and carry fragments of the target genes inserted between two convergent T7 promoters. Silencing of the VFL3-A subfamily genes was performed using a vector containing a fragment of the VFL3-1 gene (from base 20 to 637) linked to a fragment of the VFL3-2 gene (from base 861 to 1375). The two fragments were amplified by PCR and cloned successively into the vector. The RNAi of the VFL3-B subfamily was performed using a vector containing a VFL3-3 fragment extending from base 328 to 701 and linked to a part of the VFL3-4 gene extending from base 1124 to 1725 (see Additional file 1: Figure S1). Finally, for OFD1 silencing, a fragment extending from base 254 to 1310 was used. Amplifications were performed with high fidelity Phusion DNA polymerase using standard procedures. After cloning, the genes were entirely sequenced to ensure that no error was introduced during the amplification. Possible RNAi-off target effects were analyzed for each construct using the RNAi-off target tool available at ParameciumDB .
RNA extraction and Northern blot analysis
Total RNA was extracted using TRIzol as previously described . For Northern blots, 20 μg of total RNA were denatured and loaded on 1% agarose gels. Electrophoresis in 1× Tris–borate–EDTA transferred to Hybond N+ membranes (GE Healthcare), and hybridization with 32P-labeled DNA probes were performed as previously described . VFL3-1 probe: sequence of the VFL3-1 gene from base 861 to 1375. VFL3-2 probe: sequence of VFL3-2 from base 20 to 637. VFL3-3 probe: sequence of VFL3-3 from base 1354 to 1725. VFL3-4 probe: sequence of VFL3-4 from base 330 to 831. The sequences of these probes are specific for each gene and should not cross-react with the double strand RNA produced by the RNAi vector. Radioactive signals were quantified using ImageQuant.
nd7-1 mutant cells unable to discharge their trichocysts were transformed by microinjection into their macronucleus of filtered and concentrated plasmid DNA containing a mixture of the plasmids of interest (5 µg/µl) digested by Sfi1 and of plasmid DNA directing the expression of the ND7 wild-type gene . Transformants were first screened for their ability to discharge their trichocysts and if so, further analyzed. Microinjection was made under an inverted Nikon phase-contrast microscope, using a Narishige micromanipulation device and an Eppendorf air pressure microinjector.
Gene silencing was performed by the feeding method as previously described . Stationary wild-type cells or growing transformant cells were fed individually with the appropriate double-stranded RNA-expressing HT115 bacteria. They gave rise to independent cell lines containing a population of cells which presented a homogenous phenotype. Growing cells were transferred daily into fresh feeding medium as needed. Control cells were fed with bacteria transformed by the L4440 plasmid or a plasmid containing the ND7 gene involved the exocytosis . Generally, at least three independent experiments were conducted, and routinely, 6–9 clones obtained from each independent experiment were analyzed.
Each immunolabeling was carried out on samples of 20–50 cells as previously described . Briefly, paramecia were permeabilized 5 min in 1% TritonX100 in PHEM buffer (PIPES 60 mM, HEPES 25 mM, EGTA 10 mM, MgCl2 2 mM pH 6.9) fixed in 2% paraformaldehyde in PHEM buffer for 10 min, and rinsed 2 × 15 min in TBST/BSA (Tris 10 mM, NaCl 0.15 M, Tween20, 3% BSA, EGTA 10 mM MgCl2 2 mM). All subsequent steps were performed in this same buffer. Cells were incubated in the primary antibody for 20 min, rinsed twice and incubated in the secondary antibody for 25 min, rinsed in TBST/BSA supplemented with Hoescht 2 mg/ml and then in TBST/BSA without Hoechst.
The monoclonal anti-tubulin antibody 1D5  at a dilution 1:100; a polyclonal affinity purified PolyE anti-tubulin  at 1:7000 dilution; the monoclonal CTS32 directed against epiplasmins diluted 1:20 to decorate the epiplasm , the monoclonal anti-Myc antibody 9E10 (Sigma) at a dilution 1:1000, the polyclonal anti-ciliary rootlet serum (KD2)  at a dilution 1:800, the monoclonal anti-tubulin acetylated TEU318  culture supernatant at a dilution 1:10, and secondary antibodies labeled with Alexa Fluor 488 or 568 from Invitrogen-Molecular Probes (Eugene, OR), or Cy5 (Life Technologies) at a dilution 1:250–1:500.
Cells were observed with a Zeiss Axioskop 2 plus equipped with epifluorescence. Images were acquired with a CoolSnap camera coupled with Metavue. Confocal acquisitions were made with a Leica SP8 equipped a UV diode (line 405), and three laser diodes (lines 488, 552 and 635) for excitation and two PMT detectors. Images stacks were processed with imageJ and Photoshop. STED imaging was performed using a Leica TCS SP8 STED 3X (Leica Microsystems CMS GmbH, Mannheim, Germany). The system was equipped with a WLL ranging from 470 to 670 nm for excitation and with 3D STED lasers at 592, 660, and 775 nm. A 100× 1,4 Oil STED white objective was used to acquire the images. GFP, AF568, and Cy5 were excited at 488, 561, and 643 nm, respectively. Detection ranges were 500–550, 575–625, and 660–700 nm, respectively. A pixel size of 25 nm was used. For deconvolution, SVI Huygens was used.
Countings of basal bodies
Countings of basal body doublets were performed from projections of images acquired at the dorsal side of paramecium cells labeled by the 1D5 antibody. Automatic countings were performed by using the ImageJ software and were refined manually. Seven wild-type cells and eleven VFL3-A inactivated cells obtained from three independent RNAi knockdown experiments were used. The percentage of basal body containing abnormal striated rootlets was calculated from manual counting on ten cells obtained from three independent RNAi knockdown experiments.
For ultrastructural observations, the depleted and control cells were fixed in 1% (v/v) glutaraldehyde and 1% OsO4 (v/v) in 0.05 M cacodylate buffer, pH 7.4 for 30 min. After rinsing and dehydration in ethanol and propylene oxide series, they were embedded in Epon resin. For pre-embedding immunolocalization, the immunostaining process was carried out as described for immunofluorescence using a gold-coupled instead of fluorochrome-coupled secondary antibody. Then the cells were rinsed, fixed and embedded as previously described . For post-embedding immunolocalization, cells were fixed in 3% paraformaldehyde, 0.15% glutaraldehyde, in 0.05 M cacodylate buffer, pH 7.4 at 4 °C for 1 h. After washes and dehydration in ethanol, they were embedded in LR White (London Resin). Thin sections were collected on nickel grids and treated with 0.1 M NH4Cl in 0.1 M PBS and then saturated with 3% BSA and 0.1 M glycine in PBS. Sections were incubated with anti-GFP polyclonal antibody at room temperature for 45 mn. After several washes in PBS, a gold-labeled anti-rabbit IgG (GAR G10, Aurion) diluted 1/50 was applied for 30 mn. The grids were rinsed in PBS and distilled water, and finally stained with uranyl acetate. All ultrathin sections were contrasted with uranyl acetate and lead citrate. The sections were examined with Jeol 1400 (120 kV) or Philips CM10 transmission electron microscopes.
Identification of the VFL3 and OFD1 genes in Paramecium tetraurelia
While a single VFL3 gene is present in Chlamydomonas, four VFL3 genes (VFL3-1, VFL3-2 VFL3-3, and VFL3-4) are encoded in the P. tetraurelia genome which has undergone at least three successive whole-genome duplications (WGD) during the evolution . Phylogenetic analysis indicates that the four P. tetraurelia genes can be grouped into two families (VFL3-A and VFL3-B), the VFL3-A proteins being more closely related to the Chlamydomonas one (Additional file 2: Figure S2). Only one gene of each family is encoded in the Paramecium caudatum genome, a closely related species which shares only the most ancient of the three WGDs with P. tetraurelia . Only one OFD1 ortholog is encoded in the genome. The alignment of the conserved N-terminal part of this protein with its orthologs is presented in Additional file 3: Figure S3.
Role of OFD1 in the basal body anchoring process and distal part assembly
We tested the involvement of OFD1 in basal body anchoring by RNAi knockdown experiments using the feeding method . To silence the unique OFD1 gene, we performed inactivation by bacteria producing dsRNA covering 56% of the gene (1055/1887 nt). The efficiency of the vector to inactivate its target was tested by controlling for the effective depletion of the GFP-tagged protein (see Additional file 4: Figure S4).
In interphase cells, the basal bodies are embedded in cortical units which contain either one or two basal bodies, the distribution of these units being specific of the cell areas. During cell division, neoformed basal bodies are inserted in these units anterior to their parental basal bodies. Depending on the cell territories, several rounds of duplication can occur, generating several new basal bodies from a preexisting one. This is associated with the formation of the new cortical units in the daughter cells . An additional wave of duplication occurs in a subset of cortical units to generate units with basal body pairs. After division, parental and newly assembled basal bodies alternate within the rows and the daughter cells display a global pattern of cortical units identical to that of their mother .
OFD1 knockdown induced a cellular growth defect and cell size reduction after 1 division and cellular death after 4/5 divisions in all the analyzed clones (n > 24) (see “Methods”). To analyze the effects of the depletion on the basal bodies, inactivated cells were labeled by the 1D5 antibody which decorates the basal bodies and an anti-epiplasm antibody specific of the subcortical units. Figure 1 shows the defects induced by OFD1 depletion after one division upon silencing. Many basal bodies are detected in the cytoplasm of OFD1 depleted cells (Fig. 1a right panel). Whereas in the control cells, the basal bodies are located in the middle of the cortical units, in OFD1 depleted cells, basal bodies located between the units or lying under the cell surface are observed (Fig. 1b, c). At the ultrastructural level, whether undocked or partially docked at the membrane, basal bodies (n = 23) show defective distal structures (Fig. 2b, c). The three plates of the transition zone are often absent (65%; n = 15) or partially assembled (35%; n = 8). Rare intracytoplasmic undocked basal bodies prolongated by microtubules are also detected which could result from a premature extension of the basal bodies (Fig. 2d).
OFD1 is recruited early during basal body duplication and localizes at the transition zone proximal part after anchoring
To ascertain the localization of OFD1, we expressed a GFP-tagged protein either under the control of its own regulators or under the control of the constitutive calmodulin gene-regulatory elements. After transformation, we observed the signal on transformants displaying a wild-type growth rate and phenotype. Whatever the vectors used for driving the expression of the tagged protein, the fluorescence localized at the distal end of all basal bodies as demonstrated by labeling with 1D5 (Fig. 3a).
In Paramecium, the basal body duplication process starts at the equatorial zone and propagates towards the cell poles according to a precise spatio-temporal pattern. Consequently, in a dividing cell, several types of cortical units containing basal bodies at different stages of the duplication process coexist: unduplicated, duplicating and fully assembled. Because several rounds of duplication occur during cell division, one basal body can give rise to several new ones, yielding clusters of four to six basal bodies. Observations of dividing transformants (Fig. 3b) allowed us to confirm that, as already described for FOR20 , OFD1 is recruited early during basal body biogenesis before its anchoring at the cell surface. Indeed GFP fluorescent dots, not yet labeled by the anti-tubulin 1D5, are observed close to the preexisting basal bodies, at the position corresponding to the assembly site of the new basal body (Fig. 3b red inset). In the anterior part of the cells, a subset of basal body doublets does not follow a classical duplication pattern. In this region, the anterior basal body of each pair disassembles, its disappearance being followed by the assembly of a new basal body which replaces the disassembled one . In contrast to what is observed in the other parts of the cell (Fig. 3b: green inset vs red inset), a dual GFP labeling is detected near the posterior basal body of the pair when the new basal body assembles: one corresponding to the nascent basal body and the other one corresponding to the anchoring site of the disassembled one. This indicates that OFD1 is stably incorporated at the cell surface as previously shown for FOR20.
We ascertained the localization of OFD1 at the ultrastructural level (Fig. 3c). No morphological basal body defect was observed in transformants expressing GFP-tagged OFD1. Gold particles were detected in about 50% (n = 24) and 30% (n = 14) of basal body longitudinal sections observed after pre- and post-embedding processing, respectively. 87% (pre-embedding) and 95% (post-embedding) of the gold particles were located in the basal body distal part corresponding to the transition zone. Most of these particles localize between the axoneme and the membrane above the terminal plate which marks the limit between the basal body microtubule triplets and the axonemal doublets. This localization is similar to that observed for FOR20.
Relationships between OFD1–FOR20–Centrin 2
We then investigated the relationships between OFD1 and both FOR20 and Centrin 2 which are also required for assembly of the basal body distal part. We followed the fate of the fluorescence in GFP-OFD1 transformants after knockdown of FOR20 and Centrin 2 and carried reciprocal experiments in which OFD1 was silenced in transformants expressing GFP-FOR20 or GFP-Centrin 2 (Additional file 5: Figure S5). After two or three divisions under RNAi conditions, all basal bodies retained fluorescence in the GFP-OFD1 transformants silenced for Centrin 2 as well as in GFP-Centrin 2 transformants inactivated for OFD1. This indicates that there is no relationship between OFD1 and Centrin 2 for their recruitment at the basal bodies. In contrast, the GFP signal decreases in transformants expressing GFP-OFD1 silenced for FOR20 and in GFP-FOR20 transformants depleted of OFD1. This reveals an interdependence between OFD1 and FOR20 for their localization/stability at the basal bodies.
The VFL3-A gene family is involved in basal body anchoring
To ascertain the function of the VFL3-A and VFL3-B genes, two different RNAi vectors targeting specifically and simultaneously both genes of each family were used. As observed for OFD1, the knockdown of the VFL3-A genes induces a growth defect and size reduction after 1 division and cellular death after 4/5 divisions on all the analyzed clones (n > 30) see “Methods”. In contrast, inactivation of the VFL3-B genes does not produce any detectable phenotype. The efficiency and specificity of each vector to silence its targets was confirmed by controlling the effective messenger depletion of the four VFL3 genes in response to the inactivation (see Additional file 6: Figure S6). A very weak signal was detected with the probe corresponding to the VFL3-4 gene (VFL3-B family) in the control as well as in the silenced cells, indicating that the gene is probably expressed at a very low level. The lack of complementation of the VFL3-A genes by the expression of the VFL3-3 isoform (VFL3-B family) indicates that the two VFL3-A and VFL3-B gene subfamilies have different functions. The 63% reduction of the expression of the VFL3-3 gene does not seem sufficient to induce a cellular phenotype.
The effect of VFL3-A depletion on the basal bodies was tested by immunostaining with the 1D5 antibody and the anti-epiplasm antibody. Contrary to what was observed in OFD1 depleted cells, in very rare cases basal bodies were present in the cytoplasm after the first division upon inactivation (Fig. 4a). However, basal bodies mispositioned on the cortex or lying beneath the cell surface (Fig. 4b, c) were observed. An excess of cortical units containing basal body pairs [32 ± 6% in inactivated cells (n = 11) vs 6 ± 1% in control cells (n = 7)] as well as clusters of three basal bodies within one cortical unit, never observed in wild-type interphase cells, were also associated with VFL3-A depletion. This excess can arise either from over-duplication during the assembly of the basal body pairs, or from a defect in the segregation of the new versus old basal bodies after their assembly. Examination of the so-called invariant field, a region where basal bodies are associated in pairs in all units, allowed us to eliminate the first hypothesis (Fig. 5). In this region, cortical units containing a single basal body were observed, indicating that the depletion of VFL3-A could affect in some way new basal body assembly.
VFL3-A controls the correct assembly of the basal body proximal appendages
In Paramecium, three appendages associated with the proximal part of the basal body organize the movements and separation of the basal bodies during the duplication process : one striated rootlet and two microtubular ribbons. We thus analyzed the organization of these appendages in VFL3-A depleted cells. Double staining by KD2 (a marker of the striated rootlets) and 1D5 (a marker of basal bodies) or TEU318 (a marker of the microtubular ribbons) reveals abnormalities in the organization of these structures. The striated rootlets are disoriented relative to the cellular antero-posterior axis and pointed in all directions (Fig. 6a). More surprising is the presence of basal bodies devoid of striated rootlet (11.5 ± 2%) which generally also lack their microtubular appendages (Fig. 6b). In addition, we also observed basal bodies carrying more than one striated rootlet (3 ± 1.5%), a phenotype never observed in the wild-type strain. These alterations affect only the basal bodies assembled during the VFL3 depletion which represent approximatively 50% of the basal bodies present in the daughter cells after the inactivation. Thus, we can estimate that the number of basal bodies which form abnormal appendages during the VFL3-A depletion represents about 30% [(11.5% + 3%) × 2)] of the newly assembled basal bodies.
The ultrastructural analysis confirmed the basal body misorientation in VFL3-A inactivated cells and the defect in assembly of their striated rootlets. In the wild-type strain, a unique striated rootlet is associated with all basal bodies, all of them pointing to the same direction (Fig. 7A1). In addition to the defect in the orientation of these fibers (Fig. 7A2), more than one striated rootlet was found associated with some basal bodies in VFL3-A depleted cells (Fig. 7A3, A4). In most of the VFL3-A depleted cells observed after one division upon inactivation, the non-anchored basal bodies appeared fully assembled and harbored distal structures equivalent to the pro-transition zone present in the wild-type cells (Fig. 7B2). Observations on proliferating basal body units in cells undergoing their first division under VFL3-A depletion (Fig. 7C2, C4) revealed the presence of daughter basal bodies developing at inaccurate sites, with a mature length, but with abnormal orientation and unable to tilt up. Defect of spacing between old and new basal bodies as well as attachment of newly formed basal bodies to their parent by an abnormal filament were also observed as if the severing step was impaired.
One gene of each VFL3 family (VFL3-1 and VFL3-3) was tagged by the GFP or Myc epitope. After transformation, we observed the signal on transformants which displayed a wild-type growth rate and phenotype. A very weak fluorescence associated with the basal bodies was observed at all stages of the cell cycle in transformants expressing GFP or Myc-tagged VFL3-3 (Additional file 7: Figure S7). In contrast, the Myc-tagged VFL3-1 isoform is recruited to the cortex only during the division. No labeling is observed during the interphase (Fig. 8a), whereas double staining with the anti-Myc and 1D5 antibodies reveals a signal during the basal body duplication process (Fig. 8b). This anti-Myc staining was observed regardless of the regulatory sequences used for driving the expression of the fusion protein. Figure 8b illustrates the localization of Myc-VFL3-1 in a dividing cell at the onset of division. In the part of the cell where basal body duplication is in progress (below the dashed line), as judged by a weaker 1D5 labeling of the developing anterior basal body, the Myc labeling is observed. Once the basal bodies are duplicated (above the dashed line), the Myc labeling is no longer detected, indicating that VFL3-1 is present only during an early step of the new basal body assembly.
The localization of VFL3-1 with respect to the basal bodies was analyzed by stimulated emission depletion super-resolution microscopy (STED). The triple labeling by 1D5 (basal bodies), KD2 (striated rootlet) and anti-Myc antibodies indicates that VFL3-1 localizes between the proximal part of the striated rootlets and the microtubule wall of the basal body (Fig. 9).
Relationship between VFL3-A and Centrin 3
Depletion of VFL3-A and Centrin 3, both localized outside the basal bodies, yielded similar defects in the basal body anchoring . To explore the relationships between VFL3-A and Centrin 3 which is required for the basal body movement toward the cell surface, we followed the fate of the fluorescence in GFP-Centrin 3 transformants after knockdown of VFL3-A. After two or three divisions, the GFP signal decreases in indicating that VFL3-A is required for the recruitment of Centrin 3 (Additional file 8: Figure S8). This could indicate that the assembly of the transient ALF structure which requires Centrin 3 for assembly could also be impaired by the VFL3 depletion.
One homolog of OFD1 is encoded in the P. tetraurelia genome, whereas four isoforms clustered into two families (VFL3-A and VFL3-B) encode proteins related to VFL3. The VFL3-A subfamily is more closely related to the Chlamydomonas protein. We show that OFD1 localizes at the proximal part of the transition zone and is required for the assembly of the basal body tip and for correct basal body anchoring at the cell surface. While the recruitments of OFD1 and Centrin 2 proceed independently, the localization/stability of OFD1 and FOR20 at the basal body are interdependent. The knock down of the VFL3-B family which encodes proteins permanently associated with the basal bodies, did not induce any phenotype. In contrast, the VFL3-A proteins are shown to be required for both the correct positioning of the basal body at the cell surface and the regular organization of its appendages which mark its rotational asymmetry. Contrary to the effect of OFD1 depletion, the non-anchored basal bodies display a fully organized distal part. The VFL3-1 protein localizes transiently, at an early stage of the basal body duplication process, to the proximal region of the parental basal body, between the striated rootlet and the microtubules to which it is associated. These results suggest a role of VFL3-A as an extrinsic factor controlling the basal body rotational polarity. The anchoring defect could be due to the loss of this polarity, which specifies the sites of assembly of the different structures which guide the movement of the basal body toward the cell surface.
A role of OFD1 in motile cilia anchoring process
We have shown that OFD1, which localizes at the distal part of the ciliated and unciliated basal bodies, is recruited early during the formation of the basal body. Its loss leads to a defective transition zone and prevents basal body anchoring at the cell surface. This is in agreement with the known function of OFD1 in the assembly of the primary cilium in mammalian cells. Associated to the distal ends of both procentrioles and centrioles, it acts to build the distal appendages  that help to anchor the mother centriole to the ciliary vesicle during the intracytoplasmic ciliogenesis process. Our results do not fit observations in zebrafish which showed that OFD1 invalidation does not completely prevent axoneme extension of the Kupffer’s vesicle cilia. In this case however, it was considered that only a fraction of OFD1 protein was depleted . No clear defect of basal body docking was either observed in brain cells of mutant mice carrying a deletion of exons 4 and 5 of OFD1 . However a cDNA generated by the deletion  encoded a 106 amino acids long protein which retained the TOF/LisH/conserved domain required for centrosomal targeting . The expression of such a protein could thus partially complement the entire OFD1 loss. We cannot exclude the hypothesis that the anchoring process differs from one cellular type to another. Nevertheless, our results are consistent with the docking impairment observed in airways multiciliated cells of patients with OFD1 mutations (Thauvin-Robinet et al. 2013). They also indicate that OFD1 is required in the basal body anchoring process which involves a direct docking at the cell surface.
We localized OFD1 at the proximal part of the transition zone, between the ciliary membrane and the microtubule doublets. This result confirms the observations obtained through super-resolution analysis on primary and motile cilia . After labeling by an anti-OFD1 antibody, a ring of ninefold symmetry was observed at the interface between basal body and the cilium, suggesting that OFD1 could be associated with the microtubule doublets at this interface. Upon OFD1 depletion, we also observed rare intracytoplasmic basal bodies containing microtubules extension beyond the abnormal transition zone. This observation provides further support to the hypothesis that OFD1 could control the elongation of the microtubules doublets .
The interdependence of OFD1 and FOR20 for their recruitment/stability at the basal bodies suggests that they could form a complex within the basal bodies. This is supported by our previous work  which showed that like OFD1, FOR20 is recruited early in the assembly of basal bodies, is localized between the axoneme and the ciliary membrane of the transition zone and remains stably associated to the cortex when the basal body disassembles. There is evidence that such a complex could be formed in association with OFIP at the centrosome and pericentriolar satellite in mammalian cells . It is thus possible that this complex could also be formed during the basal body assembly, since OFIP is conserved in Paramecium. Surprisingly, whereas we previously showed that Centrin 2 is required for FOR20 recruitment, we observed that its loss does not prevent the localization of OFD1 at the basal body. One explanation for this discrepancy could be that, in the absence of Centrin 2, as previously shown for OFIP , FOR20 and OFD1 do not form a complex. In this case, the lack of both Centrin 2 and FOR20 would not influence the stability/recruitment of OFD1.
VFL3-A: a conserved role in the development of basal body appendages and basal body positioning
We show here that the expression of the VFL3-B genes does not complement the inactivation of VFL3-A, indicating that the genes have different functions. A single gene of each family, VFL3-A and VFL3-B, is encoded in the genome of P. caudatum, a Paramecium species which did not undergo the last two whole-genome duplications observed in P. tetraurelia . Two genes were also identified in the genome of Tetrahymena thermophila indicating that the VFL3 gene duplication is at least as old as the origin of the Ciliates. One can thus suppose that after duplication, one copy of VFL3 developed a function completely distinct of the ancestral gene in Ciliates. Another possibility is that several functions were originally performed by the ancestral gene and then have been independently retained in each isoform after sub-functionalization. In Paramecium, the inability of the basal bodies in the VFL3-A depleted cells to segregate, to dock at the cell surface as well as to assemble their appendages at the right position, is quite similar to the phenotypical defects observed in the VFL3 mutant of Chlamydomonas . The genes of the Paramecium VFL3-A family have thus conserved at least one common ancestral function with the VFL3 gene of Chlamydomonas. Further analyses will be required to understand the function of the VFL3-B family.
While OFD1 depletion leads to defects in basal body tip assembly, VFL3-A depleted cells display abnormalities in the organization of their associated rootlets, i.e., striated rootlets and microtubular ribbons. In addition, VFL3 is also required for the recruitment of Centrin 3 and thus, although not demonstrated here, its depletion would most probably affect the assembly of the ALF, a transient Centrin 3 dependent structure . All these appendages have been demonstrated to play a crucial role in basal body positioning, i.e., basal body segregation and tilting  and required for the basal anchoring at the cell surface These data strongly suggest that abnormalities of the system of rootlets result in basal body segregation and anchoring defects observed during the VFL3 depletion. This is in agreement with the conclusion of Wright et al.  concerning the effect of the VFL3 mutation in Chlamydomonas and with studies in Tetrahymena which indicate that the DisAp mutation which affect the length of the striated fibers cause defect in basal body orientation [44, 45].
VFL3-A an actor in the establishment of the basal body rotational asymmetry
Among the different phenotypical abnormalities which accompany the anchoring defect in VFL3-depleted cells, the ability of basal bodies to assemble several striated rootlets is the most surprising. The association of single striated rootlet per basal body is one character shared by all Ciliates. In all species, this rootlet has the same orientation with respect to the cell polarity and the direction of the ciliary beating [46, 47], and the same function of anchoring the basal bodies to the cell surface [37, 48, 49]. Together with two other rootlets, the postciliary and transverse microtubules ribbons, the striated rootlet marks the basal body rotational asymmetry revealing the specificity of each of the nine microtubule triplets . In wild-type cells, the asymmetry of the basal body and associated rootlets is strictly transmitted from one generation to the next during basal body duplication . The fact that several striated rootlets may be associated with basal bodies after VFL3-A depletion strongly suggests that this protein could be required for the acquisition of triplet specification allowing the establishment of rotational asymmetry.
The mechanisms of establishment of the rotational asymmetry are still unknown and remain an open debate: does this asymmetry result from an intrinsic property of the basal body/centriolar structure or is it externally imposed by the environment? One of these factors could be the mispositioning of the nascent basal body observed during VFL3-A depletion. The fact that such abnormalities are not observed after Centrin 3 depletion, which also generates mispositioned fully assembled basal bodies , suggests that the mispositioning per se does not induce a loss of polarity.
Two markers which localize asymmetrically within the basal body lumen in Chlamydomonas reveal an intrinsic asymmetry of the basal bodies: the acorn, a complex structure which attaches to specific triplets  and VFL1, a molecule which localizes near a subset of the microtubules triplets that faces one of its associated fibers . We show here for the first time that a protein required for the expression of this asymmetry, VFL3, localized transiently outside the basal body during assembly of the neoformed one. This result indicates that the asymmetry is at least partially imposed by an external factor. VFL3-A could be a crucial actor in the asymmetry by allowing the recruitment of external or internal proteins, such as VFL1, which induce a specification of the microtubular triplets of the basal body in Chlamydomonas . Further experiments on the respective roles of VFL3 and VFL1 in the establishment of the circumferential asymmetry of the basal bodies are needed to analyze these mechanisms.
The present study extends our knowledge on the molecular mechanisms involved in the basal body anchoring/positioning process in the multiciliated cell P. tetraurelia. We show that this process depends on the function of the OFD1 and VFL3-A protein families. Their loss prevents insertion and orientation of the basal bodies along the cellular antero-posterior axis and induces an accumulation of subcortical basal bodies. Although the phenotypical defects are quite similar, the underlying mechanisms are different and highlight the importance of the assembly of both distal and proximal structures of the basal body in the orientation and the docking process at the cell surface. In addition, our results suggest that the mispositioning of the basal bodies at the cell surface might result of the loss of the rotational asymmetry of their microtubular triplets.
anterior left filament
FOP-related protein of 20 kDa
oro-facial-digital syndrome 1 protein
FOR20 interacting protein
variable flagellar number mutant 3
variable flagellar number mutant 1
- Immuno EM:
Sorokin SP. Reconstructions of centriole formation and ciliogenesis in mammalian lungs. J Cell Sci. 1968;3:207–30.
Sorokin S. Centrioles and the formation of rudimentary cilia by fibroblasts and smooth muscle cells. J Cell Biol. 1962;15:363–77.
Hagiwara H, Shibasaki S, Ohwada N. Ciliogenesis in the human oviduct epithelium during the normal menstrual cycle. J Electron Microsc (Tokyo). 1992;41:321–9.
Park TJ, Mitchell BJ, Abitua PB, Kintner C, Wallingford JB. Dishevelled controls apical docking and planar polarization of basal bodies in ciliated epithelial cells. Nat Genet. 2008;40:871–9.
Burke MC, Li F-Q, Cyge B, Arashiro T, Brechbuhl HM, Chen X, et al. Chibby promotes ciliary vesicle formation and basal body docking during airway cell differentiation. J Cell Biol. 2014;207:123–37.
Wei Q, Ling K, Hu J. The essential roles of transition fibers in the context of cilia. Curr Opin Cell Biol. 2015;35:98–105.
Dippell RV. The development of basal bodies in paramecium. Proc Natl Acad Sci USA. 1968;61:461–8.
Aubusson-Fleury A, Lemullois M, de Loubresse NG, Laligné C, Cohen J, Rosnet O, et al. The conserved centrosomal protein FOR20 is required for assembly of the transition zone and basal body docking at the cell surface. J Cell Sci. 2012;125:4395–404.
Dute R, Kung C. Ultrastructure of the proximal region of somatic cilia in Paramecium tetraurelia. J Cell Biol. 1978;78:451–64.
Aubusson-Fleury A, Lemullois M, Bengueddach H, Abdallah S, Shi L, Cohen J, et al. Transition zone: the sequential assembly of its components parallels its dual role in basal body anchoring and ciliary function. Cilia. 2015;4:P26.
Iftode F, Cohen J, Ruiz F, Rueda AT, Chen-Shan L, Adoutte A, et al. Development of surface pattern during division in Paramecium. I. mapping of duplication and reorganization of cortical cytoskeletal structures in the wild type. Development. 1989;105:191–211.
Iftode F, Addoute A, Fleury A. The surface pattern of paramecium tetraurelia in interphase: an electron microscopic study of basal body variability, connections with associated ribbons and their epiplasmic environment. Eur J Protistol. 1996.
Jerka-Dziadosz M, Koll F, Włoga D, Gogendeau D, Garreau de Loubresse N, Ruiz F, et al. A Centrin 3-dependent, transient, appendage of the mother basal body guides the positioning of the daughter basal body in Paramecium. Protist. 2013;164:352–68.
Ruiz F, Garreau de Loubresse N, Klotz C, Beisson J, Koll F. Centrin deficiency in Paramecium affects the geometry of basal-body duplication. Curr Biol. 2005;15:2097–106.
Ferrante MI, Giorgio G, Feather SA, Bulfone A, Wright V, Ghiani M, et al. Identification of the gene for oral-facial-digital type I syndrome. Am J Hum Genet. 2001;68:569–76.
Singla V, Romaguera-Ros M, Garcia-Verdugo JM, Reiter JF. Ofd1, a human disease gene, regulates the length and distal structure of centrioles. Dev Cell. 2010;18:410–24.
Azimzadeh J, Nacry P, Christodoulidou A, Drevensek S, Camilleri C, Amiour N, et al. Arabidopsis TONNEAU1 proteins are essential for preprophase band formation and interact with centrin. Plant Cell. 2008;20:2146–59.
Sedjaï F, Acquaviva C, Chevrier V, Chauvin J-P, Coppin E, Aouane A, et al. Control of ciliogenesis by FOR20, a novel centrosome and pericentriolar satellite protein. J Cell Sci. 2010;123:2391–401.
Chevrier V, Bruel A-L, Dam TJPV, Franco B, Scalzo ML, Lembo F, et al. OFIP/KIAA0753 forms a complex with OFD1 and FOR20 at pericentriolar satellites and centrosomes and is mutated in one individual with oral-facial-digital syndrome. Hum Mol Genet. 2016;25:497–513.
Wright RL, Chojnacki B, Jarvik JW. Abnormal basal-body number, location, and orientation in a striated fiber-defective mutant of Chlamydomonas reinhardtii. J Cell Biol. 1983;96:1697–707.
Gupta GD, Coyaud É, Gonçalves J, Mojarad BA, Liu Y, Wu Q, et al. A dynamic protein interaction landscape of the human centrosome-cilium interface. Cell. 2015;163:1484–99.
Jakobsen L, Vanselow K, Skogs M, Toyoda Y, Lundberg E, Poser I, et al. Novel asymmetrically localizing components of human centrosomes identified by complementary proteomics methods. EMBO J. 2011;30:1520–35.
Skouri F, Cohen J. Genetic approach to regulated exocytosis using functional complementation in Paramecium: identification of the ND7 gene required for membrane fusion. Mol Biol Cell. 1997;8:1063–71.
Sonneborn TM. Chapter 12 methods in Paramecium research. In: Prescott DM, editor. Methods in cell biology. Cambridge: Academic Press; 1970. p. 241–339.
Timmons L, Fire A. Specific interference by ingested dsRNA. Nature. 1998;395:854.
Arnaiz O, Cain S, Cohen J, Sperling L. ParameciumDB: a community resource that integrates the Paramecium tetraurelia genome sequence with genetic data. Nucleic Acids Res. 2007;35:D439.
Aury J-M, Jaillon O, Duret L, Noel B, Jubin C, Porcel BM, et al. Global trends of whole-genome duplications revealed by the ciliate Paramecium tetraurelia. Nature. 2006;444:171–8.
Lepère G, Bétermier M, Meyer E, Duharcourt S. Maternal noncoding transcripts antagonize the targeting of DNA elimination by scanRNAs in Paramecium tetraurelia. Genes Dev. 2008;22:1501–12.
Galvani A, Sperling L. RNA interference by feeding in Paramecium. Trends Genet. 2002;18:11–2.
Beisson J, Bétermier M, Bré MH, Cohen J, Duharcourt S, Duret L, et al. Immunocytochemistry of Paramecium cytoskeletal structures. Cold Spring Harb Protoc. 2010
Wehland J, Weber K. Turnover of the carboxy-terminal tyrosine of alpha-tubulin and means of reaching elevated levels of detyrosination in living cells. J Cell Sci. 1987;88(Pt 2):185–203.
Janke C, Bulinski JC. Post-translational regulation of the microtubule cytoskeleton: mechanisms and functions. Nat Rev Mol Cell Biol. 2011;12:773–86.
Nahon P, Coffe G, Guyader H, Darmanaden-Delorme J, Jeanmaire-Wolf R, Clerot JC, et al. Identification of the epiplasmins, a new set of cortical proteins of the membrane cytoskeleton in Paramecium. J Cell Sci. 1993;104:975–90.
Sperling L, Keryer G, Ruiz F, Beisson J. Cortical morphogenesis in Paramecium: a transcellular wave of protein phosphorylation involved in ciliary rootlet disassembly. Dev Biol. 1991;148:205–18.
Callen A-M, Adoutte A, Andrew JM, Baroin-Tourancheau A, Bré M-H, Ruiz PC, et al. Isolation and characterization of libraries of monoclonal antibodies directed against various forms of tubulin in Paramecium. Biol Cell. 1994;81:95–119.
McGrath CL, Gout J-F, Doak TG, Yanagi A, Lynch M. Insights into three whole-genome duplications gleaned from the Paramecium caudatum genome sequence. Genetics. 2014;197:1417–28.
Aubusson-Fleury A, Bricheux G, Damaj R, Lemullois M, Coffe G, Donnadieu F, et al. Epiplasmins and epiplasm in paramecium: the building of a submembraneous cytoskeleton. Protist. 2013;164:451–69.
Romero MR, Torres A. Cortical development associated with conjugation of Paramecium. Development. 1993;117:1099–112.
Iftode F, Fleury-Aubusson A. Structural inheritance in Paramecium: ultrastructural evidence for basal body and associated rootlets polarity transmission through binary fission. Biol Cell. 2003;95:39–51.
Ferrante MI, Romio L, Castro S, Collins JE, Goulding DA, Stemple DL, et al. Convergent extension movements and ciliary function are mediated by ofd1, a zebrafish orthologue of the human oral-facial-digital type 1 syndrome gene. Hum Mol Genet. 2009;18:289–303.
D’Angelo A, De Angelis A, Avallone B, Piscopo I, Tammaro R, Studer M, et al. Ofd1 controls dorso-ventral patterning and axoneme elongation during embryonic brain development. PLoS ONE. 2012;7:e52937.
Ferrante MI, Zullo A, Barra A, Bimonte S, Messaddeq N, Studer M, et al. Oral-facial-digital type I protein is required for primary cilia formation and left-right axis specification. Nat Genet. 2006;38:112–7.
Lee YL, Santé J, Comerci CJ, Cyge B, Menezes LF, Li F-Q, et al. Cby1 promotes Ahi1 recruitment to a ring-shaped domain at the centriole–cilium interface and facilitates proper cilium formation and function. Mol Biol Cell. 2014;25:2919–33.
Galati DF, Bonney S, Kronenberg Z, Clarissa C, Yandell M, Elde NC, et al. DisAp-dependent striated fiber elongation is required to organize ciliary arrays. J Cell Biol. 2014;207:705–15.
Jerka-Dziadosz M, Jenkins LM, Nelsen EM, Williams NE, Jaeckel-Williams R, Frankel J. Cellular Polarity in ciliates: persistence of global polarity in a disorganized mutant of Tetrahymena thermophila that disrupts cytoskeletal organization. Dev Biol. 1995;169:644–61.
Lynn DH. The Organization and evolution of microtubular organelles in ciliated protozoa. Biol Rev. 1981;56:243–92.
Sleigh MA, Silvester NR. Anchorage functions of the basal apparatus of cilia. J Submicrosc Cytol. 1983;15:101–4.
Peck RK. The ultrastructure of the somatic cortex of Pseudomicrothorax dubius: structure and function of the epiplasm in ciliated protozoa. J Cell Sci. 1977;25:367–85.
Williams NE, Vaudaux PE, Skriver L. Cytoskeletal proteins of the cell surface in Tetrahymena I. Identification and localization of major proteins. Exp Cell Res. 1979;123:311–20.
Grain J. Le cinétosome et ses dérivés chez les ciliés. Ann Biol. 1969;8:59–97.
Beisson J, Sonneborn TM. Cytoplasmic inheritance of the organization of the cell cortex in Paramecium aurelia. Proc Natl Acad Sci USA. 1965;53:275–82.
Geimer S, Melkonian M. The ultrastructure of the Chlamydomonas reinhardtii basal apparatus: identification of an early marker of radial asymmetry inherent in the basal body. J Cell Sci. 2004;117:2663–74.
Silflow CD, LaVoie M, Tam LW, Tousey S, Sanders M, Wu W, et al. The Vfl1 Protein in Chlamydomonas localizes in a rotationally asymmetric pattern at the distal ends of the basal bodies. J Cell Biol. 2001;153:63–74.
Dereeper A, Guignon V, Blanc G, Audic S, Buffet S, Chevenet F, et al. Phylogeny.fr: robust phylogenetic analysis for the non-specialist. Nucleic Acids Res. 2008;36:W465–9.
HB, ML, AAF, and FK made the experiments and prepared the manuscript. All authors read and approved the final manuscript.
We thank Dr. Mireille Betermier for her help in RNA extraction, Northern blot, and hybridizations. We thank all the members of the laboratory for fruitful discussions and their critical reading of the manuscript. We are indebted to Leica Microsystems for making a STED microscope available. We are very grateful to Dr. Ulf Schwarz for making the STED imaging. The present work has benefited from the core facilities of Imagerie-Gif, (http://www.i2bc.paris-saclay.fr), member of IBiSA (http://www.ibisa.net), supported by “France-BioImaging” (ANR-10-INBS-04-01), and the Labex “Saclay Plant Science” (ANR-11-IDEX-0003-02).
The authors declare that they have no competing interests.
Availability of data and materials
Data are available on request.
This study was supported by the Agence Nationale de la Recherche: Grant No. ANR15-CE11-0002-01.