Meta-analysis Reveals Genome-Wide Significance at 15q13 for Nonsyndromic Clefting of Both the Lip and the Palate, and Functional Analyses Implicate GREM1 As a Plausible Causative Gene

The heterogeneous genetic architectures of orofacial clefts

Orofacial clefts (OFCs) are common, affecting 1:1000 live births. OFCs occur across a phenotypic spectrum - including cleft lip (CL), cleft lip and palate (CLP), or cleft palate (CP) - and can be further subdivided based on laterality, severity, or specific structures affected. Herein we review what is known about the genetic architecture underlying each of these subtypes, considering both shared and subtype-specific risks. While there are more known genetic similarities between CL and CLP than CP, recent research supports both shared and subtype-specific genetic risk factors within and between phenotypic classifications of OFCs. Larger sample sizes and deeper phenotyping data will be of increasing importance for the discovery of novel genetic risk factors for OFCs and various subtypes going forward.

Ezh2-dependent methylation in oral epithelia promotes secondary palatogenesis

BACKGROUND

In addition to genomic risk variants and environmental influences, increasing evidence suggests epigenetic modifications are important for orofacial development and their alterations can contribute to orofacial clefts. Ezh2 encodes a core catalytic component of the Polycomb repressive complex responsible for addition of methyl marks to Histone H3 as a mechanism of repressing target genes. The role of Ezh2 in orofacial clefts remains unknown.

AIMS

To investigate the epithelial role of Ezh2-dependent methylation in secondary palatogenesis.

METHODS

We used conditional gene-targeting methods to ablate Ezh2 in the surface ectoderm-derived oral epithelium of mouse embryos. We then performed single-cell RNA sequencing combined with immunofluorescence and RT-qPCR to investigate gene expression in conditional mutant palate. We also employed double knockout analyses of Ezh1 and Ezh2 to address if they have synergistic roles in palatogenesis.

RESULTS

We found that conditional inactivation of Ezh2 in oral epithelia results in partially penetrant cleft palate. Double knockout analyses revealed that another family member Ezh1 is dispensable in orofacial development, and it does not have synergistic roles with Ezh2 in palatogenesis. Histochemistry and single-cell RNA-seq analyses revealed dysregulation of cell cycle regulators in the palatal epithelia of Ezh2 mutant mouse embryos disrupts palatogenesis.

CONCLUSION

Ezh2-dependent histone H3K27 methylation represses expression of cell cycle regulator Cdkn1a and promotes proliferation in the epithelium of the developing palatal shelves. Loss of this regulation may perturb movement of the palatal shelves, causing a delay in palate elevation which may result in failure of the secondary palate to close altogether.

Association between IRF6, TP63, GREM1 Gene Polymorphisms and Non-Syndromic Orofacial Cleft Phenotypes in Vietnamese Population: A Case-Control and Family-Based Study

This study aims to identify potential variants in the TP63-IRF6 pathway and GREM1 for the etiology of non-syndromic orofacial cleft (NSOFC) among the Vietnamese population. By collecting 527 case-parent trios and 527 control samples, we conducted a stratified analysis based on different NSOFC phenotypes, using allelic, dominant, recessive and over-dominant models for case-control analyses, and family-based association tests for case-parent trios. Haplotype and linkage disequilibrium analyses were also conducted. IRF6 rs2235375 showed a significant association with an increased risk for non-syndromic cleft lip and palate (NSCLP) and cleft lip with or without cleft palate (NSCL/P) in the G allele, with pallele values of 0.0018 and 0.0003, respectively. Due to the recessive model (p = 0.0011) for the NSCL/P group, the reduced frequency of the GG genotype of rs2235375 was associated with a protective effect against NSCL/P. Additionally, offspring who inherited the G allele at rs2235375 had a 1.34-fold increased risk of NSCL/P compared to the C allele holders. IRF6 rs846810 and a G-G haplotype at rs2235375-rs846810 of IRF6 impacted NSCL/P, with p-values of 0.0015 and 0.0003, respectively. In conclusion, our study provided additional evidence for the association of IRF6 rs2235375 with NSCLP and NSCL/P. We also identified IRF6 rs846810 as a novel marker associated with NSCL/P, and haplotypes G-G and C-A at rs2235375-rs846810 of IRF6 associated with NSOFC.

Genome-wide meta-analyses identify five new risk loci for nonsyndromic orofacial clefts in the Chinese Han population

BACKGROUND

Nonsyndromic orofacial clefts (NSOFCs) are the most common craniofacial birth malformations in humans and are generally classified as nonsyndromic cleft lip with or without cleft palate (NSCL/P) and nonsyndromic cleft palate only (NSCPO). Genome-wide association studies (GWASs) of NSOFCs have demonstrated multiple risk loci and candidate genes; however, published risk factors are able to explain only a small fraction of the observed NSOFCs heritability.

METHODS

Here, we performed GWASs of 1615 NSCPO cases and 2340 controls, and then conducted genome-wide meta-analyses of NSOFCs, totaling 6812 NSCL/P cases, 2614 NSCPO cases, and 19,165 controls from the Chinese Han population.

RESULTS

