Imputation of orofacial clefting data identifies novel risk loci and sheds light on the genetic background of cleft lip ± cleft palate and cleft palate only

Role of MTHFR, IRF6, PAX7 and TP63 SNPs in susceptibility to non-syndromic orofacial cleft, a candidate gene study in a Portuguese population

BACKGROUND

Non-syndromic orofacial cleft (NSOC) is a complex phenotype, involving multiple genetic and environmental factors. Association studies exploring the genetic susceptibility to this prevalent oral malformation show variability of results in different populations. Using a candidate gene approach, we aimed to verify the role of four single-nucleotide polymorphisms (SNPs) in the susceptibility to NSOC in Portuguese patients.

METHODS

A total of 254 non-consanguineous individuals of Portuguese were recruited, including 120 patients with NSOC and 134 controls. About 92% of these patients had non-syndromic cleft lip with or without cleft palate (NSCL/P) and 8% had only non-syndromic cleft palate (NSCP). SNPs in the MTHFR (rs1801133), IRF6 (rs642961), PAX7 (rs742071) and TP63 (rs9332461) genes were studied, using a real-time approach with TaqMan probes. Allelic, genotypic, dominant, recessive and over-dominant models were explored using a chi-squared test. Adjusted p-value was calculated for multiple comparisons using the Benjamini-Hochberg false discovery rate (FDR).

RESULTS

All SNPs were in Hardy-Weinberg equilibrium. For MTHFR, IRF6, and PAX7 SNPs, no statistically significant difference was highlighted for any of the evaluated models. For TP63 SNP, data fitted an over-dominant model, with a protective effect for heterozygotes (OR 1.897; CI 95% [1.144-3.147]; p < .016, when comparing controls vs. cases), but significance was lost when applying adjusted p-value for multiple comparisons (4 × 5 tests).

CONCLUSION

In this Portuguese population, there was no evidence of an association between the evaluated SNPs and NSOC. For TP63 SNP, the possibility of a protective effect of heterozygotes should be further investigated.

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.

Rare variant modifier analysis identifies variants in SEC24D associated with orofacial cleft subtypes

As one of the most common structural birth defects, orofacial clefts (OFCs) have been studied for decades, and recent studies have demonstrated that there are genetic differences between the different phenotypic presentations of OFCs. However, the contribution of rare genetic variation genome-wide to different subtypes of OFCs has been understudied, with most studies focusing on common genetic variation or rare variation within targeted regions of the genome. Therefore, we used whole-genome sequencing data from the Gabriella Miller Kids First Pediatric Research Program to conduct a gene-based burden analysis to test for genetic modifiers of cleft lip (CL) vs cleft lip and palate (CLP). We found that there was a significantly increased burden of rare variants in SEC24D in CL cases compared to CLP cases (p = 6.86 [Formula: see text] 10-7). Of the 15 variants within SEC24D, 53.3% were synonymous, but overlapped a known craniofacial enhancer. We then tested whether these variants could alter predicted transcription factor binding sites (TFBS), and found that the rare alleles destroyed binding sites for 9 transcription factors (TFs), including Pax1 (p = 0.0009), and created binding sites for 23 TFs, including Pax6 (p = 6.12 [Formula: see text] 10-5) and Pax9 (p = 0.0001), which are known to be involved in normal craniofacial development, suggesting a potential mechanism by which these synonymous variants could have a functional impact. Overall, this study indicates that rare genetic variation may contribute to the phenotypic heterogeneity of OFCs and suggests that regulatory variation may also contribute and warrant further investigation in future studies of genetic variants controlling risk to OFC.

Integrative analysis of transcriptome dynamics during human craniofacial development identifies candidate disease genes

Craniofacial disorders arise in early pregnancy and are one of the most common congenital defects. To fully understand how craniofacial disorders arise, it is essential to characterize gene expression during the patterning of the craniofacial region. To address this, we performed bulk and single-cell RNA-seq on human craniofacial tissue from 4-8 weeks post conception. Comparisons to dozens of other human tissues revealed 239 genes most strongly expressed during craniofacial development. Craniofacial-biased developmental enhancers were enriched +/- 400 kb surrounding these craniofacial-biased genes. Gene co-expression analysis revealed that regulatory hubs are enriched for known disease causing genes and are resistant to mutation in the normal healthy population. Combining transcriptomic and epigenomic data we identified 539 genes likely to contribute to craniofacial disorders. While most have not been previously implicated in craniofacial disorders, we demonstrate this set of genes has increased levels of de novo mutations in orofacial clefting patients warranting further study.

Neural crest E-cadherin loss drives cleft lip/palate by epigenetic modulation via pro-inflammatory gene-environment interaction

Gene-environment interactions are believed to play a role in multifactorial phenotypes, although poorly described mechanistically. Cleft lip/palate (CLP), the most common craniofacial malformation, has been associated with both genetic and environmental factors, with little gene-environment interaction experimentally demonstrated. Here, we study CLP families harbouring CDH1/E-Cadherin variants with incomplete penetrance and we explore the association of pro-inflammatory conditions to CLP. By studying neural crest (NC) from mouse, Xenopus and humans, we show that CLP can be explained by a 2-hit model, where NC migration is impaired by a combination of genetic (CDH1 loss-of-function) and environmental (pro-inflammatory activation) factors, leading to CLP. Finally, using in vivo targeted methylation assays, we demonstrate that CDH1 hypermethylation is the major target of the pro-inflammatory response, and a direct regulator of E-cadherin levels and NC migration. These results unveil a gene-environment interaction during craniofacial development and provide a 2-hit mechanism to explain cleft lip/palate aetiology.

Rare genetic variants in SEC24D modify orofacial cleft phenotypes

As one of the most common structural birth defects, orofacial clefts (OFCs) have been studied for decades, and recent studies have demonstrated that there are genetic differences between the different phenotypic presentations of OFCs. However, the contribution of rare genetic variation genome-wide to different subtypes of OFCs has been understudied, with most studies focusing on common genetic variation or rare variation within targeted regions of the genome. Therefore, we used whole-genome sequencing data from the Gabriella Miller Kids First Pediatric Research Program to conduct a gene-based burden analysis to test for genetic modifiers of cleft lip (CL) vs cleft lip and palate (CLP). We found that there was a significantly increased burden of rare variants in SEC24D in CL cases compared to CLP cases (p=6.86×10-7). Of the 15 variants within SEC24D, 53.3% were synonymous, but overlapped a known craniofacial enhancer. We then tested whether these variants could alter predicted transcription factor binding sites (TFBS), and found that the rare alleles destroyed binding sites for 9 transcription factors (TFs), including Pax1 (p=0.0009), and created binding sites for 23 TFs, including Pax6 (p=6.12×10-5) and Pax9 (p= 0.0001), which are known to be involved in normal craniofacial development, suggesting a potential mechanism by which these synonymous variants could have a functional impact. Overall, this study demonstrates that rare genetic variation contributes to the phenotypic heterogeneity of OFCs and suggests that regulatory variation may also contribute and warrant further investigation in future studies of genetic variants controlling risk to OFC.

