# Spatiotemporal organization of cilia drives multiscale mucus swirls in model human bronchial epithelium

Mustapha-Kamel Khelloufi, Etienne Loiseau, Marc Jaeger, Nicolas Molinari, Pascal Chanez, Delphine Gras, Annie Viallat

<figure><img src="image_1.png" /></figure>

To cite this version:

Mustapha-Kamel Khelloufi, Etienne Loiseau, Marc Jaeger, Nicolas Molinari, Pascal Chanez, et al. Spatiotemporal organization of cilia drives multiscale mucus swirls in model human bronchial epithelium. Scientific Reports, 2018, 8, pp.2447. 10.1038/s41598-018-20882-4. hal-01821276

HAL Id: hal-01821276

https://hal.science/hal-01821276

Submitted on 3 Jul 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

<figure><img src="image_3.png" /></figure>

OPEN

# Spatiotemporal organization of cilia drives multiscale mucus swirls in model human bronchial epithelium

Received: 15 December 2017

Accepted: 25 January 2018

Published online: 05 February 2018

Mustapha-Kamel Khelloufi<fnref n="1" />, Etienne Loiseau<fnref n="1" />, Marc Jaeger², Nicolas Molinari³,⁴, Pascal Chanez⁵, Delphine Gras⁵ & Annie Viallat<fnref n="1" />

Mucociliary clearance is a biomechanical mechanism of airway protection. It consists of the active transport along the bronchial tree of the mucus, a fluid propelled by the coordinated beating of a myriad of cilia on the epithelial surface of the respiratory tract. The physics of mucus transport is poorly understood because it involves complex phenomena such as long-range hydrodynamic interactions, active collective ciliary motion, and the complex rheology of mucus. We propose a quantitative physical analysis of the ciliary activity and mucus transport on a large panel of human bronchial cultures from control subjects, patients with asthma and chronic obstructive pulmonary disease obtained from endobronchial biopsies. Here we report on the existence of multiple ciliary domains with sizes ranging from the tens of a micron to the centimeter, where ciliary beats present a circular orientational order. These domains are associated with circular mucus flow patterns, whose size scales with the average cilia density. In these domains, we find that the radial increase of the ciliated cell density coupled with the increase in the orientational order of ciliary beats result in a net local force proportional to the mucus velocity. We propose a phenomenological physical model that supports our results.

The continuous increase of chronic respiratory diseases worldwide, such as severe asthma, chronic obstructive pulmonary disease (COPD) and cystic fibrosis, represent a major threat for the quality of life. The pathophysiology associated to severe asthma and COPD results from an abnormal bronchial epithelium<fnref n="1" /> leading to a sharp increase in susceptibility to infection and inflammation². The bronchial epithelium is the first barrier to protect the respiratory tract via an innate mechanism called mucociliary clearance. It consists in the active transport of a sticky fluid, the mucus, via a myriad of cilia along the epithelial surface of the airways. The mucus traps inhaled pathogens and the protective role of the mucociliary clearance relies on the ability of cilia to self-organize and coordinate their beating in order to transport the mucus along the bronchial tree till its elimination by swallowing or expectoration. From a biophysical point of view, the bronchial epithelium is an active system which involves two main key players. First, ATP-driven microscopic cilia. Located at the apical surface of epithelial ciliated cells and immersed in a swollen grafted polyanionic gel called periciliary layer (PCL)³, they beat asymmetrically at their own frequency, direction and phase with their tip poking in the mucus during the forward stroke (see Fig. 1A)⁴,⁵. Second, the mucus layer, whose complex shear-thinning properties govern its flow on top of the ciliary layer. The mucus behaves like a cohesive gel in absence of external mechanical stress and like a viscous fluid under stress⁶. Despite a large corpus of clinical and biological studies, the physical mechanisms of long-range transport of mucus remain poorly understood due to the intrinsic complexity of the mucociliary system. Indeed, the force propelling the mucus is exerted by the displacement of cilia tips in the mucus layer and depends on the patterns of individual ciliary beats, on the coordination of their phase and direction⁷ and on the nonlinear viscoelastic response of the mucus. The transport of mucus also depends on the boundary conditions, specifically on the friction between the fluid PCL and the top mucus layer, resulting in a strongly coupled system. Furthermore, the system is responsive. Ciliary beat frequencies (CBF) were found to change with the viscosity of mucus⁸, suggesting

<footnote n="1">¹Aix Marseille Univ, CNRS, CINAM, Marseille, France. ²Aix Marseille Univ, CNRS, Centrale Marseille, M2P2, Marseille, France. ³Institut Montpelliérain Alexander Grothendieck, CNRS, Univ, Montpellier, France. ⁴Department of Statistics, University of Montpellier Hospitals, Montpellier, France. ⁵Aix Marseille Univ, Inserm, INRA, C2VN, Marseille, France. Correspondence and requests for materials should be addressed to A.V. (email: viallat@cinam.univ-mrs.fr)</footnote>
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Air Liquid Interface (ALI) bronchial epithelium. (A) Schematic of the ALI chamber: the basal side of the differentiated epithelium is in contact with the culture medium through a porous membrane and the apical side points towards the air. The epithelium has a pseudo stratified structure made of basal cells, goblet cells which produce the mucus and ciliated cells. During the forward stroke cilia tips poke in the mucus layer while their recovery motion takes place in the periciliary layer. (B) Average active ciliated cell density quantified for control cultures (n = 17), mild asthma (MA, n = 7), severe asthma (SA, n = 18) and chronic obstructive pulmonary disease (COPD, n = 7). Results are expressed as median with 10 and 90 percentiles.</figcaption></figure>

