Dedifferentiated early postnatal lung myofibroblasts redifferentiate in adult disease

Alveolarization ensures sufficient lung surface area for gas exchange, and during bulk alveolarization in mice (postnatal day [P] 4.5-14.5), alpha-smooth muscle actin (SMA)+ myofibroblasts accumulate, secrete elastin, and lay down alveolar septum. Herein, we delineate the dynamics of the lineage of early postnatal SMA+ myofibroblasts during and after bulk alveolarization and in response to lung injury. SMA+ lung myofibroblasts first appear at ∼ P2.5 and proliferate robustly. Lineage tracing shows that, at P14.5 and over the next few days, the vast majority of SMA+ myofibroblasts downregulate smooth muscle cell markers and undergo apoptosis. Of note, ∼8% of these dedifferentiated cells and another ∼1% of SMA+ myofibroblasts persist to adulthood. Single cell RNA sequencing analysis of the persistent SMA- cells and SMA+ myofibroblasts in the adult lung reveals distinct gene expression profiles. For instance, dedifferentiated SMA- cells exhibit higher levels of tissue remodeling genes. Most interestingly, these dedifferentiated early postnatal myofibroblasts re-express SMA upon exposure of the adult lung to hypoxia or the pro-fibrotic drug bleomycin. However, unlike during alveolarization, these cells that re-express SMA do not proliferate with hypoxia. In sum, dedifferentiated early postnatal myofibroblasts are a previously undescribed cell type in the adult lung and redifferentiate in response to injury.

Dendritic Cell - Fibroblast Crosstalk via TLR9 and AHR Signaling Drives Lung Fibrogenesis

Idiopathic pulmonary fibrosis (IPF) is characterized by progressive scarring and loss of lung function. With limited treatment options, patients succumb to the disease within 2-5 years. The molecular pathogenesis of IPF regarding the immunologic changes that occur is poorly understood. We characterize a role for non-canonical aryl-hydrocarbon receptor signaling (ncAHR) in dendritic cells (DCs) that leads to production of IL-6 and IL-17, promoting fibrosis. TLR9 signaling in myofibroblasts is shown to regulate production of TDO2 which converts tryptophan into the endogenous AHR ligand kynurenine. Mice with augmented ncAHR signaling were created by crossing floxed AHR exon-2 deletion mice (AHR Δex2 ) with mice harboring a CD11c-Cre. Bleomycin was used to study fibrotic pathogenesis. Isolated CD11c+ cells and primary fibroblasts were treated ex-vivo with relevant TLR agonists and AHR modulating compounds to study how AHR signaling influenced inflammatory cytokine production. Human datasets were also interrogated. Inhibition of all AHR signaling rescued fibrosis, however, AHR Δex2 mice treated with bleomycin developed more fibrosis and DCs from these mice were hyperinflammatory and profibrotic upon adoptive transfer. Treatment of fibrotic fibroblasts with TLR9 agonist increased expression of TDO2. Study of human samples corroborate the relevance of these findings in IPF patients. We also, for the first time, identify that AHR exon-2 floxed mice retain capacity for ncAHR signaling.

Single-cell transcriptomic analysis of human pleura reveals stromal heterogeneity and informs in vitro models of mesothelioma

The pleural lining of the thorax regulates local immunity, inflammation and repair. A variety of conditions, both benign and malignant, including pleural mesothelioma, can affect this tissue. A lack of knowledge concerning the mesothelial and stromal cells comprising the pleura has hampered the development of targeted therapies. Here, we present the first comprehensive single-cell transcriptomic atlas of the human parietal pleura and demonstrate its utility in elucidating pleural biology. We confirm the presence of known universal fibroblasts and describe novel, potentially pleural-specific, fibroblast subtypes. We also present transcriptomic characterisation of multiple in vitro models of benign and malignant mesothelial cells, and characterise these through comparison with in vivo transcriptomic data. While bulk pleural transcriptomes have been reported previously, this is the first study to provide resolution at the single-cell level. We expect our pleural cell atlas will prove invaluable to those studying pleural biology and disease. It has already enabled us to shed light on the transdifferentiation of mesothelial cells, allowing us to develop a simple method for prolonging mesothelial cell differentiation in vitro.

Peripheral Blood Single-Cell Sequencing Uncovers Common and Specific Immune Aberrations in Fibrotic Lung Diseases

Rationale and Objectives

The extent and commonality of peripheral blood immune aberrations in fibrotic interstitial lung diseases are not well characterized. In this study, we aimed to identify common and distinct immune aberrations in patients with idiopathic pulmonary fibrosis (IPF) and fibrotic hypersensitivity pneumonitis (FHP) using cutting-edge single-cell profiling technologies.

Methods

Single-cell RNA sequencing was performed on patients and healthy controls' peripheral blood and bronchoalveolar lavage samples using 10X Genomics 5' gene expression and V(D)J profiling. Cell type composition, transcriptional profiles, cellular trajectories and signaling, and T and B cell receptor repertoires were studied. The standard Seurat R pipeline was followed for cell type composition and differential gene expression analyses. Transcription factor activity was imputed using the DoRothEA-VIPER algorithm. Pseudotime analyses were conducted using Monocle3, while RNA velocity analyses were performed with Velocyto, scVelo, and CellRank. Cell-cell connectomics were assessed using the Connectome R package. V(D)J analyses were conducted using CellRanger and Immcantation frameworks. Across all analyses, disease group differences were assessed using the Wilcoxon rank-sum test.

Measurements and Main Results

