A Mec17-Myosin II Effector Axis Coordinates Microtubule Acetylation and Actin Dynamics to Control Primary Cilium Biogenesis

Primary cilia are specialized, acetylated microtubule-based signaling processes. Cilium assembly is activated by cellular quiescence and requires reconfiguration of microtubules, the actin cytoskeleton, and vesicular trafficking machinery. How these components are coordinated to activate ciliogenesis remains unknown. Here we identify the microtubule acetyltransferase Mec-17 and myosin II motors as the key effectors in primary cilium biogenesis. We found that myosin IIB (Myh10) is required for cilium formation; however, myosin IIA (Myh9) suppresses it. Myh10 binds and antagonizes Myh9 to increase actin dynamics, which facilitates the assembly of the pericentrosomal preciliary complex (PPC) that supplies materials for cilium growth. Importantly, Myh10 expression is upregulated by serum-starvation and this induction requires Mec-17, which is itself accumulated upon cellular quiescence. Pharmacological stimulation of microtubule acetylation also induces Myh10 expression and cilium formation. Thus cellular quiescence induces Mec17 to couple the production of acetylated microtubules and Myh10, whose accumulation overcomes the inhibitory role of Myh9 and initiates ciliogenesis.


Introduction
Primary cilia are evolutionarily conserved, microtubule-based organelles critical for detecting and transmitting mechanical and chemical cues. The biological functions of primary cilia have long been overlooked until the discovery of a cohort of cilia-related human developmental disorders, including Bardet-Biedl syndrome (BBS) [1], Joubert syndrome [2] and Merkel-Gruber syndrome [3].
Human genetic studies in combination with biochemical and cell biological approaches have identified the basic components and mechanisms underlying primary cilium formation and function [4,5]. When ciliogenesis is initiated upon cellular quiescence, the mother centriole translocates to the cortical plasma membrane and forms the basal body, from where the ciliary microtubules are polymerized and form the axoneme. In coordination with axoneme growth, specialized vesicles become concentrated around the basal body and provide new membranes and proteins to support cilium growth [6]. Disruption of this pericentrosomal preciliary compartment (PPC), which is enriched for Rab11 positive recycling endosomes and proteins important for membrane fusion and transport including Rab8, PCM-1, and Cep290, leads to defects in cilia formation [6][7][8]. The coordination of PPC assembly with microtubule-axoneme growth is thus critical for ciliogenesis but its molecular basis remains poorly understood.
Primary cilium formation requires the reorganization of cellular cytoskeleton, particularly microtubules which provide both structural components and intraflagellar transport [9]. One salient feature of ciliary microtubules is the prevalent acetylation on lysine (K)-40 of a-tubulin [10]. The presence of acetylated microtubules, in fact, is the most commonly used marker for primary cilia although its exact function in the cilium remains uncertain. a-tubulin acetylation is primarily controlled by the acetyltransferase Mec-17 (also termed alphaTAT1) [11,12] and the deacetylase HDAC6 [13]. Mec-17 knockdown does not eliminate primary cilium formation [14,15]; however it disrupts the normal kinetics of cilium biogenesis [12]. On the other hand, HDAC6 has been proposed to facilitate primary cilium resorption [16]. While these findings suggest a regulatory role of microtubule acetylation in primary cilium formation, how the production of acetylated microtubules is coupled to ciliogenesis is not known.
In addition to microtubules, several components of the actin cytoskeleton were recently identified as cilium regulators [7]. The analyses of these factors have revealed a general inhibitory role of the actin cytoskeleton in ciliogenesis where stable actin cytoskeleton prevents the formation of PPC enriched for Rab11positive recycling endosomes [7]. Thus, reorganization of the actin network is required for efficient delivery of membranes and materials for cilium growth. Supporting this view, microRNA-129-3p, a positive regulator for cilium formation, targets genes involved in the formation of a branched actin network [17]. The mechanism by which quiescent cells release the inhibitory brake enforced by the actin network to activate ciliogenesis remains to be characterized.
In this report, we provide evidence that non-muscle myosin IIA and IIB and the tubulin acetyltransferase Mec-17 form the central molecular circuit that controls cilium formation. We show that myosin IIB (Myh10) promotes, whereas IIA (Myh9) inhibits, ciliogenesis. The opposing activity of Myh10 and Myh9 is mediated through the actin dynamics, which in turn controls PPC assembly. We found that Myh10 expression is positively regulated by the tubulin acetyltransferase Mec-17. Importantly, both Mec-17 and Myh10 gene expression are induced by serum starvation conditions that activate cilium formation. In Mec-17deficient cells subject to serum starvation, Myh10 is not induced and ciliogenesis is deregulated. Conversely, pharmacological inhibition of HDAC6 increases microtubule acetylation, Myh10 expression and cilium formation. Our results suggest that Mec-17-dependent microtubule acetylation is coupled to the induction of Myh10, whose accumulation overcomes Myh9-dependent actin cytoskeleton stabilization and promote primary cilium assembly.