Identification of de novo variants in nonsyndromic cleft lip with/without cleft palate patients with low polygenic risk scores

BACKGROUND

Nonsyndromic cleft lip with/without cleft palate (nsCL/P) is a congenital malformation of multifactorial etiology. Research has identified >40 genome-wide significant risk loci, which explain less than 40% of nsCL/P heritability. Studies show that some of the hidden heritability is explained by rare penetrant variants.

METHODS

To identify new candidate genes, we searched for highly penetrant de novo variants (DNVs) in 50 nsCL/P patient/parent-trios with a low polygenic risk for the phenotype (discovery). We prioritized DNV-carrying candidate genes from the discovery for resequencing in independent cohorts of 1010 nsCL/P patients of diverse ethnicities and 1574 population-matched controls (replication). Segregation analyses and rare variant association in the replication cohort, in combination with additional data (genome-wide association data, expression, protein-protein-interactions), were used for final prioritization.

CONCLUSION

In the discovery step, 60 DNVs were identified in 60 genes, including a variant in the established nsCL/P risk gene CDH1. Re-sequencing of 32 prioritized genes led to the identification of 373 rare, likely pathogenic variants. Finally, MDN1 and PAXIP1 were prioritized as top candidates. Our findings demonstrate that DNV detection, including polygenic risk score analysis, is a powerful tool for identifying nsCL/P candidate genes, which can also be applied to other multifactorial congenital malformations.

Prioritization of non-coding elements involved in non-syndromic cleft lip with/without cleft palate through genome-wide analysis of de novo mutations

Non-syndromic cleft lip with/without cleft palate (nsCL/P) is a highly heritable facial disorder. To date, systematic investigations of the contribution of rare variants in non-coding regions to nsCL/P etiology are sparse. Here, we re-analyzed available whole-genome sequence (WGS) data from 211 European case-parent trios with nsCL/P and identified 13,522 de novo mutations (DNMs) in nsCL/P cases, 13,055 of which mapped to non-coding regions. We integrated these data with DNMs from a reference cohort, with results of previous genome-wide association studies (GWASs), and functional and epigenetic datasets of relevance to embryonic facial development. A significant enrichment of nsCL/P DNMs was observed at two GWAS risk loci (4q28.1 (p = 8 × 10-4) and 2p21 (p = 0.02)), suggesting a convergence of both common and rare variants at these loci. We also mapped the DNMs to 810 position weight matrices indicative of transcription factor (TF) binding, and quantified the effect of the allelic changes in silico. This revealed a nominally significant overrepresentation of DNMs (p = 0.037), and a stronger effect on binding strength, for DNMs located in the sequence of the core binding region of the TF Musculin (MSC). Notably, MSC is involved in facial muscle development, together with a set of nsCL/P genes located at GWAS loci. Supported by additional results from single-cell transcriptomic data and molecular binding assays, this suggests that variation in MSC binding sites contributes to nsCL/P etiology. Our study describes a set of approaches that can be applied to increase the added value of WGS data.

Revisiting the embryogenesis of lip and palate development

Clefts of the lip and palate (CLP), the major causes of congenital facial malformation globally, result from failure of fusion of the facial processes during embryogenesis. With a prevalence of 1 in 500-2500 live births, CLP causes major morbidity throughout life as a result of problems with facial appearance, feeding, speaking, obstructive apnoea, hearing and social adjustment and requires complex, multi-disciplinary care at considerable cost to healthcare systems worldwide. Long-term outcomes for affected individuals include increased mortality compared with their unaffected siblings. The frequent occurrence and major healthcare burden imposed by CLP highlight the importance of dissecting the molecular mechanisms driving facial development. Identification of the genetic mutations underlying syndromic forms of CLP, where CLP occurs in association with non-cleft clinical features, allied to developmental studies using appropriate animal models is central to our understanding of the molecular events underlying development of the lip and palate and, ultimately, how these are disturbed in CLP.

Association of genetic polymorphisms of VAX1, MAFB, and NTN1 with nonsyndromic cleft lip with or without cleft palate in Chinese population

Nonsyndromic cleft lip with or without palate (NSCL/P) is a common birth defect involving genetic factors. We conducted this case-control study to verify the association of ten single-nucleotide polymorphisms (SNPs) of six genes (VAX1, MAFB, PAX7, ABCA4, NTN1, and NOG) with NSCL/P in the Chinese population. The study included 249 NSCL/P patients, 62 nonsyndromic cleft palate only (NSCPO) patients and 480 controls. Three loci, namely, VAX1 rs7078160, MAFB rs11696257, and NTN1 rs4791774, were associated with NSCL/P (Bonferroni method adjusted p values were 0.020, 0.00031, and 0.030, respectively). We also found that the disease risk of individuals carrying both VAX1 rs7078160 and NTN1 rs4791774 was higher than those carrying only one of them (p = 4.50 × 10^-4 and 6.03 × 10^-3, respectively). SNPs of genes VAX1 rs7078160, MAFB rs11696257, and NTN1 rs4791774 increased NSCL/P risk in the Chinese population.

Allele-specific transcription factor binding in a cellular model of orofacial clefting

Non-syndromic cleft lip with/without cleft palate (nsCL/P) is a frequent congenital malformation with multifactorial etiology. While recent genome-wide association studies (GWAS) have identified several nsCL/P risk loci, the functional effects of the associated non-coding variants are largely unknown. Furthermore, additional risk loci remain undetected due to lack of power. As genetic variants might alter binding of transcription factors (TF), we here hypothesized that the integration of data from TF binding sites, expression analyses and nsCL/P GWAS might help to (i) identify functionally relevant variants at GWAS loci, and (ii) highlight novel risk variants that have been previously undetected. Analysing the craniofacial TF TFAP2A in human embryonic palatal mesenchyme (HEPM) cells, we identified 2845 TFAP2A ChIP-seq peaks, several of which were located near nsCL/P candidate genes (e.g. MSX1 and SPRY2). Comparison with independent data suggest that 802 of them might be specific to craniofacial development, and genes near these peaks are enriched in processes relevant to nsCL/P. Integration with nsCL/P GWAS data, however, did not show robust evidence for co-localization of common nsCL/P risk variants with TFAP2A ChIP-seq peaks. This data set represents a new resource for the analyses of craniofacial processes, and similar approaches with additional cell lines and TFs could be applied to generate further insights into nsCL/P etiology.

