Genome-Wide Association Study of Non-syndromic Orofacial Clefts in a Multiethnic Sample of Families and Controls Identifies Novel Regions

A non-syndromic orofacial cleft risk locus links tRNA splicing defects to neural crest cell pathologies

Orofacial clefts are the most common form of congenital craniofacial malformation worldwide. The etiology of these birth defects is multifactorial, involving genetic and environmental factors. However, in most cases, the underlying causes remain unexplained, precluding a molecular understanding of disease mechanisms. Here, we integrated genome-wide association data, targeted resequencing of case and control cohorts, tissue- and cell-type-specific epigenomic profiling, and genome architecture analyses to molecularly dissect a genomic locus associated with an increased risk of non-syndromic orofacial cleft. We found that common and rare risk variants associated with orofacial cleft intersect with an enhancer (e2p24.2) that is active in human embryonic craniofacial tissue. We mapped e2p24.2 long-range interactions to a topologically associated domain harboring MYCN, DDX1, and CYRIA. We found that MYCN and DDX1, but not CYRIA, are required during craniofacial development in chicken embryos. We investigated the role of DDX1, a key component of the tRNA splicing complex, in cranial neural crest cells (cNCCs). The loss of DDX1 in cNCCs resulted in the accumulation of unspliced tRNA fragments, depletion of mature intron-containing tRNAs, and ribosome stalling at codons decoded by these tRNAs. This was accompanied by defects in both global protein synthesis and cNCC migration. We further showed that the induction of tRNA fragments is sufficient to disrupt craniofacial development. Together, these results uncovered a molecular mechanism in which impaired tRNA splicing affects cNCCs and craniofacial development and positioned MYCN, DDX1, and tRNA processing defects as risk factors in the pathogenesis of orofacial clefts.

A multi-ancestry GWAS meta-analysis of facial features and its application in predicting archaic human features

Facial morphology, a complex trait influenced by genetics, holds great significance in evolutionary research. However, due to limited fossil evidence, the facial characteristics of Neanderthals and Denisovans have remained largely unknown. In this study, we conduct a large-scale multi-ethnic meta-analysis of the genome-wide association study (GWAS), including 9674 East Asians and 10,115 Europeans, quantitatively assessing 78 facial traits using 3D facial images. We identify 71 genomic loci associated with facial features, including 21 novel loci. We develop a facial polygenic score (FPS) that enables the prediction of facial features based on genetic information. Interestingly, the distribution of FPSs among populations from diverse continental groups exhibits relevant correlations with observed facial features. Furthermore, we apply the FPS to predict the facial traits of seven Neanderthals and one Denisovan using ancient DNA and align predictions with the fossil records. Our results suggest that Neanderthals and Denisovans likely share similar facial features, such as a wider but shorter nose and a wider endocanthion distance. The decreased mouth width is characterized specifically in Denisovans. The integration of genomic data and facial trait analysis provides valuable insights into the evolutionary history and adaptive changes in human facial morphology.

A Novel De Novo Missense Variant in Netrin-1 (NTN1) Associated With Chorioretinal Coloboma, Sensorineural Hearing Loss and Polydactyly

Microphthalmia, anophthalmia and coloboma (MAC) comprise a highly heterogeneous spectrum of congenital ocular malformations with an estimated incidence of 1 in 5000 to 1 in 30 000 live births. Although there is likely to be a genetic component in the majority of cases, many remain without a molecular diagnosis. Netrin-1 was previously identified as a mediator of optic fissure closure from transcriptome analyses of chick and zebrafish and was shown to cause ocular coloboma when knocked out in both mouse and zebrafish. Here, we report the first patient with chorioretinal coloboma and microphthalmia harbouring a novel heterozygous likely pathogenic NTN1 missense variant, c.1483T>A p.(Tyr495Asn), validating a conserved gene function in ocular development. In addition, the patient displayed bilateral sensorineural hearing loss which was investigated by examining the sensory hair cells of ntn1a morphant zebrafish, suggesting a role for netrin-1 in hair cell development.

Genetic-epigenetic interactions (meQTLs) in orofacial clefts etiology

Objectives

