Breakthroughs in the genetics of orofacial clefting

Role of ZFHX4 in orofacial clefting based on human genetic data and zebrafish models

Orofacial clefting (OFC) is a frequent congenital anomaly and can occur either in the context of underlying syndromes or in isolation (nonsyndromic). The two common OFC phenotypes are cleft lip with/without cleft palate (CL/P) and cleft palate only (CPO). In this study, we searched for penetrant CL/P genes, by evaluating de novo copy number variants (CNV) from an exome sequencing dataset of 50 nonsyndromic patient-parent trios. We detected a heterozygous 86 kb de novo deletion affecting exons 4-11 of ZFHX4, a gene previously associated with OFC. Genetic and phenotypic data from our in-house and the AGORA cohort (710 and 229 individuals with nonsyndromic CL/P) together with literature and database reviews demonstrate that ZFHX4 variants can lead to both nonsyndromic and syndromic forms not only of CL/P but also CPO. Expression analysis in published single-cell RNA-sequencing data (mouse embryo, zebrafish larva) at relevant time-points support an important role of Zfhx4/zfhx4 in craniofacial development. To characterize the role of zfhx4 in zebrafish craniofacial development, we knocked out/down the zebrafish orthologue. Cartilage staining of the zfhx4 CRISPR F0 knockout and morpholino knockdown at 4 days post-fertilization showed an underdeveloped and abnormally shaped ethmoid plate and cartilaginous jaw (resembling micrognathia). While there is evidence for the dominant inheritance of ZFHX4 variants in OFC, we here present a patient with a possible recessive inheritance. In conclusion, ZFHX4 has a highly heterogeneous phenotypic spectrum and variable mode of inheritance. Our data highlight that ZFHX4 should be considered in genetic testing in patients with nonsyndromic clefting.

Combining genetic and single-cell expression data reveals cell types and novel candidate genes for orofacial clefting

Non-syndromic cleft lip with/without cleft palate (nsCL/P) is one of the most common birth defects and has a multifactorial etiology. To date, over 45 loci harboring common risk variants have been identified. However, the effector genes at these loci, and the cell types that are affected by risk alleles, remain largely unknown. To address this, we combined genetic data from an nsCL/P genome-wide association study (GWAS) with single-cell RNA sequencing data obtained from the heads of unaffected human embryos. Using the recently developed single-cell disease relevance score (scDRS) approach, we identified two major cell types involved in nsCL/P development, namely the epithelium and the HAND2+ pharyngeal arches (PA). Combining scDRS with co-expression networks and differential gene expression analysis, we prioritized nsCL/P candidate genes, some of which were additionally supported by GWAS data (e.g., CTNND1, PRTG, RPL35A, RAB11FIP1, KRT19). Our results suggest that specific epithelial and PA sub-cell types are involved in nsCL/P development, and harbor a substantial fraction of the genetic risk for nsCL/P.

Predictors for long-term relapse of orthodontic treatment in patients with cleft lip and palate. A clinical follow-up study

OBJECTIVES

To identify predictors for long-term relapse of orthodontic therapy in patients with cleft lip and palate (CLP).

MATERIALS AND METHODS

Patients with uni- and bilateral non-syndromal CLP were followed up at least two years after completion of their orthodontic therapy. Plaster casts of the start of treatment (T1), after completion of treatment (T2), and at follow-up (T3) were measured using the modified Huddart Bodenham Index. Characteristics of multidisciplinary therapy were taken from the patient files. Potentially influencing factors of relapse were investigated using logistic regression analyses and Spearman correlations.

RESULTS

In total 58.07% of the included 31 patients showed a stable treatment outcome at follow-up after an average of 6.9 years. Even if relapse occurred, 61.54% of these patients still showed improvement regarding their occlusion compared to baseline. Predictors for the occurrence of relapse were the severity of dysgnathia at baseline (p = 0.039) and the extent of therapeutic change (p = 0.041). The extent of therapeutic change was additionally a predictor for the extent of post-therapeutic relapse (ρ = 0.425; p = 0.019).

CONCLUSIONS

Patients with CLP benefit from their orthodontic therapy in the long term despite an increased tendency to relapse.

CLINICAL RELEVANCE

Results of this long-term study could be used to adapt the treatment concept for patients with CLP and reinforce the significance of a patient-centered orthodontic treatment concept for affected patients.

Longitudinal Three-Dimensional Stereophotogrammetric Growth Analysis in Infants with Unilateral Cleft Lip and Palate from 3 to 12 Months of Age

This longitudinal study aimed to evaluate facial growth and soft tissue changes in infants with complete unilateral cleft lip, alveolus, and palate (CUCLAP) at ages 3, 9, and 12 months. Using 3D images of 22 CUCLAP infants, average faces and distance maps for the entire face and specific regions were created. Color-coded maps highlighted more significant soft tissue changes from 3 to 9 months than from 9 to 12 months. The first interval showed substantial growth in the entire face, particularly in the forehead, eyes, lower lip, chin, and cheeks (p < 0.001), while the second interval exhibited no significant growth. This study provides insights into facial soft tissue growth in CUCLAP infants during critical developmental stages, emphasizing substantial improvements between 3 and 9 months, mainly in the chin, lower lip, and forehead. However, uneven growth occurred in the upper lip, philtrum, and nostrils throughout both intervals, with an overall decline in growth from 9 to 12 months. These findings underscore the dynamic nature of soft tissue growth in CUCLAP patients, highlighting the need to consider these patterns in treatment planning. Future research should explore the underlying factors and develop customized treatment interventions for enhanced facial aesthetics and function in this population.

Multimodal spatiotemporal transcriptomic resolution of embryonic palate osteogenesis

The terminal differentiation of osteoblasts and subsequent formation of bone marks an important phase in palate development that leads to the separation of the oral and nasal cavities. While the morphogenetic events preceding palatal osteogenesis are well explored, major gaps remain in our understanding of the molecular mechanisms driving the formation of this bony union of the fusing palate. Through bulk, single-nucleus, and spatially resolved RNA-sequencing analyses of the developing secondary palate, we identify a shift in transcriptional programming between embryonic days 14.5 and 15.5 pinpointing the onset of osteogenesis. We define spatially restricted expression patterns of key osteogenic marker genes that are differentially expressed between these developmental timepoints. Finally, we identify genes in the palate highly expressed by palate nasal epithelial cells, also enriched within palatal osteogenic mesenchymal cells. This investigation provides a relevant framework to advance palate-specific diagnostic and therapeutic biomarker discovery.

