PCARE requires coiled coil, RP62 kinase-binding and EVH1 domain-binding motifs for ciliary expansion

Abstract Retinitis pigmentosa (RP) is a genetically heterogeneous form of inherited retinal disease that leads to progressive visual impairment. One genetic subtype of RP, RP54, has been linked to mutations in PCARE (photoreceptor cilium actin regulator). We have recently shown that PCARE recruits WASF3 to the tip of a primary cilium, and thereby activates an Arp2/3 complex which results in the remodeling of actin filaments that drives the expansion of the ciliary tip membrane. On the basis of these findings, and the lack of proper photoreceptor development in mice lacking Pcare, we postulated that PCARE plays an important role in photoreceptor outer segment disk formation. In this study, we aimed to decipher the relationship between predicted structural and function amino acid motifs within PCARE and its function. Our results show that PCARE contains a predicted helical coiled coil domain together with evolutionary conserved binding sites for photoreceptor kinase MAK (type RP62), as well as EVH1 domain-binding linear motifs. Upon deletion of the helical domain, PCARE failed to localize to the cilia. Furthermore, upon deletion of the EVH1 domain-binding motifs separately or together, co-expression of mutant protein with WASF3 resulted in smaller ciliary tip membrane expansions. Finally, inactivation of the lipid modification on the cysteine residue at amino acid position 3 also caused a moderate decrease in the sizes of ciliary tip expansions. Taken together, our data illustrate the importance of amino acid motifs and domains within PCARE in fulfilling its physiological function.


Introduction
The vertebrate retina captures, amplifies and propagates light signals using intricately connected layers of photoreceptors, retinal neurons and optic nerve cells. This mechanism relies on the ability of photoreceptors to capture light photons and convert them into electrochemical signals via the action of opsin proteins that populate the stacked, disk-shaped membranes of photoreceptor outer segments (OSs). Photoreception and transduction place these OS under high oxidative stress, and therefore damaged components of the membrane disks require constitutive renewal of the OS. This is achieved by shedding of the OS tips at the rate of ∼10% per day, which is matched by OS biogenesis at the base (1). To support the highly metabolically active OS, a photoreceptor contains a biosynthetic inner segment that is connected to the OS by a microtubule-based bridge called the connecting cilium that is largely homologous to the transition zone of a primary non-motile cilium (2,3). The connecting cilium enters the OS in the region that cradles new opsin-loaded disks, which are formed by evagination and expansion of the ciliary plasma membrane (4,5). A branched F-actin network localizes at the site of evagination; however, an exact mechanism of its involvement in disk morphogenesis remains elusive (6). Mutations in PCARE (photoreceptor cilium actin regulator) (C2orf71) have been linked to autosomal recessive retinitis pigmentosa type 54 (RP54) (7,8). PCARE has been found to localize to the evagination sites of photoreceptors where it interacts with WASF3 (WAVE3, SCAR3) protein (9). WASF3 is an actin modifier that is able to activate the Arp2/3 complex and facilitate assembly of a branched F-actin network (10,11). PCARE expression is specific to the retina; however, when PCARE is recombinantly expressed in cultured ciliated cells, it also localizes specifically to the cilium. Moreover, upon co-expression of WASF3, ciliary tip expansions appear that harbour both PCARE and WASF3 proteins. The formation of ciliary tip expansions is inhibited by F-actin poisons and siRNA-mediated downregulation of Arp2 expression (9).
This study aims to further investigate the role of PCARE in the actin dynamics of the ciliary tip by scrutinizing the structure of this novel protein. PCARE is an intrinsically disordered protein, which limits attempts to determine its molecular structure using crystallographic methods (12). However, advances in the in silico predictions of protein structures allow us a glimpse into its structural and functional motifs. Previously, we described a tentative actin-binding WH2 short linear motif (SLiM) between amino acids 597 and 615 (9). WH2 is one of several SLiMs found in proteins that are involved in the regulation of actin filament assembly (13). However, when cotransfected with WASF3 in hTERT RPE-1 and IMCD-3 cell lines, PCARE lacking this WH2 motif retained the ability to form ciliary tip expansions (data not shown). This finding, together with an absence of robust motif representation in the eukaryotic linear motif (ELM) database (13) and low conservation score, prompted us to further study the structure-function relationship within the PCARE protein. We designed a bioinformatic study to predict evolutionary conserved amino acid motifs of PCARE. Subsequently, we revealed important functional amino acid motifs by expressing the PCARE mutant constructs lacking the motifs in hTERT RPE-1 and IMCD-3 cell lines and assessing their ability to localize to the cilia and form ciliary tip expansions. This work thereby points to a potential mechanism of actin dynamics during photoreceptor disc formation in health and disease.

