Whole genome sequencing of a family with autosomal dominant features within the oculoauriculovertebral spectrum.

Background: Oculoauriculovertebral Spectrum (OAVS) encompasses a wide variety of anomalies on derivatives from the first and second pharyngeal arches including macrostomia, hemifacial microsomia, micrognathia, preauricular tags, ocular and vertebral anomalies. We present the genetic findings of a large three-generation family with multiple members affected with macrostomia, preauricular tags and uni- or bilateral ptosis following an autosomal dominant segregation pattern. Methods: We generated whole genome sequencing data for the proband, affected parent and unaffected paternal grandparent followed by Sanger sequencing on 23 family members for the top 10 candidate genes: KCND2, PDGFRA, CASP9, NCOA3, WNT10A, SIX1, MTF1, KDR/VEGFR2, LRRK1, and TRIM2. We performed parent and sibling-based transmission disequilibrium tests and burden analysis to explore segregation and burden of candidate gene mutations. Bioinformatic analyses investigated the biological connection between genes and the abnormal phenotypes. Results: Overall, rare missense mutations in SIX1, KDR/VEGFR2, and PDGFRA showed the best evidence of segregation with the OAV phenotypes in this family. When considering affection with any of the 3 OAVS phenotypes as an outcome, parent-TDTs and sib-TDTs (unadjusted p-values) found that SIX1 (p=0.025, p=0.052), followed by PDGFRA (p=0.180, p=0.069) and KDR/VEGFR2 (p=0.180, p=0.069) have the strongest associations in this family. Burden analysis via a penalized linear mixed model identified SIX1 (RC=0.87) and PDGFRA (RC=0.98) as having the strongest association with OAVS severity. Using phenotype-specific ogfrautcomes, sib-TDTs identified associations between (1) SIX1 with uni- or bilateral ptosis (p=0.049) and ear tags (p=0.01), (2) PDGFRA and KDR/VEGFR2 with ear tags (both p<0.01). Conclusion: Our study reports the genomic findings of a large family with multiple individuals affected with OAVS phenotypes with autosomal dominant inheritance. Our findings narrow down to three potential candidate genes, SIX1, PDGFRA, and KDR/VEGFR2. Among these, SIX1 has been previously associated with OAVS ear malformations and it is co-expressed with EYA1 during ear development. Attempts to strengthen the genotype-phenotype co-relation underlying the OAVS of phenotypes are essential to discover the etiological factors leading to this complex and burdensome condition as well as for family counseling and prevention efforts.

DNA methylation differences in monozygotic twins with Van der Woude syndrome.

Introduction: Van der Woude Syndrome (VWS) is an autosomal dominant disorder responsible for 2% of all syndromic orofacial clefts (OFCs) with IRF6 being the primary causal gene (70%). Cases may present with lip pits and either cleft lip, cleft lip with cleft palate, or cleft palate, with marked phenotypic discordance even among individuals carrying the same mutation. This suggests that genetic or epigenetic modifiers may play additional roles in the syndrome's etiology and variability in expression. We report the first DNA methylation profiling of 2 pairs of monozygotic twins with VWS. Our goal is to explore epigenetic contributions to VWS etiology and variable phenotypic expressivity by comparing DNAm profiles in both twin pairs. While the mutations that cause VWS in these twins are known, the additional mechanism behind their phenotypic risk and variability in expression remains unclear. Methods: We generated whole genome DNAm data for both twin pairs. Differentially methylated positions (DMPs) were selected based on: (1) a coefficient of variation in DNAm levels in unaffected individuals < 20%, and (2) intra-twin pair absolute difference in DNAm levels >5% (delta beta > | 0.05|). We then divided the DMPs in two subgroups for each twin pair for further analysis: (1) higher methylation levels in twin A (Twin A > Twin B); and (2) higher methylation levels in twin B (Twin B >Twin A). Results and Discussion: Gene ontology analysis revealed a list of enriched genes that showed significant differential DNAm, including clef-associated genes. Among the cleft-associated genes, TP63 was the most significant hit (p=7.82E-12). Both twin pairs presented differential DNAm levels in CpG sites in/near TP63 (Twin 1A > Twin 1B and Twin 2A < Twin 2B). The genes TP63 and IRF6 function in a biological regulatory loop to coordinate epithelial proliferation and differentiation in a process that is critical for palatal fusion. The effects of the causal mutations in IRF6 can be further impacted by epigenetic dysregulation of IRF6 itself, or genes in its pathway. Our data shows evidence that changes in DNAm is a plausible mechanism that can lead to markedly distinct phenotypes, even among individuals carrying the same mutation.