Genetic models and approaches to study orofacial clefts

INTRODUCTION

Orofacial clefts (OFCs) are common craniofacial birth defects with heterogeneous phenotype and etiology. Geneticists have applied nearly every available method and technology to further our understanding of the genetic architectures of OFCs.

OBJECTIVE

This review describes the evidence for a genetic etiology in OFCs, statistical genetic approaches employed to identify genetic causes, and how the results have shaped our current understanding of the genetic architectures of syndromic and nonsyndromic OFCs.

CONCLUSION

There has been rapid progress towards elucidating the genetic architectures of OFCs due to the availability of large collections of DNA samples from cases, controls, and families with OFCs and the consistent adoption of new methodologies and novel statistical approaches as they are developed. Genetic studies have identified rare and common variants influencing risk of OFCs in both Mendelian and complex forms of OFCs, blurring the distinctions traditional categories used in genetic studies and clinical medicine.

Identification of Novel Genomic Variations in Susceptibility to Nonsyndromic Cleft Lip and Palate Patients

BACKGROUND

Nonsyndromic cleft lip with or without palate (NSCL/P) is a multifactorial and common birth malformation caused by genetic and environmental factors, as well as by teratogens. Genome-wide association studies found genetic variations with modulatory effects of NSCL/P formation in Chinese and Iranian populations. We aimed to identify the susceptibility of single-nucleotide polymorphisms (SNPs) to nonsyndromic cleft lip with or without palate in the Indian population.

MATERIAL AND METHODS

The present study was conducted on NSCL/P cases and controls. Genomic DNA was extracted from peripheral blood and Axiom- Precision Medicine Research Array (PMRA) was performed. The Axiom-PMRA covers 902,527 markers and several thousand novel risk variants. Quality control-passed samples were included for candidate genetic variation identification, gene functional enrichment, and pathway and network analysis.

RESULTS

The genome-wide association study identified fourteen novel candidate gene SNPs that showed the most significant association with the risk of NSCL/P, and eight were predicted to have regulatory sequences.

CONCLUSION

The GWAS study showed novel candidate genetic variations in NSCL/P formations. These findings contribute to the understanding of genetic predisposition to nonsyndromic cleft lip with or without palate.

Pleiotropy method reveals genetic overlap between orofacial clefts at multiple novel loci from GWAS of multi-ethnic trios

Based on epidemiologic and embryologic patterns, nonsyndromic orofacial clefts- the most common craniofacial birth defects in humans- are commonly categorized into cleft lip with or without cleft palate (CL/P) and cleft palate alone (CP), which are traditionally considered to be etiologically distinct. However, some evidence of shared genetic risk in IRF6, GRHL3 and ARHGAP29 regions exists; only FOXE1 has been recognized as significantly associated with both CL/P and CP in genome-wide association studies (GWAS). We used a new statistical approach, PLACO (pleiotropic analysis under composite null), on a combined multi-ethnic GWAS of 2,771 CL/P and 611 CP case-parent trios. At the genome-wide significance threshold of 5 × 10-8, PLACO identified 1 locus in 1q32.2 (IRF6) that appears to increase risk for one OFC subgroup but decrease risk for the other. At a suggestive significance threshold of 10-6, we found 5 more loci with compelling candidate genes having opposite effects on CL/P and CP: 1p36.13 (PAX7), 3q29 (DLG1), 4p13 (LIMCH1), 4q21.1 (SHROOM3) and 17q22 (NOG). Additionally, we replicated the recognized shared locus 9q22.33 (FOXE1), and identified 2 loci in 19p13.12 (RAB8A) and 20q12 (MAFB) that appear to influence risk of both CL/P and CP in the same direction. We found locus-specific effects may vary by racial/ethnic group at these regions of genetic overlap, and failed to find evidence of sex-specific differences. We confirmed shared etiology of the two OFC subtypes comprising CL/P, and additionally found suggestive evidence of differences in their pathogenesis at 2 loci of genetic overlap. Our novel findings include 6 new loci of genetic overlap between CL/P and CP; 3 new loci between pairwise OFC subtypes; and 4 loci not previously implicated in OFCs. Our in-silico validation showed PLACO is robust to subtype-specific effects, and can achieve massive power gains over existing approaches for identifying genetic overlap between disease subtypes. In summary, we found suggestive evidence for new genetic regions and confirmed some recognized OFC genes either exerting shared risk or with opposite effects on risk to OFC subtypes.

Integrative approaches generate insights into the architecture of non-syndromic cleft lip with or without cleft palate

Non-syndromic cleft lip with or without cleft palate (nsCL/P) is a common congenital facial malformation with a multifactorial etiology. Genome-wide association studies (GWASs) have identified multiple genetic risk loci. However, functional interpretation of these loci is hampered by the underrepresentation in public resources of systematic functional maps representative of human embryonic facial development. To generate novel insights into the etiology of nsCL/P, we leveraged published GWAS data on nsCL/P as well as available chromatin modification and expression data on mid-facial development. Our analyses identified five novel risk loci, prioritized candidate target genes within associated regions, and highlighted distinct pathways. Furthermore, the results suggest the presence of distinct regulatory effects of nsCL/P risk variants throughout mid-facial development and shed light on its regulatory architecture. Our integrated data provide a platform to advance hypothesis-driven molecular investigations of nsCL/P and other human facial defects.

Identification of Novel Variants in Cleft Palate-Associated Genes in Brazilian Patients With Non-syndromic Cleft Palate Only

The identification of genetic risk factors for non-syndromic oral clefts is of great importance for better understanding the biological processes related to this heterogeneous and complex group of diseases. Herein we applied whole-exome sequencing to identify potential variants related to non-syndromic cleft palate only (NSCPO) in the multiethnic Brazilian population. Thirty NSCPO samples and 30 sex- and genetic ancestry-matched healthy controls were pooled (3 pools with 10 samples for each group) and subjected to whole-exome sequencing. After filtering, the functional affects, individually and through interactions, of the selected variants and genes were assessed by bioinformatic analyses. As a group, 399 variants in 216 genes related to palatogenesis/cleft palate, corresponding to 6.43%, were exclusively identified in the NSCPO pools. Among those genes are 99 associated with syndromes displaying cleft palate in their clinical spectrum and 92 previously related to cleft lip palate. The most significantly biological processes and pathways overrepresented in the NSCPO-identified genes were associated with the folic acid metabolism, highlighting the interaction between LDL receptor-related protein 6 (LRP6) and 5-methyltetrahydrofolate-homocysteine methyltransferase (MTR) that interconnect two large networks. This study yields novel data on characterization of specific variants and complex processes and pathways related to NSCPO, including many variants in genes of the folate/homocysteine pathway, and confirms that variants in genes related to syndromic cleft palate and cleft lip-palate may cause NSCPO.