327,990 cells from 83 samples were profiled. Overall, changes in monocytes were common to IPF and FHP, whereas lymphocytes exhibited disease-specific aberrations. Both diseases displayed enrichment of CCL3 hi /CCL4 hi CD14+ monocytes (p<2.2e-16) and S100A hi CD14+ monocytes (p<2.2e-16) versus controls. Trajectory and RNA velocity analysis suggested that pro-fibrotic macrophages observed in BAL originated from peripheral blood monocytes. Lymphocytes exhibited disease-specific aberrations, with CD8+ GZMK hi T cells and activated B cells primarily enriched in FHP patients. V(D)J analyses revealed unique T and B cell receptor complementarity-determining region 3 (CDR3) amino acid compositions (p<0.05) in FHP and significant IgA enrichment in IPF (p<5.2e-7).

Conclusions

We identified common and disease-specific immune mechanisms in IPF and FHP; S100A hi monocytes and SPP1 hi macrophages are common to IPF and FHP, whereas GMZK hi T lymphocytes and T and B cell receptor repertoires were unique in FHP. Our findings open novel strategies for the diagnosis and treatment of IPF and FHP.

A hybrid machine learning and regression method for cell type deconvolution of spatial barcoding-based transcriptomic data

Spatial barcoding-based transcriptomic (ST) data require cell type deconvolution for cellular-level downstream analysis. Here we present SDePER, a hybrid machine learning and regression method, to deconvolve ST data using reference single-cell RNA sequencing (scRNA-seq) data. SDePER uses a machine learning approach to remove the systematic difference between ST and scRNA-seq data (platform effects) explicitly and efficiently to ensure the linear relationship between ST data and cell type-specific expression profile. It also considers sparsity of cell types per capture spot and across-spots spatial correlation in cell type compositions. Based on the estimated cell type proportions, SDePER imputes cell type compositions and gene expression at unmeasured locations in a tissue map with enhanced resolution. Applications to coarse-grained simulated data and four real datasets showed that SDePER achieved more accurate and robust results than existing methods, suggesting the importance of considering platform effects, sparsity and spatial correlation in cell type deconvolution.

iDESC: identifying differential expression in single-cell RNA sequencing data with multiple subjects

BACKGROUND

Single-cell RNA sequencing (scRNA-seq) technology has enabled assessment of transcriptome-wide changes at single-cell resolution. Due to the heterogeneity in environmental exposure and genetic background across subjects, subject effect contributes to the major source of variation in scRNA-seq data with multiple subjects, which severely confounds cell type specific differential expression (DE) analysis. Moreover, dropout events are prevalent in scRNA-seq data, leading to excessive number of zeroes in the data, which further aggravates the challenge in DE analysis.

RESULTS

We developed iDESC to detect cell type specific DE genes between two groups of subjects in scRNA-seq data. iDESC uses a zero-inflated negative binomial mixed model to consider both subject effect and dropouts. The prevalence of dropout events (dropout rate) was demonstrated to be dependent on gene expression level, which is modeled by pooling information across genes. Subject effect is modeled as a random effect in the log-mean of the negative binomial component. We evaluated and compared the performance of iDESC with eleven existing DE analysis methods. Using simulated data, we demonstrated that iDESC had well-controlled type I error and higher power compared to the existing methods. Applications of those methods with well-controlled type I error to three real scRNA-seq datasets from the same tissue and disease showed that the results of iDESC achieved the best consistency between datasets and the best disease relevance.

CONCLUSIONS

iDESC was able to achieve more accurate and robust DE analysis results by separating subject effect from disease effect with consideration of dropouts to identify DE genes, suggesting the importance of considering subject effect and dropouts in the DE analysis of scRNA-seq data with multiple subjects.

VISTA (PD-1H) Is a Crucial Immune Regulator to Limit Pulmonary Fibrosis

VISTA (V domain immunoglobulin suppressor of T cell activation, also called PD-1H [programmed death-1 homolog]), a novel immune regulator expressed on myeloid and T lymphocyte lineages, is upregulated in mouse and human idiopathic pulmonary fibrosis (IPF). However, the significance of VISTA and its therapeutic potential in regulating IPF has yet to be defined. To determine the role of VISTA and its therapeutic potential in IPF, the expression profile of VISTA was evaluated from human single-cell RNA sequencing data (IPF Cell Atlas). Inflammatory response and lung fibrosis were assessed in bleomycin-induced experimental pulmonary fibrosis models in VISTA-deficient mice compared with wild-type littermates. In addition, these outcomes were evaluated after VISTA agonistic antibody treatment in the wild-type pulmonary fibrosis mice. VISTA expression was increased in lung tissue-infiltrating monocytes of patients with IPF. VISTA was induced in the myeloid population, mainly circulating monocyte-derived macrophages, during bleomycin-induced pulmonary fibrosis. Genetic ablation of VISTA drastically promoted pulmonary fibrosis, and bleomycin-induced fibroblast activation was dependent on the interaction between VISTA-expressing myeloid cells and fibroblasts. Treatment with VISTA agonistic antibody reduced fibrotic phenotypes accompanied by the suppression of lung innate immune and fibrotic mediators. In conclusion, these results suggest that VISTA upregulation in pulmonary fibrosis may be a compensatory mechanism to limit inflammation and fibrosis, and stimulation of VISTA signaling using VISTA agonists effectively limits the fibrotic innate immune landscape and consequent tissue fibrosis. Further studies are warranted to test VISTA as a novel therapeutic target for the IPF treatment.

Emergence of division of labor in tissues through cell interactions and spatial cues

Most cell types in multicellular organisms can perform multiple functions. However, not all functions can be optimally performed simultaneously by the same cells. Functions incompatible at the level of individual cells can be performed at the cell population level, where cells divide labor and specialize in different functions. Division of labor can arise due to instruction by tissue environment or through self-organization. Here, we develop a computational framework to investigate the contribution of these mechanisms to division of labor within a cell-type population. By optimizing collective cellular task performance under trade-offs, we find that distinguishable expression patterns can emerge from cell-cell interactions versus instructive signals. We propose a method to construct ligand-receptor networks between specialist cells and use it to infer division-of-labor mechanisms from single-cell RNA sequencing (RNA-seq) and spatial transcriptomics data of stromal, epithelial, and immune cells. Our framework can be used to characterize the complexity of cell interactions within tissues.

