System-level analyses of keystone genes required for mammalian tooth development

Unraveling the hidden complexity: Exploring dental tissues through single-cell transcriptional profiling

Understanding the composition and function of cells constituting tissues and organs is vital for unraveling biological processes. Single-cell analysis has allowed us to move beyond traditional methods of categorizing cell types. This innovative technology allows the transcriptional and epigenetic profiling of numerous individual cells, leading to significant insights into the development, homeostasis, and pathology of various organs and tissues in both animal models and human samples. In this review, we delve into the outcomes of major investigations using single-cell transcriptomics to decipher the cellular composition of mammalian teeth and periodontal tissues. The recent single-cell transcriptome-based studies have traced in detail the dental epithelium-ameloblast lineage and dental mesenchyme lineages in the mouse incisors and the tooth germ of both mice and humans; unraveled the microenvironment, the identity of niche cells, and cellular intricacies in the dental pulp; shed light on the molecular mechanisms orchestrating root formation; and characterized cellular dynamics of the periodontal ligament. Additionally, cellular components in dental pulps were compared between healthy and carious teeth at a single-cell level. Each section of this review contributes to a comprehensive understanding of tooth biology, offering valuable insights into developmental processes, niche cell identification, and the molecular secrets of the dental environment.

Bank vole genomics links determinate and indeterminate growth of teeth

BACKGROUND

Continuously growing teeth are an important innovation in mammalian evolution, yet genetic regulation of continuous growth by stem cells remains incompletely understood. Dental stem cells responsible for tooth crown growth are lost at the onset of tooth root formation. Genetic signaling that initiates this loss is difficult to study with the ever-growing incisor and rooted molars of mice, the most common mammalian dental model species, because signals for root formation overlap with signals that pattern tooth size and shape (i.e., cusp patterns). Bank and prairie voles (Cricetidae, Rodentia, Glires) have evolved rooted and unrooted molars while retaining similar size and shape, providing alternative models for studying roots.

RESULTS

We assembled a de novo genome of Myodes glareolus, a vole with high-crowned, rooted molars, and performed genomic and transcriptomic analyses in a broad phylogenetic context of Glires (rodents and lagomorphs) to assess differential selection and evolution in tooth forming genes. Bulk transcriptomics comparisons of embryonic molar development between bank voles and mice demonstrated overall conservation of gene expression levels, with species-specific differences corresponding to the accelerated and more extensive patterning of the vole molar. We leverage convergent evolution of unrooted molars across the clade to examine changes that may underlie the evolution of unrooted molars. We identified 15 dental genes with changing synteny relationships and six dental genes undergoing positive selection across Glires, two of which were undergoing positive selection in species with unrooted molars, Dspp and Aqp1. Decreased expression of both genes in prairie voles with unrooted molars compared to bank voles supports the presence of positive selection and may underlie differences in root formation.

CONCLUSIONS

Our results support ongoing evolution of dental genes across Glires and identify candidate genes for mechanistic studies of root formation. Comparative research using the bank vole as a model species can reveal the complex evolutionary background of convergent evolution for ever-growing molars.

Exploring the role of DNMT1 in dental papilla cell fate specification during mouse tooth germ development through integrated single-cell transcriptomics and bulk RNA sequencing

OBJECTIVES

This study aimed to investigate the regulatory mechanisms governing dental mesenchymal cell commitment during tooth development, focusing on odontoblast differentiation and the role of epigenetic regulation in this process.

METHODS

We performed single-cell RNA sequencing (scRNA-seq) of dental cells from embryonic day 14.5 (E14.5) mice to understand the heterogeneity of developing tooth germ cells. Computational analyses including gene regulatory network (GRN) assessment were conducted. We validated our findings using immunohistochemistry (IHC) and in vitro loss-of-function analyses using the DNA methyltransferase 1 (DNMT1) inhibitor Gsk-3484862 in primary dental mesenchymal cells (DMCs) isolated from E14.5 mouse tooth germs. Bulk RNA-seq of Gsk-3484862-treated DMCs was performed to identify potential downstream targets of DNMT1.

RESULTS

