An inclusive study of deleterious missense PAX9 variants using user-friendly tools reveals structural, functional alterations, as well as potential therapeutic targets.

Mutations in the PAX9 are responsible for non-syndromic tooth agenesis in humans, although their structural and functional consequences on protein phenotype, stability, and posttranslational modifications (PTMs) have not yet been adequately investigated. This in silico study focuses on retrieving the six most deleterious mutations (L21P, R26W, R28P, G51S, I87F, and K91E) of PAX9 that has been linked to severe oligodontia. Several computational algorithm methods were used to determine the deleterious effects of PAX9 mutations. Analysis of gene ontology, protein interactions, and PTMs indicated significant functional changes caused by PAX9 mutations. The structural superimposition of the wild-type and mutant PAX9 variants revealed structural changes in locations that were present in the structures of all six variations. The conserved domain analysis revealed that the areas shared by all six variations contained unique sections that lacked DNA binding or protein-protein interaction sites, suggesting prospective drug target sites for functional restoration. The protein-protein interaction network showed KDM5B as PAX9's strongest interacting partner similar to MSX1. The PAX9 protein's structural conformations, compactness, stiffness, and function may all be impacted by changes, according to MD simulations. In addition, research on cell lines and animal models may be valuable in establishing their specific roles in functional annotations.

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

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

Four Novel PAX9 Variants and the PAX9-Related Non-Syndromic Tooth Agenesis Patterns.

The purpose of this research was to investigate and identify PAX9 gene variants in four Chinese families with non-syndromic tooth agenesis. We identified pathogenic gene variants by whole-exome sequencing (WES) and Sanger sequencing and then studied the effects of these variants on function by bioinformatics analysis and in vitro experiments. Four novel PAX9 heterozygous variants were identified: two missense variants (c.191G > T (p.G64V) and c.350T > G (p.V117G)) and two frameshift variants (c.352delC (p.S119Pfs*2) and c.648_649insC(p.Y217Lfs*100)). The bioinformatics analysis showed that these variants might be pathogenic. The tertiary structure analysis showed that these four variants could cause structural damage to PAX9 proteins. In vitro functional studies demonstrated that (1) the p.Y217Lfs*100 variant greatly affects mRNA stability, thereby affecting endogenous expression; (2) the p. S119Pfs* 2 variant impairs the subcellular localization of the nuclear expression of the wild-type PAX9 protein; and (3) the four variants (p.G64V, p.V117G, p.S119Pfs*2, and p.Y217Lfs*100) all significantly affect the downstream transcriptional activity of the BMP4 gene. In addition, we summarized and analyzed tooth missing positions caused by PAX9 variants and found that the maxillary second molar (84.11%) and mandibular second molar (84.11%) were the most affected tooth positions by summarizing and analyzing the PAX9-related non-syndromic tooth agenesis positions. Our results broaden the variant spectrum of the PAX9 gene related to non-syndromic tooth agenesis and provide useful information for future genetic counseling.

Characterization of PAX9 variant P20L identified in a Japanese family with tooth agenesis.

Transcription factors PAX9 and MSX1 play crucial roles in the development of permanent teeth at the bud stage, and their loss-of-function variants have been associated with congenital tooth agenesis. We sequenced the coding regions of the PAX9 and MSX1 genes from nine patients with non-syndromic tooth agenesis, and identified a missense mutation, P20L, of PAX9 in a single familial case involving three patients in two generations. Identical mutation was previously reported by other authors, but has not been characterized in detail. The mutation was located in a highly conserved N-terminal subdomain of the paired domain and co-segregated as a heterozygote with tooth agenesis. The patients showed defects primarily in the first and second molars, which is typical for cases attributable to PAX9 mutation. Luciferase reporter assay using the 2.3-kb promoter region of BMP4 and electrophoretic mobility shift assay using the CD19-2(A-ins) sequence revealed that P20L substitution eliminated most of the transactivation activity and specific DNA binding activity of PAX9 under the experimental conditions we employed, while some residual activity of the mutant was evident in the former assay. The hypomorphic nature of the variant may explain the relatively mild phenotype in this case, as compared with other PAX9 pathogenic variants such as R26W.

