Cell Population-resolved Multi-Omics Atlas of the Developing Lung.

The lung is a vital organ that undergoes extensive morphological and functional changes during postnatal development. To disambiguate how different cell populations contribute to organ development, we performed proteomic and transcriptomic analyses of four sorted cell populations from the lung of human subjects aged 0 to 8 years-old with a focus on early life. The cell populations analyzed included epithelial, endothelial, mesenchymal, and immune cells. Our results revealed distinct molecular signatures for each of the sorted cell populations that enable the description of molecular shifts occurring in these populations during post-natal development. We confirmed that the proteome of the different cell populations was distinct regardless of age and identified functions specific to each population. We identified a series of cell population protein markers, including those located at the cell surface, that show differential expression and distribution on RNA in situ hybridization and immunofluorescence imaging. We validated the spatial distribution of AT1 and endothelial cell surface markers. Temporal analyses of the proteomes of the four populations revealed processes modulated during postnatal development and clarified the findings obtained from whole tissue proteome studies. Finally, the proteome was compared to a transcriptomics survey performed on the same lung samples to evaluate processes under post-transcriptional control.

Single-cell resolution of human airway epithelial cells exposed to bronchiolitis obliterans-associated chemicals.

Bronchiolitis obliterans (BO) is a fibrotic lung disease characterized by progressive luminal narrowing and obliteration of the small airways. In the nontransplant population, inhalation exposure to certain chemicals is associated with BO; however, the mechanisms contributing to disease induction remain poorly understood. This study's objective was to use single-cell RNA sequencing for the identification of transcriptomic signatures common to primary human airway epithelial cells after chemical exposure to BO-associated chemicals-diacetyl or nitrogen mustard-to help explain BO induction. Primary airway epithelial cells were cultured at air-liquid interface and exposed to diacetyl, nitrogen mustard, or control vapors. Cultures were dissociated and sequenced for single-cell RNA. Differential gene expression and functional pathway analyses were compared across exposures. In total, 75,663 single cells were captured and sequenced from all exposure conditions. Unbiased clustering identified 11 discrete phenotypes, including 5 basal, 2 ciliated, and 2 secretory cell clusters. With chemical exposure, the proportion of cells assigned to keratin 5+ basal cells decreased, whereas the proportion of cells aligned to secretory cell clusters increased compared with control exposures. Functional pathway analysis identified interferon signaling and antigen processing/presentation as pathways commonly upregulated after diacetyl or nitrogen mustard exposure in a ciliated cell cluster. Conversely, the response of airway basal cells differed significantly with upregulation of the unfolded protein response in diacetyl-exposed basal cells, not seen in nitrogen mustard-exposed cultures. These new insights provide early identification of airway epithelial signatures common to BO-associated chemical exposures.NEW & NOTEWORTHY Bronchiolitis obliterans (BO) is a devastating fibrotic lung disease of the small airways, or bronchioles. This original manuscript uses single-cell RNA sequencing for identifying common signatures of chemically exposed airway epithelial cells in BO induction. Chemical exposure reduced the proportion of keratin 5+ basal cells while increasing the proportion of keratin 4+ suprabasal cells. Functional pathways contributory to these shifts differed significantly across exposures. These new results highlight similarities and differences in BO induction across exposures.

Bulk RNA sequencing of human pediatric lung cell populations reveals unique transcriptomic signature associated with postnatal pulmonary development.

Postnatal lung development results in an increasingly functional organ prepared for gas exchange and pathogenic challenges. It is achieved through cellular differentiation and migration. Changes in the tissue architecture during this development process are well-documented and increasing cellular diversity associated with it are reported in recent years. Despite recent progress, transcriptomic and molecular pathways associated with human postnatal lung development are yet to be fully understood. In this study, we investigated gene expression patterns associated with healthy pediatric lung development in four major enriched cell populations (epithelial, endothelial, and nonendothelial mesenchymal cells, along with lung leukocytes) from 1-day-old to 8-yr-old organ donors with no known lung disease. For analysis, we considered the donors in four age groups [less than 30 days old neonates, 30 days to < 1 yr old infants, toddlers (1 to < 2 yr), and children 2 yr and older] and assessed differentially expressed genes (DEG). We found increasing age-associated transcriptional changes in all four major cell types in pediatric lung. Transition from neonate to infant stage showed highest number of DEG compared with the number of DEG found during infant to toddler- or toddler to older children-transitions. Profiles of differential gene expression and further pathway enrichment analyses indicate functional epithelial cell maturation and increased capability of antigen presentation and chemokine-mediated communication. Our study provides a comprehensive reference of gene expression patterns during healthy pediatric lung development that will be useful in identifying and understanding aberrant gene expression patterns associated with early life respiratory diseases.NEW & NOTEWORTHY This study presents postnatal transcriptomic changes in major cell populations in human lung, namely endothelial, epithelial, mesenchymal cells, and leukocytes. Although human postnatal lung development continues through early adulthood, our results demonstrate that greatest transcriptional changes occur in first few months of life during neonate to infant transition. These early transcriptional changes in lung parenchyma are particularly notable for functional maturation and activation of alveolar type II cell genes.

