Shared genetic risk between major orofacial cleft phenotypes in an African population

Nonsyndromic orofacial clefts (NSOFCs) represent a large proportion (70%-80%) of all OFCs. They can be broadly categorized into nonsyndromic cleft lip with or without cleft palate (NSCL/P) and nonsyndromic cleft palate only (NSCPO). Although NSCL/P and NSCPO are considered etiologically distinct, recent evidence suggests the presence of shared genetic risks. Thus, we investigated the genetic overlap between NSCL/P and NSCPO using African genome-wide association study (GWAS) data on NSOFCs. These data consist of 814 NSCL/P, 205 NSCPO cases, and 2159 unrelated controls. We generated common single-nucleotide variants (SNVs) association summary statistics separately for each phenotype (NSCL/P and NSCPO) under an additive genetic model. Subsequently, we employed the pleiotropic analysis under the composite null (PLACO) method to test for genetic overlap. Our analysis identified two loci with genome-wide significance (rs181737795 [p = 2.58E-08] and rs2221169 [p = 4.5E-08]) and one locus with marginal significance (rs187523265 [p = 5.22E-08]). Using mouse transcriptomics data and information from genetic phenotype databases, we identified MDN1, MAP3k7, KMT2A, ARCN1, and VADC2 as top candidate genes for the associated SNVs. These findings enhance our understanding of genetic variants associated with NSOFCs and identify potential candidate genes for further exploration.

Damaging Mutations in Contribute to Risk of Nonsyndromic Cleft Lip With or Without Cleft Palate

Novel or rare damaging mutations have been implicated in the developmental pathogenesis of nonsyndromic cleft lip with or without cleft palate (nsCL ± P). Thus, we investigated the human genome for high-impact mutations that could explain the risk of nsCL ± P in our cohorts.

We conducted next-generation sequencing (NGS) analysis of 130 nsCL ± P case-parent African trios to identify pathogenic variants that contribute to the risk of clefting. We replicated this analysis using whole-exome sequence data from a Brazilian nsCL ± P cohort. Computational analyses were then used to predict the mechanism by which these variants could result in increased risks for nsCL ± P.

We discovered damaging mutations within the AFDN gene, a cell adhesion molecule (CAMs) that was previously shown to contribute to cleft palate in mice. These mutations include p.Met1164Ile, p.Thr453Asn, p.Pro1638Ala, p.Arg669Gln, p.Ala1717Val, and p.Arg1596His. We also discovered a novel splicing p.Leu1588Leu mutation in this protein. Computational analysis suggests that these amino acid changes affect the interactions with other cleft-associated genes including nectins (PVRL1, PVRL2, PVRL3, and PVRL4) CDH1, CTNNA1, and CTNND1.

This is the first report on the contribution of AFDN to the risk for nsCL ± P in humans. AFDN encodes AFADIN, an important CAM that forms calcium-independent complexes with nectins 1 and 4 (encoded by the genes PVRL1 and PVRL4). This discovery shows the power of NGS analysis of multiethnic cleft samples in combination with a computational approach in the understanding of the pathogenesis of nsCL ± P.

Rare Variants Analyses Suggest Novel Cleft Genes in the African Population

Non-syndromic orofacial clefts (NSOFCs) are common birth defects with a complex etiology. While over 60 common risk loci have been identified, they explain only a small proportion of the heritability for NSOFC. Rare variants have been implicated in the missing heritability. Thus, our study aimed to identify genes enriched with nonsynonymous rare coding variants associated with NSOFCs. Our sample included 814 non-syndromic cleft lip with or without palate (NSCL/P), 205 non-syndromic cleft palate only (NSCPO), and 2150 unrelated control children from Nigeria, Ghana, and Ethiopia. We conducted a gene-based analysis separately for each phenotype using three rare-variants collapsing models: (1) protein-altering (PA), (2) missense variants only (MO); and (3) loss of function variants only (LOFO). Subsequently, we utilized relevant transcriptomics data to evaluate associated gene expression and examined their mutation constraint using the gnomeAD database. In total, 13 genes showed suggestive associations (p = E-04). Among them, eight genes (ABCB1, ALKBH8, CENPF, CSAD, EXPH5, PDZD8, SLC16A9, and TTC28) were consistently expressed in relevant mouse and human craniofacial tissues during the formation of the face, and three genes (ABCB1, TTC28, and PDZD8) showed statistically significant mutation constraint. These findings underscore the role of rare variants in identifying candidate genes for NSOFCs. Main documents (excluding the methods section) word count: 2145.