Nonsyndromic orofacial clefts (OFCs) etiology involves multiple genetic and environmental factors with over 60 identified risk loci; however, they account for only a minority of the estimated risk. Epigenetic factors such as differential DNA methylation (DNAm) are also associated with OFCs risk and can alter risk for different cleft types and modify OFCs penetrance. DNAm is a covalent addition of a methyl (CH3) group to the nucleotide cytosine that can lead to changes in expression of the targeted gene. DNAm can be affected by environmental influences and genetic variation via methylation quantitative loci (meQTLs). We hypothesize that aberrant DNAm and the resulting alterations in gene expression play a key role in the etiology of OFCs, and that certain common genetic variants that affect OFCs risk do so by influencing DNAm.

Methods

We used genotype from 10 cleft-associated SNPs and genome-wide DNA methylation data (Illumina 450K array) for 409 cases with OFCs and 456 controls and identified 23 cleft-associated meQTLs. We then used an independent cohort of 362 cleft-discordant sib pairs for replication. We used methylation-specific qPCR to measure methylation levels of each CpG site and combined genotypic and methylation data for an interaction analysis of each SNP-CpG pair using the R package MatrixeQTL in a linear model. We also performed a Paired T-test to analyze differences in DNA methylation between each member of the sibling pairs.

Results

We replicated 9 meQTLs, showing interactions between rs13041247 (MAFB) - cg18347630 (PLCG1) (P=0.04); rs227731 (NOG) - cg08592707 (PPM1E) (P=0.01); rs227731 (NOG) - cg10303698 (CUEDC1) (P=0.001); rs3758249 (FOXE1) - cg20308679 (FRZB) (P=0.04); rs8001641 (SPRY2) - cg19191560 (LGR4) (P=0.04); rs987525(8q24) - cg16561172(MYC) (P=0.00000963); rs7590268(THADA) - cg06873343 (TTYH3) (P=0.04); rs7078160 (VAX1) - cg09487139 (P=0.05); rs560426 (ABCA4/ARHGAP29) - cg25196715 (ABCA4/ARHGAP29) (P=0,03). Paired T-test showed significant differences for cg06873343 (TTYH3) (P=0.04); cg17103269 (LPIN3) (P=0.002), and cg19191560 (LGR4) (P=0.05).

Conclusions

Our results confirm previous evidence that some of the common non-coding variants detected through GWAS studies can influence the risk of OFCs via epigenetic mechanisms, such as DNAm, which can ultimately affect and regulate gene expression. Given the large prevalence of non-coding SNPs in most OFCs genome wide association studies, our findings can potentially address major knowledge gaps, like missing heritability, reduced penetrance, and variable expressivity associated with OFCs phenotypes.

Association of selected gene variants with nonsyndromic orofacial clefts in Kuwait

INTRODUCTION AND OBJECTIVES

Non-syndromic orofacial clefts (NSOFCs) are complex congenital abnormalities involving both environmental and genetic factors involved in orofacial development. This study aimed to investigate the genetic association of specific genetic variants at different CYRIA gene loci with the development of NSOFCs in Kuwait.

METHODS

Four genetic variants (rs7552, rs3758249, rs3821949, and rs3917201) at four selected gene loci (CYRIA, FOXE1, MSX1, and TGFB3) were genotyped in a total of 240 DNA samples (patients (n = 114) and random controls (n = 126)) employing TaqMan® allele discrimination assay. For each variant and its genotype, the frequencies were determined and tested for Hardy-Weinberg Equilibrium. Genotype frequencies was compared between patients and controls using Pearson's test. Logistic regression analyses were employed to test for the associations of the four selected variants with the occurrence of NSOFCSs.

RESULTS

Significant differences in the distribution of genotypes between cases and controls, rs7552, rs3821949, and rs3917201 were found to have a positive association with NSOFCs. After adjusting for gender, the GG genotype of the rs7552 variant, the AG genotype of the rs3821949 variant, and the CC genotype of the rs3917201 variant showed nearly a two-fold increased risk of NSOFC (p < 0.05).

CONCLUSION

This study reports significant findings on the contribution and modest effect of CYRIA rs7552, MSX1 rs3821949, and TGFB3 rs3917201 in the development of NSOFCs. Our findings provide further evidence on the molecular mechanism and the role of the selected genes in NSOFCs.

Multi-ancestry Genome Wide Association Study Meta-analysis of Non-syndromic Orofacial Clefts

