Recurrence risk for offspring of twins discordant for oral cleft: a population-based cohort study of the Danish 1936-2004 cleft twin cohort

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.

Wavelet Screening identifies regions highly enriched for differentially methylated loci for orofacial clefts

DNA methylation is the most widely studied epigenetic mark in humans and plays an essential role in normal biological processes as well as in disease development. More focus has recently been placed on understanding functional aspects of methylation, prompting the development of methods to investigate the relationship between heterogeneity in methylation patterns and disease risk. However, most of these methods are limited in that they use simplified models that may rely on arbitrarily chosen parameters, they can only detect differentially methylated regions (DMRs) one at a time, or they are computationally intensive. To address these shortcomings, we present a wavelet-based method called 'Wavelet Screening' (WS) that can perform an epigenome-wide association study (EWAS) of thousands of individuals on a single CPU in only a matter of hours. By detecting multiple DMRs located near each other, WS identifies more complex patterns that can differentiate between different methylation profiles. We performed an extensive set of simulations to demonstrate the robustness and high power of WS, before applying it to a previously published EWAS dataset of orofacial clefts (OFCs). WS identified 82 associated regions containing several known genes and loci for OFCs, while other findings are novel and warrant replication in other OFCs cohorts.

The Prevalence of Nonsyndromic Oral Clefts in Twins Compared to Singletons: The Association With Birth Weight

PURPOSE

To examine the prevalence of nonsyndromic oral clefts in twins compared to singletons in the United States and to evaluate the association between birth weight and nonsyndromic oral clefts.

DESIGN

A large population-based cross-sectional study was performed using the data from the US National Center for Health Statistics database in 2017.

PARTICIPANTS

Our sample consisted of 128 310 twins and 3 723 273 singletons.

METHODS

The variables collected were sociodemographic variables, environmental predictors, and clinical measures. Descriptive analysis, bivariate, and multivariate logistic regression were performed.

MAIN OUTCOME MEASURE

The main outcome variable in our study is nonsyndromic oral clefts.

RESULTS

The prevalence of nonsyndromic oral clefts was 5.22 per 10 000 in twins and 5.12 per 10 000 in singletons. Results show no significant risk of nonsyndromic oral clefts in twins compared to singletons (P = .92). There was a significant relationship between birth weight and infant diagnosed with nonsyndromic oral clefts (P = .01). Unadjusted odds ratio for birth weight was 2.52 (95% CI: 2.25-2.82). Adjusted odds for potential confounders such as mother's age, race, mother's education, gender of the infant, APGAR 5-minute score, gestational age, prenatal smoking, number of prenatal care visits, and mother's body mass index were resulted in similar but with a slightly lower odds of 2.11 (95% CI: 1.78-2.50).

CONCLUSION

Compared to singletons, twins did not have higher risk of nonsyndromic oral clefts. Infants with low birth weight were more prone to have nonsyndromic oral clefts.

Genetic Factors in Nonsyndromic Orofacial Clefts

Orofacial clefts (OFCs) are the most common congenital birth defects in humans and immediately recognized at birth. The etiology remains complex and poorly understood and seems to result from multiple genetic and environmental factors along with gene-environment interactions. It can be classified into syndromic (30%) and nonsyndromic (70%) clefts. Nonsyndromic OFCs include clefts without any additional physical or cognitive deficits. Recently, various genetic approaches, such as genome-wide association studies (GWAS), candidate gene association studies, and linkage analysis, have identified multiple genes involved in the etiology of OFCs. This article provides an insight into the multiple genes involved in the etiology of OFCs. Identification of specific genetic causes of clefts helps in a better understanding of the molecular pathogenesis of OFC. In the near future, it helps to provide a more accurate diagnosis, genetic counseling, personalized medicine for better clinical care, and prevention of OFCs.

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.

A genome-wide scan of cleft lip triads identifies parent-of-origin interaction effects between ANK3 and maternal smoking, and between ARHGEF10 and alcohol consumption

