Single-cell transcriptomics of a dynamic cell behavior in murine airways
Abstract
Despite advances in high-dimensional cellular analysis, the molecular profiling of dynamic behaviors of cells in their native environment remains a major challenge. We present a method that allows us to couple the physiological behaviors of cells in an intact murine tissue to deep molecular profiling of individual cells. This method enabled us to establish a novel molecular signature for a striking migratory cellular behavior following injury in murine airways.
Editor's evaluation
This study presents a useful combination of live cell imaging with single-cell transcriptomic analyses. This is a first step to further expanding the description of cellular heterogeneity, including cellular behavior as well as gene expression profiles. The method, with more technical improvements, will be of interest to researchers who study dynamic changes in cell morphology and gene expression.
https://doi.org/10.7554/eLife.76645.sa0Introduction
Cells in a living organism are dynamic entities, changing their characteristics over space and time and constantly interacting with the host and pathogens. The ability to obtain such information and link it to detailed molecular phenotypes of the cells would be highly useful for biomedical investigations but has been underappreciated. Here, we present a method that allows us to characterize complex physiologic behaviors of cells in an intact tissue and then perform live imaging-guided sequencing of the same cells. We validate this approach using a regeneration model of airway tissues and demonstrate how this method leads to new biological findings.
There is a pressing need for a comprehensive understanding of cellular behaviors in the lung, the site where aberrant cellular behavior has been linked to asthma (Kim et al., 2012; Park et al., 2015) pulmonary fibrosis (Fukumoto et al., 2016), and viral infections including influenza and coronaviruses (Kumar et al., 2011). Single-cell RNA-sequencing (scRNA-seq) has emerged as a precise way to define cell type and cell state, and new techniques are being developed to determine the spatial distribution of sequenced cells in tissues (Marx, 2021). However, the molecular pathways that drive the cellular behavior in situ continue to be inferred from time-lapse tissue sampling or transcriptional kinetics (La Manno et al., 2018). Moving beyond inference requires coupling visualized in situ cell behavior with deep molecular profiling of visualized cells.
Live cell imaging is an established technique for capturing morphology and cellular dynamics such as cellular migration during skin regeneration (Park et al., 2019; Park et al., 2017), but imaging in the lung remains challenging due to difficult access and the constant motion of the respiratory system. Additionally, molecular information that accompanies live imaging is largely limited to a few fluorescent reporters. Prior attempts to link deep molecular profiling with live imaging have relied on imaging dissociated cells (Lane et al., 2017; Yuan et al., 2018), cell monolayers (Hu et al., 2020), organoids (Konen et al., 2017) rather than on cell behaviors in their native tissue environment.
Our aim is to bridge the divide between two powerful methodologies - cell behavioral observation through live imaging and transcriptional profiling through single-cell sequencing, ultimately allowing the identification of a transcriptional signature that corresponds to that cell behavior. We present an important step forward toward linking dynamic cell behaviors with single-cell transcriptomics.
Results
We describe a novel approach to linking live tissue imaging with single-cell profiling (Figure 1a). In order to visualize the airway epithelium at high resolution over days, we explant a mouse trachea and secure it in a custom imaging platform without tissue submersion. This approach minimizes sample movement during imaging, maintains a constant supply of nutrients from below the explant, and preserves an air-liquid interface (ALI), which is required for the maintenance of the normal cellular architecture of the airway epithelium (You et al., 2002; Figure 1a and Figure 1—figure supplement 1a). This platform allows imaging of common and rare cell types in the airway epithelium at high resolution in their native environment (Figure 1—figure supplement 1b). The explant culture also allows an uninjured tracheal epithelium to survive with its native cellular anatomy for weeks with daily high-resolution imaging (Figure 1—figure supplement 1c).
Discernable cell behaviors have a broad time scale, ranging from milliseconds to days. Thus, we imaged a wide range of cellular behaviors, from rapid fluctuations of ciliary beating and directional mucociliary transport over milliseconds to wholescale regeneration of the airway epithelium, which occurs over days after the injury (Figure 1, c to f, ). Furthermore, this method allows single-cell level registration within tissues that are live-imaged and subsequently fixed and stained, which enables a unique comparison between live fluorescence cellular patterns and immunostains that describe cell identity and function (Figure 1—figure supplement 1d).
Remarkably, this airway imaging platform faithfully recapitulates and captures cellular dynamics of epithelial regeneration from native basal stem cells after an extensive epithelial injury induced by sulfur dioxide (SO2) (Figure 1f and Figure 2—figure supplement 1). In the first 5 days, the basal cells divide, increasing cellular density and reforming the pseudostratified epithelium. In the next 5–10 days, the epithelium differentiates, leading to the restoration of the full epithelium, including the regeneration of ciliated cells. Complete regeneration requires an air-liquid interface (Figure 1f). Overall, this murine trachea explant ALI culture retains the nearly complete 3D organization and microenvironment of the basal progenitor cells and, therefore, offers a unique model to study organ physiology and regeneration outside of the body.
Continuous time-lapse imaging of the airway epithelium for up to 80 hr after injury (Figure 2a and Figure 2—figure supplement 2a and Video 1) demonstrated changes in cellular architecture over time, including an increase in the average cellular density and epithelial thickness, without apparent phototoxicity (Figure 2—figure supplement 2a). Imaging of trachea explant controls from uninjured mice over 19 hr revealed no cellular displacement in the airway epithelium (Video 2). We examined cell movement using single-cell tracking following segmentation of cell nuclei and particle image velocimetry (PIV) of non-segmented images (Figure 2—figure supplement 2b,c). These analyses revealed a variety of regeneration cellular dynamics inaccessible without live imaging. For example, Hertwig’s Rule predicts that a cell division plane is perpendicular to the long axis of the cell during the preceding interphase (Minc and Piel, 2012). This was established in plants (Besson and Dumais, 2011) and developing simple model organisms (Aigouy et al., 2010; Concha and Adams, 1998; Tsou et al., 2003) but has never been probed in an adult regenerating tissue. We found that the long axis in most cells predicts the cell division axis, while the axis of cellular movement prior to cell division does not (Figure 1e and Figure 2—figure supplement 3a,b).
We also found a surprising degree of heterogeneity of collective cellular migration during regeneration throughout the injured airways. In regions that demonstrated rapid cellular movement after an injury, the movement peaked at 26–38 hr after SO2 injury and the speed declined significantly in most regions by 50 hr after injury (Figure 2c). There was a significant interaction between time and the mean speed, but no significant difference between the mouse and mean speed (Figure 2c). Variable migratory behavior of airway epithelial cells has been observed in cell culture models (Kim et al., 2020; Park et al., 2015) but not previously in an intact regenerating airway tissue. We found that the frequency distribution of cellular speed in different regions at 26–38 hr demonstrated large variability ranging from ‘non-mover’ regions (<1.5 μm/hr) to ‘mover’ regions (>4 μm/hr) (Figure 2d). Furthermore, videos with higher temporal resolution revealed that the cells with much slower movements in the ‘non-mover’ regions had no directional preference and oscillated in place, whereas the ‘mover’ regions with faster cellular movements were more uni-directional with packs of cells migrating in waves (Figure 2—figure supplement 3c,d and Video 3).
Distinct subsets of cells have been theorized to contribute to the regeneration process (Pardo-Saganta et al., 2015; Tadokoro et al., 2014), but it is unclear whether these heterogeneous transcriptional cell states reflect gene expression stochasticity or correlate with unique cell behaviors. To determine the molecular signatures of cells with directional movement compared to regenerating non-moving cells, we marked epithelial regions by photoconversion using the Kaede transgenic mouse model, in which cells change from green to red after exposure to violet light (Tomura et al., 2008). We photoconverted cells at 24 hr after injury, imaged every 6 hr, screened for ‘movers’ with >50 μm displacement (>3 μm/hr) at 18 hr after photoconversion, and then isolated photoconverted epithelial cells by FACS for plate-based scRNA-seq (Figure 3a and Materials and methods). Dimensionality reduction revealed that cells from a ‘moving’ region (M) and a ‘non-moving’ (NM) region cluster separately (Figure 3b). Using unsupervised clustering and cell identity signatures (Methods, Figure 3—figure supplement 1), we found that nearly all the cells in the M region are basal cells, whereas the NM region contains basal and club cells (Figure 3c). We identified the differences in gene expression between the M basal cells and NM basal cells and identified gene signatures that are enriched (FDR <0.05, likelihood-ratio test) either in the M or the NM basal cells (Figure 3c and Figure 3—figure supplement 2).
We hypothesized that the identified phenotypes are a common feature of injury-induced epithelial regeneration. We examined published data of an independent injury model (Plasschaert et al., 2018) and analyzed the prevalence of these signatures during repair after polidocanol injury. As predicted, the M basal cell signature is strongly enriched 24 hr post injury (hpi), declines at 48 hpi and 72 hpi, and returns to baseline at 1 week after injury (all p<10–16, Mann-Whitney U test, Figure 4a). Similarly, the NM signature is decreased at 24 hpi when cell migration is presumed to be active, increases at 48 hpi and 72 hpi when cell migration is presumed to be diminished, and returns to baseline at 1 week post-injury when regeneration is complete (Figure 4a). Furthermore, at 24 hpi we found that scoring basal cells using M and NM signatures segregated basal cells into two statistically distinct cell populations (Figure 4b), indicating that polidocanol regeneration is likely also characterized by these cell phenotypes. To test this possibility, we used unsupervised clustering (Methods) to define two groups of basal cells at 24 hpi and found that these two populations were indeed separately enriched for the M and NM basal cell signatures (Figure 4c), confirming the presence of distinct M and NM basal cells during polidocanol regeneration. Taken together, these findings suggest that distinct M and NM cell behaviors are conserved features of early epithelial regeneration and demonstrate that our live imaging-guided single-cell profiling approach can discover generalizable principles of tissue biology.
Discussion
The rapid progress in single-cell analysis is enabling the discovery and characterization of transcriptionally heterogeneous cells in diverse tissue contexts (Lee et al., 2021; Ståhl et al., 2016). However, these methods do not capture the dynamics of cell behaviors that often define the unique biological processes that occur in the tissues. To address this gap, we developed an approach to examine the association of molecular and behavioral phenotypes of cells in their native tissues. We established an airway organ explant culture that maintains tissue dynamics for an extended length of time and combined this platform with live imaging in order to observe distinct airway cellular behaviors at a broad time scale, spanning cell migration, cell division, and ciliary beating.
To link cell behavior to molecular analysis, we used photoconversion to mark cells for subsequent single-cell genomics analysis. We found that a subpopulation of basal stem cells migrates within the lung during early regeneration. We used recently developed single-cell RNA-sequencing approaches to establish molecular signatures for moving and nonmoving basal cells. These distinct cell signatures were identified across independent lung regeneration models, suggesting that M and NM cell behaviors are conserved cellular features of early epithelial regeneration.
There are important limitations to this approach. For example, the optical requirements of the Kaede photoconversion model limit its application for single cell photoconversion. However, a similar approach can take advantage of different methods of cell labeling and live imaging, including genetic labeling of single cells with fluorescent reporters (Figure 1—figure supplement 1b). The requirement for fluorescence is another constraint of this approach, although label-free imaging approaches may be used to visualize the behavior of individual cells in native tissues (Shah et al., 2022). Another consideration is that long-term live imaging itself can influence gene expression and cell behavior. Although phototoxicity that leads to photobleaching or dramatic changes in cell integrity can be avoided, subtle changes to gene expression or cell behavior due to live imaging may be difficult to detect (Magidson and Khodjakov, 2013).
This report focuses on cellular behavior in response to injury in the airway epithelium. Other epithelial organs such as the cornea (Park et al., 2019) and esophagus Doupé et al., 2012; Yokobori et al., 2016 have been cultured on similar platforms, and we anticipate this live imaging-guided single-cell profiling approach can be extended to other tissues to discover underlying principles of heterogeneous cellular behaviors at homeostasis and in disease. Although the explant platform lacks some elements of the in vivo microenvironment of these organs, live imaging of the explant followed by molecular analysis overcomes the costly constraints of spatial and temporal resolution during in vivo live imaging (Haase et al., 2022), the type of resolution required to link dynamic cell behaviors with single cell transcriptomics.
Materials and methods
Mice
mT-mG (stock no. 007676), nT-nG (stock no. 023035), CAGs-LSL-rtTA3 (stock no. 029617), and Col1a1-tetO-H2B-mCherry (stock no. 014602), CD11cCre (stock no. 007567), and Ascl1nGFP (stock no 012881) mice were purchased from the Jackson Laboratory. Foxj1Cre (Zhang et al., 2007), KRT5rtTA (Diamond et al., 2000), B1EGFP (Miller et al., 2005), Foxj1CreER (Rawlins and Hogan, 2008), CC10CreER (Rawlins et al., 2009), and Kaede (Tomura et al., 2008) lines were previously described. A line of Membrane-GFP (mG) mice was generated by selecting GFP-positive pups of a Foxj1Cre-mTmG male parent (with mT to mG recombination in the sperm) and backcrossing to WT background to eliminate the Cre allele. The mG line without Cre was crossed to nT-nG to generate the ‘nT-mG’ strain. Mice were maintained in an Association for Assessment and Accreditation of Laboratory Animal Care-accredited animal facility at the Massachusetts General Hospital, and procedures were performed with Institutional Animal Care and Use Committee (IACUC)-approved protocols. Mice of all strains were housed in an environment with controlled temperature and humidity, on 12 hr light-dark cycles, and fed with regular rodent’s chow.
SO2 injury
Request a detailed protocolSO2 injury model was performed as previously described (Kim et al., 2012; Pardo-Saganta et al., 2015). In brief, mice were exposed to 500 p.p.m. of SO2 for 3 hr 40 min and the trachea was collected 16–24 hr after injury for imaging and explant culture.
Tracheal explant and tissue ALI culture
Request a detailed protocolTracheas were dissected, cleared of connective tissue and adjacent organs, and opened longitudinally along the anterior tracheal wall. The tracheas were placed on ice in DMEM/F-12 Media with Primocin (InVivoGen) and 15 mM HEPES until culture. Trachea explants were then sutured onto a silicone o-ring and placed in a custom-made tissue culture dish over an inverted ALI insert secured in a 60 mm tissue culture dish by PDMS. This approach ensured stability during high-resolution imaging. The media contacted the explant from below through the ALI membrane. The dish was placed in a physiological live imaging chamber (CO2 and temperature-controlled, TokaiHit) on the stage of the two-photon microscope.
Physiological two-photon imaging
Request a detailed protocolTrachea explant imaging was performed on an Olympus FVMPE-RS multiphoton laser scanning microscope equipped with a MaiTai HPDS-O IR pulsed laser (900 nm for GFP and SHG) and INSIGHT X3-OL IR pulsed laser (1100 nm for tdTomato), using a 25 X water immersion lens (NA 1.05). Explants were imaged at time points as indicated in the Figures. For orthogonal view reconstruction, we scanned the trachea with 0.75 µm Z steps. To reimage the same trachea at high resolution at different time points, landmarks such as cartilage rings and vascularity patterns were used as fiducial marks. These fiducial marks were also used for 2D and 3D registration of different time points.
Image analysis
Request a detailed protocol4D images (x, y, z, t) were imported into MATLAB and/or ImageJ for image processing and analysis. Because the curvature of the tissue changes over time, we first normalized each 3D image to generate a flat basement membrane. As the SHG signal is maximal at the basement membrane, we computed the z height of the basement membrane across the image after applying a Gaussian blur (typical σ values: 10–25 μm in xy, 1–4 μm in z). This height was subtracted from the original 3D data to level the basement membrane. MATLAB code for flattening the 3D images is publicly available on Github at https://github.com/sheldonjjlase/matlab/, (copy archived at swh:1:rev:230ca1a3c03039ae553d88c3d3f1286404a34a85) (Kwok, 2022). Other image processing steps including brightness and contrast adjustments, background subtraction, photobleaching correction, pseudocoloring, 3D time-lapse registration, and stitching were performed using built-in functions in ImageJ.
Cilia beating was recorded by acquiring time-lapse two-photon images of the epithelial surface at 150 Hz over 200 frames using a resonant galvanometer scanner. To estimate the cilia beat frequency (CBF), we estimated the power spectral density of the fluorescence intensity fluctuations across the image using Welch’s method in MATLAB. The peak fluctuation frequency was computed for each pixel across the image corresponding to bright cilia. Mucociliary transport was measured by applying 1 μm fluorescent spherical beads to the epithelial surface and recording their displacement over time after equilibration.
To track individual cells over time-lapse imaging, images were imported into ilastik for segmentation and cell tracking (Berg et al., 2019). Briefly, pixels corresponding to nuclei were first classified using manual training and machine learning. Next, individual cells were similarly identified through manual training and machine learning algorithms to classify objects. Finally, classified cells were tracked over time using a conservation tracking algorithm. Segmented and tracked cells were then imported into MATLAB for quantitative analysis, including the computation of individual cell speeds over time.
For the automated estimation of cell speed from time-lapse imaging, we performed automated particle imaging velocimetry (PIV). Image sequences were imported into MATLAB and analyzed using the PIVlab plugin (Thielicke and Stamhuis, 2014). A direct Fourier transforms correlation with multiple passes of sizes consisting of 24 μm, 16 μm, and 10 μm was used. This generated displacement vectors arranged in a grid with 10 μm spacing for each sequential pair of images. The average cell speed for each time-point was estimated by computing the average absolute displacement estimated by PIV in each (10 × 10) μm2 region divided by the time between images. To quantify the directionality of cell movement, we computed the circular variance of the displacement vectors generated by PIV analysis.
Kaede photoconversion
Request a detailed protocolTrachea from Kaede mice (Tomura et al., 2008) were explanted 20 hr after SO2 inhalation injury, sutured onto a silicone O-ring, and secured on an inverted ALI insert in media on ice, and placed on the imaging platform of a FV3000 Olympus Laser Scanning confocal microscope. Selected regions were outlined and photoconverted using the 405 nm laser for 2 min, while both disappearances of KaedeGreen and the appearance of KaedeRed were simultaneously visualized using the 488 nm and the 561 nm lasers, respectively.
To identify regions of movement and no movement, we explanted Kaede tracheas 20 hr after SO2 injury, photoconverted distinct regions with a specific shape, and proceeded with timelapse live imaging, screening for regions with significant shape displacement over time (from epithelial movement) (Figure 3a). After defining whether a region moved or remained non-moving, we excised a trachea fragment, dissociated the fragment into single cells, and used flow-activated cell sorting (FACS) to isolate photoconverted (KaedeRed) epithelial cells. We then proceeded to single-cell RNA sequencing of cells isolated from moving and non-moving regions.
Cell dissociation and FACS
Request a detailed protocolAirway epithelial cells were dissociated using papain solution. Tracheal fragments with photoconverted regions were trimmed and incubated in papain dissociation solution and incubated at 37 °C for 2 hr. After incubation, dissociated tissues were passed through a cell strainer and centrifuged, and pelleted at 500 g for 5 min. Cell pellets were dispersed and incubated with Ovo-mucoid protease inhibitor (Worthington Biochemical, cat. no. LK003182) to inactivate residual papain activity by incubating on a rocker at 4 °C for 20 min. Cells were then pelleted and stained with EpCAM–BV421 (1:50; BD Bioscience, #563214) for 30 min in 2.5% FBS in PBS on ice. After washing, cells were sorted by fluorescence (antibody staining, Kaede-Green, and Kaede-Red) on a BD FACS Aria (BD Biosciences) using FACS Diva software, and analysis was performed using FlowJo (version 10) software.
Single cells were sorted into each well of a 96-well PCR plate containing 5 μl buffer. After sorting, the plate was sealed with a Microseal F, centrifuged at 800 g for 1 min and immediately frozen on dry ice. Plates were stored at −80°C and submitted to a core facility for cDNA library generation, amplification, and sequencing.
Single-cell sequencing and sequence analysis
Request a detailed protocolcDNA was generated from single cells in the 96-well plate using the SmartSeq v4 kits (Takara Bio) using 1/4th volume reactions dispensed using a Mantis dispenser (Formulatrix). Samples were amplified using 18 cycles of PCR. Resulting cDNA was then made into Illumina-compatible libraries using the Nextera XT kit (Illumina Inc). Libraries were sequenced on a NextSeq using a Mid Output 150 cycle kit (Illumina Inc) using 75 bp paired-end reads.
Pre-processing of plate-based scRNA-seq data
Request a detailed protocolBCL files were converted to merged, de-multiplexed FASTQ files using the Illumina Bcl2Fastq software package v.2.17.1.14. Paired-end reads were mapped to the UCSC mm10 mouse transcriptome using Bowtie (Langmead et al., 2009) with parameters ‘-q–phred33-quals -n 1 -e 99999999 l 25 -I 1 X 2000 a -m 15 S -p 6’, which allows alignment of sequences with one mismatch. Expression levels of genes were quantified as transcript-per-million (TPM) values by RSEM (Li and Dewey, 2011) v.1.2.3 in paired-end mode. For each cell, we determined the number of genes for which at least one read was mapped, and then excluded all cells with fewer than 1,000 or more than 10,000 detected genes, or less than 25% of reads mapping to the transcriptome.
To identify variable genes a logistic regression was fit to the cellular detection fraction, using the total number of transcripts per cell as a predictor. Outliers from this curve are genes that are expressed in a lower fraction of cells than would be expected given the total number of transcripts mapping to that gene, that is, cell-type or state-specific genes. We used a threshold of deviance <−0.15, producing a set of 1910 variable genes.
Dimensionality reduction by PCA and t-SNE
Request a detailed protocolWe restricted the expression matrix to the subsets of variable genes and high-quality cells noted above, and values were log2-transformed, and then centered and scaled before input to PCA, which was implemented using the R function ‘prcomp’ from the ‘stats’ package. After PCA, significant principal components were identified by inspection of the scree plot. Only scores from the first 20 PCs were used as the input to further analysis.
For visualization purposes only (and not for clustering), dimensionality was further reduced using the Barnes–Hut approximate version of t-SNE (Van Der Maaten, 2014; Figure 3b). This was implemented using the ‘Rtsne’ function from the ‘Rtsne’ R package.
To identify cell types within the data, hierarchical clustering was used using the ‘Ward.D2’ metric in the ‘hclust’ R package. Genes were filtered to epithelial cell type marker genes (Montoro et al., 2018) before input to the clustering algorithm. Pearson’s correlation was used as a distance metric. This produced three clusters, two were clearly identifiable as Basal and Club cells, based on the disjoint expression of known markers Krt5 and Scgb1a1, respectively, while the third was distinguished by much lower technical quality (an average of 2373 genes detected per cell compared to 5193 for the Basal and 5480 for the club clusters, respectively, p=0.0004, Mann-Whitney U-test). These low-quality cells were not used for DE testing.
To identify the signature of moving vs non-moving basal cells (Figure 3c) we ran differential expression tests between cells in the Basal cluster between the two conditions (moving and non-moving), and selected genes that were differentially expressed (FDR <0.05). Differential expression tests were carried out using a two-part ‘hurdle’ model to control for both technical quality and mouse-to-mouse variation. This was implemented using the R package MAST (Finak et al., 2015), and p values for differential expression were computed using the likelihood-ratio test. Multiple hypothesis testing corrections was performed by controlling the false discovery rate using the R function ‘p.adjust.’.
Re-analysis of polidocanol injury dataset
Request a detailed protocolPreviously published single-cell RNA sequencing data from mouse trachea injured using polidocanol (Plasschaert et al., 2018) was downloaded from the NCBI GEO (GSE102580). All available unique molecular identifier (UMI) counts tables from mice at 24, 48, 72, and 168 hr after injury along with uninjured controls were downloaded. Cell types were determined using the authors provided annotations. To determine the expression of migration-associated genes in the injury response, we scored the Basal cells for the set of genes (Figure 3c) both significantly up- (‘mover’) and down-regulated (‘non-mover’) (Figure 3d). Scoring cells was computed as described previously (Montoro et al., 2018). To obtain a score for a specific set of n genes in a given cell, a ‘background’ gene set was defined to control for differences in sequencing coverage and library complexity. The background gene set was selected for similarity to the genes of interest in terms of expression level. Specifically, the 10 n nearest neighbors in the 2D space defined by mean expression and detection frequency across all cells were selected. The signature score for that cell was then defined as the mean expression of the n signature genes in that cell, minus the mean expression of the 10 n background genes in that cell.
Unsupervised cluster analysis of polidocanol-injured basal cells 24 hr after injury was computed using default settings in Seurat. Briefly, variable genes were selected using the method ‘vst,’ and then PCA was computed using only these genes. Shared-nearest neighbor (SNN)-based clustering was implemented using the ‘FindClusters’ function (resolution parameter = 0.25) using the first 25 principal components as input, resulting in two clusters (Figure 4b).
Statistical analysis
Request a detailed protocolData were compared among groups using the Student’s t-test (unpaired, two-tailed) unless otherwise specified in the Figure legends. Analysis was performed with GraphPad Prism software (version 9.1.0).
Data availability
Sequencing data have been deposited in GEO under accession code GSE193954.
-
NCBI Gene Expression OmnibusID GSE193954. Single-cell transcriptomics of dynamic cell behaviors.
-
NCBI Gene Expression OmnibusID GSE102580. A single cell atlas of the airway epithelium reveals the CFTR-rich pulmonary ionocyte.
References
-
Ilastik: interactive machine learning for (BIO) image analysisNature Methods 16:1226–1232.https://doi.org/10.1038/s41592-019-0582-9
-
Conditional gene expression in the epidermis of transgenic mice using the tetracycline-regulated transactivators TTA and Rta linked to the keratin 5 promoterThe Journal of Investigative Dermatology 115:788–794.https://doi.org/10.1046/j.1523-1747.2000.00144.x
-
In vivo imaging of tracheal epithelial cells in mice during airway regenerationAmerican Journal of Respiratory Cell and Molecular Biology 47:864–868.https://doi.org/10.1165/rcmb.2012-0164OC
-
Unjamming and collective migration in MCF10A breast cancer cell linesBiochemical and Biophysical Research Communications 521:706–715.https://doi.org/10.1016/j.bbrc.2019.10.188
-
Circumventing photodamage in live-cell microscopyMethods in Cell Biology 114:545–560.https://doi.org/10.1016/B978-0-12-407761-4.00023-3
-
V-Atpase B1-subunit promoter drives expression of EGFP in intercalated cells of kidney, clear cells of epididymis and airway cells of lung in transgenic miceAmerican Journal of Physiology. Cell Physiology 288:C1134–C1144.https://doi.org/10.1152/ajpcell.00084.2004
-
Predicting division plane position and orientationTrends in Cell Biology 22:193–200.https://doi.org/10.1016/j.tcb.2012.01.003
-
Unjamming and cell shape in the asthmatic airway epitheliumNature Materials 14:1040–1048.https://doi.org/10.1038/nmat4357
-
Ciliated epithelial cell lifespan in the mouse trachea and lungAmerican Journal of Physiology. Lung Cellular and Molecular Physiology 295:L231–L234.https://doi.org/10.1152/ajplung.90209.2008
-
PIVlab – towards user-friendly, affordable and accurate digital particle image velocimetry in MatlabJournal of Open Research Software 2:1.https://doi.org/10.5334/jors.bl
-
PAR-dependent and geometry-dependent mechanisms of spindle positioningJournal of Cell Biology 160:845–855.https://doi.org/10.1083/jcb.200209079
-
Accelerating t-SNE using tree-based algorithmsJournal of Machine Learning Research: JMLR 15:3221–3245.
-
Intestinal epithelial culture under an air-liquid interface: a tool for studying human and mouse esophagiDiseases of the Esophagus 29:843–847.https://doi.org/10.1111/dote.12346
-
Growth and differentiation of mouse tracheal epithelial cells: selection of a proliferative populationAmerican Journal of Physiology. Lung Cellular and Molecular Physiology 283:L1315–L1321.https://doi.org/10.1152/ajplung.00169.2002
-
A transgenic foxj1-cre system for gene inactivation in ciliated epithelial cellsAmerican Journal of Respiratory Cell and Molecular Biology 36:515–519.https://doi.org/10.1165/rcmb.2006-0475RC
Article and author information
Author details
Funding
National Heart, Lung, and Blood Institute (P01HL120839)
- Seok-Hyun Yun
National Heart, Lung, and Blood Institute (F32HL154638)
- Daniel T Montoro
National Cancer Institute (R01EB033155)
- Seok-Hyun Yun
National Heart, Lung, and Blood Institute (K08HL124298)
- Vladimir Vinarsky
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Drs. Jin-ah Park, Jeffrey Fredberg, and Jeffrey Drazen for stimulating discussion about the air-liquid-interface tissue model and migratory behaviors of cells at the beginning of this project, which was partially funded by the National Institutes of Health (NIH) via P01HL120839. We thank all of the laboratories of Dr. Jayaraj Rajagopal for providing the cell type-specific fluorescent reporter mice, access to advanced microscopy, and critical appraisal of the project. We thank NIH R01EB033155, the HSCI-CRM Flow Cytometry and Microscopy Core Facility at the Massachusetts General Hospital, and The Bauer Core Facility at Harvard University for facilitating this project.
Ethics
Mice were maintained in an Association for Assessment and Accreditation of Laboratory Animal Care-accredited animal facility at the Massachusetts General Hospital, and procedures were performed with Institutional Animal Care and Use Committee (IACUC)-approved protocol 2009N000119.
Copyright
© 2023, Kwok, Montoro, Haber et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,812
- views
-
- 190
- downloads
-
- 3
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Developmental Biology
The morphogen FGF8 establishes graded positional cues imparting regional cellular responses via modulation of early target genes. The roles of FGF signaling and its effector genes remain poorly characterized in human experimental models mimicking early fetal telencephalic development. We used hiPSC-derived cerebral organoids as an in vitro platform to investigate the effect of FGF8 signaling on neural identity and differentiation. We found that FGF8 treatment increases cellular heterogeneity, leading to distinct telencephalic and mesencephalic-like domains that co-develop in multi-regional organoids. Within telencephalic domains, FGF8 affects the anteroposterior and dorsoventral identity of neural progenitors and the balance between GABAergic and glutamatergic neurons, thus impacting spontaneous neuronal network activity. Moreover, FGF8 efficiently modulates key regulators responsible for several human neurodevelopmental disorders. Overall, our results show that FGF8 signaling is directly involved in both regional patterning and cellular diversity in human cerebral organoids and in modulating genes associated with normal and pathological neural development.
-
- Cell Biology
- Developmental Biology
Chondrocyte columns, which are a hallmark of growth plate architecture, play a central role in bone elongation. Columns are formed by clonal expansion following rotation of the division plane, resulting in a stack of cells oriented parallel to the growth direction. In this work, we analyzed hundreds of Confetti multicolor clones in growth plates of mouse embryos using a pipeline comprising 3D imaging and algorithms for morphometric analysis. Surprisingly, analysis of the elevation angles between neighboring pairs of cells revealed that most cells did not display the typical stacking pattern associated with column formation, implying incomplete rotation of the division plane. Morphological analysis revealed that although embryonic clones were elongated, they formed clusters oriented perpendicular to the growth direction. Analysis of growth plates of postnatal mice revealed both complex columns, composed of ordered and disordered cell stacks, and small, disorganized clusters located in the outer edges. Finally, correlation between the temporal dynamics of the ratios between clusters and columns and between bone elongation and expansion suggests that clusters may promote expansion, whereas columns support elongation. Overall, our findings support the idea that modulations of division plane rotation of proliferating chondrocytes determines the formation of either clusters or columns, a multifunctional design that regulates morphogenesis throughout pre- and postnatal bone growth. Broadly, this work provides a new understanding of the cellular mechanisms underlying growth plate activity and bone elongation during development.