Perceptions and beliefs of community gatekeepers about genomic risk information in African cleft research

BACKGROUND

A fundamental ethical issue in African genomics research is how socio-cultural factors impact perspectives, acceptance, and utility of genomic information, especially in stigmatizing conditions like orofacial clefts (OFCs). Previous research has shown that gatekeepers (e.g., religious, political, family or community leaders) wield considerable influence on the decision-making capabilities of their members, including health issues. Thus, their perspectives can inform the design of engagement strategies and increase exposure to the benefits of genomics testing/research. This is especially important for Africans underrepresented in genomic research. Our study aims to investigate the perspectives of gatekeepers concerning genomic risk information (GRI) in the presence of OFCs in a sub-Saharan African cohort.

METHODS

Twenty-five focus group discussions (FGDs) consisting of 214 gatekeepers (religious, community, ethnic leaders, and traditional birth attendants) in Lagos, Nigeria, explored the opinions of participants on genomic risk information (GRI), OFC experience, and the possibility of involvement in collaborative decision-making in Lagos, Nigeria. Transcripts generated from audio recordings were coded and analyzed in NVivo using thematic analysis.

RESULTS

Three main themes-knowledge, beliefs, and willingness to act-emerged from exploring the perspective of gatekeepers about GRI in this group. We observed mixed opinions regarding the acceptance of GRI. Many participants believed their role is to guide and support members when they receive results; this is based on the level of trust their members have in them. However, participants felt they would need to be trained by medical experts to do this. Also, religious and cultural beliefs were crucial to determining participants' understanding of OFCs and the acceptance and utilization of GRI.

CONCLUSIONS

Incorporating cultural sensitivity into public engagement could help develop appropriate strategies to manage conflicting ideologies surrounding genomic information in African communities. This will allow for more widespread access to the advances in genomics research in underrepresented populations. We also recommend a synergistic relationship between community health specialists/scientists, and community leaders, including spiritual providers to better understand and utilize GRI.

Parents and Provider Perspectives on the Return of Genomic Findings for Cleft Families in Africa

BACKGROUND

Inadequate knowledge among health care providers (HCPs) and parents of affected children limits the understanding and utility of secondary genetic findings (SFs) in under-represented populations in genomics research. SFs arise from deep DNA sequencing done for research or diagnostic purposes and may burden patients and their families despite their potential health importance. This study aims to evaluate the perspective of both groups regarding SFs and their choices in the return of results from genetic testing in the context of orofacial clefts.

METHODS

Using an online survey, we evaluated the experiences of 252 HCPs and 197 parents across participating cleft clinics in Ghana and Nigeria toward the return of SFs across several domains.

RESULTS

Only 1.6% of the HCPs felt they had an expert understanding of when and how to incorporate genomic medicine into practice, while 50.0% agreed that all SFs should be returned to patients. About 95.4% of parents were willing to receive all the information from genetic testing (including SFs), while the majority cited physicians as their primary information source (64%).

CONCLUSIONS

Overall, parents and providers were aware that genetic testing could help in the clinical management of diseases. However, they cited a lack of knowledge about genomic medicine, uncertain clinical utility, and lack of available learning resources as barriers. The knowledge gained from this study will assist with developing guidelines and policies to guide providers on the return of SFs in sub-Saharan Africa and across the continent.

Whole Exome Sequencing Identifies Damaging Variants in Indonesians with Clefts

OBJECTIVES

The interaction between genomics, genetic and environmental factors have been implicated in non-syndromic orofacial cleft development. In the current study, we investigated the contributions of rare and novel genetic variants in known cleft genes using whole exome sequencing (WES) data of Indonesians with non-syndromic orofacial clefts.

DESIGN

WES was conducted on 6 individuals. Variants in their exons were called and annotated. These variants were filtered for novelty and rarity using MAF of 0 and 1%.

SETTING

Hospital in Indonesia.

PATIENTS/PARTICIPANTS

Indonesians with non-syndromic orofacial clefts.

INTERVENTIONS

Deleterious variants were prioritized. Pathogenic amino acid changes effect on protein structure and function were identified using HOPE.

MAIN OUTCOME MEASURE(S)