We identify 47 risk loci with genome-wide pmeta -value <5.0 × 10-8 , 5 risk loci (1p32.1, 3p14.1, 3p14.3, 3p21.31, and 13q22.1) of which are new. All of the 47 susceptibility loci conjointly account for 44.12% of the NSOFCs' heritability in the Chinese Han population.

CONCLUSION

Our results improve the comprehending of genetic susceptibility to NSOFCs and provide new views into the genetic etiology of craniofacial anomalies.

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

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

Translational Research Techniques for the Facial Plastic Surgeon: An Overview

The field of facial plastic and reconstructive surgery (FPRS) is an incredibly diverse, multispecialty field that seeks innovative and novel solutions for the management of physical defects on the head and neck. To aid in the advancement of medical and surgical treatments for these defects, there has been a recent emphasis on the importance of translational research. With recent technological advancements, there are now a myriad of research techniques that are widely accessible for physician and scientist use in translational research. Such techniques include integrated multiomics, advanced cell culture and microfluidic tissue models, established animal models, and emerging computer models generated using bioinformatics. This study discusses these various research techniques and how they have and can be used for research in the context of various important diseases within the field of FPRS.

Rare genetic variants in SEC24D modify orofacial cleft phenotypes

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

Genome Analysis Using Whole-Exome Sequencing of Non-Syndromic Cleft Lip and/or Palate from Malagasy Trios Identifies Variants Associated with Cilium-Related Pathways and Asian Genetic Ancestry

Background: Orofacial clefts (OFCs) are common congenital disabilities that can occur as isolated non-syndromic events or as part of Mendelian syndromes. OFC risk factors vary due to differences in regional environmental exposures, genetic variants, and ethnicities. In recent years, significant progress has been made in understanding OFCs, due to advances in sequencing and genotyping technologies. Despite these advances, very little is known about the genetic interplay in the Malagasy population. Methods: Here, we performed high-resolution whole-exome sequencing (WES) on non-syndromic cleft lip with or without palate (nCL/P) trios in the Malagasy population (78 individuals from 26 families (trios)). To integrate the impact of genetic ancestry admixture, we computed both global and local ancestries. Results: Participants demonstrated a high percentage of both African and Asian admixture. We identified damaging variants in primary cilium-mediated pathway genes WNT5B (one family), GPC4 (one family), co-occurrence in MSX1 (five families), WDR11 (one family), and tubulin stabilizer SEPTIN9 (one family). Furthermore, we identified an autosomal homozygous damaging variant in PHGDH (one family) gene that may impact metabiotic activity. Lastly, all variants were predicted to reside on local Asian genetic ancestry admixed alleles. Conclusion: Our results from examining the Malagasy genome provide limited support for the hypothesis that germline variants in primary cilia may be risk factors for nCL/P, and outline the importance of integrating local ancestry components better to understand the multi-ethnic impact on nCL/P.

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.

Gene×environment associations in orofacial clefting

This chapter reviews the evidence of gene×environment interactions (G×E) in the etiology of orofacial cleft birth defects (OFCs), specifically cleft lip (CL), cleft palate (CP), and cleft lip with or without cleft palate (CL/P). We summarize the current state of our understanding of the genetic architecture of nonsyndromic OFCs and the evidence that maternal exposures during pregnancy influence risk of OFCs. Further, we present possible candidate gene pathways for these exposures including metabolism of folates, metabolism of retinoids, retinoic acid receptor signaling, aryl hydrocarbon receptor signaling, glucocorticoid receptor signaling, and biotransformation and transport. We review genes in these pathways with prior evidence of association with OFCs, genes with evidence from prior candidate gene G×E studies, and genes identified from genome-wide searches specifically for identifying G×E. Finally, we suggest future directions for G×E research in OFCs.

Polygenic risk impacts PDGFRA mutation penetrance in non-syndromic cleft lip and palate

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is a common, severe craniofacial malformation that imposes significant medical, psychosocial and financial burdens. NSCL/P is a multifactorial disorder with genetic and environmental factors playing etiologic roles. Currently, only 25% of the genetic variation underlying NSCL/P has been identified by linkage, candidate gene and genome-wide association studies. In this study, whole-genome sequencing and genome-wide genotyping followed by polygenic risk score (PRS) and linkage analyses were used to identify the genetic etiology of NSCL/P in a large three-generation family. We identified a rare missense variant in PDGFRA (c.C2740T; p.R914W) as potentially etiologic in a gene-based association test using pVAAST (P = 1.78 × 10-4) and showed decreased penetrance. PRS analysis suggested that variant penetrance was likely modified by common NSCL/P risk variants, with lower scores found among unaffected carriers. Linkage analysis provided additional support for PRS-modified penetrance, with a 7.4-fold increase in likelihood after conditioning on PRS. Functional characterization experiments showed that the putatively causal variant was null for signaling activity in vitro; further, perturbation of pdgfra in zebrafish embryos resulted in unilateral orofacial clefting. Our findings show that a rare PDGFRA variant, modified by additional common NSCL/P risk variants, have a profound effect on NSCL/P risk. These data provide compelling evidence for multifactorial inheritance long postulated to underlie NSCL/P and may explain some unusual familial patterns.

