Abstract
Objectives
This study aimed to comprehensively characterise genetic variants of amelogenesis imperfecta in a single Korean family through whole-exome sequencing and bioinformatics analysis.
Material and methods
Thirty-one individuals of a Korean family, 9 of whom were affected and 22 unaffected by amelogenesis imperfecta, were enrolled. Whole-exome sequencing was performed on 12 saliva samples, including samples from 8 affected and 4 unaffected individuals. The possible candidate genes associated with the disease were screened by segregation analysis and variant filtering. In silico mutation impact analysis was then performed on the filtered variants based on sequence conservation and protein structure.
Results
Whole-exome sequencing data revealed an X-linked dominant, heterozygous genomic missense mutation in the mitochondrial gene holocytochrome c synthase (HCCS). We also found that HCCS is potentially related to the role of mitochondria in amelogenesis. The HCCS variant was expected to be deleterious in both evolution-based and large population-based analyses. Further, the variant was predicted to have a negative effect on catalytic function of HCCS by in silico analysis of protein structure. In addition, HCCS had significant association with amelogenesis in literature mining analysis.
Conclusions
These findings suggest new evidence for the relationship between amelogenesis and mitochondria function, which could be implicated in the pathogenesis of amelogenesis imperfecta.
Clinical relevance
The discovery of HCCS mutations and a deeper understanding of the pathogenesis of amelogenesis imperfecta could lead to finding solutions for the fundamental treatment of this disease. Furthermore, it enables dental practitioners to establish predictable prosthetic treatment plans at an early stage by early detection of amelogenesis imperfecta through personalised medicine.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
Introduction
Amelogenesis imperfecta (AI) is an encompassing term used to describe the presence of numerous abnormal forms of enamel. Many cases of AI are inherited via autosomal dominant, autosomal recessive, and X-linked mechanisms [1]. The incidence of AI varies from 1:714 to 1:14,000, depending on the region [2]. Patients with AI may have dental sensitivity, loss of occlusal vertical dimension, and difficulty maintaining oral hygiene [3]. These problems lead to disability of the masticatory function and lead to poor quality of life [4]. Mastication is important not only for food intake but also for general systemic and physical functions [5]. In addition, AI can also cause aesthetic problems and low self-esteem in individuals suffering from the disease. Affected patients show reduced chances of social interaction and reduced self-consciousness because of the characteristic aspects of their facial appearance [4].
AI can be classified into four different phenotypes. Based on these phenotypes (hypoplastic, hypomaturation, hypocalcified, and hypomaturation-hypoplastic), AI was subdivided into 15 subtypes and subsequently classified by the defining characteristics of the inherited trait. Among these phenotypes, the hypoplastic type does not achieve the normal thickness of enamel, but it is distinguished from dentin by showing normal enamel density in radiograph [6]. However, different families with the disease have very distinct traits, and it is difficult to clearly perform according to subgroup classification. In addition, although the AI phenotypes are largely classified into four types, the method of dividing the detailed phenotypes is not yet clear, and these are classified in various ways [1]. This is likely because different families may have distinct characteristics, with possibly different aetiologies. Depending on the phenotype of AI, in which part of the mechanism for generating enamel damage is caused can be explained differently. Therefore, we conducted research on a single family to avoid the inclusion of other phenotypes of AI with different aetiologies. This would make it possible to clearly investigate the cause.
AI affects the tooth enamel structure and appearance. Most of the amelogenesis processes leading to the formation of the enamel are known to be affected by associated genetic variations [7, 8]. For example, mutations in the AMELX [9], ENAM [10], MMP20 [11], KLK-4 [12], FAM83H [13], WDR72 [14], C4ORF26 [15], SLC24A4 [16], LAMB3 [17], and ITGB6 [16] genes have been found and known to cause AI (non-syndromic form). AMELX and ENAM encode extracellular matrix proteins of the developing tooth enamel, and KLK-4 and MMP20 encode proteases that help degrade organic matter from the enamel matrix during the maturation stage of amelogenesis [9,10,11,12, 16]. SLC24A4 encodes a calcium transporter that mediates calcium transport to developing enamel during tooth development [16].
Amelogenesis can be divided into four stages: presecretory, secretory, transition, and maturation. In the presecretory phase, ameloblasts begin to secrete enamel matrix proteins (EMP), and in the transition phase, the secretion of these proteins decreases, and ameloblasts are reconstituted [18]. In the maturation stage, the width and thickness of the enamel crystals increase [16]. Finally, the matrix is transformed into a mature enamel with little protein [19]. For the amelogenesis process to work normally, ameloblasts must keep contact with the secreted extracellular matrix during the secretion, transition, and maturation stages. Many of the genes reported to cause AI are either related to cell-cell matrix adhesion proteins or EMP. Depending on the cause, genetically occurring mutations produce a variety of phenotypes [20].
Genetic studies of AI have been conducted since the early 2000s. Candidate gene approaches were mainly focused on protein synthesis genes involved in the amelogenesis process. Sequence changes of the candidate genes in AI patients were identified, and the effects of each sequence variation on protein expression and structure were deduced [21, 22]. Due to the increased understanding of the amelogenesis process, new candidate genes, such as SLC4A4 (a sodium bicarbonate cotransporter), have been discovered [23]. This candidate gene approach attempted to identify the main functions of these genes through animal experiments using knockout mice having null mutations and to confirm the central role of proteins, such as amelogenin [24] and enamelysin [11]. However, because of the severe heterogeneity between the groups of families, the candidate gene approach could not account for AI pathogenesis in which these candidate gene mutations were not identified. In such cases, AI-causing mutations may lie on different genes or loci that were not considered in previous studies. Therefore, additional studies are required to expand the repository of AI-causing genes [25].
With the development of genome research and DNA chip technology, genomic information can be obtained economically. The genome-wide association study (GWAS) in 2002 was the first attempt to explore genetic factors for diseases collectively. AI-related GWAS was conducted using 14 enamel hypoplasia (resembling human AI) samples from Italian greyhounds and 45 individuals in a healthy control group. This study confirmed a specific functional mutation in ENAM and revealed its recessive inheritance pattern [26].
The lower cost of next-generation sequencing has accelerated the discovery of mutations that cause AI [16]. These findings contributed to a better understanding of amelogenesis pathogenesis, including cell-cell adhesion, cell matrix adhesion, intercellular transport, signalling regulation, and the enamel mineralisation process during amelogenesis. The regulation of genes involved in these processes plays an important role in enamel formation.
In this study, we discovered a new variant in HCCS associated with AI in Koreans via whole-exome sequencing (WES). In silico mutation analysis based on sequence conservation and protein structure revealed deleterious effect of variant on HCCS function, affecting the mitochondria function. Further, literature mining analysis showed the potential association between this variant and the pathogenesis of amelogenesis. The new variant in HCCS, which is a mitochondrial protein, will help understand the relationship between the process of AI and dysfunction of mitochondria.
Materials and methods
Subjects
A family comprising 31 members across four generations was enrolled in this study. Of the 31 members, 9 were affected, and 22 were unaffected. Written consents were obtained for using the clinical medical information from 12 individuals of the family, including 8 subjects affected by AI. Since four of the participants were minors, written consent was obtained from their guardians. Saliva samples were collected to extract DNA. We analysed the family members affected by AI who gave their informed consent; these participants had no specific disease or clinical history besides amelogenesis imperfecta. This clinical research was approved by the Institutional Review Board of Yonsei University College of Dentistry (Yonsei IRB No. 2-2018-0055, Approved on 28 January 2019). This clinical study was conducted in accordance with the Helsinki Declaration.
Clinical/Radiographic assessment and measurements
The diagnosis of AI was based on clinical and radiological evidence.
Sample collection
From each consenting individual, 2 ml of saliva was collected using the Oragene DNA Self-Collection Kit (DNA Genotek Inc., Ottawa, Canada). A preservative solution in the tube was mixed with the saliva. Samples were sent to DNA Link Inc. (Seoul, South Korea) for DNA collection, extraction, and further analysis.
Whole-exome sequencing
WES was performed on a Novaseq6000 using SureSelectXT Human All Exon V5. DNA quality was confirmed using 1% agarose gel electrophoresis and PicoGreen® dsDNA Assay (Invitrogen). The SureSelect sequencing library was prepared by complying with the manufacturer’s instructions (Agilent SureSelectXT Human All Exon V5). Genomic DNA (200 ng) in 50 μl EB buffer was fragmented to 150 bp size using a Covaris-S2 instrument (Covaris). Sequencing adapters were ligated to the DNA fragment according to the manufacturer’s protocol (Agilent). The adapter-ligated DNA was amplified by PCR. A hybridisation buffer was prepared by mixing SureSelect hyb #1, #2, #3, and #4 reagents (Agilent).
The amplified DNA fragment was concentrated to 750 ng in 3.4 μl, and SureSelect blocks #1, #2, and #3 reagents (Agilent) were added. The DNA blocking agent mixture and hybridisation buffer were incubated at 95 °C for 5 min and then at 65 °C for 10 min. An RNase block (Agilent) was added to the SureSelect oligo capture library (Agilent) and incubated at 65 °C for 2 min. After the hybridisation buffer was added, a DNA blocking agent mix was added to the capture library, and the mixture was incubated at 65 °C for 24 h in a thermal cycler. Streptavidin-coated beads (50 μl), Dynabeads MyOne Streptavidin T1 (Invitrogen), were washed three times with 200 ml of SureSelect binding buffer (Agilent) and resuspended in 200 ml of binding buffer. The hybridisation mixture was added to the bead suspension and incubated for 30 min at room temperature while mixing. Beads were washed with 500 ml of SureSelect wash buffers #1 and #2 (Agilent), and DNA was eluted with 30 µl of nuclease-free water. The captured library was amplified to add index tags.
After QPCR was performed using the SYBR Green PCR Master Mix (Applied Biosystems, Thermo Fischer Scientific), the libraries tagged with equimolar amounts in the pool were combined. The Illumina Novaseq 6000 system was used according to the protocol provided for 2×100 sequencing.
WES data segregation analysis and variant filtering
Variants were identified from the WES results. Only variants in protein-coding transcripts were considered. Pedigree analysis of the family alluded to two possible inheritance patterns: autosomal dominant and X-linked dominant. To identify variants corresponding to those inheritance scenarios, we searched for heterozygous variants in the affected family members, but absent in unaffected members. After segregation analysis, we subsequently collected variants with minor allele frequency (MAF) < 5% in the Asian genome or have no reported MAF. Next, we applied additional filters to exclude the variants annotated as LOW & MODIFIER impact in SnpEff [27] (Supplement Table 1). The remaining 4 variants were subjected to in silico mutation impact analysis.
In silico mutation impact analysis
To identify deleterious mutations that may cause pathogenicity of hypoplastic AI, we used Sorting Intolerant From Tolerant (SIFT) [28] and Polymorphism Phenotyping v2 (PolyPhen2) [29]. For the 4 missense variants, all possible non-redundant protein sequences in Ensembl GRCh37.p13 were analysed. These methods help assess mutational impacts through evolutionary analysis of protein sequences. Specifically, SIFT uses sequence conservation across multiple species with the assumption that a mutation in a highly conserved site is intolerable. PolyPhen2 utilises sequence co-evolution and protein structure to predict mutational impact. The SIFT and PolyPhen2 scores for variants of the 4 candidate genes were obtained using SnpEff [27], which uses precalculated scores in DbNsfp [30]. The candidate genes in which variants were specified with mutation impact analysis are summarised in Table 1.
In silico analysis of mutation impact on protein structure
To study the effect of mutations in Val64Met on the protein structure of HCCS, we employed the 3D structure of HCCS predicted by AlphaFold2, a state-of-the-art approach for predicting protein structure [31]. We found the interacting residues of Val64 using Protein Contacts Atlas, a tool for studying protein structures at atomic resolution utilising residue-residue interaction networks [32]. The relative solvent accessibility of the residues was quantified using Naccess calculating the atomic accessible surface established by rolling a solvent molecule over the van der Waals surface [33].
Large population-based mutation impact analysis
To measure the harmful effect of abnormal function of the candidate genes, we used loss-of-function (LoF) observed/expected upper bound fraction [34] and the probability of being LoF intolerant (pLI) [35] by employing the genome aggregation database (gnomAD) v2.1.1 dataset [34], which spans 125,748 exome sequences and 15,708 whole-genome sequences from unrelated individuals sequenced as part of various disease-specific and population genetic studies. In addition, the missense z-score, which measures negative selection on missense mutation [36] of candidate genes, was also collected from the precalculated score in gnomAD v2.1.1 dataset.
Literature-based association analysis of candidate genes
To query publication articles on specific genes and diseases, we used PubTator Central [37] and queried papers with Entrez ID of each candidate gene. Next, we queried “MESH:D00567”, which is the MESH term of AI. Lastly, we queried papers on both candidate gene ID and “MESH:D00567”. To measure the significance of co-citation rate between candidate genes and AI, we used Fisher’s exact test and considered p < 0.05 as significant.
Function relationship analysis of HCCS
DAVID v6.8 [38] was used to perform GO term enrichment analysis for HCCS. To analyse the functional terms and pathways associated with the variants, we expanded the list of genes based on functional association using STRING v11.0 [39]. We extracted the first neighbour genes of HCCS in the STRING network with a 500 confidence level. Forty neighbour genes were identified (Supplement Table 2). The following categories were used in DAVID analysis: for “Gene_Ontology”: “GOTERM_BP_DIRECT”, “GOTERM_MF_DIRECT”, and “GOTERM_CC_DIRECT”; for “pathway”: “KEGG”, “REACTOME”. Table 2 presents functions significantly enriched in HCCS at p < 0.05.
Results
Patient characteristics
A total of 12 samples were used for WES analysis, consisting of samples from 7 females and 5 males, 8 of which were AI patients. All patients had permanent dentition, except for a 9-year-old patient (Fig. 1; individual 4-8) with mixed dentition.
Clinical manifestation
According to clinical and radiographic assessment, all the patients showed hypoplastic type AI in which all parts of the enamel seem not to reach normal thickness but were hard and translucent in radiographs. Clinically, the enamel was thin, and the mesiodistal width of the tooth was small such that the adjacent teeth were not in contact with each other (Fig. 2). The dentin showed normal features in radiographs, and the enamel was either very thin or not observed (Fig. 3). There was no enamel layer in severe cases, and there was severe wear on all posterior teeth in individuals 3-3 and 3-5 from the family. The overall shapes were not normal, and the surface was rough. Unusual pits and grooves were found on the occlusal surface.
Candidate genes discovered in WES analysis
WES was performed on all the consenting individuals of the family. The AI patients from the family showed two possible inheritance patterns: autosomal dominant and X-linked dominant. A total of 20 variants were identified through segregation analysis: 14 variants were accounted for by autosomal dominant inheritance, and 6 variants are explained by X-linked dominant inheritance (Table 3 and Supplement Table 3). Four of the 20 variants met all the filtering criteria (Table 1). Three missense variants in these three genes are rare in the Asian population: rs143029488 in MROH7 with Asian MAF of 0.0105, rs7718054 in FAT2 with Asian MAF of 0.0437, and rs191165757 in HCCS with Asian MAF of 0.0017. A missense variant (ENST00000240050; c.92C>G) in MTERF2 had an undefined MAF (Table 3). For example, the variant in HCCS showed a heterozygous genotype in the affected family member 2-2, but was not present in the unaffected family member 3-6 (Table 4). All 4 variants were used in the subsequent mutation impact analysis.
In silico mutation impact analysis of the variants
For the 4 genes (MROH7, MTERF2, FAT2, and HCCS) with filtered variants, we conducted mutation impact analysis using two distinct computational tools: SIFT [28] and PolyPhen2 [29]. For the 4 missense variants, canonical Ensembl sequences were analysed (Table 1). A missense variant in MROH7 was predicted to be damaging by SIFT and possibly damaging by PolyPhen2 (Table 1). A variant in FAT2 was predicted to be benign by both tools. The variant in MTERF2 was predicted to have deleterious effects by SIFT but not by polyphen2. A missense variant (c.190G>A, p.Val64Met) in HCCS was predicted to be confidently deleterious by both tools. It implied that the variant is in evolutionarily conserved position and might be under strong selective pressure. In Fig. 4a, for example, the missense variant is located in exon 3 of HCCS, and valine at amino acid position 64 in this exon is well-conserved in mammal orthologs (Fig. 4).
Next, to measure the putative impact of abnormal function of each gene, we conducted large population-based mutation impact analysis of the 4 candidate genes. From 125,748 exome sequences and 15,708 whole-genome sequences in the genome aggregation database [34], we measured organism-level impact of each candidate gene’s inactivation using LoF observed/expected upper bound fraction [34] and the pLI [35] (Table 5). Loss-of-function observed/expected upper bound fraction (LOEUF) and pLI measure the degree of selection against predicted LoF (pLoF) variation. Low LOEUF scores and high pLI scores imply strong selection against pLoF variation in a given gene, which means the malfunction of the gene has great impact on patient survival. We found that HCCS had the most deleterious effect when its function was lost among the 4 candidate genes. HCCS had the lowest LOEUF value and the highest pLI value (0.41 and 0.89, respectively). In addition, measurement of the impact of missense mutation showed that HCCS had the highest intolerance (missense z-score [36] = 1.2). For further validation of the deleterious effect of abnormal gene function, we next checked the number of reported pathogenic variants in each candidate gene using ClinVar [40], which is a public archive of reports of the relationships among human variations and phenotypes, with supporting evidence. Among the 4 candidate genes, only HCCS had reported pathogenic variants that are related to microphthalmia with linear skin defect syndrome (MLS), which is reported to occasionally have abnormal dental enamel morphology as a symptom [41]. Therefore, we conducted further evaluation of HCCS using in silico analysis of protein structure.
Potential effect of HCCS mutation on protein 3D structure
The Val64Met mutation is predicted to have a negative impact on protein structure. In silico analysis of the tertiary structure of the HCCS revealed that the Val64 residue is located in a highly packing region in the 3D structure of HCCS (Fig. 5a, red area). The Val64 has direct interactions with seven amino acids (Fig. 5a, b, yellow region). In addition, Val64 residue is located in a deeply buried region and has 8.00 relative solvent accessibility (RSA, Fig. 5a, c). A residue with a low RSA value tends to have a small exposure to solvent and is usually located in an inner space of a protein. The Val64Met substitution, a change to a larger side chain in a low RSA position, could disrupt the protein folding. Furthermore, His154, one of the direct interacting residues of Val64, is known as a critical region for the catalytic activity of HCCS [42]. Therefore, the Val64Met substitution can affect interaction with the His154 and potentially impair the protein function.
Evaluation of the effect of genes thorough literature-based association analysis
To investigate whether there is evidence on the relationship between candidate genes and AI, we used a PubMed text-mining method to conduct an independent literature-based association analysis of the 4 candidate genes. We found that only HCCS was significantly co-cited with AI in PubMed articles (p-value = \(3.23\times {10}^{-10}\), Table 6). HCCS was co-mentioned with AI in 6 articles among 906 articles queried for AI. To test whether this number is higher than expected, we measured the enrichment odds ratio and p-values between the articles mentioning HCCS and the articles mentioning AI using Fisher’s exact test.
Function enrichment analysis of HCCS and its neighbour genes
To understand candidate genes’ functional roles in the amelogenesis process, we analysed the functional terms of HCCS and its neighbour genes in the STRING network. The analysis yielded three notable categories: GO, KEGG, and REACTOME. We discovered 32 enriched functional terms for HCCS (Table 2). The topmost overrepresented term of HCCS was “mitochondrial inner membrane”.
Discussion
This study is the first to combine WES with bioinformatic tools to report hypoplastic AI in a Korean family. A missense mutation (rs191165757) in HCCS was identified through segregation analysis of WES data and variant filtering. Interestingly, we found that the MAF for the HCCS variant (rs191165757, MAF = 0.0017) was similar to the AI incidence rate (from 1:714 to 1:14,000 [2]). The similarity between the MAF for the HCCS variant and the AI incidence rate as well as the inheritance pattern of HCCS is intriguing in that this mutation may trigger AI with only a single occurrence.
Until now, most of the variants found in AI were autosomal dominant and autosomal recessive, and the only X-linked inherited mutation was that of AMELX. AMELX is an amelogenin gene, and its defect causes enamel hypoplasia. AMELX is involved in both normal thickness enamel formation and normal mineralisation processes [9]. HCCS, a variant of which was found in this study, is a newly discovered gene by X chromosome analysis. Although HCCS is close to genomic region of known AI-related genes, AMELX, no harmful variant was detected in AMELX from affected family members. Therefore, we suppose that dysfunction of HCCS itself might be related to critical defects in dental morphology. During amelogenesis, the number of ameloblasts is one of the factors that affect the amount of enamel produced. In the transition and maturation stages of amelogenesis, 25% of ameloblasts undergo apoptosis [16]. Hobson noted that HCCS protein is involved in regulating apoptosis and cell necrosis. Problems in the apoptosis process due to HCCS mutation could be a factor contributing to enamel defect in patients [43].
Additional evidence of the potential pathogenicity of abnormal function of HCCS has been provided in a large-scale population study. We found that loss of function in HCCS had a low LOEUF score (0.41, Table 5), which means it has a fatal effect on affected individuals. This corresponds to the male lethality of MLS syndrome, which is caused by loss of function in HCCS [41]. Although the discovered variant is not protein-truncating, the amino acid sequence changed, and the variant was consistently expected to have deleterious effect through in silico mutation impact analysis (Table 1). Therefore, we suppose that the missense variant in HCCS might be related to critical defects in dental morphology.
Through in silico mutation impact analysis based on sequence conservation, the missense mutation in HCCS was predicted to be deleterious due to Val64 which is highly conserved across the species. This is supported by position of Val64 residue in 3D protein structure of HCCS predicted by AlphaFold. We observed that the mutated residue (Val64) in HCCS is predicted to have direct interaction with His154 residue. In the mitochondria, HCCS catalyse cytochrome c and heme to matured holocytochrome c, and the conserved histidine (His154) in HCCS provides the key ligand to the heme iron [42]. Therefore, it is postulated that Val64Met substitution in HCCS negatively affects the synthesis of holocytochrome c, an essential mitochondrial electron carrier and an important component of the apoptosis pathway, leading to adverse effect on the mitochondria function. It suggests that one of the unknown pathogenic mechanisms of AI is underlying mitochondria dysfunction.
By discovering defective mitochondrial genes in AI patients, this study advances efforts to elucidate the role of mitochondria in amelogenesis. We found that HCCS and its first neighbour genes in a protein-protein network were significantly enriched in functions and pathways related to mitochondria: “mitochondrial protein import pathway”, “respiratory chain complex IV assembly”, and “mitochondrial electron transport” (Table 2). There are several evidences of the functional relationship between mitochondria and amelogenesis. At the maturation stage of amelogenesis, the transport of Ca2+ and PO43− is increased, and the width and thickness of the enamel crystal are expanded [44]. In this process, ameloblasts alternate between a ruffle-ended (RE) border and a smooth-ended (SE) border [19]. In some reports, the number of mitochondria in ameloblasts increased from the secretory stage to the maturation stage. In addition, mitochondria themselves also increase in size during maturation and have an efficiency that provides a lot of energy to cells [45]. Cytochrome oxidase (CO) involved in oxidative phosphorylation reflects the functional activity of mitochondria. The proportion of CO-reactive mitochondria is significantly higher in SE ameloblasts than in RE ameloblasts. This indicates that a higher energy level is required to convert from SE to RE or that SE ameloblasts may still play an unknown, highly energy-demanding role [46]. Costiniti reported an increase in oxidative phosphorylation, a measure that quantifies the oxygen consumption rate of mitochondria at the maturation stage, indicating an increase in energy demand [47]. Enamel mineralisation depends on the transport system of Ca2+ and PO43− and other ions by ion transport, channels, and pumps [19, 44, 48]. During the crystallisation of the enamel in the maturation stage, pH change adjustments occur together to buffer acidification due to proton release [49]. Besides, in maturation stage ameloblasts, mitochondria accumulate near the apical end of the RE border and show a different localisation from the secretory ameloblasts [47]. This mitochondrial accumulation can contribute to the movement of ions required for mineralisation [50].
Even though the new findings of this study are from a limited condition, the association between genetic variation and cellular function can provide worthy information for the further study of AI to understand the physiology of amelogenesis and pathogenesis of AI. In particular, this study could help identify genes that play a role in the maturation stage of amelogenesis. The discovery of HCCS revealed a genetic link to the role of mitochondria in the process of amelogenesis. In this context, this study could help uncover the mineralisation process of the maturation stage ameloblasts. Since AI has several different subtypes, it is worthy to investigate one specific family due to the idiopathic origin of the disease. For further studies with a large number of affected individuals from different families, it is important to perform proper subtype grouping after careful consideration of clinical and radiographical examinations.
References
Aldred MJ, Savarirayan R, Crawford PJ (2003) Amelogenesis imperfecta: a classification and catalogue for the 21st century. Oral Dis 9:19–23
Backman B, Holm AK (1986) Amelogenesis imperfecta: prevalence and incidence in a northern Swedish county. Community Dent Oral Epidemiol 14:43–47
Seow WK (1993) Clinical diagnosis and management strategies of amelogenesis imperfectavariants. Pediatr Dent 15:384–393
Coffield KD, Phillips C, Brady M, Roberts MW, Strauss RP, Wright JT (2005) The psychosocial impact of developmental dental defects in people with hereditary amelogenesis imperfecta. J Am Dent Assoc 136:620–630
Nakata M (1998) Masticatory function and its effects on general health. Int Dent J 48:540–548. https://doi.org/10.1111/j.1875-595x.1998.tb00489.x
Witkop CJ Jr (1988) Amelogenesis imperfecta, dentinogenesis imperfecta and dentin dysplasia revisited: problems in classification. J Oral Pathol 17:547–553
Thesleff I, Vaahtokari A, Kettunen P, Aberg T (1995) Epithelial-mesenchymal signaling during tooth development. Connect Tissue Res 32:9–15
Mina M, Kollar EJ (1987) The induction of odontogenesis in non-dental mesenchyme combined with early murine mandibular arch epithelium. Arch Oral Biol 32:123–127
Aldred MJ, Crawford PJ, Roberts E, Thomas NS (1992) Identification of a nonsense mutation in the amelogenin gene (AMELX) in a family with X-linked amelogenesis imperfecta (AIH1). Hum Genet 90:413–416. https://doi.org/10.1007/bf00220469
Rajpar MH, Harley K, Laing C, Davies RM, Dixon MJ (2001) Mutation of the gene encoding the enamel-specific protein, enamelin, causes autosomal-dominant amelogenesis imperfecta. Hum Mol Genet 10:1673–1677. https://doi.org/10.1093/hmg/10.16.1673
Caterina JJ, Skobe Z, Shi J, Ding Y, Simmer JP, Birkedal-Hansen H, Bartlett JD (2002) Enamelysin (matrix metalloproteinase 20)-deficient mice display an amelogenesis imperfecta phenotype. J Biol Chem 277:49598–49604. https://doi.org/10.1074/jbc.M209100200
Wang SK, Hu Y, Simmer JP et al (2013) Novel KLK4 and MMP20 mutations discovered by whole-exome sequencing. J Dent Res 92:266–271. https://doi.org/10.1177/0022034513475626
Kim JW, Lee SK, Lee ZH, Park JC, Lee KE, Lee MH, Park JT, Seo BM, Hu JC, Simmer JP (2008) FAM83H mutations in families with autosomal-dominant hypocalcified amelogenesis imperfecta. Am J Hum Genet 82:489–494. https://doi.org/10.1016/j.ajhg.2007.09.020
El-Sayed W, Parry DA, Shore RC et al (2009) Mutations in the beta propeller WDR72 cause autosomal-recessive hypomaturation amelogenesis imperfecta. Am J Hum Genet 85:699–705. https://doi.org/10.1016/j.ajhg.2009.09.014
Parry DA, Brookes SJ, Logan CV et al (2012) Mutations in C4orf26, encoding a peptide with in vitro hydroxyapatite crystal nucleation and growth activity, cause amelogenesis imperfecta. Am J Hum Genet 91:565–571. https://doi.org/10.1016/j.ajhg.2012.07.020
Smith PJA, Antanaviciute A, Kirkham J, Brookes SJ, Inglehearn CF, Mighell AJ (2017) Amelogenesis imperfecta; genes, proteins, and pathways. Front Physiol 8:435. https://doi.org/10.3389/fphys.2017.00435
Poulter JA, El-Sayed W, Shore RC, Kirkham J, Inglehearn CF, Mighell AJ (2014) Whole-exome sequencing, without prior linkage, identifies a mutation in LAMB3 as a cause of dominant hypoplastic amelogenesis imperfecta. Eur J Hum Genet 22:132–135. https://doi.org/10.1038/ejhg.2013.76
Reith EJ (1970) The stages of amelogenesis as observed in molar teeth of young rats. J Ultrastruct Res 30:111–151. https://doi.org/10.1016/s0022-5320(70)90068-7
Smith, (1998) Cellular and chemical events during enamel maturation. Crit Rev Oral Biol Med 9:128–161. https://doi.org/10.1177/10454411980090020101
Gadhia K, McDonald S, Arkutu N, Malik K (2012) Amelogenesis imperfecta: an introduction. Br Dent J 212:377–379. https://doi.org/10.1038/sj.bdj.2012.314
Kim JW, Simmer JP, Lin BP, Seymen F, Bartlett JD, Hu JC (2006) Mutational analysis of candidate genes in 24 amelogenesis imperfecta families. Eur J Oral Sci 114 Suppl 1:3–12; discussion 39–41, 379. https://doi.org/10.1111/j.1600-0722.2006.00278.x
Mardh CK, Backman B, Holmgren G, Hu JC, Simmer JP, Forsman-Semb K (2002) A nonsense mutation in the enamelin gene causes local hypoplastic autosomal dominant amelogenesis imperfecta (AIH2). Hum Mol Genet 11:1069–1074. https://doi.org/10.1093/hmg/11.9.1069
Urzua B, Ortega-Pinto A, Morales-Bozo I, Rojas-Alcayaga G, Cifuentes V (2011) Defining a new candidate gene for amelogenesis imperfecta: from molecular genetics to biochemistry. Biochem Genet 49:104–121. https://doi.org/10.1007/s10528-010-9392-6
Gibson CW, Yuan ZA, Hall B et al (2001) Amelogenin-deficient mice display an amelogenesis imperfecta phenotype. J Biol Chem 276:31871–31875. https://doi.org/10.1074/jbc.M104624200
Wright JT, Torain M, Long K, Seow K, Crawford P, Aldred MJ, Hart PS, Hart TC (2011) Amelogenesis imperfecta: genotype-phenotype studies in 71 families. Cells Tissues Organs 194:279–283. https://doi.org/10.1159/000324339
Gandolfi B, Liu H, Griffioen L, Pedersen NC (2013) Simple recessive mutation in ENAM is associated with amelogenesis imperfecta in Italian Greyhounds. Anim Genet 44:569–578. https://doi.org/10.1111/age.12043
Cingolani P, Platts A, le Wang L, Coon M, Nguyen T, Wang L, Land SJ, Lu X, Ruden DM (2012) A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly (Austin) 6:80–92. https://doi.org/10.4161/fly.19695
Vaser R, Adusumalli S, Leng SN, Sikic M, Ng PC (2016) SIFT missense predictions for genomes. Nat Protoc 11:1–9. https://doi.org/10.1038/nprot.2015.123
Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, Kondrashov AS, Sunyaev SR (2010) A method and server for predicting damaging missense mutations. Nat Methods 7:248–249. https://doi.org/10.1038/nmeth0410-248
Liu X, Li C, Mou C, Dong Y, Tu Y (2020) dbNSFP v4: a comprehensive database of transcript-specific functional predictions and annotations for human nonsynonymous and splice-site SNVs. Genome Med 12:103. https://doi.org/10.1186/s13073-020-00803-9
Jumper J, Evans R, Pritzel A et al (2021) Highly accurate protein structure prediction with AlphaFold. Nature 596:583–589. https://doi.org/10.1038/s41586-021-03819-2
Kayikci M, Venkatakrishnan AJ, Scott-Brown J, Ravarani CNJ, Flock T, Babu MM (2018) Visualization and analysis of non-covalent contacts using the Protein Contacts Atlas. Nat Struct Mol Biol 25:185–194. https://doi.org/10.1038/s41594-017-0019-z
Lee B, Richards FM (1971) The interpretation of protein structures: estimation of static accessibility. J Mol Biol 55:379–400. https://doi.org/10.1016/0022-2836(71)90324-x
Karczewski KJ, Francioli LC, Tiao G et al (2020) The mutational constraint spectrum quantified from variation in 141,456 humans. Nature 581:434–443. https://doi.org/10.1038/s41586-020-2308-7
Lek M, Karczewski KJ, Minikel EV et al (2016) Analysis of protein-coding genetic variation in 60,706 humans. Nature 536:285–291. https://doi.org/10.1038/nature19057
Samocha KE, Robinson EB, Sanders SJ et al (2014) A framework for the interpretation of de novo mutation in human disease. Nat Genet 46:944–950. https://doi.org/10.1038/ng.3050
Wei CH, Allot A, Leaman R, Lu Z (2019) PubTator central: automated concept annotation for biomedical full text articles. Nucleic Acids Res 47:W587-w593. https://doi.org/10.1093/nar/gkz389
Dennis G Jr, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA (2003) DAVID: database for annotation, visualization, and integrated discovery. Genome Biol 4:P3
Szklarczyk D, Gable AL, Lyon D et al (2019) STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res 47:D607-d613. https://doi.org/10.1093/nar/gky1131
Landrum MJ, Chitipiralla S, Brown GR et al (2020) ClinVar: improvements to accessing data. Nucleic Acids Res 48:D835-d844. https://doi.org/10.1093/nar/gkz972
Prakash SK, Cormier TA, McCall AE, Garcia JJ, Sierra R, Haupt B, Zoghbi HY, Van Den Veyver IB (2002) Loss of holocytochrome c-type synthetase causes the male lethality of X-linked dominant microphthalmia with linear skin defects (MLS) syndrome. Hum Mol Genet 11:3237–3248. https://doi.org/10.1093/hmg/11.25.3237
San Francisco B, Bretsnyder EC, Kranz RG (2013) Human mitochondrial holocytochrome c synthase’s heme binding, maturation determinants, and complex formation with cytochrome c. Proc Natl Acad Sci USA 110:E788-797. https://doi.org/10.1073/pnas.1213897109
Hobson GM, Gibson CW, Aragon M, Yuan ZA, Davis-Williams A, Banser L, Kirkham J, Brook AH (2009) A large X-chromosomal deletion is associated with microphthalmia with linear skin defects (MLS) and amelogenesis imperfecta (XAI). Am J Med Genet A 149a:1698–1705. https://doi.org/10.1002/ajmg.a.32968
Lacruz RS (2017) Enamel: molecular identity of its transepithelial ion transport system. Cell Calcium 65:1–7. https://doi.org/10.1016/j.ceca.2017.03.006
Assis GF, Ribeiro DA, Campos PD, Cestari TM, Taga R (2003) Comparative stereologic study between secretory and maturation ameloblasts in rat incisors. J Appl Oral Sci 11:144–149. https://doi.org/10.1590/s1678-77572003000200011
Ohshima H, Maeda T, Takano Y (1998) Cytochrome oxidase activity in the enamel organ during amelogenesis in rat incisors. Anat Rec 252:519–531. https://doi.org/10.1002/(sici)1097-0185(199812)252:4%3c519::Aid-ar3%3e3.0.Co;2-i
Costiniti V, Bomfim GH, Li Y, Mitaishvili E, Ye ZW, Zhang J, Townsend DM, Giacomello M, Lacruz RS (2020) Mitochondrial function in enamel development. Front Physiol 11:538. https://doi.org/10.3389/fphys.2020.00538
Lacruz RS, Smith CE, Bringas P Jr, Chen YB, Smith SM, Snead ML, Kurtz I, Hacia JG, Hubbard MJ, Paine ML (2012) Identification of novel candidate genes involved in mineralization of dental enamel by genome-wide transcript profiling. J Cell Physiol 227:2264–2275. https://doi.org/10.1002/jcp.22965
Lacruz RS, Nanci A, Kurtz I, Wright JT, Paine ML (2010) Regulation of pH during amelogenesis. Calcif Tissue Int 86:91–103. https://doi.org/10.1007/s00223-009-9326-7
Eckstein M, Aulestia FJ, Nurbaeva MK, Lacruz RS (2018) Altered Ca(2+) signaling in enamelopathies. Biochim Biophys Acta Mol Cell Res 1865:1778–1785. https://doi.org/10.1016/j.bbamcr.2018.04.013
Funding
This work was supported in part by the National Research Foundation of Korea grants (2021R1F1A104718511, 2020R1A6A1A03047902, and 2021R1A2B5B01001903).
Author information
Authors and Affiliations
Contributions
H.C. and K.L. contributed equally by conceiving and designing the study. H.C. conducted clinical observation, collection of clinical data, and interpretation of clinical symptoms. K.L. performed statistical segregation analysis and gene functional analysis. D.K. performed mutation impact analysis. H.C., K.L., S.K., and J.H.L. interpreted analysis results. H.C., K.L., D.K., S.K., and J.H.L. drafted the manuscript.
Corresponding authors
Ethics declarations
Ethics approval
This clinical research was approved by the Institutional Review Board of Yonsei University College of Dentistry (Yonsei IRB No. 2-2018-0055, approved on 28 January 2019). This clinical study was conducted in accordance with the Helsinki Declaration.
Consent to participate
Informed consent was obtained from all individual participants or legal guardians included in the study.
Conflict of interest
Hyejin Choi declares that she has no conflict of interest. Kwanghwan Lee declares that he has no conflict of interest. Donghyo Kim declares that he has no conflict of interest. Sanguk Kim declares that he has no conflict of interest. Jae Hoon Lee declares that he has no conflict of interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Choi, H., Lee, K., Kim, D. et al. The implication of holocytochrome c synthase mutation in Korean familial hypoplastic amelogenesis imperfecta. Clin Oral Invest 26, 4487–4498 (2022). https://doi.org/10.1007/s00784-022-04413-0
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00784-022-04413-0