Non-syndromic orofacial clefts (NSOC) are common craniofacial birth defects, and result from both genetic and environmental factors. NSOC include three major sub-phenotypes: non-syndromic cleft lip with palate (NSCLP), non-syndromic cleft lip only (NSCLO) and non-syndromic cleft palate only (NSCPO), NSCLP and NSCLO are also sometimes grouped as non-syndromic cleft lip with or without cleft palate (NSCL/P) based on epidemiology. Currently known loci only explain a limited proportion of the heritability of NSOC. Further, differences in genetic susceptibility among the sub-phenotypes are poorly characterized. We performed a multi-ancestry GWAS meta-analysis on 44,094 individuals (9,381 cases, 28,510 controls, 2042 case-parent trios and 18 multiplex pedigrees) of East Asian, European, Latin and South American, and African ancestry for both NSOC and subtypes. We identified 50 loci, including 11 novel loci: four loci ( CALD1 , SHH , NRG1 and LINC00320 ) associated with both NSOC and NSCL/P, two loci ( NTRK1 and RUNX1 ) only associated with NSOC, four loci ( HMGCR , PRICKLE1 , SOX9 and MYH9 ) only associated with NSCL/P and one locus ( ALX1 ) specifically associated with NSCLO. Five of the novel loci are located in regions containing genes associated with syndromic orofacial clefts ( SHH , NTRK1, CALD1, ALX1 and SOX9 ); seven of the novel loci are located in regions containing genes-implicated in craniofacial development ( HMGCR, SHH, PRICKLE1, ALX1, SOX9, RUNX1, MYH9 ). Genetic correlation and colocalization analyses revealed an overlap between signals associated with NSCLO, NSCPO and NSCLP, but there were also notable differences, emphasizing the complexity of common and distinct genetic processes affecting lip and palate development.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSION

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

Functional analysis of ESRP1/2 gene variants and CTNND1 isoforms in orofacial cleft pathogenesis

Orofacial cleft (OFC) is a common human congenital anomaly. Epithelial-specific RNA splicing regulators ESRP1 and ESRP2 regulate craniofacial morphogenesis and their disruption result in OFC in zebrafish, mouse and humans. Using esrp1/2 mutant zebrafish and murine Py2T cell line models, we functionally tested the pathogenicity of human ESRP1/2 gene variants. We found that many variants predicted by in silico methods to be pathogenic were functionally benign. Esrp1 also regulates the alternative splicing of Ctnnd1 and these genes are co-expressed in the embryonic and oral epithelium. In fact, over-expression of ctnnd1 is sufficient to rescue morphogenesis of epithelial-derived structures in esrp1/2 zebrafish mutants. Additionally, we identified 13 CTNND1 variants from genome sequencing of OFC cohorts, confirming CTNND1 as a key gene in human OFC. This work highlights the importance of functional assessment of human gene variants and demonstrates the critical requirement of Esrp-Ctnnd1 acting in the embryonic epithelium to regulate palatogenesis.

Genetic associations and parent-of-origin effects of PVRL1 in non-syndromic cleft lip with or without cleft palate across multiple ethnic populations

OBJECTIVES

This study investigated the associations of PVRL1 gene variants with non-syndromic cleft lip with or without cleft palate (NSCL/P) by evaluating transmission distortion and parent-of-origin (POO) effects in multiple ethnic populations.

METHODS

We conducted allelic and genotypic transmission disequilibrium tests (TDT) on 10 single-nucleotide variants (SNVs) in PVRL1 using data from 142 Korean families with an affected child. POO effects were analyzed using the POO likelihood ratio test, comparing transmission rates of maternally and paternally inherited alleles. To assess generalizability and ethnic heterogeneity, we compared results from Korean families with data from the Center for Craniofacial and Dental Genetics, which included 2,226 individuals from 497 European and 245 Asian trios.

RESULTS

TDT analysis identified significant over-transmission of the rs7940667 (G361V) C allele in Korean families (p=0.007), a finding replicated in both Asian (p=6.5×10-7) and European families (p=1.6×10-10). Eight SNVs showed strong TDT evidence in larger Asian and European datasets after multiple comparison corrections (p<0.0073). Of these, 4 SNVs (rs7940667, rs7103685, rs7129848, and rs4409845) showed particularly robust association (p<5×10-8). POO analysis revealed significant maternal over-transmission of the rs10790330-A allele in Korean families (p=0.044). This finding was replicated in European families (p=9.0×10-4). Additionally, 3 other SNVs, rs7129848 (p=0.001) and the linked SNVs rs3935406 and rs10892434 (p=0.025), exhibited maternal over-transmission in the validation datasets.

CONCLUSIONS