Identification of novel susceptibility genes for non-syndromic cleft lip with or without cleft palate using NGS-based multigene panel testing

For non-syndromic cleft lip with or without cleft palate (ns-CL/P), the proportion of heritability explained by the known risk loci is estimated to be about 30% and is captured mainly by common variants identified in genome-wide association studies. To contribute to the explanation of the "missing heritability" problem for orofacial clefts, a candidate gene approach was taken to investigate the potential role of rare and private variants in the ns-CL/P risk. Using the next-generation sequencing technology, the coding sequence of a set of 423 candidate genes was analysed in 135 patients from the Polish population. After stringent multistage filtering, 37 rare coding and splicing variants of 28 genes were identified. 35% of these genetic alternations that may play a role of genetic modifiers influencing an individual's risk were detected in genes not previously associated with the ns-CL/P susceptibility, including COL11A1, COL17A1, DLX1, EFTUD2, FGF4, FGF8, FLNB, JAG1, NOTCH2, SHH, WNT5A and WNT9A. Significant enrichment of rare alleles in ns-CL/P patients compared with controls was also demonstrated for ARHGAP29, CHD7, COL17A1, FGF12, GAD1 and SATB2. In addition, analysis of panoramic radiographs of patients with identified predisposing variants may support the hypothesis of a common genetic link between orofacial clefts and dental abnormalities. In conclusion, our study has confirmed that rare coding variants might contribute to the genetic architecture of ns-CL/P. Since only single predisposing variants were identified in novel cleft susceptibility genes, future research will be required to confirm and fully understand their role in the aetiology of ns-CL/P.

Association of soluble epoxide hydrolase 2 gene with the risk of non-syndromic cleft lip with or without cleft palate in western Han Chinese population

OBJECTIVES

This study aimed to explore the associations between soluble epoxide hydrolase 2 gene (EPHX2) variants and non-syndromic cleft lip with or without cleft palate (NSCL/P) in Chinese Han population.

METHODS

We recruited 159 NSCL/P cases from Chinese Han population and carried out targeted resequencing using the whole genome sequencing data of 542 healthy Chinese individuals from Novegene internal database as controls. We classified EPHX2 variants as common or rare according to their minor allele frequency and performed an association analysis for common variations and a burden analysis for rare variations.

RESULTS

The lowest P-value in NSCL/P was observed at rs57699806 (P=0.000 13, OR=2.849 and 95% CI: 1.691-4.800), followed by rs4732723 (P=0.006 50, OR=0.662 and 95%CI: 0.491-0.892), rs7829267 (P=0.009 20, OR=1.496 and 95%CI: 1.117-2.005), rs721619 (P=0.011 00, OR=1.474 and 95%CI: 1.098-1.980), and rs7816586 (P=0.040 00, OR=1.310 and 95%CI: 1.015-1.691). The odds ratios suggested the C allele at rs4732723 as a protective factor for NSCL/P and the reference alleles at other single nucleotide polymorphisms (SNPs) as the risk factors for NSCL/P. Burden analysis showed no statistical significance (P>0.05).

CONCLUSIONS

Through targeted resequencing, this study identified five SNPs named rs57699806, rs4732723, rs7829267, rs721619, and rs7816586 around the region of EPHX2 gene associated with NSCL/P in Chinese Han population. Four SNPs of rs57699806, rs4732723, rs7829267, and rs7816586 were first identified.

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

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

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

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

Exploring the Molecular Mechanism of lncRNA-miRNA-mRNA Networks in Non-Syndromic Cleft Lip with or without Cleft Palate

Background

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is a common craniofacial birth defect. Growing evidence has demonstrated the competing endogenous RNA (ceRNA) hypothesis has played a role in the pathogenesis of NSCL/P. Here, we identified the important lncRNAs in NSCL/P and constructed a ceRNA regulatory network to predict their underlying functional mechanism.

Methods

Total RNA isolated from the peripheral blood samples were analyzed by the Human Clariom D Affymetrix platform and differentially expressed genes (DEGs) were identified. Using the limma package in R software, DEGs in the expression profile of GSE42589 were identified from Gene Expression Omnibus (GEO) database. Co-differentially expressed lncRNAs (co-DElncRNAs) were used to predict the microRNAs that may bind to them. Co-differentially expressed mRNAs (co-DEmRNAs) were subjected to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. The hub genes were screened using the cytohubba plug-in in Cytoscape. A ceRNA network was built to investigate the molecular mechanism underlying the etiology of NSCL/P. The expression levels of lncRNAs, miRNAs, and mRNAs in the network were assessed by quantitative real-time polymerase chain reaction (qRT-PCR).

Results

We found 116 DElncRNAs and 2955 DEmRNAs from the GSE42589 dataset, and 2626 DElncRNAs and 2771 DEmRNAs from the Human Clariom D gene chip. A network of co-DEmRNAs containing 3712 edges and 621 nodes were identified by PPI analysis. A ceRNA regulatory network comprising lncRNA USP17L6P, hsa-miR-449c-5p, and MYC was established. qRT-PCR results revealed significantly lower expression levels of lncRNA USP17L6P and c-Myc in NSCL/P tissues, while the expression level of hsa-miR-449c-5p was higher as compared to control samples (p < 0.05).