G-quadruplex formation enhances splicing efficiency of PAX9 intron 1.

G-quadruplexes are secondary structures present in DNA and RNA molecules, which are formed by stacking of G-quartets (i.e., interaction of four guanines (G-tracts) bounded by Hoogsteen hydrogen bonding). Human PAX9 intron 1 has a putative G-quadruplex-forming region located near exon 1, which is present in all known sequenced placental mammals. Using circular dichroism (CD) analysis and CD melting, we showed that these sequences are able to form highly stable quadruplex structures. Due to the proximity of the quadruplex structure to exon-intron boundary, we used a validated double-reporter splicing assay and qPCR to analyze its role on splicing efficiency. The human quadruplex was shown to have a key role on splicing efficiency of PAX9 intron 1, as a mutation that abolished quadruplex formation decreased dramatically the splicing efficiency of human PAX9 intron 1. The less stable, rat quadruplex had a less efficient splicing when compared to human sequences. Additionally, the treatment with 360A, a strong ligand that stabilizes quadruplex structures, further increased splicing efficiency of human PAX9 intron 1. Altogether, these results provide evidences that G-quadruplex structures are involved in splicing efficiency of PAX9 intron 1.

Mutational analysis of AXIN2, MSX1, and PAX9 in two Mexican oligodontia families.

The genes for axin inhibition protein 2 (AXIN2), msh homeobox 1 (MSX1), and paired box gene 9 (PAX9) are involved in tooth root formation and tooth development. Mutations of the AXIN2, MSX1, and PAX9 genes are associated with non-syndromic oligodontia. In this study, we investigated phenotype and AXIN2, MSX1, and PAX9 gene variations in two Mexican families with non-syndromic oligodontia. Individuals from two families underwent clinical examinations, including an intra-oral examination and panoramic radiograph. Retrospective data were reviewed, and peripheral blood samples were collected. The exons and exon-intronic boundaries of the AXIN2, MSX1, and PAX9 genes were sequenced and analyzed. Protein and messenger RNA structures were predicted using bioinformative software programs. Clinical and oral examinations revealed isolated non-syndromic oligodontia in the two Mexican families. The average number of missing teeth was 12. The sequence analysis of exons and exon-intronic regions of AXIN2, MSX1, and PAX9 revealed 11 single-nucleotide polymorphisms (SNPs), including seven in AXIN2, two in MSX1, and three in PAX9. One novel SNP of MSX1, c.476T>G (Leu159Arg), was found in all of the studied patients in the families. MSX1 Leu159Arg and PAX9 Ala240Pro change protein and messenger RNA structures. Our findings suggested that a combined reduction of MSX1 and PAX9 gene dosages increased the risk for oligodontia in the Mexican families, as in vivo investigation has indicated that interaction between Msx1 and Pax9 is required for tooth development.

Pathogenic mechanisms of tooth agenesis linked to paired domain mutations in human PAX9.

Mutations in the paired-domain transcription factor PAX9 are associated with non-syndromic tooth agenesis that preferentially affects posterior dentition. Of the 18 mutations identified to date, eight are phenotypically well-characterized missense mutations within the DNA-binding paired domain. We determined the structural and functional consequences of these paired domain missense mutations and correlated our findings with the associated dental phenotype variations. In vitro testing included subcellular localization, protein-protein interactions between MSX1 and mutant PAX9 proteins, binding of PAX9 mutants to a DNA consensus site and transcriptional activation from the Pax9 effector promoters Bmp4 and Msx1 with and without MSX1 as co-activator. All mutant PAX9 proteins were localized in the nucleus of transfected cells and physically interacted with MSX1 protein. Three of the mutants retained the ability to bind the consensus paired domain recognition sequence; the others were unable or only partly able to interact with this DNA fragment and also showed a similarly impaired capability for activation of transcription from the Msx1 and Bmp4 promoters. For seven of the eight mutants, the degree of loss of DNA-binding and promoter activation correlated quite well with the severity of the tooth agenesis pattern seen in vivo. One of the mutants however showed neither reduction in DNA-binding nor decrease in transactivation; instead, a loss of responsiveness to synergism with MSX1 in target promoter activation and a dominant negative effect when expressed together with wild-type PAX9 could be observed. Our structure-based studies, which modeled DNA binding and subdomain stability, were able to predict functional consequences quite reliably.

