Six NSCL/P Loci Show Associations With Normal-Range Craniofacial Variation

A common cis-regulatory variant impacts normal-range and disease-associated human facial shape through regulation of PKDCC during chondrogenesis

Genome-wide association studies (GWAS) identified thousands of genetic variants linked to phenotypic traits and disease risk. However, mechanistic understanding of how GWAS variants influence complex morphological traits and can, in certain cases, simultaneously confer normal-range phenotypic variation and disease predisposition, is still largely lacking. Here, we focus on rs6740960, a single nucleotide polymorphism (SNP) at the 2p21 locus, which in GWAS studies has been associated both with normal-range variation in jaw shape and with an increased risk of non-syndromic orofacial clefting. Using in vitro derived embryonic cell types relevant for human facial morphogenesis, we show that this SNP resides in an enhancer that regulates chondrocytic expression of PKDCC - a gene encoding a tyrosine kinase involved in chondrogenesis and skeletal development. In agreement, we demonstrate that the rs6740960 SNP is sufficient to confer chondrocyte-specific differences in PKDCC expression. By deploying dense landmark morphometric analysis of skull elements in mice, we show that changes in Pkdcc dosage are associated with quantitative changes in the maxilla, mandible, and palatine bone shape that are concordant with the facial phenotypes and disease predisposition seen in humans. We further demonstrate that the frequency of the rs6740960 variant strongly deviated among different human populations, and that the activity of its cognate enhancer diverged in hominids. Our study provides a mechanistic explanation of how a common SNP can mediate normal-range and disease-associated morphological variation, with implications for the evolution of human facial features.

A comprehensive consolidation of data on the relationship between IRF6 polymorphisms and non-syndromic cleft lip/palate susceptibility: From 79 case-control studies

BACKGROUND

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is a prevalent craniofacial birth defect on a global scale. A number of candidate genes have been identified as having an impact on NSCL/P. However, the association between interferon regulatory factor 6 (IRF6) polymorphisms and NSCL/P has yielded inconsistent results, prompting the need for a meta-analysis to obtain more accurate estimates.

METHODS

We conducted a thorough screening of all relevant articles published up until November 15, 2023, in online bibliographic databases. The statistical analysis of the collected data was performed using the Comprehensive Meta-Analysis (Version 4.0) software.

RESULTS

A total of 79 case-control studies, comprising 14,003 cases and 19,905 controls, were included in our analysis. The combined data indicated that the IRF6 rs642961 and rs2235371 polymorphisms were associated with an increased risk of NSCL/P in the overall population. However, no significant association was found between the rs2013162 and rs2235375 polymorphisms and the risk of NSCL/P in the overall population. Furthermore, subgroup analyses revealed significant correlations between the IRF6 rs642961, rs2235371, and rs2235375 polymorphisms and the risk of NSCL/P based on ethnic background and country of origin. Nevertheless, the rs2013162 polymorphism plays a protective role in Caucasians and mixed populations.

CONCLUSIONS

Our collective data indicates a significant association between the rs642961 and rs2235371 polymorphisms and the risk of NSCL/P in the overall population. The rs2235375 polymorphism could influence the susceptibility to NSCL/P based on ethnic background. Meanwhile, the rs2013162 polymorphism provides protective effects in Caucasian, mixed populations, and the Brazilian population.

Rs28446116 in PTCH1 is associated with non-syndromic cleft lip with or without palate in the Ningxia population, China

OBJECTIVES

To investigate the association between PTCH1 single nucleotide polymorphism(SNP) and non-syndromic cleft lip with or without palate (NSCL/P) in the Ningxia Hui Autonomous region and predict the function of single nucleotide polymorphism through bioinformatics analysis.

DESIGN

A case-control analysis of 31 single nucleotide polymorphism locus alleles on PTCH1 gene (504 cases and 455 controls) was performed to explore the association between PTCH1 gene polymorphisms and non-syndromic cleft lip with or without palate in Ningxia region. Transcription factors, 3D single nucleotide polymorphism and other related information of single nucleotide polymorphism loci with statistical significance were screened by the case-control experiments, and then analyzed the corresponding transcription factors through the NCBI database.