CD4+ T Cell Tolerance to Tissue-Restricted Self Antigens Is Mediated by Antigen-Specific Regulatory T Cells Rather Than Deletion.

Deletion of self-antigen-specific T cells during thymic development provides protection from autoimmunity. However, it is unclear how efficiently this occurs for tissue-restricted self antigens, or how immune tolerance is maintained for self-antigen-specific T cells that routinely escape deletion. Here we show that endogenous CD4+ T cells with specificity for a set of tissue-restricted self antigens were not deleted at all. For pancreatic self antigen, this resulted in an absence of steady-state tolerance, while for the lung and intestine, tolerance was maintained by the enhanced presence of thymically-derived antigen-specific Foxp3+ regulatory T (Treg) cells. Unlike deletional tolerance, Treg cell-mediated tolerance was broken by successive antigen challenges. These findings reveal that for some tissue-restricted self antigens, tolerance relies entirely on nondeletional mechanisms that are less durable than T cell deletion. This might explain why autoimmunity is often tissue-specific, and it offers a rationale for cancer vaccine strategies targeting tissue-restricted tumor antigens.

Gene Expression Risk Scores for COVID-19 Illness Severity.

BACKGROUND: The correlates of coronavirus disease 2019 (COVID-19) illness severity following infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are incompletely understood. METHODS: We assessed peripheral blood gene expression in 53 adults with confirmed SARS-CoV-2 infection clinically adjudicated as having mild, moderate, or severe disease. Supervised principal components analysis was used to build a weighted gene expression risk score (WGERS) to discriminate between severe and nonsevere COVID-19. RESULTS: Gene expression patterns in participants with mild and moderate illness were similar, but significantly different from severe illness. When comparing severe versus nonsevere illness, we identified >4000 genes differentially expressed (false discovery rate < 0.05). Biological pathways increased in severe COVID-19 were associated with platelet activation and coagulation, and those significantly decreased with T-cell signaling and differentiation. A WGERS based on 18 genes distinguished severe illness in our training cohort (cross-validated receiver operating characteristic-area under the curve [ROC-AUC] = 0.98), and need for intensive care in an independent cohort (ROC-AUC = 0.85). Dichotomizing the WGERS yielded 100% sensitivity and 85% specificity for classifying severe illness in our training cohort, and 84% sensitivity and 74% specificity for defining the need for intensive care in the validation cohort. CONCLUSIONS: These data suggest that gene expression classifiers may provide clinical utility as predictors of COVID-19 illness severity.

New insights into the natural history of bronchopulmonary dysplasia from proteomics and multiplexed immunohistochemistry.