Rational engineering of lung alveolar epithelium

Engineered whole lungs may one day expand therapeutic options for patients with end-stage lung disease. However, the feasibility of ex vivo lung regeneration remains limited by the inability to recapitulate mature, functional alveolar epithelium. Here, we modulate multimodal components of the alveolar epithelial type 2 cell (AEC2) niche in decellularized lung scaffolds in order to guide AEC2 behavior for epithelial regeneration. First, endothelial cells coordinate with fibroblasts, in the presence of soluble growth and maturation factors, to promote alveolar scaffold population with surfactant-secreting AEC2s. Subsequent withdrawal of Wnt and FGF agonism synergizes with tidal-magnitude mechanical strain to induce the differentiation of AEC2s to squamous type 1 AECs (AEC1s) in cultured alveoli, in situ. These results outline a rational strategy to engineer an epithelium of AEC2s and AEC1s contained within epithelial-mesenchymal-endothelial alveolar-like units, and highlight the critical interplay amongst cellular, biochemical, and mechanical niche cues within the reconstituting alveolus.

A Multi-omic Analysis of the Human Lung Reveals Distinct Cell Specific Aging and Senescence Molecular Programs

Age is a major risk factor for lung disease. To understand the mechanisms underlying this association, we characterized the changing cellular, genomic, transcriptional, and epigenetic landscape of lung aging using bulk and single-cell RNAseq (scRNAseq) data. Our analysis revealed age-associated gene networks that reflected hallmarks of aging, including mitochondrial dysfunction, inflammation, and cellular senescence. Cell type deconvolution revealed age-associated changes in the cellular composition of the lung: decreased alveolar epithelial cells and increased fibroblasts and endothelial cells. In the alveolar microenvironment, aging is characterized by decreased AT2B cells and reduced surfactant production, a finding that was validated by scRNAseq and IHC. We showed that a previously reported senescence signature, SenMayo, captures cells expressing canonical senescence markers. SenMayo signature also identified cell-type specific senescence-associated co-expression modules that have distinct molecular functions, including ECM regulation, cell signaling, and damage response pathways. Analysis of somatic mutations showed that burden was highest in lymphocytes and endothelial cells and was associated with high expression of senescence signature. Finally, aging and senescence gene expression modules were associated with differentially methylated regions, with inflammatory markers such as IL1B, IL6R, and TNF being significantly regulated with age. Our findings provide new insights into the mechanisms underlying lung aging and may have implications for the development of interventions to prevent or treat age-related lung diseases.

microRNA-33 deficiency in macrophages enhances autophagy, improves mitochondrial homeostasis, and protects against lung fibrosis

Idiopathic pulmonary fibrosis (IPF) is a progressive and ultimately fatal disease. Recent findings have shown a marked metabolic reprogramming associated with changes in mitochondrial homeostasis and autophagy during pulmonary fibrosis. The microRNA-33 (miR-33) family of microRNAs (miRNAs) encoded within the introns of sterol regulatory element binding protein (SREBP) genes are master regulators of sterol and fatty acid (FA) metabolism. miR-33 controls macrophage immunometabolic response and enhances mitochondrial biogenesis, FA oxidation, and cholesterol efflux. Here, we show that miR-33 levels are increased in bronchoalveolar lavage (BAL) cells isolated from patients with IPF compared with healthy controls. We demonstrate that specific genetic ablation of miR-33 in macrophages protects against bleomycin-induced pulmonary fibrosis. The absence of miR-33 in macrophages improves mitochondrial homeostasis and increases autophagy while decreasing inflammatory response after bleomycin injury. Notably, pharmacological inhibition of miR-33 in macrophages via administration of anti-miR-33 peptide nucleic acids (PNA-33) attenuates fibrosis in different in vivo and ex vivo mice and human models of pulmonary fibrosis. These studies elucidate a major role of miR-33 in macrophages in the regulation of pulmonary fibrosis and uncover a potentially novel therapeutic approach to treat this disease.

Lung Cell Atlases in Health and Disease

The human lung cellular portfolio, traditionally characterized by cellular morphology and individual markers, is highly diverse, with over 40 cell types and a complex branching structure highly adapted for agile airflow and gas exchange. While constant during adulthood, lung cellular content changes in response to exposure, injury, and infection. Some changes are temporary, but others are persistent, leading to structural changes and progressive lung disease. The recent advance of single-cell profiling technologies allows an unprecedented level of detail and scale to cellular measurements, leading to the rise of comprehensive cell atlas styles of reporting. In this review, we chronical the rise of cell atlases and explore their contributions to human lung biology in health and disease.

Single-cell RNA-seq uncovers cellular heterogeneity and provides a signature for paediatric sleep apnoea

BACKGROUND

Obstructive sleep apnoea (OSA) is a highly prevalent disease and a major cause of systemic inflammation leading to neurocognitive, behavioural, metabolic and cardiovascular dysfunction in children and adults. However, the impact of OSA on the heterogeneity of circulating immune cells remains to be determined.

METHODS

We applied single-cell transcriptomics analysis (scRNA-seq) to identify OSA-induced changes in transcriptional landscape in peripheral blood mononuclear cell (PBMC) composition, which uncovered severity-dependent differences in several cell lineages. Furthermore, a machine-learning approach was used to combine scRNAs-seq cell-specific markers with those differentially expressed in OSA.

RESULTS