that cilia might act as motile mechanosensitive organelles able to adapt their CBF with the drag force they experience in the mucus layer. Cilia are also able to transiently change their beat directions, as observed in the brain ventricles of mice⁹. *In-vivo* and *ex-vivo* explorations of the spatiotemporal coordination of ciliary beats and its coupling to mucociliary clearance has been little experimentally studied because suitable methods like inhalation of radiolabeled markers¹⁰ to assess mucus velocity or methods directly using *ex-vivo* preparation of trachea are invasive. Though providing valuable information, mucociliary clearance rates measured in animal models are not fully relevant for humans, especially in mice models¹¹⁻¹⁵. The development of well-differentiated human airway cultures at the air-liquid interface (ALI) allows overcoming previous limitations¹⁶. This culture model retains the mixed mucus-secretory and ciliated surface phenotype expressed *in vivo*¹⁷. Seminal works on cystic fibrosis¹⁸⁻²¹ reported that mucus could be spontaneously transported as a rigid disc in a circular pattern over the whole culture chamber for more than one week. The associated ciliary activity was however not investigated. It has been recently investigated in the ventricular system of mice brain, where highly organized directional patterns of cilia modules were discovered, leading to spatiotemporal modulated complex flow patterns of the cerebrospinal fluid, including swirls⁹. The spatiotemporal organization of ciliary activity and its interplay with mucus transport in the mucociliary clearance are still to be understood in health and disease, in particular, dyskinesia and abnormal cilia beat patterns have already been reported on bronchial brushings in severe asthma²². Here, we propose a quantitative physical analysis of ciliary activity and mucus transport on a large panel of human bronchial ALI cultures from control subjects, patients with mild asthma, severe asthma and COPD obtained from endobronchial biopsies. We report on the existence of multiple ciliary domains, whose sizes ranges from tens of microns to centimeters, where ciliary beats present a circular orientational order. These domains are associated with a circular mucus flow patterns, whose size scales with the average cilia density. We quantify the high degree of organization of ciliary activity on these domains and propose a phenomenological model which links mucus velocity to ciliary organization. Finally, we discuss our findings in the context of severe asthma and COPD.

## Results

We first controlled, using confocal microscopy, that ALI cultures exhibit the pseudo-stratified structure (see Fig. S1A) typical of bronchial epithelium, with basal cells, goblet cells, and ciliated cells (see the schematic Fig. 1A). The observed average ciliary length of 6.8 µm is in agreement with values reported in the literature²³. For
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Cilia orientation map in a swirl. (A) The image is a standard deviation projection over 500 images acquired at 50 fps of a the central part of a domain where ciliary beats have a circular orientational order (mild asthma culture). The spots are the tufts of cilia and the trails are due to dead cells/fragments embedded in the mucus. The color codes for the orientation of both mucus flow and ciliary beats. The white lines demarcate the main areas with a preferential orientation indicated by the arrows. On average there is a correlation between ciliary beat orientations and mucus flow. Scale bar, 50 µm. On the left, orientation scale. (B and C) Zoom on the cilia. Bright field image in (B), the white arrow indicates the direction of mucus transport. Corresponding standard deviation projection in (C). Scale bars, 20 µm.</figcaption></figure>

patients suffering of severe asthma, epithelial defects such as loosely connected adjacent cells were often observed (Fig. S1B). We quantified the ciliary activity of ALI cultures from 17 control subjects, 7 patients with mild asthma, 18 patients with severe asthma and 7 patients with COPD.

**Average active ciliated cell density.** The average active ciliated-cell density, $v$, was quantified by measuring the fraction of the epithelial surface covered by beating cilia. We found that $v$ ranges from 0 to 0.70 with a large variability among different donors and an average value of 0.16 for controls. No significant change of $v$ is observed after elimination of the mucus layer at the surface of the epithelium. These values are in agreement with those previously reported on ALI cultures²⁴,²⁵, but lower than those observed in the trachea, wherein the percentage of ciliated cell is about 50%²⁶. $v$ is not altered in mild asthma ($v$=0.142, non significant) but strongly decreases to 0.055 ($p$ = 0.0032) and 0.0395 ($p$ = 0.0168) in severe asthma and COPD respectively (Fig. 1B). This strong decrease supports previous studies performed on epithelial cells collected from endobronchial biopsies and observed in solution²². It also supports that severe asthma and COPD are epithelial diseases¹, associated with morphological abnormalities of the bronchial epithelium.