Bioinformatics analysis of PCARE protein showed structural and amino acid motifs potentially involved in interaction with actin and proteins involved in actin remodelling
PCARE is a retina-specific protein consisting of 1289 amino acids with a predicted molecular weight of 140 kDa. On the basis of the prediction using the DISOPRED and AlphaFold servers, the PCARE protein is mostly intrinsically disordered with a likely structured region between amino acids 171 and 338 (14,15). This region corresponds to a helical domain framed with two long helices forming an antiparallel coiled coil predicted with both the PSIPRED workbench and AlphaFold (16). The amino acid residue analyses showed no other strongly predicted structural domains (Supplementary Material, Fig. S1). Bioinformatic analysis revealed that the glycine residue at amino acid position 2 (p.Gly2) and the cysteine residue at amino acid position 3 (p.Cys3) may undergo myristoylation and palmitoylation lipid modifications, respectively (7). In the middle part of the protein, two target sites for the MAK (RP62) kinase were predicted (amino acid positions 513-516 and 577-580) and, in a proline-rich region, three enabled/VASP homology 1 domain (EVH1)-binding motifs (amino acids positions 805-809, 830-834 and 1056-1060).

PCARE mutants show stable protein expression comparable with wild-type
We modified recombinant PCARE cDNA to express a series of HA-tagged PCARE protein variants containing disrupted amino acid motifs (Fig. 1). To prevent lipid modification on the second glycine and third cysteine residues, the cDNA sequence was changed to encode the catalytically inert amino acid residue alanine at these positions. All other predicted amino acid motifs were disrupted by deleting the entire motif. We observed that the mutant PCARE constructs were stably expressed when transfected into HEK293T cells and their levels of expression were comparable to that of the recombinant wildtype protein (Fig. 2).

Disturbance of predicted PCARE motifs differentially affects ciliary characteristics
When transfected into hTERT RPE-1 cells, wild-type PCARE has been shown to localize along the cilia (9). We repeated this experiment with wild-type and mutant PCARE isoforms and observed ciliary localization for each of them except for the mutant harbouring the helix domain deletion (PCARE helix) where the protein localized into the cytosol only (Fig. 3).
When co-transfected with WASF3, PCARE induces the formation of ciliary tip expansions. We classified three ciliary phenotypes: type 1, where PCARE and WASF3 co-localize into the cilia and form expansions of the ciliary tip; type 2, where the proteins co-localize to the cilia but form no expansions and type 3, where cilia also appear normal but PCARE and WASF3 retain within the cytosol and/or nucleus (Fig. 4, Supplementary Material, Fig. S3). With wild-type PCARE, 50% of the colocalization events resulted in the formation of ciliary expansions of type 1 per total double-positive cilia (n = 271 number of double-positive cilia counted, SD = 7.8%). Mutant constructs formed less of type 1 expansions, for example, PCARE EVH1a, b and c formed 39%, 37% and 48% of type 1 expansions, respectively (n = 71, SD = 8.8%, n = 70, SD = 8.5% and n = 62, SD = 10.5%), whereas PCARE RP62a and PCARE RP62b (n = 71, SD = 10.0% and n = 66, SD = 9.6%) formed 48% and 42% of expansion, respectively ( Fig. 4A-C). We observed the second most prominent decrease in type 1 expansions in PCARE EVH1a-c (14%, n = 61, SD = 5.2%). PCARE helix (n = 34) failed to localize into the cilia and thus did not form expansions at all. Furthermore, we quantified the size of the ciliary tip expansions formed upon co-transfection of PCARE mutant constructs with WASF3 in hTERT RPE-1 cells. We observed the mean expansion size of wild-type PCARE to correlate with previous observations (0.97 μm 2 , n = 117, SD = 0.4 μm 2 ) (9). Our results showed that most mutated constructs formed smaller expansions than the wildtype. From those, the PCARE EVH1a and b constructs formed significantly smaller expansions with means at 0.82 μm 2 and 0.75 μm 2 , respectively (n = 73, SD = 19 μm 2 , P = 0.001 and n = 37, SD = 0.2 μm 2 , P < 0.0001). The tripledeletion PCARE EVH1a-c showed the smallest mean expansion size of 0.71 μm 2 (n = 10, SD = 0.22 μm 2 , P = 0.006). Also, a removal of lipid modification from p.Cys3 might result in a decreased expansion size (size = 0.84 μm 2 , n = 36, SD = 0.28 μm 2 , P = 0.035) (Fig. 4D). At the N-terminus of PCARE, the glycine residue at amino acid position 2 (p.Gly2), and the cysteine residue at amino acid position 3 (p.Cys3) exhibit tentative palmitoylation and myristoylation lipid modifications, respectively. As predicted with DISOPRED server, a helical structural domain is located between amino acids 171 and 338. Using Jalview, the following conserved SLiMs were predicted: two RP62 kinase-binding motifs, and three EVH1 domainbinding motifs. To generate PCARE mutant proteins, the p.Gly2 and p.Cys3 residues were substituted by the inactive amino acid alanine residue, whereas predicted protein domains that are annotated in this schematic were deleted via site-directed mutagenesis, including one mutant in which all three EVH1 domain-binding motifs were deleted.
The absence of lipid modification on p.Gly2 or deletions of RP62 candidate phosphomotifs did not exhibit a statistically significant difference in ciliary expansion size ( Table 1).