Conclusion

The identified lncRNAs and the established ceRNA regulatory network provide novel insight into the pathogenesis of NSCL/P, therefore hold great promise in NSCL/P management in clinical practice.

Genetic models and approaches to study orofacial clefts

INTRODUCTION

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

OBJECTIVE

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

CONCLUSION

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

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

BACKGROUND

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

MATERIAL AND METHODS

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

RESULTS

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

CONCLUSION

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

Genetic architecture of orbital telorism

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

Family based and case-control designs reveal an association of TFAP2A in nonsyndromic cleft lip only among Vietnamese population

Aims

Dozens of causative genes and their mechanisms of nonsyndromic cleft lip with or without cleft palate (NSCL/P) were revealed through genome‐wide association and linkage studies. Results were, however, not always replicated in different populations or methodologies. This study used case–control and family based approaches to investigate the etiology of NSCL/P and its two subtypes: nonsyndromic cleft lip only (NSCLO) and nonsyndromic cleft lip and palate (NSCLP) among the Vietnamese population.

Methods

Two hundred and seventeen NSCL/P case‐parent trios (one affected child and two parents), including 105 NSCLO and 112 NSCLP were involved for a family based design; and 273 ethnic and region‐matched healthy controls with no cleft history in their families were recruited for a case–control design. Three SNPs consisting of TFAP2A (rs1675414 and rs303048) and 8q24 (rs987525) were genotyped using the TaqMan SNP genotyping assay.

Results

TFAP2A rs1675414 was associated with NSCLO, replicated by both case‐control and family based tests. Other SNPs yielded no evidence of susceptibility to NSCL/P or two subtypes.

Conclusion

The current investigation suggests an intriguing role of TFAP2A in the etiology of NSCLO among the Vietnamese population.

This study used case‐control and family‐based approaches to investigate the etiology of NSCL/P and its two subtypes: nonsyndromic cleft lip only (NSCLO), nonsyndromic cleft lip and palate (NSCLP) among Vietnamese population. TFAP2A rs1675414 was associated with NSCLO, replicated by both case‐control and family‐based tests.

Association analysis of SNPs in GRHL3, FAF1, and KCNJ2 with NSCPO sub-phenotypes in Han Chinese

OBJECTIVES

Non-syndromic cleft palate only (NSCPO) is a common congenital deformity with complex etiologies. GRHL3, FAF1, and KCNJ2 have been reported to be involved in the pathogenesis of NSCPO. Up till now, there have been no replication studies based on large Han Chinese. Therefore, this study aimed to investigate associations between GRHL3, FAF1, KCNJ2, and NSCPO sub-phenotypes patients in Han Chinese.

MATERIALS AND METHODS

Firstly, we selected 2 SNPs based on previous literatures: FAF1 (rs3827730) and GRHL3 (rs41268753). Also, we selected 8 tagSNPs in GRHL3 (rs557811, rs609352, rs10903078, rs6659209, rs12401714, rs12568599, rs3887581, rs12024148) and 2 tagSNPs in KCNJ2 (rs75855040 and rs236514). Afterward, we evaluated these SNPs among 1668 NSCPO patients and 1811 normal controls from Han Chinese. Following data were analyzed by PLINK and Haploview program.

RESULTS

Association analysis under additive model showed that allele A at rs12568599 in GRHL3 gene is significantly associated with NSCPO (p = 0.0034, OR = 1.38 and 95%CI: 1.11-1.72) and its sub-phenotype incomplete cleft palate (ICP) (p = 0.0039, OR = 1.4 and 95%CI: 1.11-1.75), and it could increase the risk of both NSCPO and ICP.

CONCLUSIONS

This study firstly found that rs12568599 in GRHL3 is associated with NSCPO and ICP in Han Chinese, indicating that sub-phenotypes of NSCPO have different etiologies.

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.

FAT4 identified as a potential modifier of orofacial cleft laterality

Orofacial clefts (OFCs) are common (1 in 700 births) congenital malformations that include a cleft lip (CL) and cleft lip and palate (CLP). These OFC subtypes are also heterogeneous themselves, with the CL occurring on the left, right, or both sides of the upper lip. Unilateral CL and CLP have a 2:1 bias towards left-sided clefts, suggesting a nonrandom process. Here, we performed a study of left- and right-sided clefts within the CL and CLP subtypes to better understand the genetic factors controlling cleft laterality. We conducted genome-wide modifier analyses by comparing cases that had right unilateral CL (RCL; N = 130), left unilateral CL (LCL; N = 216), right unilateral CLP (RCLP; N = 416), or left unilateral CLP (LCLP; N = 638), and identified a candidate region on 4q28, 400 kb downstream from FAT4, that approached genome-wide significance for LCL versus RCL (p = 8.4 × 10^-8 ). Consistent with its potential involvement as a genetic modifier of CL, we found that Fat4 exhibits a specific domain of expression in the mesenchyme of the medial nasal processes that form the median upper lip. Overall, these results suggest that the epidemiological similarities in left- to right-sided clefts in CL and CLP are not reflected in the genetic association results.

