Genes and Pathways Associated with Skeletal Sagittal Malocclusions: A Systematic Review

Skeletal class II and III malocclusions are craniofacial disorders that negatively impact people’s quality of life worldwide. Unfortunately, the growth patterns of skeletal malocclusions and their clinical correction prognoses are difficult to predict largely due to lack of knowledge of their precise etiology. Inspired by the strong inheritance pattern of a specific type of skeletal malocclusion, previous genome-wide association studies (GWAS) were reanalyzed, resulting in the identification of 19 skeletal class II malocclusion-associated and 53 skeletal class III malocclusion-associated genes. Functional enrichment of these genes created a signal pathway atlas in which most of the genes were associated with bone and cartilage growth and development, as expected, while some were characterized by functions related to skeletal muscle maturation and construction. Interestingly, several genes and enriched pathways are involved in both skeletal class II and III malocclusions, indicating the key regulatory effects of these genes and pathways in craniofacial development. There is no doubt that further investigation is necessary to validate these recognized genes’ and pathways’ specific function(s) related to maxillary and mandibular development. In summary, this systematic review provides initial insight on developing novel gene-based treatment strategies for skeletal malocclusions and paves the path for precision medicine where dental care providers can make an accurate prediction of the craniofacial growth of an individual patient based on his/her genetic profile.


Introduction
Among the three dimensions of craniofacial structure, the sagittal dimension is the primary focus and most crucial plane in orthodontic diagnosis and treatment planning [1]. Skeletal class II and III malocclusions are disorders with a sagittal discrepancy between the maxilla and the mandible [2,3]: patients diagnosed with skeletal class II malocclusion exhibit an anterior position of the maxilla in comparison with the mandible, which could result from mandibular retrognathism (reduced lower jaw growth and size) and/or maxillary prognathism (excessive upper jaw growth and size) [2,4]. In contrast, skeletal class III malocclusion, characterized by a posterior maxilla position related to the mandible, could Int. J. Mol. Sci. 2021, 22, 13037 2 of 17 be attributed to mandibular prognathism, maxillary retrognathism, or a combination of the two [3,5]. Both skeletal class II and III malocclusions significantly diminish patients' ability to chew food effectively [6]. In addition, patients with skeletal class III malocclusion have a relatively high level of abrupted, discontinuous masticatory motions and an impaired ability to swallow a bolus [7]. Skeletal class II and III malocclusions may also deleteriously affect the functions of other components of the digestive system [8]. Moreover, skeletal class II and III malocclusions have been linked to a diversity of maladies that concern oral health, such as myogenic temporomandibular disorders [9], and even extend to other parts of the body, such as myofascial pain [10] and gastroesophageal reflux disease [8]. Furthermore, from a cosmetic standpoint, skeletal class II and III malocclusions are known to be visually unappealing traits [11,12], which may also harm patients' psychological well-being and decrease their perception of their own worth. In contrast, surgical correction could markedly increase their self-esteem and decrease anxieties about their visual appearance [11]. Together, skeletal class II and III malocclusions significantly negatively impact people's quality of life [13]. Strikingly, approximately a quarter of the general population suffers from one or the other of these types of malocclusions [14]. Indeed, the global incidence of skeletal class II malocclusion is estimated at nearly 20% and that of skeletal class III malocclusion at 6% [14]. Therefore, treating these types of malocclusions and thereby rectifying patients' dentition would tremendously benefit individuals and society.
To date, orthognathic surgical correction is still the last-resort treatment of severe skeletal class II and III malocclusions to achieve balanced anatomic and functional jaw relationships and improve facial appearance [15]. However, as an invasive procedure, orthognathic surgery is associated with various, sometimes severe, risks, including impaired sensation, infection, bleeding, tissue injury, and alterations in the osteosynthesis systems [16]. More importantly, to achieve proper diagnosis and avoid additional surgeries to correct malocclusion relapse due to post-surgery growth, most orthognathic surgeries can only be achieved when patients are skeletally matured [17]. Furthermore, the longstanding expectation of surgery during childhood and adolescence will worsen patients' physical, social, and psychological health. In order to avoid surgical operations and the consequences associated with such treatments, a large amount of effort has been dedicated to developing growth modification appliances for skeletal class II and III malocclusion correction, especially those that can be delivered around the puberty stage to achieve the normal dentition without surgery [18][19][20][21]. Unfortunately, despite the tremendous efforts that have been devoted to developing a diversity of treatment technologies [11,[18][19][20][21][22][23][24][25], dental care providers are not yet able to predict whether an orthodontic functional appliance will successfully correct a patient's skeletal discrepancy or if that patient will not respond to early orthopedic correction and require surgery to properly rectify their malocclusion later. In other words, people's expected development of skeletal malocclusion and their clinical correction prognoses are currently difficult to predict, which can be largely attributed to the present lack of knowledge of the precise etiology of skeletal class II and III malocclusions.
Since a strong inheritance pattern of a specific type of skeletal malocclusion running in families has been observed in various races [14,26], there is no doubt that genomic alteration(s) play essential roles in skeletal malocclusion development. Here, previous genome-wide association studies (GWAS) are summarized and reanalyzed to identify skeletal class II and III malocclusion-associated genes and signal pathways. Therefore, by discussing the genetic factors that may regulate the abnormal maxillary and mandibular growth that cause skeletal class II and III malocclusions, this systematic review aims to provide initial insight into developing novel gene-based treatment strategies for skeletal malocclusions and paves the path for precision medicine where dental care providers can make an accurate prediction of the craniofacial growth of an individual patient based on his/her genetic profile.

