Integrated analysis of noncoding RNAs and mRNAs reveals their potential roles in the biological activities of the growth hormone receptor.

Abstract Accumulating evidence has indicated that noncoding RNAs (ncRNAs) have important regulatory potential in various biological processes. The molecular mechanisms by which growth hormone receptor (GHR) deficiency protects against age-related pathologies, reduces the incidence and delays the occurrence of fatal neoplasms are unclear. The aim of this study was to investigate miRNA, lncRNA and mRNA expression profiles and the potential functional roles of these RNA molecules in GHR knockout (GHR-KO) mice. Microarray expression profiles of miRNAs, lncRNAs and mRNAs were determined in wild type control mice and in GHR-KO mice. Differential expression, pathway and gene network analyses were developed to identify the possible biological roles of functional RNA molecules. Compared to wild type control mice, 1695 lncRNAs, 914 mRNAs and 9 miRNAs were upregulated and 1747 lncRNAs, 786 mRNAs and 21 miRNAs were downregulated in female GHR-KO mice. Moreover, 1265 lncRNAs, 724 mRNAs and 41 miRNAs were upregulated and 1377 lncRNAs, 765 mRNAs and 16 miRNAs were downregulated in male GHR-KO mice compared to wild type mice. Co-expression analysis of mRNAs, lncRNAs, and miRNAs showed that mRNAs including Hemxi2, Ero1Ib, 4933434i20RIK, Pde7a and Lgals1, lncRNAs including ASMM9PARTA014848, EL605414-P1, ASMM9PARTA051724, ASMM9PARTA045378 and ASMM9PARTA049185, and miRNAs including miR-188-3p, miR-690, miR-709 and miR-710 are situated at the core position of a three-dimensional lncRNA–mRNA–miRNA regulatory network. KEGG analysis showed that the most significantly regulated pathway was steroid hormone biosynthesis. We identified a set of lncRNAs, miRNAs and mRNAs that were aberrantly expressed in GHR-KO mice. Our results provide a foundation and an expansive view of the biological activities of the GHR.