scRNA-seq analysis revealed diverse cell populations within the tooth germs, including epithelial, mesenchymal, immune, and muscle cells. Using single-cell regulatory network inference and clustering (SCENIC), we identified Dnmt1 as a key regulator of early odontoblast development. IHC analysis showed the ubiquitous expression of DNMT1 in the dental papilla and epithelium. Bulk RNA-seq of cultured DMCs showed that Gsk-3484862 treatment upregulated odontoblast-related genes, whereas genes associated with cell division and the cell cycle were downregulated. Integrated analysis of bulk RNA-seq data with scRNA-seq SCENIC profiles was used to identify the potential Dnmt1 target genes.

CONCLUSIONS

Dnmt1 may negatively affect odontoblast commitment and differentiation during tooth development. These findings contribute to a better understanding of the molecular mechanisms underlying tooth development and future development of hard-tissue regenerative therapies.

Vole genomics links determinate and indeterminate growth of teeth

Continuously growing teeth are an important innovation in mammalian evolution, yet genetic regulation of continuous growth by stem cells remains incompletely understood. Dental stem cells responsible for tooth crown growth are lost at the onset of tooth root formation. Genetic signaling that initiates this loss is difficult to study with the ever-growing incisor and rooted molars of mice, the most common mammalian dental model species, because signals for root formation overlap with signals that pattern tooth size and shape (i.e., cusp patterns). Different species of voles (Cricetidae, Rodentia, Glires) have evolved rooted and unrooted molars that have similar size and shape, providing alternative models for studying roots. We assembled a de novo genome of Myodes glareolus, a vole with high-crowned, rooted molars, and performed genomic and transcriptomic analyses in a broad phylogenetic context of Glires (rodents and lagomorphs) to assess differential selection and evolution in tooth forming genes. We identified 15 dental genes with changing synteny relationships and six dental genes undergoing positive selection across Glires, two of which were undergoing positive selection in species with unrooted molars, Dspp and Aqp1. Decreased expression of both genes in prairie voles with unrooted molars compared to bank voles supports the presence of positive selection and may underlie differences in root formation. Bulk transcriptomics analyses of embryonic molar development in bank voles also demonstrated conserved patterns of dental gene expression compared to mice, with species-specific variation likely related to developmental timing and morphological differences between mouse and vole molars. Our results support ongoing evolution of dental genes across Glires, revealing the complex evolutionary background of convergent evolution for ever-growing molars.

Is there a loophole in Dollo's law? A DevoEvo perspective on irreversibility (of felid dentition)

There is a longstanding interest in whether the loss of complex characters is reversible (so-called "Dollo's law"). Reevolution has been suggested for numerous traits but among the first was Kurtén, who proposed that the presence of the second lower molar (M2 ) of the Eurasian lynx (Lynx lynx) was a violation of Dollo's law because all other Felids lack M2 . While an early and often cited example for the reevolution of a complex trait, Kurtén and Werdelin used an ad hoc parsimony argument to support their. Here I revisit the evidence that M2 reevolved lynx using explicit parsimony and maximum likelihood models of character evolution and find strong evidence that Kurtén and Werdelin were correct-M2 reevolved in E. lynx. Next, I explore the developmental mechanisms which may explain this violation of Dollo's law and suggest that the reevolution of lost complex traits may arise from the reevolution of cis-regulatory elements and protein-protein interactions, which have a longer half-life after silencing that protein coding genes. Finally, I present a developmental model to explain the reevolution M2 in E. lynx, which suggest that the developmental programs required for the establishment of serially homologous characters may never really be lost so long as a single instance of the character remains-thus the gain and loss and regain of serially homologous characters, such mammalian molars, may be developmentally and evolutionarily "simple."

Transcriptomic landscape of early hair follicle and epidermal development

Morphogenesis of ectodermal organs, such as hair, tooth, and mammary gland, starts with the formation of local epithelial thickenings, or placodes, but it remains to be determined how distinct cell types and differentiation programs are established during ontogeny. Here, we use bulk and single-cell transcriptomics and pseudotime modeling to address these questions in developing hair follicles and epidermis and produce a comprehensive transcriptomic profile of cellular populations in the hair placode and interplacodal epithelium. We report previously unknown cell populations and marker genes, including early suprabasal and genuine interfollicular basal markers, and propose the identity of suprabasal progenitors. By uncovering four different hair placode cell populations organized in three spatially distinct areas, with fine gene expression gradients between them, we posit early biases in cell fate establishment. This work is accompanied by a readily accessible online tool to stimulate further research on skin appendages and their progenitors.