Background: Although both genetic and environmental factors have been reported to influence the risk of isolated cleft lip with or without cleft palate (CL/P), the exact mechanisms behind CL/P are still largely unaccounted for. We recently developed new methods to identify parent-of-origin (PoO) interactions with environmental exposures (PoOxE) and now apply them to data from a genome-wide association study (GWAS) of families with children born with isolated CL/P. Methods: Genotypes from 1594 complete triads and 314 dyads (1908 nuclear families in total) with CL/P were available for the current analyses. Of these families, 1024 were Asian, 825 were European and 59 had other ancestries. After quality control, 341,191 SNPs remained from the original 569,244. The exposures were maternal cigarette smoking, use of alcohol, and use of vitamin supplements in the periconceptional period. Our new methodology detects if PoO effects are different across environmental strata and is implemented in the R-package Haplin. Results: Among Europeans, there was evidence of a PoOxSmoke effect for ANK3 with three SNPs (rs3793861, q=0.20, p=2.6e-6; rs7087489, q=0.20, p=3.1e-6; rs4310561, q=0.67, p=4.0e-5) and a PoOxAlcohol effect for ARHGEF10 with two SNPs (rs2294035, q=0.32, p=2.9e-6; rs4876274, q=0.76, p=1.3e-5). Conclusion: Our results indicate that the detected PoOxE effects have a plausible biological basis, and thus warrant replication in other independent cleft samples. Our demonstration of the feasibility of identifying complex interactions between relevant environmental exposures and PoO effects offers new avenues for future research aimed at unravelling the complex etiology of cleft lip defects.

A genome-wide scan of cleft lip triads identifies parent-of-origin interaction effects between ANK3 and maternal smoking, and between ARHGEF10 and alcohol consumption

Background: Although both genetic and environmental factors have been reported to influence the risk of isolated cleft lip with or without cleft palate (CL/P), the exact mechanisms behind CL/P are still largely unaccounted for. We recently developed new methods to identify parent-of-origin (PoO) interactions with environmental exposures (PoOxE) and applied them to families with children born with isolated cleft palate only. Here, we used the same genome-wide association study (GWAS) dataset and methodology to screen for PoOxE effects in the larger sample of CL/P triads. Methods: Genotypes from 1594 complete triads and 314 dyads (1908 nuclear families in total) with CL/P were available for the current analyses. Of these families, 1024 were Asian, 825 were European and 59 had other ancestries. After quality control, 341,191 SNPs remained from the original 569,244. The exposures were maternal cigarette smoking, use of alcohol, and use of vitamin supplements in the periconceptional period. The methodology applied in the analyses is implemented in the R-package Haplin. Results: Among Europeans, there was evidence of a PoOxSmoke effect for ANK3 with three SNPs (rs3793861, q=0.20, p=2.6e-6; rs7087489, q=0.20, p=3.1e-6; rs4310561, q=0.67, p=4.0e-5) and a PoOxAlcohol effect for ARHGEF10 with two SNPs (rs2294035, q=0.32, p=2.9e-6; rs4876274, q=0.76, p=1.3e-5). Conclusion: Our results indicate that the detected PoOxE effects have a plausible biological basis, and thus warrant replication in other independent cleft samples. Our demonstration of the feasibility of identifying complex interactions between relevant environmental exposures and PoO effects offers new avenues for future research aimed at unravelling  the complex etiology of cleft lip defects.

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.

Genome-wide analysis of parent-of-origin interaction effects with environmental exposure (PoOxE): An application to European and Asian cleft palate trios

Cleft palate only is a common birth defect with high heritability. Only a small fraction of this heritability is explained by the genetic variants identified so far, underscoring the need to investigate other disease mechanisms, such as gene-environment (GxE) interactions and parent-of-origin (PoO) effects. Furthermore, PoO effects may vary across exposure levels (PoOxE effects). Such variation is the focus of this study. We upgraded the R-package Haplin to enable direct tests of PoOxE effects at the genome-wide level. From a previous GWAS, we had genotypes for 550 case-parent trios, of mainly European and Asian ancestry, and data on three maternal exposures (smoking, alcohol, and vitamins). Data were analyzed for Europeans and Asians separately, and also for all ethnicities combined. To account for multiple testing, a false discovery rate method was used, where q-values were generated from the p-values. In the Europeans-only analyses, interactions with maternal smoking yielded the lowest q-values. Two SNPs in the 'Interactor of little elongation complex ELL subunit 1' (ICE1) gene had a q-value of 0.14, and five of the 20 most significant SNPs were in the 'N-acetylated alpha-linked acidic dipeptidase-like 2' (NAALADL2) gene. No evidence of PoOxE effects was found in the other analyses. The connections to ICE1 and NAALADL2 are novel and warrant further investigation. More generally, the new methodology presented here is easily applicable to other traits and exposures in which a family-based study design has been implemented.

