Using RNA-Seq to Identify Reference Genes of the Transition from Brown to White Adipose Tissue in Goats

Identification and validation of stable reference genes for expression profiling of target genes in diverse ovine tissues

Quantitative PCR (qPCR) is a widely-used technique for quantifying the expression of target genes across various tissues, as well as under different pathological and physiological conditions. One of the challenges associated with this method is the need to identify optimal reference genes (RGs) that maintain consistent expression levels under diverse experimental settings, thereby ensuring accurate biological interpretation. In this study, we conducted a thorough analysis of 18 candidate RGs (ACTB, BACH1, B2M, GAPDH, HMBS, HPRT1, PGK1, PPIA, PPIB, RPLP0, RPL19, RPS9, RPS15, RPS28, SDHA, TBP, UXT, and YWHAZ) across 10 ovine tissues (muscle, skin, kidney, liver, intestine, rumen, lung, testis, heart, and spleen) obtained from five individual sheep. We aimed to identify genes with stable expression across these tissues. A literature-based survey helped us shortlist candidate genes representing various functional classes from multiple livestock species. We employed four algorithms: geNorm, NormFinder, BestKeeper, and Delta Ct (ΔCt), to rank these genes based on their stability. A consistent trend in the rankings was observed across these different algorithms. RefFinder was then used for a comprehensive ranking, integrating the outputs from the various methods. ACTB, PPIB, BACH1, and B2M emerged as the most stable RGs, while RPS9, RPS15, and PGK1 displayed variable expression. We validated our findings through qPCR analysis of four target genes (ACTN2, CRYAB, DLK1, and TRIM54) in the skin samples from two different sheep breeds. Based on these results, we recommend ACTB, PPIB, BACH1, and B2M as reliable internal control genes for qPCR experiments involving diverse ovine tissues.

Screening the optimal housekeeping genes (HKGs) of placenta tissues by RNA-sequence and qRT-PCR throughout gestation in goat (Capra Hircus)

Selection of stable housekeeping genes (HKGs) is very important for accurate calculation of relative expression levels of target genes by quantitative real-time polymerase chain reaction (qRT-PCR). At present, the appropriate HKGs have not been identified in placental tissues throughout the pregnancy of the goat. In our study, 20 HKGs were tentatively selected from RNA-seq data and previous reports. The cycle threshold (Ct) of HKGs was determined by qRT-PCR in trophoblast membrane and cotyledon villus collected from 38 Dazu Black goats on gestation days of 20, 25, 30, 45, 60, 90, 120, and 150 (birth). The expression stability of the HKGs was analyzed by geNorm, Normfinder, Bestkeeper and Delta Ct algorithms, and comprehensively evaluated by ReFinder and ComprFinder. In addition, the optimal HKGs were further verified by placenta-specific genes (SPP1, VEGFA and PAG6). The 16 candidate HKGs (except POP4, TBP, RNF10, UBC) showed a qualified Ct value, less than 28. Among them, YWHAZ, EIF3K and PPIB showed the most stable expression in placental tissues during early, mid-late pregnancy and postpartum, but the least stable expression was B2M at early and mid-late stage, and PPIB at postpartum. After comprehensive analysis, RPLP0, EIF3K and YWHAZ were found to be the most stable placental HKGs throughout pregnancy. The classical HKGs, ACTB, GAPDH and 18S RNA have unstable expressions and even ranked at the bottom of the list from comprehensive index, suggesting an inappropriate for target gene normalization. Taken together, our study confirmed that YWHAZ, EIF3K, HMBS and RPLP0 may be the optimal HKGs in goat placenta at different stage of pregnancy, which provided a valuable reference of HKGs on functional gene expression detection for further research on placenta development and growth in ruminants.

Implications of MCU complex in metabolic diseases

Metabolic diseases are considered the primary culprit for physical and mental health of individuals. Although the diagnosis of these diseases is relatively easy, more effective and convenient potent drugs are still being explored. Ca2+ across the inner mitochondrial membrane is a vital intracellular messenger that regulates energy metabolism and cellular Ca2+ homeostasis and is involved in cell death. Mitochondria rely on a selective mitochondrial Ca2+ unidirectional transport complex (MCU complex) in their inner membrane for Ca2+ uptake. We found that the channel contains several subunits and undergoes dramatic transformations in various pathological processes, especially in metabolic diseases. In this way, we believe that the MCU complex becomes a target with significant potential for these diseases. However, there is no review linking the two factors, thus hindering the possibility of new drug production. Here, we highlight the connection between MCU complex-related Ca2+ transport and the pathophysiology of metabolic diseases, adding understanding and insight at the molecular level to provide new insights for targeting MCU to reverse metabolism-related diseases.