Rare and Common Variants Conferring Risk of Tooth Agenesis.

We present association results from a large genome-wide association study of tooth agenesis (TA) as well as selective TA, including 1,944 subjects with congenitally missing teeth, excluding third molars, and 338,554 controls, all of European ancestry. We also tested the association of previously identified risk variants, for timing of tooth eruption and orofacial clefts, with TA. We report associations between TA and 9 novel risk variants. Five of these variants associate with selective TA, including a variant conferring risk of orofacial clefts. These results contribute to a deeper understanding of the genetic architecture of tooth development and disease. The few variants previously associated with TA were uncovered through candidate gene studies guided by mouse knockouts. Knowing the etiology and clinical features of TA is important for planning oral rehabilitation that often involves an interdisciplinary approach.

A Population-Based Study of Effects of Genetic Loci on Orofacial Clefts.

Prior genome-wide association studies for oral clefts have focused on clinic-based samples with unclear generalizability. Prior samples were also small for investigating effects by cleft type and exclusively studied isolated clefts (those occurring without other birth defects). We estimated the effects of 17 top loci on cleft types in both isolated and nonisolated cases in the largest consortium to date of European-descent population-based studies. Our analytic approach focused on a mother-child dyad case-control design, but it also allowed analyzing mother-only or child-only genotypes to maximize power. Our total sample included 1,875 cases with isolated clefts, 459 cases with nonisolated clefts, and 3,749 controls. After correcting for multiple testing, we observed significant associations between fetal single-nucleotide polymorphisms (SNPs) at IRF6, PAX7, 8q21.3, 8q24, KIAA1598-VAX1, and MAFB and isolated cleft lip only (CLO) and cleft lip and palate (CLP). Significant associations were observed between isolated CLO and fetal SNPs near TPM1 and NOG1 and between CLP and fetal SNPs at ABCA4-ARHGAP29, THADA, FOXE1, and SPRY2. Overall, effects were similar for isolated CLO and CLP, except for ABCA4-ARHGAP29. A protective effect was observed for the fetal NOG1 SNP on cleft palate only, opposite in direction to the effect on CLO. For most fetal SNPs, a dose-response allelic effect was observed. No evidence of parent-of-origin or maternal genome effects was observed. Overall, effect direction and magnitude were similar between isolated and nonisolated clefts, suggesting that several loci are modifiers of cleft risk in both isolated and nonisolated forms. Our results provide reliable estimates of the effects of top loci on risks of oral clefts in a population of European descent.

Dental Decay Phenotype in Nonsyndromic Orofacial Clefting.

Although children with oral clefts have a higher risk for dental anomalies when compared with the general population, prior studies have shown conflicting results regarding their dental decay risk. Also, few studies have assessed dental decay risk in unaffected relatives of children with clefts. Thus, the question of increased risk of dental decay in individuals with oral clefts or their unaffected relatives is still open for empirical investigation. This study characterizes dental decay in the largest international cohort to date of children with nonsyndromic clefts and their relatives, as compared with controls, and it addresses whether families with oral clefts have a significantly increased risk for dental decay versus the general population. A total of 3,326 subjects were included: 639 case probands, 1,549 unaffected relatives, and 1,138 controls. Decay was identified from in-person dental examinations or intraoral photographs. Case-control differences were tested with regression analysis. No significant differences were shown in percentage decayed and filled teeth and decayed teeth in the primary dentition (dft, dt) and permanent dentition (DFT, DT) in cases versus controls. In the cleft region, no significant differences were seen in primary or permanent decay (dt, DT) when compared with controls. No difference was found with regard to cleft type and percentage dft, dt, DFT, and DT in case probands. Nonsignificant differences were found in unaffected siblings and parents versus controls (primary and permanent dentitions). Collectively, these findings indicate that individuals with nonsyndromic oral clefts and their families do not have a higher dental decay risk as compared with the general population. These results suggest that either genetic or environmental factors underlying a higher susceptibility for dental anomalies do not increase caries risk or that the seemingly higher risk for dental decay associated with increased dental anomalies in case probands may be superseded by possible greater access to dental care.

Psychosocial well-being of parents of children with oral clefts.

