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.

Multiple Plasticity Regulators Reveal Targets Specifying an Induced Predatory Form in Nematodes.

The ability to translate a single genome into multiple phenotypes, or developmental plasticity, defines how phenotype derives from more than just genes. However, to study the evolutionary targets of plasticity and their evolutionary fates, we need to understand how genetic regulators of plasticity control downstream gene expression. Here, we have identified a transcriptional response specific to polyphenism (i.e., discrete plasticity) in the nematode Pristionchus pacificus. This species produces alternative resource-use morphs-microbivorous and predatory forms, differing in the form of their teeth, a morphological novelty-as influenced by resource availability. Transcriptional profiles common to multiple polyphenism-controlling genes in P. pacificus reveal a suite of environmentally sensitive loci, or ultimate target genes, that make up an induced developmental response. Additionally, in vitro assays show that one polyphenism regulator, the nuclear receptor NHR-40, physically binds to promoters with putative HNF4α (the nuclear receptor class including NHR-40) binding sites, suggesting this receptor may directly regulate genes that describe alternative morphs. Among differentially expressed genes were morph-limited genes, highlighting factors with putative "on-off" function in plasticity regulation. Further, predatory morph-biased genes included candidates-namely, all four P. pacificus homologs of Hsp70, which have HNF4α motifs-whose natural variation in expression matches phenotypic differences among P. pacificus wild isolates. In summary, our study links polyphenism regulatory loci to the transcription producing alternative forms of a morphological novelty. Consequently, our findings establish a platform for determining how specific regulators of morph-biased genes may influence selection on plastic phenotypes.

A fern AINTEGUMENTA gene mirrors BABY BOOM in promoting apogamy in Ceratopteris richardii.

Asexual reproduction is widespread in land plants, including ferns where 10% of all species are obligate asexuals. In these ferns, apogamous sporophytes are generated directly from gametophytes, bypassing fertilization. In the model fern Ceratopteris richardii, a sexual species, apogamy can be induced by culture on high sugar media. BABY BOOM (BBM) genes in angiosperms are known to promote somatic embryogenesis, which like apogamy produce sporophytes without fertilization. Here, a Brassica napus BBM (BnBBM) was used to investigate genetic similarity between apogamy in ferns and somatic embryogenesis in angiosperms. A C. richardii transcriptome was constructed from which one AINTEGUMENTA-LIKE unigene, CrANT, was identified. Whole mount in situ hybridization showed that CrANT is expressed in sperm and fertilized eggs. Phylogenetic analysis grouped CrANT with other non-seed-plant ANT genes to the euANT clade but in a branch separate from BBM genes. Overexpression of CrANT or BnBBM promotes apogamy in C. richardii without sugar supplement. CrANT knockdown gametophytes responded weakly to sugar for apogamy promotion. Theses results suggest some genetic conservation between apogamy and somatic embryogenesis and that such asexual reproduction may be ancient.

A Ceratopteris EXCESS MICROSPOROCYTES1 suppresses reproductive transition in the fern vegetative leaves.

Land plant sexual reproduction involves the transition of cells from somatic to reproductive identity during post-embryonic development. In Arabidopsis, the leucine-rich repeat receptor-like kinase EXCESS MICROSPOROCYTES1 (EXS/EMS1) restricts the number of sporogenous cells during the transition from diploid tissue to haploid spore production by promoting the formation of the tapetum cell layer within the anther. Although all land plants studied contain EMS1 genes, its function is unknown beyond a few angiosperms. In the model fern Ceratopteris (Ceratopteris richardii), we discovered an EMS1 homolog (CrEMS1) that functions to suppress formation of reproductive structures on vegetative leaves of the fern sporophyte, a role not found in angiosperms. Suppression of CrEMS1 by RNAi did not affect sporogenesis on reproductive leaves but did affect antheridium production of the fern gametophyte. Expression patterns of CrEMS1 across developmental stages suggest threshold levels of CrEMS1 control the specification of reproductive organs during both generations of the fern. Additional EMS1 homologs present in the fern genome suggest a dynamic role of EMS1 receptors in the evolution of reproductive development in vascular plants.

An integrated cell atlas of the lung in health and disease.

