Candidate gene analysis of tooth agenesis identifies novel mutations in six genes and suggests significant role for WNT and EDA signaling and allele combinations

Orodental malformations associated with human MSX1 sequence variants

BACKGROUND

MSX1 sequence variants have been known to cause human tooth agenesis (TA) with or without orofacial clefts. However, their roles during the whole processes of tooth development are not fully understood. This study aimed to characterize a 4-membered family with TA carrying a novel MSX1 pathogenic variant and investigate the disease mechanism.

METHODS

The authors conducted whole exome analysis to define the disease-causing sequence variant. They performed microcomputed tomography, morphometric analyses, transcriptome profiling, and molecular characterization to study the affected teeth and the gene variant.

RESULTS

The authors identified an MSX1 pathogenic variant, p.Glu232∗, in affected family members with TA and concomitant orodental anomalies, namely, prominent maxillary labial frenum, central incisor diastema, median maxillary anterior alveolar cleft, tooth fusion, mandibular molar dysmorphology, thin dentin layer, and slender dental roots. MSX1-defective teeth were not apparently microdontic but had thin dentin layers. The mandibular molars showed a homeotic transformation to maxillary counterparts. Genes involved in extracellular matrix organization and dentinogenesis, such as DMP1 and MMP20, were downregulated in dental pulp tissues of MSX1-defective teeth. The p.Glu232∗-truncated MSX1 properly localized to the nucleus but partially lost its transactivation ability. Analyzing reported cases indicated that truncation sequence variants within the homeobox domain of MSX1 caused a more severe TA phenotype than those outside of the homeobox domain, probably due to dominant negativity compared with haploinsufficiency.

CONCLUSIONS

This study provides in vivo evidence that MSX1 contributes to developmental processes of various orodental tissues in humans.

PRACTICAL IMPLICATIONS

Clinically, hypertrophic labial frenum, incisor diastema, and median maxillary anterior alveolar cleft might be considered diagnostic for MSX1-associated TA.

Digenic inheritance accounts for phenotypic variability in amelogenesis imperfecta

Amelogenesis imperfecta (AI) represents a group of clinically and genetically heterogeneous disorders that affect enamel formation and mineralization. Although AI is commonly considered a monogenic disorder, digenic inheritance is rarely reported. In this study, we recruited two nonconsanguineous Chinese families exhibiting diverse phenotypes of enamel defects among affected family members. Digenic variants were discovered in both probands. In family 1, the proband inherited a paternal frameshift variant in LAMA3 (NM_198129.4:c.3712dup) and a maternal deletion encompassing the entire AMELX gene. This resulted in a combined hypoplastic and hypomineralized AI phenotype, which was distinct from the parents' manifestations. In family 2, whole-exome sequencing analysis revealed the proband carried a maternal heterozygous splicing variant in COL17A1 (NC_000010.11 (NM_000494.3): c.4156 + 2dup) and compound heterozygous variants in RELT (paternal: NM_032871.4:c.260A > T; maternal: NM_032871.4:c.521 T > G). These genetic changes caused the abundant irregular enamel defects observed in the proband, whereas other affected family members carrying heterozygous variants in both COL17A1 and RELT displayed only horizontal grooves as their phenotype. The pathogenicity of the novel COL17A1 splice site variant was confirmed through RT-PCR and minigene assay. This study enhances our understanding by highlighting the potential association between the co-occurrence of variants in two genes and variable phenotypes observed in AI patients.

[The tooth: A marker of developmental abnormalities]

Tooth formation results from specific epithelial-mesenchymal interactions, which summarize a number of developmental processes. Tooth anomalies may thus reflect subclinical diseases of the kidney, bone and more broadly of the mineral metabolism, skin or nervous system. Odontogenesis starts from the 3rd week of intrauterine life by the odontogenic orientation of epithelial cells by a first PITX2 signal. The second phase is the acquisition of the number, shape, and position of teeth. It depends on multiple transcription and growth factors (BMP, FGF, SHH, WNT). These ecto-mesenchymal interactions guide cell migration, proliferation, apoptosis and differentiation ending in the formation of the specific dental mineralized tissues. Thus, any alteration will have consequences on the tooth structure or shape. Resulting manifestations will have to be considered in the patient phenotype and the multidisciplinary care, but also may contribute to identify the altered genetic circuity.

Characterization of novel MSX1 variants causally associated with non-syndromic oligodontia in Chinese families

BACKGROUND