Rare and novel variants in known cleft genes were filtered and their deleteriousness were predicted using polyphen, SIFT and CADD.

RESULTS

We identified rare (MAF <1%) deleterious variants in 4 craniofacial genes namely MMACHC (rs371937044, MAF = 0.00011). SOS1 (rs190222208, MAF = 0.00045), TULP4 (rs199583035, MAF = 0.067), and MTHFD1L (rs143492706, MAF = 0.0044). MMACHC has a mouse knockout model with facial cleft and failure of palatal fusion. The individual with variant in MMACHC presented with nsCPO.

CONCLUSIONS

Our study provides additional evidence for the role of TULP4, SOS1, MTHFD1L and MMACHC genes in nsOFC development. This is the first time MMACHC is implicated in nsOFC development in humans.

Clinically actionable secondary findings in 130 triads from sub-Saharan African families with non-syndromic orofacial clefts

INTRODUCTION

The frequency and implications of secondary findings (SFs) from genomic testing data have been extensively researched. However, little is known about the frequency or reporting of SFs in Africans, who are underrepresented in large-scale population genomic studies. The availability of data from the first whole-genome sequencing for orofacial clefts in an African population motivated this investigation.

METHODS

In total, 130 case-parent trios were analyzed for SFs within the ACMG SFv.3.0 list genes. Additionally, we filtered for four more genes (HBB, HSD32B, G6PD and ACADM).

RESULTS

We identified 246 unique variants in 55 genes; five variants in four genes were classified as pathogenic or likely pathogenic (P/LP). The P/LP variants were seen in 2.3% (9/390) of the subjects, a frequency higher than ~1% reported for diverse ethnicities. On the ACMG list, pathogenic variants were observed in PRKAG (p. Glu183Lys). Variants in the PALB2 (p. Glu159Ter), RYR1 (p. Arg2163Leu) and LDLR (p. Asn564Ser) genes were predicted to be LP.

CONCLUSION

This study provides information on the frequency and pathogenicity of SFs in an African cohort. Early risk detection will help reduce disease burden and contribute to efforts to increase knowledge of the distribution and impact of actionable genomic variants in diverse populations.

Novel IRF6 variant in orofacial cleft patients from Durban, South Africa

BACKGROUND

To date, there are over 320 variants identified in the IRF6 gene that cause Van der Woude syndrome or popliteal pterygium syndrome. We sequenced this gene in a South African orofacial cleft cohort to identify the causal IRF6 variants in our population.

METHOD

Saliva samples from 100 patients with syndromic and non-syndromic CL ± P were collected. Patients were recruited from the cleft clinics at two public, tertiary hospitals in Durban, South Africa (SA), namely Inkosi Albert Luthuli Central Hospital (IALCH) and KwaZulu-Natal Children's Hospital (KZNCH). We prospectively sequenced the exons of IRF6 in 100 orofacial cleft cases, and where possible, we also sequenced the parents of the individuals to determine the segregation pattern.

RESULTS

Two variants were identified; one novel (p.Cys114Tyr) and one known (p.Arg84His) missense variant in IRF6 gene were identified. The patient with the p.Cys114Tyr variant was non-syndromic with no clinical VWS phenotype expected of individuals with IRF6 coding variants, and the patient with the p.Arg84His had phenotypic features of popliteal pterygium syndrome. The p.Arg84His variant segregated in the family, with the father also being affected.

CONCLUSIONS

This study provides evidence that IRF6 variants are found in the South African population. Genetic counselling is essential for affected families, particularly in the absence of a known clinical phenotype since it helps with the plans for future pregnancies.

Variant analyses of candidate genes in orofacial clefts in multi-ethnic populations

OBJECTIVES

Cleft lip with/without cleft palate and cleft palate only is congenital birth defects where the upper lip and/or palate fail to fuse properly during embryonic facial development. Affecting ~1.2/1000 live births worldwide, these orofacial clefts impose significant social and financial burdens on affected individuals and their families. Orofacial clefts have a complex etiology resulting from genetic variants combined with environmental covariates. Recent genome-wide association studies and whole-exome sequencing for orofacial clefts identified significant genetic associations and variants in several genes. Of these, we investigated the role of common/rare variants in SHH, RORA, MRPL53, ACVR1, and GDF11.

MATERIALS AND METHODS