RESULTS

The case-control study showed that 5 of the 31 single nucleotide polymorphism loci rs357564 (P = 0.0233), rs1805155 (P = 0.0371), rs28446116 (P = 0.0408), rs2282041 (P = 0.0439), rs56119276 (P = 0.0256) had statistically significant differences in allele frequencies between the case and control groups. Bioinformatics analysis revealed that EP300 and RUNX3, among the transcription factors associated with rs28446116, may be associated with the development of non-syndromic cleft lip with or without palate.

CONCLUSION

PTCH1 gene may be associated with the occurrence of non-syndromic cleft lip with or without palate in the Ningxia region, which may be related to the role of EP300 and RUNX3 in the development of cleft lip and palate.

The heritability of vocal tract structures estimated from structural MRI in a large cohort of Dutch twins

While language is expressed in multiple modalities, including sign, writing, or whistles, speech is arguably the most common. The human vocal tract is capable of producing the bewildering diversity of the 7000 or so currently spoken languages, but relatively little is known about its genetic bases, especially in what concerns normal variation. Here, we capitalize on five cohorts totaling 632 Dutch twins with structural magnetic resonance imaging (MRI) data. Two raters placed clearly defined (semi)landmarks on each MRI scan, from which we derived 146 measures capturing the dimensions and shape of various vocal tract structures, but also aspects of the head and face. We used Genetic Covariance Structure Modeling to estimate the additive genetic, common environmental or non-additive genetic, and unique environmental components, while controlling for various confounds and for any systematic differences between the two raters. We found high heritability, h², for aspects of the skull and face, the mandible, the anteroposterior (horizontal) dimension of the vocal tract, and the position of the hyoid bone. These findings extend the existing literature, and open new perspectives for understanding the complex interplay between genetics, environment, and culture that shape our vocal tracts, and which may help explain cross-linguistic differences in phonetics and phonology.

What’s Shape Got to Do With It? Examining the Relationship Between Facial Shape and Orofacial Clefting

Nonsyndromic orofacial clefts belong to a class of congenital malformations characterized by a complex and multifactorial etiology. During early facial development, multiple factors can disrupt fusion leading to a cleft; this includes the shape of the embryonic face. The face shape hypothesis (FSH) of orofacial clefting emerged in the 1960s, influenced by morphological differences observed within affected families, comparative studies of mouse models, and advances in modeling genetic liability for complex traits in populations. For the past five decades, studies have documented changes in the shape or spatial arrangement of facial prominences in embryonic mice and altered post-natal facial shape in individuals at elevated risk for orofacial clefting due to their family history. Moreover, recent studies showing how genes that impact facial shape in humans and mice are providing clues about the genetic basis of orofacial clefting. In this review, I discuss the origins of the FSH, provide an overview of the supporting evidence, and discuss ways in which the FSH can inform our understanding of orofacial clefting.

Decoding the Human Face: Challenges and Progress in Understanding the Genetics of Craniofacial Morphology

Variations in the form of the human face, which plays a role in our individual identities and societal interactions, have fascinated scientists and artists alike. Here, we review our current understanding of the genetics underlying variation in craniofacial morphology and disease-associated dysmorphology, synthesizing decades of progress on Mendelian syndromes in addition to more recent results from genome-wide association studies of human facial shape and disease risk. We also discuss the various approaches used to phenotype and quantify facial shape, which are of particular importance due to the complex, multipartite nature of the craniofacial form. We close by discussing how experimental studies have contributed and will further contribute to our understanding of human genetic variation and then proposing future directions and applications for the field.

Evaluation of genome-wide association signals for nonsyndromic cleft lip with or without cleft palate in a multiethnic Brazilian population

OBJECTIVE

To evaluate previous nonsyndromic cleft lip with or without cleft palate (NSCL±P) associated signals in 4p16.2, 8p11.23, 12q13.13, 12q13.2 and 17q21.32 in a multiethnic Brazilian cohort.

