Journal of the College of Physicians and Surgeons Pakistan
ISSN: 1022-386X (PRINT)
ISSN: 1681-7168 (ONLINE)
Affiliations
doi: 10.29271/jcpsp.2022.05.570ABSTRACT
Objective: To discover novel ectodysplasin‐A (EDA) and wingless‐type MMTV integration site family, member 10A (WNT10A) mutations in tooth agenesis (TA) patients.
Study Design: Case series.
Place and Duration of Study: Guanghua School of Stomatology, Guangzhou, China, from March 2018 to August 2020.
Methodology: EDA and WNT10A were analysed in eleven TA families by PCR and Sanger sequencing. Bioinformatics and structure modelling analyses were performed after identifying different variants, to predict the resulting conformational alterations in WNT10A and EDA.
Results: Two novel mutations (c.796C>A (p.L266I), c.769G>A (p.G257R)) in EDA and two reported mutations (c.637G>A (p.G213S), c.511C>T (p.R171C))in WNT 10A were detected. Combined with the 3D structural analysis, we discovered a correlation between alterations in hydrogen bond formation and the observed phenotypes, potentially affecting protein binding.
Conclusions: The mutations were predicted to be pathogenic through bioinformatics analyses. In addition, by identifying novel mutations, our knowledge regarding the TA spectrum and tooth development was considerably expanded.
Key Words: Anodontia, EDA, WNT 10A, Whole exome sequencing, Odontogenesis.
INTRODUCTION
Congenital tooth agenesis (TA) is a common developmental abnormality. The congenital loss of teeth might cause masticatory, speech, esthetic, and psychological problems, which put a heavy burden on the patients and associated societies. (From preliminary studies, the reported prevalence of this disease ranges from 2.2% to 10.1%.1,2 Whether the patient has other systemic diseases, the agenesis is classified into nonsyndromic tooth agenesis(NSTA) and syndromic tooth agenesis(STA).3
TA can be inherited through an X-linked, autosomal recessive mode or an autosomal dominant trait. Morphogenesis and tooth growth are contingent on the interrelationship between morphogenesis and the epithelium in the course of embryonic development.
During the development of the tooth, over 300 signalling molecules, transcription factors and growth factors are involved in various signalling pathways including the FGF pathway, the SHH pathway, WNT/β-catenin pathway, along with the NF-κB pathway.4,5
Associations have been identified between TA and mutations in ectodysplasin‐A (EDA) and wingless-type MMTV integration site family, member 10A (WNT10A).4,6 Nevertheless, despite records of genomic variants in affected populations, many unidentified genomic mutations causing TA have not been identified to date. Furthermore, the signaling pathways and molecular mechanisms underlying these disorders have not been fully clarified. Therefore, the aim of this study was to discover novel EDA and WNT10A genomic mutations and reveal the genetic basis of TA patients. This expands the TA spectrum and provide a new theoretical basis for the clinical prevention and treatment of TA.
METHODOLOGY
Pakistan with X-linked hypohidrotic ectodermal dysplasia (XLHED). All cases were referred to the Hospital of Stomatology, Sun Yat-Sen University from March 2018 to August 2020 and eleven unrelated Chinese families carrying TA were identified. The absence of permanent teeth was not due to extraction or injuries and was confirmed by all probands. A detailed clinical examination of the systemic and dentition conditions confirmed the diagnosis of NSTA in nine patients and STA in two patients. When needed, dentists, dermatologists, and other doctors thoroughly evaluated the systemic conditions, including the nails, hair, sweat and teeth of the patients.
Criteria with two of the three following systems was considered criteria for XLHED: hypohidrosis, sparse hair, and oligodontia and included.7 Criteria was excluded once the diagnosis of another type of ectodermal dysplasia was found.8
This study was approved by the Ethical Review Committee at the Guanghua School and Hospital of Stomatology at Sun Yat-Sen University. TA patients were recruited, and basic information on TA and its genetic pathogenicity was provided. Related examination and experimental methods were also performed. Their treatment was guaranteed, and convenience was provided if they failed to meet our requirements after the examination. Written informed consent was obtained from the participants, and the recommendations of the Declaration of Helsinki were followed.
From each individual, 4 ml of peripheral blood was collected. Genomic DNA was extracted using the QiaAmp Kit (Qiagen, Düsseldorf, Germany). Applying Oligo 7.0, primers were designed covering the flanking intronic sequences of WNT10A and EDA along with the exons. In previous work, EDA primer sequences were used.9 The products were sequenced using an ABI 3730XL genetic analyser. Sequence Scanner Software v1.0 was utilised to analyse the sequencing results. The nomenclature of mutation, with +1 corresponding to the A of the ATG translation initiation codon of the reference sequences NM_001399.4 (EDA) and NM_025216.2 (WNT10A), was applied.9
Earlier studies described the technique for the bioinformatics analyses performed here.6,9 The Human Protein Reference Database (http://hprd.org/query) provided the domain information for WNT10A and EDA. Applying CLUSTAL X (1.83), the human EDA amino acid sequence (ENST00000374552) was briefly compared to those of the rhesus (ENSMMUT00000024953), cattle (ENSBTAT00000016649), mouse (ENSMUST000001137 79), and chicken (ENSGALT00000007137). PROVEAN, Ploy Phen2, Mutation Taster, and SIFT enabled the prediction of the pathogenic influence of the undiscovered mutations.10-12
Swiss PBD Viewer was used to model the structure of the wild-type and mutant WNT10A and EDA proteins. The structures of WNT-8 (PDB ID 4F0A; X-ray, resolution 3.25 Å) and EDA (PDB ID 1RJ7; X-ray, resolution 2.3 Å) were employed for homology modelling. Figures of three-dimensional structures were built using PyMol v2.4 (The PyMOL Molecular Graphics System, Version 2.4 Schrödinger, LLC., Cambridge, MA, USA).
RESULTS
Congenital absence of teeth excluding other symptoms was observed in three NSTA families, carrying normally shaped teeth. Afflictions other than hypotrichosis, hypohidrosis and hypodontia, were found in two STA patients. Table I shows the numbers of missing permanent teeth. The panoramic radiographs confirmed the diagnosis of TA (Figure 1A–E). Since permanent teeth could have replaced some primary teeth, the numbers of primary teeth of all individuals were not evaluated.
By analysing EDA in XLHED Family 1, we identified one undiscovered mutation c.796C>A (p.Leu266Ile) in the patient and the mother. His mother was an asymptomatic heterozygous carrier of the mutation, while the father was wild-type at c.796. Moreover, a mutation c.637G>A (p.Gly213Ser) was found during our screening of WNT10A. His father carried the mutation with no symptoms and his mother was wild-type at c.637 (Figure 1).
Figure 1: Pedigrees and corresponding mutation of XLHED Families 1 and 2 and NSTA Families 1–3. (A) In XLHED Family 1, an EDA mutation c.796C>A (p.Leu266Ile) was discovered, and a known WNT10A mutation c.637G>A(p.Gly213Ser) was found. (B) One novel EDA mutation c.769G>A (p.Gly257Arg) was discovered in XLHED Family 1. (C) In NSTA Family 1, the patient had the heterozygous WNT10A mutation c.511C>T (p.Arg171Cys). (D)Moreover, for NSTA Family 2, one covered homogeneous WNT10A mutation c.637G>A (p.Gly213Ser) was confirmed within the gene of the patient. (E) For NSTA Family 3, the patient was confirmed to carry the WNT10A mutation c.637G>A (p.Gly213Ser)). Black arrows indicate the probands. Each individual for whom blood samples were obtained is indicated with ‘+’. Arrows indicate the mutations. XLHED: X‐linked hypohidrotic ectodermal dysplasia. NSTA: nonsyndromic tooth agenesis.
Table I: Clinical data and corresponding mutations of EDA and WNT10A genes.
Family Number |
Proband |
Gender & Age |
Corresponding gene |
Nucleotide Varients |
Genotype |
Predicted Protein Changes |
Number of absent Teeth # |
XLHED family 1 |
II:1 |
Male,4Y |
EDA |
c.796G>A |
Hemizygous |
p.Leu266Ile |
22 |
WNT10A |
c.637C>T |
Heterozygous |
p.Gly213Ser |
||||
XLHED family 2 |
II:1 |
Male,8Y |
EDA |
c.769C>A |
Hemizygous |
p.Gly257Arg |
21 |
NSTA family 1 |
II:1 |
Male,12Y |
WNT10A |
c.511C>T |
Heterozygous |
p.Arg171Cys |
2 |
NSTA family2 |
II:2 |
Female,7Y |
WNT10A |
c.637G>A |
Homozygous |
p.Gly213Ser |
11 |
NSTA Family3 |
II:1 |
Female,7Y |
WNT10A |
c.637G>A |
Heterozygous |
p.Gly213Ser |
7 |
#The third molars and deciduous teeth were not included. Bold type: undiscovered mutation; Genotype: Homozygous, Heterozygous or Hemizygous. XLHED: X‐linked hypohidrotic ectodermal dysplasia. NSTA: nonsyndromic tooth agenesis. |
Figure 2: Tumour necrosis factor homology domain of wild-type EDA protein (2-A), Leu266Ile EDA protein (2-B) and Gly257Arg EDA protein (2-C). The three trimers were similar in morphology and structure but slightly different in primary structures. There were five more amino acid residues in the Leu266Ile EDA protein (B) and Gly257Arg EDA protein (2-C) than in the wild-type EDA protein (2-A): GLY242, THR243, ARG244, GLU245 and ASN246. These five amino acids had close hydrogen bond interactions with the surrounding residues. In addition, in the Leu266Ile EDA protein (2-B), the original LEU266 was replaced by ILE.
By sequencing the EDA genes in the XLHED Family 2, we uncovered a new mutation c.769G>A (p.Gly257Arg), shared by the proband and his mother. His father was wild-type at c.769.
Both EDA and WNT10A genes were examined in the patients of three NSTA families. A WNT10A mutation c.511C>T (p.Arg171Cys) was discovered within NSTA Family 1. His mother carried the mutation with no symptoms. Blood specimens of the father were not available. For Family 2, the covered homozygous WNT10A mutation c.637G>A (p.Gly213Ser) was found in the genome of the proband. Both the father and the mother were heterozygous at c.637, manifesting no symptoms. Within NSTA Family 3, we discovered the covered WNT10A mutation c.637G>A (p.Gly213Ser) in the proband as well as her mother (Figure 1).
As mentioned before, EDA variants c.796C>A were discovered in XLHED Family 1 and c.769G>A in XLHED Family 2. Bioinformatic analyses and experiments were performed to determine whether they were pathogenic.
The EDA mutation c.796C>A (p.Leu266Ile) found in XLHED Family 1 was not published in Pubmed, the Human Gene Mutation Database (HGMD; public version) or 1000 Genomes. A leucine was substituted for an isoleucine because of the variant at the position. It was by PolyPhen2 predicted to be ‘possibly damaging’. Moreover, the prediction of Mutation Taster was ‘disease-causing’ since the mutation transferred the amino acid sequence and changed the splice site.
Another mutation, c.769G>A (p.Gly257Arg), was also unavailable in Pubmed, the Human Gene Mutation Database or 1000 Genomes. A glycine was replaced with an arginine caused by the variant at the position. PolyPhen2 predicted it as being ‘probably damaging’. A status of ‘disease-causing’ was predicted by Mutation Taster because it altered the amino acid sequence. Moreover, the mutation may obstruct the link between the 2nd and 3rd exon areas, resulting in protein features changes.
Structural modelling of the tumour necrosis factor homology domain of the wild-type EDA protein, Leu266Ile EDA protein and Gly257Arg EDA protein (Figure 2) indicated similar morphologies and structures but slightly different primary structures. There were five more amino acid residues in the Leu266Ile EDA protein (Figure 2-B) and Gly257Arg EDA protein (Figure 2-C) than in the wild-type EDA protein (Figure 2-A): GLY242, THR243, ARG244, GLU245 and ASN246. These five amino acids had close hydrogen bond interactions with the surrounding residues. In addition, in the Leu266Ile EDA protein, the original LEU266 was replaced by ILE.
DISCUSSION
This study presents synthetic genetic research on NSTA and XLHED. Genetic mutations in five families carrying NSTA or XLHED were discovered along with Sanger sequencing. EDA mutations c.796C>A (p.Leu266Ile) as well as c.769G>A (p.Gly257Arg) were novel. Structural modelling and bioinformatics analyses were performed to analyse the pathogenicity of the mutations.
Notably, in the XLHED Family 1, simultaneous EDA and WNT10A mutations resulted in STA. In the XLHED Family 1, the novel mutation c.796C>A (p.Leu266Ile) was present in the patient along with the mother. In addition, a previously identified mutation c.637G>A (p.Gly213Ser) was found in the patient and his father through our screening of WNT10A. Both parents were asymptomatic. In previous research, the WNT10A mutation c.637G>A (p.Gly213Ser) was identified and analysed, and it was considered to contribute to TA.13 Therefore, whether the newly discovered mutation can lead to XLHED is worthy of detailed discussion. Conversely, the corresponding mutation is the most significant pathogenic factor in XLHED, since EDA results in more than half of TA cases.14 Moreover, according to the bioinformatics analyses, including cross-species comparison, PolyPhen2, and structural modelling, c.796C>A (p.Leu266Ile) is likely to be pathogenic. The phenomenon in which the mother had no symptoms could be explained by female carriers of EDA mutations presenting with a normal or very mild phenotype, since the EDA gene is located on chromosome X.15,16 In contrast, the phenotype of the patient seemed more severe than the previous single WNT10A mutation c.637G>A (p.Gly213Ser) patient reported previously.13 Clinical manifestations further support this conjecture. Previous studies have indicated that mutations in different genes could also lead to one disease in one patient.17,18
EDA and WNT10A mutations can influence each other resulting in TA.19 In this patient, both EDA and WNT10A mutations were found, and as mentioned earlier in the article, analyses indicated that the EDA and WNT10A mutations may affect protein function. Therefore, it was conjectured that the interaction between EDA and WNT10A mutations might result in XLHED.
The same WNT10A mutation c.637G>A(p.Gly213Ser) was found in XLHED Family 1, NSTA family2 and NSTA family 3. In NSTA family 2, the proband was homozygous and showed TA. Both the father and the mother were heterozygous at c.637, showing no symptoms.
In XLHED Family 1, the proband and his father carried the same mutation gene, and were both heterozygous. While the proband showed TA, his father was asymptomatic. This phenomenon also occurred in the proband and her mother in NSTA Family 3. This genetic phenomenon might be explained as irregular dominance and was reported in TA before.19,20 Other gene factors may be involved in regulating WNT10A c.637 gene expression during tooth development. Individuals with different genetic backgrounds could be an important cause of irregular dominant inheritance. Environmental factors might also be one of the reasons for the irregular dominance. Nevertheless, further research is needed in the future to reveal the mechanism in a larger sample.
The proband and his mother shared another new EDA mutation c.769G>A (p.Gly257Arg) in XLHED Family 2. His father was wild-type at c.769. For EDA mutation c.769G>A, a glycine was replaced with an arginine, caused by the mutation at the position. In the previous report, the same site but different nucleotide mutation of EDA c.769G>C(p.G257R) was found in the patient with NSTA.21 It is speculated that tooth development at this site is more sensitive to the mutation of EDA, providing a more sufficient basis for c.769G>A to be a pathogenic mutation. This can be explained as the mutation of EDA at c.769, affecting the protein and playing an important role in tooth development. PolyPhen2 and Mutation Taster predicted it to be disease-causing as it altered the amino acid sequence. This conclusion is also supported by the comparative structural modelling of the EDA protein.
In recent years, there have been many reports about novel mutations involved EDA-EDAR-NF-κB and WNT/beta-catenin signaling pathways in TA.22,23 A close relationship is likely to exist between EDA and WNT pathways. Previous studies showed that the WNT and EDA pathways had a sequential interdependency during hair follicle morphogenesis and development.24 Moreover, a cascade of WNT and EDA signaling was reported to be essential for touch dome Merkel cell development.25 Therefore, there might be an interaction between EDA and WNT10A mutations, which together result in abnormalities in odontogenesis.
CONCLUSIONS
This research confirmed two novel EDA mutations c.769G>A (p.Gly257Arg) and c.796C>A (p.Leu266Ile) and two discovered WNT10A mutations c.511C>T (p.Arg171Cys) and c.637G>A (p.Gly213Ser) among eleven TA patients. This study extrapolated the two undiscovered mutations as pathogenic representing a contribution to precision medicine.
ETHICAL APPROVAL:
This study was approved by the Ethical Review Committee at the Guanghua School and Hospital of Stomatology at Sun Yat-Sen University.
PATIENT’S CONSENT:
Informed consent was obtained from patients to publish the data concerning this case.
COMPETING INTERESTS:
The authors declared no conflict of interest.
ACKNOWLEDGEMENTS:
This project was funded by the National Natural Sciences Foundation of China (81873711 and 81974146) and the Science and the Fundamental Research Funds of the Central Universities, Sun Yat-sen University (19ykpy86), and China Postdoctoral Science Foundation (2020M673023).
AUTHORS’ CONTRIBUTION:
WX: Conceptualisation, investigation, methodology, writing-original draft, and writing-review and editing.
BZ: Data curation, investigation, validation, and writing-original draft.
Pei Li: Data curation, and methodology.
Duoling Xu: Data curation.
DY: Data curation, investigation, methodology, writing-original draft, and writing-review and editing.
Wei Zhao: Conceptualisation, funding acquisition, supervision, writing-original draft, writing-review and editing.
All authors approved the final version of the manuscript to be published.
REFERENCES