Integrated analysis of microRNAs, circular RNAs, long non-coding RNAs, and mRNAs revealed competing endogenous RNA networks involved in brown adipose tissue whitening in rabbits.

BACKGROUND: The brown adipose tissue (BAT) is a target for treating obesity. BAT losses thermogenic capacity and gains a "white adipose tissue-like" phenotype ("BAT whitening") under thermoneutral environments, which is a potential factor causing a low curative effect in BAT-related obesity treatments. Circular RNAs (circRNAs) and long non-coding RNAs (lncRNAs) can act as competing endogenous RNAs (ceRNA) to mRNAs and function in various processes by sponging shared microRNAs (miRNAs). However, the roles of circRNA- and lncRNA-related ceRNA networks in regulating BAT whitening remain litter known. RESULTS: In this study, BATs were collected from rabbits at day0 (D0), D15, D85, and 2 years (Y2). MiRNA-seq was performed to investigate miRNA changes during BAT whitening. Then, a combined analysis of circRNA-seq and whole-transcriptome sequencing was used for circRNA assembly and quantification during BAT whitening. Our data showed that 1187 miRNAs and 6204 circRNAs were expressed in the samples, and many of which were identified as significantly changed during BAT whitening. Target prediction showed that D0-selective miRNAs were significantly enriched in the Ras, MAPK, and PI3K-Akt signaling pathways, and Y2-selective miRNAs were predicted to be involved in cell proliferation. The cyclization of several circRNAs could form novel response elements of key thermogenesis miRNAs at the back-splicing junction (BSJ) sites, and in combination with a dual-luciferase reporter assay confirmed the binding between the BSJ site of novel_circ_0013792 and ocu-miR-378-5p. CircRNAs and lncRNAs have high cooperativity in sponging miRNAs during BAT whitening. Both circRNA-miRNA-mRNA and lncRNA-miRNA-mRNA triple networks were significantly involved in immune response-associated biological processes. The D15-selective circRNA might promote BAT whitening by increasing the expression of IDH2. The Y2-selective circRNA-related ceRNA network and lncRNA-related ceRNA network might regulate the formation of the WAT-like phenotype of BAT via MAPK and Ras signaling pathways, respectively. CONCLUSIONS: Our work systematically revealed ceRNA networks during BAT whitening in rabbits and might provide new insight into BAT-based obesity treatments.

Genome-wide identification and characterization of long non-coding RNAs during differentiation of visceral preadipocytes in rabbit.

Emerging evidence suggests that long non-coding RNAs (lncRNAs) are critical regulators of diverse biological processes, including adipogenesis. Despite being considered an ideal animal model for studying adipogenesis, little is known about the roles of lncRNAs in the regulation of rabbit preadipocyte differentiation. In the present study, visceral preadipocytes isolated from newborn rabbits were cultured in vitro and induced for differentiation, and global lncRNA expression profiles of adipocytes collected at days 0, 3, and 9 of differentiation were analyzed by RNA-seq. A total of 2066 lncRNAs were identified from nine RNA-seq libraries. Compared to protein-coding transcripts, lncRNA transcripts exhibited characteristics of a longer length and lower expression level. Furthermore, 486 and 357 differentially expressed (DE) lncRNAs were identified when comparing day 3 vs. day 0 and day 9 vs. day 3, respectively. Target genes of DE lncRNAs were predicted by the cis-regulating approach. Prediction of functions revealed that DE lncRNAs when comparing day 3 vs. day 0 were involved in gene ontology (GO) terms of developmental growth, growth, developmental cell growth, and stem cell proliferation, and involved in Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of PI3K-Akt signaling pathway, fatty acid biosynthesis, and the insulin signaling pathway. The DE lncRNAs when comparing day 9 vs. day 3 were involved in GO terms that associated with epigenetic modification and were involved in the KEGG pathway of cAMP signaling pathway. This study provides further insight into the regulatory function of lncRNAs in rabbit visceral adipose and facilitates a better understanding of different stages of preadipocyte differentiation.

Genome-wide identification and characterization of long non-coding RNAs during postnatal development of rabbit adipose tissue.

BACKGROUND: The rabbit is widely used as an important experimental model for biomedical research, and shows low adipose tissue deposition during growth. Long non-coding RNAs (lncRNAs) are associated with adipose growth, but little is known about the function of lncRNAs in the rabbit adipose tissue. METHODS: Deep RNA-sequencing and comprehensive bioinformatics analyses were used to characterize the lncRNAs of rabbit visceral adipose tissue (VAT) at 35, 85 and 120 days after birth. Differentially expressed (DE) lncRNAs were identified at the three growth stages by DESeq. The cis and trans prediction ways predicted the target genes of the DE lncRNAs. To explore the function of lncRNAs, Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed on the candidate genes. RESULTS: A total of 991,157,544 clean reads were generated after RNA-Seq of the three growth stages, of which, 30,353 and 107 differentially expressed (DE) lncRNAs were identified. Compared to the protein-coding transcripts, the rabbit lncRNAs shared some characteristics such as shorter length and fewer exons. Cis and trans target gene prediction revealed, 43 and 64 DE lncRNAs respectively, corresponding to 72 and 20 protein-coding genes. GO enrichment and KEGG pathway analyses revealed that the candidate DE lncRNA target genes were involved in oxidative phosphorylation, glyoxylate and dicarboxylate metabolism, and other adipose growth-related pathways. Six DE lncRNAs were randomly selected and validated by q-PCR. CONCLUSIONS: This study is the first to profile the potentially functional lncRNAs in the adipose tissue growth in rabbits, and contributes to our understanding of mammalian adipogenesis.