DNA Methylation Variation Is Identified in Monozygotic Twins Discordant for Non-syndromic Cleft Lip and Palate

Non-syndromic cleft lip with or without cleft palate (NSCLP) is the most common craniofacial birth defect. The etiology of NSCLP is complex with multiple genes and environmental factors playing causal roles. Although studies have identified numerous genetic markers associated with NSCLP, the role of epigenetic variation remains relatively unexplored. Because of their identical DNA sequences, monozygotic (MZ) twins discordant for NSCLP are an ideal model for examining the potential contribution of DNA methylation to non-syndromic orofacial clefting. In this study, we compared the patterns of whole genome DNA methylation in six MZ twin pairs discordant for NSCLP. Differentially methylated positions (DMPs) and regions (DMRs) were identified in NSCLP candidate genes, including differential methylation in MAFB and ZEB2 in two independent MZ twin pairs. In addition to DNA methylation differences in NSCLP candidate genes, we found common differential methylation in genes belonging to the Hippo signaling pathway, implicating this mechanosensory pathway in the etiology of NSCLP. The results of this novel approach using MZ twins discordant for NSCLP suggests that differential methylation is one mechanism contributing to NSCLP, meriting future studies on the role of DNA methylation in familial and sporadic NSCLP.

Symmetry and fluctuation of cell movements in neural crest-derived facial mesenchyme

In the face, symmetry is established when bilateral streams of neural crest cells leave the neural tube at the same time, follow identical migration routes and then give rise to the facial prominences. However, developmental instability exists, particularly surrounding the steps of lip fusion. The causes of instability are unknown but inability to cope with developmental fluctuations are a likely cause of congenital malformations, such as non-syndromic orofacial clefts. Here, we tracked cell movements over time in the frontonasal mass, which forms the facial midline and participates in lip fusion, using live-cell imaging of chick embryos. Our mathematical examination of cell velocity vectors uncovered temporal fluctuations in several parameters, including order/disorder, symmetry/asymmetry and divergence/convergence. We found that treatment with a Rho GTPase inhibitor completely disrupted the temporal fluctuations in all measures and blocked morphogenesis. Thus, we discovered that genetic control of symmetry extends to mesenchymal cell movements and that these movements are of the type that could be perturbed in asymmetrical malformations, such as non-syndromic cleft lip.

This article has an associated ‘The people behind the papers’ interview.

Highlighted Article: Live imaging of the chick embryo face followed by mathematical analysis of mesenchymal cell tracks captures novel fluctuations between states of order/disorder as well as symmetry/asymmetry, revealing developmental instabilities that are part of normal morphogenesis.

The PAX1 locus at 20p11 is a potential genetic modifier for bilateral cleft lip

Nonsyndromic orofacial clefts (OFCs) are a common birth defect and are phenotypically heterogenous in the structure affected by the cleft - cleft lip (CL) and cleft lip and palate (CLP) - as well as other features, such as the severity of the cleft. Here, we focus on bilateral and unilateral clefts as one dimension of OFC severity, because the genetic architecture of these subtypes is not well understood. We tested for subtype-specific genetic associations in 44 bilateral CL (BCL) cases, 434 unilateral CL (UCL) cases, 530 bilateral CLP cases (BCLP), 1123 unilateral CLP (UCLP) cases, and unrelated controls (N = 1626), using a mixed-model approach. While no novel loci were found, the genetic architecture of UCL was distinct compared to BCL, with 44.03% of suggestive loci having different effects between the two subtypes. To further understand the subtype-specific genetic risk factors, we performed a genome-wide scan for modifiers and found a significant modifier locus on 20p11 (p=7.53×10-9), 300kb downstream of PAX1, that associated with higher odds of BCL vs. UCL, and replicated in an independent cohort (p=0.0018) with no effect in BCLP (p>0.05). We further found that this locus was associated with normal human nasal shape. Taken together, these results suggest bilateral and unilateral clefts may have different genetic architectures. Moreover, our results suggest BCL, the rarest form of OFC, may be genetically distinct from the other OFC subtypes. This expands our understanding of modifiers for OFC subtypes and further elucidates the genetic mechanisms behind the phenotypic heterogeneity in OFCs.

Clinical characterization of 266 patients and family members with cleft lip and/or palate with associated malformations and syndromes

Objectives

To clinically characterize patients and family members with cleft lip and/or palate (CL/P) and associated congenital malformations or syndromes and propose possible inheritance patterns.

Materials and methods

An observational study of patients with CL/P, including medical and family history and intra- and extra-oral examination of their family members, was performed.

Results

Two hundred sixty-six patients, 1257 family members, and 42 pedigrees were included in the study. The distribution of patients according to the cleft type was 57.9% with CLP, 25.2% with cleft palate (CPO), and 12.8% with cleft lip with/without alveolus (CL/A). Seventy-four (27.8%) patients had associated malformations, and 24 (9.2%) a syndrome. The skeletal (27.7%), cardiovascular (19.3%) systems, and eyes (22.9%) were most commonly affected. Pierre Robin Sequence (7 patients) and van der Woude (4) were the most common syndromes. The majority of patients with CPO (19/24) had an associate syndrome. The families had an average of 2.45 affected members.

Conclusion

Individual and interfamilial phenotypic variability in patients with CL/P makes the understanding of etiopathogenesis challenging.

Clinical relevance

The overall prevalence of individuals with CL/P and their pedigrees with associated malformations and syndromes emphasize the need for early identification, interdisciplinary, and long-term planning.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00784-021-03863-2.

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.

Non-syndromic Cleft Palate: An Overview on Human Genetic and Environmental Risk Factors

The epithelial and mesenchymal cells involved in early embryonic facial development are guided by complex regulatory mechanisms. Any factor perturbing the growth, approach and fusion of the frontonasal and maxillary processes could result in orofacial clefts that represent the most common craniofacial malformations in humans. The rarest and, probably for this reason, the least studied form of cleft involves only the secondary palate, which is posterior to the incisive foramen. The etiology of cleft palate only is multifactorial and involves both genetic and environmental risk factors. The intention of this review is to give the reader an overview of the efforts made by researchers to shed light on the underlying causes of this birth defect. Most of the scientific papers suggesting potential environmental and genetic causes of non-syndromic cleft palate are summarized in this review, including genome-wide association and gene–environment interaction studies.

Genome-wide Analyses Identify a Novel Risk Locus for Nonsyndromic Cleft Palate