We sequenced these five genes in 1255 multi-ethnic cleft lip with/without palate and cleft palate only samples in order to find variants that may provide potential explanations for the missing heritability of orofacial clefts. Rare and novel variants were further analyzed using in silico predictive tools.

RESULTS

Ninteen total variants of interest were found, with variant types including stop-gain, missense, synonymous, intronic, and splice-site variants. Of these, 3 novel missense variants were found, one in SHH, one in RORA, and one in GDF11.

CONCLUSION

This study provides evidence that variants in SHH, RORA, MRPL53, ACVR1, and GDF11 may contribute to risk of orofacial clefts in various populations.

Novel GRHL3 Variants in a South African Cohort With Cleft Lip and Palate

OBJECTIVE

The etiology of cleft palate (CP) is poorly understood compared with that of cleft lip with or without palate (CL ± P). Recently, variants in Grainyhead like transcription factor 3 (GRHL3) were reported to be associated with a risk for CP in European and some African populations including Nigeria, Ghana, and Ethiopia. In order to identify genetic variants that may further explain the etiology of CP, we sequenced GRHL3 in a South African population to determine if rare variants in GRHL3 are associated with the presence of syndromic or nonsyndromic CP.

DESIGN

We sequenced the exons of GRHL3 in 100 cases and where possible, we sequenced the parents of the individuals to determine the segregation pattern and presence of de novo variants.

SETTING

The cleft clinics from 2 public, tertiary hospitals in Durban, South Africa (SA), namely Inkosi Albert Luthuli Central Hospital and KwaZulu-Natal Children's Hospital.

PATIENTS, PARTICIPANTS

One hundred patients with CL ± P and their parents.

INTERVENTIONS

Saliva samples were collected.

MAIN OUTCOME MEASURES

To ascertain the genetic variants in the GRHL3 gene in patients with CL ± P in SA.

RESULTS

Five variants in GRHL3 were observed; 3 were novel and 2 were known variants. The novel variants were intronic variants (c.1062 + 77A>G and c.627 + 1G>A) and missense variant (p.Asp169Gly).

CONCLUSIONS

This study provides further evidence that variants in GRHL3 contribute to the risk of nonsyndromic CP in African populations, specifically, in the South African population.

Investigating the relationship between cancer and orofacial clefts using GWAS significant loci for cancers: A case-control and case-triad study

Background

Several population-based case-control studies have reported concurrent presentation of cancer and congenital malformations. Many associations have been made between oral clefting and cancers, though some of these results are conflicting. Some studies have reported an increased risk of cancer among 1st-degree relatives of cleft cases and vice versa, and also an excess risk of cancers of the breast, lung, and brain among those with oral clefts. This study aimed to determine if the genetic polymorphisms found in some cancers are also associated with orofacial cleft in an African cohort.

Methods

The study was a case-control and case-triad study in which cases were 400 individuals clinically diagnosed with non-syndromic cleft lip and/or palate (CL/P), while controls were 450 individuals without CL/P. Samples were obtained from three African countries while DNA extraction, PCR, and genotyping were carried out at the University of Iowa, US. Eleven SNPs in genes coding for SWI/SNF subunits and 13 GWAS significant SNPs for cancers associated with orofacial cleft were selected. Case-control analysis, transmission disequilibrium test (TDT), and DFAM to combine the parent-offspring trio data and unrelated case/control data in a single analysis were carried out using PLINK.

Results

For the case-control analyses that included all the clefts and for the CLP subtype, none of the SNPs were statistically significant. Statistically increased risk for the following SNPs rs34775372 (p = 0.02; OR = 1.54, CI:1.07–2.22), rs55658222 (p = 0.009; OR = 2.64, CI:1.28–5.45) and rs72728755 (p = 0.02; OR=2.27, CI:1.17–4.45) was observed with the CL only sub-group. None of these were significant after Bonferoni correction. In the TDT analyses, a significantly reduced risk with rs10941679 (p = 0.003; OR = 0.43, CI:0.24–0.75) was observed and this was significant after Bonferroni correction. The rs10941679 was also significant (p = 0.003) in the DFAM analyses as well even after Bonferroni correction.

Conclusion

The results from this study represent an important starting point for understanding the concurrent presentation of some cancers in orofacial clefts, and cancer risks in cleft patients. The associations observed warrant further investigation in a larger cohort and will set the stage for a more mechanistic approach toward understanding the risk for cancers in families with clefts.