Our findings provide robust evidence supporting the associations of PVRL1 variants with NSCL/P susceptibility. Further research is necessary to explore the potential clinical applications of these findings.

DeepFace: Deep-learning-based framework to contextualize orofacial-cleft-related variants during human embryonic craniofacial development

Orofacial clefts (OFCs) are among the most common human congenital birth defects. Previous multiethnic studies have identified dozens of associated loci for both cleft lip with or without cleft palate (CL/P) and cleft palate alone (CP). Although several nearby genes have been highlighted, the "casual" variants are largely unknown. Here, we developed DeepFace, a convolutional neural network model, to assess the functional impact of variants by SNP activity difference (SAD) scores. The DeepFace model is trained with 204 epigenomic assays from crucial human embryonic craniofacial developmental stages of post-conception week (pcw) 4 to pcw 10. The Pearson correlation coefficient between the predicted and actual values for 12 epigenetic features achieved a median range of 0.50-0.83. Specifically, our model revealed that SNPs significantly associated with OFCs tended to exhibit higher SAD scores across various variant categories compared to less related groups, indicating a context-specific impact of OFC-related SNPs. Notably, we identified six SNPs with a significant linear relationship to SAD scores throughout developmental progression, suggesting that these SNPs could play a temporal regulatory role. Furthermore, our cell-type specificity analysis pinpointed the trophoblast cell as having the highest enrichment of risk signals associated with OFCs. Overall, DeepFace can harness distal regulatory signals from extensive epigenomic assays, offering new perspectives for prioritizing OFC variants using contextualized functional genomic features. We expect DeepFace to be instrumental in accessing and predicting the regulatory roles of variants associated with OFCs, and the model can be extended to study other complex diseases or traits.

Functional analysis of ESRP1/2 gene variants and CTNND1 isoforms in orofacial cleft pathogenesis

Orofacial cleft (OFC) is a common human congenital anomaly. Epithelial-specific RNA splicing regulators ESRP1 and ESRP2 regulate craniofacial morphogenesis and their disruption result in OFC in zebrafish, mouse and humans. Using esrp1/2 mutant zebrafish and murine Py2T cell line models, we functionally tested the pathogenicity of human ESRP1/2 gene variants. We found that many variants predicted by in silico methods to be pathogenic were functionally benign. Esrp1 also regulates the alternative splicing of Ctnnd1 and these genes are co-expressed in the embryonic and oral epithelium. In fact, over-expression of ctnnd1 is sufficient to rescue morphogenesis of epithelial-derived structures in esrp1/2 zebrafish mutants. Additionally, we identified 13 CTNND1 variants from genome sequencing of OFC cohorts, confirming CTNND1 as a key gene in human OFC. This work highlights the importance of functional assessment of human gene variants and demonstrates the critical requirement of Esrp-Ctnnd1 acting in the embryonic epithelium to regulate palatogenesis.

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.

Updates in Cleft Care

Cleft lip and/or palate is a congenital malformation with a wide range of presentations, and its effective treatment necessitates sustained, comprehensive care across an affected child's life. Early diagnosis, ideally through prenatal imaging or immediately postbirth, is paramount. Access to longitudinal care and long-term follow-up with a multidisciplinary approach, led by the recommendations of the American Cleft Palate Association, is the best way to ensure optimal outcomes. Multiple specialties including plastic surgery, otolaryngology, speech therapy, orthodontists, psychologists, and audiologists all may be indicated in the care of the child. Primary repair of the lip, nose, and palate are generally conducted during infancy. Postoperative care demands meticulous oversight to detect potential complications. If necessary, revisional surgeries should be performed before the child begin primary school. As the child matures, secondary procedures like alveolar bone grafting and orthognathic surgery may be requisite. The landscape of cleft care has undergone significant transformation since early surgical correction, with treatment plans now tailored to the specific type and severity of the cleft. The purpose of this text is to outline the current standards of care in children born with cleft lip and/or palate and to highlight ongoing advancements in the field.

A Cross-Sectional Study of the Dental Arch Relationship and Palatal Morphology after Cleft Surgery in Italian Children with Unilateral Cleft and Lip Palate