The 3 major subphenotypes observed in patients with nonsyndromic orofacial clefts (NSOFCs) are nonsyndromic cleft lip only (NSCLO), nonsyndromic cleft lip with palate (NSCLP), and nonsyndromic cleft palate only (NSCPO). However, the genetic architecture underlying NSCPO is largely unknown. Here we performed a 2-stage genome-wide association study (GWAS) on NSCPO and replication analyses of selected variants in other NSOFCs from the Chinese Han population. We identified a novel locus (15q24.3) and a known locus (1q32.2) where variants in or near the gene reached genome-wide significance (2.80 × 10^-13 < P < 1.72 × 10^-08) in a test for association with NSCPO in a case-control design. Although a variant from 15q24.3 was found to be significantly associated with both NSCPO and NSCLP, the direction of estimated effects on risk were opposite. Our functional annotation of the risk alleles within 15q24.3 coupled with previously established roles of the candidate genes within identified risk loci in periderm development, embryonic patterning, and/or regulation of cellular processes supports their involvement in palate development and the pathogenesis of cleft palate. Our study advances the understanding of the genetic basis of NSOFCs and provides novel insights into the pathogenesis of NSCPO.

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.

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.

Non-Syndromic Cleft Lip with or without Cleft Palate: Genome-Wide Association Study in Europeans Identifies a Suggestive Risk Locus at 16p12.1 and Supports SH3PXD2A as a Clefting Susceptibility Gene

Non-syndromic cleft lip with or without cleft palate (nsCL/P) ranks among the most common human congenital malformations, and has a multifactorial background in which both exogenous and genetic risk factors act in concert. The present report describes a genome-wide association study (GWAS) involving a total of 285 nsCL/P patients and 1212 controls from the Netherlands and Belgium. Twenty of the 40 previously reported nsC/LP susceptibility loci were replicated, which underlined the validity of this sample. SNV-based analysis of the data identified an as yet unreported suggestive locus at chromosome 16p12.1 (p-value of the lead SNV: 4.17 × 10^-7). This association was replicated in two of three patient/control replication series (Central European and Yemeni). Gene analysis of the GWAS data prioritized SH3PXD2A at chromosome 10q24.33 as a candidate gene for nsCL/P. To date, support for this gene as a cleft gene has been restricted to data from zebrafish and a knockout mouse model. The present GWAS was the first to implicate SH3PXD2A in non-syndromic cleft formation in humans. In summary, although performed in a relatively small sample, the present GWAS generated novel insights into nsCL/P etiology.

The functional variant of NTN1 contributes to the risk of nonsyndromic cleft lip with or without cleft palate

Previous genome-wide association study of nonsyndromic cleft lip with or without cleft palate (NSCL/P) identified a susceptible variant (rs4791774). We hypothesized that the functional single nucleotide polymorphism (SNP) may be in linkage disequilibrium with this lead SNP. The potential functional SNP (rs4791331) was identified by bioinformatic analysis. A case–control study with 891 orofacial cleft cases and 830 controls was designed to investigate its association with orofacial cleft. The allele-specific DNA-protein binding preference was predicted by JASPAR database. Cell proliferation, cycle and apoptosis, luciferase activity and netrin-1 (NTN1) expression were examined after transfection with the rs4791331 C/T vector in HEK-293 and HEPM cell lines. Forty-six lip tissues of NSCL/P patients were collected to detect NTN1 expression. ntn1a knockout zebrafish models were generated by CRISPR/Cas9 and observed with micro-CT. In the case–control study, the rs4791331-T allele was associated with an increased risk of nonsyndromic orofacial cleft (OR = 1.41, 95% CI = 1.19–1.68), as well as the subgroups cleft lip only (OR = 1.46, 95% CI = 1.14–1.87) and cleft lip and palate (OR = 1.58, 95% CI = 1.27–1.96). The T allele of rs4791331 exhibited anti-apoptotic effects and promoted cell cycle progression at the G1/S transition. Decreased enhancer activity and reduced NTN1 expression following transfection of the T allele were observed. Carriers of the CT/TT genotypes showed significantly lower expression of NTN1 than CC carriers. The ntn1a^−/− zebrafish showed relatively wider intermaxillary fissures. These results indicate that rs4791331 (C > T) disrupted motif binding and led to abnormal expression of NTN1, which may be involved in the development of NSCL/P.

Novel genetic loci affecting facial shape variation in humans

The human face represents a combined set of highly heritable phenotypes, but knowledge on its genetic architecture remains limited, despite the relevance for various fields. A series of genome-wide association studies on 78 facial shape phenotypes quantified from 3-dimensional facial images of 10,115 Europeans identified 24 genetic loci reaching study-wide suggestive association (p < 5 × 10-8), among which 17 were previously unreported. A follow-up multi-ethnic study in additional 7917 individuals confirmed 10 loci including six unreported ones (padjusted < 2.1 × 10-3). A global map of derived polygenic face scores assembled facial features in major continental groups consistent with anthropological knowledge. Analyses of epigenomic datasets from cranial neural crest cells revealed abundant cis-regulatory activities at the face-associated genetic loci. Luciferase reporter assays in neural crest progenitor cells highlighted enhancer activities of several face-associated DNA variants. These results substantially advance our understanding of the genetic basis underlying human facial variation and provide candidates for future in-vivo functional studies.

Genetic factors define CPO and CLO subtypes of nonsyndromicorofacial cleft

Nonsyndromic orofacial cleft (NSOFC) is a severe birth defect that occurs early in embryonic development and includes the subtypes cleft palate only (CPO), cleft lip only (CLO) and cleft lip with cleft palate (CLP). Given a lack of specific genetic factor analysis for CPO and CLO, the present study aimed to dissect the landscape of genetic factors underlying the pathogenesis of these two subtypes using 6,986 cases and 10,165 controls. By combining a genome-wide association study (GWAS) for specific subtypes of CPO and CLO, as well as functional gene network and ontology pathway analysis, we identified 18 genes/loci that surpassed genome-wide significance (P < 5 × 10⁻⁸) responsible for NSOFC, including nine for CPO, seven for CLO, two for both conditions and four that contribute to the CLP subtype. Among these 18 genes/loci, 14 are novel and identified in this study and 12 contain developmental transcription factors (TFs), suggesting that TFs are the key factors for the pathogenesis of NSOFC subtypes. Interestingly, we observed an opposite effect of the genetic variants in the IRF6 gene for CPO and CLO. Moreover, the gene expression dosage effect of IRF6 with two different alleles at the same single-nucleotide polymorphism (SNP) plays important roles in driving CPO or CLO. In addition, PAX9 is a key TF for CPO. Our findings define subtypes of NSOFC using genetic factors and their functional ontologies and provide a clue to improve their diagnosis and treatment in the future.

