Understanding the Collective Dynamics of Motile Cilia in Human Airways

Eukaryotic organisms rely on the coordinated beating of motile cilia for a multitude of fundamental reasons. In smaller organisms, such as Paramecium and the single cell alga Chlamydomonas reinhardtii, it is a matter of propulsion, to swim towards a higher concentration of nutrients or away from damaging environments. Larger organisms use instead the coordinated motion of cilia to push fluid along an epithelium: examples common to mammals are the circulation of cerebrospinal fluid in the brain, the transport of ovules in the fallopian tubes, and breaking the left/right symmetry in the embryo. Another notable example, and one that is central to this thesis, is mucociliary clearance in human airways: A carpet of motile cilia helps keeping the cell surface free from pathogens and foreign particles by constantly evacuating from lungs, bronchi, and trachea a barrier of mucus. The question of how motile cilia interact with one another to beat in a coordinated fashion is an open and pressing one, with immediate implications for the medical community. In order for the fluid propulsion to be effective, the motion of cilia needs to be phase-locked across significant distances, in the form of travelling waves (``metachronal waves''). It is still not known how this long-range coordination emerges from local rules, as there is no central node regulating the coordination among cilia. In the first part of this thesis I will focus on studying the coordination in carpets of cilia with a top-down approach, by proposing, implementing, and applying a new method of analysing microscope videos of ciliated epithelia. Chapter 1 provides the reader with an introduction on motile cilia and flagella, treating their structure and motion and reporting the different open questions currently tackled by the scientific community, with particular interest in the coordination mechanisms of cilia and the mucociliary clearance apparatus. Chapter 2 introduces Differential Dynamic Microscopy (DDM), a powerful and versatile image analysis tool that bridges the gap between spectroscopy and microscopy by allowing to perform scattering experiments on a microscope. The most interesting aspects of DDM for this work are that it can be applied to microscope videos where it is not possible to resolve individual objects in the field of view, and it requires no user input. These two characteristics make DDM a perfect candidate for analysing several hundred microscope videos of weakly scattering filaments such as cilia. In Chapter 3 I will present how it is possible to employ DDM to extract a wealth of often-overlooked information from videos of ciliated epithelia: DDM can successfully probe the ciliary beat frequency (CBF) in a sample, measure the direction of beating of the cilia, and detect metachronal waves and read their direction and wavelength. In vitro ciliated epithelia however often do not show perfect coordination or alignment among cilia. For the analysis of these samples, where the metachronal coordination might not be evident, we developed a new approach, called multiscale DDM (multiDDM), to measure a coordination length scale, a characteristic length of the system over which the coordination between cilia is lost. The new technique of multiDDM is employed in Chapter 4 to study how the coordination among cilia changes as a response to changes in the rheology of the mucous layer. In particular, we show that cilia beating under a thick, gel-like mucus layer show a larger coordination length scale, as if the mucus acted as an elastic raft effectively coupling cilia over long distances. This is corroborated by the coordination length scale being larger in samples from patients affected by Cystic Fibrosis than in healthy samples, and much shorter when the mucus layer is washed and cilia therefore beat in a near-Newtonian fluid. We then show how it is possible to employ multiDDM to measure the effectiveness of drugs in recovering, in CF samples, a coordination length scale typical of a healthy phenotype. In the second part I will focus instead on the single cilium scale, showing how we can attempt to link the beating pattern of cilia to numerical simulations studying synchronisation in a model system. In particular in Chapter 5 I will describe our approach to quantitatively describe the beating pattern of single cilia obtained from human airway cells of either healthy individuals or patients affected by Primary Ciliary Dyskinesia. Our description of the beating pattern, and the selection of a few meaningful, summary parameters, are then shown to be accurate enough to discriminate between different mutations within Primary Ciliary Dyskinesia. In Chapter 6 instead I report the results obtained by coarse-graining the ciliary beat pattern into a model system consisting of two ``rotors''. The rotors are simulated colloidal particles driven along closed trajectories while leaving their phase free. In my study, the trajectories followed by the rotors are analytical fits of experimental trajectories of the centre of drag of real cilia. The rotors, that are coupled only via hydrodynamics interactions, are seen to phase-lock, and the shape of the trajectory they are driven along is seen to influence the steady state of the system