Dental arch relationships (DARs) and palatal morphology (PM) were evaluated in in non-syndromic complete unilateral cleft lip and palate (UCLP) Italian patients after surgery. Pre- and postnatal factors affecting the results were investigated. Sixty-six children with UCLP (40 boys and 26 girls, with a mean age of 10.1 ± 2.9 years), predominantly Caucasian (77%), were consecutively enrolled in this cross-sectional study. Twenty children had received a one-stage protocol consisting of an early periosteal palate surgical repair and lip closure and forty-six were submitted to a staged surgical protocol with delayed palate repair (DPR). A single clinician collected data on their medical history and carried out a dental examination. The DAR and PM were graded on dental casts according to the Eurocran index and dichotomised as favourable and unfavourable based on the treatment outcome. Multiple logistic regression analyses demonstrated that female sex (OR = 6.08, 95% CI: 1.47-25.23, p = 0.013), DPR (OR = 4.77, 95% CI: 1.14-19.93, p = 0.032) and the use of a neonatal plate (OR = 4.68, 95% CI: 1.27-17.16, p = 0.020) increased the odds of having favourable DAR, while only DPR (OR = 9.76, 95% CI: 2.40-39.71, p = 0.001) was significantly associated with a favourable PM. Based on these findings, only DPR had a significantly favourable effect on both DAR and DM in Italian children with complete UCLP.

High throughput miRNA sequencing and bioinformatics analysis identify the mesenchymal cell proliferation and apoptosis related miRNAs during fetal mice palate development

BACKGROUND

Palatogenesis requires a precise spatiotemporal regulation of gene expression. Recent studies indicate that microRNAs (miRNAs) are key factors in normal palatogenesis. The present study aimed to explain the regulatory mechanisms of miRNAs during palate development.

METHODS

Pregnant ICR mice were choose at embryonic day 10.5 (E10.5). Hemotoxylin and eosin (H&E) staining was used to observe the morphological changes during the development of palatal process at embryonic day (E)13.5, E14.0, E14.5, E15.0 and E15.5. The fetal palatal tissues were collected at E13.5, E14.0, E14.5 and E15.0 to explore miRNA expression and function by high throughput sequencing and bioinformatic analysis. Mfuzz cluster analysis was used to look for miRNAs related to the fetal mice palate formation. The target genes of miRNAs were predicted by miRWalk. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was performed base on target genes. The mesenchymal cell proliferation and apoptosis related miRNAs-genes networks were predicted and constructed using miRWalk and Cytoscape software. The expression of mesenchymal cell proliferation and apoptosis related miRNAs at the E13.5, E14.0, E14.5, and E15.0 was detected by a quantitative real-time PCR (RT-qPCR) assay.

RESULTS

H&E staining found that the palatal process grows vertically along the sides of the tongue at E13.5, the position of the tongue begins to descend and the bilateral palatal processes rise above the tongue at E14.0, the palatal process grows horizontally at E14.5, there is palatal contact fusion at E15.0, and the palatal suture disappeared at E15.5. Nine clusters of miRNA expression changes were identified in the fetal mice palate formation progression, including two reducing trends, two rising trends and five disordered trends. Next, the heatmap showed the miRNA expression from Clusters 4, 6, 9, 12 in the E13.5, E14.0, E14.5 and E15.0 groups. GO functional and KEGG pathway enrichment analysis found target genes of miRNAs in clusters involved in regulation of mesenchymal phenotype and the mitogen-activated protein kinase (MAPK) signaling pathway. Next, mesenchymal phenotype related miRNA-genes networks were constructed. The heatmap showing that the mesenchymal phenotype related miRNA expression of Clusters 4, 6, 9 and 12 at E13.5, E14.0, E14.5 and E15.0. Furthermore, the mesenchymal cell proliferation and apoptosis related miRNA-gene networks were identified in Clusters 6 and 12, including mmu-miR-504-3p-Hnf1b, etc. The expression level of mesenchymal cell proliferation and apoptosis related miRNAs at the E13.5, E14.0, E14.5, and E15.0 was verified by a RT-qPCR assay.

CONCLUSIONS

For the first time, we identified that clear dynamic miRNA expression during palate development. Furthermore, we demonstrated that mesenchymal cell proliferation and apoptosis related miRNAs, genes and the MAPK signaling pathway are important during fetal mice palate development.

Genetics and orofacial clefts: a clinical perspective

Orofacial clefts (OFCs) are the most common congenital craniofacial anomaly seen in humans. Most OFCs are sporadic and isolated - these are thought to be multifactorial in origin. Chromosomal and monogenic variants account for the syndromic forms and for some of the non-syndromic inherited forms. This review discusses the importance of genetic testing and the current clinical strategy to deliver a genomics service that is of direct benefit to patients and their families.

