Defining human mesenchymal and epithelial heterogeneity in response to oral inflammatory disease

  1. Ana J Caetano
  2. Val Yianni
  3. Ana Volponi
  4. Veronica Booth
  5. Eleanor M D'Agostino
  6. Paul Sharpe  Is a corresponding author
  1. Centre for Craniofacial and Regenerative Biology, Faculty of Dentistry, Oral & Craniofacial Sciences, King’s College London, United Kingdom
  2. Department of Periodontology, Faculty of Dentistry, Oral & Craniofacial Sciences, King’s College London, United Kingdom
  3. Unilever R&D, Colworth Science Park, Sharnbrook, Bedfordshire, United Kingdom
5 figures, 1 table and 4 additional files

Figures

Figure 1 with 3 supplements
Single-cell Atlas of Gingiva Biopsies from Healthy Individuals and Periodontitis Patients.

(A) Overview of the experimental workflow. All samples were processed immediately after clinical surgery. (B) scRNA-seq data obtained from healthy and periodontitis cells (n = 12,411) from four donors illustrated by UMAP coloured by cell-type annotation (nHealthy = 4639, nMild = 4401, nSevere = 3367). (C) Heatmap of the mean expression of the most differentially expressed marker genes for each cluster identified. (D) Haematoxylin and eosin staining of gingival sections from healthy, mild, and severe patient samples showing increasing changes in tissue architecture with loss of epithelial rete ridges definition and infiltration of leukocytes. (E) Changes in tissue composition in periodontitis showing UMAP of progressive diseased states from healthy, mild, and severely diseased donors.

Figure 1—figure supplement 1
Flow Cytometry Gating Strategies on Human Gingival Cells.

Cells were gated based on size using standard SSC-A and FSC-A parameters. Doublets were excluded using SSC-A and SSC-W parameters. Live cells were selected as cells dimly fluorescing in DAPI. Samples were analysed on BD FACS Aria III Fusion machine.

Figure 1—figure supplement 2
Single-cell profiling of healthy human gingiva datasets using 10x Chromium.

(A) UMAP visualisation of human gingiva clusters from healthy human donors. (B) Scatter plots showing differential expressed genes across the two healthy samples. Panels A-B, n = 2 individuals.

Figure 1—figure supplement 3
Single-cell profiling of healthy and disease human gingiva using 10x Chromium.

(A) UMAP illustration of scRNA-seq data obtained from healthy and periodontitis cells (n = 12,411) from four donors coloured by condition. (B) Cell subset distributions across conditions. Shown are percentage of cells (y axis) in each cell subset (bars) that are derived from each healthy (orange), mild (green), and severe (blue) samples. X axis represents cluster ID. (C) Feature Plot showing the expression of lineage marker genes used for cell-type classification.

Figure 2 with 1 supplement
Cellular and molecular map of the stromal gingival compartment in health and disease identifies subpopulations with potential role in disease progression.

(A) UMAP plot of gingival stromal cells. Single cells coloured by cluster annotation. (B) UMAP plot of stromal cells during disease progression. (C) Heatmap showing subset-specific markers. (D) GO enrichment terms for S1 (myofibroblast) and S3 (pericyte). -log adjusted p-value shown (dotted line corresponds to FDR = 0.05). (E) Immunofluorescence staining showing COLVI, ACTA2 expression throughout disease progression. Scale bars, 100 μm. n = 3 patient samples/condition. Feature plots showing COLVI and ACTA2 expression across clusters and conditions. (F) UMAP annotated with CytoTRACE analysis to predict stromal stem populations. Transcriptional diversity is used here to predict maturation states. (G) Immunohistochemistry staining showing AEBP1+ cells tissue distribution. Scale bar, 100 μm. n = 6 patient samples.

Figure 2—figure supplement 1
Re-clustering of human stromal gingival cells in health and disease, Related to Figure 2.

(A) Feature Plots showing the expression of individual genes used for cell-type assignment of different stromal subsets. (B) GO enrichment terms for the different stromal subsets. -log adjusted p-value shown (dotted line corresponds to FDR = 0.05).

Figure 3 with 1 supplement
Cellular and molecular map of the epithelial gingival compartment in health and disease.

(A) UMAP plot of human gingival epithelial cells. Single cells coloured by cluster annotation. (B) UMAP plot of epithelial cells during disease progression. (C) Heatmap showing subset-specific markers. (D) GO enrichment terms for E0, E6, and E8 with -log adjusted p-value shown (dotted line corresponds to FDR = 0.05). (E) Dot plot showing top predicted ligands expressed by epithelial cells that modulate the E0 (stem) compartment. (F) Expression of KI67 and SOX2 in human healthy tissue. KI67 marks proliferative cells (cluster E5), and SOX2 marks an epithelial stem cell compartment (cluster E0). Scale bars = 100 μm (A, B). Scale bars, 50 μm (A’, B’). n = 4 patient samples/condition.

Figure 3—figure supplement 1
Re-clustering of human epithelial gingival cells in health and disease, Related to Figure 3.

(A) Feature Plots showing the expression of individual genes used for cell-type assignment of different epithelial subsets. (B) GO enrichment terms for the different epithelial subsets. -log adjusted p-value shown (dotted line corresponds to FDR = 0.05).

