Search

The process of skeletal muscle development is intricate and involves the regulation of a diverse array of genes. Accurate gene expression profiles are crucial for studying muscle development, making it essential to choose the right reference genes for real-time quantitative PCR (RT-qPCR). In the present study, eight candidate reference genes were identified from our previous transcriptome sequencing analysis of caprine skeletal muscle satellite cells (MuSCs), and two traditional reference genes (ACTB and GAPDH) were assessed. The quantitative levels of the candidate reference genes were determined through the RT-qPCR technique, while the stability of their expression was evaluated utilizing the GeNorm, NormFinder, BestKeeper, and RefFinder programs. Furthermore, the chosen reference genes were utilized for the normalization of the gene expression levels of PCNA and Myf5. It was determined that conventional reference genes, including ACTB and GAPDH, were not appropriate for normalizing target gene expression. Conversely, RPL14 and RPS15A, identified through RNA sequencing analysis, exhibited minimal variability and were identified as the optimal reference genes for normalizing gene expression during the proliferation and differentiation of goat MuSCs. Our research offers a validated panel of optimal reference genes for the detection of differentially expressed genes in goat muscle satellite cells using RT-qPCR.

Skeletal myogenesis is a complicated biological event that involves a succession of tightly controlled gene expressions. In order to identify novel regulators of this process, we performed mRNA-Seq studies of goat skeletal muscle satellite cells (MuSCs) cultured under proliferation (GM) and differentiation (DM1/DM5) conditions. A total of 19,871 goat genes were expressed during these stages, 198 of which represented novel transcripts. Notably, in pairwise comparisons at the different stages, 2551 differentially expressed genes (DEGs) were identified (p < 0.05), including 1560 in GM vs. DM1, 1597 in GM vs. DM5, and 959 in DM1 vs. DM5 DEGs. The time-series expression profile analysis clustered the DEGs into eight gene groups, three of which had significantly upregulated and downregulated patterns (p < 0.05). Functional enrichment analysis showed that DEGs were enriched for essential biological processes such as muscle structure development, muscle contraction, muscle cell development, striated muscle cell differentiation, and myofibril assembly, and were involved in pathways such as the MAPK, Wnt and PPAR signaling pathways. Moreover, the expression of eight DEGs (MYL2, DES, MYOG, FAP, PLK2, ADAM, WWC1, and PRDX1) was validated. These findings offer novel insights into the transcriptional regulation of skeletal myogenesis in goats.

As the largest chamber of the ruminant stomach, the rumen not only serves as the principal absorptive surface and nutrient transport pathway from the lumen into the animal, but also plays an important short-chain fatty acid (SCFA) metabolic role in addition to protective functions. Accurate characterization of the gene expression profiles of genes of interest is essential to the exploration of the intrinsic regulatory mechanisms of rumen development in goats. Thus, the selection of suitable reference genes (RGs) is an important prerequisite for real-time quantitative PCR (RT-qPCR). In the present study, 16 candidate RGs were identified from our previous transcriptome sequencing of caprine rumen tissues. The quantitative expressions of the candidate RGs were measured using the RT-qPCR method, and the expression stability of the RGs was assessed using the geNorm, NormFinder, and BestKeeper programs. GeNorm analysis showed that the M values were less than 0.5 for all the RGs except GAPT4 , indicating that they were stably expressed in the rumen tissues throughout development. RPS4X and RPS6 were the two most stable RGs. Furthermore, the expressions of two randomly selected target genes ( IGF 1 and TOP2A ), normalized by the selected most stable RGs ( RPS4X and RPS6 ), were consistent with the results of RNA sequencing, while the use of GAPDH and ACTB as RGs resulted in altered profiles. Overall, RPS4X and RPS6 showed the highest expression stability and the lowest coefficients of variation, and could be used as the optimal reference combination for quantifying gene expression in rumen tissues via RT-qPCR analysis.

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.

The rumen is an important digestive organ of ruminants. From the fetal to adult stage, the morphology, structure and function of the rumen change significantly. However, the knowledge of the intrinsic genetic regulation of these changes is still limited. We previously reported a genome-wide expression profile of miRNAs in pre-natal goat rumens. In this study, we combined and analyzed the transcriptomes of rumen miRNAs during pre-natal (E60 and E135) and post-natal (D30 and D150) stages. A total of 66 differentially expressed miRNAs (DEMs) were identified in the rumen tissues from D30 and D150 goats. Of these, 17 DEMs were consistently highly expressed in the rumens at the pre-weaning stages (E60, E135 and D30), while down-regulated at D150. Noteworthy, annotation analysis revealed that the target genes regulated by the DEMs were mainly enriched in MAPK signaling pathway, Jak-STAT signaling pathway and Ras signaling pathway. Interestingly, the expression of miR-148a-3p was significantly high in the embryonic stage and down-regulated at D150. The potential binding sites of miR-148a-3p in the 3'-UTR of QKI were predicted by the TargetScan and verified by the dual luciferase report assay. The co-localization of miR-148a-3p and QKI through in situ hybridization was observed in the rumen tissues but not in the intestinal tracts. Moreover, the expression of miR-148a-3p in the epithelium was significantly higher than that in the other layers of the rumen, suggesting that miR-148a-3p is involved in the development of the rumen epithelial cells by targeting QKI . Subsequently, miR-148a-3p inhibitor was found to induce the proliferation of GES-1 cells. Taken together, our study identified DEMs involved in the development of the rumen and provides insights into the regulation mechanism of rumen development in goats.

Brown adipose tissues have unique non-shivering thermogenesis functions, can be found in newborn ruminate animals, and then are gradually replaced by white adipose tissues in adulthood. For the purpose of exploring the intrinsic mechanism underlying the conversion process from brown (BAT) to white adipose tissue (WAT), it is necessary to utilize Quantitative PCR (qPCR) to study gene expression profiling. In this study, we identified reference genes that were consistently expressed during the transformation from goat BAT to WAT using RNA-seq data. Then, twelve genes were evaluated as candidate reference genes for qPCR in goat perirenal adipose tissue using three tools (geNorm, Normfinder, and BestKeeper). In addition, the selected reference genes were used to normalize the gene expression of PGC-1α and GPAT4 . It was found that traditional reference genes, such as GAPDH , RPLP0 , HPRT1 , and PPIA were not suitable for target gene normalization. In contrast, CTNNB , PFDN5 , and EIF3M , selected from RNA sequencing data, showed the least variation and were recommended as the best reference genes during the transformation from BAT to WAT.