Identification of putative regulatory single-nucleotide variants in NTN1 gene associated with NSCL/P

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is a common polygenetic disease. Although genome-wide association studies (GWAS) identified NTN1 gene as a high-priority candidate of NSCL/P, the comprehensive genetic architecture of NTN1 weren't yet known. Thus, this study aimed to determine full-scale genetic variants of NTN1 for NSCL/P in Chinese Han people. Initially, targeted sequencing of NTN1 gene was performed on 159 NSCL/P patients to identify susceptible single nucleotide polymorphisms (SNPs) associated with NSCL/P. Then, association analysis and burden analysis were separately used to validate the common variants and rare variants identified among large size of samples (1608 NSCL/P cases and 2255 controls). Additionally, NSCL/P subtype association analysis was applied to elucidate the etiology discrepancy of non-syndromic cleft lip with palate (NSCLP) and non-syndromic cleft lip only (NSCLO). Lastly, bioinformatics analysis was performed to annotate and prioritize candidate variants. We found 15 NSCL/P-associated SNPs including rs4791774 (P = 1.10E-08, OR = 1.467, 95% CI: 1.286~1.673) and rs9788972 (P = 1.28E-07, OR = 1.398, 95% CI : 1.235~1.584) originally detected by previous GWASs in Chinese Han ancestry. Four NSCLO risk-associated SNPs and eight specific NSCLP associated SNPs were found. Three SNPs (rs4791331, rs4791774 and rs9900753) were predicted to locate at regulatory region of NTN1. Our study validated the association between NTN1 gene and pathogenesis of NSCL/P and reinforced the hypothesis that NSCLP have a different etiology from NSCLO. We also identified three putative regulatory SNPs in NTN1 gene.

The genetic determinants of oral diseases in Africa: The gaps should be filled

Oral diseases are a major health concern and are among the most prevalent diseases globally. This problem is becoming more prominent in the rapidly growing populations of Africa. It is well documented that Africa exhibits the most diverse genetic make-up in the world. However, little work has been conducted to understand the genetic basis of oral diseases in Africans. Oral health is often neglected and receives low prioritisation from funders and governments. The genetic determinants of highly prevalent oral diseases such as dental caries and periodontal disease, and regionally prevalent conditions such as oral cancer and NOMA, are largely under-researched areas despite numerous articles alluding to a high burden of these diseases in African populations. Therefore, this review aims to shed light on the significant gaps in research on the genetic and genomic aspects of oral diseases in African populations and highlights the urgent need for evidence-based dentistry, in tandem with the development of the dentist/scientist workforce.

Decoding the Human Face: Progress and Challenges 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.

Genome-wide association study of multiethnic nonsyndromic orofacial cleft families identifies novel loci specific to family and phenotypic subtypes

Nonsyndromic orofacial clefts (OFCs) are among the most common craniofacial birth defects worldwide, and known to exhibit phenotypic and genetic heterogeneity. Cleft lip plus cleft palate (CLP) and cleft lip only (CL) are commonly combined together as one phenotype (CL/P), separately from cleft palate alone. In comparison, our study analyzes CL and CLP separately. A sample of 2218 CL and CLP cases, 4537 unaffected relatives of cases, and 2673 pure controls with no family history of OFC were selected from the Pittsburgh Orofacial Cleft (Pitt-OFC) multiethnic study.genome-wide association studies were run for seven specific phenotypes created based on the cleft type(s) observed within these families, as well as the combined CL/P phenotype. Five novel genome-wide significant associations, 3q29 (rs62284390), 5p13.2 (rs609659), 7q22.1 (rs6465810), 19p13.3 (rs628271), and 20q13.33 (rs2427238), and nine associations (p ≤ 1.0E-05) within previously confirmed OFC loci-PAX7, IRF6, FAM49A, DCAF4L2, 8q24.21, ARID3B, NTN1, TANC2 and the WNT9B:WNT3 gene cluster-were observed. We also found that single nucleotide polymorphisms within a subset of the associated loci, both previously known and novel, differ substantially in terms of their effects across cleft- or family-specific phenotypes, indicating not only etiologic differences between CL and CLP, but also genetic heterogeneity within each of the two OFC subtypes.

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

OBJECTIVE

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

DESIGN

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

RESULTS

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

CONCLUSION

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

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.