Abstract
Long non-coding RNAs (lncRNAs) are thought to play important roles in non-syndromic orofacial clefts (NSOFC). Clinical diagnosis was categorized as either non-syndromic cleft lip with or without cleft palate (NSCL/P), or non-syndromic cleft palate only (NSCPO). Tissues excised from the trimmed wound edge were reserved as experimental samples; adjacent normal control was used as a positive control, and tissue from healthy individuals was used as a blank control. Target lncRNAs in the collected tissues were identified using microarrays and quantitative reverse transcription PCR (RT-qPCR). Immunohistochemical (IHC) staining and RT-qPCR were used to verify the target mRNAs. Pathway, gene ontology (GO) enrichment, and TargetScan predictions were employed to construct competing endogenous RNA networks (ceRNA networks) and explore their potential functions. RNA-Seq revealed 24 upregulated and 43 downregulated lncRNAs; MALAT1 and NEAT1 were screened and validated using RT-qPCR. Common NSOFC risk factors were positively correlated with MALAT1 and NEAT1 expression. Bioinformatics predicted four ceRNA networks; GO enrichment focused on their potential functions. RT-qPCR and IHC data were consistent with respect to expression levels of proteins and the mRNAs that encode them. As MALAT1 and NEAT1 are associated with the severity of NSOFC, they represent potential therapeutic targets and prognostic biomarkers.