Integrated analysis of lncRNA and mRNA expression profiles in the submandibular glands of DIO mice

OBJECTIVE Obesity contributes to the dysfunction of salivary gland. To explore the specific underlying mechanism for obesity-induced hyposalivation, a model for high-fat diet-induced obese (DIO) mice were constructed to analyze long non-coding RNAs (lncRNAs) and messenger RNAs (mRNAs) expression profiles. METHODS The DIO group and control group were fed a diet containing 60 kcal% fat and a normal chow diet for 16 weeks respectively. Microarray analyses were performed to detect the expression profiles of lncRNA and mRNA in submandibular gland tissues from control group mice and DIO mice. Gene ontology, kyoto encyclopedia of genes and genomes, protein-protein interaction, coding-non-coding gene co-expression, transcription factors and competing endogenous RNA analyses were performed to examine the function of differentially expressed genes. RESULTS Microarray analyses identified that 624 lncRNAs, along with 297 mRNAs were differentially expressed. Bioinformatic analyses revealed that "complement and coagulation cascades," "glutathione metabolism," "cysteine and methionine metabolism," and "estrogen signaling pathway" were significantly associated with candidate lncRNAs. Transcription factors analysis on candidate lncRNAs revealed several genes such as tribbles pseudokinase 3 may play regulatory roles. CONCLUSIONS Our results revealed the expression profiles of lncRNAs and mRNAs and provided new insights into the mechanism of obesity-induced hyposalivation using bioinformatic analyses.