Transfection of PCARE mutant constructs into stable IMCD-3 expressing WASF3 confirms decreased ability to form tip expansions
To minimize the variability of the double-transfection experimental design, we also generated IMCD-3 cell lines stably expressing WASF3 that allowed for a standardized level of WASF3 expression between samples (Supplementary Material, Fig. S2). In the IMCD-3 cells, the average percentage of cilia with type 1 phenotype increased by 12% while type 3 decreased by 9%, whereas the average SD of expansion sizes decreased by 16%. Most of the results observed in IMCD-3 cells stably expressing WASF3 followed the trends of hTERT RPE1 cell results (Fig. 5). For example, removal of lipid modification from the cysteine residue at amino acid position 3 also showed a decrease inthe ciliary tip expansion sizes Figure 2. PCARE mutants are expressed and stable. Wild-type and mutant PCARE constructs were transfected in HEK293T cells. Western blot analysis revealed stable expression of all PCARE constructs. Wild-type PCARE protein has a molecular weight of 140 kDa, the molecular weight of mutant PCARE proteins is approximately the same, except for the mutant carrying the helix deletion which is smaller. Tubulin was loaded for control of the total amount of protein. In contrast to hTERT-RPE1, no effect on mean expansion size was observed for PCARE EVH1b (0.77 μm 2 , SD = 0.24 μm 2 , P = 0.052), but deletion of each MAK/RP62 When counted in the IHC slide and quantified, with wild-type PCARE normalized to hundred percent, all PCARE mutants generated less ciliary tip expansions than the wild-type (type 1 in yellow, type 2 in green, type 3 in red). (B) hTERT RPE-1 cells were transfected with 3×HA-PCARE (wild-type or mutant) and 3×FLAG-WASF3 constructs in equal amounts. Three distinct types of the ciliary tip expansions can be distinguished: type 1, where PCARE and WASF3 co-localize into the cilia and form expansions of the ciliary tip, type 2, where the proteins co-localize to the cilia but formed no expansions and, finally, type 3, where the proteins are retained in the cytosol (3×HA-PCARE in green, WASF3 in red, ARL13B in magenta, DAPI in blue). Untransfected cells (Utr) are shown as a negative control. Scale bars 5 μm. (C) The sizes of the ciliary tip expansions of type 1 were counted and quantified using an automated Fiji script. The triple-deletion PCARE EVH1a-c and single-deletions of PCARE EVH1a and PCARE EVH1b showed a significant decrease in the size of the expansion tip ( * * * P ≤ 0.001, * * P ≤ 0.01, * P ≤ 0.05, red line shows mean and SD). Furthermore, for PCARE EVH1c, a decrease of expansion size was also observed (mean = 0.75 μm 2 , n = 115, SD = 0.26, P = 0.014) ( Table 1).