BACKGROUND: Parents of children with oral clefts may be impacted psychosocially in several ways, but empirical evidence remains relatively sparse. The aim of this study was to identify predictors of psychosocial well-being of parents of affected children. METHODS: The study included a total sample of 287 parents (171 mothers and 116 fathers) of children with oral clefts. Parents completed validated psychosocial instruments to measure social avoidance and distress, fear of negative evaluation scale, self-esteem and interpersonal support. Regression analysis was used to evaluate how selected child, parent and household characteristics relate to psychosocial outcomes focusing on child's cleft type, sex and age, differences between mothers and fathers, marital status and household income. RESULTS: Fathers had higher self-esteem than mothers (P = 0.01) and lower concern of being negatively judged by others (P < 0.0001) but also had lower perception of having someone to talk to about their problems (P = 0.01). High household income was associated with greater self-esteem and perception of social support (<0.05). Parents of male affected children had greater perception of social support than parents of female affected children (P = 0.04). No significant differences in parental psychosocial status measures were found by cleft type. Similarly, there is little evidence of changes with child age, except for an increase in parental distress and decline in self-esteem during mid-adolescence (age 15-17 years). CONCLUSIONS: The results indicate that mothers and fathers of children with oral clefts may differ in their psychosocial adjustment and that mothers may overall experience more psychosocial problems than fathers. Also, parents from less wealthy households may be at greater risk. Parental psychosocial status should be considered in holistic family-based treatment approaches to reduce burden on affected families and improve their well-being.

Spectrum of Dental Phenotypes in Nonsyndromic Orofacial Clefting.

Children with oral clefts show a wide range of dental anomalies, adding complexity to understanding the phenotypic spectrum of orofacial clefting. The evidence is mixed, however, on whether the prevalence of dental anomalies is elevated in unaffected relatives and is mostly based on small samples. In the largest international cohort to date of children with nonsyndromic clefts, their relatives, and controls, this study characterizes the spectrum of cleft-related dental anomalies and evaluates whether families with clefting have a significantly higher risk for such anomalies compared with the general population. A total of 3,811 individuals were included: 660 cases with clefts, 1,922 unaffected relatives, and 1,229 controls. Dental anomalies were identified from in-person dental exams or intraoral photographs, and case-control differences were tested using χ(2) statistics. Cases had higher rates of dental anomalies in the maxillary arch than did controls for primary (21% vs. 4%, P = 3 × 10(-8)) and permanent dentitions (51% vs. 8%, P = 4 × 10(-62)) but not in the mandible. Dental anomalies were more prevalent in cleft lip with cleft palate than other cleft types. More anomalies were seen in the ipsilateral side of the cleft. Agenesis and tooth displacements were the most common dental anomalies found in case probands for primary and permanent dentitions. Compared with controls, unaffected siblings (10% vs. 2%, P = 0.003) and parents (13% vs. 7%, P = 0.001) showed a trend for increased anomalies of the maxillary permanent dentition. Yet, these differences were nonsignificant after multiple-testing correction, suggesting genetic heterogeneity in some families carrying susceptibility to both overt clefts and dental anomalies. Collectively, the findings suggest that most affected families do not have higher genetic risk for dental anomalies than the general population and that the higher prevalence of anomalies in cases is primarily a physical consequence of the cleft and surgical interventions.

Phenotype-genotype correlations of facial width and height proportions in patients with Class II malocclusion.

OBJECTIVES: To characterize soft-tissue facial height and width variation in Class II malocclusion and test for correlations with genes HMGA2, AJUBA, and ADK. SETTING AND SAMPLE POPULATION: Nine facial proportions were estimated from 2D frontal repose photographs of 330 Caucasian adults with Class II malocclusion. MATERIAL AND METHODS: After adjustments for age and gender, the facial proportions were submitted to a principal component analyses (PCA). The most meaningful phenotypic variations were correlated with SNPs rs7924176 (ADK), rs17101923 (HMGA2), and rs997154 (AJUBA) genotyped in 106 individuals. RESULTS: Principal component analyses resulted in four principal components (PCs), which explained 75% of total variation. PC1 captured variation in the intercanthus distance and explained 28% of total variation. PC2 explained 21% of the variations in facial taper and facial index. PC3 explained 14% and reflected variations in the vertical dimension of the lower face. PC4 explained 12% and captured variations in distance between the eyes, width of the commissures, and the length of the superior aspect of the lower face height corresponding to the vertical dimension of the philtrum of the upper lip. A suggestive association (p<0.05) was observed between PC4 and rs997154 corroborating the role of AJUBA in variation of facial dimensions. CONCLUSION: 2D frontal photographs can be used to derive quantitative measures of soft-tissue phenotypes that are of clinical relevance. The methods described are suitable for discovery and replication of associations between genotypes and malocclusion phenotypes.

Genetics of the dentofacial variation in human malocclusion.