**Frequency and orientational order of Ciliary beats.** The average ciliary beat frequency (CBF) measured on control subject is 18.3 Hz with a dispersion of ±3.2. A typical frequency map is shown on Fig. S2 for the ciliary domain of Fig. 2B and C and Movie S1. More details on CBF in health and diseases are given in SI 1 and Fig. S2. Steric interactions between neighboring cilia impose a directional order at the scale of one ciliated cell or adjacent ciliated cells. A typical example for an isolated group of neighboring cells is shown in Movie S2, where two main beat directions are observed. Surprisingly, at a larger scale and up to the scale of the whole culture chamber, we almost systematically observed the existence of spatial domains on which ciliary beats present an striking circular orientational order. The radius of these domains ranges between 10 microns (a few cells) to the whole culture chamber (R = 6 mm). Several swirls can be found within a culture chamber, randomly distributed
in the central or peripheral zones of the cell culture. (see movie S3). A typical example of an orientation map of ciliary beats is shown in Fig. 2A on an epithelium with $v$ equal to 0.2, where the circular domain spans the whole culture chamber ($R = 6$ mm). The tangential orientational order of the ciliary beats, illustrated on a disc in the central zone is visible in Fig. 2A. The orientation map of ciliary beats at a peripheral zone (Fig. 2B and Movie S1) of the circular domain is shown in Fig. 2C. Although ciliated cells may be distant from one another and cilia do not form a dense mat, an orientational order of ciliary beats is clearly visible. This suggests that additional mechanisms to steric interactions are involved in the coordination of beat direction. It should be noted that not all cilia are beating in the same direction (see cilia in purple, for example). We will now show that these circular orientational patterns of ciliary beats generate multiscale flows of the fluid located at the epithelial surface.

**Rotational surface fluid flows.** To study the ability of cilia to generate local flows, we performed an apical washing to remove the native mucus and replaced it by a dilute solution of 1-µm beads in PBS. Typical bead trajectories are shown in Fig. 3A and B and movie S4. In areas lacking of ciliated cells the beads remain immobile, while elsewhere many beads present a net motion over hundreds of micrometers. Two sort of trajectories are illustrated in Fig. 3B: a linear or pseudo-linear bead motion associated with a linear spatial repartition of cilia, and, strikingly, a circular motion on domains of radii ranging from 10 µm (a few ciliated cells) to hundreds of microns depending on the culture chamber. Such a transport over several hundreds of microns reveals that on these circular trajectories the local force that applies on every bead is purely tangential of the whole cellular culture. Several clockwise and counterclockwise swirls can coexist in a single chamber, together with linear flows indicating that these circular domains are independent from one another and from the geometry of the chamber. These swirls are observed at the surface of epithelial domains on which ciliary beats present a circular orientational order, thus showing that these swirls are directly generated by the self-organization of ciliary beats. Indeed, the local forces controlling the micron-size displacement of the latex beads must be remarkably organized at the macroscopic scale of the bead trajectories to ensure their circular displacement. This point will be detailed thereafter.

It should be noted that before performing an apical washing, the samples where latex beads rotate on radii of several hundreds of microns, exhibited a macroscopic rotational motion of native mucus spanning the whole chamber (see Fig. 3C and movie S3). Such a motion has been reported by Matsui et al.¹⁸. Similarly to latex bead transport, mucus was observed to rotate clockwise or counterclockwise with an equal probability, thus revealing a phenomenon of stochastic origin In many cases the angular velocity of mucus is constant (see Fig. S3), thus showing that it behaves like a rigid body. Interestingly, the mucus velocity at the edge of the chamber is of the order of 200 µm.s⁻¹, similar to the range of physiological bronchial mucus velocities¹⁴,²⁷,²⁸.

**The average density of cilia sets the size of the swirls.** The size of the largest circular domains of bead/mucus transport in each culture chamber scales, via a power law, with $v$ over three orders of magnitude (see Fig. 4). For circular domains spanning the whole chamber (see Fig. 3C and Movie S5), $v$ is in the range 20–30%, indicating that a dense ciliated epithelium is not necessary to drive the mucus flow at the centimeter scale. In addition, physiological data reported in the literature on mucus transport in the trachea fall on the same curve. This master curve points out that the average ciliated cell density directly governs the size of epithelial surface on which a directional flow transport is observed, or in other words, the distance over which a surface fluid (beads in a Newtonian fluid or mucus) can be transported along the epithelium. This curve may be one key to understand mucus congestion. In particular, this curve shows that low $v$ values observed for patients with severe asthma and COPD displayed in Fig. 1B might be the main cause of altered mucociliary clearance.