Myh10 Is Required for Cilium Formation
To investigate the role of the actin network in ciliogenesis, we focused on nonmuscle myosin II motors, which are known to regulate actin network dynamics and actinomyosin-microtubule crosstalk [18,19]. We adopted the serum starvation protocol to induce cilium formation in RPE-Mchr1 GFP cells, a retina pigment epithelial cell line expressing GFP-fused melanin-concentrating hormone receptor 1 (Mchr1) protein [20]. To study the role of myosin IIA (Myh9) and IIB (Myh10) in cilium formation, we first designed siRNA duplexes specifically targeting Myh10 or Myh9 mRNA. Western blot analysis and Q-PCR confirmed the effective knockdown of Myh10 and Myh9 expression by siRNAs (S1A Figure). When RPE-Mchr1 GFP cells were subjected to serum deprivation for 24 hours to induce ciliogenesis, Myh9 siRNA-treated cells were able to form cilia as efficiently as control siRNA group (Fig. 1A). In stark contrast, two different Myh10 siRNA duplexes both dramatically inhibited cilium formation (Fig. 1A). Knocking down Myh10, but not Myh9, in mouse IMCD3 cells also potently inhibited cilium formation in IMCD3 cells ( Fig. 1B and Fig. S1B). Moreover, re-introducing a wild type human Myh10 restored cilium formation in Myh10-knockdown (KD)-IMCD3 cells (Fig. 1C). These findings show that Myh10 is required for efficient ciliogenesis. It is also worthwhile to point out that Myh10 knockdown does not stop cells from proliferating before reaching confluency in our experiments, suggesting the effect of Myh10 on ciliogenesis is not an indirect consequence of altered cell cycle.
Cilia are important for sonic hedgehog signaling [21]. To confirm that Myh10 is required for cilia-mediated signal transduction, we stimulated cells with SAG [21], a small molecule agonist of Smoothened, to induce sonic hedgehog signaling in serum-starved cells. We measured Gli-1 mRNA induction, a downstream target of hedgehog signaling, by RT-PCR. As shown in Fig. 1D, Myh10 (red bar) but not Myh9 (green bar) knockdown strongly inhibited Gli-1 induction upon SAG stimulation, indicating that Myh10 is required for cilium-mediated sonic hedgehog signal transduction. These experiments collectively demonstrated that Myh10 is required for proper ciliogenesis.