Although GWAS have discovered 43 genes/loci associated with NSOFC, most previous studies used mixed samples of CL/P subtypes rather than CPO or CLO separately. Our findings define the CPO and CLO subtypes using genetic factors and their functional ontologies based on CPO and CLO GWAS data. In this study, we identified 18 genes/loci that contribute to CPO, CLO or CLP. Fourteen of them are novel and identified, and 12 contain developmental transcription factors (TFs), suggesting that TFs are the key factors for the pathogenesis of NSOFC subtypes. We observed an opposite effect in the strongest associated locus IRF6 for CPO and CLO; this information was omitted by previous CL/P GWAS. Furthermore, we reveal that the gene expression dosage of IRF6 plays important roles in driving CPO or CLO. In addition, we found that PAX9 is a strong genetic factor for CPO. These results suggest that transcription factors are the key genetic reason for the pathogenesis of NSOFC subtypes.

High-Resolution Epigenomic Atlas of Human Embryonic Craniofacial Development

Defects in patterning during human embryonic development frequently result in craniofacial abnormalities. The gene regulatory programs that build the craniofacial complex are likely controlled by information located between genes and within intronic sequences. However, systematic identification of regulatory sequences important for forming the human face has not been performed. Here, we describe comprehensive epigenomic annotations from human embryonic craniofacial tissues and systematic comparisons with multiple tissues and cell types. We identified thousands of tissue-specific craniofacial regulatory sequences and likely causal regions for rare craniofacial abnormalities. We demonstrate significant enrichment of common variants associated with orofacial clefting in enhancers active early in embryonic development, while those associated with normal facial variation are enriched near the end of the embryonic period. These data are provided in easily accessible formats for both craniofacial researchers and clinicians to aid future experimental design and interpretation of noncoding variation in those affected by craniofacial abnormalities.

A systematic genetic analysis and visualization of phenotypic heterogeneity among orofacial cleft GWAS signals

Phenotypic heterogeneity is a hallmark of complex traits, and genetic studies of such traits may focus on them as a single diagnostic entity or by analyzing specific components. For example, in orofacial clefting (OFC), three subtypes-cleft lip (CL), cleft lip and palate (CLP), and cleft palate (CP) have been studied separately and in combination. To further dissect the genetic architecture of OFCs and how a given associated locus may be contributing to distinct subtypes of a trait we developed a framework for quantifying and interpreting evidence of subtype-specific or shared genetic effects in complex traits. We applied this technique to create a "cleft map" of the association of 30 genetic loci with three OFC subtypes. In addition to new associations, we found loci with subtype-specific effects (e.g., GRHL3 [CP], WNT5A [CLP]), as well as loci associated with two or all three subtypes. We cross-referenced these results with mouse craniofacial gene expression datasets, which identified additional promising candidate genes. However, we found no strong correlation between OFC subtypes and expression patterns. In aggregate, the cleft map revealed that neither subtype-specific nor shared genetic effects operate in isolation in OFC architecture. Our approach can be easily applied to any complex trait with distinct phenotypic subgroups.

Three GLI2 mutations combined potentially underlie non-syndromic cleft lip with or without cleft palate in a Chinese pedigree

Background

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is the most common craniofacial birth defect. Its etiology is complex and it has a lifelong influence on affected individuals. Despite many studies, the pathogenic gene alleles are not completely clear. Here, we recruited a Chinese NSCL/P family and explored the candidate causative variants in this pedigree.

Methods

We performed whole‐exome sequencing on two patients and two unaffected subjects of this family. Variants were screened based on bioinformatics analysis to identify the potential etiological alleles. Species conservation analysis, mutation function prediction, and homology protein modeling were also performed to preliminarily evaluate the influence of the mutations.

Results

We identified three rare mutations that are located on a single chromatid (c.2684C > T_p.Ala895Val, c.4350G > T_p.Gln1450His, and c.4622C > A_p.Ser1541Tyr) in GLI2 as candidate causative variants. All of these three mutations were predicted to be deleterious, and they affect amino acids that are conserved in many species. The mutation c.2684C > T was predicted to affect the structure of the GLI2 protein.

Conclusion

Our results further demonstrate that GLI2 variants play a role in the pathogenesis of NSCL/P, and the three rare missense mutations combined are probably the potential disease‐causing variants in this family.

p63 establishes epithelial enhancers at critical craniofacial development genes

The transcription factor p63 is a key mediator of epidermal development. Point mutations in p63 in patients lead to developmental defects, including orofacial clefting. To date, knowledge on how pivotal the role of p63 is in human craniofacial development is limited. Using an inducible transdifferentiation model, combined with epigenomic sequencing and multicohort meta-analysis of genome-wide association studies data, we show that p63 establishes enhancers at craniofacial development genes to modulate their transcription. Disease-specific substitution mutation in the DNA binding domain or sterile alpha motif protein interaction domain of p63, respectively, eliminates or reduces establishment of these enhancers. We show that enhancers established by p63 are highly enriched for single-nucleotide polymorphisms associated with nonsyndromic cleft lip ± cleft palate (CL/P). These orthogonal approaches indicate a strong molecular link between p63 enhancer function and CL/P, illuminating molecular mechanisms underlying this developmental defect and revealing vital regulatory elements and new candidate causative genes.

A comparison of DNA methylation in newborn blood samples from infants with and without orofacial clefts

Background

Isolated orofacial clefts are among the most common congenital birth defects. Although the underlying biological mechanisms remain largely unknown, clefts are thought to be complex disorders influenced by genetic, environmental, and potentially epigenetic factors.

Methods

In blood samples from 2- to 3-day-old infants (n = 747) collected in a nationwide population-based study of orofacial clefts in Norway, we measured DNA methylation profiles for more than 450,000 CpGs and then conducted epigenome-wide association analyses (EWAS). We tested methylation profile difference at each CpG between controls (n = 436) and each of the cleft subtypes (92 cleft lip only, CLO; 84 cleft palate only, CPO; 132 cleft lip and palate, CLP). We also compared controls to various combinations of case groups and compared case subtypes to each other. Finally, using the EWAS results, we searched for larger differentially methylated regions (DMRs) associated with orofacial clefts.

Results

In EWAS comparing controls to individual cleft subtypes, we found no significant associations at a Bonferroni P value threshold of 10⁻⁷. After pooling case groups, we found two significantly differentially methylated CpGs: cg09696939 near gene BICC1 is associated with CLO+CLP (P = 9.58 × 10⁻⁸); cg26985354 in gene CLASRP is associated with CPO+CLP (P = 7.38 × 10⁻⁸). In DMR analysis, we identified a total of 56 significant regions when comparing controls to individual cleft subtypes (10 for CLO, 6 for CPO, 41 for CLP). Only one DMR is shared among the three cleft groups. In combined case group analysis, we found 26 DMRs for CLP+CLO, 31 for CLP+CPO, and 37 when all subtypes are combined. Finally, in case-case comparisons of subtypes, we identified 10 DMRs when comparing CLP to CPO, 9 in CLP compared to CLO, and 13 in CLP compared to CPO.