Cleft palate in fetuses: feasibility of early diagnosis by Crystal and Realistic Vue rendering 3D ultrasound technology in the first trimester

Objectives

This study aimed to evaluate the feasibility of direct visualization of a normal fetal palate and detect cleft palate in the first trimester with a novel three-dimensional ultrasound (3D US) technique, Crystal and Realistic Vue (CRV) rendering technology.

Methods

Two-dimensional (2D) images and 3D volumes of healthy and cleft palate fetuses at 11-13⁺⁶ weeks were obtained prospectively. 2D ultrasound views included the coronal view of the retronasal triangle and the midsagittal view of the face. 3D-CRV views were analyzed by multiplanar mode display. The pregnancy outcomes of all fetuses were determined during the follow-up period.

Results

In our study, 124 fetuses were recruited, including 100 healthy fetuses and 24 cleft palate fetuses. The cleft palate with lip was observed in 23 fetuses (bilateral in 15, unilateral in 6, median in 2), and one cleft palate was only found in the abnormal group. The bilateral (n = 12) and median (n = 2) cleft palates with lips and the cleft palate alone (n = 1) were associated with other anatomical or chromosomal abnormalities, and one unilateral cleft palate with cleft lip had concomitant NT thickening. In the cleft palate fetus group, 16 fetuses suffered intrauterine death, which was associated with other structural or chromosomal abnormalities in 14 fetuses, seven cases were terminated after consultation, and one was delivered at term. The coronal view of the retronasal triangle and the midsagittal view was easily obtained in all fetuses. 3D-CRV images of palatal parts were clearly obtained in all cases. Unilateral, bilateral, and median cleft palates with cleft lips were visually demonstrated and classified by the 3D-CRV technique.

Conclusion

It is feasible to identify the palate by 3D-CRV in the first trimester in both healthy and cleft palate fetuses. Together with 2D ultrasonography as a complementary diagnostic tool, 3D-CRV is helpful in classifying the cleft palate with a reasonable degree of certainty.

Integrated spatiotemporal transcriptomic resolution of embryonic palate osteogenesis

The differentiation of osteoblasts and the subsequent formation of bone marks an important terminal phase in palate formation that leads to the separation of the oral and nasal cavities. While the developmental events that precede palatal osteogenesis are well explored, major gaps remain in our understanding of the molecular mechanisms that lead to the bony union of fusing palatal shelves. Herein, the timeline of osteogenic transcriptional programming is unveiled in the embryonic palate by way of integrated bulk, single-cell, and spatially resolved RNA-seq analyses. We define spatially restricted expression patterns of key marker genes, both regulatory and structural, that are differentially expressed during palatal fusion, including the identification of several novel genes ( Deup1, Dynlrb2, Lrrc23 ) spatially restricted in expression to the palate, providing a relevant framework for future studies that identify new candidate genes for cleft palate anomalies in humans as well as the timing of mammalian embryonic palatal osteogenesis.

Identification of Novel Risk Variants of Non-Syndromic Cleft Palate by Targeted Gene Panel Sequencing

Non-syndromic cleft palate (ns-CP) has a genetically heterogeneous aetiology. Numerous studies have suggested a crucial role of rare coding variants in characterizing the unrevealed component of genetic variation in ns-CP called the "missing heritability". Therefore, this study aimed to detect low-frequency variants that are implicated in ns-CP aetiology in the Polish population. For this purpose, coding regions of 423 genes associated with orofacial cleft anomalies and/or involved with facial development were screened in 38 ns-CP patients using the next-generation sequencing technology. After multistage selection and prioritisation, eight novel and four known rare variants that may influence an individual's risk of ns-CP were identified. Among detected alternations, seven were located in novel candidate genes for ns-CP, including COL17A1 (c.2435-1G>A), DLG1 (c.1586G>C, p.Glu562Asp), NHS (c.568G>C, p.Val190Leu-de novo variant), NOTCH2 (c.1997A>G, p.Tyr666Cys), TBX18 (c.647A>T, p.His225Leu), VAX1 (c.400G>A, p.Ala134Thr) and WNT5B (c.716G>T, p.Arg239Leu). The remaining risk variants were identified within genes previously linked to ns-CP, confirming their contribution to this anomaly. This list included ARHGAP29 (c.1706G>A, p.Arg569Gln), FLNB (c.3605A>G, Tyr1202Cys), IRF6 (224A>G, p.Asp75Gly-de novo variant), LRP6 (c.481C>A, p.Pro161Thr) and TP63 (c.353A>T, p.Asn118Ile). In summary, this study provides further insights into the genetic components contributing to ns-CP aetiology and identifies novel susceptibility genes for this craniofacial anomaly.

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.

Analysis of candidate genes for cleft lip ± cleft palate using murine single-cell expression data

Introduction: Cleft lip ± cleft palate (CL/P) is one of the most common birth defects. Although research has identified multiple genetic risk loci for different types of CL/P (i.e., syndromic or non-syndromic forms), determining the respective causal genes and understanding the relevant functional networks remain challenging. The recent introduction of single-cell RNA sequencing (scRNA-seq) has provided novel opportunities to study gene expression patterns at cellular resolution. The aims of our study were to: (i) aggregate available scRNA-seq data from embryonic mice and provide this as a resource for the craniofacial community; and (ii) demonstrate the value of these data in terms of the investigation of the gene expression patterns of CL/P candidate genes. Methods and Results: First, two published scRNA-seq data sets from embryonic mice were re-processed, i.e., data representing the murine time period of craniofacial development: (i) facial data from embryonic day (E) E11.5; and (ii) whole embryo data from E9.5-E13.5 from the Mouse Organogenesis Cell Atlas (MOCA). Marker gene expression analyses demonstrated that at E11.5, the facial data were a high-resolution representation of the MOCA data. Using CL/P candidate gene lists, distinct groups of genes with specific expression patterns were identified. Among others we identified that a co-expression network including Irf6, Grhl3 and Tfap2a in the periderm, while it was limited to Irf6 and Tfap2a in palatal epithelia, cells of the ectodermal surface, and basal cells at the fusion zone. The analyses also demonstrated that additional CL/P candidate genes (e.g., Tpm1, Arid3b, Ctnnd1, and Wnt3) were exclusively expressed in Irf6+ facial epithelial cells (i.e., as opposed to Irf6- epithelial cells). The MOCA data set was finally used to investigate differences in expression profiles for candidate genes underlying different types of CL/P. These analyses showed that syndromic CL/P genes (syCL/P) were expressed in significantly more cell types than non-syndromic CL/P candidate genes (nsCL/P). Discussion: The present study illustrates how scRNA-seq data can empower research on craniofacial development and disease.

