Cited by 197

The DOI matches its Crossref record; 11 of the 11 references it cites were found in a registry or an index.

Guided construction of single cell reference for human and mouse lung

Guo, Minzhe, Morley, Michael P., Jiang, Cheng, Wu, Yixin, Li, Guangyuan, Du, Yina, Zhao, Shuyang, Wagner, Andrew, Cakar, Adnan Cihan, Kouril, Michal, Jin, Kang, Gaddis, Nathan, Kitzmiller, Joseph A., Stewart, Kathleen, Basil, Maria C., Lin, Susan M., Ying, Yun, Babu, Apoorva, Wikenheiser-Brokamp, Kathryn A., Mun, Kyu Shik, Naren, Anjaparavanda P., Clair, Geremy, Adkins, Joshua N., Pryhuber, Gloria S., Misra, Ravi S., Aronow, Bruce J., Tickle, Timothy L., Salomonis, Nathan, Sun, Xin, Morrisey, Edward E., Whitsett, Jeffrey A., NHLBI LungMAP Consortium, Lin, Sara, Xu, Yan·2023·Nature Communications

Published version of a preprint

Abstract

Accurate cell type identification is a key and rate-limiting step in single-cell data analysis. Single-cell references with comprehensive cell types, reproducible and functionally validated cell identities, and common nomenclatures are much needed by the research community for automated cell type annotation, data integration, and data sharing. Here, we develop a computational pipeline utilizing the LungMAP CellCards as a dictionary to consolidate single-cell transcriptomic datasets of 104 human lungs and 17 mouse lung samples to construct LungMAP single-cell reference (CellRef) for both normal human and mouse lungs. CellRefs define 48 human and 40 mouse lung cell types catalogued from diverse anatomic locations and developmental time points. We demonstrate the accuracy and stability of LungMAP CellRefs and their utility for automated cell type annotation of both normal and diseased lungs using multiple independent methods and testing data. We develop user-friendly web interfaces for easy access and maximal utilization of the LungMAP CellRefs.

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1038/s41467-023-40173-5

Cites

11 of 11 references verified0 citation errors · 0 not checkable

Statistics

Figures