scRNA-seq demonstrated OSA-induced heterogeneity in cellular composition and enabled the identification of previously undescribed cell types in PBMCs. We identified a molecular signature consisting of 32 genes, which distinguished OSA patients from various controls with high precision (area under the curve 0.96) and accuracy (93% positive predictive value and 95% negative predictive value) in an independent PBMC bulk RNA expression dataset.

CONCLUSION

OSA deregulates systemic immune function and displays a molecular signature that can be assessed in standard cellular RNA without the need for pre-analytical cell separation, thereby making the assay amenable to application in a molecular diagnostic setting.

CD38 Mediates Lung Fibrosis by Promoting Alveolar Epithelial Cell Aging

Rationale: A prevailing paradigm recognizes idiopathic pulmonary fibrosis (IPF) originating from various alveolar epithelial cell (AEC) injuries, and there is a growing appreciation of AEC aging as a key driver of the pathogenesis. Despite this progress, it is incompletely understood what main factor(s) contribute to the worsened alveolar epithelial aging in lung fibrosis. It remains a challenge how to dampen AEC aging and thereby mitigate the disease progression. Objectives: To determine the role of AEC CD38 (cluster of differentiation 38) in promoting cellular aging and lung fibrosis. Methods: We used single-cell RNA sequencing, real-time PCR, flow cytometry, and Western blotting. Measurements and Main Results: We discovered a pivotal role of CD38, a cardinal nicotinamide adenine dinucleotide (NAD) hydrolase, in AEC aging and its promotion of lung fibrosis. We found increased CD38 expression in IPF lungs that inversely correlated with the lung functions of patients. CD38 was primarily located in the AECs of human lung parenchyma and was markedly induced in IPF AECs. Similarly, CD38 expression was elevated in the AECs of fibrotic lungs of young mice and further augmented in those of old mice, which was in accordance with a worsened AEC aging phenotype and an aggravated lung fibrosis in the old animals. Mechanistically, we found that CD38 elevation downregulated intracellular NAD, which likely led to the aging promoting impairment of the NAD-dependent cellular and molecular activities. Furthermore, we demonstrated that genetic and pharmacological inactivation of CD38 improved these NAD dependent events and ameliorated bleomycin-induced lung fibrosis. Conclusions: Our study suggests targeting alveolar CD38 as a novel and effective therapeutic strategy to treat this pathology.

Microenvironmental sensing by fibroblasts controls macrophage population size

Animal tissues comprise diverse cell types. However, the mechanisms controlling the number of each cell type within tissue compartments remain poorly understood. Here, we report that different cell types utilize distinct strategies to control population numbers. Proliferation of fibroblasts, stromal cells important for tissue integrity, is limited by space availability. In contrast, proliferation of macrophages, innate immune cells involved in defense, repair, and homeostasis, is constrained by growth factor availability. Examination of density-dependent gene expression in fibroblasts revealed that Hippo and TGF-β target genes are both regulated by cell density. We found YAP1, the transcriptional coactivator of the Hippo signaling pathway, directly regulates expression of Csf1, the lineage-specific growth factor for macrophages, through an enhancer of Csf1 that is specifically active in fibroblasts. Activation of YAP1 in fibroblasts elevates Csf1 expression and is sufficient to increase the number of macrophages at steady state. Our data also suggest that expression programs in fibroblasts that change with density may result from sensing of mechanical force through actin-dependent mechanisms. Altogether, we demonstrate that two different modes of population control are connected and coordinated to regulate cell numbers of distinct cell types. Sensing of the tissue environment may serve as a general strategy to control tissue composition.

Computation and visualization of cell-cell signaling topologies in single-cell systems data using Connectome

Single-cell RNA-sequencing data has revolutionized our ability to understand of the patterns of cell-cell and ligand-receptor connectivity that influence the function of tissues and organs. However, the quantification and visualization of these patterns in a way that informs tissue biology are major computational and epistemological challenges. Here, we present Connectome, a software package for R which facilitates rapid calculation and interactive exploration of cell-cell signaling network topologies contained in single-cell RNA-sequencing data. Connectome can be used with any reference set of known ligand-receptor mechanisms. It has built-in functionality to facilitate differential and comparative connectomics, in which signaling networks are compared between tissue systems. Connectome focuses on computational and graphical tools designed to analyze and explore cell-cell connectivity patterns across disparate single-cell datasets and reveal biologic insight. We present approaches to quantify focused network topologies and discuss some of the biologic theory leading to their design.

Characterization of the COPD alveolar niche using single-cell RNA sequencing

Chronic obstructive pulmonary disease (COPD) is a leading cause of death worldwide, however our understanding of cell specific mechanisms underlying COPD pathobiology remains incomplete. Here, we analyze single-cell RNA sequencing profiles of explanted lung tissue from subjects with advanced COPD or control lungs, and we validate findings using single-cell RNA sequencing of lungs from mice exposed to 10 months of cigarette smoke, RNA sequencing of isolated human alveolar epithelial cells, functional in vitro models, and in situ hybridization and immunostaining of human lung tissue samples. We identify a subpopulation of alveolar epithelial type II cells with transcriptional evidence for aberrant cellular metabolism and reduced cellular stress tolerance in COPD. Using transcriptomic network analyses, we predict capillary endothelial cells are inflamed in COPD, particularly through increased CXCL-motif chemokine signaling. Finally, we detect a high-metallothionein expressing macrophage subpopulation enriched in advanced COPD. Collectively, these findings highlight cell-specific mechanisms involved in the pathobiology of advanced COPD.

A Pulmonary Vascular Model From Endothelialized Whole Organ Scaffolds