Orofacial clefts lead to increased pro-inflammatory cytokine levels on neonatal oral mucosa

Orofacial clefts (OFC) are frequent congenital malformations characterized by insufficient separation of oral and nasal cavities and require presurgical infant orthopedics and surgical interventions within the first year of life. Wound healing disorders and higher prevalence of gingivitis and plaque levels are well-known challenges in treatment of children with OFC. However, oral inflammatory mediators were not investigated after birth using non-invasive sampling methods so far. In order to investigate the impact of OFC on oral cytokine levels, we collected tongue smear samples from 15 neonates with OFC and 17 control neonates at two time points (T), T0 at first consultation after birth, and T1, 4 to 5 weeks later. The samples were analyzed using multiplex immunoassay. Overall, we found significantly increased cytokine levels (TNF, IL-1β/-2/-6/-8/-10) in tongue smear samples from neonates with OFC compared to controls, especially at T0. The increase was even more pronounced in neonates with a higher cleft severity. Further, we detected a significant positive correlation between cleft severity score and distinct pro-inflammatory mediators (GM-CSF, IL-1β, IL-6, IL-8) at T0. Further, we found that breast-milk (bottle) feeding was associated with lower levels of pro-inflammatory cytokines (IL-6/-8) in neonates with OFC compared to formula-fed neonates. Our study demonstrated that neonates with OFC, especially with high cleft severity, are characterized by markedly increased inflammatory mediators in tongue smear samples within the first weeks of life potentially presenting a risk for oral inflammatory diseases. Therefore, an inflammatory monitoring of neonates with (severe) OFC and the encouragement of mother to breast-milk (bottle) feed might be advisable after birth and/or prior to cleft surgery.

Association of gestational hypertension and preeclampsia with nonsyndromic orofacial clefts in China: a large prospective cohort study

BACKGROUND

The associations between hypertensive disorders of pregnancy and nonsyndromic orofacial clefts (NSOFCs) are not consistent or based on case-control study design. We hypothesize that OFCs and NSOFCs are associated with hypertensive disease in pregnancy.

METHODS

Data were collected from the Project for Neural Tube Defects Prevention (1993-1996), a large population-based cohort study conducted in two southern provinces of China. We used a system to record all births after 20 complete gestational weeks, including live births and stillbirths, and all structural congenital anomalies regardless of gestational week. A total of 200 215 singleton live births without other external birth defects were finally included.

RESULTS

The incidence of NSOFCs was 20.2 per 10 000 for the whole population, and 20.5 and 39.2 per 10 000 for women with gestational hypertension and preeclampsia, respectively. Compared with the nonpreeclampsia group, preeclampsia was associated with an increased risk of NSOFCs [adjusted risk ratio (RR) = 2.02, 95% confidence interval (CI): 1.27-3.20], cleft lip with or without cleft palate (CL/P) [adjusted RR = 2.24, 95% CI 1.37-3.65], and cleft lip and palate (CLP) [adjusted RR = 2.60, 95% CI 1.45-4.67] but not cleft lip only (CLO) [adjusted RR = 1.66, 95% CI 0.68-4.07] or cleft palate only (CPO) [adjusted RR = 1.09, 95% CI 0.27-4.45]. No associations were observed between gestational hypertension and any types of NSOFCs.

CONCLUSION

Our study supported that among hypertensive disorders of pregnancy, only preeclampsia increased the risk of NSOFCs and its subtypes (CL/P and CLP).

[Association study between haplotypes of WNT signaling pathway genes and nonsyndromic oral clefts among Chinese Han populations]

OBJECTIVE

To explore whether WNT signaling pathway genes were associated with non-syndromic oral clefts (NSOC) based on haplotypes analyses among 1 008 Chinese NSOC case-parent trios.

METHODS

The genome-wide association study (GWAS) data of 806 Chinese non-syndromic cleft lip with or without cleft palate (NSCL/P) trios and 202 Chinese non-syndromic cleft palate (NSCP) case-parent trios were drawn from the International Consortium to Identify Genes and Interactions Controlling Oral Clefts (ICOCs) study GWAS data set, whose Chinese study population were recruited from four provinces in China, namely Taiwan, Shandong, Hubei, and Sichuan provinces. The process of DNA genotyping was conducted by the Center for Inherited Disease Research in the Johns Hopkins University, using Illumina Human610-Quad v.1_B Bead Chip. The method of sliding windows was used to determine the haplotypes for analyses, including 2 SNPs haplotypes and 3 SNPs haplotypes. Haplotypes with a frequency lower than 1% were excluded for further analyses. To further assess the association between haplotypes and NSOC risks, and the transmission disequilibrium test (TDT) was performed. The Bonferroni method was adopted to correct multiple tests in the study, with which the threshold of statistical significance level was set as P < 0.05 divided by the number of tests, e.g P < 3.47×10^-4 in the current stu-dy. All the statistical analyses were performed by using plink (v1.07).

RESULTS

After quality control, a total of 144 single nucleotide polymorphisms (SNPs) mapped in seven genes in WNT signaling pathway were included for the analyses among the 806 Chinese NSCL/P trios and 202 Chinese NSCP trios. A total of 1 042 haplotypes with frequency higher than 1% were included for NSCL/P analyses and another 1 057 haplotypes with frequency higher than 1% were included for NSCP analyses. Results from the TDT analyses showed that a total of 69 haplotypes were nominally associated with the NSCL/P risk among Chinese (P < 0.05). Another 34 haplotypes showed nominal significant association with the NSCP risk among Chinese (P < 0.05). However, none of these haplotypes reached pre-defined statistical significance level after Bonferroni correction (P>3.47×10^-4).

CONCLUSION

This study failed to observe any statistically significant associations between haplotypes of seven WNT signaling pathway genes and the risk of NSOC among Chinese. Further studies are warranted to replicate the findings here.

Target sequencing reveals the association between variants in VAX1 and NSCL/P in Chinese population

OBJECTIVE

A significant genetic association between rs7078160 in VAX1 and NSCL/P has been established through genome-wide association studies (GWAS), and we previously replicated the association in the Chinese population. The critical issue in the post-GWAS era is to identify functional variations that have a real impact on disease in the susceptible regions highlighted by GWAS. This study aimed to elucidate functional variants in VAX1 fully.

MATERIALS AND METHODS

Firstly, target sequencing was performed on 159 NSCL/P patients, followed by association analysis to discover disease-associated single nucleotide polymorphisms (SNPs); we then replicated the findings using a larger sample (1626 cases, 2255 controls) and investigated how candidate SNPs affect disease occurrence using extensive annotation databases. Additionally, we compared the genetic profiles of NSCL/P subtypes.