**Experimental analysis of the local forces exerted on swirls.** The circular transport of mucus at a constant angular velocity on a disk requires a fine regulation over the entire disk, in terms of direction and magnitude, of the forces that locally propel the mucus. Indeed, along a diameter, the local forces should have the same tangential direction and an opposite direction on either side of the center. The magnitude of local forces must also vary along the radius to drive a transport at a constant angular velocity, that is with a curvilinear velocity which increases linearly with the distance from the center of the swirl. In this section we quantify how such a fine regulation is provided by the ciliated epithelium. We analyse, within a swirl, the spatial and dynamic self-organization of the cilia underneath the mucus for 6 domains from 5 different culture chambers. We performed the experiments on commercial cultures from healthy human samples (MucilAir from Epithelix, Switzerland). Most commercial Epithelix cultures are characterized by a high $v$ (more than 65%) and thus, exhibit mucus rotation over the whole culture chamber (the radius of the chamber is 3.25 mm). We studied two such cultures. Three other cultures were provided by Epithelix at a younger stage, with epithelial cells still being differentiated. On reception, the number of ciliated cells in the culture was small and $v$ was therefore low. After elimination of mucus and 1-µm beads addition, no significant bead motion was observed. After a few days, upon increase of $v$, large rotational flow domains were observed at the epithelial surface, where patches of mucus rotated at constant angular velocity during more than 24 hours (Fig. S3 and Movie S3). The existence of local rotational flows and epithelial domains with an orientational order of ciliary beats, both on cultures prepared in our laboratory and on commercial ones, clearly show the general and robust nature of these phenomena. To decipher the physical parameters that determine the mucus propelling force on a rotational flow domain, we consider that each active cilium exerts a force on the surface fluid when the cilium tip penetrates the fluid. This force has the direction of the ciliary beat and its useful component is the tangential one, which circularly propels the surface fluid (see Fig. 5A). At first approximation we assume that the force exerted on the surface fluid by an area element of the epithelium is proportional to $v$ on this area element. By decomposing the rotational domain into concentric annuli, the force per unit area on an annulus located at a distance $r$ from the center of the domain is therefore proportional to $v_r < \sin \alpha >_r$, with $\alpha$ the angle between a beat direction and the radial direction of the circular domain, $v_r$ and $< \sin \alpha >_r$ are $v$ and the average of $\sin \alpha$ on the considered annulus.
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Cilia driven transport. (A) Transport of 1-µm fluorescent beads at the apical surface of a washed epithelium. The image is a standard deviation projection over 50 images. Scale bar, 200 µm. (B) Local transport. Trajectories are obtained by tracking 1-µm beads. Scale bar, 20 µm. (C) Rotational transport of mucus spanning the whole chamber (same sample than in A). The image is a standard deviation projection over 50 images. Scale bar, 200 µm.</figcaption></figure>

We quantified both $v_r$ and $\langle sin\alpha \rangle_r$ for different values of r on each rotational domain. Their respective radial evolutions are plotted in Fig. 5B and C. Both the local ciliated density and the average tangential component of ciliary beats increase with the distance to the center of the domain. Moreover, the product $v_r \langle sin\alpha \rangle_r$, which is proportional to the mucus driving force, increases linearly with r on each localized domains (A to C) and on the central part ($r < 200$ µm) of the two entire chamber-spanning domains D and E (see Fig. 5D). This remarkable radial increase of the driving force ensures a stable circular movement of the mucus. These results show that there is a fine spatial organization of the epithelium in terms of ciliary beat direction and cilia density. The adequate local cilia density could be regulated during the initial process of differentiation or by direct transdifferentiation without proliferation among epithelial cells²⁹. This transdifferentiation might be triggered by mechanotransduction due to the resistive mechanical force exerted by the periciliary layer and by the viscoelastic mucus on beating cilia. Another way to spatially regulate cilia density is the re-arrangement of the ciliated cells within the
<figure><img src="image_7.png" /><figcaption>Figure 4. The size of swirls scales with $\nu$. The radius of the swirls is plotted versus the average $\nu$ for control (red circles), mild asthma (black squares) and severe asthma (blue diamonds) ALI cultures. The data for the trachea (green triangles) are from the literature$^{26,35}$. The size of the swirls scales with $\nu$ over 3 orders of magnitude, from R ~ 10 $\mu$m up to the radius of the ALI culture (R = 6 mm). Interestingly, the data for the trachea fall on the same curve.</figcaption></figure>

epithelium. Indeed we observed over hours the displacement of the cells with an average drift of 22 $\mu$m/h. This re-arrangement could lead in changes in cell orientations and planar cell polarity, which might drive changes in beating direction$^7$. In addition, the radial increase in orientational order might result from steric constraints between neighboring cilia due to the radial increase of $\nu$.

Interestingly, the force per unit area does not vanish at the center of a domain. Its value is expected to be directly linked to the resistive friction exerted by the periciliary layer on the surface mucus layer. Finally, as both the local force per unit area and the tangential velocity of the surface fluid vary linearly with r, they also vary linearly with each other as shown in Fig. 5E. Consequently, the mucus is transported at constant angular velocity. It is therefore not radially sheared in these rotational domains, thus limiting viscous dissipation. This is energetically favorable but requires a high degree of self-organization of spatial cilia density and temporal ciliary beats. The physical interpretation of this linear relationship is not straightforward. In order to better understand the origin of this minimal force that has to be applied to generate a surface-fluid flow and the linear velocity-force relationship, we developed a semi-quantitative phenomenological model that we now detail.