Myh10 Antagonizes Myh9 Activity in Cilium Formation
To further ascertain the role of Myh10 in cilium formation, we determined whether an increase of Myh10 expression is sufficient to promote the growth of , Myh10 siRNA #1 and #2 (middle panels) and Myh9 siRNA duplex (right panel) and then serum starved for 36 hours. Cells were fixed by methanol for indirect fluorescence staining using rabbit anti-glu-tubulin (pseudo colored with green) and mouse anti-c-tubulin (red). Nuclei cilia. To this end, we stably expressed Myh10-GFP in IMCD3 cells. Myh9-GFP [22] was included as a control. We then labeled acetylated tubulin to examine cilium formation in these two stable cell lines. Myh10-GFP IMCD3 cells form cilia at a modestly higher percentage than control cells ( Fig. 2A, left panel, middle row, and right panel, red bar); however, cilium length in Myh10-GFP expressing cells is significantly longer (Fig. 2B). This data indicate that increasing Myh10 expression stimulates cilia elongation. Unexpectedly, cilium formation was dramatically inhibited in Myh9-GFP expressing IMCD3 cells ( Fig. 2A, left panel, bottom row, and right panel, green bar). Thus Myh9, in contrast to Myh10, appears to be a negative regulator of cilium formation.
As overexpression of Myh9 and Myh10 caused opposite phenotypes on cilia morphology, we tested if Myh10 promotes cilium formation by antagonizing Myh9 activity, as these two myosin motors can form hetero-oligomers (Fig. 2F). To test this possibility, we knocked down both Myh9 and Myh10 expression by siRNA in RPE-Mchr1 GFP cells and assessed cilium formation. In contrast to Myh10 single knockdown, where few cilia were observed, concurrent knockdown of Myh9 restored cilium formation in more than 40% of the cells (Fig. 2C). Further supporting the conclusion that Myh9 is epistatic to Myh10, application of a myosin II-specific inhibitor, Blebbistatin, also markedly restored cilium formation in Myh10-knockdown cells (Fig. 2D). We noted that the cilia formed under these conditions were generally shorter than those in wild type control cells, which suggests incomplete rescue of Myh10 function by Myh9 inhibition. These data indicate that Myh9 activity is, at least partially, responsible for cilium inhibition caused by Myh10 knockdown. We next asked whether Myh9 activity might function to inhibit ciliogenesis in the presence of high serum when ciliogenesis is typically low. As shown in Fig. 2E, knocking down Myh9 expression indeed significantly increased ciliogenesis. These experiments indicate that Myh9 suppresses ciliogenesis under non-permissive conditions (e.g. high concentration of serum) and Myh10 can counteract this activity.
were counter stained with DAPI. Scale bar: 10 mm. Percentage of ciliated cells (acetylated tubulin staining longer than 1 mm) were quantified under confocal microscopy (right panel). Error bars represent standard deviations (s.d.) from triple replicates. ***, t-test p,0.001. (B) Myh10 is required for ciliogenesis in IMCD3 cells. IMCD3 cells were transfected with siRNA duplexes targeting Myh9 or Myh10 (left panel) and were allow to grow confluent for 72 hours after transfection. Cells were methanol fixed and immunostained with mouse anti-acetylated tubulin (red) or rabbit anti-c-tubulin (green) to label cilia and centrosomes respectively. Yellow insets were magnified to allow better visualization of representative primary cilium images in each knockdown group (left panel, bottom). Percentage of ciliated IMCD3 cells from each knockdown group was quantified under the confocal microscopy (right panel). Error bars represent standard deviations (s.d.) from triple replicates. ***, t-test p,0.001.

Myh10 Antagonizes Myh9 and Increases Actin Dynamics
Myh9 is known to stabilize actin bundles and the filament network [18]. We therefore asked whether Myh10 might oppose Myh9 and increase actin network dynamics. To test this possibility, we first performed Phalloidin staining in ARPE-19 cells (an immortalized human retina pigment epithelia cell line) to examine the influence of Myh10 knockdown on the F-actin network. As reported previously, Myh10 knockdown did not disrupt stress fiber formation (S2A Figure, arrow heads), which is controlled by Myh9 [18]. Myh10 knockdown, however, dramatically increased F-actin staining intensity at the cell edge and the cell cortex ( Fig. 3A and S2B Figure), suggesting the stabilization of the F-actin network. To quantitatively measure actin dynamics, we performed a fluorescence recovery after photobleaching (FRAP) assay using b-actin-GFP in ARPE-19 cells. We found that in contrast to Myh9 knockdown, which accelerated the initial phase of fluorescence recovery after photobleaching, Myh10 knockdown dramatically blocked the recovery (Fig. 3B, S2C Figure). This result indicates that loss of Myh10 leads to a more stable actin network.
The areas with enhanced F-actin staining in Myh10-knockdown cells were also enriched for p34/Arc, a marker for the branched actin network (S2D Figure). Branched actin network has been proposed to prevent the concentration of pericentrosomal recycling endosomes that provide membranes for primary cilium growth [7,17]. We therefore examined whether Myh10 regulates pericentrosomal recycling endosome cluster formation. As shown in Fig. 3C, Myh10 knockdown dramatically reduced the percentage of cells with Rab11-positive recycling endosomes clusters, consistent with the function of Myh10 in maintaining actin dynamics. Taken together, Myh10, contrary to Myh9, functions to maintain a dynamic actin network and supports pericentrosomal recycling endosome cluster formation.