DESIGN

The single nucleotide polymorphisms (SNPs) rs34246903 in 4p16.2, rs13317 in 8p11.23 (FGFR1, fibroblast growth factor receptor 1), rs3741442 in 12q13.13, rs705704 in 12q13.2 and rs4968247 in 17q21.32 were genotyped with TaqMan allelic discrimination assays in a case-control sample including 801 NSCL±P patients [233 nonsyndromic cleft lip (NSCLO) and 568 nonsyndromic cleft lip and palate (NSCLP)] and 881 healthy controls. Multiple logistic regression analyses, considering sex and genomic ancestry as covariates, were conducted, and the p value was adjusted with Bonferroni multiple correction testing (p ≤ 0.01).

RESULTS

Although several associations were identified, those that resisted the multiple correction testing involved the alleles and genotypes of rs34246903 and rs13317. The NSCLO group had a lower frequency of the minor C allele of rs34246903 compared to controls, giving an odds ratio (OR) of 0.74 [95% confidence interval (CI): 0.59-0.93, p = 0.01]. The rs34246903 CC genotype (homozygous) and the recessive model revealed significant protective associations with NSCLO, yielding ORs of 0.50 (95% CI: 0.29-0.85, p = 0.005) and 0.55 (95% CI: 0.33-0.93, p = 0.01) respectively. The presence of C variant allele of rs13317 (OR: 0.81, 95% CI: 0.69-0.96, p = 0.01) as well the TC genotype (OR: 0.77, 95% CI: 0.62-0.94, p = 0.01) and the dominant model (OR: 0.77, 95% CI: 0.63-0.94, p = 0.009) showed significant associations with reduced risk of NSCL±P.

CONCLUSION

Our study is the first to support the association of rs34246903 (4p16.2) with NSCLO and rs13317 within FGFR1 with NSCL±P in the highly admixed Brazilian population. Further studies are needed to determine the functionality of those SNPs or to identify the causal markers in linkage disequilibrium with those susceptibility markers.

Association between PTCH1 and RAD54B Single-Nucleotide Polymorphisms and Non-syndromic Orofacial Clefts in the Northeast Population of Iran

Background

Non-Syndromic Cleft Lip with or without cleft Palate (NSCL/P) is a common developmental disorder of the head and neck with a multifactorial etiology. The current study aimed to evaluate the potential association of PTCH1 (rs10512248) and RAD54B (rs12681366) polymorphisms with NSCL/P in the Northeast Iranian population.

Methods

In the present study, blood samples were taken from 122 subjects with NSCL/P and 161 healthy controls. Polymerase Chain Reaction (PCR) followed by Restriction Fragment Length Polymorphism (RFLP) were used to conduct genotyping of single-nucleotide polymorphisms.

Results

Although differences were observed between cases and controls in rs10512248 and rs12681366, our data did not support a significant association of these polymorphisms with NSCL/P in our population.

Conclusion

Our findings suggest that polymorphisms of rs10512248 and rs12681366 may not be potential risk factors for NSCL/P in the Northeast Iranian population due to the multifactorial and multiethnicity characteristics of some genes.

Genetic architecture of orbital telorism

The interocular distance, or orbital telorism, is a distinctive craniofacial trait that also serves as a clinically informative measure. While its extremes, hypo- and hypertelorism, have been linked to monogenic disorders and are often syndromic, little is known about the genetic determinants of interocular distance within the general population. We derived orbital telorism measures from cranial magnetic resonance imaging by calculating the distance between the eyeballs’ centre of gravity, which showed a good reproducibility with an intraclass correlation coefficient of 0.991 (95% confidence interval 0.985–0.994). Heritability estimates were 76% (standard error = 12%) with a family-based method (N = 364) and 39% (standard error = 2.4%) with a single nucleotide polymorphism-based method (N = 34 130) and were unaffected by adjustment for height (model II) and intracranial volume (model III) or head width (model IV). Genome-wide association studies in 34 130 European individuals identified 56 significantly associated genomic loci (P < 5 × 10⁻⁸) across four different models of which 46 were novel for facial morphology, and overall these findings replicated in an independent sample (N = 10 115) with telorism-related horizontal facial distance measures. Genes located nearby these 56 identified genetic loci were 4.9-fold enriched for Mendelian hypotelorism and hypertelorism genes, underlining their biological relevance. This study provides novel insights into the genetic architecture underlying interocular distance in particular, and the face in general, and explores its potential for applications in a clinical setting.