RESULTS

In this study, 6 SNPs in VAX1 were identified to be associated with NSCL/P in the Western Han Chinese population. Five of them were predicted to influence transcriptional factor biding ability and were expression quantitative trait loci (eQTLs) of nearby genes in multiple tissues.

CONCLUSION

The previously reported association between rs7078160 and NSCL/P was successfully replicated. Moreover, our findings firstly revealed that 5 SNPs in VAX1 are associated with NSCL/P in the Western Han Chinese population.

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.

Extending the allelic spectrum at noncoding risk loci of orofacial clefting

Genome-wide association studies (GWAS) have generated unprecedented insights into the genetic etiology of orofacial clefting (OFC). The moderate effect sizes of associated noncoding risk variants and limited access to disease-relevant tissue represent considerable challenges for biological interpretation of genetic findings. As rare variants with stronger effect sizes are likely to also contribute to OFC, an alternative approach to delineate pathogenic mechanisms is to identify private mutations and/or an increased burden of rare variants in associated regions. This report describes a framework for targeted resequencing at selected noncoding risk loci contributing to nonsyndromic cleft lip with/without cleft palate (nsCL/P), the most frequent OFC subtype. Based on GWAS data, we selected three risk loci and identified candidate regulatory regions (CRRs) through the integration of credible SNP information, epigenetic data from relevant cells/tissues, and conservation scores. The CRRs (total 57 kb) were resequenced in a multiethnic study population (1061 patients; 1591 controls), using single-molecule molecular inversion probe technology. Combining evidence from in silico variant annotation, pedigree- and burden analyses, we identified 16 likely deleterious rare variants that represent new candidates for functional studies in nsCL/P. Our framework is scalable and represents a promising approach to the investigation of additional congenital malformations with multifactorial etiology.

Study on the role of methylation in nonsyndromic cleft lip with or without cleft palate using a monozygotic twin model

OBJECTIVE

The research on the etiology of nonsyndromic cleft lip with or without cleft palate(NSCL/P) is challenging, and DNA methylation has an impact on the formation of cleft lip and palate.

SUBJECTS

In this study, one of a pair of monozygotic twins (T1) had nonsyndromic cleft lip (NSCL), and one of a pair of monozygotic twins (T2) had nonsyndromic cleft lip and palate (NSCLP). We determined the methylation profiles of more than 850,000 CpGs in the DNA of the blood samples from the two pairs of monozygotic twins.

RESULT

Methylation data indicated that 1184 differentially methylated CpG sites were found in the T1 group (651 hypermethylated and 533 hypomethylated) and 8099 differentially methylated CpG sites in the T2 group (1713 hypermethylated and 6386 hypomethylated) compared with the healthy twin.The common difference was 107 methylation sites.GO enrichment analysis showed that regulation of smooth muscle cell migration and actin cytoskeleton reorganization were the most prominent classes.KEGG pathway enrichment analysis showed that the TGF-β signaling pathway, Notch signaling pathway and Wnt signaling pathway are relevant to the formation of NSCL/P.Two selected genes (NTN1 and PLEKHA7) are involved in the formation of NSCL/P.

CONCLUSION

These findings provide some support for the hypothesis that abnormal DNA methylation may influence the formation of clefts.

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.

[Evaluating the effect of WNT pathway genes considering interactions on the risk of non-syndromic oral clefts among Chinese populations]

OBJECTIVE

In this study, we used genome-wide association study (GWAS) data to explore whether WNT pathway genes were associated with non-syndromic oral clefts (NSOC) considering gene-gene interaction and gene-environment interaction.

METHODS

We conducted the analysis using 806 non-syndromic cleft lip with or without cleft palate (NSCL/P) case-parent trios and 202 non-syndromic cleft palate (NSCP) case-parent trios among Chinese populations selected from an international consortium established for a GWAS of non-syndromic oral clefts. Genotype data and maternal environmental exposures were collected through DNA samples and questionnaires. Conditional Logistic regression models were adopted to explore gene-gene interaction and gene-environment in teraction using trio package in R software. The threshold of significance level was set as 3.47×10^-4 using Bonferroni correction.

RESULTS

A total of 144 single nucleotide polymorphisms (SNPs) in seven genes passed the quality control process in NSCL/P trios and NSCP trios, respectively. Totally six pairs of SNPs interactions showed statistically significant SNP-SNP interaction (P < 3.47×10^-4) after Bonferroni correction, which were rs7618735 (WNT5A) and rs10848543 (WNT5B), rs631948 (WNT11) and rs556874 (WNT5A), and rs631948 (WNT11) and rs472631 (WNT5A) among NSCL/P trios; rs589149 (WNT11) and rs4765834 (WNT5B), rs1402704 (WNT11) and rs358792 (WNT5A), and rs1402704 (WNT11) and rs358793 (WNT5A) among NSCP trios, respectively. In addition, no significant result was found for gene-environment interaction analysis in both of the NSCL/P trios and NSCP trios.

CONCLUSION

Though this study failed to detect significant association based on gene-environment interactions of seven WNT pathway genes and the risk of NSOC, WNT pathway genes may influence the risk of NSOC through potential gene-gene interaction.

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.

A comprehensive analysis of AHRR gene as a candidate for cleft lip with or without cleft palate

Cleft lip and palate (CL/P) is among the most common congenital malformations and affects 1 in 700 newborns. CL/P is caused by genetic and environmental factors (maternal smoking, alcohol or drug use and others). Many genes and loci were associated with cleft lip/palate but the amount of heterogeneity justifies identifying new causal genes and variants. AHRR (Aryl-Hydrocarbon Receptor Repressor) gene has recently been related to CL/P however, few functional studies analyze the genotypephenotype interaction of AHRR with CL/P. Several studies associate the molecular pathway of AHRR to CL/P which indicates this gene as a functional candidate in CL/P etiology.

METHODS

Systematic Literature Review was performed using PUBMED database with the keywords cleft lip, cleft palate, orofacial cleft, AHRR and synonyms. SLR resulted in 37 included articles.

RESULTS

AHRR is a positional and functional candidate gene for CL/P. In silico analysis detected interactions with other genes previously associated to CL/P like ARNT and CYP1A1. AHRR protein regulates cellular toxicity through TCDD mediated AHR pathway. Exposure to TCDD in animal embryos is AHR mediated and lead to cleft palate due to palate fusion failure and post fusion rupture. AHRR regulates cellular growth and differentiation, fundamental to lip and palatogenesis. AHRR decreases carcinogenesis and recently a higher tumor risk has been described in CL/P patients and families. AHRR is also a smoking biomarker due to changed methylation sites found in smokers DNA although folate intake may partially revert these methylation alterations. This corroborates the role of maternal smoking and lack of folate supplementation as risk factors for CL/P.