Gene expression detection in developing mouse tissue using in situ hybridization and µCT imaging

High resolution and noninvasiveness have made soft-tissue X-ray microtomography (µCT) a widely applicable three-dimensional (3D) imaging method in studies of morphology and development. However, scarcity of molecular probes to visualize gene activity with µCT has remained a challenge. Here, we apply horseradish peroxidase-assisted reduction of silver and catalytic gold enhancement of the silver deposit to in situ hybridization in order to detect gene expression in developing tissues with µCT (here called GECT, gene expression CT). We show that GECT detects expression patterns of collagen type II alpha 1 and sonic hedgehog in developing mouse tissues comparably with an alkaline phosphatase-based detection method. After detection, expression patterns are visualized with laboratory µCT, demonstrating that GECT is compatible with varying levels of gene expression and varying sizes of expression regions. Additionally, we show that the method is compatible with prior phosphotungstic acid staining, a conventional contrast staining approach in µCT imaging of soft tissues. Overall, GECT is a method that can be integrated with existing laboratory routines to obtain spatially accurate 3D detection of gene expression.

Clinical, genetic, epidemiologic, evolutionary, and functional delineation of TSPEAR-related autosomal recessive ectodermal dysplasia 14

TSPEAR variants cause autosomal recessive ectodermal dysplasia (ARED) 14. The function of TSPEAR is unknown. The clinical features, the mutation spectrum, and the underlying mechanisms of ARED14 are poorly understood. Combining data from new and previously published individuals established that ARED14 is primarily characterized by dental anomalies such as conical tooth cusps and hypodontia, like those seen in individuals with WNT10A-related odontoonychodermal dysplasia. AlphaFold-predicted structure-based analysis showed that most of the pathogenic TSPEAR missense variants likely destabilize the β-propeller of the protein. Analysis of 100000 Genomes Project (100KGP) data revealed multiple founder TSPEAR variants across different populations. Mutational and recombination clock analyses demonstrated that non-Finnish European founder variants likely originated around the end of the last ice age, a period of major climatic transition. Analysis of gnomAD data showed that the non-Finnish European population TSPEAR gene-carrier rate is ∼1/140, making it one of the commonest AREDs. Phylogenetic and AlphaFold structural analyses showed that TSPEAR is an ortholog of drosophila Closca, an extracellular matrix-dependent signaling regulator. We, therefore, hypothesized that TSPEAR could have a role in enamel knot, a structure that coordinates patterning of developing tooth cusps. Analysis of mouse single-cell RNA sequencing (scRNA-seq) data revealed highly restricted expression of Tspear in clusters representing enamel knots. A tspeara -/-;tspearb -/- double-knockout zebrafish model recapitulated the clinical features of ARED14 and fin regeneration abnormalities of wnt10a knockout fish, thus suggesting interaction between tspear and wnt10a. In summary, we provide insights into the role of TSPEAR in ectodermal development and the evolutionary history, epidemiology, mechanisms, and consequences of its loss of function variants.

A Genome-wide association study of premolar agenesis in a chinese population

OBJECTIVE

Premolar agenesis is a common subtype of tooth agenesis. Although a genome-wide study (GWAS) has identified some variants involved in tooth agenesis in Europeans, the genetic mutation related to premolar agenesis in the Chinese population remains unclear.

MATERIALS AND METHODS

We present a GWAS in 218 premolar agenesis cases and 1,222 controls using the Illumina Infinium^® Global Screening Array. 5,585,618 single nucleotide polymorphisms (SNPs) were used for tests of associations with premolar agenesis.

RESULTS

Four independent SNPs on chromosome 2 were identified as susceptibility loci, including rs147680216, rs79743039, rs60540881, and rs6738629. The genome-wide significant SNP rs147680216 (p = 6.09 × 10^-9 ) was predicted to change the structure of the WNT10A protein and interact with hedgehog signaling pathway components. Meta-analysis showed that the rs147680216 A allele significantly increased the risk of tooth agenesis (p = 0.000). The other three SNPs with nominal significance are novel susceptibility loci. Of them, rs6738629 (p = 5.40 × 10^-6 ) acts as a potential transcriptional regulator of GCC2, a gene playing a putative role in dental and craniofacial development.