Genome scans of facial features in East Africans and cross-population comparisons reveal novel associations

Facial morphology is highly variable, both within and among human populations, and a sizable portion of this variation is attributable to genetics. Previous genome scans have revealed more than 100 genetic loci associated with different aspects of normal-range facial variation. Most of these loci have been detected in Europeans, with few studies focusing on other ancestral groups. Consequently, the degree to which facial traits share a common genetic basis across diverse sets of humans remains largely unknown. We therefore investigated the genetic basis of facial morphology in an East African cohort. We applied an open-ended data-driven phenotyping approach to a sample of 2,595 3D facial images collected on Tanzanian children. This approach segments the face into hierarchically arranged, multivariate features that capture the shape variation after adjusting for age, sex, height, weight, facial size and population stratification. Genome scans of these multivariate shape phenotypes revealed significant (p < 2.5 × 10⁻⁸) signals at 20 loci, which were enriched for active chromatin elements in human cranial neural crest cells and embryonic craniofacial tissue, consistent with an early developmental origin of the facial variation. Two of these associations were in highly conserved regions showing craniofacial-specific enhancer activity during embryological development (5q31.1 and 12q21.31). Six of the 20 loci surpassed a stricter threshold accounting for multiple phenotypes with study-wide significance (p < 6.25 × 10⁻¹⁰). Cross-population comparisons indicated 10 association signals were shared with Europeans (seven sharing the same associated SNP), and facilitated fine-mapping of causal variants at previously reported loci. Taken together, these results may point to both shared and population-specific components to the genetic architecture of facial variation.

Genetic factors play an important role in shaping human facial features. Over the last decade, studies have identified numerous genes associated with various facial traits. The vast majority of these studies have focused on European or Asian populations, while African populations have been underrepresented. Increasing the diversity of these analyses can reveal novel associations and cross-population analyses can help deepen our understanding of known genetic associations. We therefore performed a genome scan of 3D facial features in African children from Tanzania and then compared our results to Europeans. We found 20 regions of the genome associated with facial shape in Tanzanian children, 10 of which were also present in Europeans, indicating evidence for a partly shared genetic basis for human facial shape across populations. In addition, about half of the genetic associations observed in Tanzanians were not present in Europeans, and some of the shared signals differed between populations in the specific genetic variants associated or specific facial traits affected. These results shed light on the shared and population-specific genetic contributors to normal-range facial variation.

The Intersection of the Genetic Architectures of Orofacial Clefts and Normal Facial Variation

Unaffected relatives of individuals with non-syndromic cleft lip with or without cleft palate (NSCL/P) show distinctive facial features. The presence of this facial endophenotype is potentially an expression of underlying genetic susceptibility to NSCL/P in the larger unselected population. To explore this hypothesis, we first partitioned the face into 63 partially overlapping regions representing global-to-local facial morphology and then defined endophenotypic traits by contrasting the 3D facial images from 264 unaffected parents of individuals with NSCL/P versus 3,171 controls. We observed distinct facial features between parents and controls across 59 global-to-local facial segments at nominal significance (p ≤ 0.05) and 52 segments at Bonferroni corrected significance (p < 1.2 × 10-3), respectively. Next, we quantified these distinct facial features as univariate traits in another dataset of 8,246 unaffected European individuals and performed a genome-wide association study. We identified 29 independent genetic loci that were associated (p < 5 × 10-8) with at least one of the tested endophenotypic traits, and nine genetic loci also passed the study-wide threshold (p < 8.47 × 10-10). Of the 29 loci, 22 were in proximity of loci previously associated with normal facial variation, 18 were near genes that show strong evidence in orofacial clefting (OFC), and another 10 showed some evidence in OFC. Additionally, polygenic risk scores for NSCL/P showed associations with the endophenotypic traits. This study thus supports the hypothesis of a shared genetic architecture of normal facial development and OFC.