Bronchopulmonary dysplasia (BPD) is a disease of prematurity related to the arrest of normal lung development. The objective of this study was to better understand how proteome modulation and cell-type shifts are noted in BPD pathology. Pediatric human donors aged 1-3 yr were classified based on history of prematurity and histopathology consistent with "healed" BPD (hBPD, n = 3) and "established" BPD (eBPD, n = 3) compared with respective full-term born (n = 6) age-matched term controls. Proteins were quantified by tandem mass spectroscopy with selected Western blot validations. Multiplexed immunofluorescence (MxIF) microscopy was performed on lung sections to enumerate cell types. Protein abundances and MxIF cell frequencies were compared among groups using ANOVA. Cell type and ontology enrichment were performed using an in-house tool and/or EnrichR. Proteomics detected 5,746 unique proteins, 186 upregulated and 534 downregulated, in eBPD versus control with fewer proteins differentially abundant in hBPD as compared with age-matched term controls. Cell-type enrichment suggested a loss of alveolar type I, alveolar type II, endothelial/capillary, and lymphatics, and an increase in smooth muscle and fibroblasts consistent with MxIF. Histochemistry and Western analysis also supported predictions of upregulated ferroptosis in eBPD versus control. Finally, several extracellular matrix components mapping to angiogenesis signaling pathways were altered in eBPD. Despite clear parsing by protein abundance, comparative MxIF analysis confirms phenotypic variability in BPD. This work provides the first demonstration of tandem mass spectrometry and multiplexed molecular analysis of human lung tissue for critical elucidation of BPD trajectory-defining factors into early childhood.NEW & NOTEWORTHY We provide new insights into the natural history of bronchopulmonary dysplasia in donor human lungs after the neonatal intensive care unit hospitalization. This study provides new insights into how the proteome and histopathology of BPD changes in early childhood, uncovering novel pathways for future study.

Proteomic Analysis of Human Lung Development.

Rationale: The current understanding of human lung development derives mostly from animal studies. Although transcript-level studies have analyzed human donor tissue to identify genes expressed during normal human lung development, protein-level analysis that would enable the generation of new hypotheses on the processes involved in pulmonary development are lacking. Objectives: To define the temporal dynamic of protein expression during human lung development. Methods: We performed proteomics analysis of human lungs at 10 distinct times from birth to 8 years to identify the molecular networks mediating postnatal lung maturation. Measurements and Main Results: We identified 8,938 proteins providing a comprehensive view of the developing human lung proteome. The analysis of the data supports the existence of distinct molecular substages of alveolar development and predicted the age of independent human lung samples, and extensive remodeling of the lung proteome occurred during postnatal development. Evidence of post-transcriptional control was identified in early postnatal development. An extensive extracellular matrix remodeling was supported by changes in the proteome during alveologenesis. The concept of maturation of the immune system as an inherent part of normal lung development was substantiated by flow cytometry and transcriptomics. Conclusions: This study provides the first in-depth characterization of the human lung proteome during development, providing a unique proteomic resource freely accessible at Lungmap.net. The data support the extensive remodeling of the lung proteome during development, the existence of molecular substages of alveologenesis, and evidence of post-transcriptional control in early postnatal development.

CX3CR1 Engagement by Respiratory Syncytial Virus Leads to Induction of Nucleolin and Dysregulation of Cilia-related Genes.

Respiratory syncytial virus (RSV) contains a conserved CX3C motif on the ectodomain of the G-protein. The motif has been indicated as facilitating attachment of the virus to the host initiating infection via the human CX3CR1 receptor. The natural CX3CR1 ligand, CX3CL1, has been shown to induce signaling pathways resulting in transcriptional changes in the host cells. We hypothesize that binding of RSV to CX3CR1 via CX3C leads to transcriptional changes in host epithelial cells. Using transcriptomic analysis, the effect of CX3CR1 engagement by RSV was investigated. Normal human bronchial epithelial (NHBE) cells were infected with RSV virus containing either wildtype G-protein, or a mutant virus containing a CX4C mutation in the G-protein. RNA sequencing was performed on mock and 4-days-post-infected cultures. NHBE cultures were also treated with purified recombinant wild-type A2 G-protein. Here we report that RSV infection resulted in significant changes in the levels 766 transcripts. Many nuclear associated proteins were upregulated in the WT group, including nucleolin. Alternatively, cilia-associated genes, including CC2D2A and CFAP221 (PCDP1), were downregulated. The addition of recombinant G-protein to the culture lead to the suppression of cilia-related genes while also inducing nucleolin. Mutation of the CX3C motif (CX4C) reversed these effects on transcription decreasing nucleolin induction and lessening the suppression of cilia-related transcripts in culture. Furthermore, immunohistochemical staining demonstrated decreases in in ciliated cells and altered morphology. Therefore, it appears that engagement of CX3CR1 leads to induction of genes necessary for RSV entry as well as dysregulation of genes associated with cilia function.ImportanceRespiratory Syncytial Virus (RSV) has an enormous impact on infants and the elderly including increased fatality rates and potential for causing lifelong lung problems. Humans become infected with RSV through the inhalation of viral particles exhaled from an infected individual. These virus particles contain specific proteins that the virus uses to attach to human ciliated lung epithelial cells, initiating infection. Two viral proteins, G-protein and F-protein, have been shown to bind to human CX3CR1and nucleolin, respectively. Here we show that the G-protein induces nucleolin and suppresses gene transcripts specific to ciliated cells. Furthermore, we show that mutation of the CX3C-motif on the G-protein, CX4C, reverses these transcriptional changes.

