{'fraction_non_alphanumeric': 0.03230429086144776, 'fraction_numerical': 0.018015067147068457, 'mean_word_length': 5.026400197384653, 'pattern_counts': {'":': 0, '<': 3, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Actin and microtubules drive differential aspects of planar cell polarity in multiciliated cells

Actin dynamics are required for proper cilia spacing, global coordination of cilia polarity, and coordination of metachronic cilia beating, whereas cytoplasmic microtubule dynamics are required for local coordination of polarity between neighboring cilia.


Introduction
Ciliated epithelia are found throughout nature, performing loco motor functions in a variety of aquatic organisms and providing directed fluid flow in a variety of developmental and physio logical contexts. The loss of directed fluid flow in humans can result in hydrocephaly, situs inversus, infertility, and respiratory dysfunction (Zariwala et al., 2007). The epithelium of Xenopus laevis embryos is covered with multiciliated cells, which generate a robust flow oriented from anterior to poste rior. This system has proven useful for dissecting the role that planar cell polarity (PCP) signaling and hydrodynamic forces play in orienting motile cilia (Mitchell et al., 2007(Mitchell et al., , 2009Park et al., 2008).
Individual multiciliated cells have dozens of cilia that beat in a directed manner termed rotational polarity (Wallingford, 2010). In Xenopus embryos, cilia rotational polarity is estab lished during a discrete developmental window (stages 23-29). First, non-cellautonomous PCP signaling cues instruct the ori entation of ciliated cells, as reflected by a bias in the orientation of cilia beating (Park et al., 2008;Mitchell et al., 2009). This bias initiates a weak flow that establishes a positive feedback loop in which cilia respond to the prevailing hydrodynamic forces by refining their orientation until they are precisely aligned (Mitchell et al., 2007). In the ciliated epithelium of the mouse ependyma, this response to flow requires the PCP signal ing molecule Van Gogh-like 2 and presumably acts through downstream mechanisms in common with the PCP pathway (Guirao et al., 2010).
Cytoskeletal rearrangements have been well characterized to be involved in generating cell polarity (Eaton, 1997;Fanto and McNeill, 2004;Vladar et al., 2009). Asymmetric accumula tion of Frizzled in Drosophila melanogaster wing cells is driven by differential trafficking along microtubules, indicating that polarized microtubules are an integral component of early steps in cell polarity (Shimada et al., 2006;Harumoto et al., 2010). Additionally, RhoAdriven modulation of actin is downstream of PCP signaling in a variety of cell polarity events (Fanto et al., 2000;Strutt, 2001;Marlow et al., 2002). In multiciliated cells, Disheveled (Dvl) and active RhoA accumulate at the base of cilia and are required to coordinate cilia orientation (Park et al., 2008;Mitchell et al., 2009;Hirota et al., 2010). Detailed EM analysis has identified interactions between basal bodies and both actin and microtubules at the apical surface of multiciliated P lanar cell polarization represents the ability of cells to orient within the plane of a tissue orthogonal to the apical basal axis. The proper polarized function of multiciliated cells requires the coordination of cilia spacing and cilia polarity as well as the timing of cilia beating during metachronal synchrony. The planar cell polarity pathway and hydrodynamic forces have been shown to instruct cilia polarity. In this paper, we show how intracellular effectors interpret polarity to organize cellular morphology in accordance with asymmetric cellular function. We observe that both cellular actin and microtubule networks undergo drastic reorganization, providing differential roles during the polarized organization of cilia. Using computational angular correlation analysis of cilia orientation, we report a graded cellular organization downstream of cell polarity cues. Actin dynamics are required for proper cilia spacing, global coordination of cilia polarity, and coordination of metachronic cilia beating, whereas cytoplasmic microtubule dynamics are required for local coordination of polarity between neighboring cilia. intact (Fig. 2 A9). In contrast, the subapical pool of actin, bridg ing neighboring cilia, fails to organize in a polarized manner and is almost completely absent (Fig. 2 A).
To determine whether this subapical pool of actin is involved in orienting cilia, we treated embryos with Cyto D during the developmental window in which cilia orientation is established (stages 23-29) and determined flow velocity, beat frequency, and cilia orientation. Cyto D treatment reduces flow velocity but does not affect either gross cilia morphology or ciliary beat fre quency (Fig. 2, A, D, and E). Consistent with a reduction in flow velocity, Cyto D treatment completely blocks the establishment of coordinated cilia polarity, as reflected by a high circular SD (CSD) of 73 and a short mean vector length (vector length = 1/variance; cells (Gordon, 1982;Lemullois et al., 1988;Sandoz et al., 1988;Chailley et al., 1989a,b). Additionally, actin and actin interacting proteins have been implicated in basal body apical migration and docking as well as ciliogenesis (BoisvieuxUlrich et al., 1990;Pan et al., 2007;Bershteyn et al., 2010;Kim et al., 2010;Ravanelli and Klingensmith, 2011). Here, we focus on the cellular events occurring after ciliogenesis to determine whether cytoskeletal interactions are actively involved in the polarized organization of cilia within a cell.

Results and discussion
Two distinct pools of actin form at the apical surface of multiciliated cells Cilia nucleate from a ninefold symmetric microtubulebased structure termed the basal body, which can be reliably marked with fluorescently tagged centrin. Projecting asymmetrically off of the basal body is the basal foot, which projects in the direction of cilia beating, and the striated rootlet, which projects in the opposite direction (Steinman, 1968;Park et al., 2008). Fluores cently tagged CLAMP has been previously proposed to mark the striated rootlet because of the fact that CLAMP localization proj ects away from the basal body toward the nucleus in the opposite direction of cilia beating. We first analyzed the localization of actin at the base of cilia in relation to both centrin and CLAMP and observed two distinct but interconnected pools of actin at the apical surface of ciliated cells (Fig. 1, A and B). As pre viously reported, a meshworklike actin cap forms at the apical membrane in the same plane as basal bodies such that the ciliary axoneme emerges from the cell through holes in this meshwork ( Fig. 1 A; Pan et al., 2007). Additionally, we observed a sec ond organized pool of actin that localizes subapically, 0.5 µm below the cell surface ( Fig. 1 B). In polarized late cells (stage 29), this second pool of actin consists of short filaments that align to connect neighboring cilia ( Fig. 1 C). These actin filaments link the basal body of one cilium and the distal tip of the striated root let of the posterior neighboring cilium along the axis of cilia ori entation ( Fig. 1 E and Video 1). Early cells (stage 23) that have not established rotational polarity still have an organized apical actin meshwork; however, the subapical pool of actin is largely absent and disorganized ( Fig. 1 D). Thus, this second pool con nects basal bodies and rootlets to create a highly organized matrix reflecting the polarization state of the cell (Fig. 1 F). This sug gests that distinct pools of actin are differentially regulated and may play distinct roles in ciliated cells.

Subapical actin is essential for the establishment of global but not local coordination of cilia polarity
The stability of actin filaments varies considerably depend ing on cellular function, allowing for differential targeting of distinct pools of actin with drugs that affect stability (Seagull and Heath, 1980;Gottlieb et al., 1993). We performed a dose response analysis of the actindepolymerizing drug cytochala sin D (Cyto D) and identified a dose (10 µM) that specifically targeted the subapical pool of actin. In Cyto D-treated embryos, we observed that the apical actin meshwork appears largely Importantly, Cyto D-treated cells maintain their initial polarized bias, indicating that overall cell polarity is not lost but rather that there is a disruption in downstream cellular organiza tion (arrow direction in Fig. 3 C). In an effort to quantify this cel lular organization, we developed software that calculates angular correlation of cilia polarity as a function of distance between each basal body pair in the cell. The comparison of nearest neighbor angular correlations with farthest neighbor angular correlations allows us to quantify cellular organization and to distinguish between local coordination and global cellwide coordination. Interestingly, whereas cellwide coordination is completely lost, Cyto D-treated embryos maintain local coordination within individual basal body clusters, as indicated by a high nearest . This CSD is significantly higher than polarized wildtype late cells (CSD = 20, P < 0.001) but is similar to un polarized wildtype early cells (CSD = 79).
In addition to the loss of cilia polarity, we observed that cells treated with Cyto D have defective basal body spacing. Instead of distributing uniformly along the cell surface, basal bodies form clusters or long filamentlike structures separated by large areas devoid of basal bodies (Fig. 3 C). We quanti fied this defect by scoring the distance of each individual basal body relative to its nearest neighbor (Fig. 3 F). We calculated the mean basal body nearest neighbor distribution for wildtype cells to be 0.97 µm, which is significantly decreased to 0.77 µm (P < 0.001) when subapical actin polymerization is perturbed. foot (Gordon, 1982;Lemullois et al., 1988;Sandoz et al., 1988;Chailley et al., 1989a). However, the function of these inter actions and the consequence of interfering with them remain un known. Because loss of actin dynamics cannot fully explain the generation of refined rotational polarity in multiciliated cells, we expanded our analysis to include microtubules. To test whether microtubule dynamics are actively involved in establishing rota tional polarity, we treated embryos with various concentrations of the microtubuledepolymerizing drug nocodazole (Noc). We determined that treatment with 1 µM Noc between stages 23 and 29 did not affect either gross ciliary morphology (Fig. 2, B and C) or beat frequency (Fig. 2 E). However, using the plusend microtubuletracking protein CLIP170, we determined that this neighbor correlation (Fig. 3 G). This is in stark contrast to early wildtype cells that have a low level of local coordination (Fig. 3 G). This result suggests that cells with perturbed actin dy namics are still able to locally coordinate cilia orientation (high nearest neighbor angular correlation) but fail to integrate polar ity information over the entire cell (high CSD and low farthest neighbor angular correlation).
Microtubule dynamics are essential for the establishment of local coordination of cilia polarity EM studies have described elaborate interactions between basal bodies and cytoplasmic microtubules, specifically at the basal in the presence of multiple microtubuleorganizing centers, we visualized microtubules with a triple GFP-tagged version of the ensconsin microtubulebinding domain (EMTB; Bulinski et al., 2001). EMTB3XGFP localizes to a highly organized micro tubule network that forms within the plane of basal bodies in multiciliated cells (Fig. 4). Microtubules converge in bright foci in close proximity to individual basal bodies (Fig. 4, A and C-D) at the proximal end of the striated rootlet (Fig. 4, B-D). This is consistent with EM data showing microtubules enriched at the basal foot, a basal body appendage that projects 180° opposite to the striated rootlet Park et al., 2008). In late polarized cells, microtubules consistently converge at the posterior ventral side of the nearest basal body aligned with the direction of cilia beating (Fig. 4 D, CSD = 28). In early cells, microtubules still converge in close proximity to basal bodies, but the position relative to the anterior-posterior axis is random ized (Fig. 4 C, CSD = 76). The CSD values are similar to those obtained using centrin/CLAMP (Fig. 3, A and B), indicating that concentration of Noc was sufficient to inhibit cytoplasmic micro tubule dynamics (Fig. 2, B' and C'; Perez et al., 1999).
Furthermore, Noctreated embryos fail to establish robust flow but maintain a weak directed flow, similar to wildtype early embryos (Fig. 2 D). Most notably, these embryos also fail to estab lish rotational polarity (CSD of 85, P < 0.001; Fig. 3, D and E). In contrast to Cyto D, Noc treatment does not affect basal body distri bution on the cell surface (mean distance = 0.96 µm; Fig. 3 F) but instead prevents even local coordination of cilia orientation, as in dicated by a low nearest neighbor angular correlation (Fig. 3 G). Collectively, these data suggest that both cytoplasmic actin and microtubules contribute differentially to the establishment of rota tional polarity. Specifically, actin is required for proper basal body spacing and global cellular organization, whereas microtubules are required for local cellular organization.
Basal bodies are equivalent to centrioles, which serve as the primary microtubuleorganizing center of the cell. To de termine how ciliated cells organize cytoplasmic microtubules led us to hypothesize that these links serve as a physical lever that facilitates propagation of the metachronal wave (Fig. 5 D). To test this possibility, we treated polarized stage 29 embryos with 10 µM Cyto D, which specifically perturbs the actin bridges (Fig. 5 A) but does not perturb the already established cilia ori entation (Fig. 5 B), indicating that the subapical actin pool is not involved in maintaining cilia polarity. To test whether these bridges were involved in propagating the metachronal wave, we captured highspeed videos of mature ciliated cells either mock treated or treated with Cyto D (Videos 2 and 3). Metachrony can be ob served in kymographs taken above multiciliated cells that capture the passing of the effective stroke from multiple cilia as well as a pause during the recovery stroke, resulting in a discrete pattern of lines and pauses (Fig. 5 C). Significantly, in mocktreated em bryos, ciliated cells maintained a highly coordinated beat pattern (Fig. 5 C and Video 2). In contrast, Cyto D-treated embryos lost coordination of cilia beating (Fig. 5 C and Video 3), indicating a loss of the metachronal wave. Because polarity and beat fre quency are unperturbed in these embryos, we propose that this defect is caused by a loss in the physical connection linking the basal body of one cilium to the rootlet of its posterior neighbor.
The high degree of polar, spatial, and temporal coordina tion required for multiciliated cells to generate robust directed flow offers a unique view into the complexity of cell polarity. We are able to use this complexity to separate out the roles that actin and microtubules have in translating polarity cues into cel lular organization. We show that both actin and microtubules form elaborate networks at the apical surface of multiciliated cells linking neighboring cilia via their basal bodies and rootlets (Fig. 5 D). Organization of these networks reflects the polariza tion state of the cell, and, most notably, perturbation of either of these networks leads to defects in cilia polarity. Interestingly, in the absence of cytoplasmic cytoskeletal dynamics, cilia are still capable of beating with normal beat frequencies but are unable to translate the generated hydrodynamic forces into rotational polarity. Importantly, these cells maintain their polarized bias, indicating that these defects are downstream of PCP signaling. Consistent with this, disruption of PCP cell-cell signaling leads to a loss of cell polarity but not cilia organization, whereas dis ruption of the PCP effector Dvl leads to a loss in both cell polar ity and cilia organization (Mitchell et al., 2009).
Our data indicate that in multiciliated cells, there is an in tricate interplay between actin and microtubules that is required to coordinate cell polarity. Microtubules link individual basal bodies to their neighbors via connections at their basal foot (op posite the striated rootlet) and promote local coordination of cilia orientation. In contrast, subapical actin links basal bodies to the striated rootlets of their nearest neighbor along the axis of cilia beating. These connections appear to serve multiple func tions. Cells without subapical actin maintain a significant level of cellular organization, but that organization is no longer in formed by global polarity cues.
The mouse ependyma has a distinct component of cell polarity termed translational polarity, in which the cilia are asymmetrically enriched on one side of the cell. This relies on signaling through a primary cilium and occurs in a myosin IIdependent manner Mirzadeh et al., 2010). microtubule organization reflects the polarity state of the cell, which is consistent with their role in coordinating cilia polarity.

Actin bridges are required for intracellular propagation of the metachronal wave
To generate robust flow, cilia beating must be coordinated not only in direction but also in timing. Cilia at the front edge of the cell initiate the effective stroke first and are subsequently fol lowed by more posterior cilia in a process called metachronal synchrony. The physical connection of neighboring cilia via short actin filaments linking their basal bodies to the posterior neigh boring striated rootlet along the anterior-posterior axis of the cell using a 60× oil Plan-Apo objective with a 1.4 NA. Beating of fluorescently labeled cilia was imaged using resonance-scanning confocal microscopy at 240 frames per second over a 5-s time period. Images were acquired using Elements software (Nikon). To determine flow velocity, we used a digital camera (165FC and DFC295; Leica) to visualize the displacement of 10 µm FluoSpheres polystyrene microspheres (yellow-green fluorescent 505/515; Invitrogen) along the skin of live embryos. Flow videos were acquired using Application Suite software (Leica), and flow velocities were scored using Elements software. High-speed videos of cilia beating were captured at 6,688 frames per second on a Phantom camera (7.2; Vision Research) mounted on a microscope (BX51; Olympus) with a 100× objective. All microscopy was performed at room temperature.
Image processing and quantifications Images were deconvolved using AutoQuant Blind deconvolution in Elements software. Cilia orientation was scored as previously described (Mitchell et al., 2007(Mitchell et al., , 2009Park et al., 2008). In brief, individual rootlet orientation was scored manually by measuring the angle of orientation of the rootlets relative to the anterior-posterior axis of the embryo using Elements software. Excel (Microsoft) and Oriana 2.0 software were used for graphical representation of the datasets and to calculate mean vector length and CSD. Cilia beat frequency was scored using kymographs to capture the effective stroke of an individual cilium over a 5-s time period. Basal body spacing was calculated by determining the absolute position of individual basal bodies using Elements software and by calculating the distance of each basal body to its nearest neighbor using a custom macro in Excel.

Angular correlation analysis
Angular correlation analysis was performed using a custom program. After determining absolute coordinates for basal bodies and distal tips of associated rootlets, unit vectors (length = 1) were constructed that pointed from the basal body to the rootlet tip of each cilia in a ciliated cell. Let A = (A x , A y ) and B = (B x , B y ) denote two such vectors. The angular correlation is the dot product A • B = A x B x + A y B y = cos(), in which  is the angle between the vectors if their tails (basal bodies) were at the same point. The angular correlation (dot product) varies between 1 and +1. The distance between the basal bodies was computed and gives the separation.
Online supplemental material Video 1 shows complex interactions between actin, basal bodies, and rootlets. Video 2 shows metachronal synchrony of cilia beating. Video 3 shows loss of metachronal synchrony after treatment with Cyto D. Online supplemental material is available at http://www.jcb.org/cgi/content/full/ jcb.201106110/DC1.
We give thanks to Lance Davidson and Bill Bement for reagents and to Jennifer Stubbs, Chad Pearson, Ezster Vladar, and John Wallingford for critical reading and discussions. This work was supported by grants from the Parker B. Francis Foundation to B.J. Mitchell, the University of California Irvine Center for Complex Biological Systems, and the National Institutes of Health/National Institute of General Medical Sciences to B.J. Mitchell and C.C. Yu. In contrast, Xenopus multiciliated cells, similar to the mamma lian respiratory tract, lack translational polarity and have cilia that are spaced evenly over the apical surface of the cell. Our data suggest that this spacing requires interactions between neighboring cilia via subapical actin filaments. Surprisingly, interactions between neighboring basal bodies via cytoplasmic microtubules appear dispensable for proper spacing.
Finally, subapical actin connections contribute to the syn chronization of cilia beating. This synchrony has been shown to require both mechanosensory regulation and a conformational switch in the outer dynein arms (Rompolas et al., 2010). Con siderable modeling efforts have suggested that this coordination is derived from hydrodynamic coupling (Gueron et al., 1997;Guirao and Joanny, 2007;Yang et al., 2008). We hypothesize that, in addition to purely hydrodynamic forces, physical links in the form of actin bridges propagate a mechanical signal that is required for the generation of a metachronal wave.

Plasmids and mRNA
The basal body marker centrin (UniGene accession no. Xl.50473) was cloned into pCS2+ and was C-terminally fused with either GFP or RFP, and the rootlet marker CLAMP (UniGene accession no. Xl.26316) was cloned into pCS2+ and was N-terminally fused with GFP, RFP, or mCherry (Park et al., 2008). Cilia were labeled with GFP-Xl.16654, which was amplified by PCR from Xenopus cDNA using the primers 59-GCCTCGA-GATGTTCAAAAACACCTTCCA-39 and 59-GCTCTAGATTGCTTGGCTTT-ATTCTGTA-39 and cloned into pCS2-GFPN3 with XhoI and XbaI. To make pCS2-tag BFPN1, TagBFP (Evogen, Inc.) was amplified by PCR, cut with BamHI and BglII, and cloned into the BamHI site of pCS2+. BFP-centrin was generated by amplifying centrin from pCS2-centrin-RFP and by cloning into pCS2-tag BFPN1. Plasmids were linearized with NotI, and mRNA was generated by in vitro transcription using SP6 RNA polymerase as previously described (Sive et al., 2000). The pCS2-EMTB-3XGFP was a gift from W. Bement (University of Wisconsin, Madison, WI), and the pCS2-CLIP170 was a gift from L. Davidson (University of Pittsburgh, Pittsburgh, PA). mRNA injection and drug treatments Xenopus embryos were obtained by in vitro fertilization using standard protocols (Sive et al., 2000) approved by the Northwestern University Institutional Animal Care and Use Committee. Purified mRNA was injected into two-or four-cell stage embryos into all blastomeres. Embryos were incubated at 16-25°C in 0.1× Marc's Modified Ringer's with gentamicin until they reached the desired developmental stage. Noc (Sigma-Aldrich) and Cyto D (Sigma-Aldrich) were used at final concentrations of 1 µM and 10 µM, respectively. Optimal concentrations were determined from dilution series using the lowest drug concentration that showed significant effects on actin or microtubule organization. Final dilutions were obtained by diluting 1,000× stocks in 0.1× Marc's Modified Ringer's with gentamicin. When embryos were treated at stage 23, they were incubated in the presence of the drug overnight at 16°C until they reached stage 29. When embryos were treated at stage 29, they were incubated in the presence of the drug for 12 h at 16°C.

Microscopy
All immunofluorescence, as well as measurements of cilia beat frequency, was performed on a laser-scanning confocal microscope (A1R; Nikon)