Cytoskeleton polarity is essential in determining orientational order in basal bodies of multi-ciliated cells
11 Cytoskeleton polarity is essential in determining orientational order in basal bodies of multi-ciliated cells
Toshinori Namba and Shuji Ishihara Graduate School of Arts and Sciences, The University of Tokyo, Komaba, 153-8902 Tokyo, Japan; Universal Biology Institute, The University of Tokyo, Komaba, 153-8902 Tokyo, Japan; * Corresponding author E-mail: [email protected].
Abstract
In multi-ciliated cells, directed and synchronous ciliary beating in the apical membrane occurs through appropriate configuration of basal bodies (BBs, roots of cilia). Although it has been experimentally shown that the position and orientation of BBs are coordinated by apical cytoskeletons (CSKs), such as microtubules (MTs), and planar cell polarity (PCP), the underlying mechanism for achieving the patterning of BBs is not yet understood. In this study, we propose that polarity in bundles of apical MTs play a crucial role in the patterning of BBs. First, the necessity of the polarity was discussed by theoretical consideration on the symmetry of the system. The existence of the polarity was investigated by measuring relative angles between the MTs and BBs using published experimental data. Next, a mathematical model for BB patterning was derived by combining the polarity and self-organizational ability of CSKs. In the model, BBs were treated as finite-size particles in the medium of CSKs and excluded volume effects between BBs and CSKs were taken into account. The model reproduces the various experimental observations, including normal and drug-treated phenotypes. Our model with polarity provides a coherent and testable mechanism for apical BB pattern formation. We have also discussed the implication of our study on cell chirality.
Author Summary
Synchronous and directed ciliary beating in trachea allows transport and ejection of virus and dust from the body. This ciliary function depends on the coordinated configuration of basal bodies (root of cilia) in apical cell membrane. However, the mechanism for their formation remains unknown. In this study, we show that the polarity in apical microtubule bundles plays a significant role in the organization of basal bodies. A mathematical model incorporating polarity has been formulated which provides a coherent explanation and is able to reproduce experimental observations. We have clarified both necessity (‘why polarity is required for pattern formation’) and sufficiency (‘how polarity works for pattern formation’) of cytoskeleton polarity for correct pattering of basal bodies with verification by experimental data. This model further leads us to a possible mechanism for cellular chirality.
Introduction
Synchronous and coordinated ciliary beating underlies various biological functions in organisms ranging from ciliates to mammals [1]. In the mammalian trachea, synchronous ciliary beating of multi-ciliated cells produces mucociliary transport, that is the flow of mucus to the oral side of the epithelial cell surface, ejecting viruses and dust to the outside of the body [2]. For efficient mucociliary transport, hundreds of cilia covering the apical membrane of a cell beat synchronously in the same direction. This organized directionality in ciliary beating is achieved by basal bodies (BBs), which are roots of the cilia. These BBs point in the same orientation (Fig 1A), which is defined by the relative position of the basal foot (BF), a major appendage associated with BBs. In addition, previous studies have revealed that BBs are regularly aligned in the apical membrane of the mature epithelial cells in mouse trachea [3]. The positional and orientational order of BBs is crucial for proper ciliary beating, as reported in many other systems where disturbance in ciliary organization resulted in dysfunction of organs, such as inner ear hearing loss [4], infertility [5], hydrocephalus [6], and situs inversus [7]. Therefore, elucidation of the mechanism of ciliary organization on the apical membrane is of paramount biological and clinical importance. However, it is largely unknown how BBs are coordinated through the maturation of the tracheal multi-ciliated cells.
Fig 1. Schematics of the positional and orientational orders of BBs in the trachea. (A) The mammalian trachea sends mucus to the mouth on the surface of tracheal epithelial cells (red arrow). Tracheal epithelia is composed of multi-ciliated (green) and mucous cells (red). The mucus flow is produced by beating of multiple cilia in the apical membrane of the multi-ciliated cells. The orientation and positioning of cilia are determined by BF (red corn) and BBs (green rod), respectively. BF is an MT-organizing center [13,41]. (B) Typical pattern of BB distribution. During maturation of multi-ciliated cells, the pattern of BBs is changed from cluster to alignment patterns [3]. (C) Polar MTs can orient BBs to the definite direction. (i) Binding of the BF to the neighboring MT bundles is not sufficient to realize unidirectional orientation of BBs [14]. (ii) BBs can orient to same direction if polarity exists in the MT bundles (denoted as p ) and BBs can recognize it. The polarity is also shown by blue arrows in panel A. (iii) TEM images of BBs and detected MTs (adapted from [3]). BBs contact MTs at a biased angle. (D) Histograms of the angle of BFs measured using TEM (From (I) Movie S8 in [9], (II) Fig 6 in [9] (III) Video 6 in [3], and (IV) Video 7 in [3]). Long-term live cell imaging of cultured tracheal cells revealed that maturation of BBs in the apical membrane proceeds though several stages, characterized by a typical pattern in BB distribution (Fig 1B) [3]. Clusters of BBs are observed in the early stages after differentiation to multi-ciliated cells. The clusters are then scattered, followed by formation of short and branched arrays of several BBs termed as ‘partial alignment’ [3]. Through connection and reconnection of these arrays, BBs finally become regularly aligned over the apical membrane. The number of BBs remains nearly unchanged during the maturation period. BBs become increasingly oriented in the same direction indicated by the relative position of the BF, concurrent with the alignment of BBs. Quantification of positional and orientational orders of BBs showed that these two ordering processes in position and orientation were highly correlated [3], suggesting that both are coordinated by common factors. Cytoskeletons (CSK) in the apical membrane were shown to be involved in the regulation of BB distribution. During maturation of cells, microtubules (MTs) and 10 nm-filaments, which are sparsely distributed in the early stage, grow to a dense mesh network and surround the BBs [8]. Most MT bundles run in parallel between the aligned BB array, which might contribute to the positioning formation of the BB merely by invading and separating the array of BBs [9]. Further, since BF directly binds to MTs as confirmed by transmission electron microscopy (TEM) for which -tubulin in BF might mediate the association [41], it was suggested that the orientation of the BBs is also regulated by MTs [9]. Remarkably, inhibition of MT polymerization by nocodazole treatment results in disruption of BB alignment pattern and formation of BB cluster pattern , with disruption in the orientation of BB as well. Removal of nocodazole by washing results in recovery of BB positional order [3]. In summary, there is sufficient experimental evidence to show that CSKs are responsible for the determination of position as well as orientation of BBs. However, the underlying mechanism explaining these observations remains elusive. In particular, the attachment of BF to MTs might be important for the determination of BB orientation [41], but its mechanistic details have not been identified as yet. In the present study, we have generated a mathematical model to explain the ordering of position and orientation in BBs. Our aim here is not to propose a detailed model that can explain all aspects of the experimental findings, but rather to provide a possible model to explore the role of MTs for the observed BB patterning. Firstly, the symmetry of the system has been taken into consideration to show that determination of orientation in BBs requires polarity in MT bundles. The presence of polarity is supported by a quantitative relationship between BB orientation and CSK bundles, determined from available experimental data [3,9]. This was followed by derivation of a mathematical model based on the characteristics of CSKs and their interaction with BBs. ‘Active matter’ formalism was employed by which the generic form of system equations is obtained by taking into account the symmetry of the system [10,11]. This kind of formulation was conducted for studying pattern formation on apical membranes for a system without BBs [12]. Our early study on tracheal cells also employed this scheme and produced positional order of BBs [3], but did not consider the orientation of BBs and the CSK-BB interaction. Using the derived model, we have investigated the conditions for BBs to align and point in the same direction, which further confirms that MT polarity is crucial for the establishment of BB pattern. The information transmission from planar cell polarity (PCP) to the BB orientation has also been discussed. In the Discussion, we have summarized our findings and challenges that need to be addressed in future. We have also considered the implication of our findings in determining a possible mechanism for cellular chirality in determining the relationship between PCP and CSK polarity. Results
Polarity in cytoskeleton guides orientation of BBs
During the maturation period of multi-ciliated cells, BBs establish regularity in their position and orientation. The patterning of BB arrays follows the formation of the CSK matrix, particularly the MT bundles forming on the cytosolic apical membrane [3,9]. The BBs orient in the same direction at the same time [3]. Since the BF, which is an appendage of BBs, has a component binding to MTs as confirmed by TEM images and mutant strains in Odf2 [9,13,41], it was suggested that the connection of BF with neighboring MT bundles is a determinant of BB orientation [14]. However, this binding between BF and MTs allows BBs to associate with either of the two neighboring MT bundles and hence they cannot distinguish between the two opposite directions (Fig 1C (i)). Therefore, how BBs can determine their definite direction remains a mystery. A probable hypothesis to fill this knowledge gap is the presence of polarity in the MT bundles. It may be assumed that there is a component in CSK that provides the direction which the BBs are able to interpret. The presence of polarity therefore enables BBs to distinguish the direction by breaking the symmetry in the relationship between BBs and MT bundles (Fig 1C (ii)). To assess this hypothesis, we focused on the contact angles between BBs and MT bundles, since the angles are apparently neither perpendicular nor uniformly distributed but are specifically biased. We measured the angles between MTs and the orientation of the associated BBs from published TEM images (Fig 1C (iii)) [3,9]. The contact angles were biased to the positive side as shown in Fig
1D and Table 1 (n = 4). This observation indicates a polar nature in MTs and chirality in BF, which enables the connection of BBs to MTs within a specific range of angles.
Table 1. Measured angles of BFs and its numbers in Fig 1D.
TEM images mean S.D. numbers Reference I o ±10.5 o
10 Movie S8 in [9] II o ±8.6 o
9 Fig 6 in [9]
III o ±20.6 o
23 Video 6 in [3] IV o ±19.0 o
36 Video 7 in [3] Mean and standard deviation of contact angle of BBs to MTs. Each sample is adopted from TEM images shown in Reference. Mean and standard deviation are defined according to directional statistics [15].
Mathematical model for the interaction between CSK and BB
A mathematical model was generated based on the polarity of MTs, demonstrating that polarity is a necessary and sufficient condition for the reproduction of the observed regular positional and orientational orders of BBs on the apical membrane. The following processes were taken into account while generating the model (Fig 2); (i) CSKs form bundles on the apical membrane that sustain the intrinsic polarity. (ii) CSKs and BBs exclude each other. (iii) BBs repel each other. (iv) Polymerization and depolymerization of CSKs. (v) BBs orient in the direction determined by CSK polarity. The processes (i)-(iii) are modeled by a free energy function F , a generic form of which is as follows [10-12,16]. 𝐹 = 𝐹 𝑐𝑠𝑘 + 𝐹 𝑐𝜙 + 𝐹 𝜙𝜙 (1) 𝐹 𝑐𝑠𝑘 = ∫ [𝑓 𝑐 (𝑐) + 𝐷 [ ∇𝑐 ] 𝑡 ( 𝐼 + 𝛼 𝑄 𝑄 ) ∇𝑐 − β P ( 𝑐 ) | 𝒑 | + | 𝒑 | + 𝐾2 | ∇𝒑 | ] 𝑑𝑟 (2) 𝐹 𝑐𝜙 = 𝜆 𝑐𝜙 ∑ ∫ 𝜙 𝑖 𝑐 𝑑𝑟 𝑖 (3) 𝐹 𝜙𝜙 = 𝜆 𝜙𝜙 ∑ ∫ 𝜙 𝑖 𝜙 𝑗 𝑑𝑟 [𝑖,𝑗] (4) Here, c ( t , x ) and p ( t , x ) are continuum field variables representing concentration and polarity of CSKs at time t in a two-dimensional space x . 𝐹 𝑐𝑠𝑘 is the free energy of bundling and polarity of CSK. We employed a simple Landau energy form 𝑓 c (𝑐) = −(𝛼 /2)(c − c ) + (𝛼 /4)(𝑐 − 𝑐 ) , according to which, MT accumulation is induced to form bundles. In the second term, I is the identity matrix and 𝑄 = 𝒑⨂𝒑 − (1 2⁄ )Tr[𝒑⨂𝒑]𝐼 is a nematic tensor depending on the polarity of CSKs. This term represents ‘surface tension’ of the MT bundles, in which, polarity p prefers parallel alignment to the CSK concentration boundary with a positive value of 𝛼 𝑄 , i.e., p is perpendicular to the cell concentration gradient ∇𝑐 (Fig 2A(i)). The term 𝐹 𝑐𝑠𝑘 describes the formation of CSK polarity. By setting 𝛽 𝑃 (𝑐) = 𝛽(𝑐 − 𝑐 𝑐 ) , the polarity increases with MT concentration exceeding the critical value c c , up to an amplitude of |𝒑| = √𝛽 𝑃 (𝑐) . The last term represents the tendency of the polarity to align with each other, where one constant approximation is employed. Taken together, F csk provides the free energy that expresses process (i), where MT concentration and polarity are coupled via nematic tensor Q . The simulation results without BBs has been shown in the methods section and S1 Fig. Fig 2. Model for positional and orientational orders of BBs. (A) Illustration of the processes considered in the model. (i) CSK filaments are condensed and bundled in parallel and their polarity orient in the same direction (blue arrow, p ). (ii, iii) One BB repels the CSK bundle and other BBs. (iv) Polymerization and depolymerization of CSKs. (v) The orientation of BB, 𝜃 , is determined by CSK polarity. (B) Typical time series of BB pattern development ( 𝛼 𝑄 = 2.0 ). Green indicates CSK concentration, and blue arrow indicates the polarity of CSK. Gray disk is a BB and the line segment in the disk indicates the orientation of BB. Patten of BBs reaches “alignment”. (C) Time series when the coupling between polarity and CSK concentration is weak ( 𝛼 𝑄 = 0.2 ). Eventual pattern is “partial alignment”. Bars, 0.5 µm. 𝐹 𝑐𝜙 (Eq. 3) and 𝐹 𝜙𝜙 (Eq. 4) are the free energies representing processes (ii) and (iii), respectively. The variable 𝜙 i (𝑡, 𝒙) represents the region occupied by the i -th BB, a circular disk with radius R centered at X i . We set 𝜙 i = [1 + tanh(𝑣(𝑅 − |𝒓 − 𝑿 𝒊 |))]] 2⁄ by which 𝜙 i = 1 for the internal region of the BB, and 𝜙 i = 0 otherwise [17]. 𝑣 controls the steepness of the boundary of the BB surface and is set to be sufficiently large. 𝐹 𝑐𝜙 is the energy term for the excluded volume effect between BBs and CSKs, which was not considered in the previous model [3]. 𝐹 𝜙𝜙 is the energy term representing the excluded volume effect among BBs, to avoid overlap between different BBs, summed over all pairs of BBs. The dynamics of the system is derived from the free energy F as follows. 𝜕 𝑡 𝑐 = 𝛼 𝐽 ∇ ( 𝛿𝐹𝛿𝑐 ) + 𝑘 𝑝 𝑐𝑐+𝐾 𝑀 − 𝑘 𝑑 𝑐 (5) 𝜕 𝑡 𝒑 = −𝛾 𝑃 𝛿𝐹𝛿𝒑 (6) 𝜕 𝑡 𝑿 𝒊 = −𝛾 𝐵𝐵 𝜕𝐹𝜕𝑿 𝒊 + 𝜉 𝑖 (7) In Eq. (5), time evolution of CSK concentration c is derived from the energy relaxation with mass conservation, where 𝛼 𝐽 is the kinetic coefficient. Polymerization and depolymerization (process iv) are represented by 𝑘 𝑝 and 𝑘 𝑑 being maximum polymerization rate and depolymerization rate, respectively, with 𝐾 𝑀 being a constant. The mean concentration of CSKs in the system is determined by the balance between these rates. We assumed 𝑘 𝑝 is rescaled by 𝑘 𝑑 to keep the MT concentration almost constant (see Table 2). The time evolution rule of polarity, 𝒑 , and the center of the i -th BB, 𝑿 𝑖 , are determined by simple relaxation of the free energy given by the derivative of F . In Eq. (7), the center of the i -th BB moves by repulsive interaction with CSKs and other BBs. Stochastic noise is added to the evolution of X i , where 𝜉 𝑖 is a white Gaussian random variable obtained from the normal distribution with zero-mean and standard deviation 𝜎 𝐵𝐵 Δ𝑡 at each time step, where Δ𝑡 is the discretized time step in the simulation. Finally, we assume that the polarity of CSKs guides BB orientation via the interaction between CSKs and BFs (process v). Although there is a preferred angle between the direction of CSK polarity and orientation of a BB (Fig 2A(v)), the preferred angle has been assumed to be 𝜋 2⁄ in this model. This assumption does not change the generality of the model. ∂ t 𝜃 𝑖 = −γ BF (〈𝑐 𝒑〉 𝜙 𝑖 ) ⋅ 𝒏 𝑖 + 𝜉 𝐵𝐹 (8) Here, 𝜃 𝑖 is the orientation of the i -th BB and 𝒏 𝒊 = (cos𝜃 𝑖 , sin𝜃 𝑖 ) is its director. The angle brackets 〈⋯ 〉 denote the average over the surface of the BB. Gaussian white noise term 𝜉 𝐵𝐹 is also added with standard deviation 𝜎 𝐵𝐹 Δ𝑡 . The above argument provides the closed form of the equations to be explored below. We have numerically solved non-dimensionalized equations where the units of time and length in the simulations correspond to 5 min and 0.04 m, respectively. Details of the equations, parameter values, and numerical methods have been provided in the methods section. The mathematical model reproduces BB positioning and orientation
Depending on the parameter values used in the model equations, several patterns were observed in the numerical simulations. Fig 2B shows an example of the time series obtained in the simulation when 𝛼 𝑄 , the coupling strength between the polarity and concentration of CSK in Eq. (2), is set at 𝛼 𝑄 = 2.0 (see Table 2 for the other parameters). The polarity of CSK, p, becomes oriented in the same direction as the system and the CSK concentration, c , becomes striped along the direction of the polarity. BB arrays are formed and aligned between the stripes of CSK concentration, with all BBs being oriented in the same direction as the polarity. This pattern obtained in the numerical simulation recapitulates the experimental observations in normal conditions. We refer to this pattern as ‘alignment’. Decreasing α Q retains direction of the polarity with stripes of CSK concentration often becoming branched and BB arrays becoming winding (Fig 2C, 𝛼 𝑄 = 0.2 ), which is referred to as ‘partial alignment’ [3]. This is due to the weak coupling between polarity and CSK concentration, suggesting the importance of polarity in establishing an alignment pattern. This will be further studied in the next section. Table 2. Parameters used in Fig 2B. * 𝑘 𝑝 is kept to be proportional to 𝑘 𝑑 for the MT concentration to be almost constant. Parameters
Free energy 𝛼 𝛼 𝑐 𝐷 α Q 𝒑 𝛽 Contribution from concentration to polarity (Active term) 𝑐 𝑐 Threshold of CSK concentration at which polarity can grow 𝐾 𝜆 𝑐𝜙 𝜆 𝜙𝜙 𝑣 𝑅 Dynamics 𝛼 𝐽 𝑘 𝑝 𝑘 𝑑 (𝑐 + 𝐾 𝑀 ) Rate of polymerization of CSK * 𝑘 𝑑 √0.1 Rate of depolymerization of CSK 𝐾 𝑀 𝛾 𝑃 𝛾 𝐵𝐵 𝛾 𝐵𝐹 Noise strength 𝜎 𝐵𝐵 Standard deviation of white Gaussian random variable 𝜉 𝑖 𝜎 𝐵𝐹 Standard deviation of white Gaussian random variable 𝜉 𝐵𝐹 System Parameters N L x x -axis (same as above) L y y -axis (same as above) dx Fig 3 illustrates the phase diagram of the patterns observed by the numerical simulation against 𝛼 𝑄 and the polymerization/depolymerization rate 𝑘 𝑑 . Some typical patterns are shown in the diagrams which are dependent on the parameter values. The boundaries between these patterns are not so sharp due to the initial condition dependence and difficulty in classification, particularly due to the appearance of mixed patterns. For weak polymerization rate ( 𝑘 𝑑 < 10 −1.5 ), clusters of BBs appear regardless of the value of α Q . By increasing α Q , double strand arrays of BBs are obtained. Such a pattern is occasionally observed in cultured tracheal cells [3]. On the other hand, BBs are scattered for sufficiently high polymerization rate ( 𝑘 d > 10 ) since the amounts of CSKs surrounding BBs recover quickly to the mean concentration and isolate them. Alignment and partial alignment patterns appear in the middle region of k d , depending on the value of α Q as mentioned above. Altogether, configuration of BBs is dependent on the polymerization rate and coupling strength between CSK polarity and concentration. Appearance of BB clusters at low polymerization rate is of interest since it resembles those observed in the nocodazole treated clusters [3], which will be discussed further. Fig 3. Phase diagram of BB pattern against parameters 𝛂 𝐐 and 𝒌 𝐝 . α Q is the coupling constant between CSK concentration and polarity. 𝑘 d is the polymerization/depolymerization rate of CSK. Patterns are classified into Alignment (A), Cluster (C), Scatter (S), partial alignment (PA), and Double strand (DS). Typical BBs pattern in each region is shown are also shown. Bars, 0.5 µm. Polarity is required for alignment pattern of BBs
To further explore the role of polarity in CSKs, we studied the change in the behavior of the system with change in the parameters α Q and K , the latter is bending modulus of MT bundles (Eq. 2). Our numerical simulation suggests that these two parameters must be sufficiently high for achieving regular alignment pattern of BBs, indicating importance of the polarity in CSKs for BB alignment pattern. In S2 Fig, the straightness of BB arrays is quantified by calculating 〈−cos𝜓 𝑖 〉 , which is the average of the negative cosine of the opening angle between line segments connecting the i -th BB and its two neighboring BBs. When either α Q and K are low, arrays of BBs are short or are not aligned straight, classifying them into ‘partial alignment’ or a ‘scatter pattern’. Therefore, polarity in CSKs plays two roles in our model. It does not only convey the orientational information to BBs (Eq. 8), but also enables BBs to be regularly aligned through the interaction with CSK concentration and BBs. Bundling of CSKs enhances BB alignment formation
MTs are bundled in the apical membrane, which stiffen the apical CSK network and could contribute towards the formation of aligned BB arrays. For more insights into the role of the MT bundles for BB pattern, we changed the functional form of the free energy f c ( c ) in Eq. (2) which drives the formation of MT bundles (condensation of MTs). We re-parametrized the free energy as f c ( c ) = 𝜆 W (𝑐 − 𝑐 − 𝑊) (𝑐 − 𝑐 + 𝑊) in which the parameter 𝜆 W controls the scale of the energy (i.e., ‘depth’ of the potential), while the positions of the two minima are controlled by 𝑊 . With a proper form of f c ( c ), BB arrays are well lined up as shown in S3 Fig, consistent with the preceding sections. Since in the system without BBs, the double minimum form is indispensable for inducing the phase separation (condensation) of CSKs resulting in the stripe pattern (see Methods and S1 Fig), we hypothesize 𝑊 to be positive-finite. If that were not the case, BBs would scatter and be randomly configured. This is confirmed for large λ 𝑊 and small 𝑊 where BBs fails to align (S3 Fig (iv)). However, when both λ 𝑊 and 𝑊 are small, there were unexpected observations in which arrays of BBs were formed and were aligned (S3 Fig (ii)). Since the form of f c is flattened in the parameter region, we then conducted numerical simulation by erasing the free energy (i.e., f c ( c ) = 0). We found that the aligned pattern of BBs is achieved as shown in Fig 4A. The obtained phase diagram on k d - Q plane shown in Fig 4B is similar to those obtained in Fig 3, indicating that the condensation of CSKs is not essential for the pattern of BBs. However, the region of the ‘aligned pattern’ becomes narrower, and the BB arrays are often branched in the aligned region (Fig 4B(i)) and are less straight. As shown in Fig 4C, in the model including f c ( c ) term, BBs align in straighter over a wide range of noise strengths in Eq. (7). Thus, MT bundles are useful for more robust and ordered positioning of BBs. Fig 4. Model without the condensation term and simulation for nocodazole treated cells. (A) Typical time series of BB pattern when 𝑓 c (𝑐) = 0 . Parameter values are the same as Fig 2B. (B) Phase diagram of the model without condensation. (i) and (ii) show examples of “alignment” and “partial alignment”, where the same parameters are used as corresponding points in Fig 3. (C) Robustness of BB for noise strength, σ BB . The straightness of BB arrays, 〈−cos𝜓 𝑖 〉 , was compared between models with and without 𝑓 c (𝑐) term. Gray shadow area indicates σ BB used in the other simulations. The error bars represent the standard error among four independent simulations. (D) Simulated patterns corresponding to nocodazole treated cells, which is produced by the model without condensation term and coupling with polarity (See Methods). Bars, 0.5 µm. We conclude from the above numerical results that even without the condensation effect of MTs, aligned patterns can be realized to some extent by the exclusion volume effect between BBs and CSK. This was in contrast to the results of the previous model, where exclusive interaction between BBs and CSKs was not considered and condensation of CSKs was indispensable for patterning of BBs [3]. However, it is possible that the fluctuation is suppressed and the aligned pattern of BBs is achieved more robustly through the use of MT bundles. This result is biologically relevant since the maturation of multi-ciliated cells must progress under various external disturbances [14,30].
Simulation for nocodazole treated cells
Given the role of polarity and condensation of CSKs investigated in the preceding sections, here we studied the formation of BB clusters in cells treated with nocodazole, which inhibits MT polymerization [3]. As shown in Figs 3 and 4B, BB clusters are formed at low polymerization and depolymerization ( 𝑘 𝑑 < 10 −1.5 ), consistent with the experimental result. To clarify the mechanism of cluster formation, we simplified our model to ignore polarity and condensation of MTs, since they would play only an insignificant role when MT polymerization is inhibited. The simplified model is obtained by replacing the free energy F csk in Eq. (2) with the following one; 𝐹 𝑐𝑠𝑘 = ∫ 𝐷 [ ∇𝑐 ] 𝑑𝑟 . (9) The numerical simulation of the simplified model demonstrates the formation of BB clusters (Fig 4D), similar to the experimental observation in nocodazole treated cells. This result indicates that clustering of BBs occurs by merely minimizing the quadratic term of the concentration gradient, attributed to the depletion force [18]. Since the concentration gradient is localized in the boundary of BBs, accumulation of BBs is preferable for decreasing the total length of the boundary. Meanwhile, interaction length of the resulting attractive force between BBs is too short, resulting in the appearance of multiple clusters depending on the initial position of BBs. MTs transmit directional information in PCP to BBs
All the cilia orient to the oral side of the tracheal tissue in vivo , which is coordinated by PCP. In fact, directionality of BBs is reduced by knocking out the PCP-related protein Vangl1, compared with those in wild-type (WT) [19]. Importance of PCP in cilial organization is also evident in many systems [4,6,7,20]. In
Drosophila epithelial cells, the apical MTs are aligned along the proximal-distal axis (coinciding with the PCP axis), which is involved in directed molecular transportation within the cell [21-23]. Recent experiments suggest that apical CSKs mediate the transmission of directional information from PCP to BBs and cilia, where PCP proteins are localized along the oral and lung side of the adherent membrane in each cell (Fig 5A). However, the detailed information transmission mechanisms remain unclear and may be different among different cell types and species [19,24]. Fig 5. Simulation of BBs positioning by PCP. (A) PCP proteins are localized along the oral and lung boundaries on the apical membrane. CSK concentration in the apical is low in immature cell and becomes increased during maturation. In matured cells, MTs are more distributed on the oral side (indicated by localization of DVI1/3 and Daple). (B) Polymerization rate 𝑘 p is increased (blue line) and the boundary condition is changed at t =66.7 h (gray dashed line) in the simulation. The straightness of BB arrays obtained by simulation is shown by a red line. (C) Simulation result for pattern formation of BB alignment. Black triangle represents change in boundary condition (see text). Bars, 0.5 µm. It was observed that MT concentration is significantly high along the oral side of the cell boundary membrane in a tracheal cell, which was attributed to the interaction with PCP proteins such as Frizzled, Dishevelled, and Daple which accumulate on the oral side of the cell membrane (see Fig 5A) [3,19,25]. Although the asymmetric accumulation of the cortex MTs can also be related with MT organization over the apical membrane, and can thus proceed at the same time with the pattern formation of BBs [42], here we do not address the entire process. We restrict our study to simply demonstrate a possible mechanism of how the pattern of BBs can be modulated by the PCPs. In the framework of our model, the effect of PCP is incorporated as the boundary condition as follows [3]. For the oral side boundary, concentration of CSKs is set to be high (i.e., Dirichlet condition 𝑐 = 𝑐 H is applied, see Methods for details) and direction of the polarity is restricted to be parallel along the border. Neumann conditions are applied for the variables for the other boundary. In order to demonstrate the maturation process of a multi-ciliated cell during high CSK concentration, we simulated a condition where the CSK concentration on the oral boundary is switched from Neuman to Dirichlet boundary conditions mentioned above, at time t = 66.7 h, and the parameter k p also increases with time (Fig 5B blue line; see Methods for details). The time series of the simulation result is shown in Fig 5C. Initially scattered pattern of BBs ( t = 10 h in the panel) evolve to the partial alignment pattern ( t = 50 h and 80 h), where short or branched arrays of BBs are formed by the increase in CSK concentration and BBs orient to the same direction locally (Fig 5C). BBs eventually establish a fully aligned pattern that is parallel to the oral border through connection and reconnection of BB arrays, ( t = 110h). This time evolution of BB pattern agrees with those observed in the normal maturation process of a multi-ciliated cell. Discussion
On the apical membrane of epithelial cells, CSKs such as MTs, actins, and 10 nm-filaments can perform various biological functions through the process of pattern formation attributed to the self-organizational ability of CSKs [12,16,26-28]. In the present study, we investigated positional and orientation orders of BBs in the apical membrane of multi-ciliated cell which are responsible for coordinated ciliary beating to generate mucociliary transport. Apical MTs, which are bound to the BF of BBs, play a crucial role in achieving the positional and orientational orders of BBs. Most importantly, consideration of the symmetry led us to propose that polarity should be maintained in the MT bundle filaments for properly guiding orientation in BBs. This idea is supported by the biased distribution of angles between BB orientation and MT filaments in TEM images (Fig 1D). Although more direct evidence of the presence of polarity in MT bundles is still lacking, this can be investigated and confirmed in future experiments. We propose a mathematical model for patterning of BBs, taking into consideration the polarity of MT bundles. The model was constructed based on assumptions such as CSK characteristics and volume exclusion among BBs and CSKs. Since the fundamental mechanism of how BBs align and orient is not yet known, the model was not aimed at being precise or quantitative, but rather provided a basic description of the processes for understanding the underlying mechanism following the generic formalism [10-12]. The derived model is distinct from the previous one [3] in incorporating polarity as well as the repulsive interaction between BBs and CSKs, by which it can reproduce both the positional and orientational orders of BBs in multi-ciliated cells. Analysis of the model demonstrated that polarity is crucial not only for determination of BB orientation, but also for achieving alignment pattern of BBs, the latter depending on the nematicity of the CSK determined by polarity. It was also shown that BBs can align even without the condensation effect of CSK concentration, i.e., exclusion of BBs and CSKs was sufficient for the formation of alignment pattern. However, the formation of MT bundles enabled robust formation of the alignment pattern against external disturbance. These simulation results are feasible given the experimental observation that MTs are indispensable for establishing alignment pattern in the apical membrane of multi-ciliated cells where MT bundles segment the aligned array of BBs. Positional order of BBs was disrupted by suppressing MT polymerization with BB clusters being formed, which could also be reproduced by the present model. Taken together, it was demonstrated that incorporation of CSK polarity in the model was necessary and sufficient to reproduce prior experimental observations. It should be noted that our model still depends on various assumptions and misses some crucial aspects observed in experiments. First, in our model, the effect of PCP was simplified into a fixed boundary condition (i.e., Dirichlet condition) that represents accumulated tyrosinated MTs on the oral side of cellular cortex [19]. As explained before, this accumulation of the cortical MTs itself can be regulated by both PCP proteins and apical MTs. These processes should be incorporated to study a more precise time course of patterning of BBs. Second, in our model, it was assumed that the orientation of BBs merely follows the direction of MTs (process (v) in the model assumption) for which the detailed mechanism was not taken into account. Since BF contains -tubulin [41], it will be interesting to explore whether nucleation of the apical MTs at BF has a possible impact on the determination of BB orientation to MTs. In addition, while PCP and apical CSKs are indispensable for the coordination of BBs and ciliary beating, it was reported that hydrodynamic interaction or direct physical contacts among cilia were part of a positive feedback loop for accurate BB patterning and contributed to the correct directionality in BBs and ciliary beating to enhance mucociliary transport [24,29,30]. It would be interesting to extend our model to understand how the mechanical perturbation from the synchronously beating cilia to the BBs in apical membranes enables smooth BB rearrangement and promotes their alignment pattern. Finally, we would like to discuss the implication of our study in determining cellular chirality. Presence of polarity in CSKs aligning perpendicular to the PCP axis (oral to lung axis) indicates that the cell can distinguish between left and right directions from the PCP axis (Figs 1A and 5A, blue arrows). From a theoretical point of view, our study provides a possible mechanism for implementing cellular chirality by using PCP and CSK polarity. In the biological system, chirality appears at various levels ranging from molecular [31] and cellular [32] to the tissue level [33]. Cellular chirality is responsible for body formation. For example, cell intrinsic chirality underlies left-right asymmetry of chicken cardiac looping [34,35]. Nodal cilia are well-known examples in which the cilia are tilted in a definite direction and generate the left-ward flow on the nodal tissue surface, which induces the left-right asymmetry of animal body [7,36,37]. Although the mouse trachea studied in the present work has no apparent chiral asymmetry, it is possible that the system can regulate cellular chirality through some mechanisms, which is of potential importance in generating chirality at the tissue level. We believe that our model provides useful insights into exploring the common mechanisms of cellular and tissue chirality. Methods
Mathematical model for BB pattern formation in multi-ciliated cells
Details of the model Eqs. (5-7) in the main text are as follows. ∂ t 𝑐 = α J ∇ (−𝛼 (𝑐 − 𝑐 ) + 𝛼 (𝑐 − 𝑐 ) − 𝐷 (∇ 𝑐 + 𝛼 𝑄 ∇ ⋅ (𝑄∇𝑐)) + 2𝜆 𝑐𝜙 𝑐 ∑ 𝜙 𝑖𝑖 ) +𝑘 𝑝 𝑐𝑐+𝐾 𝑀 − 𝑘 𝑑 𝑐 (10) ∂ t 𝒑 = 𝛾 𝑃 ((𝛽 𝑐 (𝑐) − |𝒑| )𝒑 + 𝐾∇ 𝒑 − α Q 𝐷 𝑀 𝑐 𝒑) (11) ∂ t 𝑿 𝒊 = 2γ BB 𝑣 ∫ [(𝜆 𝑐𝜙 𝑐 + 𝜆 𝜙𝜙 (∑ 𝜙 𝑗𝑗≠𝑖 )) 𝜙 i (1 − 𝜙 i ) 𝑿 𝒊 −𝒓|𝑿 𝒊 −𝒓| ] 𝑑𝑟 + 𝜉 𝑖 (12) Ginzburg-Landau potential 𝑓 𝑐 (𝑐) = −(𝛼
2⁄ )(c − c ) + (𝛼
4⁄ )(𝑐 − 𝑐 ) has been employed for this model which has two minima at 𝑐 ± = 𝑐 ± √ 𝛼 /𝛼 . The last term in Eq. (11) represents the reactive term to align the polarity along the CSK concentration gradient, where 𝑀 𝑐 =(∇𝑐)⨂(∇𝑐) − (1 2⁄ )Tr[(∇𝑐)⨂(∇𝑐)]𝐼 . The reactive term associated with 𝛽 𝑐 (c) was ignored for simplicity in Eq. (10). This was verified by interpreting that the emergence of polarity p is governed by active dynamics. In practice, the results are almost independent of the presence of the reactive term. In Eq. (3), we adopted 𝐹 𝑐𝜙 = 𝜆 𝑐𝜙 𝜙 𝑖 𝑐 for exclusive volume effect between i -th BB and MTs. Alternative symmetrically permissible form is to adopt 𝐹 𝑐𝜙 = 𝜆̃ 𝑐𝜙 𝜙 𝑖 𝑐 , by which we found similar patterns with appropriate choice of 𝜆̃ 𝑐𝜙 and 𝜆 𝜙𝜙 . The dynamics of CSK and BBs was numerically simulated by integrating the differential equations (8 and 10-12). The non-dimensionalized equations of CSK defined on a rectangular domain 𝐿 x × 𝐿 y were calculated. The time integration is done by the fourth-order finite difference scheme with time step ∆𝑡 = 1.0 × 10 −3 . Periodic boundary condition was employed unless stated otherwise. The system size was chosen as L x = L y = 64, for which the domain is discretized by a square grid with ∆𝑥 = 1.0 spacing. The number of BBs is N = 80 for most of the simulations. Numerical simulations were performed with system size L x = L y = 32 and N = 20 for producing the phase diagrams in Figs 3 and 4. Four independent simulations were performed for each set of parameters to check the dependency on initial conditions. Initial conditions were set as follows. The position and orientation of BBs, X i and 𝜃 𝑖 , were random, the polarity of CSK p was set at | p | = 0.1 with random direction, and the concentration of CSK is set at c = 0.2 with small value of additional noise. The values of the parameters used in the simulations are summarized in Table 2. Since exact kinetic values were not determined experimentally, we determined them so as to produce experimental patterns, using the following rationale. The units of length and time were set to be equal to 0.04 m and 5 min, respectively. Radius of a BB was chosen as R = 2.5 corresponding to 0.1 m [3,9]. System size was 𝐿 x × 𝐿 y = 64 × 64 corresponding to , approximately quarter of the apical area of single multi-ciliated cells. N = 80 since ~300 BBs exist per cell. Typical speed of BBs was about 0.02-0.03 m/ 5min [3]. Hence, 𝛾 𝐵𝐵 was chosen to satisfy 𝑑𝑋 𝑖 𝑑𝑡⁄ ~𝛾 𝐵𝐵 𝑅 = 0.03 m/ 5min. The concentration of CSKs was not known. Therefore, we normalized it and used the dimensionless concentration for c by setting parameters to ensure c - = 0 and c + to be almost unity. Amplitude of polarity | p | was also chosen to be almost unity. Another prerequisite for c , modeled using Ginzburg-Landau potential combined with the reaction term [38], was that the system must show bundling. Reaction rate k d should be smaller than the critical value evaluated by linear stability analysis, 𝑘 𝑑𝑐 = 𝛼 𝐽 𝛼 (1 + 𝐾 𝑀 /𝑐 )/(4𝐷 ) = 1.05 (see S1 Fig). Typically, we used 𝑘 𝑑−1 = 3.16 which corresponds to ~15 min in reality. This is comparable with the observed time scale of MT bundle dynamics in vitro [39]. In addition, we set 𝛼 𝐽 = 𝛾 𝑃 𝐾 following the model for actin filament in [12]. With these choice of parameters, the pattern formation of BBs took 600~1200 time in the simulations (Figs 2 and 5), corresponding to 2~4 days. This is consistent with actual developing time and experimental observation in cultured tracheal cells [3]. We should note that the parameter values themselves, such as equilibrium concentration of MTs ( k p and k d ), can be gradually changed during the developmental stage, which could govern the time scale of patterning ultimately. We set 𝑘 d = 1.0 × 10 −2 to obtain the results in Fig 4D, where polarity is ignored and the dynamics of CSK and BBs are simulated by using Eq. (12) and ∂ t 𝑐 = α J ∇ (−𝐷 ∇ 𝑐 + 2𝜆 𝑐𝜙 𝑐 ∑ 𝜙 𝑖𝑖 ) + 𝑘 𝑝 𝑐𝑐+𝐾 𝑀 − 𝑘 𝑑 𝑐 . (13) We set the system size as N = 30, L x = 30, and L y = 45 to simulate the model coupled with PCP information shown in Fig 5. As briefly mentioned in the main text, Neumann boundary condition was initially employed, followed by time step t = 66.7 h. The boundary condition of the oral side (top boundary in Fig 5C) was changed to Dirichlet where CSK concentration was fixed at c = c H to mimic the experimental observation that the distribution of MTs is dense in the opposite side of the cell boundary where PCP protein Vangl1 is accumulated (i.e., lung side) [3,19]. c H was chosen as 𝑐 H = 𝑐 + = 𝑐 + √𝛼 /𝛼 . Polymerization rate k p was set to be an increasing function of time as 𝑘 p 𝑘 𝑑 ⁄ = 0.6 + 0.1 tanh(1.0 × 10 −3 𝑡 − 0.5) (Fig 5B, blue line), since MTs and intermediate filaments increase during maturation of cells [8]. Model without BBs
It is useful to summarize the case where BBs are absent, to understand the behavior of the model. The model is composed of phase-separation dynamics with first order chemical reaction [40] and polarity dynamics [10]. The phase separation is driven by Landau-Ginzaburg potential 𝑓 c (𝑐) . We have numerically solved Eqs. (10 and 11) without the terms representing interaction with the BBs. As shown in S1 Fig, the pattern of CSK concentration is dependent on 𝛼 Q and the polymerization/depolymerization parameter k d [40]. For lower chemical reaction rate (i.e., smaller k d ), the CSK domain concentration is larger. CSK concentration shows stripe pattern in the middle range of k d . The range of aligned stripe is wider when the coupling constant with polarity, 𝛼 Q , is larger, while it is very narrow at 𝛼 Q = 0 . In addition, the stripe pattern is less sharp for smaller 𝛼 Q . The CSK concentration becomes uniform for larger value of k d . In summary, the model shows aligned robust bundles of CSKs for proper value of k d and higher value of 𝛼 Q . The model without BBs satisfies the minimal prerequisite to express the bundling of CSKs in our model. Acknowledgments
We thank Prof. Sachiko Tsukita and her laboratory members for helpful discussions, especially A. Tamura, T. Yano, S. Konishi, D. Taniguchi, E. Herawati, S. Nakayama, H. Kanoh, for helpful discussions and showing us their original data. We also thank Dr. D. Taniguchi for his useful advices. Reference Sengupta P. Cilia and sensory signaling: the journey from “animalcules” to human disease. PLoS Biol. 2017; 15(4): e2002240. doi: 10.1371/journal.pbio.2002240 2.
Fulford GR, Blake JR. Muco-ciliary transport in the lung. J Theor Biol. 1986; 121(4): 381-402. doi: 10.1016/S0022-5193(86)80098-4 3.
Herawati E, Taniguchi D, Kanoh H, Tateishi K, Ishihara S, Tsukita S. Multiciliated cell basal bodies align in stereotypical patterns coordinated by the apical cytoskeleton. J Cell Biol. 2016; 214(5): 571-586. doi: 10.1083/jcb.201601023 4.
Siletti K, Tarchini B, Hudspeth AJ. Daple coordinates organ-wide and cell-intrinsic polarity to pattern inner-ear hair bundles. Proc Natl Acad Sci USA. 2017; 114(52): E11170-E11179. doi: 10.1073/pnas.1716522115 5.
Pedersen H. Absence of dynein arms in endometrial cilia: cause of infertility?. Acta Obstet Gynecol Scand. 1983; 62(6): 625-627. doi: 10.3109/00016348309156261 6.
Ohata S, Nakatani J, Herranz-Pérez V, Cheng J, Belinson H, Inubushi T, et al. Loss of Dishevelleds disrupts planar polarity in ependymal motile cilia and results in hydrocephalus. Neuron. 2014; 83(3): 558-571. doi: 10.1016/j.neuron.2014.06.022 7.
Hirokawa N, Tanaka Y, Okada Y, Takeda S. Nodal flow and the generation of left-right asymmetry. Cell. 2006; 125(1): 33-45. doi: 10.1016/j.cell.2006.03.002 8.
Tateishi K, Nishida T, Inoue K, Tsukita S. Three-dimensional organization of layered apical cytoskeletal networks associated with mouse airway tissue development. Sci Rep. 2017; 7: 43783. doi: 10.1038/srep43783 9.
Kunimoto K, Yamazaki Y, Nishida T, Shinohara K, Ishikawa H, Hasegawa T, et al. Coordinated ciliary beating requires Odf2-mediated polarization of basal bodies via basal feet. Cell. 2012; 148(1-2): 189-200. doi: 10.1016/j.cell.2011.10.052 10.
Marchetti MC, Joanny JF, Ramaswamy S, Liverpool TB, Prost J, Rao M, et al. Hydrodynamics of soft active matter. Rev Mod Phys. 2013; 85(3): 1143. doi: 10.1103/RevModPhys.85.1143 11.
Prost J, Jülicher F, Joanny JF. Active gel physics. Nat Phys. 2015; 11(2): 111. doi: 10.1038/nphys3224 12.
Gowrishankar K, Ghosh S, Saha S, Rumamol C, Mayor S, Rao M. Active remodeling of cortical actin regulates spatiotemporal organization of cell surface molecules. Cell. 2012; 149(6): 1353-1367. doi:
Clare DK, Magescas J, Piolot T, Dumoux M, Vesque C, Pichard E, et al. Basal foot MTOC organizes pillar MTs required for coordination of beating cilia. Nat Commun. 2014; 5: 4888. doi: Werner ME, Hwang P, Huisman F, Taborek P, Clare CY, Mitchell BJ. Actin and microtubules drive differential aspects of planar cell polarity in multiciliated cells. J Cell Biol. 2011; 195(1): 19-26. doi:
Berens P. CircStat: a MATLAB toolbox for circular statistics. J Stat Softw. 2009; 31(10): 1-21. doi:
Chaudhuri A, Bhattacharya B, Gowrishankar K, Mayor S, Rao M. Spatiotemporal regulation of chemical reactions by active cytoskeletal remodeling. Proc Natl Acad Sci USA. 2011; 108(36): 14825-14830. doi:
Yamamoto R, Nakayama Y, Kim K. A smooth interface method for simulating liquid crystal colloid dispersions. J Phys Condens Matter. 2004; 16(19): S1945. doi:
Asakura S, Oosawa F. On interaction between two bodies immersed in a solution of macromolecules. J Chem Phys. 1954; 22(7): 1255-1256. doi:
Vladar EK, Bayly RD, Sangoram AM, Scott MP, Axelrod JD. Microtubules enable the planar cell polarity of airway cilia. Curr Biol. 2012; 22(23): 2203-2212. Doi:
Devenport D. The cell biology of planar cell polarity. J cell Biol. 2014; 207(2): 171-179. doi:
Shimada Y, Yonemura S, Ohkura H, Strutt D, Uemura T. Polarized transport of Frizzled along the planar microtubule arrays in Drosophila wing epithelium. Dev Cell. 2006; 10(2): 209-222. doi:
Matis M, Russler-Germain DA, Hu Q, Tomlin CJ, Axelrod JD. Microtubules provide directional information for core PCP function. Elife. 2014; 3: e02893. doi: 10.7554/eLife.02893 23.
Butler MT, Wallingford JB. Planar cell polarity in development and disease. Nat Rev Mol Cell Biol. 2017; 18(6): 375. doi:
Werner ME, Mitchell BJ. Planar cell polarity: microtubules make the connection with cilia. Curr Biol. 2012; 22(23): R1001-R1004. doi:
Takagishi M, Sawada M, Ohata S, Asai N, Enomoto A, Takahashi K, et al. Daple coordinates planar polarized microtubule dynamics in ependymal cells and contributes to hydrocephalus. Cell Rep. 2017; 20(4): 960-972. doi:
Surrey T, Nédélec F, Leibler S, Karsenti E. Physical properties determining self-organization of motors and microtubules. Science. 2001; 292(5519): 1167-1171. doi:
Nogucci H, Ishihara S. Collective dynamics of active filament complexes. Phys Rev E. 2016; 93(5): 052406. doi: Vignaud T, Blanchoin L, Théry M. Directed cytoskeleton self-organization. Trends Cell Biol. 2012; 22(12): 671-682. doi:
Mitchell B, Jacobs R, Li J, Chien S, Kintner C. A positive feedback mechanism governs the polarity and motion of motile cilia. Nature. 2007; 447(7140): 97. doi:
Guirao B, Meunier A, Mortaud S, Aguilar A, Corsi JM, Strehl L, et al. Coupling between hydrodynamic forces and planar cell polarity orients mammalian motile cilia.
Nat Cell Biol. 2010; 12(4): 341. doi:
Namba T, Nishikawa M, Shibata T. The relation of signal transduction to the sensitivity and dynamic range of bacterial chemotaxis. Biophys J. 2012; 103(6): 1390-1399. doi:
Tee YH, Shemesh T, Thiagarajan V, Hariadi RF, Anderson KL, Page C, et al. Cellular chirality arising from the self-organization of the actin cytoskeleton. Nat Cell Biol. 2015; 17(4): 445. doi:
Inaki M, Liu J, Matsuno K. Cell chirality: its origin and roles in left–right asymmetric development. Philos Trans R Soc Lond B Biol Sci. 2016; 371(1710): 20150403. doi:
Ray P, Chin AS, Worley KE, Fan J, Kaur G, Wu M, et al. Intrinsic cellular chirality regulates left–right symmetry breaking during cardiac looping. Proc Natl Acad Sci USA. 2018; 115(50): E11568-E11577. doi:
Chin AS, Worley KE, Ray P, Kaur G, Fan J, Wan LQ. Epithelial cell chirality revealed by three-dimensional spontaneous rotation. Proc Natl Acad Sci USA. 2018; 115(48): 12188-12193. doi:
Nonaka S, Tanaka Y, Okada Y, Takeda S, Harada A, Kanai Y, et al. Randomization of left–right asymmetry due to loss of nodal cilia generating leftward flow of extraembryonic fluid in mice lacking KIF3B motor protein. Cell. 1998; 95(6): 829-837. doi:
Yoshiba S, Shiratori H, Kuo IY, Kawasumi A, Shinohara K, Nonaka S, et al. Cilia at the node of mouse embryos sense fluid flow for left-right determination via Pkd2. Science. 2012; 338(6104): 226-231. doi:
Schmid F. Physical mechanisms of micro-and nanodomain formation in multicomponent lipid membranes. Biochim Biophys Acta Biomembr. 2017; 1859(4): 509-528. doi:
Torisawa T, Taniguchi D, Ishihara S, Oiwa K. Spontaneous formation of a globally connected contractile network in a microtubule-motor system. Biophys J. 2016; 111(2): 373-385. doi: Huberman BA. Striations in chemical reactions. J Chem Phys. 1976; 65(5): 2013-2019. doi:
Turk E, Wills AA, Kwon T, Sedzinski J, Wallingford JB, Stearns T. Zeta-tubulin is a member of a conserved tubulin module and is a component of the centriolar basal foot in multiciliated cells. Current Biology. 2015; 25(16): 2177-2183. doi: 10.1016/j.cub.2015.06.063 42.
Boutin C, Labedan P, Dimidschstein J, Richard F, Cremer H, André P, et al. A dual role for planar cell polarity genes in ciliated cells. Proc Natl Acad Sci USA. 2014; 111(30): E3129-E3138. doi: 10.1073/pnas.1404988111 Supporting information
S1 Fig. Patterns of CSKs without BBs.
The patterns are dependent on α Q and 𝑘 d ; Striped pattern appears for finite value of α Q (striped patterns which are not sharp may appear at α Q = 0 in a narrow parameter region). CSKs concentration becomes uniform at very high values of polymerization/depolymerization rate k d . Bars, 0.5 µm. 9 S2 Fig. BB positional order depends on the elasticity parameter, K . Color indicates the straightness of BB quantified by 〈−cos𝜓 𝑖 〉 . Inset figures show typical pattern of BBs. Gray dashed line ( 𝐾 = 3.0 ) corresponds to the line connecting between symbol A and PA in Fig 3. Bars, 0.5 µm.
S3 Fig. Dependence of BB pattern on the form of condensation term 𝒇 𝐜 (𝒄) . The form of 𝑓 c is expressed as 𝑓 c (𝑐) = 𝜆 W (𝑐 − 𝑐 − 𝑊) (𝑐 − 𝑐 + 𝑊) . Left: phase diagram against 𝜆 W and 𝑊 . Color indicates the straightness of BB ( 〈−cos𝜓 i 〉 ). Right: Four representative patterns are shown (i-iv) with corresponding black squares in the left panel. Bars, 0.5 µm. S1 Video.
Model simulation of BBs corresponding to Fig 2B. BBs reach “alignment” pattern. Bars, 0.5 µm.
S2 Video.
Model simulation of BBs corresponding to Fig 2C. BBs reach “partial alignment” pattern. Bars, 0.5 µm.
S3 Video.
Model simulation of BBs without the condensation term, corresponding to Fig 4A. Bars, 0.5 µm.
S4 Video.
Model simulation for nocodazole treated cells, corresponding to Fig 4D. Bars, 0.5 µm.