A GWAS in Latin Americans identifies novel face shape loci, implicating VPS13B and a Denisovan introgressed region in facial variation

To characterize the genetic basis of facial features in Latin Americans, we performed a genome-wide association study (GWAS) of more than 6000 individuals using 59 landmark-based measurements from two-dimensional profile photographs and ~9,000,000 genotyped or imputed single-nucleotide polymorphisms. We detected significant association of 32 traits with at least 1 (and up to 6) of 32 different genomic regions, more than doubling the number of robustly associated face morphology loci reported until now (from 11 to 23). These GWAS hits are strongly enriched in regulatory sequences active specifically during craniofacial development. The associated region in 1p12 includes a tract of archaic adaptive introgression, with a Denisovan haplotype common in Native Americans affecting particularly lip thickness. Among the nine previously unidentified face morphology loci we identified is the VPS13B gene region, and we show that variants in this region also affect midfacial morphology in mice.

Msx1 deficiency interacts with hypoxia and induces a morphogenetic regulation during mouse lip development

Nonsyndromic clefts of the lip and palate are common birth defects resulting from gene-gene and gene-environment interactions. Mutations in human MSX1 have been linked to orofacial clefting and we show here that Msx1 deficiency causes a growth defect of the medial nasal process (Mnp) in mouse embryos. Although this defect alone does not disrupt lip formation, Msx1-deficient embryos develop a cleft lip when the mother is transiently exposed to reduced oxygen levels or to phenytoin, a drug known to cause embryonic hypoxia. In the absence of interacting environmental factors, the Mnp growth defect caused by Msx1 deficiency is modified by a Pax9-dependent 'morphogenetic regulation', which modulates Mnp shape, rescues lip formation and involves a localized abrogation of Bmp4-mediated repression of Pax9 Analyses of GWAS data revealed a genome-wide significant association of a Gene Ontology morphogenesis term (including assigned roles for MSX1, MSX2, PAX9, BMP4 and GREM1) specifically for nonsyndromic cleft lip with cleft palate. Our data indicate that MSX1 mutations could increase the risk for cleft lip formation by interacting with an impaired morphogenetic regulation that adjusts Mnp shape, or through interactions that inhibit Mnp growth.

Sources of variation in the 3dMDface and Vectra H1 3D facial imaging systems

As technology advances and collaborations grow, our ability to finely quantify and explore morphological variation in 3D structures can enable important discoveries and insights into clinical, evolutionary, and genetic questions. However, it is critical to explore and understand the relative contribution of potential sources of error to the structures under study. In this study, we isolated the level of error in 3D facial images attributable to four sources, using the 3dMDface and Vectra H1 camera systems. When the two camera systems are used separately to image human participants, this analysis finds an upper bound of error potentially introduced by the use of the 3dMDface or Vectra H1 camera systems, in conjunction with the MeshMonk registration toolbox, at 0.44 mm and 0.40 mm, respectively. For studies using both camera systems, this upper bound increases to 0.85 mm, on average, and there are systematic differences in the representation of the eyelids, nostrils, and mouth by the two camera systems. Our results highlight the need for careful assessment of potential sources of error in 3D images, both in terms of magnitude and position, especially when dealing with very small measurements or performing many tests.

Craniofacial Analysis May Indicate Co-Occurrence of Skeletal Malocclusions and Associated Risks in Development of Cleft Lip and Palate