Testing the face shape hypothesis in twins discordant for nonsyndromic orofacial clefting

Nonsyndromic orofacial clefts (OFCs) are complex traits characterized by multifactorial inheritance and wide phenotypic variability. Numerous studies have shown subtle differences in the faces of unaffected relatives from cleft families compared to controls, the implication being that such outward differences are an incomplete expression reflecting an underlying genetic predisposition. Twins discordant for OFCs provide a unique opportunity to further test this idea, as the unaffected co-twin shares on average 50% (for dizygotic twins) and 100% (for monozygotic twins) of the genetic risk factors as the affected twin. We used 3D surface imaging and spatially-dense morphometry to compare facial shape in a sample of 44 unaffected co-twins and age- and sex-matched unaffected controls (n = 241). Unaffected co-twins showed statistically significant differences in the midface, lateral upper face, and forehead regions, compared to controls. Furthermore, co-twins were characterized by a distinct pattern of midfacial retrusion, broader upper faces, and greater protrusion of the mandible and brow ridges. This same general facial pattern was shown in both unaffected monozygotic and dizygotic co-twin subsets. These results provide additional support that altered facial shape is a phenotypic marker for OFC susceptibility.

Exploring Subclinical Phenotypic Features in Twin Pairs Discordant for Cleft Lip and Palate

OBJECTIVE

Monozygotic twins of an individual with an orofacial cleft have a significantly elevated risk for orofacial cleft compared with the general population, but still the concordance rate for orofacial cleft in monozygotic twins is about 40% to 50%. The goal of this study was to determine whether unaffected cotwins have an increased frequency of orbicularis oris muscle defects, a subclinical form of orofacial cleft. The presence of such defects may reduce the overall rate of discordance.

METHOD

A total of 63 discordant monozygotic and dizygotic twin pairs, 262 unaffected nontwin siblings, and 543 controls with no history of orofacial clefts were assessed for orbicularis oris defects by high-resolution ultrasound. Frequencies were compared by the Fisher exact test.

RESULTS

Unaffected cotwins from discordant monozygotic pairs had a higher frequency of defects (12.5%) than the other test groups (6.38% to 6.99%), but the difference was not statistically significant (P = .74).

CONCLUSIONS

In this study, orbicularis oris defects were not statistically significantly more common among the unaffected twins from orofacial cleft discordant twin pairs. The trends in the results warrant future studies with larger sample sizes and additional subclinical phenotypes.

From Bench to Bedside and Back: Improving Diagnosis and Treatment of Craniofacial Malformations Utilizing Animal Models

Craniofacial anomalies are among the most common birth defects and are associated with increased mortality and, in many cases, the need for lifelong treatment. Over the past few decades, dramatic advances in the surgical and medical care of these patients have led to marked improvements in patient outcomes. However, none of the treatments currently in clinical use address the underlying molecular causes of these disorders. Fortunately, the field of craniofacial developmental biology provides a strong foundation for improved diagnosis and for therapies that target the genetic causes of birth defects. In this chapter, we discuss recent advances in our understanding of the embryology of craniofacial conditions, and we focus on the use of animal models to guide rational therapies anchored in genetics and biochemistry.

Regulatory variant in FZD6 gene contributes to nonsyndromic cleft lip and palate in an African-American family

Nonsyndromic cleft lip with or without cleft palate (NSCLP) is a common birth defect affecting 135,000 newborns worldwide each year. While a multifactorial etiology has been suggested as the cause, despite decades of research, the genetic underpinnings of NSCLP remain largely unexplained. In our previous genome-wide linkage study of a large NSCLP African-American family, we identified a candidate locus at 8q21.3-24.12 (LOD = 2.98). This region contained four genes, Frizzled-6 (FZD6), Matrilin-2 (MATN2), Odd-skipped related 2 (OSR2) and Solute Carrier Family 25, Member 32 (SLC25A32). FZD6 was located under the maximum linkage peak. In this study, we sequenced the coding and noncoding regions of these genes in two affected family members, and identified a rare variant in intron 1 of FZD6 (rs138557689; c.-153 + 432A>C). The variant C allele segregated with NSCLP in this family, through affected and unaffected individuals, and was found in one other NSCLP African-American family. Functional assays showed that this allele creates an allele-specific protein-binding site and decreases promoter activity. We also observed that loss and gain of fzd6 in zebrafish contributes to craniofacial anomalies. FZD6 regulates the WNT signaling pathway, which is involved in craniofacial development, including midfacial formation and upper labial fusion. We hypothesize, therefore, that alteration in FZD6 expression contributes to NSCLP in this family by perturbing the WNT signaling pathway.

