Comparative Transcriptome Profiling of mRNA and lncRNA Related to Tail Adipose Tissues of Sheep

Non-Coding RNAs in Regulating Fat Deposition in Farm Animals

Fat deposition represents a crucial feature in the expenditure of physical energy and affects the meat quality of farm animals. It is regulated by multiple genes and regulators. Of them, non-coding RNAs (ncRNAs) play a critical role in modulating the fat deposition process. As well as being an important protein source, farm animals can be used as medical models, so many researchers worldwide have explored their mechanism of fat deposition. This article summarizes the transcription factors, regulatory genes, and signaling pathways involved in the molecular regulation process of fat deposition; outlines the progress of researching the roles of microRNAs (miRNAs), long ncRNAs (lncRNAs), and circular RNAs (circRNAs) in fat deposition in common farm animals including pigs, cattle, sheep, ducks, and chickens; and identifies scientific problems in the field that must be further investigated. It has been demonstrated that ncRNAs play a critical role in regulating the fat deposition process and have great potential in improving meat quality traits.

Understanding lncRNAs: key regulators of myogenesis and lipogenesis in farm animals

Long non-coding RNAs (lncRNAs) are RNA molecules exceeding 200 nucleotides in length. Recent studies have demonstrated their involvement in regulating gene expression and various biological processes. Among these, myogenesis and lipogenesis are particularly important because of their direct effects on muscle development and fat deposition in farm animals. These processes are crucial for determining meat quality, growth rates, and overall economic value in animal husbandry. Although the specific mechanisms through which lncRNAs influence these pathways are still under investigation, further research into their roles in muscle and fat development is crucial for optimizing farm animal breeding strategies. Here, we review the characteristics of lncRNAs, including their biogenesis, localization, and structures, with a particular focus on their association with myogenesis and adipogenesis. This review seeks to establish a theoretical foundation for enhancing farm animal production. In particular, focusing on lncRNAs may reveal how these molecules can enhance the economic traits of farm animals, thereby contributing to the optimization of farm animal breeding processes.

Ovine LncRSFD1 Mined from RNA-Seq: Identification, Expression Profile, Promotion of Preadipocyte Differentiation, Promoter Activity, and Its Polymorphisms Related to Phenotypic Traits

Tail fat is essential for sheep survival in extreme environments, yet its significance is often overlooked, leading to the decline of fat-tailed breeds. This study identified a novel lncRNA, lncRSFD1 (TCONS_00054953), through transcriptome sequencing, showing differential expression in the tail adipose tissues of Lanzhou Fat-Tailed (LFT) sheep and Tibetan (TS) sheep. Highly expressed in adipose tissues, lncRSFD1 inhibits preadipocyte proliferation and promotes 3T3-L1 differentiation, suggesting its role in regulating fat deposition. Located in both the cytoplasm and nucleus, lncRSFD1 targets the neighboring gene PDE4DIP and may function as a molecular sponge for conserved miRNAs, including oar-miR-30a-3p, oar-miR-329b-5p, and oar-miR-431, which are known to influence fat and muscle-related physiological processes. Moreover, the core promoter of lncRSFD1 (-2607 bp to -1776 bp) harbors four SNPs (g.-2429G>A, g.-2030T>C, g.-2016C>T, g.-2015G>A) significantly associated with growth traits such as body height in Guiqian Semi-Fine Wool (GSFW) sheep. These findings suggest lncRSFD1 plays a key role in fat deposition and growth regulation, offering new insights into the molecular mechanisms of lncRNAs in sheep. It provides a potential target for genetic improvement and molecular breeding to enhance fat deposition and adaptability in sheep breeds.

Lipidomics analysis of adipose depots at differently aged Sunit sheep

The objective of this study was to optimize the utilization of deposited fat in Sunit sheep, with a focus on dietary nutrition. This study also elucidated variations in lipid metabolism among subcutaneous fat (SF), perirenal fat (PF), and tail fat (TF) in sheep of different ages using non-targeted lipidomic techniques. In total, 173 different lipids were identified, of which triacylglycerol (TG) and phosphatidylcholine (PC) were prominent. The relative intensity of TG was highest at 6 months of age in three adipose depots. Glycerophospholipids (PLs) were expressed at peak levels in TF and SF at 18 months of age. Pathway analysis revealed that biosynthesis of unsaturated fatty acids, linoleic acid metabolism, glycerophospholipid metabolism, and fatty acid biosynthesis were the main pathways involved in the metabolism of adipose depots. These findings provide a comprehensive reference for the metabolic characteristics and pathways of adipose tissue in sheep and the utilization of its by-products.

A meta-analysis approach for annotation and identification of lncRNAs controlling perirenal fat deposition in suckling lambs

Long non-coding RNAs (lncRNAs) are being studied in farm animals due to their association with traits of economic interest, such as fat deposition. Based on the analysis of perirenal fat transcriptomes, this research explored the relevance of these regulatory elements to fat deposition in suckling lambs. To that end, meta-analysis techniques have been implemented to efficiently characterize and detect differentially expressed transcripts from two different RNA-seq datasets, one including samples of two sheep breeds that differ in fat deposition features, Churra and Assaf (n = 14), and one generated from Assaf suckling lambs with different fat deposition levels (n = 8). The joint analysis of the 22 perirenal fat RNA-seq samples with the FEELnc software allowed the detection of 3953 novel lncRNAs. After the meta-analysis, 251 differentially expressed genes were identified, 21 of which were novel lncRNAs. Additionally, a co-expression analysis revealed that, in suckling lambs, lncRNAs may play a role in controlling angiogenesis and thermogenesis, processes highlighted in relation to high and low fat deposition levels, respectively. Overall, while providing information that could be applied for the improvement of suckling lamb carcass traits, this study offers insights into the biology of perirenal fat deposition regulation in mammals.

Structural variant landscapes reveal convergent signatures of evolution in sheep and goats

BACKGROUND

Sheep and goats have undergone domestication and improvement to produce similar phenotypes, which have been greatly impacted by structural variants (SVs). Here, we report a high-quality chromosome-level reference genome of Asiatic mouflon, and implement a comprehensive analysis of SVs in 897 genomes of worldwide wild and domestic populations of sheep and goats to reveal genetic signatures underlying convergent evolution.

RESULTS

We characterize the SV landscapes in terms of genetic diversity, chromosomal distribution and their links with genes, QTLs and transposable elements, and examine their impacts on regulatory elements. We identify several novel SVs and annotate corresponding genes (e.g., BMPR1B, BMPR2, RALYL, COL21A1, and LRP1B) associated with important production traits such as fertility, meat and milk production, and wool/hair fineness. We detect signatures of selection involving the parallel evolution of orthologous SV-associated genes during domestication, local environmental adaptation, and improvement. In particular, we find that fecundity traits experienced convergent selection targeting the gene BMPR1B, with the DEL00067921 deletion explaining ~10.4% of the phenotypic variation observed in goats.

CONCLUSIONS

Our results provide new insights into the convergent evolution of SVs and serve as a rich resource for the future improvement of sheep, goats, and related livestock.

Deciphering the role of alternative splicing as a potential regulator in fat-tail development of sheep: a comprehensive RNA-seq based study