Transcriptomic Analysis Reveals Fibroblast Growth Factor 11 (FGF11) Role in Brown Adipocytes in Thermogenic Regulation of Goats

Brown adipose tissue (BAT) is the main site of adaptive thermogenesis, generates heat to maintain body temperature upon cold exposure, and protects against obesity by promoting energy expenditure. RNA-seq analysis revealed that FGF11 is enriched in BAT. However, the functions and regulatory mechanisms of FGF11 in BAT thermogenesis are still limited. In this study, we found that FGF11 was significantly enriched in goat BAT compared with white adipose tissue (WAT). Gain- and loss-of-function experiments revealed that FGF11 promoted differentiation and thermogenesis in brown adipocytes. However, FGF11 had no effect on white adipocyte differentiation. Furthermore, FGF11 promoted the expression of the UCP1 protein and an EBF2 element was responsible for UCP1 promoter activity. Additionally, FGF11 induced UCP1 gene expression through promoting EBF2 binding to the UCP1 promoter. These results revealed that FGF11 promotes differentiation and thermogenesis in brown adipocytes but not in white adipocytes of goats. These findings provide evidence for FGF11 and transcription factor regulatory functions in controlling brown adipose thermogenesis of goats.

Screening of Candidate Housekeeping Genes in Uterus Caruncle by RNA-Sequence and qPCR Analyses in Different Stages of Goat (Capra hircus)

The uterus is a critical pregnancy organ for mammals. The normal growth and development of ruminant uterus caruncles are crucial to maintain gestation and fetal health in goats. Quantitative real-time polymerase chain reaction (qRT-PCR) is a reliable tool to study gene expression profiling for exploring the intrinsic mechanism underlying the conversion process of uterus caruncle tissue. However, the candidate housekeeping genes (HKGs) are required for normalizing the expression of function genes. In our study, 22 HKGs were selected from analyzing transcriptome data at non-pregnancy and pregnancy processes and previous reports about HKGs in goat tissues. We assessed them for expression suitability in 24 samples from uterus tissues at 15 non-pregnant days (Stage 1), early (Stage 2), and medium-later pregnant days (Stage 3). The expression stability of these genes was evaluated by using geNorm, Normfinder, Bestkeeper, and Delta Ct algorithms and, comprehensively, by ReFinder. In addition, the most and least stable HKGs were used to normalize the target genes expression of SPP1, VEGFA, and PAG8. It was found that traditional reference genes, such as ACTB and GAPDH, were not suitable for target gene normalization. In contrast, PPIB selected from RNA sequencing data and EIF3K selected from previous references showed the least variation and were recommended as the best HKGs during the nonpregnant stage and the whole stages of goat uterus caruncle tissue, respectively. It is the first time the HKGs genes in uterus during the non-pregnant day and throughout the total pregnancy have been explored. These findings found suitable HKGs in uterus caruncle tissues at various stages of non-pregnancy and pregnancy; these can be useful for gene expression studies to reveal the molecular mechanisms of uterus development in goats.

A Novel LncRNA MSTRG.310246.1 Promotes Differentiation and Thermogenesis in Goat Brown Adipocytes

Brown adipose tissue (BAT) plays a critical role in maintaining the body temperature in newborn lamb due to its unique non-shivering thermogenesis. Previous studies have found that BAT thermogenesis is regulated by several long non-coding RNAs (lncRNAs). Here, we identified a novel lncRNA, MSTRG.310246.1, which was enriched in BAT. MSTRG.310246.1 was localized in both the nuclear and cytoplasmic compartments. In addition, MSTRG.310246.1 expression was upregulated during brown adipocyte differentiation. Overexpression of MSTRG.310246.1 increased the differentiation and thermogenesis of goat brown adipocytes. On the contrary, the knockdown of MSTRG.310246.1 inhibited the differentiation and thermogenesis of goat brown adipocytes. However, MSTRG.310246.1 had no effect on goat white adipocyte differentiation and thermogenesis. Our results show that MSTRG.310246.1 is a BAT-enriched LncRNA that improves the differentiation and thermogenesis of goat brown adipocytes.

Transcriptome-based selection and validation of optimal reference genes in perirenal adipose developing of goat (Capra hircus)