CONCLUSION

Our GWAS indicates that rs147680216 and additional three novel susceptibility loci on chromosome 2 are associated with the risk of premolar agenesis in the Chinese population.

Domestic cat embryos reveal unique transcriptomes of developing incisor, canine, and premolar teeth

Division of the dentition into morphologically distinct classes of teeth (incisors, canines, premolars, and molars) and the acquisition of tribosphenic molars facilitated precise occlusion between the teeth early in mammal evolution. Despite the evolutionary and ecological importance of distinct classes of teeth with unique cusp, crest, and basin morphologies, relatively little is known about the genetic basis for the development of different tooth classes within the embryo. Here we investigated genetic differences between developing deciduous incisor, canine, and premolar teeth in the domestic cat (Felis catus), which we propose to be a new model for tooth development. We examined differences in both developmental timing and crown morphology between the three tooth classes. Using RNA sequencing of early bell stage tooth germs, we showed that each of the three deciduous tooth classes possess a unique transcriptional profile. Three notable groups of genes emerged from our differential expression analysis; genes involved in the extracellular matrix (ECM), Wnt pathway signaling, and members of multiple homeobox gene families (Lhx, Dlx, Alx, and Nkx). Our results suggest that ECM genes may play a previously under-appreciated role in shaping the surface of the tooth crown during development. Differential regulation of these genes likely underlies differences in tooth crown shape and size, although subtle temporal differences in development between the tooth germs could also be responsible. This study provides foundational data for future experiments to examine the function of these candidate genes in tooth development to directly test their potential effects on crown morphology.

Integration of multimodal data in the developing tooth reveals candidate regulatory loci driving human odontogenic phenotypes

Human odontogenic aberrations such as abnormal tooth number and delayed tooth eruption can occur as a symptom of rare syndromes or, more commonly, as nonsyndromic phenotypes. These phenotypes can require extensive and expensive dental treatment, posing a significant burden. While many dental phenotypes are heritable, most nonsyndromic cases have not been linked to causal genes. We demonstrate the novel finding that common sequence variants associated with human odontogenic phenotypes are enriched in developmental craniofacial enhancers conserved between human and mouse. However, the bulk nature of these samples obscures if this finding is due to the tooth itself or the surrounding tissues. We therefore sought to identify enhancers specifically active in the tooth anlagen and quantify their contribution to the observed genetic enrichments. We systematically identified 22,001 conserved enhancers active in E13.5 mouse incisors using ChIP-seq and machine learning pipelines and demonstrated biologically relevant enrichments in putative target genes, transcription factor binding motifs, and in vivo activity. Multi-tissue comparisons of human and mouse enhancers revealed that these putative tooth enhancers had the strongest enrichment of odontogenic phenotype-associated variants, suggesting a role for dysregulation of tooth developmental enhancers in human dental phenotypes. The large number of these regions genome-wide necessitated prioritization of enhancer loci for future investigations. As enhancers modulate gene expression, we prioritized regions based on enhancers' putative target genes. We predicted these target genes and prioritized loci by integrating chromatin state, bulk gene expression and coexpression, GWAS variants, and cell type resolved gene expression to generate a prioritized list of putative odontogenic phenotype-driving loci active in the developing tooth. These genomic regions are of particular interest for downstream experiments determining the role of specific dental enhancer:gene pairs in odontogenesis.

An epithelial signalling centre in sharks supports homology of tooth morphogenesis in vertebrates