Methods
This systematic review followed the 2020 Preferred Reporting Items for Systemic Reviews and Meta-Analyses (PRISMA) guidelines [27], and the article search was conducted in June 2021. The PubMed and Google Scholar databases were searched using the following keywords: "gene", "genetic", "GWAS", "mandible", "mandibular", "prognathism", "retrognathism", "skeletal class II", and "skeletal class III". No publication date limits were set. Reviews were excluded from consideration to avoid double-counting of genes, and only studies conducted on humans and written in English were included. Studies were evaluated independently by two authors (E. G. and C. L.) and were selected for use in the current systematic review if both authors determined that the studies were eligible ( Figure 1). It is important to note that the genes found in the databases pertained to skeletal class II and III malocclusions, with some research specifically focused on a given gene's association with the maxilla or mandible only. Because gene expression in the mandible and the maxilla can contribute to skeletal class II and III malocclusions, both potential sources of malocclusion were considered in this review.
Skeletal classes II or III malocclusion-associated genes were then separately entered into the Reactome database for pathway enrichment [28]. The pathways associated with each skeletal malocclusion with a p-value and a false discovery rate (FDR) less than 0.05 were recorded, and the top 20 pathways with the lowest FDR values were further discussed. Pathway functional annotations were based on PubMed focusing on bone, cartilage, osteoblasts, osteoclasts, and osteocytes. Functions of genes not recognized in the Reactome database were annotated by the Uniprot database [29]. make an accurate prediction of the craniofacial growth of an individual patient based on his/her genetic profile.

Methods
This systematic review followed the 2020 Preferred Reporting Items for Systemic Reviews and Meta-Analyses (PRISMA) guidelines [27], and the article search was conducted in June 2021. The PubMed and Google Scholar databases were searched using the following keywords: "gene", "genetic", "GWAS", "mandible", "mandibular", "prognathism", "retrognathism", "skeletal class II", and "skeletal class III". No publication date limits were set. Reviews were excluded from consideration to avoid double-counting of genes, and only studies conducted on humans and written in English were included. Studies were evaluated independently by two authors (E. G. and C. L.) and were selected for use in the current systematic review if both authors determined that the studies were eligible ( Figure 1). It is important to note that the genes found in the databases pertained to skeletal class II and III malocclusions, with some research specifically focused on a given gene's association with the maxilla or mandible only. Because gene expression in the mandible and the maxilla can contribute to skeletal class II and III malocclusions, both potential sources of malocclusion were considered in this review.
Skeletal classes II or III malocclusion-associated genes were then separately entered into the Reactome database for pathway enrichment [28]. The pathways associated with each skeletal malocclusion with a p-value and a false discovery rate (FDR) less than 0.05 were recorded, and the top 20 pathways with the lowest FDR values were further discussed. Pathway functional annotations were based on PubMed focusing on bone, cartilage, osteoblasts, osteoclasts, and osteocytes. Functions of genes not recognized in the Reactome database were annotated by the Uniprot database [29].

Gene Identification
Using the PubMed and Google Scholar databases, 15 human studies targeting skeletal class II malocclusion were collected. Nineteen skeletal class II malocclusion-associated

