Nonsyndromic cleft lip and/or cleft palate (NSCL/P) is an isolated phenotype of orofacial clefts with skewed sex ratio in prevalence. This study aims to identify differentially expressed genes (DEGs) and microRNAs (DEMs) of NSCL/P by integrated bioinformatics analysis, revealing mechanisms for sexual dimorphism in prevalence.
Materials and Methods:
First, we downloaded the expression profile data from Gene Expression Omnibus database to identify DEGs and DEMs. Second, Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses performed DEGs’ functions. Then, clustered DEGs were identified through protein–protein interaction networks. Combining clustered DEGs with key genes searched in GeneCards enlarged NSCL/P-related genes. Moreover, the genes were linked by transcription factors (TFs). Subsequently, connected by the above TFs, DEMs and genes were used to establish the miRNA-TF-messenger RNA (mRNA) regulatory networks.
Results:
The DEGs in sex-ignored group, female-only group, and male-only group were obtained, respectively. Among the DEMs, miR-378 was downregulated in females but upregulated in males. In female-only group, the miRNA-TF-mRNA regulatory networks showed miR-378-SP1-POLE2/CDK6/EZR regulatory axis was found to be key candidates of NSCL/P.
Conclusions:
Our findings suggest that different expression of miR-378 is consistent with the skewed sex ratio in the prevalence of NSCL/P.
AbramyanJRichmanJM. Recent insights into the morphological diversity in the amniote primary and secondary palates. Dev Dyn. 2015;244(12):1457–1468. doi:10.1002/dvdy.24338
BeatyTHMurrayJCMarazitaMLMungerRGRuczinskiIHetmanskiJBLiangKYWuTMurrayTFallinMD, et al.A genome-wide association study of cleft lip with and without cleft palate identifies risk variants near MAFB and ABCA4. Nat Genet. 2010;42(6):525–529. doi:10.1038/ng.580
4.
BeatyTHRuczinskiIMurrayJCMarazitaMLMungerRGHetmanskiJBMurrayTRedettRJFallinMDLiangKY, et al.Evidence for gene-environment interaction in a genome wide study of nonsyndromic cleft palate. Genet Epidemiol. 2011;35(6):469–478. doi:10.1002/gepi.20595
5.
ChenGLiMXWangHXHongJWShenJYWangQShiQMGeXDingZZhangJP, et al.Identification of key genes in cleft lip with or without cleft palate regulated by miR-199a-5p. Int J Pediatr Otorhinolaryngol. 2018;111:128–137. doi:10.1016/j.ijporl.2018.06.005
6.
ChouCHShresthaSYangCDChangNWLinYLLiaoKWHuangWCSunTHTuSJLeeWH, et al.miRTarBase update 2018: a resource for experimentally validated microRNA-target interactions. Nucleic Acids Res. 2018;46(D1):D296–D302. doi:10.1093/nar/gkx1067
7.
CuraFPalmieriAGirardiACarinciFMorselliPGNouriNPezzettiFScapoliLMartinelliM. Possible effect of SNAIL family transcriptional repressor 1 polymorphisms in non-syndromic cleft lip with or without cleft palate. Clin Oral Investig. 2018;22(7):2535–2541. doi:10.1007/s00784-018-2350-0
8.
DhulipalaVCWelshonsWVReddyCS. Inhibition of human embryonic palatal mesenchymal cell cycle by secalonic acid D: a probable mechanism of its cleft palate induction. Orthod Craniofac Res. 2004;7(4):227–236.
9.
DixonMJMarazitaMLBeatyTHMurrayJC. Cleft lip and palate: understanding genetic and environmental influences. Nat Rev Genet. 2011;12(3):167–178. doi:10.1038/nrg2933
10.
FonsecaRFde CarvalhoFMPolettaFAMontanerDDopazoJMerebJCMoreiraMAMSeuanezHNVieiraARCastillaEEOrioliIM. Family-based genome-wide association study in Patagonia confirms the association of the DMD locus and cleft lip and palate. Eur J Oral Sci. 2015;123(5):381–384. doi:10.1111/eos.12212
11.
GrosselMJHindsPW. Beyond the cell cycle: a new role for Cdk6 in differentiation. J Cell Biochem. 2006;97(3):485–493. doi:10.1002/jcb.20712
12.
GuzziPHDi MartinoMTTagliaferriPTassonePCannataroM. Analysis of miRNA, mRNA, and TF interactions through network-based methods. EURASIP J Bioinform Syst Biol. 2015;2015:4. doi:10.1186/s13637-015-0023-8
13.
JugessurASkareOLieRTWilcoxAJChristensenKChristiansenLNguyenTTMurrayJCGjessingHK. X-linked genes and risk of orofacial clefts: evidence from two population-based studies in Scandinavia. PLoS One. 2012;7(6):e39240. doi:10.1371/journal.pone.0039240
14.
KobayashiGSAlviziLSunagaDYFrancis-WestPKutaAAlmadaBVFerreiraSGde Andrade-LimaLCBuenoDFRaposo-AmaralCE, et al.Susceptibility to DNA damage as a molecular mechanism for non-syndromic cleft lip and palate. PLoS One. 2013;8(6):e65677. doi:10.1371/journal.pone.0065677
15.
LiYAsaharaHPatelVSZhouSLinnS. Purification, cDNA cloning, and gene mapping of the small subunit of human DNA polymerase epsilon. J Biol Chem. 1997;272(51):32337–32344. doi:10.1074/jbc.272.51.32337
16.
LiuFKorcM. Cdk4/6 inhibition induces epithelial-mesenchymal transition and enhances invasiveness in pancreatic cancer cells. Mol Cancer Ther. 2012;11(10):2138–2148. doi:10.1158/1535-7163.MCT-12-0562
17.
LiuJYeXWuFX. Characterizing dynamic regulatory programs in mouse lung development and their potential association with tumourigenesis via miRNA-TF-mRNA circuits. BMC Syst Biol. 2013;7(suppl 2):S11. doi:10.1186/1752-0509-7-S2-S11
18.
MachicoaneMde FrutosCAFinkJRocancourtMLombardiYGarelSPielMEchardA. SLK-dependent activation of ERMs controls LGN-NuMA localization and spindle orientation. J Cell Biol. 2014;205(6):791–799. doi:10.1083/jcb.201401049
19.
MasottiCBritoLANicaACLudwigKUNunesKSavastanoCPMalcherCFerreiraSGKobayashiGSBuenoDF, et al.MRPL53, a new candidate gene for orofacial clefting, identified using an eQTL approach. J Dent Res. 2018;97(1):33–40. doi:10.1177/0022034517735805
20.
MosseyPACatillaEE. Global Registry and Database on Craniofacial Anomalies: Report of a WHO Registry Meeting on Craniofacial Anomalies; 2003.
21.
MosseyPALittleJMungerRGDixonMJShawWC. Cleft lip and palate. Lancet. 2009;374(9703):1773–1785. doi:10.1016/S0140-6736(09)60695-4
22.
RahimovFJugessurAMurrayJC. Genetics of nonsyndromic orofacial clefts. Cleft Palate Craniofac J. 2012;49(1):73–91. doi:10.1597/10-178
23.
RebhanMChalifa-CaspiVPriluskyJLancetD. GeneCards: integrating information about genes, proteins and diseases. Trends Genet. 1997;13(4):163.
24.
RitchieMEPhipsonBWuDHuYLawCWShiWSmyth GK. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47. doi:10.1093/nar/gkv007
25.
ShiHZWangDSunXNShengL. MicroRNA-378 acts as a prognosis marker and inhibits cell migration, invasion and epithelial-mesenchymal transition in human glioma by targeting IRG1. Eur Rev Med Pharmacol Sci. 2018;22(12):3837–3846. doi:10.26355/eurrev_201806_15268
26.
ShuXShuSChengH. Genome-wide mRNA-seq profiling reveals that LEF1 and SMAD3 regulate epithelial-mesenchymal transition through the hippo signaling pathway during palatal fusion. Genet Test Mol Biomarkers. 2019;23(3):197–203. doi:10.1089/gtmb.2018.0221
27.
SkareOGjessingHKGjerdevikMHaalandOARomanowskaJLieRTJugessurA. A new approach to chromosome-wide analysis of X-linked markers identifies new associations in Asian and European case-parent triads of orofacial clefts. PLoS One. 2017;12(9):e0183772. doi:10.1371/journal.pone.0183772
28.
SlikKKurkiSKorpelaTCarpenOKorkeilaESundstromJ. Ezrin expression combined with MSI status in prognostication of stage II colorectal cancer. PLoS One. 2017;12(9):e0185436. doi:10.1371/journal.pone.0185436
29.
SunDBaurSHayED. Epithelial-mesenchymal transformation is the mechanism for fusion of the craniofacial primordia involved in morphogenesis of the chicken lip. Dev Biol. 2000;228(2):337–349. doi:10.1006/dbio.2000.9946
30.
TakigawaTShiotaK. Terminal differentiation of palatal medial edge epithelial cells in vitro is not necessarily dependent on palatal shelf contact and midline epithelial seam formation. Int J Dev Biol. 2004;48(4):307–317. doi:10.1387/ijdb.041840tt
31.
VaquerizasJMKummerfeldSKTeichmannSALuscombeNM. A census of human transcription factors: function, expression and evolution. Nat Rev Genet. 2009;10(4):252–263. doi:10.1038/nrg2538
32.
WangHQiuTShiJLiangJWangYQuanLZhangYZhangQTaoK. Gene expression profiling analysis contributes to understanding the association between non-syndromic cleft lip and palate, and cancer. Mol Med Rep. 2016;13(3):2110–2116. doi:10.3892/mmr.2016.4802
33.
WatkinsSEMeyerREStraussRPAylsworthAS. Classification, epidemiology, and genetics of orofacial clefts. Clin Plast Surg. 2014;41(2):149–163. doi:10.1016/j.cps.2013.12.003
34.
WehbyGLCassellCH. The impact of orofacial clefts on quality of life and healthcare use and costs. Oral Dis. 2010;16(1):3–10. doi:10.1111/j.1601-0825.2009.01588.x
35.
WiseASShiMWeinbergCR. Family-based multi-SNP X chromosome analysis using parent information. Front Genet. 2016;7:20. doi:10.3389/fgene.2016.00020
36.
WyszynskiDF. Cleft Lip and Palate: From Origin to Treatment. Oxford University Press; 2002.
37.
XuZLieRTWilcoxAJSaugstadODTaylorJA. A comparison of DNA methylation in newborn blood samples from infants with and without orofacial clefts. Clin Epigenetics. 2019;11(1):40. doi:10.1186/s13148-019-0638-9
38.
YuGWangLGHanYHeQY. Cluster profiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–287. doi:10.1089/omi.2011.0118
39.
ZengMZhuLLiLKangC. miR-378 suppresses the proliferation, migration and invasion of colon cancer cells by inhibiting SDAD1. Cell Mol Biol Lett. 2017;22:12. doi:10.1186/s11658-017-0041-5
40.
ZhangXMGuoLChiMHSunHMChenXW. Identification of active miRNA and transcription factor regulatory pathways in human obesity-related inflammation. BMC Bioinformatics. 2015;16(1):76. doi:10.1186/s12859-015-0512-5
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.