Development of tooth shape is regulated by the enamel knot signalling centre, at least in mammals. Fgf signalling regulates differential proliferation between the enamel knot and adjacent dental epithelia during tooth development, leading to formation of the dental cusp. The presence of an enamel knot in non-mammalian vertebrates is debated given differences in signalling. Here, we show the conservation and restriction of fgf3, fgf10, and shh to the sites of future dental cusps in the shark (Scyliorhinus canicula), whilst also highlighting striking differences between the shark and mouse. We reveal shifts in tooth size, shape, and cusp number following small molecule perturbations of canonical Wnt signalling. Resulting tooth phenotypes mirror observed effects in mammals, where canonical Wnt has been implicated as an upstream regulator of enamel knot signalling. In silico modelling of shark dental morphogenesis demonstrates how subtle changes in activatory and inhibitory signals can alter tooth shape, resembling developmental phenotypes and cusp shapes observed following experimental Wnt perturbation. Our results support the functional conservation of an enamel knot-like signalling centre throughout vertebrates and suggest that varied tooth types from sharks to mammals follow a similar developmental bauplan. Lineage-specific differences in signalling are not sufficient in refuting homology of this signalling centre, which is likely older than teeth themselves.

Spatiotemporal Changes in Transcriptome of Odontogenic and Non-odontogenic Regions in the Dental Arch of Mus musculus

Over the past 40 years, studies on tooth regeneration have been conducted. These studies comprised two main flows: some focused on epithelial-mesenchymal interaction in the odontogenic region, whereas others focused on creating a supernumerary tooth in the non-odontogenic region. Recently, the scope of the research has moved from conventional gene modification and molecular therapy to genome and transcriptome sequencing analyses. However, these sequencing data have been produced only in the odontogenic region. We provide RNA-Seq data of not only the odontogenic region but also the non-odontogenic region, which loses tooth-forming capacity during development and remains a rudiment. Sequencing data were collected from mouse embryos at three different stages of tooth development. These data will expand our understanding of tooth development and will help in designing developmental and regenerative studies from a new perspective.

Chromosomal neighbourhoods allow identification of organ specific changes in gene expression

Although most genes share their chromosomal neighbourhood with other genes, distribution of genes has not been explored in the context of individual organ development; the common focus of developmental biology studies. Because developmental processes are often associated with initially subtle changes in gene expression, here we explored whether neighbouring genes are informative in the identification of differentially expressed genes. First, we quantified the chromosomal neighbourhood patterns of genes having related functional roles in the mammalian genome. Although the majority of protein coding genes have at least five neighbours within 1 Mb window around each gene, very few of these neighbours regulate development of the same organ. Analyses of transcriptomes of developing mouse molar teeth revealed that whereas expression of genes regulating tooth development changes, their neighbouring genes show no marked changes, irrespective of their level of expression. Finally, we test whether inclusion of gene neighbourhood in the analyses of differential expression could provide additional benefits. For the analyses, we developed an algorithm, called DELocal that identifies differentially expressed genes by comparing their expression changes to changes in adjacent genes in their chromosomal regions. Our results show that DELocal removes detection bias towards large changes in expression, thereby allowing identification of even subtle changes in development. Future studies, including the detection of differential expression, may benefit from, and further characterize the significance of gene-gene neighbour relationships.

Development of organs is typically associated with small and hard to detect changes in gene expression. Here we examined how often genes regulating specific organs are neighbours to each other in the genome, and whether this gene neighbourhood helps in the detection of changes in gene expression. We found that genes regulating individual organ development are very rarely close to each other in the mouse and human genomes. We built an algorithm, called DELocal, to detect changes in gene expression that incorporates information about neighbouring genes. Using transcriptomes of developing mouse molar teeth containing gene expression profiles of thousands of genes, we show how genes regulating tooth development are ranked high by DELocal even if their expression level changes are subtle. We propose that developmental biology studies can benefit from gene neighbourhood analyses in the detection of differential expression and identification of organ specific genes.

Developmental influence on evolutionary rates and the origin of placental mammal tooth complexity

Interactions during development among genes, cells, and tissues can favor the more frequent generation of some trait variants compared with others. This developmental bias has often been considered to constrain adaptation, but its exact influence on evolution is poorly understood. Using computer simulations of development, we provide evidence that molecules promoting the formation of mammalian tooth cusps could help accelerate tooth complexity evolution. Only relatively small developmental changes were needed to derive the more complex, rectangular upper molar typical of early placental mammals from the simpler triangular ancestral pattern. Development may therefore have enabled the relatively fast divergence of the early placental molar dentition.

