Dynamics of transcriptome and chromatin accessibility revealed sequential regulation of potential transcription factors during the brown adipose tissue whitening in rabbits.

Brown adipose tissue (BAT) represents a valuable target for treating obesity in humans. BAT losses of thermogenic capacity and gains a "white adipose tissue-like (WAT-like)" phenotype (BAT whitening) under thermoneutral environments, which could lead to potential low therapy responsiveness in BAT-based obesity treatments. However, the epigenetic mechanisms of BAT whitening remain largely unknown. In this study, BATs were collected from rabbits at day0 (D0), D15, D85, and 2 years (Y2). RNA-sequencing (RNA-seq) and the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) were performed to investigate transcriptome and chromatin accessibility of BATs at the four whitening stages, respectively. Our data showed that many genes and chromatin accessible regions (refer to as "peaks") were identified as significantly changed during BAT whitening in rabbits. The BAT-selective genes downregulated while WAT-selective genes upregulated from D0 to Y2, and the de novo lipogenesis-related genes reached the highest expression levels at D85. Both the highly expressed genes and accessible regions in Y2 were significantly enriched in immune response-related signal pathways. Analysis of different relationships between peaks and their nearby genes found an increased proportion of the synchronous changes between chromatin accessibility and gene expression during BAT whitening. The synergistic changes between the chromatin accessibility of promoter and the gene expression were found in the key adipose genes. The upregulated genes which contained increased peaks were significantly enriched in the PI3K-Akt signaling pathway, steroid biosynthesis, TGF-beta signaling pathway, osteoclast differentiation, and dilated cardiomyopathy. Moreover, the footprinting analysis suggested that sequential regulation of potential transcription factors (TFs) mediated the loss of thermogenic phenotype and the gain of a WAT-like phenotype of BAT. In conclusion, our study provided the transcriptional and epigenetic frameworks for understanding BAT whitening in rabbits for the first time and might facilitate potential insights 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.

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.

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.

De novo transcriptome assembly and identification of genes associated with feed conversion ratio and breast muscle yield in domestic ducks.

Breast muscle yield and feed conversion efficiency are the major breeding aims in duck breeding. Understanding the role of specific transcripts in the muscle and small intestine might lead to the elucidation of interrelated biological processes. In this study, we obtained jejunum and breast muscle samples from two strains of Peking ducks that were sorted by feed conversion ratio (FCR) and breast muscle percentage into two-tailed populations. Ten RNA-Seq libraries were developed from the pooled samples and sequenced using the Hiseq2000 platform. We created a reference duck transcript database using de novo assembly methods, which included 16 663 irredundant contigs with an N50 length of 1530 bp. This new duck reference cDNA dataset significantly improved the mapping rate for RNA-Seq data, from 50% to 70%. Mapping and annotation were followed by Gene Ontology analysis, which showed that numerous genes were differentially expressed between the low and high FCR groups. The differentially expressed genes in the jejunum were enriched in biological processes related to immune response and immune response activation, whereas those in the breast muscle were significantly enriched in biological processes related to muscle cell differentiation and organ development. We identified new candidate genes, that is, PCK1, for improving the FCR and breast muscle yield of ducks and obtained much better reference duck transcripts. This study suggested that de novo assembly is essential when applying transcriptome analysis to a species with an incomplete genome.