Soft pinning of liquid domains on topographical hemispherical caps

The role of lipid composition as a regulator or mediator of processes that take place in biological membranes is a very topical question, and important insights can be gained by studying in-vitro model lipid mixture systems. A particular question is the coupling of local curvature to the local phases in membranes of mixed composition. Working with an experimental system of giant unilamellar vesicles of ternary composition, the curvature is imposed by approaching the membrane to a topographically (on the micron scale) patterned surface. Performing experiments, we show that domains of the more disordered phase localise preferentially to regions of higher curvature. We characterise and discuss the strength of this "caging" behaviour. In future, the setup we discuss here could prove useful as a platform to localise domains rich in membrane proteins, or to promote the onset of biochemical processes at specific locations. Finally, we note that the methods developed here could have also applications in bio-sensing, as a similar -but metal coated- topography can sustain plasmonic resonances

Code supporting "Phenotyping ciliary dynamics and coordination in response to CFTR-modulators in Cystic Fibrosis respiratory epithelial cells"

MATLAB code to run multiDDM analysis. See the associated README file for further information

Quantitative high-speed video profiling discriminates between <i>DNAH11</i> and <i>HYDIN</i> variants of primary ciliary dyskinesia

Primary ciliary dyskinesia (PCD) is a genetically heterogeneous condition. The number of genes associated with primary ciliary dyskinesia is rising, and the link between genotype and clinical phenotype is largely unknown (1). Recent advances in molecular genotyping have helped establish the association between genetic defects and aberrant cilia ultrastructure in PCD, as detected by transmission electron microscopy. However, up to 30% of PCD cases do not show any specific ultrastructural defects and thus cannot be detected using this method. A classic example of this occurs in PCD cases caused by mutations in the DNAH11 (dynein axonemal heavy chain 1) protein, where aberrant ciliary beating can be detected via high-speed video microscopy (HSVM) analysis, but no ultrastructural defect is observed (2, 3). Even in cases where ultrastructural defects do exist, they can be difficult to detect using standard transmission electron microscopy techniques, and therefore patients who harbor such subtle defects are at risk of being misdiagnosed. PCD caused by mutations in the protein HYDIN, for example, display normal ciliary beat frequency and rarely exhibit abnormal ciliary transmission electron microscopy, yet mucociliary clearance is compromised

Behavioral fingerprints predict insecticide and anthelmintic mode of action

Abstract Novel invertebrate‐killing compounds are required in agriculture and medicine to overcome resistance to existing treatments. Because insecticides and anthelmintics are discovered in phenotypic screens, a crucial step in the discovery process is determining the mode of action of hits. Visible whole‐organism symptoms are combined with molecular and physiological data to determine mode of action. However, manual symptomology is laborious and requires symptoms that are strong enough to see by eye. Here, we use high‐throughput imaging and quantitative phenotyping to measure Caenorhabditis elegans behavioral responses to compounds and train a classifier that predicts mode of action with an accuracy of 88% for a set of ten common modes of action. We also classify compounds within each mode of action to discover substructure that is not captured in broad mode‐of‐action labels. High‐throughput imaging and automated phenotyping could therefore accelerate mode‐of‐action discovery in invertebrate‐targeting compound development and help to refine mode‐of‐action categories

Behavioral fingerprints predict insecticide and anthelmintic mode of action