Conclusions

We identified two individual CpGs and multiple DMRs that differ between controls and cleft case subtypes. Although we find some evidence for the possible role of DNA methylation in etiology of orofacial clefts, our study does not support previous reports of widespread differences in blood DNA methylation between babies with and without facial clefts.

Electronic supplementary material

The online version of this article (10.1186/s13148-019-0638-9) contains supplementary material, which is available to authorized users.

Association of low-frequency genetic variants in regulatory regions with nonsyndromic orofacial clefts

Genome-wide scans have shown that common risk alleles for orofacial clefts (OFC) tend to be located in noncoding regulatory elements and cumulatively explain only part of the heritability of OFCs. Low-frequency variants may account for some of the "missing" heritability. Therefore, we scanned low-frequency variants located within putative craniofacial enhancers to identify novel OFC risk variants and implicate new regulatory elements in OFC pathogenesis. Analyses were performed in a multiethnic sample of 1,995 cases of cleft lip with or without cleft palate (CL/P), 221 cases with cleft palate (CP) only, and 1,576 unaffected controls. One hundred and nineteen putative craniofacial enhancers identified from ChIP-Seq studies in craniofacial tissues or cell lines contained multiple low-frequency (0.01-1%) variants, which we genotyped in participants using a custom Illumina panel. Two complementary statistical approaches, sequence kernel association test and combined multivariate and collapsing, were used to test association of the aggregated low-frequency variants across each enhancer region with CL/P and CP. We discovered a significant association between CP and a branchial arch enhancer near FOXP1 (mm60; p-value = .0002). Additionally, we observed a suggestive association between CL/P and a forebrain enhancer near FOXE1 (hs1717; p-value = .001). These findings suggest that low-frequency variants in craniofacial enhancer regions contribute to the complex etiology of nonsyndromic OFCs.

Associations of genetic variants in endocytic trafficking of epidermal growth factor receptor super pathway with risk of nonsyndromic cleft lip with or without cleft palate

BACKGROUND

The genetic etiology of nonsyndromic cleft lip with or without cleft palate (NSCL/P) has not been fully clarified to date. Epidermal growth factor receptor (EGFR) was reportedly involved in its biological establishment and regulation of cell migration during the embryonic stage.

METHODS

We selected a super pathway of endocytic trafficking of EGFR and investigated the associations of single-nucleotide polymorphisms (SNPs) in the super pathway with the risk of NSCL/P by analyzing our published genome-wide association study (GWAS) data from 504 NSCL/P individuals and 455 controls. After the false discovery rate (FDR) control, we conducted linkage disequilibrium (LD) analyses and conditional regression analyses to obtain independent lead SNPs. We performed LD analyses between the lead SNPs and the reported SNPs to find novel ones from our study. We annotated the lead SNPs and investigated their mapped genes in silico.

RESULTS

A total of 82 SNPs showed a statistical association with the risk of NSCL/P after FDR control. They contained three reported SNPs which were g.117068049G>A (rs7078160), g.117086783C>G (rs10886040), and g.117101266G>T (rs17095681). Four independent lead SNPs were obtained, including g.116979803 T>C (rs1905539) and g.117037960A>G (rs7902502) at 10q25.3, g.35720163G>C (rs75656820) at 17q12, and g.156864512G>A (rs1800877) at 1q23.1. Three of them were in low LD (r²  < 0.5) with the reported SNPs except g.117037960A>G (rs7902502), so these three were newly identified. Lead SNPs were mapped to three genes: SHTN1, AP2B1, and NTRK1. The three genes were relatively more highly expressed in the human craniofacial region and in the proximal maxillary location during the craniofacial development stage of the embryonic mouse.

CONCLUSION

Our results suggested that SHTN1, AP2B1, and NTRK1 might be associated with the development of NSCL/P.

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

Objectives: Orofacial clefting is one of the most prevalent craniofacial malformations. Previous research has demonstrated that unaffected relatives of patients with non-syndromic cleft lip with/without cleft palate (NSCL/P) show distinctive facial features, which can be an expression of underlying NSCL/P susceptibility genes. These results support the hypothesis that genes involved in the occurrence of a cleft also play a role in normal craniofacial development. In this study, we investigated the influence of genetic variants associated with NSCL/P on normal-range variation in facial shape.

Methods: A literature review of genome wide association studies (GWAS) investigating the genetic etiology of NSCL/P was performed, resulting in a list of 75 single nucleotide polymorphisms (SNPs) located in 38 genetic loci. Genotype data were available for 65 of these selected SNPs in three datasets with a combined sample size of 7,418 participants of European ancestry, whose 3D facial images were also available. The effect of each SNP was tested using a multivariate canonical correlation analysis (CCA) against 63 hierarchically-constructed facial segments in each of the three datasets and meta-analyzed. This allowed for the investigation of associations between SNPs known to be involved in NSCL/P and normal-range facial shape variations in a global-to-local perspective, without preselecting specific facial shape features or characteristics.

Results: Six NSCL/P SNPs showed significant associations with variation in normal-range facial morphology. rs6740960 showed significant effects in the chin area (p = 3.71 × 10⁻²⁸). This SNP lies in a non-coding area. Another SNP, rs227731 near the NOG gene, showed a significant effect in the philtrum area (p = 1.96 × 10⁻¹⁶). Three SNPs showed significant effects on the shape of the nose. rs742071 (p = 8.71 × 10⁻¹⁴), rs34246903 (p = 6.87 × 10⁻¹²), and rs10512248 (p = 8.4 × 10⁻⁹). Respectively, these SNPs are annotated to PAX7, MSX1, and PTCH1. Finally, rs7590268, an intron variant of THADA, showed an effect in the shape of the supraorbital ridge (p = 3.84 × 10⁻⁷).

Conclusions: This study provides additional evidence NSCL/P-associated genetic variants influence normal-range craniofacial morphology, with significant effects observed for the chin, the nose, the supraorbital ridges and the philtrum area.

Investigating the shared genetics of non-syndromic cleft lip/palate and facial morphology