Elevated diversity in loci linked to facial morphology is consistent with the hypothesis that individual facial recognition is important across hominoids

OBJECTIVES

The ability to use visual signals to identify individuals is an important feature of primate social groups, including humans. Sheehan and Nachman (2014) showed that loci linked to facial morphology had elevated levels of diversity and interpreted this as evidence that the human face is under frequency-dependent selection to enhance individual recognition (Nature Communications 5). In our study, we tested whether this pattern is found in non-human ape species, to help understand whether individual recognition might also play a role in species other than humans.

MATERIALS AND METHODS

We examined levels of genetic diversity in an available population genomic dataset of humans, chimpanzees, bonobos, gorillas, and orangutans for three sets of loci, (1) loci linked to facial morphology, (2) loci linked to height, and (3) neutrally evolving regions. We tested whether loci linked to facial morphology were more variable than loci linked to height or neutrally evolving loci in each of these species.

RESULTS

We found significantly elevated diversity in loci linked to facial morphology in chimpanzees, gorillas, and Sumatran and Bornean orangutans.

DISCUSSION

Our findings closely parallel those of Sheehan and Nachman and are consistent with the idea that selection for facial diversity and individual recognition has not only shaped the evolution of the human face, but it has similarly shaped the evolution of most of our closest primate relatives. We also discuss alternative hypotheses for this pattern.

Evaluating shared genetic influences on nonsyndromic cleft lip/palate and oropharyngeal neoplasms

It has been hypothesised that nonsyndromic cleft lip/palate (nsCL/P) and cancer may share aetiological risk factors. Population studies have found inconsistent evidence for increased incidence of cancer in nsCL/P cases, but several genes (e.g., CDH1, AXIN2) have been implicated in the aetiologies of both phenotypes. We aimed to evaluate shared genetic aetiology between nsCL/P and oral cavity/oropharyngeal cancers (OC/OPC), which affect similar anatomical regions. Using a primary sample of 5,048 OC/OPC cases and 5,450 controls of European ancestry and a replication sample of 750 cases and 336,319 controls from UK Biobank, we estimate genetic overlap using nsCL/P polygenic risk scores (PRS) with Mendelian randomization analyses performed to evaluate potential causal mechanisms. In the primary sample, we found strong evidence for an association between a nsCL/P PRS and increased odds of OC/OPC (per standard deviation increase in score, odds ratio [OR]: 1.09; 95% confidence interval [CI]: 1.04, 1.13; p = .000053). Although confidence intervals overlapped with the primary estimate, we did not find confirmatory evidence of an association between the PRS and OC/OPC in UK Biobank (OR 1.02; 95% CI: 0.95, 1.10; p = .55). Mendelian randomization analyses provided evidence that major nsCL/P risk variants are unlikely to influence OC/OPC. Our findings suggest possible shared genetic influences on nsCL/P and OC/OPC.

Association of Genetic Polymorphisms of GREM1 Gene with Susceptibility to Non-Syndromic Cleft Lip with or without Cleft Palate in an Iranian Population

Background: Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is common congenital birth anomaly with multifactorial etiology. The GREM1 gene has been proposed to play a role in oral clefts development.Objective: The aim of the present study was to evaluate the correlation between GREM1 polymorphisms and the risk of NSCL/P in an Iranian population.Methods: Genotyping of rs7162202, rs12915554, rs3743105, rs1129456, and rs10318 polymorphisms of GREM1 gene in 150 NSCL/P and 152 healthy subjects was determined by the PCR-RFLP or T-ARMS-PCR.Results: The findings showed that the rs12915554 variant significantly increased the risk of NSCL/P in heterozygous (OR = 4.20, 95%CI = 2.46-7.11, p < 0.0001, AC vs AA), and allele (OR = 3.17, 95%CI = 2.00-5.08, p < 0.0001, C vs A) genetic models. The rs3743105 polymorphism was correlated with reduced risk of NSCL/P in heterozygous (OR = 0.49, 95%CI = 0.29-0.83, p = 0.008, AG vs GG) and dominant (OR = 0.54, 95%CI = 0.33-0.89, p = 0.018, GA + AA vs GG) genetic models. The rs1129456 variant was positively associated with the risk of NSCL/P in heterozygous (OR = 2.91, 95%CI = 1.12-7.38, p = 0.028, AT vs AA) and allele (OR = 2.80, 95%CI = 2.80-6.95, p = 0.031, T vs C). The rs10318 polymorphism significantly reduced NSCL/P risk in homozygous (OR = 0.20, 95%CI = 0.06-0.67, p = 0.013, TT vs CC), dominant (OR = 0.57, 95%CI = 0.36-0.91, p = 0.019, CT + CC vs CC), recessive (OR = 0.24, 95%CI = 0.07-0.76, p = 0.031, TT vs CT + CC), and allele (OR = 0.57, 95%CI = 0.38-0.84, p = 0.005, T vs C). No correlation was observed between rs7162202 polymorphism and NSCL/P.Conclusion: The findings support that GREM1 polymorphisms are involved in NSCL/P susceptibility in an Iranian population.

