Nonsyndromic cleft lip and/or cleft palate (NSCL/P) have multifactorial etiology where genetic factors, gene–environment interactions, stochastic factors, gene–gene interactions, and parent-of-origin effects (POEs) play cardinal roles. POEs arise when the parental origin of alleles differentially impacts the phenotype of the offspring. The aim of this study was to identify POEs that can increase risk for NSCL/P in humans using a genome-wide dataset.
Methods
The samples (174 case-parent trios from Ghana, Ethiopia, and Nigeria) included in this study were from the African only genome wide association studies (GWAS) that was published in 2019. Genotyping of individual DNA using over 2 million multiethnic and African ancestry-specific single-nucleotide polymorphisms from the Illumina Multi-Ethnic Genotyping Array v2 15070954 A2 (genome build GRCh37/hg19) was done at the Center for Inherited Diseases Research. After quality control checks, PLINK was employed to carry out POE analysis employing the pooled subphenotypes of NSCL/P.
Results
We observed possible hints of POEs at a cluster of genes at a 1 mega base pair window at the major histocompatibility complex class 1 locus on chromosome 6, as well as at other loci encompassing candidate genes such as ASB18, ANKEF1, AGAP1, GABRD, HHAT, CCT7, DNMT3A, EPHA7, FOXO3, lncRNAs, microRNA, antisense RNAs, ZNRD1, ZFAT, and ZBTB16.
Conclusion
Findings from our study suggest that some loci may increase the risk for NSCL/P through POEs. Additional studies are required to confirm these suggestive loci in NSCL/P etiology.
BoylesALWilcoxAJTaylorJAMeyerKFredriksenAUelandPMDrevonCAVollsetSELieRT. Folate and one-carbon metabolism gene polymorphisms and their associations with oral facial clefts. Am J Med Genet A. 2008;146(4):440-449.
2.
ButaliAMosseyPAAdeyemoWLEsheteMAGowansLJJBuschTDJainDYuWHuanLLaurieCA, et al. Genomic analyses in African populations identify novel risk loci for cleft palate. Hum Mol Genet. 2019;28(6):1038-1051.
3.
ButaliASuzukiSCooperMEMansillaAMCuencoKLeslieEJSuzukiYNiimiTYamamotoMAyangaG, et al. Replication of genome wide association identified candidate genes confirm the role of common and rare variants in PAX7 and VAX1 in the etiology of non-syndromic CL(P). Am J Med Genet A. 2013;161(5): 965-972.
4.
CarlsonJCAnandDButaliABuxoCJChristensenKDeleyiannisFHechtJTMorenoLMOrioliIMPadillaC, et al. A systematic genetic analysis and visualization of phenotypic heterogeneity among orofacial cleft GWAS signals. Genet Epidemiol. 2019;43(6):704-716.
5.
ConomosMPThorntonT. GENESIS: gENetic EStimation and inference in structured samples (GENESIS): statistical methods for analyzing genetic data from samples with population structure and/or relatedness. R package version 2.4.0, 2016. https://bioc.ism.ac.jp/packages/3.2/bioc/html/GENESIS.html.
6.
DuanSJHuangNZhangBHShiJYHeSMaJYuQQShiBJiaZL. New insights from GWAS for the cleft palate among Han Chinese population. Med Oral Patol Oral Cir Bucal. 2017;22(2):e219-e227.
7.
GargPLudwigKUBöhmerACRubiniMSteegers-TheunissenRMosseyPAMangoldESharpAJ. Genome-wide analysis of parent-of-origin effects in non-syndromic orofacial clefts. Euro J Hum Genet. 2014;22(6):822-830.
8.
GjerdevikMHaalandOARomanowskaJLieRTJugessurAGjessingHK. Parent-of-origin-environment interactions in case parent triads with or without independent controls. Ann Hum Genet. 2018;82(2):60-73.
9.
GogartenSMBhangaleTConomosMPLaurieCAMcHughCPPainterIZhengXCrosslinDRLevineDLumleyT, et al. GWAS tools: an R/bioconductor package for quality control and analysis of genome-wide association studies. Bioinformatics2012;28(24):3329-3331.
10.
HaalandØAJugessurAGjerdevikMRomanowskaJShiMBeatyTH, Marazita ML, Murray JC, Wilcox AJ, Liet RT, et al. Genome-wide analysis of parent-of-origin interaction effects with environmental exposure (PoOxE): an application to European and Asian cleft palate trios. PLoS ONE2017;12(9):e0184358.
11.
HaalandØARomanowskaJGjerdevikMLieRTGjessingHKJugessurA. A genome-wide scan of cleft lip triads identifies parent-of-origin interaction effects between ANK3 and maternal smoking, and between ARHGEF10 and alcohol consumption. F1000Res.2019;8:960.
12.
HuangLJiaZShiYDuQShiJWangZMouYWangQZhangBWangQ, et al. Genetic factors define CPO and CLO subtypes of nonsyndromic orofacial cleft. PLoS Genet. 2019;15(10):e1008357.
13.
JungSHLeeAYParkJWBaekSHKimYH. Investigation of parental transmission of RUNX2 single nucleotide polymorphism and its association with nonsyndromic cleft lip with or without palate. Cleft Palate–Craniofac. J2014;51(2): 234-239.
14.
KhanMFJLittleJMosseyPASteegers-TheunissenRPMBonsiMAndreasiRBRubiniM. Association between a common missense variant in LOXL3 gene and the risk of non-syndromic cleft palate. Congenital Anom.2018;58(4):136-140.
15.
LaurieCCDohenyKFMirelDBPughEWBierutLJBhangaleTBoehmFCaporasoNECornelisMCEdenbergHJ, et al. Quality control and quality assurance in genotypic data for genome-wide association studies. Genet Epidemiol. 2010;34(6):591-602.
16.
LawsonHACheverudJMWolfJB. Genomic imprinting and parent-of-origin effects on complex traits. Nat Rev Genet. 2013;14(9): 609-617.
17.
LonjouCZhangWCollinsATapperWJElahiEManiatisNMortonNE. Linkage disequilibrium in human populations. Proc Natl Acad Sci USA. 2003;100(13):8036.
18.
MailiLLetraASilvaRBuchananEPMullikenJBGreivesMRTeichgraeberJFBlackwellSJUmmerRWeberR,et al.PBX-WNT-P63-IRF6 pathway in nonsyndromic cleft lip and palate. Birth Defects Res. 2020;112(3):234-244.
19.
Moreno UribeLMFominaTMungerRGRomittiPAJenkinsMMGjessingHKGjerdevikMChristensenKWilcoxAJMurrayJC, et al. A population-based study of effects of genetic loci on orofacial clefts. J Dent Res. 2017;96(11):1322-1329.
20.
MosseyPAModellB. Epidemiology of oral clefts 2012: an international perspective. Front Oral Biol. 2012;16:1-18.
21.
NideyNMoreno UribeLMMarazitaMMWehbyGL. Psychosocial well-being of parents of children with oral clefts. Child Care Health Dev. 2016;42(1):42-50.
22.
NideyNWehbyG. Barriers to health care for children with orofacial clefts: a systematic literature review and recommendations for research priorities. Oral Health Dent Stud. 2019;2(1):2.
23.
OseniGODeeptiJMosseyPABuschTDGowansLJJEsheteMAAdeyemoWLLaurieCALaurieCCOwaisA, et al. Identification of paternal uniparental disomy on chromosome 22 and a de-novo deletion on chromosome 18 in individuals with orofacial clefts. Mol Genet Genomic Med. 2018;6(6):924-932.
24.
ParkBYSullJWParkJYJeeSHBeatyTH. Differential parental transmission of markers in BCL3 among Korean cleft case-parent trios. J Prev Med Public Health. 2009;42(1):1-4.
25.
PetersJ. The role of genomic imprinting in biology and disease: an expanding view. Nat Rev Genet. 2014;15:517-530.
26.
PurcellSNealeBTodd-BrownKThomasLFerreiraMARBenderDMallerJSklarPde BakkerPIWDalyMJ,et al.PLINK: a toolset for whole-genome association and population-based linkage analysis. Am J Hum Genet. 2007;81(3):559-575.
27.
RajuGTLakkakulaBVKSMurthyJKannanMAPaulSFD. Transmission analysis of TGFB1 gene polymorphisms in non-syndromic cleft lip with or without cleft palate. Int J Ped Otorhinolaryn. 2017;100(2017):14e17.
28.
ReutterHBirnbaumSMendeMLausterCSchmidtGHenschkeHSaffarMMartiniMLausterRSchiefkeF, et al. TGFB3 displays parent-of-origin effects among central Europeans with nonsyndromic cleft lip and palate. J Hum Genet. 2008;53(7):656-661.
29.
ShiMMurrayJCMarazitaMLMungerRGRuczinskiIHetmanskiJBWuTMurrayTRedettRJWilcoxAJ, et al. Genome wide study of maternal and parent-of-origin effects on the etiology of orofacial clefts. Am J Med Genet A. 2012;158A(4):784-794.
30.
SkareØLieRTHaalandØAGjerdevikMRomanowskaJGjessingHKJugessurA. Analysis of parent-of-origin effects on the X chromosome in Asian and European orofacial cleft triads identifies associations with DMD, FGF13, EGFL6, and additional loci at Xp22.2. Front Genet. 2018;9:25.
31.
SuazoJSantosJLJaraLBlancoR. Parent-of-origin effects for MSX1 in a Chilean population with nonsyndromic cleft lip/palate. Am J Med Genet Part A. 2010;152A(8):2011-2016.
32.
SullJWLiangKYHetmanskiJBFallinMDIngersollRGParkJWu-ChouYHChenPKChongSSCheahF, et al. Differential parental transmission of markers in RUNX2 among cleft case-parent trios from four populations. Genet Epidemiol. 2008;32(6):505-512.
33.
SullJWLiangKYHetmanskiJBFallinMDIngersollRGParkJWu-ChouYHChenPKChongSSCheahF, et al. Maternal transmission effects of the PAX genes among cleft case–parent trios from four populations. Euro J Hum Genet. 2009b;17(6):831-839.
34.
SullJWLiangKYHetmanskiJBWuTFallinMDIngersollRGParkJWWu-ChouYHChenPKChongSS, et al. Evidence that TGFA influences risk to cleft lip with/without cleft palate through unconventional genetic mechanisms. Hum Genet. 2009a; 126(3):385-394.
35.
ThomasD. Gene–environment-wide association studies: emerging approaches. Nat Rev Genet. 2010;11(4):259-272.
36.
van RooijIALMLudwigKUWelzenbachJIshorstNThonissenMGaleslootTEOngkosuwitoEBergéSJAldhoraeKRojas-MartinezA, et al. 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. Genes (Basel).2019;10(12):1023.
37.
WangHHetmanskiJBRuczinskiILiangKYFallinMDRedettRJRaymondGVWu-ChouY-HChenPK-TYeowV, et al. ROR2 Gene is associated with risk of non-syndromic cleft palate in an Asian population. Chin Med J (Engl). 2012;125(3):476-480.
38.
WeiWHHemaniGHaleyCS. Detecting epistasis in human complex traits. Nat Rev Genet. 2014;15(11):722-733.
39.
WuDWangMWangXYinNSongTLiHZhangFZhangHYeZYuJ, et al. Maternal transmission effect of a PDGF-C SNP on nonsyndromic cleft lip with or without palate from a Chinese population. PLoS ONE2012;7(9):e46477.
40.
YinBShiBJiaZL. Associations among PRDM16 polymorphisms, environmental exposure factors during mother’s pregnancy, and nonsyndromic cleft lip with or without cleft palate. Hua Xi Kou Qiang Yi Xue Za Zhi. 2018;36(5):503-507.
41.
ZengYAmadorCXiaCMarioniRSproulDWalkerRMMorrisSWBretherickACanela-XandriOBoutinTS, et al. Parent of origin genetic effects on methylation in humans are common and influence complex trait variation. Nature Comm. 2019; 10(1):1383.
42.
ZhengXLevineDShenJGogartenSMLaurieCWeirBS. A high-performance computing toolset for relatedness and principal component analysis of SNP data. Bioinformatics2012;28(24):3326-3328.
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