CONCLUSION

This research identified the importance of AHRR in dioxin response and demonstrated an example of genetic and environmental interaction, indispensable in the development of many complex diseases.

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.

First trimester sonographic diagnosis of orofacial defects. Review of literature

Ultrasound has been used since the 1950s as a useful tool for the screening of several pregnancy abnormalities. The National Institute for Excellence in Health and Care (NICE) guidelines for prenatal control recommend its routine use between 12 and 20 weeks of gestational age, given that during the first trimester, a series of very frequent markers that determine a high risk of fetal anomalies can be evaluated. Among these markers, the most frequently studied are: increased nuchal translucency, the absence of nasal bones, increased tricuspid regurgitation, and altered flow in the venous duct. There is also a new ultrasound technique consisting of the evaluation of the retronasal triangle view, which captures the coronal plane of the face in which the primary palate and the frontal process of the maxillary are simultaneously visualized, presenting high sensitivity and specificity for malformations such as oropalatine clefts, malformations of the nasal bones, and micrognathia. The purpose of this article is to make a comprehensive review of first trimester sonographic diagnosis of orofacial defects.

[Cleft lip and palate]

BACKGROUND

Cleft lip and palate (CLP) represents a group of malformations of unknown etiology but similar phenotypes. This implies consequences for the diagnostics, therapy, prevention, prognosis and risk estimation.

OBJECTIVE

Definition of CLP subtypes and the embryonic development, clarification of correlations and differences between entities using epidemiological data, overview of the present state of genetic analyses, correlation to syndromes, sequences and associations and resulting consequences for clinical practice.

MATERIAL AND METHODS

Update on embryological development of the face, summary of epidemiological and genetic studies and considerations on pedopathological and forensic aspects.

RESULTS

Syndromic and non-syndromic CLP exhibit different and highly variable etiologies, therapeutic needs and prognosis. A thorough understanding is mandatory to distinguish between the different subgroups. In addition to specific aspects of CLP for the pediatric (forensic) pathologist this article provides an overall view of the topic which aims to help understand these malformations.

Associations Between MTR A2756G, MTRR A66G, and TCN2 C776G Polymorphisms and Risk of Nonsyndromic Cleft Lip With or Without Cleft Palate: A Meta-Analysis

OBJECTIVE

We conducted a meta-analysis to investigate the associations of methionine synthase (MTR) A2756G, methionine synthase reductase (MTRR) A66G, and transcobalamin 2 (TCN2) C776G gene polymorphisms with nonsyndromic cleft lip with or without cleft palate (NSCL/P).

MATERIALS AND METHODS

The PubMed, Web of Science, Embase, and Wiley Online Library databases and the China Biomedical Literature Service System (SinoMed) were searched for relevant articles to explore the associations between the MTR A2756G, MTRR A66G, and TCN2 C776G polymorphisms and the risk of NSCL/P. We performed overall comparisons and stratified analyses according to the ethnicity, type of NSCL/P, and Hardy-Weinberg equilibrium (HWE) of the control group. Pooled odds ratios (ORs) and 95% confidence intervals (CIs) were applied to estimate the associations of these gene polymorphisms with NSCL/P risk using fixed-effects or random-effects models incorporating five genetic models.

RESULTS

Ultimately, 12 articles were included in this study. The pooled results did not reveal a significant association of the MTR A2756G polymorphism with NSCL/P risk (G vs. A: OR = 0.95, 95% CI = 0.82-1.11, p = 0.55). Similar results were observed for the MTRR A66G polymorphism (G vs. A: OR = 0.99, 95% CI = 0.82-1.18, p = 0.72) and the TCN2 C776G polymorphism (G vs. C: OR = 0.95, 95% CI = 0.86-1.06, p = 0.37).

CONCLUSION

In summary, the MTR A2756G, MTRR A66G, and TCN2 C776G polymorphisms might not be associated with NSCL/P risk.

IRF6 rs2235375 single nucleotide polymorphism is associated with isolated non-syndromic cleft palate but not with cleft lip with or without palate in South Indian population

Introduction

Transcription factors are very diverse family of proteins involved in activating or repressing the transcription of a gene at a given time. Several studies using animal models demonstrated the role of transcription factor genes in craniofacial development.

Objective

We aimed to investigate the association of IRF6 intron-6 polymorphism in the non-syndromic cleft lip with or without palate in a South Indian population.

Methods

173 unrelated nonsyndromic cleft lip with or without cleft palate patients and 176 controls without clefts patients were genotyped for IRF6 rs2235375 variant by allele-specific amplification using the KASPar single nucleotide polymorphism genotyping system. The association between interferon regulatory factor-6 gene intron-6 dbSNP208032210:g.G>C (rs2235375) single nucleotide polymorphism and non-syndromic cleft lip with or without palate risk was investigated by chi-square test.

Results

There were significant differences in genotype or allele frequencies of rs2235375 single nucleotide polymorphism between controls and cases with non-syndromic cleft lip with or without palate. IRF6 rs2235375 variant was significantly associated with increased risk of non-syndromic cleft lip with or without palate in co-dominant, dominant (OR: 1.19; 95% CI 1.03–2.51; p = 0.034) and allelic models (OR: 1.40; 95% CI 1.04–1.90; p = 0.028). When subset analysis was applied significantly increased risk was observed in cleft palate only group (OR dominant: 4.33; 95% CI 1.44–12.97; p = 0.005).

Conclusion

These results suggest that IRF6 rs2235375 SNP play a major role in the pathogenesis and risk of developing non-syndromic cleft lip with or without palate.

Association between a common missense variant in LOXL3 gene and the risk of non-syndromic cleft palate