Gene Identification
Using the PubMed and Google Scholar databases, 15 human studies targeting skeletal class II malocclusion were collected. Nineteen skeletal class II malocclusion-associated genes were identified, among which five genes (actinin alpha 3 (ACTN3), ADAM metallopeptidase with thrombospondin type 1 motif 9 (ADAMTS9), fibroblast growth factor receptor

Pathway Enrichment
A total of 17 of these 19 skeletal class II malocclusion-associated genes were noted by the Reactome database [28] and yielded 129 pathways (Supplemental Table S1). Excitingly, 78 of the 129 enriched pathways had a p-value and an FDR less than 0.05 (Supplemental Table S1), representing strong confidence in the enrichment results from the skeletal class II malocclusion data. Based on the FDR, the top 20 enriched pathways ( Table 2) were selected for further discussion regarding their functions related to bone and cartilage growth and development.

The Skeletal Class II Malocclusion-Associated Genes Not Recognized by the Reactome Database
Two genes, MSX1 and MYO1H, that are highly associated with skeletal class II malocclusion (appearing twice in the database search) were not indexed in the Reactome database. Consequently, the Uniprot database [29] was utilized to gain insights into the functions of MSX1 and MYO1H instead (Table 3). Table 3. Functional annotation of the skeletal class II malocclusion-associated genes not recognized by the Reactome database.

MSX1
activation of meiosis, anterior/posterior pattern specification, BMP signaling pathway involved in heart development, bone morphogenesis, cardiac conduction system development, cartilage morphogenesis, cell morphogenesis, cellular response to nicotine, embryonic forelimb morphogenesis, embryonic hindlimb morphogenesis, embryonic morphogenesis, embryonic nail plate morphogenesis, epithelial to mesenchymal transition involved in endocardial cushion formation, face morphogenesis, in utero embryonic development, mammary gland epithelium development, mesenchymal cell proliferation, midbrain development, middle ear morphogenesis, muscle organ development, negative regulation of apoptotic process, negative regulation of cell growth, negative regulation of cell population proliferation, negative regulation of striated muscle cell differentiation, negative regulation of transcription regulatory region DNA binding, odontogenesis of dentin-containing tooth, pituitary gland development, positive regulation of BMP signaling pathway, positive regulation of DNA damage response, signal transduction by p53 class mediator, positive regulation of intrinsic apoptotic signaling pathway by p53 class mediator, positive regulation of mesenchymal cell apoptotic process, protein localization to nucleus, protein stabilization, regulation of odontogenesis, regulation of transcription by RNA polymerase II, central roof of mouth development, signal transduction involved in regulation of gene expression, stem cell differentiation

MYO1H
actin filament organization, vesicle transport along actin filament

Pathway Enrichment
Forty-four skeletal class III malocclusion-associated genes were recognized in the Reactome database, resulting in 425 enriched pathways (Supplemental Table S2). The enrichment results yielded 129 pathways with a p-value and an FDR less than 0.05 (Supplemental Table S2). Among them, the top 20 enriched pathways were selected for further discussion (Table 5). Nine identified skeletal class III-associated genes, including myosin heavy chain 1 (MYH1), chromosome 1 open reading frame 167 (C1orf167), homeobox C cluster (HOXC), NBPF member 8 (NBPF8), NBPF member 9 (NBPF9), transcription factor 21 (TCF21), MYO1H, SSX2IP, and calneuron 1 (CALN1), were not found in the Reactome database. Notably, although 8 of these 9 genes appeared only once in previous GWAS studies, MYO1H was examined by two independent studies ( Table 4). As above, the Uniprot database [29] was utilized to functionally annotate these skeletal class III malocclusion-associated genes that had not been indexed in the Reactome database (Table 6). However, functional annotation of C1orf167, NBPF8, and NBPF9 failed in the Uniprot database, suggesting that knowledge of these genes is largely lacking. Table 6. Functional annotation of the skeletal class III malocclusion-associated genes not recognized by the Reactome database.

MYH1 muscle contraction
C1orf167 - cell adhesion, centrosome cycle, cilium assembly, intraciliary transport involved in cilium assembly, regulation of cell motility, regulation of Rac protein signal transduction

CALN1
Negatively regulates Golgi-to-plasma membrane trafficking by interacting with PI4KB and inhibiting its activity. May play a role in the physiology of neurons and is potentially important in memory and learning.

Single Nucleotide Polymorphism (SNP) Analysis of The Genes Involved in Both Skeletal Class II and III Malocclusions
It is worth noting that four genes (COL1A1, FGFR2, MATN1, and MYO1H), were associated with both class II and class III skeletal malocclusions, indicating that these genes may play key regulatory roles in craniofacial growth, especially in the sagittal dimension. One possibility is that for each of these four genes, different SNPs may lead to a more active isoform or higher expression in one class malocclusion and a more inert isoform or less expression in the other. However, available GWAS studies [30,34,35,38,40,43,[46][47][48]53,54] revealed that certain SNPs, including rs2249492 (COL1A1), rs11200014 (FGFR2), rs2162540 (FGFR2), rs10850110 (MYO1H), and rs3825393 (MYO1H), were shared by skeletal class II and class III malocclusions (Table 7), which cannot be explained by the hypothesis mentioned above. Surprisingly, most of the reported skeletal class II or class III malocclusion-associated SNPs are located in the introns of these four genes [74-81], indicating the importance of the noncoding regions of these DNA fragments in craniofacial development and skeletal malocclusion.  [85]. Moreover, FGFR2 mutations have previously been linked to bone development and growth diseases, such as Apert syndrome [86], a genetic syndrome characterized by untimely early fusion of the skull bones during development that displays skeletal class III malocclusion [87,88]. Noticeably, Apert syndrome can arise from the S252W mutation in FGFR2 [88]. In this study, enrichment results showed that FGFR2 is involved in 19 of the top 20 skeletal class III malocclusion-associated pathways (Table 5), echoing the assertion that this receptor and the associated FGFR pathway play a vital role in the sagittal disharmony of the maxillomandibular complex. Meanwhile, FGFR2 was also found to be linked to skeletal class II malocclusion (Table 1) and appeared in 17 of the top 20 skeletal class II malocclusion-associated pathways ( Table 2), suggesting that FGFR2's influence is not limited to skeletal class III malocclusion but extends to skeletal class II malocclusion as well.
Regarding FGFR2's downstream signal transduction, "Phospholipase C-mediated cascade; FGFR2" is one of the enriched pathways shared by skeletal class II and III malocclusions (Tables 2 and 5). In this pathway, phospholipase C-γ (PLCγ) is a substrate of FGFR and other receptors with tyrosine kinase activity [95]. Particularly, a previous study revealed that FGFR2 is responsible for PLCγ2 signaling activation in rat osteoblasts [96]. Furthermore, it is worth noting that PLCγ2 also promotes bone resorption by upregulating the expression of nuclear factor of activated T-cells, cytoplasmic 1 (NFATc1), a transcription factor that plays a central role in promoting osteoclastogenesis [97]. Consequently, impairment of the PLCγ2 pathway depresses osteoclastogenesis [98]. In addition, NFATc1 is a skeletal class III malocclusion-associated gene (Table 4). Thus, the FGFR2→PLCγ2→NFATc1 signal axis may be one of the essential bone growth and development regulating pathways that modulate osteoblast and osteoclast activities and thus contribute to skeletal malocclusion establishment. Meanwhile, since more FGFR2 SNPs were reported in skeletal class III malocclusion than in skeletal class II malocclusion, the influence of the FGFR2→PLCγ2→NFATc1 cassette may be more common and notable for skeletal class III malocclusion progression; however, further investigation targeting a population of greater diversity should be conducted to test this hypothesis. Downregulated by FGFR2 [99], the "PI3K Cascade" is another skeletal class II and III malocclusion-associated pathway (Tables 2 and 5). Previous studies showed that inhibition of PI3K/p70 S6K cascades increases osteoblastic differentiation induced by bone morphogenetic protein 2 (BMP2) [100] while decreasing osteoclast activity [101]. These results suggest that the "PI3K Cascade" may govern bone growth and development concerning both osteoblasts and osteoclasts, while further research is needed to parse out how this pathway is associated with skeletal class II and class III malocclusions.
Interestingly, all reported FGFR2 SNPs are located in introns ( Table 7). As opposed to exons, which directly encode protein sequences, introns integrally regulate gene expression [102]. Some FGFR2 SNPs, such as rs10736303, rs1078806, and rs2981578, were only detected in skeletal class III malocclusion (Table 7) [35], suggesting that these SNPs may be responsible for mandibular prognathism and/or maxillary retrognathism. Surprisingly, two FGFR2 SNPs from the same intron, rs11200014 and rs2162540, are linked to both skeletal class II and III malocclusions (Table 7) [34,35]. It is difficult to explain the influence of these two FGFR2 SNPs if FGFR2 is considered as functioning in an isolated manner. One possible explanation is that after post-transcriptional modification, the respective introns are removed from the pre-mRNA and serve as miRNAs that regulate the expression of other skeletal malocclusion-associated genes, while the SNPs significantly alter their regulation targets and/or their effectiveness. Alternatively, these two intronic SNPs may make FGFR2 expression more sensitive to other influences and thus increase the risk of abnormal skeletal malocclusions. These additional influences may not be limited to other growth factors or cytokines involved in bone and cartilage growth and development. It is also possible that these two intronic FGFR2 SNPs indicate a scenario in which the affected bone is more vulnerable and sensitive to mechanical stimulations, such as those from the attached muscle tissue, as proposed by functional matrix theory [103,104].

Insulin Receptor Cascade
"Insulin receptor signaling cascade", which was identified by pathway enrichment of the genes associated with skeletal class II and class III malocclusions (Tables 2 and 5), is known to promote osteoblast differentiation and osteocalcin (OCN) secretion [105][106][107]. Moreover, deletion of the Insulin-like growth factor 1 receptor (IGF1R) gene in mouse osteoblasts led to a striking decrease in cancellous bone volume, connectivity, and trabecular number, accompanied by an increase in trabecular spacing, while the rate of mineralization of osteoid was significantly decreased [108]. This correlation suggests that the "Insulin receptor signaling cascade" is essential for coupling matrix biosynthesis to sustained min-eralization, which is particularly important during the pubertal growth spurt when rapid bone formation and consolidation are required [108].
Both insulin receptor substrate 1 (IRS1) and insulin receptor substrate 2 (IRS2) are involved in the signal transduction from IGF1R to its downstream mediator RAC-alpha serine/threonine-protein kinase (Akt1) [109]. Previous studies demonstrate that IRS1 in osteoblasts is indispensable for maintaining bone turnover [110]. In contrast, IRS2 is more predominantly required for bone formation over bone resorption [111]. Indeed, Akune et al. found that IRS2 signaling is not essential for osteoclastic cells' differentiation, functioning, or survival [111]. Therefore, balanced IRS1/IRS2 signal transduction may contribute to the foundation of normal bone development and growth, while breaking the equilibrium of the "Insulin receptor signaling cascade" may trigger the formation of skeletal abnormalities, including skeletal class II and III malocclusions.
Because they are both connected to the PI3K→Akt axis, a complex interaction exists between FGF-related signal transduction and the "Insulin receptor signaling cascade" [99,109,112]. For example, not only is IRS1 involved in the "Insulin receptor signaling cascade" [109], but it is also a part of the "PI3K cascade" [113], and thus its influence may also extend to osteoclasts. In addition, stimulating the "Insulin receptor signaling cascade" could activate the parathyroid hormone (PTH) type 1 receptor (PTH1R) in osteoblasts/osteocytes to enhance the osteoblast-to-osteocyte transition [114]. Thus, the "Insulin receptor signaling cascade" may act as an intermediate coordinator of various signal pathways to integrate functions of theirs that are related to bone growth and development regulation and may warrant further investigation regarding its association with craniofacial sagittal discrepancy.

Runt-Related Transcription Factor 2 (RUNX2) and Notch Receptor 3 (NOTCH3)
As a well-known, predominant regulator of bone and cartilage growth and development [115,116], runt-related transcription factor 2 (RUNX2) promotes the activity of FGFR2 and FGFR3, leading to the proliferation of pre-osteoblastic cells [117]. Considering that RUNX2 promotes the differentiation of secondary chondrocytes and the formation of the cartilage of mandibular condyles [118], genetic variations that are more likely to have a negative influence on the RUNX2 gene may more frequently lead to mandibular retrognathism and thus induce skeletal class II malocclusion. This hypothesis is supported by previous animals studies that have shown that RUNX2-deficiency in mice resulted in a lack of mandibular condylar cartilage and mandibular bone [119]. In the current study, the RUNX2 gene was identified as the seed for three skeletal class II malocclusion-associated enriched pathways, as expected (Table 2).
On the other hand, the NOTCH signaling pathway was associated with skeletal class III malocclusion in the current study (Table 5), which is in agreement with previous research that demonstrates that a reduction in NOTCH3 signaling promoted osteogenic differentiation of mesenchymal stem stems located in a certain part of the mandible [120].
Together, the class-specific pathway identification of RUNX2 and NOTCH3 establishes confidence in the current skeletal malocclusion-associated enrichment.

Functional Matrix Theory
Excitingly, some skeletal class II and III malocclusion-associated genes identified in the current study have muscle-related functional annotations (Tables 1 and 4). One possible explanation for this is functional matrix theory, also called functional matrix hypothesis [103,104].
Functional matrix theory proposes that the growth and development of bone, including qualities such as attained length and width, are largely influenced by other body components [103,104]. For example, changes in the skeleton are not fully promoted by the genetic makeup of the bones but also rely on other biological components, particularly the skeletal muscles that belong to the periosteal category of the functional matrix [103].
Skeletal muscle activity could generate endogenous electrical fields that might orchestrate bone growth and development through various mechanisms [103]. On the other hand, the epigenetic event of muscle contraction may extend to bone cells' genomes and thus modulate bone growth and development [121]. Therefore, as demonstrated in the current study (Tables 1 and 4), although genes that regulate skeletal muscle generation and activities may not be directly involved in bone and cartilage growth and development, a number of them nevertheless play a role in the growth and development of the maxilla and mandible, influencing characteristics such as length, and consequently are involved in the development of skeletal class II and III malocclusions. On the other hand, SNPs found in genes associated with bone and cartilage growth and development, either in introns or exons, may cause the gene isoform to be more sensitive to intracellular mechanotransduction and/or more vulnerable to muscle contraction-initiated epigenetic modification, and thus contribute to skeletal malocclusion determination.

Myosin Heavy Chain (MYH) Genes
Several MYH genes were previously associated with both types of malocclusions (Tables 1 and 4). In particular, higher expression levels of MYH3, MYH6, and MYH7 were detected in patients with mandibular retrognathism than in those with mandibular prognathism, while no significant difference was found between these two populations in regard to MYH1, MYH2, and MYH8 expression levels [44]. Thus, it is safe to say that the MYH proteins could be future research targets for craniofacial skeletal sagittal growth prediction and modification. At the same time, further investigations are warranted to determine whether, which of, and how the MYH molecules are involved in skeletal malocclusion establishment and if and what characters were overlooked in these previous investigations.

Other Muscle Function Genes
Histone deacetylase 4 (HDAC4) was found to be more highly expressed in the masseter muscle of patients with skeletal class III malocclusion than in patients with skeletal class II malocclusion [51]. In addition, two independent studies revealed higher expression of lysine acetyltransferase 6B (KAT6B) in the masseter muscle of skeletal class III malocclusion patients than in the same tissue type belonging to skeletal class II malocclusion patients, thereby indicating that KAT6B may be associated with mandibular prognathism [51,52]. Clearly, these genes could also be candidates for predicting patients' expected development of skeletal malocclusion and their clinical correction prognoses in the future.
It is important to note that GWAS can only indicate the association between identified genes and a targeted clinical condition, while functional validation at the molecular level is a prerequisite for forming any conclusions. As stated above, current knowledge of the functions of the skeletal class II and III-associated genes in the musculoskeletal system is extensive. However, the specific impacts of these genes on the craniofacial region are largely unknown. We hope that the current review and the conducted pathway enrichment can provide clues about the potential relationships among the identified genes and identify which genes and pathways can be prioritized for functional validation in future studies. We would also like to emphasize that the current review focused on the sagittal dimension, which is a crucial dimension in orthodontic diagnosis and treatment planning. Because craniofacial structure develops three-dimensionally, further evaluation of the genes involved in the vertical and transverse growth and development of the craniofacial structure is also needed.

Conclusions
By reviewing and reanalyzing available human studies, 19 genes were found to be associated with skeletal class II malocclusion and 53 genes with skeletal class III malocclusion. Using the Reactome database, most of these genes were enriched in pathways related to bone and cartilage regulation and growth, as expected. Interestingly, multiple muscle-related genes and pathways have also been identified, which aligns with functional matrix theory in that the muscles surrounding the jaws may also contribute to the development of these two types of malocclusions. In addition, several SNPs were associated with both skeletal class II and III malocclusions, with the majority of these SNPs being located in the introns of their respective genes. It is possible that these SNPs result in a more sensitive isoform of the respective genes in response to factors associated with the surrounding microenvironment, such as mechanotransductive stimulation and/or pre-and post-translational regulations that might include epigenetic modification and miRNA-level interaction. Although multiple confirmation studies must be conducted, we hope our current study can provide a relatively panoramic view of the current understanding of skeletal class II and III malocclusions at the genetic level and thus potentially help dental care providers make more accurate craniofacial sagittal growth predictions and treatment strategy selections based on the patients' genetic profiles in the future.

Conflicts of Interest:
The authors declare no conflict of interest.