Brown adipose tissue (BAT) is mainly present in young mammals and is important for maintaining body temperature in neonatal mammals because of its ability to produce non-shivering thermogenesis. There is usually a large amount of BAT around the kidneys of newborn kids, but the BAT gradually "whiting" after birth. Screening and validating appropriate reference genes is a prerequisite for further studying the mechanism of goat brown adipose tissue "whiting" during the early stages. In this study, the expression stability of 17 candidate reference genes: 12 COPS8, SAP18, IGF2R, PARL, SNRNP200, ACTG1, CLTA, GANAB, GABARAP, PCBP2, CTSB, and CD151) selected based on previous transcriptome data as new candidate reference genes, 3 (PFDN5, CTNNB1, and EIF3M) recommended in previous studies, and 2 traditional reference genes (ACTB and GAPDH) was evaluated. Real-time quantitative PCR (RT-qPCR) technology was used to detect the expression level of candidate reference genes during goat BAT "whiting". Four algorithms: Normfinder, geNorm, ΔCt method, and BestKeeper, and two comprehensive algorithms: ComprFinder and RefFinder, were used to analyze the stability of each candidate reference genes. GABARAP, CLTA, GAPDH, and ACTB were identified as the most stable reference genes, while CTNNB1, CTSB, and EIF3M were the least stable. Moreover, two randomly selected target genes IDH2 and RBP4, were effectively normalized using the selected most stable reference genes. These findings collectively suggest that GABARAP, CLTA, GAPDH, and ACTB are relatively stable reference genes that can potentially be used for the development of perirenal fat in goats.

Maternal L-carnitine supplementation promotes brown adipose tissue thermogenesis of newborn goats after cold exposure

Brown adipose tissue (BAT) is an important component of energy expenditure and necessary to maintain body temperature for newborn mammals. In the previous study, we found that L-carnitine was enriched in BAT and promoted BAT adipogenesis and thermogenesis in goat brown adipocytes. However, whether dietary L-carnitine regulates BAT heat production and energy expenditure in lambs remains unclear. In this study, maternal L-carnitine supplementation elevated the rectal temperature, as well as the expression of UCP1 and mitochondrial DNA content to promote BAT thermogenesis in newborn goats. Moreover, maternal L-carnitine supplementation increased the levels of triglycerides (TG), non-esterified fatty acids (NEFA), and lactate in plasma, as well as the content of lipid droplet and glycogen in BAT of newborn goats. Lipidomic analysis showed that maternal L-carnitine supplementation remodeled the lipid composition of BAT in newborn goats. L-carnitine significantly increased the levels of TG and diglyceride (DG) and decreased the levels of glycerophospholipids and sphingolipids in BAT. Further studies showed that L-carnitine promoted TG and glycogen deposition in brown adipocytes through AMPKα. Our results indicate that maternal L-carnitine supplementation promotes BAT development and thermogenesis in newborn goats and provides new evidence for newborn goats to maintain body temperature in response to cold exposure.

Analysis of stability of reference genes for qPCR in bovine preadipocytes during proliferation and differentiation in vitro

The stability of internal reference genes is crucial to the reliability of gene expression results using real-time fluorescence quantitative PCR (qRT-PCR). Inappropriate reference genes may lead to inaccurate results or even wrong conclusions. This study aims to identify stable reference genes for analyzing the expression of proliferation-related and differentiation-inducing genes in bovine primary preadipocytes (BPPs) in vitro. In this study, the stability of 16 candidate internal reference genes (GAPDH, ACTB, PPIA, LRP10, HPRT1, YWHAZ, B2M, TBP, EIF3K, RPS9, UXT, 18S rRNA, RPLP0, MARVELD, EMD and RPS15A) for qRT-PCR at proliferation and differentiation stages of BPPs was investigated by three different algorithms (geNorm, NormFinder and BestKeeper). The expression of two marker genes, PCNA and LPL, was used to determine the validity of the candidate reference genes (RGs) at the proliferation and differentiation stages, respectively. The results showed that GAPDH and RPS15A were the most stable RGs in the proliferation of bovine primary preadipocyte, while PPIA was the least stable internal reference gene. RPLP0 and EIF3K were the most stable RGs in the differentiation induction of bovine primary preadipocyte, while GAPDH was the least stable internal reference gene. This study of RGs laid the foundation for subsequent research into the mechanism of proliferation and differentiation of BPPs in vitro using qRT-PCR.

RNA-Seq reveals miRNA role in thermogenic regulation in brown adipose tissues of goats

BACKGROUND

MicroRNAs (miRNAs) are a family of short non-coding RNA molecules and play important roles in various biological processes. However, knowledge of the expression profiles and function of miRNAs on the regulation of brown adipose tissue (BAT) thermogenesis remains largely unknown.