Single-cell technologies have transformed our understanding of human tissues. Yet, studies typically capture only a limited number of donors and disagree on cell type definitions. Integrating many single-cell datasets can address these limitations of individual studies and capture the variability present in the population. Here we present the integrated Human Lung Cell Atlas (HLCA), combining 49 datasets of the human respiratory system into a single atlas spanning over 2.4 million cells from 486 individuals. The HLCA presents a consensus cell type re-annotation with matching marker genes, including annotations of rare and previously undescribed cell types. Leveraging the number and diversity of individuals in the HLCA, we identify gene modules that are associated with demographic covariates such as age, sex and body mass index, as well as gene modules changing expression along the proximal-to-distal axis of the bronchial tree. Mapping new data to the HLCA enables rapid data annotation and interpretation. Using the HLCA as a reference for the study of disease, we identify shared cell states across multiple lung diseases, including SPP1+ profibrotic monocyte-derived macrophages in COVID-19, pulmonary fibrosis and lung carcinoma. Overall, the HLCA serves as an example for the development and use of large-scale, cross-dataset organ atlases within the Human Cell Atlas.

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.

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 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.

Single-cell RNA sequencing reveals profibrotic roles of distinct epithelial and mesenchymal lineages in pulmonary fibrosis.

Pulmonary fibrosis (PF) is a form of chronic lung disease characterized by pathologic epithelial remodeling and accumulation of extracellular matrix (ECM). To comprehensively define the cell types, mechanisms, and mediators driving fibrotic remodeling in lungs with PF, we performed single-cell RNA sequencing of single-cell suspensions from 10 nonfibrotic control and 20 PF lungs. Analysis of 114,396 cells identified 31 distinct cell subsets/states. We report that a remarkable shift in epithelial cell phenotypes occurs in the peripheral lung in PF and identify several previously unrecognized epithelial cell phenotypes, including a KRT5- /KRT17 + pathologic, ECM-producing epithelial cell population that was highly enriched in PF lungs. Multiple fibroblast subtypes were observed to contribute to ECM expansion in a spatially discrete manner. Together, these data provide high-resolution insights into the complexity and plasticity of the distal lung epithelium in human disease and indicate a diversity of epithelial and mesenchymal cells contribute to pathologic lung fibrosis.

A sulfotransferase dosage-dependently regulates mouthpart polyphenism in the nematode Pristionchus pacificus.

Polyphenism, the extreme form of developmental plasticity, is the ability of a genotype to produce discrete morphologies matched to alternative environments. Because polyphenism is likely to be under switch-like molecular control, a comparative genetic approach could reveal the molecular targets of plasticity evolution. Here we report that the lineage-specific sulfotransferase SEUD-1, which responds to environmental cues, dosage-dependently regulates polyphenism of mouthparts in the nematode Pristionchus pacificus. SEUD-1 is expressed in cells producing dimorphic morphologies, thereby integrating an intercellular signalling mechanism at its ultimate target. Additionally, multiple alterations of seud-1 support it as a potential target for plasticity evolution. First, a recent duplication of seud-1 in a sister species reveals a direct correlation between genomic dosage and polyphenism threshold. Second, inbreeding to produce divergent polyphenism thresholds resulted in changes in transcriptional dosage of seud-1. Our study thus offers a genetic explanation for how plastic responses evolve.

Author Correction: A sulfotransferase dosage-dependently regulates mouthpart polyphenism in the nematode Pristionchus pacificus.

The original version of this Article contained errors in Figure 4. In panel a, the x axis labels of bars 6-11 were incorrect, as depicted in the associated Author Correction. These errors have been corrected in both the PDF and HTML versions of the Article.

Transient and stable transformation of Ceratopteris richardii gametophytes.

BACKGROUND: Ferns, being vascular yet seedless, present unparalleled opportunities to investigate important questions regarding the evolution and development of land plants. Ceratopteris richardii, a diploid, homosporous fern has been advanced as a model fern system; however, the tenuous ability to transform the genome of this fern greatly limited its usefulness as a model organism. Here we report a simple and reliable Agrobacterium-mediated method for generating transient and stable transformants of mature C. richardii gametophytes. RESULTS: Transformation success was achieved by enzyme treatment that partially digested the cell walls of mature gametophytes to facilitate Agrobacteria infection. Co-incubation of Agrobacteria with enzymatically treated gametophytes was sufficient to generate transient transformants at a frequency of nearly 90% under optimal conditions. Stable transformation was achieved at a rate of nearly 3% by regenerating entire gametophytes from single transformed cells from T0 gametophytes on selective media. CONCLUSIONS: This transformation method will allow for the immediate observation of phenotypes in the haploid gametophytes of transformed plants, as well as the generation of stably transformed C. richardii lines for further analysis. Transformation capability will greatly facilitate gene functional studies in C. richardii, more fully realizing the potential of this model fern species. These protocols may be adapted to other plant species that are recalcitrant to Agrobacterium-mediated transformation.