Establishment of inclusive single-cell transcriptome atlases from mouse and human tooth as powerful resource for dental research

Single-cell and spatially resolved interactomics of tooth-associated keratinocytes in periodontitis

Periodontitis affects billions of people worldwide. To address relationships of periodontal niche cell types and microbes in periodontitis, we generated an integrated single-cell RNA sequencing (scRNAseq) atlas of human periodontium (34-sample, 105918-cell), including sulcular and junctional keratinocytes (SK/JKs). SK/JKs displayed altered differentiation states and were enriched for effector cytokines in periodontitis. Single-cell metagenomics revealed 37 bacterial species with cell-specific tropism. Fluorescence in situ hybridization detected intracellular 16 S and mRNA signals of multiple species and correlated with SK/JK proinflammatory phenotypes in situ. Cell-cell communication analysis predicted keratinocyte-specific innate and adaptive immune interactions. Highly multiplexed immunofluorescence (33-antibody) revealed peri-epithelial immune foci, with innate cells often spatially constrained around JKs. Spatial phenotyping revealed immunosuppressed JK-microniches and SK-localized tertiary lymphoid structures in periodontitis. Here, we demonstrate impacts on and predicted interactomics of SK and JK cells in health and periodontitis, which requires further investigation to support precision periodontal interventions in states of chronic inflammation.

Spatiotemporal cell landscape of human embryonic tooth development

Understanding the cellular composition and trajectory of human tooth development is valuable for dentistry and stem cell engineering research. Previous single-cell studies have focused on mature human teeth and developing mouse teeth, but the cell landscape of human embryonic dental development is still unknown. In this study, tooth germ tissues were collected from aborted foetus (17-24 weeks) for single-cell RNA sequence and spatial transcriptome analysis. The cells were classified into seven subclusters of epithelium, and seven clusters of mesenchyme, as well as other cell types such as Schwann cell precursor and pericyte. For epithelium, the stratum intermedium branch and the ameloblast branch diverged from the same set of outer enamel-inner enamel-ALCAM+ epithelial cell lineage, but their spatial distribution of two branches was not clearly distinct. This trajectory received spatially adjacent regulation signals from mesenchyme and pericyte, including JAG1 and APP. The differentiation of pulp cell and pre-odontoblast showed four waves of temporally distinct gene expression, which involved regulation networks of LHX9, DLX5 and SP7, and these genes were regulated by upstream ligands such as the BMP family. This provides a reference landscape for the research on early human tooth development, covering different spatial structures and developmental periods.

The temporal protein signature analyses of developing human deciduous molar tooth germ

The tooth serves as an exemplary model for developmental studies, encompassing epithelial-mesenchymal transition and cell differentiation. The essential factors and pathways identified in tooth development will help understand the natural development process and the malformations of mineralized tissues such as skeleton. The time-dependent proteomic changes were investigated through the proteomics of healthy human molars during embryonic stages, ranging from the cap-to-early bell stage. A comprehensive analysis revealed 713 differentially expressed proteins (DEPs) exhibiting five distinct temporal expression patterns. Through the application of weighted gene co-expression network analysis (WGCNA), 24 potential driver proteins of tooth development were screened, including CHID1, RAP1GDS1, HAPLN3, AKAP12, WLS, GSS, DDAH1, CLSTN1, AFM, RBP1, AGO1, SET, HMGB2, HMGB1, ANP32A, SPON1, FREM1, C8B, PRPS2, FCHO2, PPP1R12A, GPALPP1, U2AF2, and RCC2. Then, the proteomics and transcriptomics expression patterns of these proteins were further compared, complemented by single-cell RNA-sequencing (scRNA-seq). In summary, this study not only offers a wealth of information regarding the molecular intricacies of human embryonic epithelial and mesenchymal cell differentiation but also serves as an invaluable resource for future mechanistic inquiries into tooth development.

Profiles of Wnt pathway gene expression during tooth morphogenesis