RESULTS

In this study, we found that brown adipose tissue (BAT) existed within the perirenal fat at 1 day after birth (D1) and transferred into white adipose tissue (WAT) at 30 days after birth (D30) by UCP1 protein expression and immunohistochemistry analysis. After that, we performed RNA sequencing on six libraries of goat BAT and WAT. A total of 238 known miRNAs and 1834 goat novel miRNAs were identified. Moreover, 395 differentially expressed miRNAs including 167 up-regulated and 228 down-regulated miRNAs were obtained in BAT. For the known BAT enriched miRNA, 30 miRNAs were enriched in goat BAT but not in mouse BAT. In addition, miR-433 was enriched in goat BAT but not in mouse BAT. Gain- and loss-of-function experiments reveal that miR-433 reduced the lipid accumulation of brown adipocytes and decreased the expression of BAT marker and mitochondrial related genes. However, miR-433 had no effect on lipid accumulation and thermogenesis in white adipocytes. In addition, miR-433 inhibited the expression of MAPK8 by targeting to the 3'UTR of MAPK8 gene. These data demonstrate that miR-433 acts as a negative regulator in controlling brown adipocytes differentiation and thermogenesis.

CONCLUSION

The present study provides a detailed miRNAs expression landscape in BAT and WAT. Furthermore, we found that miR-433, which was highly expressed on BAT had a negative regulatory function on the thermogenesis and adipogenesis in goat brown adipocytes. This study provides evidence for understanding the role of miRNAs in regulating BAT thermogenesis and energy expenditure in goats.

Integrated application of metabolomics and RNA-seq reveals thermogenic regulation in goat brown adipose tissues

Brown adipose tissue (BAT) plays an important role on no shivering thermogenesis during cold exposure to maintain animal body temperature and energy homeostasis. However, knowledge of the cellular transition from white adipose tissue (WAT) to BAT is still limited. In this study, we provided a comprehensive metabolomics and transcriptional signatures of goat BAT and WAT. A total of 157 metabolites were significantly changed, including 81 upregulated and 76 downregulated metabolites. In addition, we identified the citric acid cycle, fatty acid elongation, and degradation pathways as coordinately activated in BAT. Interestingly, five unsaturated fatty acids (Eicosadienoic Acid, C20:2; γ-Linolenic acid, C20:3; Arachidonic Acid, C20:4; Adrenic acid, C22:4; Docosahexaenoic acid, C22:6), Succinate, L-carnitine, and L-palmitoyl-carnitine were found to be abundant in BAT. Furthermore, L-carnitine, an intermediate of fatty acid degradation, is required for goat brown adipocyte differentiation and thermogenesis through activating AMPK pathway. However, L-carnitine decreased lipid accumulation through inducing lipolysis and thermogenesis in white adipocytes. These results revealed that there are the significant alterations in transcriptomic and metabolomic profiles between goat WAT and BAT, which may contribute to better understanding the roles of metabolites in BAT thermogenesis process.

Dynamic Expression Profiles of Circular RNAs during Brown to White Adipose Tissue Transformation in Goats (Capra hircus)

Simple Summary

In our study, we launched RNA-seq in order to investigate the potential functions of circRNA during brown adipose tissue (BAT) to white adipose tissue (WAT) transformation. As a result, 6610 circRNAs and 61 differentially expressed circRNAs (DEcircRNAs) were identified. Moreover, 65 miRNAs were detected that could potentially interact with DEcircRNAs. The present study provides a detailed circRNA expression landscape and evidence for circRNA functions in the transformation from BAT to WAT.

Abstract

Adipose tissues are mainly divided into brown adipose tissue (BAT) and white adipose tissue (WAT). WAT mainly functions to buffer excess calories, whereas BAT plays a role in the non-shivering thermogenesis to maintain body temperature and energy balance. Moreover, circRNAs play important roles in various biological processes. However, knowledge of the expression profile and function of circRNAs from BAT to WAT remains largely unknown. In this study, a total of 6610 unique circRNAs were identified in the perirenal adipose tissues of 1-day, 30-days, and 1-year goats. Functional annotation revealed that host genes of circRNAs were involved in some BAT-related pathways, such as the thyroid hormone signaling pathway, MAPK signaling pathway, and VEGF signaling pathway. Furthermore, a total of 61 DEcircRNAs were detected across three stages. Additionally, five selected circRNAs were validated by RNase R assay, qPCR, and Sanger sequencing. Finally, the circRNA–miRNA network was constructed between the DEcircRNAs and their miRNA binding sites.