SNPs at TP63 gene was specifically associated with right-side cleft lip in Han Chinese population

OBJECTIVES

Non-syndromic cleft lip with or without palate is one of the most common birth malformations. TP63 and GREM1 were recently reported to be associated with NSCL/P. However, there were few studies focused on their associations in non-syndromic cleft lip only (NSCLO).

DESIGN

Initial screening and replication in large cohorts were used to locate the susceptible SNPs of TP63 and GREM1. Firstly, variations were screened among 192 NSCLO cases by the Sanger sequencing. Then, we selected five associated SNPs in initial screening phase and replicated among 1,006 NSCLO cases and 1,823 normal controls.

RESULTS

Initial chi-square test showed that rs7653848, rs7624324, rs6790167, and rs1345186 in TP63 and rs2280738 in GREM1 achieved statistical significance (p < .05); the subsequent replication analysis showed that rs1345186 was specifically significant in right-side cleft lip (RCL; p = .017, OR = 1.33, and 95% CI: 1.05-1.69).

CONCLUSION

This study firstly used the subphenotype of cleft lip samples to verify the association between TP63 and GREM1, which indicated that TP63 is a promising susceptible gene for RCL in Chinese population. And further confirmed the different etiology in the right-sided cleft lip, left-sided cleft lip, and bilateral cleft lip of NSCLO. This will give new reference for the future research and genetic counseling.

Genetics and signaling mechanisms of orofacial clefts

Craniofacial development involves several complex tissue movements including several fusion processes to form the frontonasal and maxillary structures, including the upper lip and palate. Each of these movements are controlled by many different factors that are tightly regulated by several integral morphogenetic signaling pathways. Subject to both genetic and environmental influences, interruption at nearly any stage can disrupt lip, nasal, or palate fusion and result in a cleft. Here, we discuss many of the genetic risk factors that may contribute to the presentation of orofacial clefts in patients, and several of the key signaling pathways and underlying cellular mechanisms that control lip and palate formation, as identified primarily through investigating equivalent processes in animal models, are examined.

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.

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

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

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

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

Likely Pathogenic Variants in One Third of Non-Syndromic Discontinuous Cleft Lip and Palate Patients

Oral clefts are composed of cleft of the lip, cleft of the lip and palate, or cleft of the palate, and they are associated with a wide range of expression and severity. When cleft of the palate is associated with cleft of the lip with preservation of the primary palate, it defines an atypical phenotype called discontinuous cleft. Although this phenotype may represent 5% of clefts of the lip and/or palate (CLP), it is rarely specifically referred to and its pathophysiology is unknown. We conducted whole exome sequencing (WES) and apply a candidate gene approach to non-syndromic discontinuous CLP individuals in order to identify genes and deleterious variants that could underlie this phenotype. We discovered loss-of-function variants in two out of the seven individuals, implicating FGFR1 and DLG1 genes, which represents almost one third of this cohort. Whole exome sequencing of clinically well-defined subgroups of CLP, such as discontinuous cleft, is a relevant approach to study CLP etiopathogenesis. It could facilitate more accurate clinical, epidemiological and fundamental research, ultimately resulting in better diagnosis and care of CLP patients. Non-syndromic discontinuous cleft lip and palate seems to have a strong genetic basis.

High-Resolution Epigenomic Atlas of Human Embryonic Craniofacial Development

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

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

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

PAX7 nucleotide variants and the risk of non-syndromic orofacial clefts in the Polish population

OBJECTIVE

The etiology of non-syndromic cleft lip with or without cleft palate (nsCL/P) is multifactorial, heterogeneous, and still not completely understood. The aim of the present study was to examine the associations between common and rare PAX7 nucleotide variants and the risk of this common congenital anomaly in a Polish population.

SUBJECTS AND METHODS

Eight top nsCL/P-associated PAX7 variants identified in our cleft genome-wide association study (GWAS) were selected for replication analysis in an independent group of patients and controls (n = 247 and n = 445, respectively). In addition, mutation screening of the PAX7 protein-coding region was conducted.

RESULTS

Analysis of the pooled data from the GWAS and replication study confirmed that common PAX7 nucleotide variants are significantly associated with the increased risk of nsCL/P. The strongest individual variant was rs1339062 (c.586 + 15617T > C) with a p-value = 2.47E-05 (OR = 1.4, 95%CI: 1.20-1.64). Sequencing analysis identified a novel synonymous PAX7 substitution (c.87G > A, p.Val29Val) in a single patient with nsCLP. This transition located in the early exonic position was predicted to disrupt potential splice enhancer elements.

CONCLUSION

Our study confirmed that PAX7 is a strong candidate gene for nsCL/P. Nucleotide variants of this gene contribute to the etiology of nsCL/P in the homogenous Polish population.

Two novel genes TOX3 and COL21A1 in large extended Malay families with nonsyndromic cleft lip and/or palate

BACKGROUND

Nonsyndromic cleft lip and/or palate is one of the most common human birth defects worldwide that affects the lip and/or palate. The incidence of clefts varies among populations through ethnic, race, or geographical differences. The focus on Malay nonsyndromic cleft lip and/or palate (NSCL/P) is because of a scarce report on genetic study in relation to this deformity in Malaysia. We are interested to discuss about the genes that are susceptible to cause orofacial cleft formation in the family.