Although research on alternative splicing (AS) has been widely conducted in mammals, no study has investigated the splicing profiles of genes involved in fat-tail formation in sheep. Here, for the first time, a comprehensive study was designed to investigate the profile of AS events and their involvement in fat-tail development of sheep. In total, 45 RNA-Seq samples related to seven different studies, which have compared the fat-tailed vs thin-tailed sheep breeds, were analyzed. Two independent tools, rMATS and Whippet, along with a set of stringent filters were applied to identify differential AS (DAS) events between the breeds per each study. Only DAS events that were detected by both tools as well as in at least three datasets with the same ΔPSI trend (percent spliced in), were considered as the final high-confidence set of DAS genes. Final results revealed 130 DAS skipped exon events (69 negative and 61 positive ΔPSI) belonged to 124 genes. Functional enrichment analysis highlighted the importance of the genes in the underlying molecular mechanisms of fat metabolism. Moreover, protein-protein interaction network analysis revealed that DAS genes are significantly connected. Of DAS genes, five transcription factors were found that were enriched in the biological process associated with lipid metabolism like "Fat Cell Differentiation". Further investigations of the findings along with a comprehensive literature review provided a reliable list of candidate genes that may potentially contribute to fat-tail formation including HSD11B1, SIRT2, STRN3 and TCF7L2. Based on the results, it can be stated that the AS patterns may have evolved, during the evolution of sheep breeds, as another layer of regulation to contribute to biological complexity by reprogramming the gene regulatory networks. This study provided the theoretical basis of the molecular mechanisms behind the sheep fat-tail development in terms of AS.

A Single-Cell Atlas of an Early Mongolian Sheep Embryo

Cell types have been established during organogenesis based on early mouse embryos. However, our understanding of cell types and molecular mechanisms in the early embryo development of Mongolian sheep has been hampered. This study presents the first comprehensive single-cell transcriptomic characterization at E16 in Ujumqin sheep and Hulunbuir short-tailed sheep. Thirteen major cell types were identified at E16 in Ujumqin sheep, and eight major cell types were identified at E16 in Hulunbuir short-tailed sheep. Function enrichment analysis showed that several pathways were significantly enriched in the TGF-beta signaling pathway, the Hippo signaling pathway, the platelet activation pathway, the riboflavin metabolism pathway, the Wnt signaling pathway, regulation of the actin cytoskeleton, and the insulin signaling pathway in the notochord cluster. Glutathione metabolism, glyoxylate, and dicarboxylate metabolism, the citrate cycle, thyroid hormone synthesis, pyruvate metabolism, cysteine and methionine metabolism, thermogenesis, and the VEGF signaling pathway were significantly enriched in the spinal cord cluster. Steroid biosynthesis, riboflavin metabolism, the cell cycle, the Hippo signaling pathway, the Hedgehog signaling pathway, the FoxO signaling pathway, the JAK-STAT signaling pathway, and the Wnt signaling pathway were significantly enriched in the paraxial mesoderm cluster. The notochord cluster, spinal cord cluster, and paraxial mesoderm cluster were found to be highly associated with tail development. Pseudo-time analysis demonstrated that the mesenchyme can translate to the notochord in Ujumqin sheep. Molecular assays revealed that the Hippo signaling pathway was enriched in Ujumqin sheep. This comprehensive single-cell map revealed previously unrecognized signaling pathways that will further our understanding of the mechanism of short-tailed sheep formation.

Meta-analysis of RNA-Seq datasets highlights novel genes/pathways involved in fat deposition in fat-tail of sheep

Introduction

Fat-tail in sheep is considered as an important energy reservoir to provide energy as a survival buffer during harsh challenges. However, fat-tail is losing its importance in modern sheep industry systems and thin-tailed breeds are more desirable. Using comparative transcriptome analysis to compare fat-tail tissue between fat- and thin-tailed sheep breeds provides a valuable approach to study the complex genetic factors associated with fat-tail development. However, transcriptomic studies often suffer from issues with reproducibility, which can be improved by integrating multiple studies based on a meta-analysis.

Methods

Hence, for the first time, an RNA-Seq meta-analysis on sheep fat-tail transcriptomes was performed using six publicly available datasets.

Results and discussion

A total of 500 genes (221 up-regulated, 279 down-regulated) were identified as differentially expressed genes (DEGs). A jackknife sensitivity analysis confirmed the robustness of the DEGs. Moreover, QTL and functional enrichment analysis reinforced the importance of the DEGs in the underlying molecular mechanisms of fat deposition. Protein-protein interactions (PPIs) network analysis revealed the functional interactions among the DEGs and the subsequent sub-network analysis led to identify six functional sub-networks. According to the results of the network analysis, down-regulated DEGs in green and pink sub-networks (like collagen subunits IV, V, and VI, integrins 1 and 2, SCD, SCD5, ELOVL6, ACLY, SLC27A2, and LPIN1) may impair lipolysis or fatty acid oxidation and cause fat accumulation in tail. On the other hand, up-regulated DEGs, especially those are presented in green and pink sub-networks (like IL6, RBP4, LEPR, PAI-1, EPHX1, HSD11B1, and FMO2), might contribute to a network controlling fat accumulation in the tail of sheep breed through mediating adipogenesis and fatty acid biosynthesis. Our results highlighted a set of known and novel genes/pathways associated with fat-tail development, which could improve the understanding of molecular mechanisms behind fat deposition in sheep fat-tail.

MiRNA-Seq reveals key MicroRNAs involved in fat metabolism of sheep liver

There is a genetic difference between Hu sheep (short/fat-tailed sheep) and Tibetan sheep (short/thin-tailed sheep) in tail type, because of fat metabolism. Previous studies have mainly focused directly on sheep tail fat, which is not the main organ of fat metabolism. The function of miRNAs in sheep liver fat metabolism has not been thoroughly elucidated. In this study, miRNA-Seq was used to identify miRNAs in the liver tissue of three Hu sheep (short/fat-tailed sheep) and three Tibetan sheep (short/thin-tailed sheep) to characterize the differences in fat metabolism of sheep. In our study, Hu sheep was in a control group, we identified 11 differentially expressed miRNAs (DE miRNAs), including six up-regulated miRNAs and five down-regulated miRNAs. Miranda and RNAhybrid were used to predict the target genes of DE miRNAs, obtaining 3,404 target genes. A total of 115 and 67 GO terms as well as 54 and 5 KEGG pathways were significantly (padj < 0.05) enriched for predicted 3,109 target genes of up-regulated and 295 target genes of down-regulated miRNAs, respectively. oar-miR-432 was one of the most up-regulated miRNAs between Hu sheep and Tibetan sheep. And SIRT1 is one of the potential target genes of oar-miR-432. Furthermore, functional validation using the dual-luciferase reporter assay indicated that the up-regulated miRNA; oar-miR-432 potentially targeted sirtuin 1 (SIRT1) expression. Then, the oar-miR-432 mimic transfected into preadipocytes resulted in inhibited expression of SIRT1. This is the first time reported that the expression of SIRT1 gene was regulated by oar-miR-432 in fat metabolism of sheep liver. These results could provide a meaningful theoretical basis for studying the fat metabolism of sheep.

Whole-body adipose tissue multi-omic analyses in sheep reveal molecular mechanisms underlying local adaptation to extreme environments