Discussion
Mutations in PCARE are causative for retinal ciliopathy RP54 (7,8). Recently, we observed that PCARE localizes at the initiation site of new OS disks within photoreceptors and proposed its role in protein transport and actin remodelling that is required to form OS disks (9). We identified candidate SLiM motifs of PCARE that may be relevant for interactions with actin or actin modifiers and generated stable PCARE protein mutants with deletion or inactivation of these motifs. Our data strongly suggest the importance of the alpha-helix domain for PCARE ciliary trafficking and EVH1 domain-binding motifs and myristoylation on the third cysteine for modification of the actin filament network in the ciliary tip. The aim of this study was to investigate the structurefunction relationship within PCARE with regard to its ability to regulate actin dynamics. To address this, we performed a bioinformatics analysis of PCARE protein structure to find evolutionary conserved structural and functional amino acid motifs. From the list of predicted SLiMs, we selected the motifs that are potentially involved in interaction with actin or with other proteins involved in actin remodelling. We identified two target sites for a photoreceptor ciliumassociated MAK/RP62 that phosphorylates serine and threonine residues (17). Dysfunction of MAK kinase is linked to structural defects in photoreceptor cilia and recessive retinitis pigmentosa (18)(19)(20). MAK knockout mice retina shows elongated cilia, rhodopsin accumulation and photoreceptor degeneration (21). Furthermore, we identified three EVH1 domain-binding SLiMs that are specific to proteins that are involved in actin-based cytoskeleton reorganization, actin dynamics modulation and actin motility (22)(23)(24). We also modified amino acids (p.Gly2 and p.Cys3) at the N-terminus of PCARE, because of their known role as a lipid modification sites in transport and membrane localization of proteins (17). Finally, PCARE contains one predicted alpha-helical domain, one of the most common secondary structure motifs found within proteins (18). To access the function of each of these motifs, we generated PCARE mutant constructs with deletion or inactivation of the motifs of interest on the basis of wild-type PCARE cDNA.
To determine the inf luence of the predicted amino acid motifs on the PCARE stability, we expressed the PCARE mutant constructs in HEK293T cells and assessed PCARE expression levels. All mutant proteins were expressed at levels comparable to wild-type PCARE. Furthermore, to elucidate the function of the selected motifs in PCARE ciliary trafficking, we investigated the localization pattern of PCARE mutants in ciliated hTERT RPE-1 and IMCD-3 cells. Previously, we showed that, when separated into arbitrarily chosen fragments, the most Nterminal fragment of PCARE, which harbours the alphahelical coiled coil domain, localizes to the basal body and axoneme of the cilia and is sufficient to transport WASF3 into the cilia to form ciliary tip expansions (9). Interestingly, the most C-terminal fragment of PCARE alone was also capable of localization into the ciliary axoneme.
Here, we showed that full-size PCARE lacking the helix domain failed to localize to the cilia. This finding hints towards a crucial role that the minimal structural organization of PCARE has in its ciliary localization and further work is required to establish the exact role of the coiled coil domain in PCARE trafficking. Moreover, myristoylation and palmitoylation were predicted at the N-terminus of PCARE, modifications that are known to play a role in ciliary trafficking in Trypanosoma brucei and in eukaryotic cells (19,20). N-myristoylation on p.Gly2 is a protein modification where 14-carbon saturated fatty acid myristate is covalently attached to the N-terminal Gly (exposed after removal of p.Met1) of the target protein catalyzed by N-terminal myristoyltransferases. On p.Cys3, palmitoylation where 16-carbon fatty acid palmitate is covalently attached to a cysteine residue is predicted, a process that is catalyzed by the Nterminal palmitoyl acyltransferases. Myristoylation at the N-terminus also facilitates palmitoylation (25,26). The lipid-modified protein is directed to a membrane by the action of N-myristoylation, but palmitoylation is required for stable anchoring to the phospholipid bilayer (27). In X-linked retinitis pigmentosa protein 2 (RP2), myristoylation alone was not sufficient to traffic amino acid residues 1 to 16 RP2 to the cilia (20). We deactivated palmitoylation on Gly2 or myristoylation on Cys3 by Gly-to-Ala or Cys-to-Ala substitutions, but observed no effect on the ciliary trafficking. Further analysis of human retinal-organoid derived, or murine purified PCARE protein, for example by mass spectrometry, would be required to determine if p.Gly2 or p.Cys3 is respectively myristoylated or palmitoylated in vivo. Taken together, our data show that all PCARE mutants except for PCARE helix localized to the cilia, suggesting redundancy or potentially undiscovered amino acid motifs to play a role in ciliary trafficking of PCARE.
Next, we accessed the frequency and size of ciliary tip expansions generated by PCARE mutant proteins when co-expressed with WASF3, either via double transfection of both proteins in hTERT RPE-1 cells, or by a single transfection of PCARE in IMCD-3 cells stably expressing WASF3. All mutants showed a decreased frequency of ciliary tip expansions compared with the wild-type PCARE. In IMCD-3 cells stably expressing WASF3, there is a general increase in the ability to form ciliary tip expansions for all protein variants compared with double-transfected hTERT RPE-1. Compared with the wild-type, both PCARE (p.Gly2Ala), where the myristoylation was inactivated and PCARE (p.Cys3Ala), where the palmitoylation was inactivated, showed a lower number of ciliary tip expansions of type 1. When quantifying the sizes of the ciliary tip expansions, the deactivation of myristoylation had no statistically significant effect but the deactivation of palmitoylation seemed to have negatively affected the ability of PCARE to form ciliary tip expansions. Potentially, the palmitoylation-affected mutant impaired membrane anchoring and, thus, fewer units of the protein complex may reach the ciliary tip leading to fewer and smaller bulges.
Deletion of the EVH1 domain-binding motifs, both sequentially and concurrently, also resulted in fewer and smaller ciliary tip expansions in both hTERT RPE-1 and IMCD-3 cell lines. Analysis of the sequential deletions showed that the deletion of the most evolutionary conserved EVH1b motif affects the tip expansions most strongly in both cell lines. Analysis of concurrent EVH1 domain-binding SLiM deletions showed the strongest effect when all three EVH1 domain-binding motifs were deleted simultaneously. Although this effect is stronger than that of the individual deletions, more research should be done to explore a potential synergistic effect of the EVH1 deletions. Putative EVH1 domain-binding motifs would allow PCARE to interact with proteins that have the EVH1 (WH1, RanBP1-WASP; Pfam:PF00568) domain that is present in species ranging from yeast to mammals. Such proteins are involved in a wide range of signalling, nuclear transport and cytoskeletal events. Many EVH1-containing proteins are actin interactors and participate in actin cytoskeleton organization. The EVH1 domain-binding motifs of PCARE consist of five amino acid stretches containing three or more proline residues and resemble the EVH1 domain-binding motifs in other proteins such as the Ena/VASP family proteins, which are thought to assist in actin filament formation by profilinactin complexes to their proline-rich binding partners (28). For example, EVH1-binding motifs might facilitate interaction with Protein-enabled homolog (ENAH) and Vasodilator stimulated phosphoprotein (VASP) PCARE interactor proteins that we have identified in yeast two-hybrid screening and tandem affinity purification studies (9). A very recent structural study has shown that PCARE binds with a high affinity to the EVH1 domain of the protein ENAH, which as noted we had previously identified as a PCARE interactor, that plays a role in actin polymerization and cytoskeletal remodelling (29). Taken together, these studies establish that PCARE acts through ENAH to regulate actin filament assembly in the eye photoreceptor cilium.
MAK/RP62 kinase candidate phosphorylation site mutants PCARE RP62a and PCARE RP62b showed a similar frequency of expansion formation to the wildtype. Possibly, that although MAK/RP62 kinase is crucial for ciliary length and retinal health, it does not directly affect actin network polymerization via PCARE. However, for PCARE RP62a and PCARE RP62b, we observed a subtle but significant difference between hTERT RPE-1 and IMCD-3 cell lines in the ciliary tip expansion size.
When comparing the results obtained in IMCD-3 and hTERT RPE-1 cells, it is clear that although most of the findings were evident in both cell lines, the immortalized cell line models might provide contradictive readouts. Being derived from mouse kidney or human retinal epithelial cells, these models cannot fully represent the function of photoreceptor cells. Therefore, further confirmation of PCARE protein's functions is required with model systems that more faithfully represent the human retina, for example in human iPSC-derived retinal organoids (30). The recent advances in these techniques allow for precise editing of the gene of interest to induce mutations via the CRISPR-Cas9 editing machinery, which can be further studied at a different level of complexity such as protein expression, morphology or even function (20).
To conclude, all PCARE mutants that localized into the cilia formed ciliary tip expansions, suggesting that PCARE contains a redundancy of amino acid motifs via which it can interact with actin modulators and that it may interact with multiple actin-modulators. Our data also pinpoint the prominent role of proteins with EVH1 domains and hint towards the importance of palmitoylation lipid modification in the PCARE-actin interactions. This work presents a first glance at the importance of amino acid modules in PCARE and expands upon their structure and function in relationship to actin network modification, paving the way for improved understanding of mechanisms that underlie PCARE-associated retinal disease, as well as the formation of the OS disks of photoreceptors.