To investigate possible association between functional common variants in the lysyl oxidase like 3 gene and non-syndromic cleft palate we selected a common missense variant p.Ile615Phe (rs17010021), which was predicted to have a probably damaging effect on the lysyl oxidase like 3 enzyme. We genotyped 258 non-syndromic cleft palate case-parent triads of European origin and tested genetic association using the transmission disequilibrium test and log-linear regression analyses of genotypic relative risks and of parent-of-origin effects. The observed genotype frequency in parents was in Hardy-Weinberg equilibrium. Compared with wild-type Ile/Ile homozygotes, the relative risks for Phe/Phe homozygote infants was 6.87 (P value 3.0 × 10^-3 ), while that for Ile/Phe heterozygotes was not significant. Assuming an autosomal recessive model, the relative risks for Phe/Phe genotype resulted 10.54 (P value 2.9 × 10^-5 ), with a 3.6% population attributable risk. No parent-of-origin effect was observed. The identification in lysyl oxidase like 3 of a missense variant which under a recessive model associates with 10-fold increased risk of non-syndromic cleft palate supports the hypothesis that the genetic etiology of this congenital anomaly includes relatively uncommon recessive variants with moderate penetrance and located in genes which are also involved in syndromes that include cleft palate as part of the phenotype. Our findings require functional validation and replication in a larger independent genetic association study.

Association Between DLX4 Polymorphisms and Nonsyndromic Orofacial Clefts in a Chinese Han Population

OBJECTIVE

Distal-less 4 ( DLX4) was recently identified as the causative gene for a syndromic form of cleft lip with or without cleft palate, and further biological analyses have established the importance of Dlx4 gene in craniofacial development, which suggested DLX4 as a promising candidate to further investigate any possible association between DLX4 polymorphisms and risk to nonsyndromic orofacial clefts (NSOFCs).

DESIGN

Single-nucleotide polymorphisms (SNPs) with minor allele frequency >5% in the Han Chinese population which locate in the 5' flanking region, 5'/3'-untranslated region, or coding region with nonsynonymous changes in DLX4 were selected. Four SNPs (rs58769681, rs1058562, rs1058564, and rs8066341) were thus included in the following genotyping using the TaqMan 5'-exonuclease allelic discrimination assay in a case-control cohort with 1522 individuals.

RESULTS

None of SNPs were associated with NSOFCat the allele and genotype levels in general and stratified single-marker analysis, including genotypic distributions under different modes of inheritance. In linkage disequilibrium (LD) analysis, we found strong LD ( r² > 0.8) between any 2 of the SNPs, respectively. Further haplotyping identified haplotypes C-C (formed by rs1058564 and rs1058562) and C-C-A (formed by rs1058564, rs1058562, and rs58769681) which reached the significance threshold ( P < .05); nevertheless, none of them survived the multiple comparison correction.

CONCLUSIONS

Our findings indicated the hypothesis that DLX4 variants contributing to NSOFC risk should be interpreted with caution. Further replications in diverse ethnic origins and larger cohorts are still warranted.

Nonsyndromic cleft palate: An association study at GWAS candidate loci in a multiethnic sample

BACKGROUND

Nonsyndromic cleft palate only (nsCPO) is a common and multifactorial form of orofacial clefting. In contrast to successes achieved for the other common form of orofacial clefting, that is, nonsyndromic cleft lip with/without cleft palate (nsCL/P), genome wide association studies (GWAS) of nsCPO have identified only one genome wide significant locus. Aim of the present study was to investigate whether common variants contribute to nsCPO and, if so, to identify novel risk loci.

METHODS

We genotyped 33 SNPs at 27 candidate loci from 2 previously published nsCPO GWAS in an independent multiethnic sample. It included: (i) a family-based sample of European ancestry (n = 212); and (ii) two case/control samples of Central European (n = 94/339) and Arabian ancestry (n = 38/231), respectively. A separate association analysis was performed for each genotyped dataset, and meta-analyses were performed.

RESULTS

After association analysis and meta-analyses, none of the 33 SNPs showed genome-wide significance. Two variants showed nominally significant association in the imputed GWAS dataset and exhibited a further decrease in p-value in a European and an overall meta-analysis including imputed GWAS data, respectively (rs395572: P_MetaEU  = 3.16 × 10^-4 ; rs6809420: P_MetaAll  = 2.80 × 10^-4 ).

CONCLUSION

Our findings suggest that there is a limited contribution of common variants to nsCPO. However, the individual effect sizes might be too small for detection of further associations in the present sample sizes. Rare variants may play a more substantial role in nsCPO than in nsCL/P, for which GWAS of smaller sample sizes have identified genome-wide significant loci. Whole-exome/genome sequencing studies of nsCPO are now warranted.

A Genome-Wide Search for Gene-Environment Effects in Isolated Cleft Lip with or without Cleft Palate Triads Points to an Interaction between Maternal Periconceptional Vitamin Use and Variants in ESRRG

Background: It is widely accepted that cleft lip with or without cleft palate (CL/P) results from the complex interplay between multiple genetic and environmental factors. However, a robust investigation of these gene-environment (GxE) interactions at a genome-wide level is still lacking for isolated CL/P. Materials and Methods: We used our R-package Haplin to perform a genome-wide search for GxE effects in isolated CL/P. From a previously published GWAS, genotypes and information on maternal periconceptional cigarette smoking, alcohol intake, and vitamin use were available on 1908 isolated CL/P triads of predominantly European or Asian ancestry. A GxE effect is present if the relative risk estimates for gene-effects in the offspring are different across exposure strata. We tested this using the relative risk ratio (RRR). Besides analyzing all ethnicities combined ("pooled analysis"), separate analyses were conducted on Europeans and Asians to investigate ethnicity-specific effects. To control for multiple testing, q-values were calculated from the p-values. Results: We identified significant GxVitamin interactions with three SNPs in "Estrogen-related receptor gamma" (ESRRG) in the pooled analysis. The RRRs (95% confidence intervals) were 0.56 (0.45-0.69) with rs1339221 (q = 0.011), 0.57 (0.46-0.70) with rs11117745 (q = 0.011), and 0.62 (0.50-0.76) with rs2099557 (q = 0.037). The associations were stronger when these SNPs were analyzed as haplotypes composed of two-SNP and three-SNP combinations. The strongest effect was with the "t-t-t" haplotype of the rs1339221-rs11117745-rs2099557 combination [RRR = 0.50 (0.40-0.64)], suggesting that the effects observed with the other SNP combinations, including those in the single-SNP analyses, were mainly driven by this haplotype. Although there were potential GxVitamin effects with rs17734557 and rs1316471 and GxAlcohol effects with rs9653456 and rs921876 in the European sample, respectively, none of the SNPs was located in or near genes with strong links to orofacial clefts. GxAlcohol and GxSmoke effects were not assessed in the Asian sample because of a lack of observations for these exposures. Discussion/Conclusion: We identified significant interactions between vitamin use and variants in ESRRG in the pooled analysis. These GxE effects are novel and warrant further investigations to elucidate their roles in orofacial clefting. If validated, they could provide prospects for exploring the impact of estrogens and vitamins on clefting, with potential translational applications.