Unbiased cell–cell interaction analysis and its effect in the stromal microenvironment.

(A) Schematic representation of the NicheNet analysis of upstream ligand–receptor pairs and stromal target genes inducing DE genes in periodontitis. Created with BioRender.com. (B) Dot plots depicting which gingival cell populations express top-ranked ligands contributing to the transcriptional response observed from health to mild disease and from mild to severe in the stromal compartment. (C) Schematic representation of the NicheNet analysis of epithelial-mesenchymal crosstalk in mild disease. Created with BioRender.com. (D) Top predicted epithelial ligands driving the stromal inflammatory response and dot plot showing which epithelial subpopulation express these ligands. (E) Ligand-receptor heatmap of potential receptors expressed by stromal cells associated with each epithelial ligand. (F) Ligand-target heatmap of stromal and perivascular target genes of the identified epithelial ligands. (G) Dot plot confirming upregulation of the identified stromal target genes in disease. (H) UMAPs feature plots mapping the identified epithelial ligands and target genes to the respective target genes expressed by stromal cells.

Figure 5 with 1 supplement
Periodontitis induces an IgG plasma B-cell signature in human gingiva.

(A) UMAP analysis of human B cells identifying follicular and IgG plasma B cells split by condition. (B) UMAP expression plots of human B-cell subset markers. Cells coloured by normalised expression of indicated genes. (C) CD79A in human gingival tissue across health and disease. Scale bars, 100 μm. n = 3 patient samples/condition. (D) Gene enrichment analysis of IgG plasma B cells. -log adjusted p-value shown (dotted line corresponds to FDR = 0.05).

Figure 5—figure supplement 1
Re-clustering of human B and T cells in health and disease.

(A) UMAP plot of B cells in health and disease. (B) Heatmap showing subset-specific markers. (C) UMAP plot of B cells during disease progression. (D) Feature Plots showing the expression of individual genes used for assignment of different B-cell subsets. (E) GO enrichment terms for IgG plasma B cells with -log adjusted p-value shown (dotted line corresponds to FDR = 0.05). (F) UMAP plot of T cells in health and disease. (G) UMAP plot of T cells during disease progression. (H) Feature Plots showing the expression of CD4 and CD8 T-cell subsets. (I) UMAP showing number of detected genes.

Tables

Key resources table
Reagent type
(species) or
resource
DesignationSource or
reference
IdentifiersAdditional
information
Biological sample (Human)Human gingival biopsiesPeriodontology department, King’s College London
Antibodyanti-COLVI (Rabbit monoclonal)AbcamCat #ab182744, RRID:AB_2847919IHC (1:500)
Antibodyanti-ACTA2 (Mouse monoclonal)AbcamCat #ab7817, RRID:AB_262054IHC (1:200)
antibodyanti-MCAM (Rabbit monoclonal)AbcamCat #ab75769, RRID:AB_2143375IHC (1:100)
Antibodyanti-KI67 (Rabbit polyclonal)AbcamCat #ab15580, RRID:AB_443209IHC (1:100)
Antibodyanti-SOX2
(Rabbit monoclonal)
AbcamCat# ab92494, RRID:AB_10585428IHC (1:100)
Antibodyanti-CD79A (Rabbit monoclonal)AbcamCat# ab79414,
RRID:AB_2260147
IHC (1:100)
Antibodyanti-AEBP1 (Rabbit polyclonal)Atlas AntibodiesCat# HPA064970, RRID:AB_2685394IHC (1:200)
Commercial assay or kitImmPRESS Excel Staining Kit, Anti-Rabbit IgVector LaboratoriesCat# MP-7601, RRID:AB_2336533
Chemical compound, drugUltraPure BSA (50 mg/mL)ThermoFisher ScientificCat# AM26180.04%
Commercial assay or kitChromium Single Cell 3’ Library and Gel Bead Kit v310X GenomicsCat# PN-1000092
Commercial assay or kitChromium Single Cell B Chip Kit10X GenomicsCat# PN-1000074
Commercial assay or kitWhole Skin Dissociation Kit, humanMiltenyi BiotecCat# 130-101-540
Software, algorithmCellRanger Version 410X GenomicsRRID:SCR_017344
Software, algorithmSeurat Version 3.0R BioconductorRRID:SCR_007322
https://satijalab.org/seurat/
Software, algorithmEnrichrChen et al., 2013RRID:SCR_001575
Software, algorithmCytoTRACER Bioconductorhttps://cytotrace.stanford.edu
Software, algorithmNicheNetGitHubhttps://github.com/saeyslab/nichenetr
OtherGRCh38CellRanger, 10X Genomicshttps://support.10xgenomics.com/single-cell-gene-expression/software/downloads/latest
OtherDAPI stainInvitrogenD1306, RRID:AB_2629482(1 ug/mL)

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Ana J Caetano
  2. Val Yianni
  3. Ana Volponi
  4. Veronica Booth
  5. Eleanor M D'Agostino
  6. Paul Sharpe
(2021)
Defining human mesenchymal and epithelial heterogeneity in response to oral inflammatory disease
eLife 10:e62810.
https://doi.org/10.7554/eLife.62810