**Phenomenological model.** The linear relationship between velocity and force results from two opposite actions on the surface mucus layer, the propulsion exerted by cilia and the resistive friction exerted by the periciliary layer. Here we propose a model that describes the mucociliary system as a periciliary layer with cilia tufts separated by empty zones, where cilia beat near a no-slip wall. Above, the mucus is considered as a yield-stress fluid (Fig. 6). We hypothesize that the yield-stress is reached only when cilia tips penetrate into mucus. Everywhere else mucus behaves like a rigid body and has a constant velocity $V_{SF}$. Within a circular domain, the propulsive force per unit area exerted in an annulus at a distance r from the center writes as:

$$ \frac{\varphi n_c f_c}{A_{cc}} \nu_r < \sin \alpha >_r, \qquad (1) $$

where $\varphi$ is the time fraction during which the propulsion occurs, $n_c$ is the average number of cilia per ciliated cell, $A_{cc}$ is the area of a ciliated cell and $f_c$ is the force exerted by one cilium to the surface mucus layer. $f_c$ results from the motion of the cilium tip in mucus and is expected to be proportional to the CBF and the fluid viscosity (see SI 2). The resistive friction force arises during the recovery of the cilia to their initial position. They drag the periciliary fluid backward in the ciliated zones and in the neighboring non-ciliated zones, whereas the top of the periciliary layer is dragged forward by mucus. This results in a shear layer in the upper part of the periciliary layer as shown in Fig. 6. The stress in this layer in an annulus at distance r from the center of the circular domain is $\nu_r(1 - \varphi)\eta_{pc}^{\frac{(V_{SF}(r) + V_R)}{e}}$, in ciliated zones and $(1 - \nu_r)\eta_{pc}^{\frac{(V_{SF}(r) + V_R)}{e}}$ in neighboring non-ciliated zones, where e is the thickness of the sheared layer, $\eta_{pc}$ is the viscosity of the periciliary fluid, and $V_R$ is the velocity of the recovering cilia tips. The shear in the lower part of the periciliary layer is lower (small variation in speed on a large fluid thickness). The total stress is $\tau = (1 - \varphi\nu_r)\eta_{pc}^{\frac{(V_{SF}(r) + V_R)}{e}} \approx \eta_{pc}^{\frac{(V_{SF}(r) + V_R)}{e}}$ since $\varphi\nu_r \ll 1$. The balance between propulsive and resistive force yields an affine relation of the force with the surface fluid velocity, in agreement with experimental results:

$$ \frac{\varphi n_c f_c}{A_{cc}} \nu_r < \sin \alpha >_r = \eta_{pc} \frac{(V_{SF}(r) + V_R)}{e}. \qquad (2) $$
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Driving force on 6 swirls. (A) The size and the center of a swirl are determined on an image, which results from the projection of the maximum intensity of the acquired movie. 1-µm fluorescent beads are added on the mucus to facilitate the visualization. The standard deviation projection of the movie results in an image on which we can quantify *v* and the ciliary beat directions on annuli. *α* is the angle between the ciliary beat direction of a ciliary tuft and the radius passing through the tuft. (B) The evolution of *v* measured on an annulus at distance *r* from the center of the swirl, *v*<sub>r</sub>, is plotted versus *r* for 6 different swirls. (C) The evolution of Sin *α* averaged on an annulus at distance *r* from the center of the swirl, <Sin *α*>, is plotted versus *r* for the swirls analyzed in (B) (D) The propulsive force, proportional to the product of *v*<sub>r</sub> with <Sin *α*> increases radially in an affine way. (E) The curvilinear mucus velocity increases linearly with the propulsive force.</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Schematic of the mucociliary coupled system. Geometry of the mucociliary system considered for the phenomenological model. $V_R$ is the recovery ciliary velocity, $V_{SF}$ is the surface fluid velocity, $e$ is the thickness of the sheared layer.</figcaption></figure>

According to Eq. 2, the constant friction term observed at $r=0$ is directly linked to the backward cilia velocity. Therefore, the slope variations observed in Fig. 5E may be due to the viscosity of the periciliary layer, the number of cilia per ciliated cell or the ciliary beat pattern motion. The slope of the variations of $V_{SF}$ versus $v_r \sin\alpha > r$, $\frac{e}{\eta_{pc}} \frac{\varphi^n f_c}{A_{cc}}$, can be roughly estimated to $4 \times 10^{-5} \text{ m.s}^{-1} (40 \mu\text{m.s}^{-1})$, by using typical numerical values: $e \approx 10^{-6} \text{ m}$, $\eta_{pc} \approx 510^{-3} \text{ Pa.s}$, $\varphi = 0.1$, $n_c \approx 160$, $A_{cc} \approx 810^{-11} \text{ m}^2$ where $f_c$ has been estimated to $1 \text{ pN}$, as detailed in SI 2. This value is in the range on the ones experimentally observed in Fig. 5E. This simple phenomenological model yields numerical predictions in agreement with experiments and it captures the linear increase of the propulsive force with the distance to the domain center and the linear relationship between local force and mucus velocity. It shows that a 'rigid-like' motion at constant angular velocity is generated on the epithelial surface, enabling mucus to be transported without undergoing shearing forces.