A systems genomics approach uncovers molecular associates of RSV severity.

Respiratory syncytial virus (RSV) infection results in millions of hospitalizations and thousands of deaths each year. Variations in the adaptive and innate immune response appear to be associated with RSV severity. To investigate the host response to RSV infection in infants, we performed a systems-level study of RSV pathophysiology, incorporating high-throughput measurements of the peripheral innate and adaptive immune systems and the airway epithelium and microbiota. We implemented a novel multi-omic data integration method based on multilayered principal component analysis, penalized regression, and feature weight back-propagation, which enabled us to identify cellular pathways associated with RSV severity. In both airway and immune cells, we found an association between RSV severity and activation of pathways controlling Th17 and acute phase response signaling, as well as inhibition of B cell receptor signaling. Dysregulation of both the humoral and mucosal response to RSV may play a critical role in determining illness severity.

Prenatal histological, cellular, and molecular anomalies in trisomy 21 lung.

Down syndrome (DS), also known as trisomy 21 (T21), is the most common human chromosomal anomaly. Although DS can affect many organ systems, lung and heart disease are the leading causes of death. An abundance of existing data suggests that lung abnormalities originate postnatally in DS. However, a single report of branching insufficiency in DS has inferred a potential prenatal origin. The histology of T21 fetal lungs (n = 15) was assessed by an experienced pathologist. Spatial differences in cellular phenotypes were examined using immunohistochemistry (IHC). Comprehensive gene expression in prenatal T21 lungs (n = 19), and age-matched controls (n = 19), was performed using high-throughput RNA sequencing (RNAseq) and validated by RT-qPCR. Histopathological abnormalities were observed in approximately half of T21 prenatal lung samples analyzed, which included dilated terminal airways/acinar tubules, dilated lymphatics, and arterial wall thickening. IHC for Ki67 revealed significant reductions in epithelial and mesenchymal cell proliferation, predominantly in tissues displaying pathology. IHC demonstrated that airway smooth muscle was reduced and discontinuous in the proximal airway in conjunction with reduced SOX2. RNAseq identified 118 genes significantly dysregulated (FDR < 0.05) in T21 lung when unadjusted and 316 genes when adjusted for age. Ontology analysis showed that IFN pathway genes were appreciably upregulated, whereas complement and coagulation cascades and extracellular matrix pathway genes were downregulated. RT-qPCR confirmed the changes in genes associated with these pathways in prenatal T21 lungs. Our data demonstrate that specific histological, cellular, and molecular abnormalities occur prenatally in different compartments of human T21 lung, which could be representative of premature stage progression. © 2021 The Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Temporal Dysbiosis of Infant Nasal Microbiota Relative to Respiratory Syncytial Virus Infection.

BACKGROUND: Respiratory syncytial virus (RSV) is a leading cause of infant respiratory disease. Infant airway microbiota has been associated with respiratory disease risk and severity. The extent to which interactions between RSV and microbiota occur in the airway, and their impact on respiratory disease susceptibility and severity, are unknown. METHODS: We carried out 16S rRNA microbiota profiling of infants in the first year of life from (1) a cross-sectional cohort of 89 RSV-infected infants sampled during illness and 102 matched healthy controls, and (2) a matched longitudinal cohort of 12 infants who developed RSV infection and 12 who did not, sampled before, during, and after infection. RESULTS: We identified 12 taxa significantly associated with RSV infection. All 12 taxa were differentially abundant during infection, with 8 associated with disease severity. Nasal microbiota composition was more discriminative of healthy vs infected than of disease severity. CONCLUSIONS: Our findings elucidate the chronology of nasal microbiota dysbiosis and suggest an altered developmental trajectory associated with RSV infection. Microbial temporal dynamics reveal indicators of disease risk, correlates of illness and severity, and impact of RSV infection on microbiota composition.

Airway Gene Expression Correlates of Respiratory Syncytial Virus Disease Severity and Microbiome Composition in Infants.