The development of an in vitro system for the study of lung vascular disease is critical to understanding human pathologies. Conventional culture systems fail to fully recapitulate native microenvironmental conditions and are typically limited in their ability to represent human pathophysiology for the study of disease and drug mechanisms. Whole organ decellularization provides a means to developing a construct that recapitulates structural, mechanical, and biological features of a complete vascular structure. Here, we developed a culture protocol to improve endothelial cell coverage in whole lung scaffolds and used single-cell RNA-sequencing analysis to explore the impact of decellularized whole lung scaffolds on endothelial phenotypes and functions in a biomimetic bioreactor system. Intriguingly, we found that the phenotype and functional signals of primary pulmonary microvascular revert back—at least partially—toward native lung endothelium. Additionally, human induced pluripotent stem cell-derived endothelium cultured in decellularized lung systems start to gain various native human endothelial phenotypes. Vascular barrier function was partially restored, while small capillaries remained patent in endothelial cell-repopulated lungs. To evaluate the ability of the engineered endothelium to modulate permeability in response to exogenous stimuli, lipopolysaccharide (LPS) was introduced into repopulated lungs to simulate acute lung injury. After LPS treatment, proinflammatory signals were significantly increased and the vascular barrier was impaired. Taken together, these results demonstrate a novel platform that recapitulates some pulmonary microvascular functions and phenotypes at a whole organ level. This development may help pave the way for using the whole organ engineering approach to model vascular diseases.

Transcriptomics of bronchoalveolar lavage cells identifies new molecular endotypes of sarcoidosis

BACKGROUND

Sarcoidosis is a multisystem granulomatous disease of unknown origin with a variable and often unpredictable course and pattern of organ involvement. In this study we sought to identify specific bronchoalveolar lavage (BAL) cell gene expression patterns indicative of distinct disease phenotypic traits.

METHODS

RNA sequencing by Ion Torrent Proton was performed on BAL cells obtained from 215 well-characterised patients with pulmonary sarcoidosis enrolled in the multicentre Genomic Research in Alpha-1 Antitrypsin Deficiency and Sarcoidosis (GRADS) study. Weighted gene co-expression network analysis and nonparametric statistics were used to analyse genome-wide BAL transcriptome. Validation of results was performed using a microarray expression dataset of an independent sarcoidosis cohort (Freiburg, Germany; n=50).

RESULTS

Our supervised analysis found associations between distinct transcriptional programmes and major pulmonary phenotypic manifestations of sarcoidosis including T-helper type 1 (Th1) and Th17 pathways associated with hilar lymphadenopathy, transforming growth factor-β1 (TGFB1) and mechanistic target of rapamycin (MTOR) signalling with parenchymal involvement, and interleukin (IL)-7 and IL-2 with airway involvement. Our unsupervised analysis revealed gene modules that uncovered four potential sarcoidosis endotypes including hilar lymphadenopathy with increased acute T-cell immune response; extraocular organ involvement with PI3K activation pathways; chronic and multiorgan disease with increased immune response pathways; and multiorgan involvement, with increased IL-1 and IL-18 immune and inflammatory responses. We validated the occurrence of these endotypes using gene expression, pulmonary function tests and cell differentials from Freiburg.

CONCLUSION

Taken together, our results identify BAL gene expression programmes that characterise major pulmonary sarcoidosis phenotypes and suggest the presence of distinct disease molecular endotypes.

Integrated Single-Cell Atlas of Endothelial Cells of the Human Lung

BACKGROUND

Cellular diversity of the lung endothelium has not been systematically characterized in humans. We provide a reference atlas of human lung endothelial cells (ECs) to facilitate a better understanding of the phenotypic diversity and composition of cells comprising the lung endothelium.

METHODS

We reprocessed human control single-cell RNA sequencing (scRNAseq) data from 6 datasets. EC populations were characterized through iterative clustering with subsequent differential expression analysis. Marker genes were validated by fluorescent microscopy and in situ hybridization. scRNAseq of primary lung ECs cultured in vitro was performed. The signaling network between different lung cell types was studied. For cross-species analysis or disease relevance, we applied the same methods to scRNAseq data obtained from mouse lungs or from human lungs with pulmonary hypertension.

RESULTS

Six lung scRNAseq datasets were reanalyzed and annotated to identify >15 000 vascular EC cells from 73 individuals. Differential expression analysis of EC revealed signatures corresponding to endothelial lineage, including panendothelial, panvascular, and subpopulation-specific marker gene sets. Beyond the broad cellular categories of lymphatic, capillary, arterial, and venous ECs, we found previously indistinguishable subpopulations; among venous EC, we identified 2 previously indistinguishable populations: pulmonary-venous ECs (COL15A1neg) localized to the lung parenchyma and systemic-venous ECs (COL15A1pos) localized to the airways and the visceral pleura; among capillary ECs, we confirmed their subclassification into recently discovered aerocytes characterized by EDNRB, SOSTDC1, and TBX2 and general capillary EC. We confirmed that all 6 endothelial cell types, including the systemic-venous ECs and aerocytes, are present in mice and identified endothelial marker genes conserved in humans and mice. Ligand-receptor connectome analysis revealed important homeostatic crosstalk of EC with other lung resident cell types. scRNAseq of commercially available primary lung ECs demonstrated a loss of their native lung phenotype in culture. scRNAseq revealed that endothelial diversity is maintained in pulmonary hypertension. Our article is accompanied by an online data mining tool (www.LungEndothelialCellAtlas.com).

CONCLUSIONS

Our integrated analysis provides a comprehensive and well-crafted reference atlas of ECs in the normal lung and confirms and describes in detail previously unrecognized endothelial populations across a large number of humans and mice.

Chronic lung diseases are associated with gene expression programs favoring SARS-CoV-2 entry and severity