The fat tail of sheep is an important organ that has evolved to adapt to extreme environments. However, the genetic mechanisms underlying the fat tail phenotype remain poorly understood. Here, we characterize transcriptome and lipidome profiles and morphological changes in 250 adipose tissues from two thin-tailed and three fat-tailed sheep populations in summer and winter. We implement whole-genome selective sweep tests to identify genetic variants related to fat-tails. We identify a set of functional genes that show differential expression in the tail fat of fat-tailed and thin-tailed sheep in summer and winter. These genes are significantly enriched in pathways, such as lipid metabolism, extracellular matrix (ECM) remodeling, molecular transport, and inflammatory response. In contrast to thin-tailed sheep, tail fat from fat-tailed sheep show slighter changes in adipocyte size, ECM remodeling, and lipid metabolism, and had less inflammation in response to seasonal changes, indicating improved homeostasis. Whole-genome selective sweep tests identify genes involved in preadipocyte commitment (e.g., BMP2, PDGFD) and terminal adipogenic differentiation (e.g., VEGFA), which could contribute to enhanced adipocyte hyperplasia. Altogether, we establish a model of regulatory networks regulating adipose homeostasis in sheep tails. These findings improve our understanding of how adipose homeostasis is maintained, in response to extreme environments in animals.

Analyses of lncRNAs, circRNAs, and the Interactions between ncRNAs and mRNAs in Goat Submandibular Glands Reveal Their Potential Function in Immune Regulation

As part of one of the main ruminants, goat salivary glands hardly secrete digestive enzymes, but play an important role in immunity. The immune function of goat salivary glands significantly changes with age, while the expression profile and specific function of non-coding RNA during this process are unknown. In this study, transcriptome sequencing was performed on submandibular gland (SMG) tissues of 1-month-old, 12-month-old, and 24-month-old goats, revealing the expression patterns of lncRNA and circRNA at different ages. A total of 369 lncRNAs and 1699 circRNAs were found to be differentially expressed. Functional enrichment analyses showed that the lncRNA regulated target mRNAs and circRNA host genes were significantly enriched in immune-related GO terms and pathways. CeRNA network analysis showed that the key differentially expressed circRNAs and lncRNAs mainly regulate the key immune-related genes ITGB2, LCP2, PTPRC, SYK, and ZAP70 through competitive binding with miR-141-x, miR-29-y, and chi-miR-29b-3p, thereby affecting the natural killer cell-mediated cytotoxicity pathway, the T cell receptor signaling pathway, and other immune-related pathways. It should be noted that the expression of key circRNAs, lncRNAs, and key immune-related genes in goat SMGs decreased significantly with the growth of the goat. This is the first reporting of lncRNAs, circRNAs, and ceRNA network regulation in goat SMGs. Our study contributes to the knowledge of changes in the expression of non-coding RNAs during SMG development in goats and provides new insights into the relationship between non-coding RNAs and salivary gland immune function in goats.

The landscape of the long non-coding RNAs and circular RNAs of the abdominal fat tissues in the chicken lines divergently selected for fatness

BACKGROUND

Excessive deposition of abdominal fat poses serious problems in broilers owing to rapid growth. Recently, the evolution of the existing knowledge on long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) have established their indispensable roles in multiple physiological metabolic processes, including adipogenesis and fat deposition. However, not much has been explored on their profiles in the abdominal fat tissues of broilers to date. In the study, we aimed to characterize the vital candidates of lncRNAs and circRNAs and their underlying regulations for abdominal fat deposition in broilers.

RESULTS

The present study sequenced the lncRNAs and circRNAs expression profiles in the abdominal fat tissues isolated from 7-week-old broilers, who were divergently selected for their fatness. It identified a total of 3359 lncRNAs and 176 circRNAs, demonstrating differential expressed (DE) 30 lncRNAs and 17 circRNAs between the fat- and lean-line broilers (|log2FC| ≥ 1, P < 0.05). Subsequently, the 20 cis-targets and 48 trans-targets of the candidate DE lncRNAs were identified for depositing abdominal fat by adjacent gene analysis and co-expression analysis, respectively. In addition, the functional enrichment analysis showed the DE lncRNAs targets and DE circRNAs host genes to be mainly involved in the cellular processes, amino/fatty acid metabolism, and immune inflammation-related pathways and GO terms. Finally, the vital 16 DE lncRNAs located in cytoplasm and specifically expressed in fat/lean line and their targets were used to construct the lncRNA-miRNA-mRNA competing endogenous RNA (ceRNA) regulatory network, comprising 7 DE lncRNAs, 28 miRNAs, 11 DE mRNAs. Notably, three lncRNAs including XR_001468036.2, XR_003077610.1 and XR_001466431.2 with the most connected degrees might play hub regulatory roles in abdominal fat deposition of broilers.

CONCLUSIONS

This study characterized the whole expression difference of lncRNAs and circRNAs between the two lines broilers with divergently ability of abdominal fat. The vital candidate DE lncRNAs/circRNAs and ceRNA regulations were identified related to the deposition of abdominal fat in chicken. These results might further improve our understanding of regulating the non-coding RNAs in obesity.

A 7-nt nucleotide sequence variant within the sheep KDM3B gene affects female reproduction traits

Lysine demethylase 3B (KDM3B) gene is a histone demethylase, demonstrating specific demethylation of the histone H3 lysine 9. It was detected as a sheep reproductive candidate gene by genome-wide scans, and related studies also showed its significance in female reproductive process. However, rare study researched its polymorphism. Herein, we hypothesized that the polymorphisms of KDM3B gene were associated with sheep reproduction traits. A 7-nt nucleotide sequence variant (rs1088697156) within KDM3B gene was identified in a total of 888 individuals, including the Australian White (AUW) sheep and Lanzhou Fat-tailed (LFT) sheep. II (insertion/insertion) and ID (insertion/deletion) genotypes of 7-nt variant were detected, which were at Hardy-Weinberg equilibrium (HWE) in detected breeds. Association analysis illustrated the 7-nt variant was significantly associated with the litter size, duration of pregnancy, live lamb number, live lamb rate, stillbirth number, stillbirth rate of average and different parity (P < 0.05) in AUW sheep. Moreover, 'ID' was the dominant genotype with excellent consistency in reproductive traits. It is instrumental to select individuals with ID genotype for improving the sheep reproduction traits. These findings suggest that the 7-nt variant within KDM3B gene can be used as a candidate marker of reproduction traits for sheep breeding improvement by marker-assisted selection.

MicroRNA and circular RNA profiling in the deposited fat tissue of Sunite sheep

As the most typical deposited fat, tail fat is an important energy reservoir for sheep adapted to harsh environments and plays an important role as a raw material in daily life. However, the regulatory mechanisms of microRNA (miRNA) and circular RNA (circRNA) in tail fat development remain unclear. In this study, we characterized the miRNA and circRNA expression profiles in the tail fat of sheep at the ages of 6, 18, and 30 months. We identified 219 differentially expressed (DE) miRNAs (including 12 novel miRNAs), which exhibited a major tendency to be downregulated, and 198 DE circRNAs, which exhibited a tendency to be upregulated. Target gene prediction analysis was performed for the DE miRNAs. Functional analysis revealed that their target genes were mainly involved in cellular interactions, while the host genes of DE circRNAs were implicated in lipid and fatty acid metabolism. Subsequently, we established a competing endogenous RNA (ceRNA) network based on the negative regulatory relationship between miRNAs and target genes. The network revealed that upregulated miRNAs play a leading role in the development of tail fat. Finally, the ceRNA relationship network with oar-miR-27a_R-1 and oar-miR-29a as the core was validated, suggesting possible involvement of these interactions in tail fat development. In summary, DE miRNAs were negatively correlated with DE circRNAs during sheep tail fat development. The multiple ceRNA regulatory network dominated by upregulated DE miRNAs may play a key role in this developmental process.