BACKGROUND: Respiratory syncytial virus (RSV) is the leading cause of severe respiratory disease in infants. The causes and correlates of severe illness in the majority of infants are poorly defined. METHODS: We recruited a cohort of RSV-infected infants and simultaneously assayed the molecular status of their airways and the presence of airway microbiota. We used rigorous statistical approaches to identify gene expression patterns associated with disease severity and microbiota composition, separately and in combination. RESULTS: We measured comprehensive airway gene expression patterns in 106 infants with primary RSV infection. We identified an airway gene expression signature of severe illness dominated by excessive chemokine expression. We also found an association between Haemophilus influenzae, disease severity, and airway lymphocyte accumulation. Exploring the time of onset of clinical symptoms revealed acute activation of interferon signaling following RSV infection in infants with mild or moderate illness, which was absent in subjects with severe illness. CONCLUSIONS: Our data reveal that airway gene expression patterns distinguish mild/moderate from severe illness. Furthermore, our data identify biomarkers that may be therapeutic targets or useful for measuring efficacy of intervention responses.

Lung disease manifestations in Down syndrome.

Down syndrome (DS) is one of the most prevalent chromosomal abnormalities worldwide, affecting 1 in 700 live births. Although multiple organ systems are affected by the chromosomal defects, respiratory failure and lung disease are the leading causes of morbidity and mortality observed in DS. Manifestations of DS in the respiratory system encompass the entire lung starting from the nasopharynx to the trachea/upper airways to the lower airways and alveolar spaces, as well as vascular and lymphatic defects. Most of our knowledge on respiratory illness in persons with DS arises from pediatric studies; however, many of these disorders present early in infancy, supporting developmental mechanisms. In this review, we will focus on the different lung phenotypes in DS, as well as the genetic and molecular pathways that may be contributing to these complications during development.

Rare Pulmonary Connective Tissue Type Mast Cells Regulate Lung Endothelial Cell Angiogenesis.

Within the human lung, mast cells typically reside adjacent to the conducting airway and assume a mucosal phenotype (MCT). In rare pathologic conditions, connective tissue phenotype mast cells (MCTCs) can be found in the lung parenchyma. MCTCs accumulate in the lungs of infants with severe bronchopulmonary dysplasia, a chronic lung disease associated with preterm birth, which is characterized by pulmonary vascular dysmorphia. The human mast cell line (LUVA) was used to model MCTCs or MCTs. The ability of MCTCs to affect vascular organization during fetal lung development was tested in mouse lung explant cultures. The effect of MCTCs on in vitro tube formation and barrier function was studied using primary fetal human pulmonary microvascular endothelial cells. The mechanistic role of MCTC proteases was tested using inhibitors. MCTCLUVA but not MCTLUVA was associated with vascular dysmorphia in lung explants. In vitro, the addition of MCTCLUVA potentiated fetal human pulmonary microvascular endothelial cell interactions, inhibited tube stability, and disrupted endothelial cell junctions. Protease inhibitors ameliorated the ability of MCTCLUVA to alter endothelial cell angiogenic activities in vitro and ex vivo. These data indicate that MCTCs may directly contribute to disrupted angiogenesis in bronchopulmonary dysplasia. A better understanding of factors that regulate mast cell subtype and their different effector functions is essential.

Collagen VI Deficiency Results in Structural Abnormalities in the Mouse Lung.

Collagen VI (COL6) is known for its role in a spectrum of congenital muscular dystrophies, which are often accompanied by respiratory dysfunction. However, little is known regarding the function of COL6 in the lung. We confirmed the presence of COL6 throughout the basement membrane region of mouse lung tissue. Lung structure and organization were studied in a previously described Col6a1-/- mouse, which does not produce detectable COL6 in the lung. The Col6a1-/- mouse displayed histopathologic alveolar and airway abnormalities. The airspaces of Col6a1-/- lungs appeared simplified, with larger (29%; P < 0.01) and fewer (31%; P < 0.001) alveoli. These airspace abnormalities included reduced isolectin B4+ alveolar capillaries and surfactant protein C-positive alveolar epithelial type-II cells. Alterations in lung function consistent with these histopathologic changes were evident. Col6a1-/- mice also displayed multiple airway changes, including increased branching (59%; P < 0.001), increased mucosal thickness (34%; P < 0.001), and increased epithelial cell density (13%; P < 0.001). Comprehensive transcriptome analysis revealed that the loss of COL6 is associated with reductions in integrin-paxillin-phosphatidylinositol 3-kinase signaling in vivo. In vitro, COL6 promoted steady-state phosphorylated paxillin levels and reduced cell density (16% to 28%; P < 0.05) at confluence. Inhibition of phosphatidylinositol 3-kinase, or its downstream effectors, resulted in increased cell density to a level similar to that seen on matrices lacking COL6.