Evaluating LINE-1 methylation in cleft lip tissues and its association with early pregnancy exposures

Aim: To pilot investigation of methylation of long interspersed nucleotide element-1 in lip tissues from infants with nonsyndromic cleft lip, and its association with maternal periconceptional exposures. Methods: The lateral and medial sides of the cleft lips of 23 affected infants were analyzed for long interspersed nucleotide element-1 methylation by bisulfite conversion and pyrosequencing. Results: The medial side showed 1.8% higher methylation compared with the lateral side; p = 0.031, particularly in male infants (2.7% difference; p = 0.011) or when the mothers did not take folic acid during periconceptional period (2.4% difference; p = 0.011). These results were not statistically significant when Bonferroni adjustment was used. Conclusion: The observed differences in DNA methylation, although nonsignificant after correction for multiple comparisons, suggest that differential regulation of the two sides may impact lip fusion and warrant larger-scale replication.

Syndromale und nichtsyndromale orofaziale Spalten

Orofaziale Spalten sind die zweithäufigste angeborene Fehlbildung. Die beiden häufigen Untergruppen sind die Formen der Lippen-Kiefer-Gaumen-Spalte (LKGS) und der reinen Gaumenspalte (GS). Beide können im Zuge zahlreicher Syndrome auftreten, LKGS sind häufiger nichtsyndromal (ca. 70 %) als GS (ca. 50 %). Nichtsyndromale Spalten sind multifaktoriell und bergen relativ geringe Wiederholungsrisiken für Angehörige. Manche syndromalen Spaltformen beruhen auf Chromosomenaberrationen, die meisten sind jedoch monogen. Das autosomal-dominante Van-der-Woude-Syndrom ist eine wichtige Differenzialdiagnose zu nichtsyndromalen Spaltformen, die sich nur durch Unterlippengrübchen von diesen unterscheidet, aber höhere Risiken für eine Spaltbildung bei Nachkommen birgt. Die Pränataldiagnostik orofazialer Spalten erfolgt sonographisch. Erfahrene Untersucher können LKGS in der 20. Schwangerschaftswoche diagnostizieren; die Detektion von GS ist schwieriger.

Investigation of dominant and recessive inheritance models in genome-wide association studies data of nonsyndromic cleft lip with or without cleft palate

Nonsyndromic cleft lip with or without cleft palate (nsCL/P) is one of the most common congenital malformation worldwide, and its etiology involves both genetic and environmental factors. Recent genome-wide and targeted genetic studies of nsCL/P have identified numerous genetic risk loci, under the hypothesis of a multiplicative mode of inheritance. The present study investigated whether novel nsCL/P risk loci could be identified by analyzing dominant/recessive genetic effects in single nucleotide polymorphism (SNP) data from genome-wide association studies. For this purpose, a genome-wide investigation of dominant/recessive common SNP effects was performed in our previously published meta-analysis data set. Twenty-four loci were identified as candidate regions. In a subsequent association analysis in an independent study cohort of 224 nsCL/P patients and 986 controls of European descent, none of the loci could be replicated. Therefore, our strategy of identifying novel loci by applying different genetic models did not yield any novel findings, suggesting that recessive/dominant common variation only make a limited contribution to nsCL/P in Europeans. However, we cannot rule out that such effects are present at some of the loci that have previously been identified, or are present in different populations.

A new approach to chromosome-wide analysis of X-linked markers identifies new associations in Asian and European case-parent triads of orofacial clefts

Background

GWAS discoveries on the X-chromosome are underrepresented in the literature primarily because the analytical tools that have been applied were originally designed for autosomal markers. Our objective here is to employ a new robust and flexible tool for chromosome-wide analysis of X-linked markers in complex traits. Orofacial clefts are good candidates for such analysis because of the consistently observed excess of females with cleft palate only (CPO) and excess of males with cleft lip with or without cleft palate (CL/P).

Methods

Genotypes for 14,486 X-chromosome SNPs in 1,291 Asian and 1,118 European isolated cleft triads were available from a previously published GWAS. The R-package HAPLIN enables genome-wide–level analyses as well as statistical power simulations for a range of biologic scenarios. We analyzed isolated CL/P and isolated CPO for each ethnicity in HAPLIN, using a sliding-window approach to haplotype analysis and two different statistical models, with and without X-inactivation in females.

Results

There was a larger number of associations in the Asian versus the European sample, and similar to previous reports that have analyzed the same GWAS dataset using different methods, we identified associations with EFNB1/PJA1 and DMD. In addition, new associations were detected with several other genes, among which KLHL4, TBX22, CPXCR1 and BCOR were noteworthy because of their roles in clefting syndromes. A few of the associations were only detected by one particular X-inactivation model, whereas a few others were only detected in one sex.

Discussion/Conclusion

We found new support for the involvement of X-linked variants in isolated clefts. The associations were specific for ethnicity, sex and model parameterization, highlighting the need for flexible tools that are capable of detecting and estimating such effects. Further efforts are needed to verify and elucidate the potential roles of EFNB1/PJA1, KLHL4, TBX22, CPXCR1 and BCOR in isolated clefts.

Candidate Genes for Nonsyndromic Cleft Palate Detected by Exome Sequencing

Nonsyndromic cleft palate only (nsCPO) is a facial malformation that has a livebirth prevalence of 1 in 2,500. Research suggests that the etiology of nsCPO is multifactorial, with a clear genetic component. To date, genome-wide association studies have identified only 1 conclusive common variant for nsCPO, that is, a missense variant in the gene grainyhead-like-3 ( GRHL3). Thus, the underlying genetic causes of nsCPO remain largely unknown. The present study aimed at identifying rare variants that might contribute to nsCPO risk, via whole-exome sequencing (WES), in multiply affected Central European nsCPO pedigrees. WES was performed in 2 affected first-degree relatives from each family. Variants shared between both individuals were analyzed for their potential deleterious nature and a low frequency in the general population. Genes carrying promising variants were annotated for 1) reported associations with facial development, 2) multiple occurrence of variants, and 3) expression in mouse embryonic palatal shelves. This strategy resulted in the identification of a set of 26 candidate genes that were resequenced in 132 independent nsCPO cases and 623 independent controls of 2 different ethnicities, using molecular inversion probes. No rare loss-of-function mutation was identified in either WES or resequencing step. However, we identified 2 or more missense variants predicted to be deleterious in each of 3 genes ( ACACB, PTPRS, MIB1) in individuals from independent families. In addition, the analyses identified a novel variant in GRHL3 in 1 patient and a variant in CREBBP in 2 siblings. Both genes underlie different syndromic forms of CPO. A plausible hypothesis is that the apparently nonsyndromic clefts in these 3 patients might represent hypomorphic forms of the respective syndromes. In summary, the present study identified rare variants that might contribute to nsCPO risk and suggests candidate genes for further investigation.