METHODS

Genome-wide linkage analysis was carried out on eight large extended families of NSCL/P with the total of 91 individuals among Malay population using microarray platform. Based on linkage analyses findings, copy number variation (CNV) of LPHN2, SATB2, PVRL3, COL21A1, and TOX3 were identified in four large extended families that showed linkage evidence using quantitative polymerase chain reaction (qPCR) as for a validation purpose. Copy number calculated (CNC) for each genes were determined with Applied Biosystems CopyCallerTM Software v2.0. Normal CNC of the target sequence expected was set at two.

RESULTS

Genome-wide linkage analysis had discovered several genes including TOX3 and COL21A1 in four different loci 4p15.2-p16.1, 6p11.2-p12.3, 14q13-q21, and 16q12.1. There was significant decreased, p < 0.05 of SATB2, COL21A1, and TOX3 copy number in extended families compared to the normal controls.

CONCLUSION

Novel linkage evidence and significant low copy number of COL21A1 and TOX3 in NSCLP family was confirmed. These genes increased the risks toward NSCLP formation in that family traits.

Genomic analyses in African populations identify novel risk loci for cleft palate

Orofacial clefts are common developmental disorders that pose significant clinical, economical and psychological problems. We conducted genome-wide association analyses for cleft palate only (CPO) and cleft lip with or without palate (CL/P) with ~17 million markers in sub-Saharan Africans. After replication and combined analyses, we identified novel loci for CPO at or near genome-wide significance on chromosomes 2 (near CTNNA2) and 19 (near SULT2A1). In situ hybridization of Sult2a1 in mice showed expression of SULT2A1 in mesenchymal cells in palate, palatal rugae and palatal epithelium in the fused palate. The previously reported 8q24 was the most significant locus for CL/P in our study, and we replicated several previously reported loci including PAX7 and VAX1.

Understanding the participation of GREM1 polymorphisms in nonsyndromic cleft lip with or without cleft palate in the Brazilian population

BACKGROUND

GREM1, which encodes Gremlin 1, an antagonist of bone morphogenic proteins with effects on proliferation and apoptosis, has been considered a candidate gene for nonsyndromic cleft lip with or without cleft palate (NSCL±P). In this study, we investigated potential associations of single nucleotide polymorphisms (SNP) in GREM1 and NSCL±P risk in the Brazilian population. Additionally, SNP-SNP interactions of GREM1 with previously reported rs1880646 variant in NTN1 (netrin 1), a gene also responsible for apoptotic phenotypes were verified.

METHODS

Applying Taqman allelic discrimination assays, we evaluated the variants rs16969681, rs16969816, rs16969862, and rs1258763 in 325 case-parent trios and in 1,588 isolated samples in a case-control study. Allelic and genotypic analyses, as well as interaction tests assessing gene-environmental factor (GxE) and SNP-SNP interaction with rs1880646 variant in NTN1, were performed based on logistic regression analysis adjusted for the effects of gender and genomic ancestry proportions.

RESULTS

The risk alleles of all SNP were undertransmitted in NSCL±P trios, though the case-control analysis confirmed only the association with rs16969862 alleles (OR: 0.78, 95% CI: 0.63-0.96, p = .02). The GxE interaction analysis revealed a significant interaction between maternal environmental contact with agrotoxics and rs16969816 (OR: 0.25, 95% CI: 0.08-0.74, p = .01), and pairwise interaction test with NTN1 rs1880646 yielded significant p values in the 1,000 permutation test for rs16969681, rs16969816, and rs16969862.

CONCLUSION

The GREM1 is involved in the etiology of NSCL±P in the Brazilian population and reveal that the interaction between GREM1 and NTN1 may be related with the pathogenesis of this common craniofacial malformation.

Birth Prevalence of Orofacial Clefts in Kuwait From Hospital-Based Registration

Introduction:

Cleft lip and palate (CLAP) are the most common craniofacial anomalies and birth defects globally. Despite the fact that a tertiary care registry of clefts has existed in Kuwait since 2008, to date there is no published data regarding the prevalence of orofacial clefts in this population.

Objective:

To tabulate the pattern of orofacial clefts from tertiary care center registration during 2009 through 2014 and to estimate the prevalence and trend using population-based records.

Methodology:

Data from all CLAP cases (born in Kuwait) registered in the central cleft center registry of the Al-Amiri hospital, Kuwait City, Kuwait, from January 2009 to December 2014 were obtained. Data regarding the type, severity, gender as well as nationality, parental consanguinity, and associated syndrome were obtained from medical records. Birth prevalence was tabulated against the population statistics for the period obtained from the central department of statistics.

Result:

A total of 202 CLAP patients were recorded in the study period with a mean birth prevalence of 0.57 per 1000 live births (95% confidence interval [CI] .57 ± .23). The registry recorded 108 (53.2%) males and 94 (47.8%) females. Children born to Kuwaitis represented 53.7% of cases while those born to non-Kuwaitis represented 45.3%. The most common oral cleft was CLAP (47.3%), followed by cleft palate (30.5%), cleft lip (20.2%), and other facial clefts (2%). Other congenital anomalies were recognized in 33% of all cases. There were no statistically significant differences in oral cleft prevalence across gender or nationality.