MicroRNA-17-92, a direct Ap-2α transcriptional target, modulates T-box factor activity in orofacial clefting

Among the most common human congenital anomalies, cleft lip and palate (CL/P) affects up to 1 in 700 live births. MicroRNA (miR)s are small, non-coding RNAs that repress gene expression post-transcriptionally. The miR-17-92 cluster encodes six miRs that have been implicated in human cancers and heart development. We discovered that miR-17-92 mutant embryos had severe craniofacial phenotypes, including incompletely penetrant CL/P and mandibular hypoplasia. Embryos that were compound mutant for miR-17-92 and the related miR-106b-25 cluster had completely penetrant CL/P. Expression of Tbx1 and Tbx3, the DiGeorge/velo-cardio-facial (DGS) and Ulnar-mammary syndrome (UMS) disease genes, was expanded in miR-17-92 mutant craniofacial structures. Both Tbx1 and Tbx3 had functional miR seed sequences that mediated gene repression. Analysis of miR-17-92 regulatory regions uncovered conserved and functional AP-2α recognition elements that directed miR-17-92 expression. Together, our data indicate that miR-17-92 modulates expression of critical T-box transcriptional regulators during midface development and is itself a target of Bmp-signaling and the craniofacial pioneer factor AP-2α. Our data are the first genetic evidence that an individual miR or miR cluster is functionally important in mammalian CL/P.

Genome-wide approaches (GWA) in oral and craniofacial diseases research

Oral Diseases (2012) Underlying molecular genetic mechanisms of diseases can be deciphered with unbiased strategies using recently developed technologies enabling genome-wide scale investigations. These technologies have been applied in scanning for genetic variations, gene expression profiles, and epigenetic changes for oral and craniofacial diseases. However, these approaches as applied to oral and craniofacial conditions are in the initial stages, and challenges remain to be overcome, including analysis of high throughput data and their interpretation. Here, we review methodology and studies using genome-wide approaches in oral and craniofacial diseases and suggest future directions.

Whole-blood 3-hydroxyisovalerylcarnitine as a risk factor for orofacial clefts

BACKGROUND/PURPOSE

In mice, biotin deficiency is one of the most potent clefting factors. Increased 3-hydroxyisovalerylcarnitine (C5OH) is regarded as a biomarker of biotin deficiency. This retrospective study was undertaken to determine whether increased C5OH in newborns is associated with orofacial clefts.

MATERIALS AND METHODS

Seventy newborns with non-syndromic cleft lip with or without cleft palate and 140 control newborns without congenital anomalies were investigated. Whole-blood C5OH concentrations were measured using tandem mass spectrometry.

RESULTS

The median (interquartile range, IQR) concentrations of C5OH in patients with clefts and controls were 0.16 (0.13-0.22)μmoll(-1) and 0.17 (0.13-0.20)μmoll(-1), respectively (p=0.90). The receiver operating characteristic analysis did not find out cut-off values for C5OH discriminating between cases and controls.

CONCLUSION

There appears to be no association between biotin deficiency, as indexed by an increase of C5OH, and orofacial clefts in the investigated group of patients.

The evolution of human genetic studies of cleft lip and cleft palate

Orofacial clefts (OFCs)--primarily cleft lip and cleft palate--are among the most common birth defects in all populations worldwide, and have notable population, ethnicity, and gender differences in birth prevalence. Interest in these birth defects goes back centuries, as does formal scientific interest; scientists often used OFCs as examples or evidence during paradigm shifts in human genetics, and have also used virtually every new method of human genetic analysis to deepen our understanding of OFC. This review traces the evolution of human genetic investigations of OFC, highlights the specific insights gained about OFC through the years, and culminates in a review of recent key OFC genetic findings resulting from the powerful tools of the genomics era. Notably, OFC represents a major success for genome-wide approaches, and the field is poised for further breakthroughs in the near future.