Bioinformatic prediction of structural and functional amino acid motifs within PCARE protein
The human PCARE SEQUENCE (UniProt:A6NGG8) was submitted to the Eukaryotic Linear Motif resource (http://elm.eu.org) server (13). Among the candidate matches were three proline-rich EVH1-binding SLiMs (805-809, 830-834, 1056-1060: http://elm.eu.org/elms/ LIG_EVH1_1.html). A multiple sequence alignment of vertebrate PCARE sequences was prepared and visualized with Jalview (31). All vertebrate PCARE sequences had multiple EVH1 motifs, with some variation in number and positioning, as often occurs with SLiMs. Since these candidates are in regions of predicted native disorder, amino acid conservation is driven by functional, not structural selection: Thus, these motifs are strong candidates for actin-regulatory function via EVH1 domain binding (32). By visual inspection, we also noted two conserved RPx[ST] motifs (513-516, 577-580) that are candidate phosphorylation sites (33) for the MAK (RP62) ciliary protein kinase (UniProt:P20794) required for long-term photoreceptor survival (33). All vertebrate PCARE sequences begin with the residues MetGlyCys indicative of N-terminal myristoylation and prenylation. Structural context was provided by native disorder predictors IUPRED (34) and DISOPRED (14) and the PSIPRED 4.0 (http://bioinf.cs.ucl.ac.uk/psipred/) workbench (16). More recently, the structural context was re-evaluated with AlphaFold (35), supporting the disorder predictions and providing high confidence assignment for the long helices (171 and 338) as an anti-parallel coiled coil pair (see Supplementary Material, Fig. S1).

Generation of mutant PCARE constructs
Human C2orf71/PCARE [NCBI GeneID:388939] full-length cDNA was obtained by PCR using Human Retina Marathon ® -Ready cDNA (Clontech) and cloned into Gateway™ pDONR201 vector (Thermo Fisher) as previously described (9). PCARE mutant constructs lacking functional and structural motifs were constructed by performing site-directed mutagenesis on pDONR with cloned PCARE cDNA construct with primer pairs that introduced deletions of nucleotide coding for the amino acid motifs via PCR amplification with Q5 High-fidelity DNA polymerase (New England Biolabs, M0491L). PCR reactions were performed using 0.2 ml PCR tubes in a final volume of 50 μl including the following components: 50 ng of template DNA, 0.5 μM of each primer (forward and reverse), 200 μM of each dNTPs, 1× Q5 Buffer and 1.6 U of Q5 High-fidelity DNA polymerase. 30 cycles of amplification were carried out under the following conditions: initial denaturation at 98 • C for 4 min, denaturation at 98 • C for 30 s, annealing at 50 • C for 30 s and extension at 72 • C for 3 min, following by final extension at 72 • C for 7 min. Primer sequences are listed in Supplementary Material, Table S1. Parental plasmid DNA was eliminated by DpnI treatment (1 μl of 20 U/μl) (New England Biolabs, R0176S) of a 10-μl aliquot of the linear amplification reaction for 1 at 37 • C. DNA was then directly transformed to competent Escherichia coli DH5α cells and purified with NucleoSpin Plasmid EasyPure (BIOKÉ, MN740727). Constructs were checked for deletions by amplification of the 200 nucleotide target region around the deletion with forward and reverse primers with Taq polymerase (Roche, 11 647 687 001). Thirty-five cycles of amplification were carried out under the following conditions: initial denaturation at 94 • C for 5 min, denaturation at 94 • C for 30 s, annealing at 58 • C for 30 s and extension at 72 • C for 20 s, following by final extension at 72 • C for 5 min. At the completion of the PCR, fragments were separated on 2.5% agarose gels containing 0.5 μg/ml ethidium bromide and further verified by Sanger sequencing. The entire coding region of constructs containing the modification of interest was verified by Sanger sequencing. The obtained entry clones were used to generate expression clones containing 3×HA tag using Gateway Technology (Life Technologies) and according to the manufacturer's instructions.

Statistical analysis
Statistical analyses were performed in Excel and Graphpad/Prism 6. An ALPACA tool (version 1.53c) was used for quantification of sizes of the ciliary tip expansions (38). Error bars represent standard deviations. An unpaired t test with Welch's correction was used to compare each mutant cilium tip expansions' sizes to the wild-type. N values are stated in the figures and represent cilia. Pvalues are stated in the figures and symbols indicate the following P-values: ns, P > 0.05; * P ≤ 0.05; * * P ≤ 0.01; * * * P ≤ 0.001.

Supplementary Material
Supplementary Material is available at HMG online.