MSX1 (OMIM #142983) is crucial to normal dental development, and variants in MSX1 are associated with dental anomalies. The objective of this study was to characterize the pathogenicity of novel MSX1 variants in Chinese families with non-syndromic oligodontia (NSO).

METHODS

Genomic DNA was extracted from individuals representing 35 families with non-syndromic oligodontia and was analyzed by Sanger sequencing and whole-exome sequencing. Pathogenic variants were screened via analyses involving PolyPhen-2, Sorting-Intolerant from Tolerant, and MutationTaster, and conservative analysis of variants. Patterns of MSX1-related NSO were analyzed. MSX1 structural changes suggested functional consequences in vitro.

RESULTS

Three previously unreported MSX1 heterozygous variants were identified: one insertion variant (c.576_577insTAG; p.Gln193*) and two missense variants (c. 871T>C; p.Tyr291His and c. 644A>C; p.Gln215Pro). Immunofluorescence analysis revealed abnormal subcellular localization of the p.Gln193* MSX1 variant. In addition, we found that these MSX1 variants likely lead to the loss of second premolars.

CONCLUSION

Three novel MSX1 variants were identified in Chinese Han families with NSO, expanding the MSX1 variant spectrum and presenting a genetic origin for the pathogenesis detected in patients and their families.

A Missense Mutation in the Collagen Triple Helix of EDA Is Associated with X-Linked Recessive Hypohidrotic Ectodermal Dysplasia in Fleckvieh Cattle

Mutations within the ectodysplasin A (EDA) gene have been associated with congenital hypotrichosis and anodontia (HAD/XHED) in humans, mice, dogs and cattle. We identified a three-generation family of Fleckvieh cattle with male calves exhibiting clinical and histopathological signs consistent with an X-linked recessive HAD (XHED). Whole genome and Sanger sequencing of cDNA showed a perfect association of the missense mutation g.85716041G>A (ss2019497443, rs1114816375) within the EDA gene with all three cases following an X-linked recessive inheritance, but normal EDAR and EDARADD. This mutation causes an exchange of glycine (G) with arginine (R) at amino acid position 227 (p.227G>R) in the second collagen triple helix repeat domain of EDA. The EDA variant was associated with a significant reduction and underdevelopment of hair follicles along with a reduced outgrowth of hairs, a complete loss of seromucous nasolabial and mucous tracheal and bronchial glands and a malformation of and reduction in number of teeth. Thermostability of EDA G227R was reduced, consistent with a relatively mild hair and tooth phenotype. However, incisors and canines were more severely affected in one of the calves, which correlated with the presence of a homozygous missense mutation of RNF111 (g.51306765T>G), a putative candidate gene possibly associated with tooth number in EDA-deficient Fleckvieh calves.

Phenotypic characteristics of taurodontism and a novel WNT10A variant in non-syndromic oligodontia family

OBJECTIVE

Variants in wingless-type MMTV integration site family member 10A (WNT10A) have been proposed to be the most common cause of non-syndromic oligodontia (NSO). The goal of the present study was to identify the novel WNT10A variants in Chinese families with NSO.

DESIGN

Clinical data were collected from 39 families with oligodontia admitted to the Hospital of Stomatology Hebei Medical University (China) from 2016 to 2022. Whole-exome sequencing (WES) and Sanger sequencing were performed to identify WNT10A variants in three families with non-syndromic oligodontia. Amino acid conservation analysis and protein conformational analysis were conducted for the WNT10A variant. Genotype-phenotype analysis was performed on the previously reported WNT10A variants related to NSO.

RESULTS

We found a novel heterozygous WNT10A variant c.1127 G>A (p.Cys376Tyr) and two reported heterozygous variants c.460 C>A (p.Leu154Met) and c.511 C>T (p.Arg171Cys). Structural modeling showed that the novel WNT10A variant was located in a highly conserved domain, which led to structural damage of WNT10A protein. In addition, we found that the phenotype of the WNT10A variants affected the maxillary second premolars, followed by the mandibular second premolars, and rarely affected the maxillary central incisor. Herein, it is the first time to report that NSO patients with WNT10A monoallele mutation carry taurodontism phenotype and 6.1% prevalence of taurodontism in WNT10A-related NSO patients.

CONCLUSIONS

Our results demonstrated that the novel variant c.1127 G>A (p.Cys376Tyr) of WNT10A causes NSO. The present study expanded the known variation spectrum of WNT10A and provided valuable information for genetic counseling of families.

Identification of potential key variants in mandibular premolar hypodontia through whole-exome sequencing

Determining genotype-phenotype correlations in patients with hypodontia is important for understanding disease pathogenesis, although only a few studies have elucidated it. We aimed to identify genetic variants linked to non-syndromic bilateral mandibular second premolar hypodontia in a Korean population for the first time by specifying the phenotype of hypodontia. Twenty unrelated individuals with non-syndromic bilateral mandibular second premolar hypodontia were enrolled for whole-exome sequencing. Using a tooth agenesis gene set panel consisting of 112 genes based on literature, potential candidate variants were screened through variant filtering and prioritization. We identified 13 candidate variants in 12 genes, including a stop-gain variant (c.4750C>T) in LAMA3. Through the functional enrichment analysis of the prioritized genes, several terms related to tooth development were enriched in a protein-protein interaction network of candidate genes for mandibular premolar hypodontia. The hypodontia group also had approximately 2-fold as many mutated variants in all four genes related to these key terms, which are CDH1, ITGB4, LAMA3, LAMB3, as those in the 100 healthy control group individuals. The relationship between enriched terms and pathways and mandibular premolar hypodontia was also investigated. In addition, we identified some known oligodontia variants in patients with hypodontia, strengthening the possibility of synergistic effects in other genes. This genetic investigation may be a worthwhile preliminary attempt to reveal the pathogenesis of tooth agenesis and sets a background for future studies.

PAX9 mutations and genetic synergism in familial tooth agenesis

Familial tooth agenesis (FTA) is one of the most common craniofacial anomalies in humans. Loss-of-function mutations in PAX9 and WNT10A have been known to cause FTA with various expressivity. In this study, we identified five FTA kindreds with novel PAX9 disease-causing mutations: p.(Glu7Lys), p.(Val83Leu), p.(Pro118Ser), p.(Ser197Argfs*23), and c.771+4A>G. Concomitant PAX9 and WNT10A pathogenic variants found in two probands with severe phenotypes suggested an effect of mutational synergism. All overexpressed PAX9s showed proper nuclear localization, excepting the p.(Pro118Ser) mutant. Various missense mutations caused differential loss of PAX9 transcriptional ability. PAX9 overexpression in dental pulp cells upregulated LEF1 and AXIN2 expression, indicating a positive regulatory role for PAX9 in canonical Wnt signaling. Analyzing 176 cases with 63 different mutations, we observed a distinct pattern of tooth agenesis for PAX9-associated FTA: Maxillary teeth are in general more frequently affected than mandibular ones. Along with all second molars, maxillary bicuspids and first molars are mostly involved, while maxillary lateral incisors and mandibular bicuspids are relatively less affected. Genotypically, missense mutations are associated with fewer missing teeth than frameshift and nonsense variants. This study significantly expands the phenotypic and genotypic spectrums of PAX9-associated disorders and reveals a molecular mechanism of genetic synergism underlying FTA variable expressivity.

The EDA/EDAR/NF-κB pathway in non-syndromic tooth agenesis: A genetic perspective

Non-syndromic tooth agenesis (NSTA) is one of the most common dental developmental malformations affected by genetic factors predominantly. Among all 36 candidate genes reported in NSTA individuals, EDA, EDAR, and EDARADD play essential roles in ectodermal organ development. As members of the EDA/EDAR/NF-κB signaling pathway, mutations in these genes have been implicated in the pathogenesis of NSTA, as well as hypohidrotic ectodermal dysplasia (HED), a rare genetic disorder that affects multiple ectodermal structures, including teeth. This review provides an overview of the current knowledge on the genetic basis of NSTA, with a focus on the pathogenic effects of the EDA/EDAR/NF-κB signaling pathway and the role of EDA, EDAR, and EDARADD mutations in developmental tooth defects. We also discuss the phenotypic overlap and genetic differences between NSTA and HED. Ultimately, this review highlights the importance of genetic analysis in diagnosing and managing NSTA and related ectodermal disorders, and the need for ongoing research to improve our understanding of these conditions.

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.

Genotype-phenotype pattern analysis of pathogenic PAX9 variants in Chinese Han families with non-syndromic oligodontia

Background: Non-syndromic oligodontia is characterized by the absence of six or more permanent teeth, excluding third molars, and can have aesthetic, masticatory, and psychological consequences. Previous studies have shown that PAX9 is associated with autosomal dominant forms of oligodontia but the precise molecular mechanisms are still unknown.

Methods: Whole-exome and Sanger sequencing were performed on a cohort of approximately 28 probands with NSO, for mutation analysis. Bioinformatic analysis was performed on the potential variants. Immunofluorescence assay, western blotting, and qPCR were used to explore the preliminary functional impact of the variant PAX9 proteins. We reviewed PAX9-related NSO articles in PubMed to analyze the genotype-phenotype correlations.

Results: We identified three novel PAX9 variants in Chinese Han families: c.152G>T (p.Gly51Val), c.239delC (p.Thr82Profs*3), and c.409C>T (q.Gln137Ter). In addition, two previously reported missense variants were identified: c.140G>C (p.Arg47Pro) and c.146C>T (p.Ser49Leu) (reference sequence NM_006194.4). Structural modeling revealed that all missense variants were located in the highly conserved paired domain. The other variants led to premature termination of the protein, causing structural impairment of the PAX9 protein. Immunofluorescence assay showed abnormal subcellular localizations of the missense variants (R47P, S49L, and G51V). In human dental pulp stem cells, western blotting and qPCR showed decreased expression of PAX9 variants (c.140G>C, p.R47P, and c.152G>T, p.G51V) compared with the wild-type group at both the transcription and translation levels. A review of published papers identified 64 PAX9 variants related to NSO and found that the most dominant feature was the high incidence of missing upper second molars, first molars, second premolars, and lower second molars.

Conclusion: Three novel PAX9 variants were identified in Chinese Han families with NSO. These results extend the variant spectrum of PAX9 and provide a foundation for genetic diagnosis and counseling.

Next generation sequencing panel target genes: possible diagnostic tool for ectodermal dysplasia related diseases

BACKGROUND

Ectodermal dysplasias (EDs) are a large and complex group of disorders affecting the ectoderm-derived organs; the clinical and genetic heterogeneity of these conditions renders an accurate diagnosis more challenging. The aim of this study is to demonstrate the clinical utility of a targeted resequencing panel through enhancing the molecular and clinical diagnosis of EDs. Given the recent developments in gene and protein-based therapies for X-linked hypohidrotic ectodermal dysplasia, there is a re-emerging interest in identifying the genetic basis of EDs and the respective phenotypic presentations, in an aim to facilitate potential treatments for affected families.

METHODS

We assessed seventeen individuals, from three unrelated families, who presented with diverse phenotypes suggestive of ED. An extensive multidisciplinary clinical evaluation was performed followed by a targeted exome resequencing panel (including genes that are known to cause EDs). MiSeqTM data software was used, variants with Qscore >30 were accepted.

RESULTS

Three different previously reported hemizygous EDA mutations were found in the families. However, a complete genotype-phenotype correlation could not be established, neither in our patients nor in the previously reported patients.

CONCLUSIONS

Targeted exome resequencing can provide a rapid and accurate diagnosis of EDs, while further contributing to the existing ED genetic data. Moreover, the identification of the disease-causing mutation in an affected family is crucial for proper genetic counseling and the establishment of a genotype-phenotype correlation which will subsequently provide the affected individuals with a more suitable treatment plan.

A novel EDAR variant identified in non-syndromic tooth agenesis: Insights from molecular dynamics

OBJECTIVE

This study aims to investigate a novel pathogenic variant in a Chinese family of non-syndromic tooth agenesis (NSTA) and study the impact of the variant on related protein and pathway.

DESIGN

One NSTA family was collected. Whole exome sequencing and Sanger sequencing were performed on the proband with NSTA and his 5 family members. The pathogenic influence of the mutant is evaluated by bioinformatics analyses including evolutionary conservation analysis and secondary structure prediction. Molecular dynamics (MD) simulations and binding free energy calculations were then performed to explore changes in the tertiary structure and binding ability of the protein.

RESULTS

We found a novel missense ectodysplasin A receptor (EDAR) variant (c .1292 T > G; p.Ile431Arg) in all affected family members. The results of bioinformatics analyses revealed that the EDAR had harmful changes after mutation. MD simulations and the binding free energy calculations results showed that the mutant EDAR protein and EDAR/ectodysplasin-A receptor-associated adapter (EDARADD) complex displayed tertiary structural change, and EDAR possessed a lower affinity to EDARADD after mutation.

CONCLUSIONS

We found a novel EDAR variant (c.1292 T > G; p.Ile431Arg) in one NSTA family, which affects the binding of EDAR and EDARADD.

Eda controls the size of the enamel knot during incisor development

Ectodysplasin (Eda) plays important roles in both shaping the developing tooth and establishing the number of teeth within the tooth row. Sonic hedgehog (Shh) has been shown to act downstream of Eda and is involved in the initiation of tooth development. Eda-/- mice possess hypoplastic and hypomineralized incisors and show changes in tooth number in the molar region. In the present study we used 3D reconstruction combined with expression analysis, cell lineage tracing experiments, and western blot analysis in order to investigate the formation of the incisor germs in Eda-/- mice. We show that a lack of functional Eda protein during early stages of incisor tooth germ development had minimal impact on development of the early expression of Shh in the incisor, a region proposed to mark formation of a rudimental incisor placode and act as an initiating signalling centre. In contrast, deficiency of Eda protein had a later impact on expression of Shh in the primary enamel knot of the functional tooth. Eda-/- mice had a smaller region where Shh was expressed, and a reduced contribution from Shh descendant cells. The reduction in the enamel knot led to the formation of an abnormal enamel organ creating a hypoplastic functional incisor. Eda therefore appears to influence the spatial formation of the successional signalling centres during odontogenesis.

Tooth number abnormality: from bench to bedside

Tooth number abnormality is one of the most common dental developmental diseases, which includes both tooth agenesis and supernumerary teeth. Tooth development is regulated by numerous developmental signals, such as the well-known Wnt, BMP, FGF, Shh and Eda pathways, which mediate the ongoing complex interactions between epithelium and mesenchyme. Abnormal expression of these crutial signalling during this process may eventually lead to the development of anomalies in tooth number; however, the underlying mechanisms remain elusive. In this review, we summarized the major process of tooth development, the latest progress of mechanism studies and newly reported clinical investigations of tooth number abnormality. In addition, potential treatment approaches for tooth number abnormality based on developmental biology are also discussed. This review not only provides a reference for the diagnosis and treatment of tooth number abnormality in clinical practice but also facilitates the translation of basic research to the clinical application.

Novel PAX9 compound heterozygous variants in a Chinese family with non-syndromic oligodontia and genotype-phenotype analysis of PAX9 variants

OBJECTIVE

Studies have reported that >91.9% of non-syndromic tooth agenesis cases are caused by seven pathogenic genes. To report novel heterozygous PAX9 variants in a Chinese family with non-syndromic oligodontia and summarize the reported genotype-phenotype relationship of PAX9 variants.

METHODOLOGY

We recruited 28 patients with non-syndromic oligodontia who were admitted to the Hospital of Stomatology Hebei Medical University (China) from 2018 to 2021. Peripheral blood was collected from the probands and their core family members for whole-exome sequencing (WES) and variants were verified by Sanger sequencing. Bioinformatics tools were used to predict the pathogenicity of the variants. SWISS-MODEL homology modeling was used to analyze the three-dimensional structural changes of variant proteins. We also analyzed the genotype-phenotype relationships of PAX9 variants.

RESULTS

We identified novel compound heterozygous PAX9 variants (reference sequence NM_001372076.1) in a Chinese family with non-syndromic oligodontia: a new missense variant c.1010C>A (p.T337K) in exon 4 and a new frameshift variant c.330_331insGT (p.D113Afs*9) in exon 2, which was identified as the pathogenic variant in this family. This discovery expands the known variant spectrum of PAX9; then, we summarized the phenotypes of non-syndromic oligodontia with PAX9 variants.

CONCLUSION

We found that PAX9 variants commonly lead to loss of the second molars.

An in vitro and computational validation of a novel loss-of-functional mutation in PAX9 associated with non-syndromic tooth agenesis

Congenital tooth agenesis (CTA) is one of the most common craniofacial anomalies. Its frequency varies among different population depending upon the genetic heterogeneity. CTA could be of familial or sporadic and syndromic or non-syndromic. Five major genes are found to be associated with non-syndromic CTA, namely PAX9, MSX1, EDA1, AXIN2, and WNT10A. Very few studies have been carried out so far on CTA on this Indian population making this study unique and important. This study was initiated to identify potential pathogenic variant associated with congenital tooth agenesis in an India family with molar tooth agenesis. CTA was investigated and a novel c.336C > G variation was identified in the exon 3 of PAX9, leading to substitution of evolutionary conserved Cys with Trp at 112th amino acid position located at the functionally significant DNA-binding paired domain region. Functional analysis revealed that p.Cys112Trp mutation did not prevent the nuclear localization although mutant protein had higher cytoplasmic retention. EMSA using e5 probe revealed that mutant protein was unable to bind with the paired-domain-binding site. Subsequently, GST pull-down assay revealed lower binding activity of the mutant protein with its known interactor MSX1. These in vitro results were consistent with the computational results. The in vitro and computational observations altogether suggest that c.336C > G (p.Cys112Trp) variation leads to loss of function of PAX9 leading to CTA in this family.

Rethinking the Genetic Etiology of Nonsyndromic Tooth Agenesis

PURPOSE OF REVIEW

Genetic studies in humans and animal models have improved our understanding of the role of numerous genes in the etiology of nonsyndromic tooth agenesis (TA). The purpose of this review is to discuss recently identified genes potentially contributing to TA.

RECENT FINDINGS

Despite research progress, understanding the genetic factors underlying nonsyndromic TA has been challenging given the genetic heterogeneity, variable expressivity, and incomplete penetrance of putatively pathogenic variants often observed associated with the condition. Next-generation sequencing technologies have provided a platform for novel gene and variant discoveries and informed paradigm-shifting concepts in the etiology of TA. This review summarizes the current knowledge on genes and pathways related to nonsyndromic TA with a focus on recently identified genes/variants. Evidence suggesting possible multi-locus variation in TA is also presented.

WNT10A and RUNX2 mutations associated with non-syndromic tooth agenesis

The goal of this study was to examine the prevalence of WNT10A and RUNX2 mutations and assess their potential impact on the phenotype of non-syndromic tooth agenesis. The study included 30 participants with non-syndromic tooth agenesis, divided into hypodontia (n = 24) and oligodontia forms (n = 6), and 42 unaffected family members. Genomic DNA from buccal epithelial cells was used for polymerase chain reaction amplification of functionally important exons of the WNT10A and RUNX2 genes. Direct sequencing reactions were performed to confirm the presence of mutations. The trend of increasing prevalence of WNT10A mutations and a slight increase in the prevalence of RUNX2 mutations were revealed in tooth agenesis cases compared to unaffected family members. There was a higher prevalence of hypodontia than oligodontia, increased frequency of females over males with missing teeth, and a wide phenotypic variability was observed in individuals and families analyzed. The common missense mutations (p.Phe228Ile, p.Arg113Cys, p.Asp217Asn, and p.Gly165Arg) and c.114-56T>C in the WNT10A gene and in-frame-deletion/insertions (11A, 24Q, 30Q), synonymous variant c.240G>A, and 424-33dupC in the RUNX2 gene were identified. These findings highlight an important role of WNT10A and RUNX2 mutations in the genetic etiology of non-syndromic tooth agenesis.

Novel Candidate Genes for Non-Syndromic Tooth Agenesis Identified Using Targeted Next-Generation Sequencing

Non-syndromic tooth agenesis (ns-TA) is one of the most common dental anomalies characterized by the congenital absence of at least one permanent tooth (excluding third molars). Regarding the essential role of genetic factors in ns-TA aetiology, the present study aimed to identify novel pathogenic variants underlying hypodontia and oligodontia. In a group of 65 ns-TA patients and 127 healthy individuals from the genetically homogenous Polish population, the coding sequences of 423 candidate genes were screened using targeted next-generation sequencing. Pathogenic and likely pathogenic variants were identified in 37 (56.92%) patients, including eight nucleotide alternations of genes not previously implicated in ns-TA (CHD7, CREBBP, EVC, LEF1, ROR2, TBX22 and TP63). However, since only single variants were detected, future research is required to confirm and fully understand their role in the aetiology of ns-TA. Additionally, our results support the importance of already known ns-TA candidate genes (AXIN2, EDA, EDAR, IRF6, LAMA3, LRP6, MSX1, PAX9 and WNT10A) and provide additional evidence that ns-TA might be an oligogenic condition involving the cumulative effect of rare variants in two or more distinct genes.

WNT10A variants: following the pattern of inheritance in tooth agenesis and self-reported family history of cancer

OBJECTIVES

The aim of this study was the analysis of WNT10A variants in seven families of probands with various forms of tooth agenesis and self-reported family history of cancer.

MATERIALS AND METHODS

We enrolled 60 young subjects (aged 13 to 17) from the Czech Republic with various forms of tooth agenesis. Dental phenotypes were assessed using Planmeca ProMax 3D (Planmeca Oy, Finland) with Planmeca Romexis software (version 2.9.2) together with oral examinations. After screening PAX9, MSX1, EDA, EDAR, AXIN2 and WNT10A genes on the Illumina MiSeq platform (Illumina, USA), we further analyzed the evolutionarily highly conserved WNT10A gene by capillary sequencing in the seven families.

RESULTS

All the detected variants were heterozygous or compound heterozygous with various levels of phenotypic expression. The most severe phenotype (oligodontia) was found in a proband who was compound heterozygous for the previously identified WNT10A variant p.Phe228Ile and a newly discovered c.748G > A variant (p.Gly250Arg) of WNT10A. The newly identified variant causes substitution of hydrophobic glycine for hydrophilic arginine.

CONCLUSIONS

We suggest that the amino acid changes in otherwise highly conserved sequences significantly affect the dental phenotype. No relationship between the presence of WNT10A variants and a risk of cancer has been found.

CLINICAL RELEVANCE

Screening of PAX9, MSX1, EDA, EDAR, AXIN2 and WNT10A genes in hope to elucidate the pattern of inheritance in families.

An Educational Assessment of Evidence Used for Variant Classification: A Report of the Association for Molecular Pathology

The Association for Molecular Pathology Variant Interpretation Testing Among Laboratories (VITAL) Working Group convened to evaluate the Standards and Guidelines for the Interpretation of Sequence Variants implementation into clinical practice, identify problematic classification rules, and define implementation challenges. Variants and associated clinical information were provided to volunteer respondents. Participant variant classifications were compared with intended consensus-derived classifications of the Working Group. The 24 variant challenges received 1379 responses; 1119 agreed with the intended response (81%; 95% CI, 79% to 83%). Agreement ranged from 44% to 100%, with 16 challenges (67%; 47% to 82%) reaching consensus (≥80% agreement). Participant classifications were also compared to a calculated interpretation of the ACMG Guidelines using the participant-reported criteria as input. The 24 variant challenges had 1368 responses with specific evidence provided and 1121 (82%; 80% to 84%) agreed with the calculated interpretation. Agreement for challenges ranged from 63% to 98%; 15 (63%; 43% to 79%) reaching consensus. Among 81 individual participants, 32 (40%; 30% to 50%) reached agreement with at least 80% of the intended classifications and 42 (52%; 41% to 62%) with the calculated classifications. This study demonstrated that although variant classification remains challenging, published guidelines are being utilized and adapted to improve variant calling consensus. This study identified situations where clarifications are warranted and provides a model for competency assessment.

Patterns of molar agenesis associated with p.P20L and p.R77Q variants in PAX9

Nonsyndromic tooth agenesis is associated with variants in several genes. There are numerous genotype-phenotype publications involving many patients and kindreds. Here, we identified six Thai individuals in two families with nonsyndromic tooth agenesis, performed exome sequencing, and conducted functional experiments. Family 1 had four affected members carrying the heterozygous PAX9 variant, c.59C>T (p.Pro20Leu). The p.Pro20Leu was previously reported in two families having four and three affected members. These seven cases and Proband-1 had agenesis of at least three third molars. Family 2 comprised two affected members with agenesis of all 12 molars. Both individuals were heterozygous for c.230G>A (p.Arg77Gln) in PAX9, which has not been reported previously. This variant is predicted to be damaging, evolutionarily conserved, and resides in the PAX9 linking peptide. The BMP4 RNA levels in Proband-1's leukocytes were not significantly different from those in the controls, whereas BMP4 levels observed in Proband-2 were significantly increased. Moreover, the p.Arg77Gln variant demonstrated nuclear localization similar to the wild-type but resulted in significantly impaired transactivation of BMP4, a PAX9 downstream gene. In conclusion, we demonstrate that the PAX9 p.Pro20Leu is highly associated with absent third molars, while the novel PAX9 p.Arg77Gln impairs BMP4 transactivation and is associated with total molar agenesis.

Ectodysplasin pathogenic variants affecting the furin-cleavage site and unusual clinical features define X-linked hypohidrotic ectodermal dysplasia in India

Hypohidrotic ectodermal dysplasia (HED) is a rare genetic disorder caused by mutational inactivation of a developmental pathway responsible for generation of tissues of ectodermal origin. The X-linked form accounts for the majority of HED cases and is caused by Ectodysplasin (EDA) pathogenic variants. We performed a combined analysis of 29 X-linked hypohidrotic ectodermal dysplasia (XLHED) families (including 12 from our previous studies). In addition to the classical triad of symptoms including loss (or reduction) of ectodermal structures, such as hair, teeth, and sweat glands, we detected additional HED-related clinical features including facial dysmorphism and hyperpigmentation in several patients. Interestingly, global developmental delay was identified as an unusual clinical symptom in many patients. More importantly, we identified 22 causal pathogenic variants that included 15 missense, four small in-dels, and one nonsense, splice site, and large deletion each. Interestingly, we detected 12 unique (India-specific) pathogenic variants. Of the 29 XLHED families analyzed, 11 (38%) harbored pathogenic variant localized to the furin cleavage site. A comparison with HGMD revealed significant differences in the frequency of missense pathogenic variants; involvement of specific exons and/or protein domains and transition/transversion ratios. A significantly higher proportion of missense pathogenic variants (33%) localized to the EDA furin cleavage when compared to HGMD (7%), of which p.R155C, p.R156C, and p.R156H were detected in three families each. Therefore, the first comprehensive analysis of XLHED from India has revealed several unique features including unusual clinical symptoms and high frequency of furin cleavage site pathogenic variants.

Intertwined Signaling Pathways Governing Tooth Development: A Give-and-Take Between Canonical Wnt and Shh

Teeth play essential roles in life. Their development relies on reciprocal interactions between the ectoderm-derived dental epithelium and the underlying neural crest-originated mesenchyme. This odontogenic process serves as a prototype model for the development of ectodermal appendages. In the mouse, developing teeth go through distinct morphological phases that are tightly controlled by epithelial signaling centers. Crucial molecular regulators of odontogenesis include the evolutionarily conserved Wnt, BMP, FGF and sonic hedgehog (Shh) pathways. These signaling modules do not act on their own, but are closely intertwined during tooth development, thereby outlining the path to be taken by specific cell populations including the resident dental stem cells. Recently, pivotal Wnt-Shh interaction and feedback loops have been uncovered during odontogenesis, showing conservation in other developing ectodermal appendages. This review provides an integrated overview of the interplay between canonical Wnt and Shh throughout mouse tooth formation stages, extending from the initiation of dental placode to the fully formed adult tooth.

Analyses of oligodontia phenotypes and genetic etiologies

Oligodontia is the congenital absence of six or more teeth and comprises the more severe forms of tooth agenesis. Many genes have been implicated in the etiology of tooth agenesis, which is highly variable in its clinical presentation. The purpose of this study was to identify associations between genetic mutations and clinical features of oligodontia patients. An online systematic search of papers published from January 1992 to June 2021 identified 381 oligodontia cases meeting the eligibility criteria of causative gene mutation, phenotype description, and radiographic records. Additionally, ten families with oligodontia were recruited and their genetic etiologies were determined by whole-exome sequence analyses. We identified a novel mutation in WNT10A (c.99_105dup) and eight previously reported mutations in WNT10A (c.433 G > A; c.682 T > A; c.318 C > G; c.511.C > T; c.321 C > A), EDAR (c.581 C > T), and LRP6 (c.1003 C > T, c.2747 G > T). Collectively, 20 different causative genes were implicated among those 393 cases with oligodontia. For each causative gene, the mean number of missing teeth per case and the frequency of teeth missing at each position were calculated. Genotype-phenotype correlation analysis indicated that molars agenesis is more likely linked to PAX9 mutations, mandibular first premolar agenesis is least associated with PAX9 mutations. Mandibular incisors and maxillary lateral incisor agenesis are most closely linked to EDA mutations.

Association between tooth agenesis and cancer: a systematic review

The congenital absence of multiple teeth may share the same genetic background of the development of some types of cancer. Objective: This systematic review aimed to investigate the possible association between dental agenesis and cancer, and the perspective of agenesis as an early predictor for cancer risk. Methodology: The electronic databases PubMed, Scopus, Web of Science, Cochrane Library, LILACS, and OpenGrey were searched and the risk of bias was evaluated using the Newcastle-Ottawa tool. The GRADE tool was used to evaluate the certainty of the evidence. Results: Six studies met the eligibility criteria. A positive co-occurrence between ovarian cancer and hypodontia was found in two articles. Three studies evaluated the association between dental agenesis and colorectal cancer and only one showed common genes for these conditions. One paper found individuals with hypodontia had a higher risk of family history of cancer. Five studies had a fair quality and one a good quality. The certainty of evidence was classified as very low. Conclusion: Notwithstanding the limited scientific evidence, there may be a possible association between dental agenesis and cancer due to genes involved in both conditions. Agenesis of multiple teeth could be an early indicator of cancer risk. Nevertheless, studies with a better level of evidence are needed to confirm this possible association.

A novel EDAR missense mutation identified by whole-exome sequencing with non-syndromic tooth agenesis in a Chinese family

Background

Causative variants in genes of the EDA/EDAR/NF‐κB pathway, such as EDA and EDARADD, have been widely identified in patients with non‐syndromic tooth agenesis (NSTA). However, few cases of NSTA are due to ectodysplasin‐A receptor (EDAR) variants. In this study, we investigated NSTA‐associated variants in Chinese families.

Methods

Peripheral blood samples were collected from the family members of 24 individuals with NSTA for DNA extraction. The coding region of the EDA gene of the 24 probands was amplified by PCR and sequenced to investigate new variants. Whole‐exome sequencing and Sanger sequencing were then performed for probands without EDA variants detected by PCR.

Results

A novel missense variant EDAR c.338G>A (p.(Cys113Tyr)) was identified in one family. In addition, three known EDA variants (c.865C>T, c.866G>A, and c.1013C>T) were identified in three families. Genotype–phenotype correlation analysis of EDAR gene mutation showed that NSTA patients were most likely to lose the maxillary lateral incisors and the maxillary central incisors were the least affected. The phenotype of mutations at codon 289 of EDA in NSTA affected patients was characterized by lateral incisors loss, rarely affecting the maxillary first molars.

Conclusion

A novel EDAR missense variant c.338G>A (p.(Cys113Tyr)) was identified in a family with NSTA, extending the mutation spectrum of the EDAR gene. Genotype–phenotype correlation analyses of EDAR and EDA mutations could help to improve disease status prediction in NSTA families.

Novel MSX1 variants identified in families with nonsyndromic oligodontia

The goal of this study was to identify MSX1 gene variants in multiple Chinese families with nonsyndromic oligodontia and analyse the functional influence of these variants. Whole-exome sequencing (WES) and Sanger sequencing were performed to identify the causal gene variants in five families with nonsyndromic oligodontia, and a series of bioinformatics databases were used for variant confirmation and functional prediction. Phenotypic characterization of the members of these families was described, and an in vitro analysis was performed for functional evaluation. Five novel MSX1 heterozygous variants were identified: three missense variants [c.662A>C (p.Q221P), c.670C>T (p.R224C), and c.809C>T (p.S270L)], one nonsense variant [c.364G>T (p.G122*)], and one frameshift variant [c.277delG (p.A93Rfs*67)]. Preliminary in vitro studies demonstrated that the subcellular localization of MSX1 was abnormal with the p.Q221P, p.R224C, p.G122*, and p.A93Rfs*67 variants compared to the wild type. Three variants (p.Q221P, p.G122*, and p.A93Rfs*67) were classified as pathogenic or likely pathogenic, while p.S270L and p.R224C were of uncertain significance in the current data. Moreover, we summarized and analysed the MSX1-related tooth agenesis positions and found that the type and variant locus were not related to the severity of tooth loss. Our results expand the variant spectrum of nonsyndromic oligodontia and provide valuable information for genetic counselling.

Phenotypic and Genotypic Features of Thai Patients With Nonsyndromic Tooth Agenesis and WNT10A Variants

Tooth agenesis is one of the most common orodental anomalies that demonstrate phenotypic and genotypic heterogeneity with a prevalence of 2.5%–7%. Mutations in WNT10A have been proposed to be the most common cause of nonsyndromic tooth agenesis (NSTA). The aim of this study was to characterize the dental features and genetic variants of NSTA in a Thai population. We recruited 13 unrelated patients with NSTA who attended the Faculty of Dentistry, Chulalongkorn University, Thailand, from 2017 to 2019. All 13 underwent whole exome sequencing that identified likely pathogenic genetic variants, all in WNT10A, in five patients. All five patients had second premolar agenesis, while three also had absent or peg-shaped upper lateral incisors. Patient 1 possessed a novel heterozygous duplication, c.916_918dupAAC (p.Asn306dup) in WNT10A. Patients 2 and 3 harbored a heterozygous and homozygous c.637G > A (p.Gly213Ser) in WNT10A, respectively. Patients 4 possessed a heterozygous c.511C > T (p.Arg171Cys) in WNT10A. Patient 5 harbored a homozygous c.511C > T (p.Arg171Cys) in WNT10A and a novel heterozygous c.413A > T (p.Asn138Ile) in EDARADD, suggesting digenic inheritance. We recruited another 18 family members of these five patients. Out of 23 participants, homozygous WNT10A variants were identified in 2 patients and heterozygous variants in 17 individuals. Both homozygous patients had NSTA. Eight out of 17 heterozygous individuals (8/17) had NSTA or a peg-shaped lateral incisor, indicating a 47% penetrance of the heterozygous variants or 53% (10/19) penetrance of either homozygous or heterozygous variants in WNT10A. The frequencies of the c.511C > T in our in-house 1,876 Thai exome database, Asian populations, and non-Asian populations were 0.016, 0.005–0.033, and 0.001, respectively; while those of the c.637G > A were 0.016, 0.004–0.029, and 0.000, respectively. In conclusion, our study reports two novel variants with one each in WNT10A and EDARADD, expanding the genotypic spectra of NSTA. Second premolar agenesis is a common phenotype in affected individuals with variants in WNT10A; however, its penetrance is incomplete. Lastly, the different frequencies of WNT10A variants, c.511C > T and c.637G > A, in diverse populations might contribute to the prevalence range of NSTA between continents.

Variants of the ectodysplasin A1 receptor gene underlying homozygous cases of autosomal recessive hypohidrotic ectodermal dysplasia

Hypohidrotic ectodermal dysplasia (HED) is a rare genetic condition resulting from defective development of ectodermal derivatives, such as hair, teeth, and sweat glands. Autosomal recessive (AR) forms of HED may be caused by pathogenic variants of the ectodysplasin A1 receptor (EDAR) gene that encodes a receptor involved in the NF-κB signaling pathway. Here, we describe three cases of AR-HED in families of Turkish, Austrian, and German-American origin (with or without known consanguinity). In these cases, two out-of-frame deletions and a pathogenic missense variant of EDAR were found to be disease-causing due to reduced availability of the respective messenger RNA or impaired interaction of the encoded protein with its binding partner leading to diminished signal transduction. The same missense variant, c.1258C>T (p.Arg420Trp), has actually been reported to be restricted to the Icelandic population and to be associated with non-syndromic tooth agenesis but not HED. As our patient has no known relationship to Icelandic individuals and displays a rather severe HED phenotype, we suggest that EDAR-Arg420Trp is a more widespread variant, possibly with variable clinical expressivity.

Characterisation of a second gain of function EDAR variant, encoding EDAR380R, in East Asia

Ectodysplasin A1 receptor (EDAR) is a TNF receptor family member with roles in the development and growth of hair, teeth and glands. A derived allele of EDAR, single-nucleotide variant rs3827760, encodes EDAR:p.(Val370Ala), a receptor with more potent signalling effects than the ancestral EDAR370Val. This allele of rs3827760 is at very high frequency in modern East Asian and Native American populations as a result of ancient positive selection and has been associated with straighter, thicker hair fibres, alteration of tooth and ear shape, reduced chin protrusion and increased fingertip sweat gland density. Here we report the characterisation of another SNV in EDAR, rs146567337, encoding EDAR:p.(Ser380Arg). The derived allele of this SNV is at its highest global frequency, of up to 5%, in populations of southern China, Vietnam, the Philippines, Malaysia and Indonesia. Using haplotype analyses, we find that the rs3827760 and rs146567337 SNVs arose on distinct haplotypes and that rs146567337 does not show the same signs of positive selection as rs3827760. From functional studies in cultured cells, we find that EDAR:p.(Ser380Arg) displays increased EDAR signalling output, at a similar level to that of EDAR:p.(Val370Ala). The existence of a second SNV with partly overlapping geographic distribution, the same in vitro functional effect and similar evolutionary age as the derived allele of rs3827760, but of independent origin and not exhibiting the same signs of strong selection, suggests a northern focus of positive selection on EDAR function in East Asia.

Novel EDAR mutation in tooth agenesis and variable associated features

Tooth agenesis (TA) is the developmental absence of one or more permanent teeth. We report on 10 members of a Pakistani family afflicted with TA with variable associated features inherited in autosomal dominant fashion with full penetrance. The malformation is bilateral in the majority of cases, and hallmark feature is the absence of lateral and central incisors and canines whereas first and second premolars are involved less often. Affected individuals also have pronounced variable features associated with TA such as diastema between central incisors, overgrown labial frenum, peg-shaped lower incisors, delayed exfoliation, over-erupted upper incisors and malocclusion but have no other signs of ectodermal dysplasia. Through linkage analysis coupled with exome sequencing, we identified novel nonsense variant EDAR c.1302G>A, p.(Trp434*). The variant is deduced to create a premature termination codon that leads to the deletion of the 15 C-terminal residues. Heterozygous EDAR variants most commonly cause hypohydrotic ectodermal dysplasia, but recently one nonsense and 10 missense variants have been reported in nonsyndromic TA, some with few mild features of hypohydrotic ectodermal dysplasia. The phenotype in the family we present, the largest with EDAR-related TA reported to date, is highly variable and without any signs of ectodermal dysplasia.

Two novel mutations in MSX1 causing oligodontia

Tooth agenesis is one of the most common developmental anomalies in humans and can affect dental occlusion and speech pronunciation. Research has identified an association between mutations in MSX1, PAX9, EDA, AXIN2, WNT10A, WNT10B and LRP6 and human tooth agenesis. Two unrelated individuals with non-syndromic tooth agenesis and their families were enrolled in this study. Using Sanger sequencing of the candidate genes, we identified two novel mutations: a missense mutation c.572 T>C and a frameshift mutation c.590_594 dup TGTCC, which were both detected in the homeodomain of MSX1. After identifying the mutations, structural modeling and bioinformatics analysis were used to predict the resulting conformational changes in the MSX1 homeodomain. Combined with 3D-structural analysis of other MSX1 mutations, we propose that there is a correlation between the observed phenotypes and alterations in hydrogen bond formation, thereby potentially affecting protein binding.

EDA, EDAR, EDARADD and WNT10A allelic variants in patients with ectodermal derivative impairment in the Spanish population

BACKGROUND

Ectodermal dysplasias (ED) are a group of genetic conditions affecting the development and/or homeostasis of two or more ectodermal derivatives. An attenuated phenotype is considered a non-syndromic trait when the patient is affected by only one impaired ectodermal structure, such as in non-syndromic tooth agenesis (NSTA) disorder. Hypohidrotic ectodermal dysplasia (HED) is the most highly represented ED. X-linked hypohidrotic ectodermal dysplasia (XLHED) is the most common subtype, with an incidence of 1/50,000-100,000 males, and is associated with the EDA gene (Xq12-q13.1); the dominant and recessive subtypes involve the EDAR (2q13) and EDARADD (1q42.3) genes, respectively. The WNT10A gene (2q35) is associated more frequently with NSTA. Our goal was to determine the mutational spectrum in a cohort of 72 Spanish patients affected by one or more ectodermal derivative impairments referred to as HED (63/72) or NSTA (9 /72) to establish the prevalence of the allelic variants of the four most frequently associated genes. Sanger sequencing of the EDA, EDAR, EDARADD and WNT10A genes and multiplex ligation-dependent probe amplification (MLPA) were performed.

RESULTS

A total of 61 children and 11 adults, comprising 50 males and 22 females, were included. The average ages were 5.4 and 40.2 years for children and adults, respectively. A molecular basis was identified in 51/72 patients, including 47/63 HED patients, for whom EDA was the most frequently involved gene, and 4/9 NSTA patients, most of whom had variants of WNT10A. Among all the patients, 37/51 had variants of EDA, 8/51 had variants of the WNT10A gene, 4/51 had variants of EDAR and 5/51 had variants of EDARADD. In 42/51 of cases, the variants were inherited according to an X-linked pattern (27/42), with the remaining showing an autosomal dominant (10/42) or autosomal recessive (5/42) pattern. Among the NSTA patients, 3/9 carried pathogenic variants of WNT10A and 1/9 carried EDA variants. A total of 60 variants were detected in 51 patients, 46 of which were different, and out of these 46 variants, 12 were novel.

CONCLUSIONS

This is the only molecular study conducted to date in the Spanish population affected by ED. The EDA, EDAR, EDARADD and WNT10A genes constitute the molecular basis in 70.8% of patients with a 74.6% yield in HED and 44.4% in NSTA. Twelve novel variants were identified. The WNT10A gene has been confirmed as the second molecular candidate that has been identified and accounts for one-half of non-EDA patients and one-third of NSTA patients. Further studies using next generation sequencing (NGS) will help to identify other contributory genes in the remaining uncharacterized Spanish patients.

Homeobox code model of heterodont tooth in mammals revised

Heterodonty is one of the hallmarks of mammals. It has been suggested that, homeobox genes, differentially expressed in the ectomesenchyme of the jaw primordium along the distal-proximal axis, would determine the tooth classes (homeobox code model) based on mouse studies. Because mouse has highly specialized tooth pattern lacking canine and premolars (dental formula: 1003/1003, for upper and lower jaws, respectively), it is unclear if the suggested model could be applied for mammals with all tooth classes, including human. We thus compared the homeobox code gene expressions in various mammals, such as opossum (5134/4134), ferret (3131/3132), as well as mouse. We found that Msx1 and BarX1 expression domains in the jaw primordium of the opossum and ferret embryos show a large overlap, but such overlap is small in mouse. Detailed analyses of gene expressions and subsequent morphogenesis of tooth germ in the opossum indicated that the Msx1/BarX1 double-positive domain will correspond to the premolar region, and Alx3-negative/Msx1-positive/BarX1-negative domain will correspond to canine. This study therefore provides a significant update of the homeobox code model in the mammalian heterodonty. We also show that the modulation of FGF-mediated Msx1 activation contributes to the variation in the proximal Msx1 expression among species.

From understanding tooth development to bioengineering of teeth

Remarkable breakthroughs in the fields of developmental biology and stem cell biology during the last 15 yr have led to a new level of understanding regarding how teeth develop and how stem cells can be programmed. As a result, the possibilities of growing new teeth and of tooth bioengineering have been explored. Currently, a great deal is known about how signaling molecules and genes regulate tooth development, and modern research using transgenic mouse models has demonstrated that it is possible to induce the formation of new teeth by tinkering with the signaling networks that govern early tooth development. A breakthrough in stem cell biology in 2006 opened up the possibility that a patient's own cells can be programmed to develop into pluripotent stem cells and used for building new tissues and organs. At present, active research in numerous laboratories around the world addresses the question of how to program the stem and progenitor cells to develop into tooth-specific cell types. Taken together, the remarkable progress in developmental and stem cell biology is now feeding hopes of growing new teeth in the dental clinic in the not-too-distant future.

Germline variant in MSX1 identified in a Dutch family with clustering of Barrett's esophagus and esophageal adenocarcinoma

The vast majority of esophageal adenocarcinoma cases are sporadic and caused by somatic mutations. However, over the last decades several families have been identified with clustering of Barrett’s esophagus and esophageal adenocarcinoma. This observation suggests that one or more hereditary factors may play a role in the initiation of Barrett’s esophagus and esophageal adenocarcinoma in these families. A Dutch family with clustering of Barrett’s esophagus and esophageal adenocarcinoma was identified. Normal DNA obtained from the proband diagnosed with Barrett’s esophagus was analyzed with SNP array and exome sequencing. A custom-made panel consisting of potential germline variants was verified in the normal DNA of the affected family members. In addition, the respective tumors were analyzed for somatic loss of the wild type allele or the presence of an inactivating somatic mutation in the wild type allele. Exome sequencing revealed 244 candidate variants in the normal DNA of the proband, of which 212 variants were verified successfully. After the normal DNA of the affected family members was analyzed for the presence of the 212 potential germline variants and subsequently the respective tumors, only one potential germline variant in MSX1 (chr4: 4861985 T > G, c.359T > G, p.V120G, NM_002448) showed loss of the wild type allele in the tumor DNAs of the affected family members. A germline variant in MSX1 was identified in a Dutch family with clustering of Barrett’s esophagus and esophageal adenocarcinoma. This finding indicates that the germline defect in MSX1 may be associated with Barrett’s esophagus and cancer in this particular family.

Persistent Wnt/β-catenin signaling in mouse epithelium induces the ectopic Dspp expression in cheek mesenchyme

Tooth development is accomplished by a series of epithelial-mesenchyme interactions. Epithelial Wnt/β-catenin signaling is sufficient to initiate tooth development by activating Shh, Bmps, Fgfs and Wnts in dental epithelium, which in turn, triggered the expression of odontogenic genes in the underlying mesenchyme. Although constitutive activation of Wnt/β-catenin signaling in oral ectoderm resulted in the continuous tooth formation throughout the life span, if the epithelial Wnt/β-catenin signaling could induce the mesenchyme other than oral mesenchyme still required to be elucidated. In this study, we found that in the K14-cre; Ctnnb1^ex3f mice, the markers of dental epithelium, such as Pitx2, Shh, Bmp2, Fgf4, and Fgf8, were not only activated in the oral ectoderm, but also in the cheek epithelium. Surprisingly, the underlying cheek mesenchymal cells were elongated and expressed Dspp. Further investigations detected that the expression of Msx1 and Runx2 extended from oral to cheek mesenchyme. These findings suggested that epithelial Wnt/β-catenin signaling was capable of inducing Dspp expression in non-dental mesenchyme. Moreover, Dspp expression in the K14-cre; Ctnnb1^ex3f oral mesenchyme was activated earlier than that in the wild type littermates. In contrast, although the elongated oral epithelial cells were detected in the K14-cre; Ctnnb1^ex3f mice, the Amelogenin expression was suppressed. The differential effects of the persistent epithelial Wnt/β-catenin signaling on ameloblast and odontoblast differentiation might result from the altered BMP signaling. In summary, our findings suggested that the epithelial Wnt/β-catenin signaling could induce craniofacial mesenchyme into odontogenic program and promote odontoblast differentiation.

WNT10A variants in relation to nonsyndromic hypodontia in eastern Slovak population

Nonsyndromic hypodontia is a congenital absence of less than six permanent teeth, with a most common subtype maxillary lateral incisor agenesis (MLIA). Mutations in several genes have been described in severe tooth agenesis. The aim of this study was to search for the variants in wingless-type MMTV-integration site family member (WNT10A), paired box 9 (PAX9) and axis inhibitor 2 (AXIN2) genes, and investigate their potential role in the pathogenesis of non-syndromic hypodontia. Clinical examination and panoramic radiograph were performed in the cohort of 60 unrelated Slovak patients of Caucasian origin with nonsyndromic hypodontia including 37 MLIA cases and 48 healthy controls. Genomic DNA was isolated from buccal swabs and Sanger sequencing of WNT10A, PAX9 and AXIN2 was performed. Altogether, we identified 23 single-nucleotide variants, of which five were novel. We have found three rare nonsynonymous variants in WNT10A (p.Gly165Arg; p.Gly213Ser and p.Phe228Ile) in eight (13.33%) of 60 patients. Analysis showed potentially damaged WNT10A variant p.Phe228Ile predominantly occurred only in MLIA patients, and with a dominant form of tooth agenesis (odds ratio (OR_dom) = 9.841; P = 0.045; 95% confidence interval (CI) 0.492-196.701;OR_rec = 0.773; P = 1.000; 95% CI 0.015-39.877). In addition, the WNT10A variant p.Phe228Ile showed a trend associated with familial nonsyndromic hypodontia (P = 0.024; OR= 1.20; 95% CI 0.97-1.48). After Bonferroni correction, these effects remained with borderline tendencies. Using a 3D WNT10A protein model, we demonstrated that the variant Phe228Ile changes the proteinsecondary structure. In PAX9 and AXIN2, common variants were detected. Our findings suggest that the identified WNT10A variant p.Phe228Ile could represent risk for the inherited nonsyndromic hypodontia underlying MLIA. However, further study in different populations is required.

KDF1 is a novel candidate gene of non-syndromic tooth agenesis

OBJECTIVE

Tooth agenesis (TA) is featured by congenital loss of teeth, and can be divided into two subtypes, non-syndromic TA (NSTA) and syndromic TA (STA). Although 12 candidate genes of NSTA have been revealed, the genetic basis of NSTA needs to be further studied. We noticed an overlap of candidate genes between NSTA and STA, and hypothesized that some candidate genes of STA may be new candidate genes of NSTA.

METHODS

Sanger sequencing, whole exome sequencing, bioinformatics analyses and immunohistochemical staining were performed to reveal the genetic basis of the patients in a family with NSTA.

RESULTS

No pathogenic mutation was found in the 12 candidate genes of NSTA. We screened the variants of 76 STA candidate genes and identified a novel pathogenic mutation c.G908C (p.R303 P) in Keratinocyte Differentiation Factor 1 (KDF1). This mutation was cosegregated with the disease in the family. Bioinformatics analyses predicted the mutation to be pathogenic. Immunohistochemical staining of kdf1 in developing tooth germs indicated that kdf1 expression is important for the development of teeth.

CONCLUSIONS

This study identified KDF1 as a novel candidate gene for NSTA. STA candidate genes may be a promising source of new NSTA genes.

Next generation sequencing reveals a novel nonsense mutation in MSX1 gene related to oligodontia

Tooth agenesis is one of the most common craniofacial disorders in humans. More than 350 genes have been associated with teeth development. In this study, we enrolled 60 child patients (age 13 to 17) with various types of tooth agenesis. Whole gene sequences of PAX9, MSX1, AXIN2, EDA, EDAR and WNT10a genes were sequenced by next generation sequencing on the Illumina MiSeq platform. We found previously undescribed heterozygous nonsense mutation g.8177G>T (c.610G>T) in MSX1 gene in one child. Mutation was verified by Sanger sequencing. Sequencing analysis was performed in other family members of the affected child. All family members carrying g.8177G>T mutation suffered from oligodontia (missing more than 6 teeth excluding third molars). Mutation g.8177G>T leads to a stop codon (p.E204X) and premature termination of Msx1 protein translation. Based on previous in vitro experiments on mutation disrupting function of Msx1 homeodomain, we assume that the heterozygous g.8177G>T nonsense mutation affects the amount and function of Msx1 protein and leads to tooth agenesis.

Identification of likely pathogenic and known variants in TSPEAR, LAMB3, BCOR, and WNT10A in four Turkish families with tooth agenesis

Tooth agenesis (TA), the failure of development of one or more permanent teeth, is a common craniofacial abnormality observed in different world populations. The genetic etiology of TA is heterogeneous; more than a dozen genes have been associated with isolated or nonsyndromic TA, and more than 80 genes with syndromic forms. In this study, we applied whole exome sequencing (WES) to identify candidate genes contributing to TA in four Turkish families. Likely pathogenic variants with a low allele frequency in the general population were identified in four disease-associated genes, including two distinct variants in TSPEAR, associated with syndromic and isolated TA in one family each; a variant in LAMB3 associated with syndromic TA in one family; and a variant in BCOR plus a disease-associated WNT10A variant in one family with syndromic TA. With the notable exception of WNT10A (Tooth agenesis, selective, 4, MIM #150400), the genotype-phenotype relationships described in the present cohort represent an expansion of the clinical spectrum associated with these genes: TSPEAR (Deafness, autosomal recessive 98, MIM #614861), LAMB3 (Amelogenesis imperfecta, type IA, MIM #104530; Epidermolysis bullosa, junctional, MIMs #226700 and #226650), and BCOR (Microphthalmia, syndromic 2, MIM #300166). We provide evidence supporting the candidacy of these genes with TA, and propose TSPEAR as a novel nonsyndromic TA gene. Our data also suggest potential multilocus genomic variation, or mutational burden, in a single family, involving the BCOR and WNT10A loci, underscoring the complexity of the genotype-phenotype relationship in the common complex trait of TA.