Patients with chronic lung disease (CLD) have an increased risk for severe coronavirus disease-19 (COVID-19) and poor outcomes. Here, we analyze the transcriptomes of 611,398 single cells isolated from healthy and CLD lungs to identify molecular characteristics of lung cells that may account for worse COVID-19 outcomes in patients with chronic lung diseases. We observe a similar cellular distribution and relative expression of SARS-CoV-2 entry factors in control and CLD lungs. CLD AT2 cells express higher levels of genes linked directly to the efficiency of viral replication and the innate immune response. Additionally, we identify basal differences in inflammatory gene expression programs that highlight how CLD alters the inflammatory microenvironment encountered upon viral exposure to the peripheral lung. Our study indicates that CLD is accompanied by changes in cell-type-specific gene expression programs that prime the lung epithelium for and influence the innate and adaptive immune responses to SARS-CoV-2 infection.

Fibroblasts positive for meflin have anti-fibrotic property in pulmonary fibrosis

The prognosis of elderly individuals with idiopathic pulmonary fibrosis (IPF) remains poor. Fibroblastic foci, in which aggregates of proliferating fibroblasts and myofibroblasts are involved, are the pathological hallmark lesions in IPF to represent focal areas of active fibrogenesis. Fibroblast heterogeneity in fibrotic lesions hampers the discovery of the pathogenesis of pulmonary fibrosis. Therefore, to determine of the pathogenesis of IPF, identification of functional fibroblasts is warranted. This study was aimed to determine the role of fibroblasts positive for meflin, identified as a potential marker for mesenchymal stromal cells, during the development of pulmonary fibrosis. We characterised meflin-positive cells in a single cell atlas established by single-cell RNA sequencing (scRNA-seq)-based profiling of 243 472 cells from 32 IPF lungs and 29 normal lung samples. scRNA-seq combined with in situ RNA hybridisation identified proliferating fibroblasts positive for meflin in fibroblastic foci, not dense fibrosis, of fibrotic lungs in IPF patients. We determined the role of fibroblasts positive for meflin using bleomycin (BLM)-induced pulmonary fibrosis. A BLM-induced lung fibrosis model for meflin-deficient mice showed that fibroblasts positive for meflin had anti-fibrotic property to prevent pulmonary fibrosis. Although transforming growth factor-β-induced fibrogenesis and cell senescence with senescence-associated secretory phenotype were exacerbated in fibroblasts via the repression or lack of meflin, these were inhibited in meflin-deficient fibroblasts with meflin reconstitution. These findings provide evidence to show the biological importance of meflin expression on fibroblasts and myofibroblasts in the active fibrotic region of pulmonary fibrosis.

Single-cell meta-analysis of SARS-CoV-2 entry genes across tissues and demographics

Angiotensin-converting enzyme 2 (ACE2) and accessory proteases (TMPRSS2 and CTSL) are needed for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) cellular entry, and their expression may shed light on viral tropism and impact across the body. We assessed the cell-type-specific expression of ACE2, TMPRSS2 and CTSL across 107 single-cell RNA-sequencing studies from different tissues. ACE2, TMPRSS2 and CTSL are coexpressed in specific subsets of respiratory epithelial cells in the nasal passages, airways and alveoli, and in cells from other organs associated with coronavirus disease 2019 (COVID-19) transmission or pathology. We performed a meta-analysis of 31 lung single-cell RNA-sequencing studies with 1,320,896 cells from 377 nasal, airway and lung parenchyma samples from 228 individuals. This revealed cell-type-specific associations of age, sex and smoking with expression levels of ACE2, TMPRSS2 and CTSL. Expression of entry factors increased with age and in males, including in airway secretory cells and alveolar type 2 cells. Expression programs shared by ACE2+TMPRSS2+ cells in nasal, lung and gut tissues included genes that may mediate viral entry, key immune functions and epithelial-macrophage cross-talk, such as genes involved in the interleukin-6, interleukin-1, tumor necrosis factor and complement pathways. Cell-type-specific expression patterns may contribute to the pathogenesis of COVID-19, and our work highlights putative molecular pathways for therapeutic intervention.

Chronic lung diseases are associated with gene expression programs favoring SARS-CoV-2 entry and severity

Patients with chronic lung disease (CLD) have an increased risk for severe coronavirus disease-19 (COVID-19) and poor outcomes. Here, we analyzed the transcriptomes of 605,904 single cells isolated from healthy and CLD lungs to identify molecular characteristics of lung cells that may account for worse COVID-19 outcomes in patients with chronic lung diseases. We observed a similar cellular distribution and relative expression of SARS-CoV-2 entry factors in control and CLD lungs. CLD epithelial cells expressed higher levels of genes linked directly to the efficiency of viral replication and innate immune response. Additionally, we identified basal differences in inflammatory gene expression programs that highlight how CLD alters the inflammatory microenvironment encountered upon viral exposure to the peripheral lung. Our study indicates that CLD is accompanied by changes in cell-type-specific gene expression programs that prime the lung epithelium for and influence the innate and adaptive immune responses to SARS-CoV-2 infection.

Single-Cell Transcriptional Archetypes of Airway Inflammation in Cystic Fibrosis