Genetics of the phenotypic evolution in sheep: a molecular look at diversity-driving genes

BACKGROUND

After domestication, the evolution of phenotypically-varied sheep breeds has generated rich biodiversity. This wide phenotypic variation arises as a result of hidden genomic changes that range from a single nucleotide to several thousands of nucleotides. Thus, it is of interest and significance to reveal and understand the genomic changes underlying the phenotypic variation of sheep breeds in order to drive selection towards economically important traits.

REVIEW

Various traits contribute to the emergence of variation in sheep phenotypic characteristics, including coat color, horns, tail, wool, ears, udder, vertebrae, among others. The genes that determine most of these phenotypic traits have been investigated, which has generated knowledge regarding the genetic determinism of several agriculturally-relevant traits in sheep. In this review, we discuss the genomic knowledge that has emerged in the past few decades regarding the phenotypic traits in sheep, and our ultimate aim is to encourage its practical application in sheep breeding. In addition, in order to expand the current understanding of the sheep genome, we shed light on research gaps that require further investigation.

CONCLUSIONS

Although significant research efforts have been conducted in the past few decades, several aspects of the sheep genome remain unexplored. For the full utilization of the current knowledge of the sheep genome, a wide practical application is still required in order to boost sheep productive performance and contribute to the generation of improved sheep breeds. The accumulated knowledge on the sheep genome will help advance and strengthen sheep breeding programs to face future challenges in the sector, such as climate change, global human population growth, and the increasing demand for products of animal origin.

Differential regulation of mRNAs and lncRNAs related to lipid metabolism in Duolang and Small Tail Han sheep

The function of long non-coding RNA (lncRNA) can be achieved through the regulation of target genes, and the deposition of fat is regulated by lncRNA. Fat has an important effect on meat quality. However, there are relatively few studies on lncRNAs in the subcutaneous adipose tissue of Duolang sheep and Small Tail Han sheep. In this study, RNA-Seq technology and bioinformatics methods were used to identify and analyze the lncRNA and mRNA in the subcutaneous adipose tissue of the two breeds of sheep. The results showed that 107 lnRNAs and 1329 mRNAs were differentially expressed. The differentially expressed genes and lncRNA target genes were significantly enriched in the biosynthesis of unsaturated fatty acids signaling pathway, fatty acid metabolism, adipocyte differentiation and other processes related to fat deposition. Among them, LOC105616076, LOC114118103, LOC105607837, LOC101116622, and LOC105603235 target FADS1, SCD, ELOVL6, HSD17B12 and HACD2, respectively. They play a key regulatory role in the biosynthesis of unsaturated fatty acids. This study lays a foundation for the study of the molecular mechanism of lncRNA on fat development, and has reference value for studying the differences in fat deposition between Duolang sheep and Small Tail Han sheep.

Transcriptome Profiling of Developing Ovine Fat Tail Tissue Reveals an Important Role for MTFP1 in Regulation of Adipogenesis

Fat-tail sheep exhibit a unique trait whereby substantial adipose tissue accumulates in the tail, a phenotype that is advantageous in many agroecological environments. In this study, we conducted histological assays, transcriptome analysis and functional assays to examine morphogenesis, characterize gene expression, and elucidate mechanisms that regulate fat tail development. We obtained the microstructure of tail before and after fat deposition, and demonstrated that measurable fat deposition occurred by the 80-day embryo (E80) stage, earlier than other tissues. Transcriptome profiling revealed 1,058 differentially expressed genes (DEGs) with six markedly different expression trends. GSEA enrichment and other downstream analyses showed important roles for genes and pathways involving in metabolism and that mitochondrial components were specifically overexpressed in the fat tail tissue of the 70-day embryo (E70). One hundred and eighty-three genes were further identified by leading edge gene analysis, among which, 17 genes have been reported in previous studies, including EEF1D, MTFP1, PPP1CA, PDGFD. Notably, the MTFP1 gene was highly correlated with the expression of other genes and with the highest enrichment score and gene expression change. Knockdown of MTFP1 in isolated adipose derived stem cells (ADSCs) inhibited cell proliferation and migration ability, besides, promoted the process of adipogenesis in vitro.

Ovine ELOVL5 and FASN genes polymorphisms and their correlations with sheep tail fat deposition

Reducing tail fat deposition can increase the economic value of a carcass and improve feed efficiency. This study aimed to explore ELOVL5 and FASN polymorphisms associated with tail fat deposition and their expression levels of sheep. Association analysis showed that ELOVL5 synonymous mutation g.62534 C > T was associated with tail width, tail fat weight, and relative tail fat weight (P < 0.05). FASN synonymous mutation g.12694 A > G was associated with tail length and width (P < 0.05). Combined effect analyses indicated significant differences between the combined genotypes and tail fat deposition. Quantitative real-time reverse transcription PCR indicated that the ELOVL5 and FASN expression levels were significantly higher in tail fat than in other tissues (P < 0.05). ELOVL5 expression levels in tail-fat tissue of big-tail sheep was significantly higher than that in small-tail sheep (P < 0.01). FASN expression levels were significantly higher in tail-fat tissue of small-tail sheep than in that of big-tail sheep (P < 0.05). During development, ELOVL5 tail fat expression increased significantly from 0 to 6 months old (P < 0.05), and FASN expression at 3 months old was significantly higher than that at 0 (minimum) and 6 months old (P < 0.05). Therefore, ELOVL5 and FASN polymorphisms could represent new candidate molecular markers and targets to reduce tail fat deposition in sheep.

Interference With ACSL1 Gene in Bovine Adipocytes: Transcriptome Profiling of mRNA and lncRNA Related to Unsaturated Fatty Acid Synthesis

The enzyme long-chain acyl-CoA synthetase 1 (ACSL1) is essential for lipid metabolism. The ACSL1 gene controls unsaturated fatty acid (UFA) synthesis as well as the formation of lipid droplets in bovine adipocytes. Here, we used RNA-Seq to determine lncRNA and mRNA that regulate UFA synthesis in bovine adipocytes using RNA interference and non-interference with ACSL1. The corresponding target genes of differentially expressed (DE) lncRNAs and the DE mRNAs were found to be enriched in lipid and FA metabolism-related pathways, according to GO and KEGG analyses. The differentially expressed lncRNA- differentially expressed mRNA (DEL-DEM) interaction network indicated that some DELs, such as TCONS_00069661, TCONS_00040771, TCONS_ 00035606, TCONS_00048301, TCONS_001309018, and TCONS_00122946, were critical for UFA synthesis. These findings assist our understanding of the regulation of UFA synthesis by lncRNAs and mRNAs in bovine adipocytes.

Transcriptome profiling of LncRNAs in sheep tail fat deposition