Conclusion:

The prevalence of oral cleft in Kuwait appears to be similar to those of other Middle Eastern populations.

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

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

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

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

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

Gene-Gene Interactions among SPRYs for Nonsyndromic Cleft Lip/Palate

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is a common birth defect with a complex genetic architecture. Gene-gene interactions have been increasingly regarded as contributing to the etiology of NSCL/P. A recent genome-wide association study revealed that a novel single-nucleotide polymorphism at SPRY1 in 4q28.1 showed a significant association with NSCL/P. In the current study, we explored the role of 3 SPRY genes in the etiology of NSCL/P by detecting gene-gene interactions: SPRY1, SPRY2, and SPRY4-with SPRY3 excluded due to its special location on the X chromosome. We selected markers in 3 SPRY genes to test for gene-gene interactions using 1,908 case-parent trios recruited from an international consortium established for a genome-wide association study of nonsyndromic oral clefts. As the trios came from populations with different ancestries, subgroup analyses were conducted among Europeans and Asians. Cordell's method based on conditional logistic regression models was applied to test for potential gene-gene interactions via the statistical package TRIO in R software. Gene-gene interaction analyses yielded 10 pairs of SNPs in Europeans and 6 pairs in Asians that achieved significance after Bonferroni correction. The significant interactions were confirmed in the 10,000-permutation tests (empirical P = 0.003 for the most significant interaction). The study identified gene-gene interactions among SPRY genes among 1,908 NSCL/P trios, which revealed the importance of potential gene-gene interactions for understanding the genetic architecture of NSCL/P. The evidence of gene-gene interactions in this study also provided clues for future biological studies to further investigate the mechanism of how SPRY genes participate in the development of NSCL/P.

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

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

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.

Application of high-resolution array platform for genome-wide copy number variation analysis in patients with nonsyndromic cleft lip and palate

BACKGROUND

Although more than 14 loci may be involved in the development of nonsyndromic cleft lip and palate (NSCLP), the etiology has not been fully elucidated due to genetic and environmental risk factor interactions. Despite advances in identifying genes associated with the NSCLP development using traditional genetic mapping strategies of candidate genes, genome-wide studies, and epidemiologic and linkage analysis, microarray techniques have become important complementary tools in the search for potential causative oral clefts genes in genetic studies. Microarray hybridization enables scanning of the whole genome and detecting copy number variants (CNVs). Although common benign CNVs are often smaller, with sizes smaller than 20 kb, here we reveal small exonic CNVs based on the importance of the encompassed genes in cleft lip and palate phenotype.

METHODS

Microarray hybridization analysis was performed in 15 individuals with NSCLP.

RESULTS

We identified 11 exonic CNVs affecting at least one exon of the candidate genes. Thirteen candidate genes (COL11A1-1p21; IRF6-1q32.3; MSX1-4p16.2; TERT-5p15.33; MIR4457-5p15.33; CLPTM1L-5p15.33; ESR1-6q25.1; GLI3-7p13; FGFR-8p11.23; TBX1-22q11.21; OFD-Xp22; PHF8-Xp11.22; and FLNA-Xq28) overlapped with the CNVs identified.

CONCLUSIONS

Considering the importance to NSCLP, the microdeletions that encompass MSX1, microduplications over TERT, MIR4457, CLPTM1L, and microduplication of PHF8 have been identified as small CNVs related to sequence variants associated with oral clefts susceptibility. Our findings represent a preliminary study on the clinical significance of small CNVs and their relationship with genes implicated in NSCLP.

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.

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

OBJECTIVE

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

PATIENTS AND METHOD

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

RESULTS

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

CONCLUSION

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

Association between a single-nucleotide polymorphism in the GREM1 gene and non-syndromic orofacial cleft in the Chinese population

BACKGROUND

Non-syndromic orofacial cleft (NSOC) is a common craniofacial deformity among newborns. The GREM1 gene is correlated with orofacial development. The aim of our study was to investigate the association between a single-nucleotide polymorphism in the GREM1 gene and this malformation in the Chinese population.

METHODS

The SNaPshot mini-sequencing technique was used to genotype the locus rs1258763 of the GREM1 gene in 331 patients with NSOC and 271 individuals in a control group.

RESULTS

For GREM1 rs1258763, there was a significant difference between the NSOC case group and control group (P = .022). Children carrying GA and GA/AA genotypes had an increased risk of NSOC (OR=1.62, 95%CI: 1.15-2.30; OR=1.52, 95%CI: 1.09-2.12). In the cleft subgroup, we found that the GREM1 rs1258763 GA genotype might contribute to the elevated risk of the cleft lip with or without cleft palate (CL/P) (P = .029). Non-significant differences were found between the cleft palate only (CPO) and control groups (P = .077).

CONCLUSION

Our findings revealed that the GREM1 polymorphism was significantly associated with the risk of NSOC in the Chinese population.

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.