Periconceptional use of vitamin A and the risk of giving birth to a child with nonsyndromic orofacial clefts-A meta-analysis

BACKGROUND

We conducted a meta-analysis of observational epidemiological studies to evaluate the association between periconceptional use of vitamin A and the risk of giving birth to a child with nonsyndromic orofacial clefts (NSOFCs).

METHODS

We carried out a systematic literature search of Embase, PubMed, Web of Science, Google Scholar, and OpenGrey from inception to June 30, 2021. Two reviewers independently evaluated the studies that met the inclusion criteria and filled out an abstraction form for each study. Study quality was assessed using the Newcastle-Ottawa Assessment Scale (NOS). Adjusted estimates were pooled with an inverse variance weighting using a random-effects model. Heterogeneity and publication bias were assessed using the Cochran's Q test and funnel plot, respectively.

RESULTS

A total of six case-control studies with moderate risk of bias were included. The pooled OR showed a 20% reduction in the risk of NSOFCs for periconceptional use of vitamin A which was not statistically significant (OR = .80; 95% CI .54-1.17, p = .25). For nonsyndromic cleft lip with or without cleft palate (NSCL/P), the studies were homogenous, and the pooled estimate showed a 13% risk reduction, which was significant (OR = .87; 95% CI .77-.99, p = .03). For nonsyndromic cleft palate only (NSCPO), the pooled estimate showed a 33% lower likelihood, which was not statistically significant (OR = .67; 95% CI .42-1.08, p = .10).

CONCLUSION

Our results suggest a possible protective effect for the periconceptional use of vitamin A on the risk of NSCL/P. This finding should be investigated further in prospective studies across multiple populations.

Genetic and Epigenetic Studies in Nonsyndromic Oral Clefts

The etiology of non-syndromic oral clefts (NSOFC) is complex with genetics, genomics, epigenetics and stochastics factors playing a role. Several approaches have been applied to understand the etiology of non-syndromic oral clefts. These include linkage, candidate gene association studies, genome-wide association studies, whole genome sequencing, copy number variations and epigenetics. In this review we shared these approaches, genes and loci reported in some studies.

Genome-wide Gene-by-Sex Interaction Studies Identify Novel Nonsyndromic Orofacial Clefts Risk Locus

Risk loci identified through genome-wide association studies have explained about 25% of the phenotypic variations in nonsyndromic orofacial clefts (nsOFCs) on the liability scale. Despite the notable sex differences in the incidences of the different cleft types, investigation of loci for sex-specific effects has been understudied. To explore the sex-specific effects in genetic etiology of nsOFCs, we conducted a genome-wide gene × sex (GxSex) interaction study in a sub-Saharan African orofacial cleft cohort. The sample included 1,019 nonsyndromic orofacial cleft cases (814 cleft lip with or without cleft palate and 205 cleft palate only) and 2,159 controls recruited from 3 sites (Ethiopia, Ghana, and Nigeria). An additive logistic model was used to examine the joint effects of the genotype and GxSex interaction. Furthermore, we examined loci with suggestive significance (P < 1E-5) in the additive model for the effect of the GxSex interaction only. We identified a novel risk locus on chromosome 8p22 with genome-wide significant joint and GxSex interaction effects (rs2720555, p_2df = 1.16E-08, p_GxSex = 1.49E-09, odds ratio [OR] = 0.44, 95% CI = 0.34 to 0.57). For males, the risk of cleft lip with or without cleft palate at this locus decreases with additional copies of the minor allele (p < 0.0001, OR = 0.60, 95% CI = 0.48 to 0.74), but the effect is reversed for females (p = 0.0004, OR = 1.36, 95% CI = 1.15 to 1.60). We replicated the female-specific effect of this locus in an independent cohort (p = 0.037, OR = 1.30, 95% CI = 1.02 to 1.65), but no significant effect was found for the males (p = 0.29, OR = 0.86, 95% CI = 0.65 to 1.14). This locus is in topologically associating domain with craniofacially expressed and enriched genes during embryonic development. Rare coding mutations of some of these genes were identified in nsOFC cohorts through whole exome sequencing analysis. Our study is additional proof that genome-wide GxSex interaction analysis provides an opportunity for novel findings of loci and genes that contribute to the risk of nsOFCs.