LncRNAs have recently received special attention due to their critical role in many important biological processes. There are few reports on its regulatory function in sheep fat deposition. In this study, two sheep populations with different tail types in Xinjiang, Bashibai sheep (fat-tailed) and the hybrid population of Bashibai sheep and wild argali (small-tailed) were selected for whole transcriptome sequencing from their tail tissues. First, 728 differentially expressed LncRNAs of tail fat between Bashibai and F2 sheep were identified by RNA-seq. Second, the tissue expression profile and relative expression difference between Bashibai and F2 sheep of 2 of 728 DE LncRNAs were analyzed by RT-PCR. LncRNA-MSTRG.24995 was highly expressed in tail fat, while lncRNA-MSTRG.36913 was highly expressed in subcutaneous fat. In addition, the expressions of LncRNA-MSTRG.24995 and LncRNA-MSTRG.36913 in tail fat of F2 sheep were significantly lower than that of Bashibai sheep, while those patterns in longissimus dorsi, quadriceps femoris and rumen were reversed. Third, the expression pattern of target genes FASN and THRSP in each tissue was similar with that of corresponding LncRNAs. The LncRNA-MSTRG.24995 directly affects tail fat deposition by FASN gene, while the LncRNA-MSTRG.36913 indirectly affects that by THRSP gene. This will help us to understand molecular mechanism of fat tail deposition from transcriptomic perspectives.

Integrative analysis of Iso-Seq and RNA-seq data reveals transcriptome complexity and differentially expressed transcripts in sheep tail fat

Background

Nowadays, both customers and producers prefer thin-tailed fat sheep. To effectively breed for this phenotype, it is important to identify candidate genes and uncover the genetic mechanism related to tail fat deposition in sheep. Accumulating evidence suggesting that post-transcriptional modification events of precursor-messenger RNA (pre-mRNA), including alternative splicing (AS) and alternative polyadenylation (APA), may regulate tail fat deposition in sheep. Differentially expressed transcripts (DETs) analysis is a way to identify candidate genes related to tail fat deposition. However, due to the technological limitation, post-transcriptional modification events in the tail fat of sheep and DETs between thin-tailed and fat-tailed sheep remains unclear.

Methods

In the present study, we applied pooled PacBio isoform sequencing (Iso-Seq) to generate transcriptomic data of tail fat tissue from six sheep (three thin-tailed sheep and three fat-tailed sheep). By comparing with reference genome, potential gene loci and novel transcripts were identified. Post-transcriptional modification events, including AS and APA, and lncRNA in sheep tail fat were uncovered using pooled Iso-Seq data. Combining Iso-Seq data with six RNA-sequencing (RNA-Seq) data, DETs between thin- and fat-tailed sheep were identified. Protein protein interaction (PPI) network, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were implemented to investigate the potential functions of DETs.

Results

In the present study, we revealed the transcriptomic complexity of the tail fat of sheep, result in 9,001 potential novel gene loci, 17,834 AS events, 5,791 APA events, and 3,764 lncRNAs. Combining Iso-Seq data with RNA-Seq data, we identified hundreds of DETs between thin- and fat-tailed sheep. Among them, 21 differentially expressed lncRNAs, such as ENSOART00020036299, ENSOART00020033641, ENSOART00020024562, ENSOART00020003848 and 9.53.1 may regulate tail fat deposition. Many novel transcripts were identified as DETs, including 15.527.13 (DGAT2), 13.624.23 (ACSS2), 11.689.28 (ACLY), 11.689.18 (ACLY), 11.689.14 (ACLY), 11.660.12 (ACLY), 22.289.6 (SCD), 22.289.3 (SCD) and 22.289.14 (SCD). Most of the identified DETs have been enriched in GO and KEGG pathways related to extracellular matrix (ECM). Our result revealed the transcriptome complexity and identified many candidate transcripts in tail fat, which could enhance the understanding of molecular mechanisms behind tail fat deposition.

Adipose gene expression profiles reveal insights into the adaptation of northern Eurasian semi-domestic reindeer (Rangifer tarandus)

Reindeer (Rangifer tarandus) are semi-domesticated animals adapted to the challenging conditions of northern Eurasia. Adipose tissues play a crucial role in northern animals by altering gene expression in their tissues to regulate energy homoeostasis and thermogenic activity. Here, we perform transcriptome profiling by RNA sequencing of adipose tissues from three different anatomical depots: metacarpal (bone marrow), perirenal, and prescapular fat in Finnish and Even reindeer (in Sakha) during spring and winter. A total of 16,212 genes are expressed in our data. Gene expression profiles in metacarpal tissue are distinct from perirenal and prescapular adipose tissues. Notably, metacarpal adipose tissue appears to have a significant role in the regulation of the energy metabolism of reindeer in spring when their nutritional condition is poor after winter. During spring, genes associated with the immune system are upregulated in the perirenal and prescapular adipose tissue. Blood and tissue parameters reflecting general physiological and metabolic status show less seasonal variation in Even reindeer than in Finnish reindeer. This study identifies candidate genes potentially involved in immune response, fat deposition, and energy metabolism and provides new information on the mechanisms by which reindeer adapt to harsh arctic conditions.

Transcriptome analysis of messenger RNA and long noncoding RNA related to different developmental stages of tail adipose tissues of sunite sheep

The tail fat of sheep is the most typical deposited fat, and it can be widely used in human daily life, such as diet, cosmetics, and industrial raw materials. To understand the potential regulatory mechanism of different growth stages of tail fat in Sunite sheep, we performed high-throughput RNA sequencing to characterize the long noncoding RNA (lncRNA) and messenger RNA (mRNA) expression profiles of the sheep tail fat at the age of 6, 18, and 30 months. A total of 223 differentially expressed genes (DEGs) and 148 differentially expressed lncRNAs were found in the tail fat of 6-, 18-, and 30-month-old sheep. Based on functional analysis, we found that fat-related DEGs were mainly expressed at 6 months of age and gradually decreased at 18 and 30 months of age. The target gene prediction analysis shows that most of the lncRNAs target more than 20 mRNAs as their transregulators. Further, we obtained several fat-related differentially expressed target genes; these target genes interact with different differentially expressed lncRNAs at various ages and play an important role in the development of tail fat. Based on the DEGs and differentially expressed lncRNAs, we established three co-expression networks for each comparison group. Finally, we concluded that the development of the sheep tail fat is more active during the early stage of growth and gradually decreases with the increase in age. The mutual regulation of lncRNAs and mRNAs may play a key role in this complex biological process.

Transcriptome Analysis Reveals Differentially Expressed Genes and Long Non-coding RNAs Associated With Fecundity in Sheep Hypothalamus With Different FecB Genotypes

Small-tailed Han sheep, with different FecB genotypes, manifest distinct ovulation rates and fecundities, which are due to differences in reproductive hormones secreted by the hypothalamic-pituitary-ovarian axis. Nevertheless, the function of the hypothalamus against a FecB mutant background on increasing ovulation rate is rarely reported. Therefore, we determined the expression profiles of hypothalamus tissue collected from six wild-type (WW) and six FecB mutant homozygous (BB) ewes at the follicular and luteal phases by whole-transcriptome sequencing. We identified 53 differentially expressed mRNAs (DEGs) and 40 differentially expressed long non-coding RNAs (DELs) between the two estrus states. Functional annotation analysis revealed that one of the DEGs, PRL, was particularly enriched in the hypothalamic function, hormone-related, and reproductive pathways. The lncRNA-target gene interaction networks and KEGG analysis in combination suggest that the lncRNAs LINC-676 and WNT3-AS cis-acting on DRD2 and WNT9B in different phases may induce gonadotropin-releasing hormone (GnRH) secretion. Furthermore, there were differences of regulatory elements and WNT gene family members involved in the follicular-luteal transition in the reproductive process between wild-type (WNT7A) and FecB mutant sheep (WNT9B). We combined the DEG and DEL data sets screened from different estrus states and genotypes. The overlap of these two sets was identified to select the mRNAs and lncRNAs that have major effects on ovulation. Among the overlapping molecules, seven DEGs and four DELs were involved in the follicular-luteal transition regulated by FecB mutation. Functional annotation analysis showed that two DEGs (FKBP5 and KITLG) were enriched in melanogenesis, oxytocin, and GnRH secretion. LINC-219386 and IGF2-AS were highly expressed in the BB ewes compared with WW ewes, modulating their target genes (DMXL2 and IGF2) to produce more GnRH during follicular development, which explains why mutated ewes produced more mature follicles. These results from expression profiling of the hypothalamus with the FecB mutation at different estrus states provide new insights into how the hypothalamus regulates ovulation under the effect of the FecB mutation.