Rationale: Cystic fibrosis (CF) is a life-shortening, multisystem hereditary disease caused by abnormal chloride transport. CF lung disease is driven by innate immune dysfunction and exaggerated inflammatory responses that contribute to tissue injury. To define the transcriptional profile of this airway immune dysfunction, we performed the first single-cell transcriptome characterization of CF sputum.Objectives: To define the transcriptional profile of sputum cells and its implication in the pathogenesis of immune function and the development of CF lung disease.Methods: We performed single-cell RNA sequencing of sputum cells from nine subjects with CF and five healthy control subjects. We applied novel computational approaches to define expression-based cell function and maturity profiles, herein called transcriptional archetypes.Measurements and Main Results: The airway immune cell repertoire shifted from alveolar macrophages in healthy control subjects to a predominance of recruited monocytes and neutrophils in CF. Recruited lung mononuclear phagocytes were abundant in CF and were separated into the following three archetypes: activated monocytes, monocyte-derived macrophages, and heat shock-activated monocytes. Neutrophils were the most prevalent in CF, with a dominant immature proinflammatory archetype. Although CF monocytes exhibited proinflammatory features, both monocytes and neutrophils showed transcriptional evidence of abnormal phagocytic and cell-survival programs.Conclusions: Our findings offer an opportunity to understand subject-specific immune dysfunction and its contribution to divergent clinical courses in CF. As we progress toward personalized applications of therapeutic and genomic developments, we hope this inflammation-profiling approach will enable further discoveries that change the natural history of CF lung disease.

Single-cell RNA-seq reveals ectopic and aberrant lung-resident cell populations in idiopathic pulmonary fibrosis

We provide a single-cell atlas of idiopathic pulmonary fibrosis (IPF), a fatal interstitial lung disease, by profiling 312,928 cells from 32 IPF, 28 smoker and nonsmoker controls, and 18 chronic obstructive pulmonary disease (COPD) lungs. Among epithelial cells enriched in IPF, we identify a previously unidentified population of aberrant basaloid cells that coexpress basal epithelial, mesenchymal, senescence, and developmental markers and are located at the edge of myofibroblast foci in the IPF lung. Among vascular endothelial cells, we identify an ectopically expanded cell population transcriptomically identical to bronchial restricted vascular endothelial cells in IPF. We confirm the presence of both populations by immunohistochemistry and independent datasets. Among stromal cells, we identify IPF myofibroblasts and invasive fibroblasts with partially overlapping cells in control and COPD lungs. Last, we confirm previous findings of profibrotic macrophage populations in the IPF lung. Our comprehensive catalog reveals the complexity and diversity of aberrant cellular populations in IPF.

Collagen-producing lung cell atlas identifies multiple subsets with distinct localization and relevance to fibrosis

Collagen-producing cells maintain the complex architecture of the lung and drive pathologic scarring in pulmonary fibrosis. Here we perform single-cell RNA-sequencing to identify all collagen-producing cells in normal and fibrotic lungs. We characterize multiple collagen-producing subpopulations with distinct anatomical localizations in different compartments of murine lungs. One subpopulation, characterized by expression of Cthrc1 (collagen triple helix repeat containing 1), emerges in fibrotic lungs and expresses the highest levels of collagens. Single-cell RNA-sequencing of human lungs, including those from idiopathic pulmonary fibrosis and scleroderma patients, demonstrate similar heterogeneity and CTHRC1-expressing fibroblasts present uniquely in fibrotic lungs. Immunostaining and in situ hybridization show that these cells are concentrated within fibroblastic foci. We purify collagen-producing subpopulations and find disease-relevant phenotypes of Cthrc1-expressing fibroblasts in in vitro and adoptive transfer experiments. Our atlas of collagen-producing cells provides a roadmap for studying the roles of these unique populations in homeostasis and pathologic fibrosis.

Platform Effects on Regeneration by Pulmonary Basal Cells as Evaluated by Single-Cell RNA Sequencing

Cell-based therapies have shown promise for treating myriad chronic pulmonary diseases through direct application of epithelial progenitors or by way of engineered tissue grafts or whole organs. To elucidate environmental effects on epithelial regenerative outcomes in vitro, here, we isolate and culture a population of pharmacologically expanded basal cells (peBCs) from rat tracheas. At peak basal marker expression, we simultaneously split peBCs into four in vitro platforms: organoid, air-liquid interface (ALI), engineered trachea, and engineered lung. Following differentiation, these samples are evaluated using single-cell RNA sequencing (scRNA-seq) and computational pipelines are developed to compare samples both globally and at the population level. A sample of native rat tracheal epithelium is also evaluated by scRNA-seq as a control for engineered epithelium. Overall, this work identifies platform-specific effects that support the use of engineered models to achieve the most physiologic differential outcomes in pulmonary epithelial regenerative applications.

Clinical management of women presenting with field effect of HPV and intraepithelial disease

Anogenital human papillomavirus is the most common sexually transmitted infection acquired through skin-to-skin contact. Most infections are cleared by an intact immune system. Persistence of these infections results in precancerous lesions and, eventually, to cancers of the cervix, vagina, vulva, and perianal area. The introduction of the prophylactic human papilloma virus (HPV) vaccinations may reduce the incidence of these infections, but the effect of these vaccinations will be seen only in the decades that follow. In the meantime, multiple therapies such as immune modulators, ablative modalities, and surgical excision are used in an attempt to treat precancerous lesions and hence prevent cancer. There is an increase in multicentric disease in young women, especially with the HIV epidemic and in women who are immune compromised. This article aims to address the challenges and management options in women who have a field effect of HPV-associated disease.

Age- and diapause-related acid and alkaline phosphatase activities in the intestine and malpighian tubules of the Colorado potato beetle, Leptinotarsa decemlineata (Say)

Specific activities for soluble (s) and membrane (m)-bound acid (ACP) and alkaline phosphatases (ALP) were determined in the midgut, hindgut, and Malpighian tubules for developing, prediapausing, and diapausing adult Colorado potato beetles, Leptinotarsa decemlineata (Say). High ACP activities were found in the hindgut and Malpighian tubules while high ALP activities were found in the Malpighian tubules. Variation in both ACP and ALP activities in each tissue reflects fluctuation in protein synthesis and secretion involved with digestion, excretion, and other unknown functions. Phosphatase activities in the tissues examined show the dynamic nature of diapause in this insect. Diapausing beetles showed increases in phosphatase activity after hormone treatments. JHA treatments increased s-ACP and m-ACP activities in all tissues but 20-HE did not increase activity in any tissue. Allatotropin tended to mimic the effects of JHA treatment. The s-ALP activity was also increased in all tissues whereas m-ALP was increased in the midgut and hindgut by JHA treatment. Malpighian tubule m-ALP activity was only increased by 20-HE treatments. Allatotropin was not as effective in increasing ALP activities as it was with ACP activities.