## Discussion

To summarize the findings, we reported, on non fully differentiated and mature cultures, on the existence of epithelial domains where ciliary beats present a stable circular orientational order. These domains reveals circular flow patterns of mucus and of micron-size beads placed at the apical surface of the epithelium after mucus removal. This generic behavior was observed in human ALI cultures both commercial and made in-house from endobronchial biopsies on controls and asthmatic patients. The size of these domains, which ranges from tens of microns to centimeters, is controlled by the average ciliated cell density via a scaling law over more than three orders of magnitude. This physical quantity governs the distance over which the mucus is efficiently transported. Noteworthy, based on this curve, we found that the epithelium from severe asthma and COPD patients does not seem to be able to generate a sufficient number of active ciliated cells to ensure a long range transport of mucus. Our detailed quantitative analysis of 6 circular flow domains showed a radial increase in ciliated cell density and in orientation order of ciliary beats on these domains as well as a linear radial increase in mucus velocity. An important and notable finding is the radial regulation of both the ciliated cell density and the ciliary beat orientation which results in the linear radial increase of the net propulsive force. Finally, we developed a phenomenological model, which captures the observed linear relationship between the mucus velocity and the propulsive force, and predicts numerical values in agreement with the experiments. The $v$ and CBF values we found are respectively significantly lower and higher than those measured *in-vivo* in the trachea. On the basis of our findings, $v$ value of about 0.2 is enough to ensure a mucus transport over the whole culture chamber of centimeter size. A higher $v$ is therefore not required for efficiently transport the mucus on these ALI cultures. An increase in CBF could be a mechanism to provide the power necessary to handle a suitable mucus velocity. Indeed, the power provided by the ciliary beats is proportional to the total number of beats per unit time and thus, to both $v$ and CBF. A high frequency could therefore compensate a low ACDD to get a mucus velocity in the culture chamber in the range observed in the trachea. Such a phenomenon of high CBF compensating the lack of ciliary coordination is indeed physiologically relevant. It has been evoked for the clearance of amniotic fluid from newborn lungs³⁰. Our study points out the fundamental role of the long-range spatial regulation of ciliated cell density and coordination of ciliary beats direction. Yet, the mechanisms underlying this regulation/coordination remain unclear. They are however crucial in chronic respiratory disease and also, in the development of the ciliated lung epithelium and in the regeneration of ciliary beat-induced flow after injury. Long-range fluid-structure hydrodynamic interactions combined with steric interactions are known to generate meso-scale coherent flow patterns, among which circular patterns, in active suspensions³¹,³². We suggest that hydrodynamic interactions in both periciliary and mucus layer may give rise to the emergence of circular flow patterns, driven by the circularly oriented drag force exerted by both the mucus and the periciliary fluid on ciliated cells. Ciliated cells could actively respond to this drag force, via a mechanotransduction mechanism, leading to spatial adaptation of ciliated cells density and of the orientation of ciliary beats. Although the mucus flow has been characterized in the present study, the flows in the periciliary layer, which may be paramount for mucus transport and ciliary activity are still poorly documented. Experimental and computational efforts on epithelial surface flows and collective motion in cilia assemblies are
today crucial to shed light on the mechanisms governing the spatiotemporal organization of ciliary activity and to take the understanding of mucociliary clearance to the next level.

**Methods**

**Endobronchial biopsy specimens and ALI cultures.** Subjects were recruited for endobronchial biopsies at the Clinique des Bronches, de l'Allergie et du Sommeil (Asistance Publique HÀ pitaux de Marseille, Marseille, France). The study protocol was approved by the local ethic committee and was part of the Bronchial Obstruction and Asthma Cohort (COBRA - Cohorte Obstruction Bronchique et Asthme) sponsored by the French National Institute of Health and Medical Research, INSERM [IDRGB 2008-A00284-51, Afssaps 2008-0113]. All subjects were informed about the nature and purpose of the study and provided written consent before enrolment. All experiments were performed in accordance with relevant guidelines and regulations.

*In-vitro* reconstituted bronchial epithelium were developed according to a previously existing protocol<sup>17,33</sup> and briefly described in SI3. Alternatively, commercial ALI cultures (MucilAir) were bought from Epithelix.

**Cilia activity.** Transmission light microscopy was performed on a Nikon Eclipse Ti microscope at 37°C and 5% CO2 under a humid atmosphere. For each subject, three separate culture chambers were examined and on each culture chamber, the central zone and two opposite side zones, each one of area 255 µm² were observed. In each zone, the ciliary motion was recorded at 250 fps for 3 seconds with a camera Photron SA4. Then, using Fiji software, the average intensity of each pixel obtained from all images of the stack of 250 × 3 images was determined and subtracted, pixel per pixel to each image of the stack, in order to enhance the visualization of mobile cilia. In each zone, the number of ciliated cells with motile cilia was counted, the fraction of the surface area covered by beating cilia was measured, and the ciliary beating frequencies were determined over 20 cycles on 5 ciliated cells randomly chosen. Cilia located on the same ciliated cell were found to beat at the same frequency. *v* and average ciliary beating frequency were determined for each patient after average of the measured quantities in the three zones and the three culture chambers.