Novel invertebrate‐killing compounds are required in agriculture and medicine to overcome resistance to existing treatments. Because insecticides and anthelmintics are discovered in phenotypic screens, a crucial step in the discovery process is determining the mode of action of hits. Visible whole‐organism symptoms are combined with molecular and physiological data to determine mode of action. However, manual symptomology is laborious and requires symptoms that are strong enough to see by eye. Here, we use high‐throughput imaging and quantitative phenotyping to measure Caenorhabditis elegans behavioral responses to compounds and train a classifier that predicts mode of action with an accuracy of 88% for a set of ten common modes of action. We also classify compounds within each mode of action to discover substructure that is not captured in broad mode‐of‐action labels. High‐throughput imaging and automated phenotyping could therefore accelerate mode‐of‐action discovery in invertebrate‐targeting compound development and help to refine mode‐of‐action categories

Single-Cell and Population Transcriptomics Reveal Pan-epithelial Remodeling in Type 2-High Asthma.

The type 2 cytokine-high asthma endotype (T2H) is characterized by IL-13-driven mucus obstruction of the airways. To further investigate this incompletely understood pathobiology, we characterize IL-13 effects on human airway epithelial cell cultures using single-cell RNA sequencing, finding that IL-13 generates a distinctive transcriptional state for each cell type. Specifically, we discover a mucus secretory program induced by IL-13 in all cell types which converts both mucus and defense secretory cells into a metaplastic state with emergent mucin production and secretion, while leading to ER stress and cell death in ciliated cells. The IL-13-remodeled epithelium secretes a pathologic, mucin-imbalanced, and innate immunity-depleted proteome that arrests mucociliary motion. Signatures of IL-13-induced cellular remodeling are mirrored by transcriptional signatures characteristic of the nasal airway epithelium within T2H versus T2-low asthmatic children. Our results reveal the epithelium-wide scope of T2H asthma and present candidate therapeutic targets for restoring normal epithelial function

Single-Cell and Population Transcriptomics Reveal Pan-epithelial Remodeling in Type 2-High Asthma.

The type 2 cytokine-high asthma endotype (T2H) is characterized by IL-13-driven mucus obstruction of the airways. To further investigate this incompletely understood pathobiology, we characterize IL-13 effects on human airway epithelial cell cultures using single-cell RNA sequencing, finding that IL-13 generates a distinctive transcriptional state for each cell type. Specifically, we discover a mucus secretory program induced by IL-13 in all cell types which converts both mucus and defense secretory cells into a metaplastic state with emergent mucin production and secretion, while leading to ER stress and cell death in ciliated cells. The IL-13-remodeled epithelium secretes a pathologic, mucin-imbalanced, and innate immunity-depleted proteome that arrests mucociliary motion. Signatures of IL-13-induced cellular remodeling are mirrored by transcriptional signatures characteristic of the nasal airway epithelium within T2H versus T2-low asthmatic children. Our results reveal the epithelium-wide scope of T2H asthma and present candidate therapeutic targets for restoring normal epithelial function

A common polymorphism in the Intelectin-1 gene influences mucus plugging in severe asthma

By incompletely understood mechanisms, type 2 (T2) inflammation present in the airways of severe asthmatics drives the formation of pathologic mucus which leads to airway mucus plugging. Here we investigate the molecular role and clinical significance of intelectin-1 (ITLN-1) in the development of pathologic airway mucus in asthma. Through analyses of human airway epithelial cells we find that ITLN1 gene expression is highly induced by interleukin-13 (IL-13) in a subset of metaplastic MUC5AC+ mucus secretory cells, and that ITLN-1 protein is a secreted component of IL-13-induced mucus. Additionally, we find ITLN-1 protein binds the C-terminus of the MUC5AC mucin and that its deletion in airway epithelial cells partially reverses IL-13-induced mucostasis. Through analysis of nasal airway epithelial brushings, we find that ITLN1 is highly expressed in T2-high asthmatics, when compared to T2-low children. Furthermore, we demonstrate that both ITLN-1 gene expression and protein levels are significantly reduced by a common genetic variant that is associated with protection from the formation of mucus plugs in T2-high asthma. This work identifies an important biomarker and targetable pathways for the treatment of mucus obstruction in asthma