Non-syndromic orofacial clefts encompass a range of morphological changes affecting the oral cavity and the craniofacial skeleton, of which the genetic and epigenetic etiologic factors remain largely unknown. The objective of this study is to explore the contribution of underlying dentofacial deformities (also known as skeletal malocclusions) in the craniofacial morphology of non-syndromic cleft lip and palate patients (nsCLP). For that purpose, geometric morphometric analysis was performed using full skull cone beam computed tomography (CBCT) images of patients with nsCLP (n = 30), normocephalic controls (n = 60), as well as to sex- and ethnicity- matched patients with an equivalent dentofacial deformity (n = 30). Our outcome measures were shape differences among the groups quantified via principal component analysis and associated principal component loadings, as well as mean shape differences quantified via a Procrustes distance among groups. According to our results, despite the shape differences among all three groups, the nsCLP group shares many morphological similarities in the maxilla and mandible with the dentofacial deformity group. Therefore, the dentoskeletal phenotype in nsCLP could be the result of the cleft and the coexisting dentofacial deformity and not simply the impact of the cleft.

Assessment of candidate genes and genetic heterogeneity in human non syndromic orofacial clefts specifically non syndromic cleft lip with or without palate

Non syndromic orofacial clefts specifically non-syndromic cleft lip/palate are one of the most common craniofacial malformation among birth defects in human having multifactorial etiology with an incidence of 1:700/1000. On the basis of association with other congenital malformations or their presence as isolated anomaly, OFC can be classified as syndromic (30%) and nonsyndromic (70%) respectively. The major cause of disease demonstrates complex interplay between genetic and environmental factors. The pathogenic mechanism of underlying factors have been provided by different genetic studies on large-scale with significant recent advances in genotyping technologies usually based on linkage or genome wide association studies (GWAS). On the basis of recent studies, new tools to identify causative genes involved in NSCL/P reported approximately more than 30 genetic risk loci that are responsible for pathogenesis of facial deformation. Despite these findings, it is still uncertain that how much of variance in NSCL/P predisposing factors can be explain by identified risk loci, as they all together accounts for only 20%-25% of NSCL/P heritability. So there is need of further findings about the problem of rare low frequency coding variants and other missing responsive factors or genetic modifiers. This review will described those potential genes and loci reported in different studies whose involvement in pathogenesis of nonsyndromic OFC has wide scientific evidence.

Polymorphic Variants of V-Maf Musculoaponeurotic Fibrosarcoma Oncogene Homolog B (rs13041247 and rs11696257) and Risk of Non-Syndromic Cleft Lip/Palate: Systematic Review and Meta-Analysis

Background: Non-syndromic cleft lip/palate (NSCL/P) has an etiology, including both genetic and environmental factors. Herein, we evaluated the association of rs13041247 and rs11696257 v-maf musculoaponeurotic fibrosarcoma oncogene homolog B (MAFB) polymorphisms with the risk of NSCL/P in a meta-analysis. Methods: The PubMed/Medline, Scopus, Cochrane Library, Web of Science, and HuGE Navigator databases were systematically searched to retrieve relevant articles published up to January 2019. The Newcastle–Ottawa scale was applied for quality evaluation of retrieved articles. The 95% confidence interval (CI) and crude odds ratio (OR) were calculated for each study using the Review Manager 5.3 software to show the association between MAFB polymorphisms and risk of NSCL/P. The comprehensive meta-analysis 2.0 software was used to calculate the publication bias. In addition, sensitivity analysis was carried out to show the stability of results. Results: Of 102 articles retrieved from the databases, 10 articles were analyzed in this meta-analysis. Ten articles, including eleven studies reporting rs13041247 MAFB polymorphism, included 3082 NSCL/P patients and 4104 controls. Three studies that reported rs11696257 MAFB polymorphism involved 845 NSCL/P patients and 927 controls. The rs11696257 MAFB polymorphism was not associated with the risk of NSCL/P, but the CC and TC genotypes of rs13041247 polymorphism were associated with the risk of NSCL/P. Nevertheless, the C allele and CC and TC genotypes were associated with a significant decline in the risk of NSCL/P in population-based studies. Conclusions: The results of this meta-analysis demonstrated that the risk of NSCL/P was related to rs13041247 polymorphism, not rs11696257 MAFB polymorphism. Well-designed studies are required to assess the interaction of MAFB and other genes with environmental factors in different ethnic groups.