**Measurement of transport properties.** We imaged the mucus layer with a Photron SA4 camera or Andor Neo sCMOS camera (×10, ×20 and ×60 objectives) and computed the mucus velocity by tracking the motion of mucus heterogeneities or latex beads (microspheres, 1 µm, red fluorescent) embedded in the mucus layer. The velocity and trajectory of latex beads on the epithelial surface after mucus removal were measured as follows. Mucus was removed by performing an apical washing with PBS supplemented with 10 mM DTT during 5 minutes, followed by 5 rinsings with PBS. A dilute suspension of latex beads in PBS was then placed on the epithelial surface in order to form a very thin layer (few microns) over the whole culture chamber surface. The motion of latex beads was then recorded with a Photron SA4 camera or Andor Neo sCMOS camera (×10, ×20 and ×60 objectives)

**Statistical analysis.** In order to statistically analyze *v* and CBF values, linear mixed-effects models were used. Dependent variables were *v* and CBF, independent variable was the patient's group (controls, mild asthma, severe asthma and COPD patients). A p-value < 0.05 is considered as significant. Statistical analyses were performed using R (version 3.1.3).

**Image processing.** We developed a homemade routine in python to compute the CBF based on Fourier analysis. The orientation maps were computed by implementing in python the algorithm described in Püspöki et al.<sup>34</sup> based on the structure tensor computation.

**References**

1. Chanez, P. Severe asthma is an epithelial disease. *European Respiratory Journal* **25**, 945–946 (2005).

2. Camner, P., Mossberg, B. & Philipson, K. Tracheobronchial clearance and chronic obstructive lung disease. *Scandinavian journal of respiratory diseases* **54**(5), 272–281 (1973).

3. Button, B. et al. A periciliary brush promotes the lung health by separating the mucus layer from airway epithelia. *Science* **337**, 937–941, https://doi.org/10.1126/science.1223012 (2012).

4. Sanderson, M. & Sleigh, M. Ciliary activity of cultured rabbit tracheal epithelium: beat pattern and metachrony. *Journal of Cell Science* **47**, 331–347 (1981).

5. Knowles, M. R. & Boucher, R. C. Mucus clearance as a primary innate defense mechanism for mammalian airways. *The Journal of clinical investigation* **109**, 571 (2002).

6. Lai, S. K., Wang, Y.-Y., Wirtz, D. & Hanes, J. Micro-and macrorheology of mucus. *Advanced drug delivery reviews* **61**, 86–100 (2009).

7. Vladar, E. K., Lee, Y. L., Stearns, T. & Axelrod, J. D. Observing planar cell polarity in multiciliated mouse airway epithelial cells. *Methods in cell biology* **127**, 37–54 (2015).

8. Sears, P. R., Yin, W.-N. & Ostrowski, L. E. Continuous mucociliary transport by primary human airway epithelial cells in *vitro*. *American Journal of Physiology-Lung Cellular and Molecular Physiology* **309**, L99–L108 (2015).

9. Faubel, R., Westendorf, C., Bodenschatz, E. & Eichele, G. Cilia-based flow network in the brain ventricles. *Science* **353**, 176–178 (2016).

10. Bennett, W. D. et al. Multisite comparison of mucociliary and cough clearance measures using standardized methods. *Journal of aerosol medicine and pulmonary drug delivery* **26**, 157–164 (2013).

11. Foster, W. M., Walters, D. M., Longphre, M., Macri, K. & Miller, L. M. Methodology for the measurement of mucociliary function in the mouse by scintigraphy. *Journal of applied physiology* **90**, 1111–1118 (2001).

12. Lay, J. C., Stang, M. R., Fisher, P. E., Yankaskas, J. R. & Bennett, W. D. Airway retention of materials of different solubility following local intrabronchial deposition in dogs. *Journal of aerosol medicine* **16**, 153–166 (2003).

13. Bhashyam, A. R. et al. A pilot study to examine the effect of chronic treatment with immunosuppressive drugs on mucociliary clearance in a vagotomized murine model. *PloS one* **7**, e45312 (2012).

14. Hoegger, M. J. et al. Assessing mucociliary transport of single particles *in vivo* shows variable speed and preference for the ventral trachea in newborn pigs. *Proceedings of the National Academy of Sciences* **111**, 2355–2360 (2014).

15. Donnelley, M. et al. High-resolution mucociliary transport measurement in live excised large animal trachea using synchrotron x-ray imaging. *Respiratory research* **18**, 95 (2017).
16. Gray, T. E., Guzman, K., Davis, C. W., Abdullah, L. H. & Nettesheim, P. Mucociliary differentiation of serially passaged normal human tracheobronchial epithelial cells. *American journal of respiratory cell and molecular biology* **14**, 104–112 (1996).