Repetitive diacetyl vapor exposure promotes ubiquitin proteasome stress and precedes bronchiolitis obliterans pathology.

Bronchiolitis obliterans (BO) is a devastating lung disease seen commonly after lung transplant, following severe respiratory tract infection or chemical inhalation exposure. Diacetyl (DA; 2,3-butanedione) is a highly reactive alpha-diketone known to cause BO when inhaled, however, the mechanisms of how inhalation exposure leads to BO development remains poorly understood. In the current work, we combined two clinically relevant models for studying the pathogenesis of DA-induced BO: (1) an in vivo rat model of repetitive DA vapor exposures with recovery and (2) an in vitro model of primary human airway epithelial cells exposed to pure DA vapors. Rats exposed to 5 consecutive days 200 parts-per-million DA 6 h per day had worsening survival, persistent hypoxemia, poor weight gain, and histologic evidence of BO 14 days after DA exposure cessation. At the end of exposure, increased expression of the ubiquitin stress protein ubiquitin-C accumulated within DA-exposed rat lung homogenates and localized primarily to the airway epithelium, the primary site of BO development. Lung proteasome activity increased concurrently with ubiquitin-C expression after DA exposure, supportive of significant proteasome stress. In primary human airway cultures, global proteomics identified 519 significantly modified proteins in DA-exposed samples relative to controls with common pathways of the ubiquitin proteasome system, endosomal reticulum transport, and response to unfolded protein pathways being upregulated and cell-cell adhesion and oxidation-reduction pathways being downregulated. Collectively, these two models suggest that diacetyl inhalation exposure causes abundant protein damage and subsequent ubiquitin proteasome stress prior to the development of chemical-induced BO pathology.

Implementing sequence-based antigenic distance calculation into immunological shape space model.

BACKGROUND: In 2009, a novel influenza vaccine was distributed worldwide to combat the H1N1 influenza "swine flu" pandemic. However, antibodies induced by the vaccine display differences in their specificity and cross-reactivity dependent on pre-existing immunity. Here, we present a computational model that can capture the effect of pre-existing immunity on influenza vaccine responses. The model predicts the region of the virus hemagglutinin (HA) protein targeted by antibodies after vaccination as well as the level of cross-reactivity induced by the vaccine. We tested our model by simulating a scenario similar to the 2009 pandemic vaccine and compared the results to antibody binding data obtained from human subjects vaccinated with the monovalent 2009 H1N1 influenza vaccine. RESULTS: We found that both specificity and cross-reactivity of the antibodies induced by the 2009 H1N1 influenza HA protein were affected by the viral strain the individual was originally exposed. Specifically, the level of antigenic relatedness between the original exposure HA antigen and the 2009 HA protein affected antigenic-site immunodominance. Moreover, antibody cross-reactivity was increased when the individual's pre-existing immunity was specific to an HA protein antigenically distinct from the 2009 pandemic strain. Comparison of simulation data with antibody binding data from human serum samples demonstrated qualitative and quantitative similarities between the model and real-life immune responses to the 2009 vaccine. CONCLUSION: We provide a novel method to evaluate expected outcomes in antibody specificity and cross-reactivity after influenza vaccination in individuals with different influenza HA antigen exposure histories. The model produced similar outcomes as what has been previously reported in humans after receiving the 2009 influenza pandemic vaccine. Our results suggest that differences in cross-reactivity after influenza vaccination should be expected in individuals with different exposure histories.

Transcriptional characterisation of human lung cells identifies novel mesenchymal lineage markers.