Effect of pyriproxyfen and photoperiod on free amino acid concentrations and proteins in the hemolymph of the Colorado potato beetle, Leptinotarsa decemlineata (Say)

A North Dakota strain of the Colorado potato beetle, Leptinotarsa decemlineata (Say), was reared under both short- (8L:16D) and long-day (17L:7D) conditions. Age-related and pyriproxyfen- (JHA-) induced changes in hemolymph free amino acids and proteins were examined. Under a short-day photoperiod, the total free amino acid concentration in the hemolymph increased gradually up to 20 days of adult life, but the long-day beetles showed marked increases during the first 10 days and then decreased afterwards. Proline, glutamine and valine were the most abundant free amino acids in both sexes of beetles held under either short- or long-day photoregims. JHA treatment of diapausing adults, held under either short- or long-day conditions after treatment, terminated diapause as indicated by re-emergence from the vermiculite, feeding, mating, changes in free amino acid levels, the disappearance of diapause protein 1 and appearance of vitellogenin in the hemolymph. Furthermore, most of the JHA-treated females held under long-day conditions also matured oocytes and oviposited, but those held under short-day conditions did not.

Ecdysteroidostatin from the house fly, Musca domestica

Ovaries from house flies maintained on sucrose secrete large amounts of ecdysteroid when they are cultured with ovarian ecdysteroidogenic hormone, OEH. However, ovarian ecdysteroid secretion is reduced by incubation with both OEH and the ovarian ecdysteroidostatin (OES). A partially purified OES fraction from a semi-preparative reverse phase HPLC C18 column caused a 98% inhibition of ovarian ecdysteroid secretion in vitro at a concentration of 0.8 equivalents per microliter. Ovaries can be activated to produce ecdysteroid in vivo by feeding diet containing protein to flies maintained on sucrose. Ecdysteroid secretion was inhibited when the in vivo stimulated ovaries were cultured with OES. This suggests that OES does not interfere with the OEH activation mechanism, but blocks ovarian ecdysteroid synthesis or release. Furthermore, OES inhibition is reversible and ecdysteroid secretion resumes when OES is removed. Musca OES could explain the decrease in ecdysteroid levels found in flies after mid-vitellogenesis. Both adult male and female abdomens contain OES, but OES was not transferred to females during mating. Evidence is presented that OES is not a trypsin modulating oostatic factor.

Phenytoin reduces the contraction of recessive dystrophic epidermolysis bullosa fibroblast populated collagen gels

Recessive dystrophic epidermolysis bullosa (RDEB) is a group of genetic disorders in which blistering occurs below the basement membrane, in many cases resulting in extensive scar formation, contractures and mitten deformities. Our aim was to compare quantitatively the contraction forces generated by normal and RDEB fibroblats and to investigate the effect of Phenytoin (5,5-diphenyl-2,4-imidazolidinedione, sodium salt; PHT). PHT is an anticonvulsant agent, that causes fibrosis as a side effect. This study utilised conventional untethered fibroblast populated collagen lattice contraction and a quantitative force measurement instrument, the culture force monitor (CFM). The RDEB cell lines were hypercontractile, generating 2.5 times the force of normal fibroblasts, though they appeared morphologically normal. In untethered collagen gels PHT (20 micrograms/ml) significantly reduced contraction of both normal and RDEB fibroblasts over 7 days. Pre-treatment of RDEB cells for 5 days also produced a 40% reduction in contraction as measured in the CFM. One suggested mechanism of PHT action is through inhibition of matrix metalloproteinase activity, but the similar effects of PHT and Colchicine (an inhibitor of microtubule polymerisation) in the CFM, indicate that it may act on contraction through disruption of microfilaments and changes to cell shape. These findings show that isolated RDEB fibroblasts retain the hypercontractile features of many of the patient's lesion sites and imply that local application of PHT may have a therapeutic effect in controlling contraction.

Mortyn Jones Memorial Lecture. Limbic regions mediating central actions of oxytocin on the milk-ejection reflex in the rat

Central oxytocin administration has a profound facilitatory effect on the patterning of the milk-ejection reflex in the lactating rat. Lesion and microinjection studies indicate that this action is, in part, mediated via a population of limbic neurones in the bed nuclei of the stria terminalis and ventrolateral septum, which have been shown to possess oxytocin receptors and to be activated by selective oxytocin-receptor agonists in vitro. In vivo electrophysiological recordings reveal that some of these neurones display cyclical activity which is highly correlated to each milk ejection, and are rapidly activated following i.c.v. administration of oxytocin, coincident with the facilitation of milk ejection activity. A hypothetical model is proposed in which this population of limbic neurones serves to gate the activity of a pacemaker which, in turn, coordinates the bursting of hypothalamic magnocellular neurones. The oxytocin innervation of these neurones and their expression of oxytocin receptors increases in the postpartum period, and the resultant enhanced sensitivity leads to a greater facilitatory response during lactation. Inhibitory opioid and noradrenergic inputs which converge on these oxytocin-sensitive neurones may function to switch off the facilitatory circuit during periods of stress. Thus, this population of limbic neurones participates in the regulation of neuroendocrine activity during lactation by providing an appropriate degree of feedback to alter the patterning of the milk-ejection reflex.