The role of long noncoding RNAs in livestock adipose tissue deposition - A review

With the development of sequencing technology, numerous, long noncoding RNAs (lncRNAs) have been discovered and annotated. Increasing evidence has shown that lncRNAs play an essential role in regulating many biological and pathological processes, especially in cancer. However, there have been few studies on the roles of lncRNAs in livestock production. In animal products, meat quality and lean percentage are vital economic traits closely related to adipose tissue deposition. However, adipose tissue accumulation is also a pivotal contributor to obesity, diabetes, atherosclerosis, and many other diseases, as demonstrated by human studies. In livestock production, the mechanism by which lncRNAs regulate adipose tissue deposition is still unclear. In addition, the phenomenon that different animal species have different adipose tissue accumulation abilities is not well understood. In this review, we summarize the characteristics of lncRNAs and their four functional archetypes and review the current knowledge about lncRNA functions in adipose tissue deposition in livestock species. This review could provide theoretical significance to explore the functional mechanisms of lncRNAs in adipose tissue accumulation in animals.

Comparative Transcriptome Analysis of Key Genes and Pathways Activated in Response to Fat Deposition in Two Sheep Breeds With Distinct Tail Phenotype

Fat tail in sheep presents a valuable energy reserve that has historically facilitated adaptation to harsh environments. However, in modern intensive and semi-intensive sheep industry systems, breeds with leaner tails are more desirable. In the present study, RNA sequencing (RNA-Seq) was applied to determine the transcriptome profiles of tail fat tissues in two Chinese sheep breeds, fat-rumped Altay sheep and thin-tailed Xinjiang fine wool (XFW) sheep, with extreme fat tail phenotype difference. Then the differentially expressed genes (DEGs) and their sequence variations were further analyzed. In total, 21,527 genes were detected, among which 3,965 displayed significant expression variations in tail fat tissues of the two sheep breeds (P < 0.05), including 707 upregulated and 3,258 downregulated genes. Gene Ontology (GO) analysis disclosed that 198 DEGs were related to fat metabolism. In Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, the majority of DEGs were significantly enriched in "adipocytokine signaling," "PPAR signaling," and "metabolic pathways" (P < 0.05); moreover, some genes were involved in multiple pathways. Among the 198 DEGs, 22 genes were markedly up- or downregulated in tail fat tissue of Altay sheep, indicating that these genes might be closely related to the fat tail trait of this breed. A total of 41,724 and 42,193 SNPs were detected in the transcriptomic data of tail fat tissues obtained from Altay and XFW sheep, respectively. The distribution of seven SNPs in the coding regions of the 22 candidate genes was further investigated in populations of three sheep breeds with distinct tail phenotypes. In particular, the g.18167532T/C (Oar_v3.1) mutation of the ATP-binding cassette transporter A1 (ABCA1) gene and g.57036072G/T (Oar_v3.1) mutation of the solute carrier family 27 member 2 (SLC27A2) gene showed significantly different distributions and were closely associated with tail phenotype (P < 0.05). The present study provides transcriptomic evidence explaining the differences in fat- and thin-tailed sheep breeds and reveals numerous DEGs and SNPs associated with tail phenotype. Our data provide a valuable theoretical basis for selection of lean-tailed sheep breeds.

Molecular mechanisms of fat deposition: IL-6 is a hub gene in fat lipolysis, comparing thin-tailed with fat-tailed sheep breeds

Tail fat content affects meat quality and varies significantly among different breeds of sheep. Ghezel (fat-tailed) and Zel (thin-tailed) are two important Iranian local sheep breeds with different patterns of fat storage. The current study presents the transcriptome characterization of tail fat using RNA sequencing in order to get a better comprehension of the molecular mechanism of lipid storage in the two mentioned sheep breeds. Seven (Zel  =  4 and Ghezel  =  3) 7-month-old male lambs were used for this experiment. The results of sequencing were analyzed with bioinformatics methods, including differentially expressed genes (DEGs) identification, functional enrichment analysis, structural classification of proteins, protein-protein interaction (PPI) and network and module analyses. Some of the DEGs, such as LIPG, SAA1, SOCS3, HIF-1 α , and especially IL-6, had a close association with lipid metabolism. Furthermore, functional enrichment analysis revealed pathways associated with fat deposition, including "fatty acid metabolism", "fatty acid biosynthesis" and "HIF-1 signaling pathway". The structural classification of proteins showed that major down-regulated DEGs in the Zel (thin-tailed) breed were classified under transporter class and that most of them belonged to the solute carrier transporter (SLC) families. In addition, DEGs under the transcription factor class with an important role in lipolysis were up-regulated in the Zel (thin-tailed) breed. Also, network analysis revealed that IL-6 and JUNB were hub genes for up-regulated PPI networks, and HMGCS1, VPS35 and VPS26A were hub genes for down-regulated PPI networks. Among the up-regulated DEGs, the IL-6 gene seems to play an important role in lipolysis of tail fat in thin-tailed sheep breeds via various pathways such as tumor necrosis factor (TNF) signaling and mitogen-activated protein kinase (MAPK) signaling pathways. Due to the probable role of the IL-6 gene in fat lipolysis and also due to the strong interaction of IL-6 with the other up-regulated DEGs, it seems that IL-6 accelerates the degradation of lipids in tail fat cells.

Quantitative proteomic analysis identified differentially expressed proteins with tail/rump fat deposition in Chinese thin- and fat-tailed lambs

Tail adipose as one of the important functional tissues can enhance hazardous environments tolerance for sheep. The objective of this study was to gain insight into the underlying development mechanisms of this trait. A quantitative analysis of protein abundance in ovine tail/rump adipose tissue was performed between Chinese local fat- (Kazakh, Hu and Lanzhou) and thin-tailed (Alpine Merino, Tibetan) sheep in the present study by using lable-free approach. Results showed that 3400 proteins were identified in the five breeds, and 804 were differentially expressed proteins, including 638 up regulated proteins and 83 down regulated proteins in the tail adipose tissues between fat- and thin-tailed sheep, and 8 clusters were distinguished for all the DEPs’ expression patterns. The differentially expressed proteins are mainly associated with metabolism pathways and peroxisome proliferator activated receptor signaling pathway. Furthermore, the proteomics results were validated by quantitative real-time PCR and Western Blot. Our research has also suggested that the up-regulated proteins ACSL1, HSD17β4, FABP4 in the tail adipose tissue might contribute to tail fat deposition by facilitating the proliferation of adipocytes and fat accumulation in tail/rump of sheep. Particularly, FABP4 highly expressed in the fat-tail will play an important role for tail fat deposition. Our study might provide a novel view to understanding fat accumulation in special parts of the body in sheep and other animals.

Identification and co-expression analysis of long noncoding RNAs and mRNAs involved in the deposition of intramuscular fat in Aohan fine-wool sheep