Pcm-1/Cep290 Localization
The requirement of Myh10 for the assembly of pericentrosomal recycling endosome pool prompted us to examine whether Myh10-dependent actin dynamics are also necessary for clustering other pericentrosomal components required for cilium formation. We examined Cep290 and PCM-1 proteins, which are co-localized to centriolar satellites and recruit Rab8 to facilitate the fusion of preciliary vesicles to the ciliary membrane [8]. In Myh10-knockdown cells pericentrosomal Cep290 staining showed a marked decrease in signal intensity (Fig. 4A) whereas Myh9 knockdown had little effect (S3A Figure). More detailed examination revealed a specific loss of Cep290 at the centriolar satellite whereas centrosome-associated Cep290 was retained (Fig. 4A, left panels, insets). PCM-1 was also dispersed from the pericentrosomal region in Myh10-knockdown cells (Fig. 4B, and S3B Figure). The overall protein levels of Cep290 and PCM-1 remained unaltered in Myh10-knockdown cells (S3C Figure). These results show that the Cep290-PCM-1 preciliary complex is mis-localized in Myh10-deficient cells.
Next, we tested if actin dynamics underlie Myh10-dependent pericentrosomal structure organization. We treated Myh10-knockdown RPE-Mchr1 GFP line with a low dosage (20 nM) of Latrunculin A, an actin depolymerization reagent, to increase actin dynamics in Myh10-knockdown cells. 24 hours after Latrunculin A (LatA) treatment, pericentrosomal PCM-1 was partially restored (Fig. 4C). Inhibiting Myh9 by the myosin II-specific inhibitor, Blebbistatin (25 mM), also significantly restored PCM-1 clusters in Myh10 knockdown cells (Fig. 4C, Blebbi panels). These data indicate that Myh10 promotes pericentrosomal concentration of the Cep290-PCM-1 complexes by opposing Myh9 and increasing dynamics of the actin cytoskeleton.  The fact that Myh10 promotes cilium elongation prompted us to examine whether Myh10 expression is regulated during cilium formation. We monitored the expression of Myh10 during cilium formation induced by serum starvation (Fig. 5A). Interestingly, Myh10 protein levels gradually increased during the early phase of cilium formation (0-12 h after starvation, Fig. 5B, C). In contrast, Myh9 protein levels were relatively constant. Real-time PCR analysis showed that Myh10 mRNA levels increased within 2 hours of serum starvation whereas Myh9 mRNA level was not induced (Fig. 5D). These data reveal that relative abundance of Myh10 vs. Mhy9 increases in response to serum starvation, a condition that would favor primary cilia formation.
We next investigated the mechanism by which Myh10 is induced prior to cilium formation. As cilia are heavily enriched for acetylated microtubules, we determined if microtubule acetyltransferase Mec-17 affects Myh10 expression. We transduced RPE-Mchr1 GFP cells with several independent Mec-17 shRNA lentiviral particles. Two (#1 and #3) independent shRNA constructs effectively reduced acetylated tubulin levels (Fig. 6A) and Mec-17 mRNA (S4A Figure). Remarkably, in both Mec-17 shRNA knockdown cell lines, Myh10 protein level was markedly reduced (Fig. 6A) in the presence of serum whereas Myh9 remained relatively unchanged. The decrease in Myh10 protein abundance was not due to protein degradation as inhibiting proteasome function by MG132 did not restore Myh10 level (data not shown). Indeed, quantitative PCR detected a significant reduction of Myh10 mRNA (,60-70%) in Mec-17-knockdown cells (S4B Figure), suggesting that Mec-17 regulates Myh10 protein expression via a transcriptional mechanism. A modest reduction of Myh9 mRNA was also observed although its protein level was not affected (S4B Figure). Importantly, Myh10 was not induced by serum starvation in Mec-17 knockdown cells (Fig. 6G, Top Panel). Collectively, these results showed that Mec-17 is required for Myh10 upregulation in response to serum starvation.
The tubulin acetyltransferase Mec-17 has been indicated to possess both catalytic-dependent and -independent functions [23]. To determine if Myh10 expression depends on Mec-17 catalytic activity, we expressed wild type and catalytic-dead (CD) D157N mutant of mouse Mec-17 [24]. As shown in Fig. 6B, wild type Mec17 restored Myh10 expression to different degrees in two Mec-17 KD lines that correlated with the extent of tublin acetylation whereas catalytic-dead (CD) D157N mutant was unable to restore tubulin acetylation or Myh10 expression. Of interest, over-expression of Mec-17 in control KD cells increased tublin acetylation and Myh10 levels (Fig. 6B, lane Control shRNA, WT). These results indicate that Mec-17 induces Myh10 via its acetyltransferase activity, likely by acetylating microtubules. Indeed, treatment with tubastatin A (TBSA), a specific inhibitor for tubulin deacetylase HDAC6 [25] significantly increased tubulin acetylation and Myh10 expression in RPE-Mchr1 GFP cells, whereas sodium butyrate (NaB), a HDAC inhibitor that does not target HDAC6, had no effect (Fig. 6C). Importantly, TBSA treatment significantly increased spontaneously formed cilia in the presence of serum, which further supports the critical role of Myh10 protein expression in cilium formation (Fig. 6D). Altogether, these  experiments support the conclusion that Mec-17 promotes Myh10 expression by increasing microtubule acetylation.
If Mec-17 is required for Myh10 expression, one would expect defective ciliogenesis in Mec-17 knockdown cells. Indeed, we found that Mec-17 knockdown caused PCM-1 dispersion (Fig. 6E) and ciliogenesis was significantly reduced in the early phase (Fig. 6F, 12 hours after serum withdrawal). Similar to Myh10 knockdown cells, the ciliogenesis defect in Mec-17 KD cells can be corrected by Blebbistatin treatment that inactivates Myh9 (Fig. 6F, red line at 12 hours). Intriguingly, cilium formation in Mec-17 knockdown cells eventually recovered at later time points (Fig. 6F, green), as it was previously reported [12]. This later phase recovery of ciliogenesis is also consistent with normal cilia formation observed in Mec-17-deficient mice [14,15]. Analysis of Myh9 protein in Mec-17-knockdown cells revealed that it was decreased in later time points after serum starvation (12-24 hours, Fig. 6G). Given the inhibitory activity of Myh9 on cilium formation, the loss of Myh9 could explain the eventual recovery of cilium formation in Mec-17 knockdown cells. Our data reveal that Mec17 regulates ciliogenesis, at least in part, through Myh10 induction.

Discussion
The formation of primary cilium requires polarized growth of the membrane material delivered by the PPC complex and a buildup of acetylated microtubules. Our study suggests that these processes are coordinated by myosin II and tubulin acetyltransferase Mec-17. Actin cytoskeleton, especially branched actin network, plays an inhibitory role in PPC complex recruitment. Evidence suggests that the disassembly of the actin network is required for cilium formation [7,17]. In this study, we have identified non-muscle myosin IIA (Myh9) and IIB (Myh10) as the key regulators of the actin cytoskeleton that controls PPC complex and cilium formation. Pericentrosomal recycling endosomes are one of the main organelles that provide a membrane source to support cilia growth [26]. Cep290-PCM-1 complex is also localized to the pericentrosomal region where it recruits Rab8 to facilitate the vesicular trafficking to the cilia [8]. Our results show that Myh10 is required for pericentrosomal recruitment of Cep290, PCM-1 as well as Rab11positive recycling endosomes. Importantly, defects in PPC recruitment in Myh10 KD cells can be partially rescued by inhibiting Myh9 or destablizing actin cytoskeleton (Fig. 4). Thus, Myh9 and Myh10 are the critical regulators of the actin cytoskeleton that controls the concentration of cilia-associated proteins and organelles for ciliogenesis. It is also important to note that Myh10 is required to maintain pericentrosomal organization of recycling endosome clusters and Cep290-PCM-1 satellites even in the presence of serum. This finding suggests that basal level of Myh10 and actin dynamics may be a pre-requisite to prime ciliogenesis by facilitating preciliary structure assembly.
Despite their homology, Myh9 and Myh10 knockdown display an opposite phenotype in ciliogenesis. Myosin IIA (Myh9) has been shown to form oligomers to stablize and promote the assembly of F-actin-based structures [18,[27][28][29]. Our analysis indicates that Myh10 opposes Myh9 activity, potentially by forming hetero-oligomers (Fig. 2F), to increase actin dynamics and stimulate cilium formation (Fig. 3). The antagonistic activity of Myh9 and Myh10 suggests that the relative abundance or activity of Myh9 and Myh10 could determine the state of cilia. Therefore, the observed induction of Myh10 upon serum starvation could potentially act as a part of the molecular switch that releases the break on cilogenesis enforced by the inhibitory actin cytoskeleton when accumulated Myh10 overcomes Myh9.
Interestingly, we found that Myh10 expression is positively regulated by the tubulin acetyltransferase Mec-17. The regulation of Myh10 requires Mec-17 acetyltransferase activity (Fig. 6B), suggesting that microtubule acetylation is coupled to Myh10 transcriptional induction. Indeed, an increase in tubulin acetylation by inhibiting HDAC6 also led to Myh10 induction and enhanced cilium formation (Fig. 6C-D). As cilia primarily consist of acetylated microtubules, this finding suggests an intriguing link between microtubule acetylation and the activation of cilium assembly. Indeed, knocking down Mec-17 prevented Myh10 induction and delayed normal ciliogenesis in quiescent cells (Fig. 6F). This microtubule acetylation-dependent regulation of gene transcription may represent a novel function of microtubule network to regulate cellular signaling events. We recently reported similar regulation in anti-inflammatory response of LPS-treated macrophage cells and suggested that this regulation may present in broader cellular processes. In macrophage anti-inflammatory response, LPS-treatment stimulated tubulin acetylation which amplified the p38 kinase cascade and activated downstream SP1-dependent transcription of IL-10 [30]. We hypothesized that similar acetylated microtubule-dependent signal amplification mechanisms may be involved in the induction of Myh10 during cell quiescenceinduced ciliogenesis but may depend on other kinases and transcription factors. Although further experiments would be required to determine how acetylated microtubules control Myh10 expression, the Mec-17-dependent Myh10 induction provides a potential mechanism for a coordinated regulation of acetylation of microtubules and delivery of ciliary membranes and proteins, two key elements for cilium assembly. Interestingly, Mec-17 mRNA level is also upregulated by serum starvation (S4C Figure). These findings suggest a simple model for ciliogenesis where Mec17, upon cellular quiescence, is activated to catalyze microtubule acetylation for building ciliary axoneme. Microtubule acetylation also induces Myh10 expression, which increases the dynamics of the actin network by counteracting Myh9, thereby assisting the assembly of PPC and enabling delivery of membrane and protein components for cilium growth. Recent genetic mouse models suggest that Mec-17 is dispensable for in vivo ciliogenesis [14,15], our data, however, indicates that Mec-17 regulates in the early steps of cilia formation and kinetics control. Whether a faster kinetics of ciliogenesis is physiological relevant in any in vivo context remains to be investigated. Nevertheless, our finding that increasing tubulin acetylation facilitates ciliogenesis may point to a potential intervention of various human diseases resulting from cilia defects and therefore may be worth further testing in animal models of ciliopathies.

Materials and Methods
Cell Lines, Constructs, Sirnas, Antibodies and Chemical Reagents

Cell Culture, Transfection and Cilia Formation Assay
RPE-Mchr1 GFP cells were maintained in DMEM (Gibco) containing 10% FBS (Fetal Bovine Serum) in a 37˚C humidified incubator with 5% CO 2 . ARPE-19 and IMCD3 cells were maintained in DMEM:F12 medium (Gibco) containing 10% FBS. Cells were passaged at 80%-90% confluence at a regular basis. To induce cilia formation, cells were allow to grow confluent on a coverslip before switching to DMEM or DMEM: F12 containing 0.2% FBS. For the time course study of ciliogenesis, confluent cell layers were first stimulated with 20% FBS for 6 hours to disassemble spontaneously formed cilia and then switched to medium containing 0.2% FBS. Transient transfections of plasmids or siRNA duplexes were performed with lipofectamine 2000 or RNAiMAX (Life Technologies) following manufacturer's guide. For stable transfection using lentivirus or retrovirus, virus was packaged with appropriate packaging plasmids (D8.9 and pCMV-VSVG for lentivirus; pCL-Ampho for retrovirus) in 293T cells. Viruses were collected 48 hours after transfection, filtered through 0.45 mm sterile cellulose acetate membrane filters. Viruses were added to cell culture medium and selected with puromycin (1 mg/ml) for a week. Survived cells were then collected or passaged for analysis. Multiple rounds of infection were performed in a need basis. For non-virus based stable transfection, plasmids were delivered into the cell by Lipofectamine 2000 transfection. Transfected cells were then split at 1:10 and selected with corresponding antibiotics (1 mg/ml puromycin or 0.5 mg/ml G418) for 3 weeks. Individual cell colonies were picked and expanded for further analysis. To maintain all stable cell lines, culture medium contained antibiotics (1 mg/ml puromycin or 0.5 mg/ml G418).

Immunofluorescence and F-Actin Staining
For immunofluorescence, cells were fixated with ice-cold methanol (for antibodies targeting acetylated-tubulin, c-tubulin, Cep290 and PCM-1) or 4% paraformaldehyde in PBS (for Rab11), washed with PBS, blocked with 10% normal goat serum, then incubated with primary antibodies diluted in 10%

Protein Interaction
To detect protein-protein interaction between Myh9 and Myh10, different amount of MSCV-HA-Flag-Myh10 plasmids (0, 1.5, 3.0, 6.0, 9.0 ug) were transfected into 293T cells using lipofectamine 2000 (Invitrogen) following the manufacturer's guide. Cells were collected 48 hours after transfection and lysed in immunoprecipitation buffer (150 mM NaCl, 20 mM pH 8.0 Tris.HCl, 0.5 mM EDTA, 0.5% NP-40). Lysate containing 1 mg total protein were incubated with anti-FLAG M2 affinity gel (Sigma #A2220) for 2 hours at 4˚C and washed 3 times with same lysis buffer. Precipitates were then mixed with protein loading buffer and boiled for western blot analysis.

Fluorescence Recovery after Photobleaching (Frap)
ARPE-19 cells were first transfected with Actin-GFP plasmid and then with control, Myh9 and Myh10 siRNA. 48 hours after siRNA transfection, cells were switched to warm phenol-free DMEM for FRAP analysis. Cells were imaged for 10 seconds, then photobleached at selected regions and imaged for 3 minutes to allow fluorescence recovery. For each knockdown group, a total of five cells were used for analysis. Intensity data from 5 cells were analyzed using FRAP profiler plugin for imageJ, averaged and plotted. The fluorescence recovery curve fitting was performed with Prism 6 software (GraphPad Software) under the exponential association function.

Real-Time Pcr Analysis of Gene Expression
Total RNA was extracted from cells using RNAeasy mini kit (Qiagen), quantified using a spectrometer. 1 ug of total RNA from each sample was reverse transcribed using Promega M-MLV Reverse Transcriptase to be used as template for realtime PCR analysis. cDNA and primers for specific genes (listed in table 1) were mixed in realtime PCR master mix (Qiagen). Real-time PCR reaction was performed on a Realplex Mastercycler (Eppendorf). Data was analyzed with Excel.

Data Processing and Statistical Analysis
All the quantifications were from three independent experiments. To quantify the percentage of ciliated cells, 100-150 randomly chosen cells were counted under the microscopy. Only cells with acetylated tubulin staining longer than 1 mm were considered as bearing mature cilia and counted as positive. For cilia length analysis, cilia lengths were determined using the line scan function in Fiji imaging processing software. Lengths of cilia were plotted as cumulative probability distribution curve. Data are demonstrated as mean ¡ s.d. Two-group hypothesis testing were performed with student t-test. Results were considered statistically significant when p,0.01.
Supporting Information S1 Figure