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.

Effects of variants in proopiomelanocortin and neuropeptide y genes on growth, carcass, and meat quality traits in rabbits.

Appetite-related neuropeptides proopiomelanocortin (POMC) and Neuropeptide Y (NPY) are essential for regulating feeding behavior and energy homeostasis. The objective of this study was to evaluate the effects of variants in POMC and NPY genes on growth, carcass and meat quality traits in rabbits. A total of six SNPs were identified for POMC (n = 2) and NPY (n = 4) genes by direct sequencing. Three SNPs were subsequently genotyped by using MassArray system (Sequenom iPLEXassay) in 235 individuals, which belong to three meat rabbit breeds, including 93 Ira rabbits; 81 Champagne rabbits and 61 Tianfu black rabbits. The SNP c.112-12G>T was in intron-exon boundaries (intron 1) of POMC gene, and the association analysis showed that individuals with TT genotype had a greater 84 d body weight (BW84), eviscerated weight and semi-eviscerated weight than those with GT genotype (p<0.05); the TT individuals were also higher than those GG in the ripe meat ratio (RMR) (p<0.05). The g.1778G>C SNP, which was in complete linkage with other three SNPs (g.1491G>A, g.1525G>T and g.1530C>T) in intron 1 of NPY gene, was significantly correlated with eviscerated slaughter percentage and semi-eviscerated slaughter percentage in rabbits, and the individuals with CC genotype had a better performance than CG genotype (p<0.05). These findings would provide primary clues for the biological roles of POMC and NPY underlying the rabbit growth-related traits.

Liver Transcriptome Changes of Hyla Rabbit in Response to Chronic Heat Stress.

Rabbit is an economically important farm animal in China and also is a widely used animal model in biological researches. Rabbits are very sensitive to the environmental conditions, therefore we investigated the liver transcriptome changes in response to chronic heat stress in the present study. Six Hyla rabbits were randomly divided into two groups: chronic heat stress (HS) and controls without heat stress (CN). Six RNA-Seq libraries totally yielded 380 million clean reads after the quality filtering. Approximately 85.07% of reads were mapped to the reference genome. After assembling transcripts and quantifying gene expression levels, we detected 51 differentially expressed genes (DEGs) between HS and CN groups with thresholds of the adjusted p-value < 0.05 and |log2(FoldChange)| > 1. Among them, 33 and 18 genes were upregulated and downregulated, respectively. Gene ontology analyses further revealed that these DEGs were mainly associated with metabolism of lipids, thyroid hormone metabolic process, and cellular modified amino acid catabolic process. The upregulated ACACB, ACLY, LSS, and CYP7A1 genes were found to be inter-related through biological processes of thioester biosynthetic process, acyl-CoA biosynthetic process, acetyl-CoA metabolic process, and others. Six DEGs were further validated by quantitative real-time PCR analysis. The results revealed the candidate genes and biological processes that will potentially be considered as important regulatory factors involved in the heat stress response in rabbits.

Association between the IRS1 and FTO genes regulates body weight in rabbits.

Insulin receptor substrate (IRS) proteins play key roles in signal transduction in insulin and insulin-like growth factor signaling to control postnatal growth. The fat mass and obesity-associated (FTO) protein also play an essential role in postnatal growth. The aim of this study was to investigate the association between the IRS1 and FTO genes and the regulation of growth traits in rabbits. A total of nine synonymous SNPs were detected in the IRS1 coding sequence using direct sequencing, and the c.189G>T and c.2574G>A SNPs from two linkage disequilibrium blocks were further genotyped for association analysis in 216 New Zealand rabbits. The association results revealed that the TT genotype of c.189G>T and the AA genotype of c.2574G>A were significantly associated with higher body weight at 70 (BW70) and 84 (BW84) days of age and with higher average daily gain (P<0.05). Linear-regression analysis revealed that the two-gene combination model of FTO c.499G>A and IRS1 c.2574G>A was associated with BW70 and BW84. The combination model of the GA genotype of FTO c.499G>A with the AA genotype of IRS1 c.2574G>A was associated with preferred values for BW70 and BW84. The performance values for the FTO c.499G>A genotypes after stratification with regard to the IRS1 c.189G>T genotypes revealed that the TT genotype of IRS1 c.189G>T reduced the FTO c.499G>A significance associated with BW70 and BW84. Together, our data indicated that the IRS1 gene was associated with growth traits in rabbits. The IRS1 and FTO combination model may be exploited to assist breeders in selecting rabbits with preferred body weight.