BACKGROUND

Intramuscular fat (IMF) content has become one of the most important indicators for measuring meat quality, and levels of IMF are affected by various genes. Long non-coding RNAs (lncRNAs) are widely expressed non-coding RNAs that play an important regulatory role in a variety of biological processes; however, research on the lncRNAs involved in sheep IMF deposition is still in its infancy. Aohan fine-wool sheep (AFWS), one of China's most important meat-hair, dual-purpose sheep breed, provides a great model for studying the role of lncRNAs in the regulation of IMF deposition. We identified lncRNAs by RNA sequencing in Longissimus thoracis et lumborum (LTL) samples of sheep at two ages: 2 months (Mth-2) and 12 months (Mth-12).

RESULTS

We identified a total of 26,247 genes and 6935 novel lncRNAs in LTL samples of sheep. Among these, 199 mRNAs and 61 lncRNAs were differentially expressed. We then compared the structural characteristics of lncRNAs and mRNAs. We obtained target genes of differentially expressed lncRNAs (DELs) and performed enrichment analyses using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG). We found that target mRNAs were enriched in metabolic processes and developmental pathways. One pathway was significantly enriched, namely tight junction. Based on the analysis of critical target genes, we obtained seven candidate lncRNAs that potentially regulated lipid deposition and constructed a lncRNA-mRNA co-expression network that included MSTRG.4051.3-FZD4, MSTRG.16157.3-ULK1, MSTRG.21053.3-PAQR3, MSTRG.19941.2-TPI1, MSTRG.12864.1-FHL1, MSTRG.2469.2-EXOC6 and MSTRG.21381.1-NCOA1. We speculated that these candidate lncRNAs might play a role by regulating the expression of target genes. We randomly selected five mRNAs and five lncRNAs to verify the accuracy of the sequencing data by qRT-PCR.

CONCLUSIONS

Our study identified the differentially expressed mRNAs and lncRNAs during intramuscular lipid deposition in Aohan fine-wool sheep. The work may widen the knowledge about the annotation of the sheep genome and provide a working basis for investigating intramuscular fat deposition in sheep.

Transcriptional Profiles of Long Non-coding RNA and mRNA in Sheep Mammary Gland During Lactation Period

Sheep milk and related products have been growing in popularity around the world in recent years. However, the sheep milk industry is limited by low milk yield, and the molecular regulators of ovine lactation remain largely unknown. To investigate the transcriptomic basis of sheep lactation, RNA-Sequencing was used to explore the expression profiles of lncRNA and mRNA of the mammary gland in Hu sheep at three key time points during the lactation stage: 5 days before the expected date of parturition perinatal period (PP), 6 days after parturition early lactation (EL), and 25 days after parturition peak lactation (PL). A total of 1111, 688, and 54 differentially expressed (DE) lncRNAs as well as 1360, 660, and 17 DE mRNAs were detected in the EL vs PP, PL vs PP, and PL vs EL comparisons, respectively. Several prominent mRNAs (e.g., CSN1S1, CSN1S2, PAEP, CSN2, CSN3, and COL3A1) and lncRNAs (e.g., LNC_018483, LNC_005678, LNC_012936, and LNC_004856) were identified. Functional enrichment analysis revealed that several DE mRNAs and target genes of DE lncRNAs were involved in lactation-related pathways, such as MAPK, PPAR, and ECM-receptor interaction. This study enhances our understanding of how transcriptomic profiles change during the lactation period and pave the way for future studies examining sheep lactation.

Analysis of Transcriptome and miRNAome in the Muscle of Bamei Pigs at Different Developmental Stages

Simple Summary

The pigs is the most popular agricultural animal in the world. Muscle growth—which has the highest economic value in pigs—can be regulated by multiple genes and involves complex regulatory mechanisms. It is necessary to understand the dynamics of muscle transcriptome during development to understand the muscle development mechanism. However, the genes and miRNAs that play regulatory roles underlying differences in the meat quality of pigs remain unclear. In the current study, qRT-PCR, miRNA-Seq, and RNA-Seq were applied to analyze and verify muscle tissues of pigs from three different developmental stages and screened genes, miRNAs and pathways related to pig muscle development. This study focused on analyzing the mechanisms of muscle development and uncover the development differences in muscle from embryo to adult.

Abstract

The growth of skeletal muscle involves complex developmental processes that play an important part in the determinization of pork quality. The investigation of skeletal muscle mRNA or miRNA profiles is especially important for finding molecular approaches to improve meat quality in pig breeding. Therefore, we studied the transcriptome (mRNA and miRNA) profiles of skeletal muscle with RNA-Seq in three developmental stages of pigs: 65-day embryonic (E65), postnatal 0 days (natal) and 10 months (adult). We found 10,035, 9050 and 4841 differentially expressed (DE) genes for natal vs. E65, adult vs. E65 and adult vs. natal, 55, 101 and 85 DE miRNA for natal vs. E65, adult vs. E65 and adult vs. natal, respectively. In addition, the target genes of DE miRNA that was in a negative correlation with the corresponding miRNA in the same comparison group were selected for enrichment analysis. Gene Ontology terms were mainly classified into developmental processes. Pathway analysis revealed enrichment in the Rap1 signaling pathway, citrate cycle and oxidative phosphorylation and carbon. Finally, RT-PCR was employed for validating the level of expression of 11 DE miRNA and 14 DEGs. The transcriptome profiles of skeletal muscle from the different developmental stages of the Bamei pigs were obtained. From these data, hundreds of DE miRNA and mRNA, and the miRNA–mRNA regulatory network can provide valuable insights into further understanding of key molecular mechanisms and improving the meat quality in pig breeding.

miR-27a is an important adipogenesis regulator associated with differential lipid accumulation between intramuscular and subcutaneous adipose tissues of sheep

Micro ribonucleic acids (miRNAs) are crucial regulators for various biological processes. Despite important function in the proliferation and differentiation of preadipocytes, miRNA studies are limited in regional differences in adipogenesis. Here, we show that miR-27a plays an important role in regulating differential lipid accumulation between intramuscular (IM) and subcutaneous (SC) adipose tissues in sheep. Invivo, we observed that miR-27a expression in IM adipose tissue is more abundant than in SC adipose tissue. However, the expression of Peroxisome Proliferator-Activated Receptor Gamma (PPARG) and retinoid X receptor alpha (RXR alpha) in IM adipose tissue was significantly lower than that in SC adipose tissue. In the ovine preadipocyte differentiation model, we found that the expression of miR-27a was significantly decreased in differentiated ovine adipocytes. Overexpression of miR-27a significantly downregulated the expression of PPARG and RXR alpha and suppressed the accumulation of triglyceride but promoted the proliferation of ovine preadipocytes. Whereas, inhibition of miR-27a suppressed preadipocyte proliferation but enhanced PPARG and RXR alpha expression and lipid droplet formation. In addition, dual-luciferase activity assays showed that RXR alpha was a direct target of miR-27a. Thus, miR-27a enhances ovine preadipocytes proliferation and inhibits ovine preadipocytes differentiation through regulating the expression of target RXR alpha. Collectively, our study demonstrates the functional importance of miR-27a in ovine adipogenesis and provides novel insights into exploring regional differences in adipogenesis.

Genome-Wide Analysis and Function Prediction of Long Noncoding RNAs in Sheep Pituitary Gland Associated with Sexual Maturation