Fig. 1 Data collection and the guided single-cell reference (CellRef) construction pipeline. A Characteristics of the collection of single cell/nucleus (sc/sn) RNA-seq datasets from normal human lung samples. B Schematic workflow for the LungMAP CellRef construction guided by using LungMAP CellCards as a cell type dictionary.
Fig. 1 Data collection and the guided single-cell reference (CellRef) construction pipeline. A Characteristics of the collection of single cell/nucleus (sc/sn) RNA-seq datasets from normal human lung samples. B Schematic workflow for the LungMAP CellRef construction guided by using LungMAP CellCards as a cell type dictionary.CC BY · PubMed Central
Fig. 2 The construction of LungMAP Human Lung CellRef. A Uniform manifold approximation and projection (UMAP) visualization of seed cells representing 48 lung cell types of normal human lung, termed LungMAP Human Lung CellRef Seed. Cells were colored by their predicted seed identities. B UMAP visualization of the complete single-cell reference for normal human lung, denoted as LungMAP Human Lung CellRef, which contains 347,970 cells from 104 donors and defines 48 cell types in normal human lung. Cells were colored by their predicted identities. C Validation of the seed cell identity using the expression of cell type selective marker genes derived from LungMAP CellCards. D Reconstruction of cell lineage relationships using hierarchical clustering analysis of cell type pseudo-bulk gene expression profiles. Color represents Pearson’s correlation value of pseudo-bulk expression profiles. Labels ending with “.Seed” represent pseudo-bulk profiles created by averaging gene expression in the cells of each cell type in the human lung CellRef Seed, while labels ending with “.CellRef” represent pseudo-bulk profiles created using gene expression of each cell type in the complete human lung CellRef.
Fig. 2 The construction of LungMAP Human Lung CellRef. A Uniform manifold approximation and projection (UMAP) visualization of seed cells representing 48 lung cell types of normal human lung, termed LungMAP Human Lung CellRef Seed. Cells were colored by their predicted seed identities. B UMAP visualization of the complete single-cell reference for normal human lung, denoted as LungMAP Human Lung CellRef, which contains 347,970 cells from 104 donors and defines 48 cell types in normal human lung. Cells were colored by their predicted identities. C Validation of the seed cell identity using the expression of cell type selective marker genes derived from LungMAP CellCards. D Reconstruction of cell lineage relationships using hierarchical clustering analysis of cell type pseudo-bulk gene expression profiles. Color represents Pearson’s correlation value of pseudo-bulk expression profiles. Labels ending with “.Seed” represent pseudo-bulk profiles created by averaging gene expression in the cells of each cell type in the human lung CellRef Seed, while labels ending with “.CellRef” represent pseudo-bulk profiles created using gene expression of each cell type in the complete human lung CellRef.CC BY · PubMed Central
Fig. 3 The construction of LungMAP Mouse Lung Development CellRef. A The developmental time points of mouse lung single-cell transcriptome data used for the guided CellRef construction. B Uniform manifold approximation and projection (UMAP) visualization of the seed cells representing 40 cell types of the developing mouse lung, termed LungMAP Mouse Lung Development CellRef Seed. Cells were colored by predicted seed identities. C UMAP visualization of CellRef for normal mouse lung development, named LungMAP Mouse Lung Development CellRef. Cells were colored by their predicted identities. D Validation of seed cell identities using expression of cell type selective marker genes. E Lineage relationships among mouse lung cell types were reconstructed using hierarchical clustering analysis using pseudo-bulk gene expression profiles. Color represents Pearson’s correlation value of pseudo-bulk expression profiles. Labels ending with “.Seed” represent pseudo-bulk profiles created by averaging gene expression in the cells of each cell type in the mouse lung CellRef Seed, while labels ending with “.CellRef” represent pseudo-bulk profiles created using gene expression of each cell type in the complete mouse lung CellRef.
Fig. 3 The construction of LungMAP Mouse Lung Development CellRef. A The developmental time points of mouse lung single-cell transcriptome data used for the guided CellRef construction. B Uniform manifold approximation and projection (UMAP) visualization of the seed cells representing 40 cell types of the developing mouse lung, termed LungMAP Mouse Lung Development CellRef Seed. Cells were colored by predicted seed identities. C UMAP visualization of CellRef for normal mouse lung development, named LungMAP Mouse Lung Development CellRef. Cells were colored by their predicted identities. D Validation of seed cell identities using expression of cell type selective marker genes. E Lineage relationships among mouse lung cell types were reconstructed using hierarchical clustering analysis using pseudo-bulk gene expression profiles. Color represents Pearson’s correlation value of pseudo-bulk expression profiles. Labels ending with “.Seed” represent pseudo-bulk profiles created by averaging gene expression in the cells of each cell type in the mouse lung CellRef Seed, while labels ending with “.CellRef” represent pseudo-bulk profiles created using gene expression of each cell type in the complete mouse lung CellRef.CC BY · PubMed Central
Fig. 4 Online interactive exploration of LungMAP CellRef Seed using Lung Gene Expression Analysis (LGEA) web portal. The LungMAP Human Lung CellRef Seed was comprised of 8080 seed cells representing 48 normal lung cell types. A The “Gene Expression Query” interface allows users to input a gene of interest (top) and visualize of the expression of the queried gene in UMAP embeddings of cells (bottom), Colors represent the seed cell identities (bottom left) or the expression of the input gene (bottom right). B Visualization of the gene expression pattern (top: expression distribution; middle: expression frequency and sensitivity; bottom: fold change and p-value of differential expression) across all cell types in the CellRef Seed. Box center lines, bounds of the box, and whiskers indicate medians, first and third quartiles, and minimum and maximum values within 1.5×IQR (interquartile range) of the box limits, respectively. P value for each cell type was determined using a nonparametric binomial test47 for single-cell RNA-seq data by comparing the expression of FOXJ1 in the cell type with its expression in all other cells in the CellRef Seed. See Fig. 4 source data table for number of cells in each cell type. C LGEA hosts comprehensive cell information related with the query cell type. D “Cell Signature Query” function retrieves signature gene expression statistics of a given cell type and bar-plot visualization of signature genes expression across all cell types in the CellRef Seed. P values were determined using a nonparametric binomial test47 for single-cell RNA-seq data by comparing gene expression in the ciliated cells (n = 200 cells) with all other cells (n = 7880 cells) in the CellRef Seed. In (A) and (B), FOXJ1 expression was shown as example. In (C) and (D), Ciliated cells were used as example.
Fig. 4 Online interactive exploration of LungMAP CellRef Seed using Lung Gene Expression Analysis (LGEA) web portal. The LungMAP Human Lung CellRef Seed was comprised of 8080 seed cells representing 48 normal lung cell types. A The “Gene Expression Query” interface allows users to input a gene of interest (top) and visualize of the expression of the queried gene in UMAP embeddings of cells (bottom), Colors represent the seed cell identities (bottom left) or the expression of the input gene (bottom right). B Visualization of the gene expression pattern (top: expression distribution; middle: expression frequency and sensitivity; bottom: fold change and p-value of differential expression) across all cell types in the CellRef Seed. Box center lines, bounds of the box, and whiskers indicate medians, first and third quartiles, and minimum and maximum values within 1.5×IQR (interquartile range) of the box limits, respectively. P value for each cell type was determined using a nonparametric binomial test47 for single-cell RNA-seq data by comparing the expression of FOXJ1 in the cell type with its expression in all other cells in the CellRef Seed. See Fig. 4 source data table for number of cells in each cell type. C LGEA hosts comprehensive cell information related with the query cell type. D “Cell Signature Query” function retrieves signature gene expression statistics of a given cell type and bar-plot visualization of signature genes expression across all cell types in the CellRef Seed. P values were determined using a nonparametric binomial test47 for single-cell RNA-seq data by comparing gene expression in the ciliated cells (n = 200 cells) with all other cells (n = 7880 cells) in the CellRef Seed. In (A) and (B), FOXJ1 expression was shown as example. In (C) and (D), Ciliated cells were used as example.CC BY · PubMed Central
Fig. 5 Cell type annotation and evaluation using the LungMAP Human Lung CellRef. A Schematic workflow of the automated cell type annotation and evaluation pipeline. B Distributions of cell type prediction scores in each test data. Prediction scores using CellRef Seed (yellow bars) are comparable to those using the complete CellRef (blue bars). Prediction scores (between 0 and 1) were calculated by the Seurat v4 MapQuery function for each cell. Box center lines, bounds of the box, and whiskers indicate medians, first and third quartiles, and minimum and maximum values within 1.5×IQR (interquartile range) of the box limits, respectively. GSM5388411: 6228 cells, GSM5388412: 8329 cells, GSM5388413: 7143 cells, GSM4504966: 8381 cells, GSM4504967: 8043 cells, GSM4035472: 5767 cells. C Consistency of cell type predictions using the CellRef Seed and CellRef in each test data. Consistency percentages (y axis) were calculated for cells in each test dataset (color) passing different thresholds of prediction scores (x axis). D–H Evaluation of automated cell type annotations for three of our test data (GSM5388411/12/13, three scRNA-seq of normal human lungs). Evaluation of the other three test data samples were shown in Supplementary Fig. 7. Basal and suprabasal cells were combined in prediction. D UMAP visualization of cells with prediction scores ≥ default cut-off (mean-1 standard deviation) and predicted annotations with at least 5 cells. Cells were colored by automated cell type annotations using the CellRef Seed as reference. Data from different donors were integrated using Seurat’s reciprocal principal components analysis (RPCA) pipeline. E Evaluation of cell type annotations using CellRef cell type markers from Supplementary Data 2. F Percentages of cell type markers (Supplementary Data 2) that are differentially expressed in their corresponding cell type predictions (n = 34 cell types) in (D). Data are shown using violin plot with dot and error bars representing mean ± SEM. G Heatmap visualization of expression of cell type specific differentially expressed genes (DEGs). H The number of DEGs for each predicted cell type. I Significantly enriched functional annotations using DEGs of the predicted AT2 cells: most enriched Gene Ontology Biological Processes (top) and ToppCell Gene Sets (bottom). Functional enrichment analysis was performed using ToppGene (https://toppgene.cchmc.org/enrichment.jsp). The minimum false discovery rate (FDR) was set to 1e−300. Please see Fig. 2 for definitions of cell type abbreviations.
Fig. 5 Cell type annotation and evaluation using the LungMAP Human Lung CellRef. A Schematic workflow of the automated cell type annotation and evaluation pipeline. B Distributions of cell type prediction scores in each test data. Prediction scores using CellRef Seed (yellow bars) are comparable to those using the complete CellRef (blue bars). Prediction scores (between 0 and 1) were calculated by the Seurat v4 MapQuery function for each cell. Box center lines, bounds of the box, and whiskers indicate medians, first and third quartiles, and minimum and maximum values within 1.5×IQR (interquartile range) of the box limits, respectively. GSM5388411: 6228 cells, GSM5388412: 8329 cells, GSM5388413: 7143 cells, GSM4504966: 8381 cells, GSM4504967: 8043 cells, GSM4035472: 5767 cells. C Consistency of cell type predictions using the CellRef Seed and CellRef in each test data. Consistency percentages (y axis) were calculated for cells in each test dataset (color) passing different thresholds of prediction scores (x axis). D–H Evaluation of automated cell type annotations for three of our test data (GSM5388411/12/13, three scRNA-seq of normal human lungs). Evaluation of the other three test data samples were shown in Supplementary Fig. 7. Basal and suprabasal cells were combined in prediction. D UMAP visualization of cells with prediction scores ≥ default cut-off (mean-1 standard deviation) and predicted annotations with at least 5 cells. Cells were colored by automated cell type annotations using the CellRef Seed as reference. Data from different donors were integrated using Seurat’s reciprocal principal components analysis (RPCA) pipeline. E Evaluation of cell type annotations using CellRef cell type markers from Supplementary Data 2. F Percentages of cell type markers (Supplementary Data 2) that are differentially expressed in their corresponding cell type predictions (n = 34 cell types) in (D). Data are shown using violin plot with dot and error bars representing mean ± SEM. G Heatmap visualization of expression of cell type specific differentially expressed genes (DEGs). H The number of DEGs for each predicted cell type. I Significantly enriched functional annotations using DEGs of the predicted AT2 cells: most enriched Gene Ontology Biological Processes (top) and ToppCell Gene Sets (bottom). Functional enrichment analysis was performed using ToppGene (https://toppgene.cchmc.org/enrichment.jsp). The minimum false discovery rate (FDR) was set to 1e−300. Please see Fig. 2 for definitions of cell type abbreviations.CC BY · PubMed Central
Fig. 6 Application of LungMAP Human Lung CellRef to disease lungs. A UMAP visualization of a published scRNA-seq of human lungs with LAM1. Cell colors represent cell identities predicted in Guo et al., 2020, including a unique disease-related cell population, named LAMCORE cells (magenta cell cluster). B UMAP visualizations of cells predicted using the CellRef Seed as reference. Basal and suprabasal cells were combined in the prediction. Prediction scores (between 0 and 1) were calculated by the Seurat v4 MapQuery function for each cell. Cells with prediction score >= the default cutoff (i.e., the mean minus 1 standard deviation value) were shown. Three singleton cell type predictions were not included. C Evaluation of cell type predictions using expression of representative CellRef marker genes. Megaka./Platelet: Megakaryocyte/Platelet. D Distributions of the cell type prediction scores in each of the original cell identities (n = 18 cell types; abbreviations were defined in Guo et al.1). The black and red horizontal line represents the mean and (1 standard deviation lower than the mean) value of the prediction scores, respectively. E–G UMAP and boxplot visualizations of application of CellRef to a published scRNA-seq of human lungs with idiopathic pulmonary fibrosis (IPF)11. E UMAP visualization of cells predicted using the CellRef Seed. Basal and suprabasal were combined, T cell subsets, and monocyte subsets were combined in the prediction. F UMAP visualization of cells colored by the prediction scores. G Left: UMAP visualization of cells colored by the original cell identities (n = 31 cell types; abbreviations were defined in Habermann et al.11). Right: boxplot visualization of the distribution of prediction scores in each of the original cell identities. The black and red horizontal line represents the mean and (1 standard deviation lower than the mean) value of the prediction scores, respectively. The disease-associated KRT5-/KRT17+ cells had prediction scores below the cutoff line. The number of data points in each boxplot in (B) and (G) can be found in Fig. 6 source data table. In (D) and (G), Box center lines, bounds of the box, and whiskers indicate medians, first and third quartiles, and minimum and maximum values within 1.5×IQR (interquartile range) of the box limits, respectively. Please see Fig. 2 for definitions of CellRef cell type abbreviations.
Fig. 6 Application of LungMAP Human Lung CellRef to disease lungs. A UMAP visualization of a published scRNA-seq of human lungs with LAM1. Cell colors represent cell identities predicted in Guo et al., 2020, including a unique disease-related cell population, named LAMCORE cells (magenta cell cluster). B UMAP visualizations of cells predicted using the CellRef Seed as reference. Basal and suprabasal cells were combined in the prediction. Prediction scores (between 0 and 1) were calculated by the Seurat v4 MapQuery function for each cell. Cells with prediction score >= the default cutoff (i.e., the mean minus 1 standard deviation value) were shown. Three singleton cell type predictions were not included. C Evaluation of cell type predictions using expression of representative CellRef marker genes. Megaka./Platelet: Megakaryocyte/Platelet. D Distributions of the cell type prediction scores in each of the original cell identities (n = 18 cell types; abbreviations were defined in Guo et al.1). The black and red horizontal line represents the mean and (1 standard deviation lower than the mean) value of the prediction scores, respectively. E–G UMAP and boxplot visualizations of application of CellRef to a published scRNA-seq of human lungs with idiopathic pulmonary fibrosis (IPF)11. E UMAP visualization of cells predicted using the CellRef Seed. Basal and suprabasal were combined, T cell subsets, and monocyte subsets were combined in the prediction. F UMAP visualization of cells colored by the prediction scores. G Left: UMAP visualization of cells colored by the original cell identities (n = 31 cell types; abbreviations were defined in Habermann et al.11). Right: boxplot visualization of the distribution of prediction scores in each of the original cell identities. The black and red horizontal line represents the mean and (1 standard deviation lower than the mean) value of the prediction scores, respectively. The disease-associated KRT5-/KRT17+ cells had prediction scores below the cutoff line. The number of data points in each boxplot in (B) and (G) can be found in Fig. 6 source data table. In (D) and (G), Box center lines, bounds of the box, and whiskers indicate medians, first and third quartiles, and minimum and maximum values within 1.5×IQR (interquartile range) of the box limits, respectively. Please see Fig. 2 for definitions of CellRef cell type abbreviations.CC BY · PubMed Central
Fig. 7 Assessment of cell type prediction accuracy of the LungMAP Human Lung CellRef. A Heatmap visualization of Pearson’s correlations of cell types between the human lung CellRef and the Human Lung Cell Atlas (HLCA)37. A pseudo-bulk profile was created for each cell type of either CellRef or HLCA by averaging each gene’s expression in the cells of the cell type. Cell types were clustered into four modules, each corresponding to one of the four major cell-lineages. Correspondences of CellRef and HLCA cell types within each of the four modules were shown based on the hierarchical clustering analysis. B, C Assessment of cell type accuracy based on marker gene expression. B Area under the receiver operating characteristic (ROC) curve (AUC) values for each of the mapped cell types (n = 42) in CellRef (orange) and HLCA (blue) calculated using the cell type selective marker genes identified from the HLCA study. Left: summary of the AUC values using violin plots. Middle: AUC values for each of the mapped cell types. Right: using CellRef AF2 (HLCA adventitial fibroblasts) as an example to show the ROC curves labeled with AUC values and 90% confidence interval. C AUCs values for each of the mapped cell types (n = 42) in the CellRef (orange) and HLCA (blue) calculated using the cell type selective marker genes identified by CellRef (Supplementary Data 5). Left: summary of the AUC values using violin plots. Middle: AUC values for each of the mapped cell types. Right: using CellRef AF2 (HLCA adventitial fibroblasts) as an example to show the ROC curves labeled with AUC values and 90% confidence interval. In both (B) and (C), the black dot and error bars represent mean ± SEM. p value represents significance of difference assessed using two-tailed paired Welch’s t test. CellRef cell type abbreviations are described in Fig. 2.
Fig. 7 Assessment of cell type prediction accuracy of the LungMAP Human Lung CellRef. A Heatmap visualization of Pearson’s correlations of cell types between the human lung CellRef and the Human Lung Cell Atlas (HLCA)37. A pseudo-bulk profile was created for each cell type of either CellRef or HLCA by averaging each gene’s expression in the cells of the cell type. Cell types were clustered into four modules, each corresponding to one of the four major cell-lineages. Correspondences of CellRef and HLCA cell types within each of the four modules were shown based on the hierarchical clustering analysis. B, C Assessment of cell type accuracy based on marker gene expression. B Area under the receiver operating characteristic (ROC) curve (AUC) values for each of the mapped cell types (n = 42) in CellRef (orange) and HLCA (blue) calculated using the cell type selective marker genes identified from the HLCA study. Left: summary of the AUC values using violin plots. Middle: AUC values for each of the mapped cell types. Right: using CellRef AF2 (HLCA adventitial fibroblasts) as an example to show the ROC curves labeled with AUC values and 90% confidence interval. C AUCs values for each of the mapped cell types (n = 42) in the CellRef (orange) and HLCA (blue) calculated using the cell type selective marker genes identified by CellRef (Supplementary Data 5). Left: summary of the AUC values using violin plots. Middle: AUC values for each of the mapped cell types. Right: using CellRef AF2 (HLCA adventitial fibroblasts) as an example to show the ROC curves labeled with AUC values and 90% confidence interval. In both (B) and (C), the black dot and error bars represent mean ± SEM. p value represents significance of difference assessed using two-tailed paired Welch’s t test. CellRef cell type abbreviations are described in Fig. 2.CC BY · PubMed Central
Fig. 8 Assessment of cell type stability of automated annotation using CellRef. A, B UMAP projection of scRNA-seq (Travaglini et al.15, n = 3 human lungs) with Azimuth projected cell type annotations using the LungMAP Human Lung CellRef Seed (A) or using the Human Lung Cell Atlas (HLCA)37 (B) as the reference. C Corresponding cell-population assignments of CellRef and HLCA (mapping percentage relative to CellRef). D Cells colored by “winning” annotations from CellRef or HLCA determined by scTriangulate based on stability assessments (shown in E) annotations. E Violin plot visualization of stability metric scores calculated using scTriangulate, including reclassification accuracy (SCCAF and reassign) or marker gene specificity (TF-IDF score), for all Azimuth assigned CellRef or HLCA cell populations (n = 42 cell populations predicted using the CellRef Seed; n = 48 cell populations predicted using HLCA) in Travaglini et al. 2020. The black dots and error bars represent mean ± SEM. p value represents significance of difference assessed using two-tailed unpaired Welch’s t test. Please see Fig. 2 for definitions of CellRef cell type abbreviations.
Fig. 8 Assessment of cell type stability of automated annotation using CellRef. A, B UMAP projection of scRNA-seq (Travaglini et al.15, n = 3 human lungs) with Azimuth projected cell type annotations using the LungMAP Human Lung CellRef Seed (A) or using the Human Lung Cell Atlas (HLCA)37 (B) as the reference. C Corresponding cell-population assignments of CellRef and HLCA (mapping percentage relative to CellRef). D Cells colored by “winning” annotations from CellRef or HLCA determined by scTriangulate based on stability assessments (shown in E) annotations. E Violin plot visualization of stability metric scores calculated using scTriangulate, including reclassification accuracy (SCCAF and reassign) or marker gene specificity (TF-IDF score), for all Azimuth assigned CellRef or HLCA cell populations (n = 42 cell populations predicted using the CellRef Seed; n = 48 cell populations predicted using HLCA) in Travaglini et al. 2020. The black dots and error bars represent mean ± SEM. p value represents significance of difference assessed using two-tailed unpaired Welch’s t test. Please see Fig. 2 for definitions of CellRef cell type abbreviations.CC BY · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Citations

Cites

Show all 11 references

Research works printed but not linked: 35

  1. Guo, M. et al. Single-cell transcriptomic analysis identifies a unique pulmonary lymphangioleiomyomatosis cell. Am. J. Respir. Crit. Care Med. 202, 1373–1387 (2020). · DOI
  2. Wang, A. et al. Single-cell multiomic profiling of human lungs reveals cell-type-specific and age-dynamic control of SARS-CoV2 host genes. Elife 9, https://doi.org/10.7554/eLife.62522 (2020). · DOI
  3. Melms, J. C. et al. A molecular single-cell lung atlas of lethal COVID-19. Nature 595, 114–119 (2021). · DOI
  4. Basil, M. C. et al. Human distal airways contain a multipotent secretory cell that can regenerate alveoli. Nature https://doi.org/10.1038/s41586-022-04552-0 (2022). · DOI
  5. Abdelaal, T. et al. A comparison of automatic cell identification methods for single-cell RNA sequencing data. Genome Biol. 20, 194 (2019). · DOI · Paper page
  6. Osumi-Sutherland, D. et al. Cell type ontologies of the Human Cell Atlas. Nat. Cell Biol. 23, 1129–1135 (2021). · DOI
  7. Miller, J. A. et al. Common cell type nomenclature for the mammalian brain. Elife 9, https://doi.org/10.7554/eLife.59928 (2020). · DOI
  8. Sun, X. et al. A census of the lung: CellCards from LungMAP. Dev. Cell 57, 112–145.e112 (2022). · DOI
  9. Habermann, A. C. et al. Single-cell RNA sequencing reveals profibrotic roles of distinct epithelial and mesenchymal lineages in pulmonary fibrosis. Sci. Adv. 6, eaba1972 (2020). · DOI · Paper page
  10. Adams, T. S. et al. Single-cell RNA-seq reveals ectopic and aberrant lung-resident cell populations in idiopathic pulmonary fibrosis. Sci. Adv. 6, eaba1983 (2020). · DOI · Paper page
  11. Deprez, M. et al. A single-cell atlas of the human healthy airways. Am. J. Respir. Crit. Care Med. 202, 1636–1645 (2020). · DOI · Paper page
  12. Travaglini, K. J. et al. A molecular cell atlas of the human lung from single-cell RNA sequencing. Nature 587, 619–625 (2020). · DOI
  13. Goldfarbmuren, K. C. et al. Dissecting the cellular specificity of smoking effects and reconstructing lineages in the human airway epithelium. Nat. Commun. 11, 2485 (2020). · DOI
  14. Li, G. et al. Decision level integration of unimodal and multimodal single cell data with scTriangulate. Nat. Commun. 14, 406 (2023). · DOI
  15. Stuart, T. et al. Comprehensive Integration of Single-Cell Data. Cell 177, 1888–1902.e1821 (2019). · DOI · Paper page
  16. Schupp, J. C. et al. Integrated single-cell atlas of endothelial cells of the human lung. Circulation 144, 286–302 (2021). · DOI · Paper page
  17. Guo, M. et al. Single cell RNA analysis identifies cellular heterogeneity and adaptive responses of the lung at birth. Nat. Commun. 10, 37 (2019). · DOI
  18. Morrisey, E. E. & Hogan, B. L. Preparing for the first breath: genetic and cellular mechanisms in lung development. Dev. Cell 18, 8–23 (2010). · DOI
  19. Rawlins, E. L., Clark, C. P., Xue, Y. & Hogan, B. L. The Id2+ distal tip lung epithelium contains individual multipotent embryonic progenitor cells. Development 136, 3741–3745 (2009). · DOI
  20. Choi, J. et al. Inflammatory signals induce AT2 cell-derived damage-associated transient progenitors that mediate alveolar regeneration. Cell Stem Cell 27, 366–382.e367 (2020). · DOI
  21. Strunz, M. et al. Alveolar regeneration through a Krt8+ transitional stem cell state that persists in human lung fibrosis. Nat. Commun. 11, 3559 (2020). · DOI
  22. Wang, G. et al. Generation of pulmonary endothelial progenitor cells for cell-based therapy using interspecies mouse-rat chimeras. Am. J. Respir. Crit. Care Med. 204, 326–338 (2021). · DOI
  23. Guo, M., Wang, H., Potter, S. S., Whitsett, J. A. & Xu, Y. SINCERA: a pipeline for single-cell RNA-seq profiling analysis. PLoS Comput. Biol. 11, e1004575 (2015). · DOI
  24. Bridges, J. P. et al. Glucocorticoid regulates mesenchymal cell differentiation required for perinatal lung morphogenesis and function. Am. J. Physiol. Lung Cell Mol. Physiol. 319, L239–L255 (2020). · DOI
  25. Ouyang, J. F., Kamaraj, U. S., Cao, E. Y. & Rackham, O. J. L. ShinyCell: simple and sharable visualisation of single-cell gene expression data. Bioinformatics, https://doi.org/10.1093/bioinformatics/btab209 (2021). · DOI
  26. Kadur Lakshminarasimha Murthy, P. et al. Human distal lung maps and lineage hierarchies reveal a bipotent progenitor. Nature 604, 111–119 (2022). · DOI
  27. Zepp, J. A. et al. Genomic, epigenomic, and biophysical cues controlling the emergence of the lung alveolus. Science 371, https://doi.org/10.1126/science.abc3172 (2021). · DOI
  28. Madissoon, E. et al. scRNA-seq assessment of the human lung, spleen, and esophagus tissue stability after cold preservation. Genome Biol. 21, 1 (2019). · DOI
  29. Sikkema, L. et al. An integrated cell atlas of the human lung in health and disease. Preprint at bioRxiv https://doi.org/10.1101/2022.03.10.483747 (2022). · DOI
  30. Miao, Z. et al. Putative cell type discovery from single-cell gene expression data. Nat. Methods, https://doi.org/10.1038/s41592-020-0825-9 (2020). · DOI
  31. DePasquale, E. A. K. et al. cellHarmony: cell-level matching and holistic comparison of single-cell transcriptomes. Nucleic Acids Res. 47, e138 (2019). · DOI · Paper page
  32. Seumois, G. & Vijayanand, P. Single-cell analysis to understand the diversity of immune cell types that drive disease pathogenesis. J. Allergy Clin. Immunol. 144, 1150–1153 (2019). · DOI
  33. Wolock, S. L., Lopez, R. & Klein, A. M. Scrublet: computational identification of cell doublets in single-cell transcriptomic data. Cell Syst. 8, 281–291.e289 (2019). · DOI · Paper page
  34. Macosko, E. Z. et al. Highly parallel genome-wide expression profiling of individual cells using nanoliter droplets. Cell 161, 1202–1214 (2015). · DOI · Paper page
  35. Lun, A. T. L. et al. EmptyDrops: distinguishing cells from empty droplets in droplet-based single-cell RNA sequencing data. Genome Biol. 20, 63 (2019). · DOI

Cited by

See all 197

Metadata

Authors
Guo, M., Morley, M. P., Jiang, C., Wu, Y., Li, G., Du, Y., Zhao, S., Wagner, A., Cakar, A. C., Kouril, M., Jin, K., Gaddis, N., Kitzmiller, J. A., Stewart, K., Basil, M. C., Lin, S. M., Ying, Y., Babu, A., Wikenheiser-Brokamp, K. A., Mun, K. S., Naren, A. P., Clair, G., Adkins, J. N., Pryhuber, G. S., Misra, R. S., Aronow, B. J., Tickle, T. L., Salomonis, N., Sun, X., Morrisey, E. E., Whitsett, J. A., NHLBI LungMAP Consortium, Lin, S., & Xu, Y.
Venue
Nature Communications
Published
2023-07-29
Volume
14
Issue
1
ISSN
2041-1723
Cited by
197
References
11
Language
en
Funders
  • U.S. Department of Health & Human Services | NIH | National Heart, Lung, and Blood Institute (U01HL122642, U01HL148856, U01HL134745, R01HL153045)
  • U.S. Department of Health & Human Services | NIH | National Institute of Diabetes and Digestive and Kidney Diseases (P30DK117467)

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review

How to cite

Cite CoChat for the checks on this page, and the paper for its own findings.

This page
CoChat. (2026). Citation check: Guided construction of single cell reference for human and mouse lung. https://app.cochat.ai/papers/10.1038/s41467-023-40173-5
The paper
Guo, M., Morley, M. P., Jiang, C., Wu, Y., Li, G., Du, Y., Zhao, S., Wagner, A., Cakar, A. C., Kouril, M., Jin, K., Gaddis, N., Kitzmiller, J. A., Stewart, K., Basil, M. C., Lin, S. M., Ying, Y., Babu, A., Wikenheiser-Brokamp, K. A., . . . Xu, Y. (2023). Guided construction of single cell reference for human and mouse lung. Nature Communications. https://doi.org/10.1038/s41467-023-40173-5