There is increasing evidence that genetic risk variants for non-syndromic cleft lip/palate (nsCL/P) are also associated with normal-range variation in facial morphology. However, previous analyses are mostly limited to candidate SNPs and findings have not been consistently replicated. Here, we used polygenic risk scores (PRS) to test for genetic overlap between nsCL/P and seven biologically relevant facial phenotypes. Where evidence was found of genetic overlap, we used bidirectional Mendelian randomization (MR) to test the hypothesis that genetic liability to nsCL/P is causally related to implicated facial phenotypes. Across 5,804 individuals of European ancestry from two studies, we found strong evidence, using PRS, of genetic overlap between nsCL/P and philtrum width; a 1 S.D. increase in nsCL/P PRS was associated with a 0.10 mm decrease in philtrum width (95% C.I. 0.054, 0.146; P = 2x10-5). Follow-up MR analyses supported a causal relationship; genetic variants for nsCL/P homogeneously cause decreased philtrum width. In addition to the primary analysis, we also identified two novel risk loci for philtrum width at 5q22.2 and 7p15.2 in our Genome-wide Association Study (GWAS) of 6,136 individuals. Our results support a liability threshold model of inheritance for nsCL/P, related to abnormalities in development of the philtrum.

Genome-wide mapping of global-to-local genetic effects on human facial shape

Genome-wide association scans of complex multipartite traits like the human face typically use preselected phenotypic measures. Here we report a data-driven approach to phenotyping facial shape at multiple levels of organization, allowing for an open-ended description of facial variation while preserving statistical power. In a sample of 2,329 persons of European ancestry, we identified 38 loci, 15 of which replicated in an independent European sample (n = 1,719). Four loci were completely new. For the others, additional support (n = 9) or pleiotropic effects (n = 2) were found in the literature, but the results reported here were further refined. All 15 replicated loci highlighted distinctive patterns of global-to-local genetic effects on facial shape and showed enrichment for active chromatin elements in human cranial neural crest cells, suggesting an early developmental origin of the facial variation captured. These results have implications for studies of facial genetics and other complex morphological traits.

Novel GREM1 Variations in Sub-Saharan African Patients With Cleft Lip and/or Cleft Palate

OBJECTIVE

Cleft lip and/or cleft palate (CL/P) are congenital anomalies of the face and have multifactorial etiology, with both environmental and genetic risk factors playing crucial roles. Though at least 40 loci have attained genomewide significant association with nonsyndromic CL/P, these loci largely reside in noncoding regions of the human genome, and subsequent resequencing studies of neighboring candidate genes have revealed only a limited number of etiologic coding variants. The present study was conducted to identify etiologic coding variants in GREM1, a locus that has been shown to be largely associated with cleft of both lip and soft palate.

PATIENTS AND METHOD

We resequenced DNA from 397 sub-Saharan Africans with CL/P and 192 controls using Sanger sequencing. Following analyses of the sequence data, we observed 2 novel coding variants in GREM1. These variants were not found in the 192 African controls and have never been previously reported in any public genetic variant database that includes more than 5000 combined African and African American controls or from the CL/P literature.

RESULTS

The novel variants include p.Pro164Ser in an individual with soft palate cleft only and p.Gly61Asp in an individual with bilateral cleft lip and palate. The proband with the p.Gly61Asp GREM1 variant is a van der Woude (VWS) case who also has an etiologic variant in IRF6 gene.

CONCLUSION

Our study demonstrated that there is low number of etiologic coding variants in GREM1, confirming earlier suggestions that variants in regulatory elements may largely account for the association between this locus and CL/P.

FOXE1 Mutation Screening in a Case with Cleft Lip, Hypothyroidism, and Thyroid Carcinoma: A New Syndrome?

A 26-year-old woman is referred to the Internal Medicine consultation due to increases in laboratory studies associated with Papillary Thyroid Carcinoma (PTC) that was confirmed by histopathological studies. Her clinical history revealed that, at 3 months of age, she was successfully treated with surgery for cleft lip (CL) and at the age of 24 years was diagnosed with hypothyroidism. Single nucleotide polymorphisms (SNPs) in FOXE1 and its promoter regions have been associated with various etiologies related to the thyroid, including orofacial clefting, specially cleft palate (CP) and CL, hypothyroidism (HT), and thyroid cancer. The association of CL, HT, and PTC might be component of a new syndrome; however FOXE1 coding region, which has been involved with these entities, has not exhibited mutations or SNPs. Further study of other genes may help in better characterization of the possible syndrome.

A Review of Systemic Corticosteroid Use in Pregnancy and the Risk of Select Pregnancy and Birth Outcomes

The evidence to date regarding corticosteroid exposure in pregnancy and select pregnancy and birth outcomes is limited and inconsistent. The authors provide a narrative review of published literature summarizing the findings for oral clefts, preterm birth, birth weight, preeclampsia, and gestational diabetes mellitus. Whenever possible, the results are limited to oral or systemic administration with a further focus on use in autoimmune disease. Although previous studies of corticosteroid exposure in pregnancy reported an increased risk of oral clefts in the offspring, more recent studies have not replicated these findings.

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.

Distinct DNA methylation profiles in subtypes of orofacial cleft

Background

Epigenetic data could help identify risk factors for orofacial clefts, either by revealing a causal role for epigenetic mechanisms in causing clefts or by capturing information about causal genetic or environmental factors. Given the evidence that different subtypes of orofacial cleft have distinct aetiologies, we explored whether children with different cleft subtypes showed distinct epigenetic profiles.

Methods

In whole-blood samples from 150 children from the Cleft Collective cohort study, we measured DNA methylation at over 450,000 sites on the genome. We then carried out epigenome-wide association studies (EWAS) to test the association between methylation at each site and cleft subtype (cleft lip only (CLO) n = 50; cleft palate only (CPO) n = 50; cleft lip and palate (CLP) n = 50). We also compared methylation in the blood to methylation in the lip or palate tissue using genome-wide data from the same 150 children and conducted an EWAS of CLO compared to CLP in lip tissue.

Results

We found four genomic regions in blood differentially methylated in CLO compared to CLP, 17 in CPO compared to CLP and 294 in CPO compared to CLO. Several regions mapped to genes that have previously been implicated in the development of orofacial clefts (for example, TBX1, COL11A2, HOXA2, PDGFRA), and over 250 associations were novel. Methylation in blood correlated with that in lip/palate at some regions. There were 14 regions differentially methylated in the lip tissue from children with CLO and CLP, with one region (near KIAA0415) showing up in both the blood and lip EWAS.

Conclusions

Our finding of distinct methylation profiles in different orofacial cleft (OFC) subtypes represents a promising first step in exploring the potential role of epigenetic modifications in the aetiology of OFCs and/or as clinically useful biomarkers of OFC subtypes.

Electronic supplementary material

The online version of this article (doi:10.1186/s13148-017-0362-2) contains supplementary material, which is available to authorized users.