RATIONALE: The lung mesenchyme gives rise to multiple distinct lineages of cells in the mature respiratory system, including smooth muscle cells of the airway and vasculature. However, a thorough understanding of the specification and mesenchymal cell diversity in the human lung is lacking. METHODS: We completed single-cell RNA sequencing analysis of fetal human lung tissues. Canonical correlation analysis, clustering, cluster marker gene identification and t-distributed stochastic neighbour embedding representation was performed in Seurat. Cell populations were annotated using ToppFun. Immunohistochemistry and in situ hybridisation were used to validate spatiotemporal gene expression patterns for key marker genes. RESULTS: We identified molecularly distinct populations representing "committed" fetal human lung endothelial cells, pericytes and smooth muscle cells. Early endothelial lineages expressed "classic" endothelial cell markers (platelet endothelial cell adhesion molecule/CD31 and claudin 5), while pericytes expressed platelet-derived growth factor receptor-ß, Thy-1 membrane glycoprotein and basement membrane molecules (collagen IV, laminin and proteoglycans). We observed a large population of "nonspecific" human lung mesenchymal progenitor cells characterised by expression of collagen I and multiple elastin fibre genes (ELN, MFAP2 and FBN1). We closely characterised the diversity of mesenchymal lineages defined by α2-smooth muscle actin (ACTA2) expression. Two cell populations, with the highest levels of ACTA2 transcriptional activity, expressed unique sets of markers associated with airway or vascular smooth muscle cells. Spatiotemporal analysis of these marker genes confirmed early and persistent spatial specification of airway (HHIP, MYLK and IGF1) and vascular (NTRK3 and MEF2C) smooth muscle cells in the developing human lung. CONCLUSION: Our data suggest that specification of distinct airway and vascular smooth muscle cell phenotypes is established early in development and can be identified using the markers we provide.

CX3CR1 as a respiratory syncytial virus receptor in pediatric human lung.

BACKGROUND: Data on the host factors that contribute to infection of young children by respiratory syncytial virus (RSV) are limited. The human chemokine receptor, CX3CR1, has recently been implicated as an RSV receptor. Here we evaluate a role for CX3CR1 in pediatric lung RSV infections. METHODS: CX3CR1 transcript levels in the upper and lower pediatric airways were assessed. Tissue localization and cell-specific expression was confirmed using in situ hybridization and immunohistochemistry. The role of CX3CR1 in RSV infection was also investigated using a novel physiological model of pediatric epithelial cells. RESULTS: Low levels of CX3CR1 transcript were often, but not always, expressed in both upper (62%) and lower airways (36%) of pediatric subjects. CX3CR1 transcript and protein expression was detected in epithelial cells of normal human pediatric lung tissues. CX3CR1 expression was readily detected on primary cultures of differentiated pediatric/infant human lung epithelial cells. RSV demonstrated preferential infection of CX3CR1-positive cells, and blocking CX3CR1/RSV interaction significantly decreased viral load. CONCLUSION: CX3CR1 is present in the airways of pediatric subjects where it may serve as a receptor for RSV infection. Furthermore, CX3CR1 appears to play a mechanistic role in mediating viral infection of pediatric airway epithelial cells in vitro.

A novel in vitro model of primary human pediatric lung epithelial cells.

BACKGROUND: Current in vitro human lung epithelial cell models derived from adult tissues may not accurately represent all attributes that define homeostatic and disease mechanisms relevant to the pediatric lung. METHODS: We report methods for growing and differentiating primary Pediatric Human Lung Epithelial (PHLE) cells from organ donor infant lung tissues. We use immunohistochemistry, flow cytometry, quantitative RT-PCR, and single cell RNA sequencing (scRNAseq) analysis to characterize the cellular and transcriptional heterogeneity of PHLE cells. RESULTS: PHLE cells can be expanded in culture up to passage 6, with a doubling time of ~4 days, and retain attributes of highly enriched epithelial cells. PHLE cells can form resistant monolayers, and undergo differentiation when placed at air-liquid interface. When grown at Air-Liquid Interface (ALI), PHLE cells expressed markers of airway epithelial cell lineages. scRNAseq suggests the cultures contained 4 main sub-phenotypes defined by expression of FOXJ1, KRT5, MUC5B, and SFTPB. These cells are available to the research community through the Developing Lung Molecular Atlas Program Human Tissue Core. CONCLUSION: Our data demonstrate that PHLE cells provide a novel in vitro human cell model that represents the pediatric airway epithelium, which can be used to study perinatal developmental and pediatric disease mechanisms.