Breeding experiments and genome-wide association analysis elucidate two genetically different forms of non-syndromic congenital cleft lip and jaw in Vorderwald × Montbéliarde cattle

Non-syndromic congenital cleft lip and jaw (CLJ) is a condition reported in Vorderwald × Montbéliarde cattle. The objective of the present study was to perform a genome-wide association study (GWAS) for 10 CLJ-affected and 50 unaffected Vorderwald × Montbéliarde cattle using the bovine Illumina high density bead chip to identify loci for this condition. Phenotypic classification of CLJ was based on a detailed recording of orofacial structures using computed tomography. A breeding experiment among CLJ-affected Vorderwald × Montbéliarde cattle and CLJ-affected Vorderwald × Montbéliarde cattle with unaffected Holsteins confirmed recessive inheritance and different loci for bilateral or left-sided versus right-sided CLJ. The GWAS for the five cases with right-sided CLJ gave a genome-wide signal on bovine chromosome (BTA) 29 at 16 Mb. For the four left-sided and one bilateral CLJ case, a genome-wide significant association was identified on BTA4 at 32 Mb. Two different loci are very likely to be involved in CLJ in Vorderwald × Montbéliarde cattle because experimental matings among affected cows and bulls with different types of CLJ did not result in CLJ-affected progeny, and in addition, two different loci were also found through GWAS and mapped on two different bovine chromosomes. Validation in 346 Vorderwald × Montbéliarde cattle for the highly associated SNPs on BTA4 and 29 gave ratios of 33/346 (0.095, BTA4) and 6/346 (0.017, BTA29) homozygous mutant genotypes. Further studies should elucidate the responsible mutations underlying the different types of CLJ in Vorderwald × Montbéliarde cattle.

The Role of Noncoding Genetic Variation in Isolated Orofacial Clefts

In the past decade, medical genetic research has generated multiple discoveries, many of which were obtained via genome-wide association studies (GWASs). A major GWAS finding is that the majority of risk variants for complex traits map to noncoding regions. This has resulted in a paradigm shift in terms of the interpretation of human genomic sequence variation, with more attention now being paid to what was previously termed “junk DNA.” Translation of genetic findings into biologically meaningful results requires 1) large-scale and cell-specific efforts to annotate non-protein–coding regions and 2) the integration of comprehensive genomic data sets. However, this represents an enormous challenge, particularly in the case of human traits that arise during embryonic development, such as orofacial clefts (OFCs). OFC is a multifactorial trait and ranks among the most common of all human congenital malformations. These 2 attributes apply in particular to its isolated forms (nonsyndromic OFC [nsOFC]). Although genetic studies (including GWASs) have yielded novel insights into the genetic architecture of nsOFC, few data are available concerning causality and affected biological pathways. Reasons for this deficiency include the complex genetic architecture at risk loci and the limited availability of functional data sets from human tissues that represent relevant embryonic sites and time points. The present review summarizes current knowledge of the role of noncoding regions in nsOFC etiology. We describe the identification of genetic risk factors for nsOFC and several of the approaches used to identify causal variants at these loci. These strategies include the use of biological and genetic information from public databases, the assessment of the full spectrum of genetic variability within 1 locus, and comprehensive in vitro and in vivo experiments. This review also highlights the role of the emerging research field “functional genomics” and its increasing contribution to our biological understanding of nsOFC.

Prevalence of Orofacial Clefts among Live Births in China: A Systematic Review and Meta-Analysis

BACKGROUND

Orofacial clefts (OFCs) are common human birth defects in China. However, studies on the prevalence of OFCs present inconsistent results. The overall prevalence and geographic distribution of OFCs are poorly described in China. Thus, we conducted a systematic review and meta-analysis to estimate the prevalence of OFCs.

METHODS

The systematic review and meta-analysis were conducted on the basis of an established protocol (PROSPERO 2015: CRD42015030198). We systematically searched for articles in four electronic databases, including Embase, PubMed, Wanfang Database, and China National Knowledge Infrastructure (CNKI) to identify relevant studies about prevalence of OFCs in China. Meta-analysis, including subgroup analysis, was conducted to estimate the pooled prevalence.

RESULTS

A total of 41 studies published between 1986 and 2015 were included in our analysis. The sample size ranged from 2,586 to 4,611,808 live births. The random-effects model of meta-analysis showed that the overall prevalence of OFCs in China was 1.4 per 1000 live births (95% confidence interval [CI], 1.1-1.7). In subgroup analysis based on geographic regions, we found that OFC prevalence in Southwest (2.3 per 1000 live births, 95% CI, 1.1-4.7) was higher than that in other regions of China. There were no significant time trends of OFCs during the study period (p-value = 0.47).

CONCLUSION

The overall prevalence of OFCs in China was 1.4 per 1000 live births. No significant secular trend of prevalence has been found in this analysis. Further studies need to be conducted to explore the etiology of OFC to better control the risk of this common birth defect. Birth Defects Research 109:1011-1019, 2017. © 2017 Wiley Periodicals, Inc.

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

Nonsyndromic cleft lip with or without cleft palate (nsCL/P) is among the most common human birth defects with multifactorial etiology. Here, we present results from a genome-wide imputation study of nsCL/P in which, after adding replication cohort data, four novel risk loci for nsCL/P are identified (at chromosomal regions 2p21, 14q22, 15q24 and 19p13). On a systematic level, we show that the association signals within this high-density dataset are enriched in functionally-relevant genomic regions that are active in both human neural crest cells (hNCC) and mouse embryonic craniofacial tissue. This enrichment is also detectable in hNCC regions primed for later activity. Using GCTA analyses, we suggest that 30% of the estimated variance in risk for nsCL/P in the European population can be attributed to common variants, with 25.5% contributed to by the 24 risk loci known to date. For each of these, we identify credible SNPs using a Bayesian refinement approach, with two loci harbouring only one probable causal variant. Finally, we demonstrate that there is no polygenic component of nsCL/P detectable that is shared with nonsyndromic cleft palate only (nsCPO). Our data suggest that, while common variants are strongly contributing to risk for nsCL/P, they do not seem to be involved in nsCPO which might be more often caused by rare deleterious variants. Our study generates novel insights into both nsCL/P and nsCPO etiology and provides a systematic framework for research into craniofacial development and malformation.