Malocclusions affect individuals worldwide, resulting in compromised function and esthetics. Understanding the etiological factors contributing to the variation in dentofacial morphology associated with malocclusions is the key to develop novel treatment approaches. Advances in dentofacial phenotyping, which is the comprehensive characterization of hard and soft tissue variation in the craniofacial complex, together with the acquisition of large-scale genomic data have started to unravel genetic mechanisms underlying facial variation. Knowledge on the genetics of human malocclusion is limited even though results attained thus far are encouraging, with promising opportunities for future research. This review summarizes the most common dentofacial variations associated with malocclusions and reviews the current knowledge of the roles of genes in the development of malocclusions. Lastly, this review will describe ways to advance malocclusion research, following examples from the expanding fields of phenomics and genomic medicine, which aim to better patient outcomes.

Candidate Gene Analyses of Skeletal Variation in Malocclusion.

This study evaluated associations between craniofacial candidate genes and skeletal variation in patients with malocclusion. Lateral cephalometric radiographs of 269 untreated adults with skeletal classes I, II, and III malocclusion were digitized with 14 landmarks. Two-dimensional coordinates were analyzed using Procrustes fit and principal component (PC) analysis to generate continuous malocclusion phenotypes. Skeletal class classifications (I, II, or III) were used as a categorical phenotype. Individuals were genotyped for 198 single-nucleotide polymorphisms (SNPs) in 71 craniofacial genes and loci. Phenotype-genotype associations were tested via multivariate linear regression for continuous phenotypes and multinomial logistic regression for skeletal malocclusion class. PC analysis resulted in 4 principal components (PCs) explaining 69% of the total skeletal facial variation. PC1 explained 32.7% of the variation and depicted vertical discrepancies ranging from skeletal deep to open bites. PC1 was associated with a SNP near PAX5 (P = 0.01). PC2 explained 21.7% and captured horizontal maxillomandibular discrepancies. PC2 was associated with SNPs upstream of SNAI3 (P = 0.0002) and MYO1H (P = 0.006). PC3 explained 8.2% and captured variation in ramus height, body length, and anterior cranial base orientation. PC3 was associated with TWIST1 (P = 0.000076). Finally, PC4 explained 6.6% and detected variation in condylar inclination as well as symphysis projection. PC4 was associated with PAX7 (P = 0.007). Furthermore, skeletal class II risk increased relative to class I with the minor alleles of SNPs in FGFR2 (odds ratio [OR] = 2.1, P = 0.004) and declined with SNPs in EDN1 (OR = 0.5, P = 0.007). Conversely, skeletal class III risk increased versus class I with SNPs in FGFR2 (OR 2.2, P = 0.005) and COL1A1 (OR = 2.1, P = 0.008) and declined with SNPs in TBX5 (OR = 0.5, P = 0.014). PAX5, SNAI3, MYO1H, TWIST1, and PAX7 are associated with craniofacial skeletal variation among patients with malocclusion, while FGFR2, EDN1, TBX5, and COL1A1 are associated with type of skeletal malocclusion.

Comparison of observer reliability of three-dimensional cephalometric landmark identification on subject images from Galileos and i-CAT cone beam CT.

OBJECTIVES: Recently, there has been increasing interest in the use of cone beam CT (CBCT) for three-dimensional cephalometric analysis and craniofacial reconstruction in orthodontic and orthognathic surgical treatment planning. However, there is a need to redefine the cephalometric landmarks in three dimensional cephalometric analysis and to demonstrate the reproducibility of landmark identification on the type of CBCT machine being used. METHODS: CBCT images of 20 subjects aged 15-25 years were selected, ten each from Galileos(®) (Sirona Dental Systems Inc., Bensheim, Germany) and Next Generation i-CAT(®) (Imaging Sciences International, Hatfield, PA). 2 observers located 18 landmarks on each subject twice using Dolphin-3D v. 11 software (Dolphin Imaging and Management Systems, Chatsworth, CA). Inter- and intraobserver reliability was assessed using Euclidean distances and linear mixed models. RESULTS: Overall, the intra- and interobserver reliability was excellent for both machines. The landmarks Gonion, Nasion, Orbitale and Anterior Nasal Spine (ANS) showed the greatest median Euclidean distances for both intra- and interobserver measurements. There were significant observer effects in the unified models for Sella, Menton and all six dental landmarks. For Sella, the distances between the measures were significantly smaller (more closely spaced) on the i-CAT machine than on the Galileos in both intra- and interobserver measurements. CONCLUSIONS: The intra- and interobserver reliability was excellent for both machines. Some of the landmarks were not as reproducible as others. Which machine produced the highest reliability depended on the landmark considered.