17. Gras, D. et al. An ex vivo model of severe asthma using reconstituted human bronchial epithelium. *Journal of Allergy and Clinical Immunology* **129**, 1259–1266 (2012).

18. Matsui, H. et al. Evidence for periciliary liquid layer depletion, not abnormal ion composition, in the pathogenesis of cystic fibrosis airways disease. *Cell* **95**, 1005–1015 (1998).

19. Tarran, R. et al. Normal and cystic fibrosis airway surface liquid homeostasis the effects of phasic shear stress and viral infections. *Journal of Biological Chemistry* **280**, 35751–35759 (2005).

20. Zhang, L. et al. Cftr delivery to 25% of surface epithelial cells restores normal rates of mucus transport to human cystic fibrosis airway epithelium. *PLoS biology* **7**, e1000155 (2009).

21. Button, B. & Boucher, R. C. of North Carolina Virtual Lung Group, U. et al. Role of mechanical stress in regulating airway surface hydration and mucus clearance rates. *Respiratory physiology & neurobiology* **163**, 189–201 (2008).

22. Thomas, B. et al. Ciliary dysfunction and ultrastructural abnormalities are features of severe asthma. *Journal of Allergy and Clinical Immunology* **126**, 722–729 (2010).

23. Chang, M. M.-J., Shih, L. & Wu, R. Pulmonary epithelium: cell types and functions. *The Pulmonary Epithelium in Health and Disease* **2**–26 (2008).

24. Parker, J. C. et al. Il-31 does not induce normal human ciliated epithelial cells to differentiate into a phenotype consistent with the pathophysiology of asthma. *Results in immunology* **2**, 104–111 (2012).

25. Adam, D. et al. Cystic fibrosis airway epithelium remodelling: involvement of inflammation. *The Journal of pathology* **235**, 408–419 (2015).

26. Raman, T. et al. Quality control in microarray assessment of gene expression in human airway epithelium. *BMC genomics* **10**, 493 (2009).

27. King, M., Gilboa, A., Meyer, F. & Silberberg, A. On the transport of mucus and its rheologic simulants in ciliated systems 1, 2. *American Review of Respiratory Disease* **110**, 740–745 (1974).

28. Hofmann, W. & Asgharian, B. The effect of lung structure on mucociliary clearance and particle retention in human and rat lungs. *Toxicological Sciences* **73**, 448–456 (2003).

29. Lafkas, D. et al. Therapeutic antibodies reveal notch control of transdifferentiation in the adult lung. *Nature* **528**, 127 (2015).

30. Francis, R. J. et al. Initiation and maturation of cilia-generated flow in newborn and postnatal mouse airway. *American Journal of Physiology-Lung Cellular and Molecular Physiology* **296**, L1067–L1075 (2009).

31. Marchetti, M. C. et al. Hydrodynamics of soft active matter. *Reviews of Modern Physics* **85**, 1143 (2013).

32. Ramaswamy, S. The mechanics and statistics of active matter. *Annual Review of Condensed Matter Physics* **1**, 323–345 (2010).

33. Vachier, I. et al. Glucocorticoid receptors in bronchial epithelial cells in asthma. *American Journal of respiratory and critical care medicine* **158**, 963–970 (1998).

34. Püspöki, Z., Storath, M., Sage, D. & Unser, M. Transforms and operators for directional bioimage analysis: a survey. In *Focus on Bio-Image Informatics*, 69–93 (Springer, 2016).

35. Staudt, M. R. et al. Smoking is associated with a loss of ciliated cells throughout the airways. In C26. TOBACCO: HOW IT HARMS US AND MORE REASONS TO QUIT, A4097–A4097 (Am Thoracic Soc, 2014).

## Acknowledgements

We thank Dr. J. Dupire and C. Garulli for preliminary experiments, Dr. I. Vachier, Dr. A. Petit, Pr. A. Bourdin, Dr. G. Massiera and Dr. J. Favier for fruitful discussions. Research was supported by the ANR MUCOCIL project, grant ANR-13-BSV5-0015 of the French Agence Nationale de la Recherche (M-K. K., E.L., P.C., D.G. and A.V.). The research leading to these results has received funding from the People Programme (Marie Curie Actions) of the European Union's Seventh Framework Programme (FP7/2007-2013) under REA grant agreement n. PCOFUND-GA-2013-609102, through the PRESTIGE programme coordinated by Campus France (E.L.). M. K. Khelloufi thanks the company MedBioMed for financial support.

## Author Contributions

A.V. and D.G. contributed equally to this work, A.V., D.G., and P.C. planned the experiment. A.V. and M.J. developed the theoretical model; D.G. and P.C. provided the epithelial cultures; M.-K.K., and E.L. performed the experiments; N.M., M.-K.K., E.L. and A.V. performed data analysis; A.V. and E.L. wrote the paper.

## Additional Information

Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-018-20882-4.

**Competing Interests:** The authors declare that they have no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_10.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2018