Development has often been viewed as a constraining force on morphological adaptation, but its precise influence, especially on evolutionary rates, is poorly understood. Placental mammals provide a classic example of adaptive radiation, but the debate around rate and drivers of early placental evolution remains contentious. A hallmark of early dental evolution in many placental lineages was a transition from a triangular upper molar to a more complex upper molar with a rectangular cusp pattern better specialized for crushing. To examine how development influenced this transition, we simulated dental evolution on “landscapes” built from different parameters of a computational model of tooth morphogenesis. Among the parameters examined, we find that increases in the number of enamel knots, the developmental precursors of the tooth cusps, were primarily influenced by increased self-regulation of the molecular activator (activation), whereas the pattern of knots resulted from changes in both activation and biases in tooth bud growth. In simulations, increased activation facilitated accelerated evolutionary increases in knot number, creating a lateral knot arrangement that evolved at least ten times on placental upper molars. Relatively small increases in activation, superimposed on an ancestral tritubercular molar growth pattern, could recreate key changes leading to a rectangular upper molar cusp pattern. Tinkering with tooth bud geometry varied the way cusps initiated along the posterolingual molar margin, suggesting that small spatial variations in ancestral molar growth may have influenced how placental lineages acquired a hypocone cusp. We suggest that development could have enabled relatively fast higher-level divergence of the placental molar dentition.

Mapping molar shapes on signaling pathways

A major challenge in evolutionary developmental biology is to understand how genetic mutations underlie phenotypic changes. In principle, selective pressures on the phenotype screen the gene pool of the population. Teeth are an excellent model for understanding evolutionary changes in the genotype-phenotype relationship since they exist throughout vertebrates. Genetically modified mice (mutants) with abnormalities in teeth have been used to explore tooth development. The relationship between signaling pathways and molar shape, however, remains elusive due to the high intrinsic complexity of tooth crowns. This hampers our understanding of the extent to which developmental factors explored in mutants explain developmental and phenotypic variation in natural species that represent the consequence of natural selection. Here we combine a novel morphometric method with two kinds of data mining techniques to extract data sets from the three-dimensional surface models of lower first molars: i) machine learning to maximize classification accuracy of 22 mutants, and ii) phylogenetic signal for 31 Murinae species. Major shape variation among mutants is explained by the number of cusps and cusp distribution on a tooth crown. The distribution of mutant mice in morphospace suggests a nonlinear relationship between the signaling pathways and molar shape variation. Comparative analysis of mutants and wild murines reveals that mutant variation overlaps naturally occurring diversity, including more ancestral and derived morphologies. However, taxa with transverse lophs are not fully covered by mutant variation, suggesting experimentally unexplored developmental factors in the evolutionary radiation of Murines.

System-level analyses of keystone genes required for mammalian tooth development

When a null mutation of a gene causes a complete developmental arrest, the gene is typically considered essential for life. Yet, in most cases, null mutations have more subtle effects on the phenotype. Here we used the phenotypic severity of mutations as a tool to examine system‐level dynamics of gene expression. We classify genes required for the normal development of the mouse molar into different categories that range from essential to subtle modification of the phenotype. Collectively, we call these the developmental keystone genes. Transcriptome profiling using microarray and RNAseq analyses of patterning stage mouse molars show highly elevated expression levels for genes essential for the progression of tooth development, a result reminiscent of essential genes in single‐cell organisms. Elevated expression levels of progression genes were also detected in developing rat molars, suggesting evolutionary conservation of this system‐level dynamics. Single‐cell RNAseq analyses of developing mouse molars reveal that even though the size of the expression domain, measured in the number of cells, is the main driver of organ‐level expression, progression genes show high cell‐level transcript abundances. Progression genes are also upregulated within their pathways, which themselves are highly expressed. In contrast, a high proportion of the genes required for normal tooth patterning are secreted ligands that are expressed in fewer cells than their receptors and intracellular components. Overall, even though expression patterns of individual genes can be highly different, conserved system‐level principles of gene expression can be detected using phenotypically defined gene categories.

The phenotypic severity of mutations on mouse teeth is used to classify genes.

Genes essential for the progression of odontogenesis are highly expressed at the organ and cell level.

Many of the genes required for normal patterning are locally expressed ligands.