Identification and functional analysis of two novel PAX9 mutations.

The paired-domain transcription factor PAX9 plays a critical role in tooth development, as heterozygous mutations in PAX9 have been shown to be associated with human tooth agenesis. In this study, we report 2 novel missense mutations, gly6arg (G6R) and ser43lys (S43K), in the paired domain of PAX9 in Chinese patients with varying degrees of nonsyndromic tooth agenesis. Excluding third molars, the individual with the G6R mutation was missing 2 mandibular incisors and a maxillary premolar, while the phenotype of individuals with the S43K mutation consisted of peg-shaped upper lateral incisors and missing molars, premolars and canines. As these 2 mutations occur at highly conserved amino acids in the PAX gene family and between different species, we further analyzed the effects of the mutations on the function of the resulting proteins. Immunofluorescence and immunoblotting studies showed that the mutations did not alter nuclear localization in mammalian cells. Gel shift and super shift assays indicate that both mutant proteins bound DNA at a lower level than the normal protein, with G6R having a greater affinity for DNA than S43K. Likewise, the G6R protein was able to transcriptionally activate a Bmp4 promoter construct to a greater extent than S43K. Our finding that the severity of tooth agenesis in the patients was correlated to the DNA-binding capacity of the mutated PAX9 9proteins supports the hypothesis that DNA binding is responsible for the genetic defect.

Natural selection and molecular evolution in primate PAX9 gene, a major determinant of tooth development.

Large differences in relation to dental size, number, and morphology among and within modern human populations and between modern humans and other primate species have been observed. Molecular studies have demonstrated that tooth development is under strict genetic control, but, the genetic basis of primate tooth variation remains unknown. The PAX9 gene, which codes for a paired domain-containing transcription factor that plays an essential role in the development of mammal dentition, has been associated with selective tooth agenesis in humans and mice, which mainly involves the posterior teeth. To determine whether this gene is polymorphic in humans, we sequenced approximately 2.1 kb of the entire four-exon region (exons 1, 2, 3 and 4; 1,026 bp) and exon-intron (1.1 kb) boundaries of 86 individuals sampled from Asian, European, and Native American populations. We provided evidence that human PAX9 polymorphisms are limited to exon 3 only and furnished details about the distribution of a mutation there in 350 Polish subjects. To investigate the pattern of selective pressure on exon 3, we sequenced ortholog regions of this exon in four species of New World monkeys and one gorilla. In addition, orthologous sequences of PAX9 available in public databases were also analyzed. Although several differences were identified between humans and other species, our findings support the view that strong purifying selection is acting on PAX9. New World and Old World primate lineages may, however, have different degrees of restriction for changes in this DNA region.

PAX9 and TGFB3 are linked to susceptibility to nonsyndromic cleft lip with or without cleft palate in the Japanese: population-based and family-based candidate gene analyses.

The prevalence of nonsyndromic cleft lip with or without cleft palate (CL/P) and cleft palate only (CPO) are believed to be higher in the Japanese than in Americans, Europeans or Africans. The purpose of this study was to investigate, in a Japanese population, relationships between CL/P or CPO and seven candidate genes (TGFB3, DLX3, PAX9, CLPTM1, TBX10, PVRL1, TBX22) that showed positive associations in other populations and are expressed in the oral/lip region in developing mice. We first searched for mutations in these genes among 112 CL/P and 16 CPO patients, and found a heterozygous missense mutation (640A > G, S214G) in exon 3 of PAX9 in two sibs with CL/P and their phenotypically normal mother from a Japanese family. A population-based case-control analysis and a family-based transmission disequilibrium test (TDT), using single nucleotide polymorphisms (SNPs), and two-SNP haplotypes of the genes, between the 112 CL/P cases with their parents and 192 controls indicated a significant association at one SNP site, IVS1 + 5321, in TGFB3 with a P-value of 0.0016. Population-based haplotyping revealed that the association was most significant for haplotype "A/A" consisting of IVS1 + 5321 and IVS1 - 1572; TDT also gave a P-value of 0.0252 in this haplotype.