Identification of competing endogenous RNA networks associated with circRNA and lncRNA in TCDD-induced cleft palate development

2,3,7,8 -tetrachlorodibenzo-p-dioxin (TCDD) is a teratogen that can induce cleft palate formation, a common birth defect. Competing endogenous RNAs (ceRNAs), including circular RNAs (circRNAs) and long non-coding RNAs (lncRNAs), indirectly regulate gene expression via sharing microRNAs (miRNAs). Nevertheless, the mechanism by which they act as ceRNAs to regulate palatal development remains to be explored in greater detail. Here, the cleft palate model of C57BL/6 N pregnant mice was constructed by gavage of TCDD (64 ug/kg) on gestation day (GD) 10.5, and the palatal shelves were taken on gestation day (GD) 14.5 for whole-transcriptome sequencing to investigate the underlying mechanisms of the roles of circRNAs and lncRNAs as ceRNAs in cleft palate. Sequencing results revealed that 293 lncRNA, 589 circRNA, 47 miRNA, and 138 messenger RNA (mRNA) were significantly dysregulated, and the cytochrome P450 (CYP) enzymes and the aryl hydrocarbon receptor (AhR) pathway play key roles in the induction of cleft palate upon exposure to TCDD. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed the function of TCDD function was mainly related to the metabolic processes of intracellular compounds, including the metabolic processes of cellular aromatic compounds and the metabolism of exogenous drugs by cytochrome P450, etc. Furthermore, quantitative reverse transcription polymerase chain reaction (qRT-PCR) indicated that the circRNA_1781/miR-30c-1–3p/PKIB and XR_380026.2/miR-1249–3p/DNAH10 ceRNA networks were hypothesized to be a hub involved in palatal development suggesting that the circRNA_1781/miR-30c-1–3p/PKIB and XR_380026.2/miR-1249–3p/DNAH10 ceRNA networks may be critical for palatogenesis, setting the foundation for the investigation of cleft palate.