Mouse and human genetic studies indicate key roles of the Wnt10a ligand in odontogenesis. Previous studies have identified effectors and regulators of the Wnt signaling pathway actively expressed during key stages of tooth morphogenesis. However, limitations in multiplexing and spatial resolution hindered a more comprehensive analysis of these signaling molecules. Here, profiling of transcriptomes using fluorescent multiplex in situ hybridization and single-cell RNA-sequencing (scRNA-seq) provide robust insight into the synchronized expression patterns of Wnt10a, Dkk1, and Sost simultaneously during tooth development. First, we identified Wnt10a transcripts restricted to the epithelium at the stage of tooth bud morphogenesis, contrasting that of Sost and Dkk1 localization to the dental mesenchyme. By embryonic day 15.5 (E15.5), a marked shift of Wnt10a expression from dental epithelium to mesenchyme was noted, while Sost and Dkk1 expression remained enriched in the mesenchyme. By postnatal day 0 (P0), co-localization patterns of Wnt10a, Dkk1, and Sost were observed in both terminally differentiating and secreting odontoblasts of molars and incisors. Interestingly, Wnt10a exhibited robust expression in fully differentiated ameloblasts at the developing cusp tip of both molars and incisors, an observation not previously noted in prior studies. At P7 and 14, after the mineralization of dentin and enamel, Wnt10a expression was limited to odontoblasts. Meanwhile, Wnt modulators showed reduced or absent signals in molars. In contrast, strong signals persisted in ameloblasts (for Wnt10a) and odontoblasts (for Wnt10a, Sost, and Dkk1) towards the proximal end of incisors, near the cervical loop. Our scRNA-seq analysis used CellChat to further contextualize Wnt pathway-mediated communication between cells by examining ligand-receptor interactions among different clusters. The co-localization pattern of Wnt10a, Dkk1, and Sost in both terminally differentiating and secreting odontoblasts of molars and incisors potentially signifies the crucial ligand-modulator interaction along the gradient of cytodifferentiation starting from each cusp tip towards the apical region. These data provide cell type-specific insight into the role of Wnt ligands and mediators during epithelial-mesenchymal interactions in odontogenesis.

Computational identification of key transcription factors for embryonic and postnatal Sox2+ dental epithelial stem cell

While many reptiles can replace their tooth throughout life, human loss the tooth replacement capability after formation of the permanent teeth. It was thought that the difference in tooth regeneration capability depends on the persistence of a specialized dental epithelial structure, the dental lamina that contains dental epithelial stem cells (DESC). Currently, we know very little about DESC such as what genes are expressed or its chromatin accessibility profile. Multiple markers of DESC have been proposed such as Sox2 and Lgr5 . Few single cell RNA-seq experiments have been performed previously, but no obvious DESC cluster was identified in these scRNA-seq datasets, possible due to that the expression level of DESC markers such as Sox2 and Lgr5 is too low or the percentage of DESC is too low in whole tooth. We utilize a mouse line Sox2-GFP to enrich Sox2+ DESC and use Smart-Seq2 protocol and ATAC-seq protocol to generate transcriptome profile and chromatin accessibility profile of P2 Sox2+ DESC. Additionally, we generate transcriptome profile and chromatin accessibility profile of E11.5 Sox2+ dental lamina cells. With transcriptome profile and chromatin accessibility profile, we systematically identify potential key transcription factors for E11.5 Sox2+ cells and P2 Sox2+ cells. We identified transcription factors including Pitx2, Id3, Pitx1, Tbx1, Trp63, Nkx2-3, Grhl3, Dlx2, Runx1, Nfix, Zfp536 , etc potentially formed the core transcriptional regulatory networks of Sox2+ DESC in both embryonic and postnatal stages.

Organoids from mouse molar and incisor as new tools to study tooth-specific biology and development

Organoid models provide powerful tools to study tissue biology and development in a dish. Presently, organoids have not yet been developed from mouse tooth. Here, we established tooth organoids (TOs) from early-postnatal mouse molar and incisor, which are long-term expandable, express dental epithelium stem cell (DESC) markers, and recapitulate key properties of the dental epithelium in a tooth-type-specific manner. TOs display in vitro differentiation capacity toward ameloblast-resembling cells, even more pronounced in assembloids in which dental mesenchymal (pulp) stem cells are combined with the organoid DESCs. Single-cell transcriptomics supports this developmental potential and reveals co-differentiation into junctional epithelium- and odontoblast-/cementoblast-like cells in the assembloids. Finally, TOs survive and show ameloblast-resembling differentiation also in vivo. The developed organoid models provide new tools to study mouse tooth-type-specific biology and development and gain deeper molecular and functional insights that may eventually help to achieve future human biological tooth repair and replacement.