Long noncoding RNA (lncRNA) plays a crucial role in the hypothalamic-pituitary-testis (HPT) axis associated with sheep reproduction. The pituitary plays a connecting role in the HPT axis. However, little is known of their expression pattern and potential roles in the pituitary gland. To explore the potential lncRNAs that regulate the male sheep pituitary development and sexual maturation, we constructed immature and mature sheep pituitary cDNA libraries (three-month-old, TM, and nine-month-old, NM, respectively, n = 3) for lncRNA and mRNA high-throughput sequencing. Firstly, the expression of lncRNA and mRNA were comparatively analyzed. 2417 known lncRNAs and 1256 new lncRNAs were identified. Then, 193 differentially expressed (DE) lncRNAs and 1407 DE mRNAs were found in the pituitary between the two groups. Moreover, mRNA-lncRNA interaction network was constructed according to the target gene prediction of lncRNA and functional enrichment analysis. Five candidate lncRNAs and their targeted genes HSD17B12, DCBLD2, PDPK1, GPX3 and DLL1 that enriched in growth and reproduction related pathways were further filtered. Lastly, the interaction of candidate lncRNA TCONS_00066406 and its targeted gene HSD17B12 were validated in in vitro of sheep pituitary cells. Our study provided a systematic presentation of lncRNAs and mRNAs in male sheep pituitary, which revealed the potential role of lncRNA in male reproduction.

Genome-Wide Identification and Characterization of Long Noncoding RNAs of Brown to White Adipose Tissue Transformation in Goats

Long noncoding RNAs (lncRNAs) play an important role in the thermogenesis and energy storage of brown adipose tissue (BAT). However, knowledge of the cellular transition from BAT to white adipose tissue (WAT) and the potential role of lncRNAs in goat adipose tissue remains largely unknown. In this study, we analyzed the transformation from BAT to WAT using histological and uncoupling protein 1 (UCP1) gene analyses. Brown adipose tissue mainly existed within the goat perirenal fat at 1 day and there was obviously a transition from BAT to WAT from 1 day to 1 year. The RNA libraries constructed from the perirenal adipose tissues of 1 day, 30 days, and 1 year goats were sequenced. A total number of 21,232 lncRNAs from perirenal fat were identified, including 5393 intronic-lncRNAs and 3546 antisense-lncRNAs. Furthermore, a total of 548 differentially expressed lncRNAs were detected across three stages (fold change ≥ 2.0, false discovery rate (FDR) < 0.05), and six lncRNAs were validated by qPCR. Furthermore, trans analysis found lncRNAs that were transcribed close to 890 protein-coding genes. Additionally, a coexpression network suggested that 4519 lncRNAs and 5212 mRNAs were potentially in trans-regulatory relationships (r > 0.95 or r < -0.95). In addition, Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses showed that the targeted genes were involved in the biosynthesis of unsaturated fatty acids, fatty acid elongation and metabolism, the citrate cycle, oxidative phosphorylation, the mitochondrial respiratory chain complex, and AMP-activated protein kinase (AMPK) signaling pathways. The present study provides a comprehensive catalog of lncRNAs involved in the transformation from BAT to WAT and provides insight into understanding the role of lncRNAs in goat brown adipogenesis.

Deep transcriptome analysis using RNA-Seq suggests novel insights into molecular aspects of fat-tail metabolism in sheep

Fat-tail content of sheep breeds is varied and the molecular mechanisms regulating fat-tail development have not been well characterized. Aiming at better identifying the important candidate genes and their functional pathways contributing to fat deposition in the tail, a comparative transcriptome analysis was performed between fat- (Lori-Bakhtiari) and thin-tailed (Zel) Iranian sheep breeds using RNA-seq. The experiment was conducted on six male lambs (three lambs per each breed) at seven months of age. Four different combinations of aligners and statistical methods including Hisat2 + edgeR, Hisat2 + DESeq2, STAR + edgeR and STAR + DESeq2 were used to identify the differentially expressed genes (DEGs). The DEGs were selected for functional enrichment analysis and protein-protein interaction (PPI) network construction. Module analysis was also conducted to mine the functional sub-networks from the PPI network. In total, 264 genes including 80 up- and 184 down-regulated genes were identified as DEGs. The RNA-Seq results were validated by Q-RT-PCR. Functional analysis of DEGs and the module analysis of PPI network demonstrated that in addition to pathways affecting lipid metabolism, a series of enriched functional terms related to "response to interleukin", "MAPK signaling pathways", "Wnt signaling pathway", "ECM-receptor interaction", "regulation of actin cytoskeleton", and "response to cAMP" might contribute to the deposition of fat in tails of sheep. Overall results using RNA-Seq analysis characterized important candidate genes involved in the fatty acid metabolism and regulation of fat deposition, suggesting novel insights into molecular aspects of fat-tail metabolism in sheep. Selected DEGs should be further investigated as potential markers associated with the fat-tail development in sheep breeds.

GROWTH AND DEVELOPMENT SYMPOSIUM: STEM AND PROGENITOR CELLS IN ANIMAL GROWTH: Long noncoding RNAs in adipogenesis and adipose development of meat animals12

Sequencing technology, especially next-generation RNA sequencing, has greatly facilitated the identification and annotation of long noncoding RNAs (lncRNAs). In mammals, a large number of lncRNAs have been identified, which regulate various biological processes. An increasing number of lncRNAs have been identified which could function as key regulators of adipogenesis (adipocyte formation), a key step of the development of adipose tissue. Because proper adipose tissue development is a key factor affecting animal growth efficiency, lean/fat ratio, and meat quality, summarizing the roles and recent advances of lncRNAs in adipogenesis is needed in order to develop strategies to effectively manage fat deposition. In this review, we updated lncRNAs contributed to the regulation of adipogenesis, focusing on their roles in fat development of farm animals.

Genetic variations in the sheep SIRT7 gene and their correlation with body size traits

As a nicotinamide adenine dinucleotide (NAD)-dependent histone deacetylase and ADP ribosyl transferase, the silent information regulator 7 (Sirtuin 7, SIRT7) plays a crucial role in regulating the differentiation of adipocytes and myoblasts, lipid metabolism, glucose metabolism, and cellular growth in mammals. It has been hypothesized that SIRT7 affects growth traits in animals; therefore, in this study, the potential insertion/deletion (indel) of genetic variations within the ovine SIRT7 gene and their correlation with sheep growth traits were explored. A total of 709 individuals from five Chinese and Mongolian sheep breeds were analyzed. Two novel indel loci of the sheep SIRT7 gene were detected and were named 5′ promoter region-insertion-7 bp (5′ promoter region-7 bp) and 3′ UTR-insertion-17 bp (3′ UTR-17 bp), respectively. In all of the sheep breeds, frequencies of the 5′ promoter region-7 bp mutation were low, whereas mutations of 3′ UTR-17 bp were high in Tong sheep and Lanzhou fat-tail sheep (LFTS). Furthermore, both indel polymorphisms had significant associations with different growth characteristics (P<0.05). Among these associations, the 3′ UTR-17 bp was highly correlated with rump width in small-tail Han sheep (STHS, rams; P<0.01), and individuals with the ID genotype had better chest depth values than those with the II genotype. In this paper, two novel indels within the sheep SIRT7 gene were identified, and genetic diversity and its effects on body size traits were explored. These findings will potentially provide useful DNA markers for the improvement of economic traits in sheep genetic breeding.