Transcriptome study digs out BMP2 involved in adipogenesis in sheep tails

Multi-omics integrated analysis reveals the molecular mechanism of tail fat deposition differences in sheep with different tail types

BACKGROUND

The accumulation of tail fat in sheep is a manifestation of adaptive evolution to the environment. Sheep with different tail types show significant differences in physiological functions and tail fat deposition. Although these differences reflect the developmental mechanism of tail fat under different gene regulation, the situation of sheep tail fat tissue at the single cell level has not been explored, and its molecular mechanism still needs to be further elucidated.

RESULTS

Here, we characterized the genomic features of sheep with different tail types, detected the transcriptomic differences in tail adipose tissue between fat-tailed and thin-tailed sheep, established a single-cell atlas of sheep tail adipose tissue, and screened potential molecular markers (SESN1, RPRD1A and RASGEF1B) that regulate differences in sheep tail fat deposition through multi-omics integrated analysis. We found that the differential mechanism of sheep tail fat deposition not only involves adipocyte differentiation and proliferation, but is also closely related to cell-specific communication networks (When adipocytes act as signal outputters, LAMININ and other signal pathways are strongly expressed in guangling large tailed sheep and hu sheep), including interactions with immune cells and tissue remodeling to drive the typing of tail fat. In addition, we revealed the differentiation trajectory of sheep tail adipocytes through pseudo-time analysis and constructed the cell communication network of sheep tail adipose tissue.

CONCLUSIONS

Our results provide insights into the molecular mechanisms of tail fat deposition in sheep with different tail types, and provide a deeper explanation for the development and functional regulation of adipocytes.

Genomic insights into the population history of fat-tailed sheep and identification of two mutations that contribute to fat tail adipogenesis

INTRODUCTION

Since their domestication, domestic sheep (Ovis aries) have been culturally and economically significant farming animals worldwide. Fat-tailed sheep serve as a unique genetic resource for understanding adipogenesis and adaptive evolution in livestock.

OBJECTIVES

Several genomic analyses have been conducted on various sheep breeds to elucidate the genome and regulation mechanism of the fat tail trait, prior genomic studies have failed to reconcile conflicting evidence about the genetic basis of tail morphology, particularly regarding the roles of PDGFD and BMP2.

METHODS

Here, we conducted whole-genome resequencing of 283 sheep, encompassing 66 domestic breeds and 5 wild ovine species, to investigate the domestication history and selection signatures of fat-tailed sheep. Additionally, we performed transcriptome sequencing on adipose tissue to identify differentially expressed genes and cellular assays to validate these results.

RESULTS

Demographic analysis revealed that domestic sheep descended from Asiatic mouflon and fat-tailed sheep began to diverge from thin-tailed sheep approximately 4.4-7.5 thousand years ago in East Asia. Chinese indigenous sheep were classified into Mongolian, Kazakh, Tibetan, and Yunnan populations. The Yunnan population may have experienced more recent genetic introgression from wild species, rather than an independent domestication event. Moreover, many potential regions associated with the fat-tailed phenotype (DDI1, PDGFD, and BMP2) were identified by selective sweep and genome-wide association analyses. Additionally, a fine-scale analysis of fat-tailed and thin-tailed sheep revealed two novel mutations: a G/A missense variant of PDGFD (Chr15: 3900312) and a C/T missense variant of BMP2 (Chr13: 48462350), both of which were significantly associated with tail adiposity. Functional validation demonstrated that mutant A-PDGFD significantly activated PFGFD expression and reduced fat deposition compared to wildtype. The C-BMP2 mutant activated BMP2 expression and promoted preadipocyte fat deposition.

CONCLUSION

Our study provides the first evidence that these genes jointly regulate fat tail development through complementary mechanisms: PDGFD promotes adipose expansion, whereas BMP2 modulates energy partitioning. These findings offer new insights into the evolutionary history of fat-tailed sheep and identify potential targets for precision breeding in small ruminants.

Genetic and Epigenetic Adaptation Mechanisms of Sheep Under Multi-Environmental Stress Environment

Sheep (Ovis aries), domesticated from wild Asian mouflon ~10,000 years ago, are an important livestock species adapted to various ecological environments. Recent advancements in high-throughput sequencing and global environmental databases have facilitated the exploration of genetic-environmental associations, uncovering the genetic and epigenetic mechanisms behind sheep's adaptation to multiple environments. Studies show that HIF-1α and EPAS1 enhance high-altitude adaptation via hypoxic stress regulation; UCP1 contributes to cold adaptation through non-shivering thermogenesis; SLC4A4 and GPX3 increase drought resistance by regulating renal water reabsorption; and SOCS2 likely plays a role in metabolic and stress response regulation. Additionally, sheep adapt to temperature, drought, and environmental stress through DNA methylation, transcriptional regulation (e.g., SOD1, GPX4), heat shock proteins (e.g., HSP70), and metabolic pathways (e.g., UCP1). These findings offer valuable insights for improving sheep breeding and genetic enhancement. This review summarizes the mechanisms of adaptation to high altitude, cold, heat, drought, and comprehensive climate stress.

Integrated multiomic profiling of tail adipose tissue highlights novel genes, lipids, and metabolites involved in tail fat deposition in sheep

BACKGROUND

Tail fat is important for fat-tailed or fat-rumped sheep to survive in harsh environments. However, the molecular mechanism underlying tail fat deposition in sheep remains unclear. In this study, we comprehensively characterized the transcriptome, untargeted lipidome, and targeted metabolome profiles of the tail adipose tissues from Large-tailed Han sheep (long fat-tailed sheep) and Hu sheep (short fat-tailed sheep).

RESULTS

We identified 183 differentially expressed genes (DEGs), 55 differential lipids (DLs) and 17 differential metabolites (DMs) in the adipose tissues of the tails from Large-tailed Han and Hu sheep. Among the 183 DEGs selected (Q values ≤ 0.05 and│Log2(FC)│≥ 0.5), 18 DEGs, such as UCP3, ELOVL7 and GDF10, were directly associated with lipid metabolism identified via Gene Ontology (GO) analysis. Some genes, such as PPP3R1A, ADRA1, and DSLC46A2, were reportedly associated with lipid metabolism. A fold change ≥ 1.2 or ≤ 0.83 and a P-value < 0.05 were set as the default threshold to select the DLs and DMs. Among the 55 DLs, 36 DLs were phosphatidylcholines and 9 DLs were phosphatidylethanolamines. The top six DLs with the greatest differences in content were LPE (20:4) (up), PC (42:10) (up), PC (42:8) (up), PC (16:1/16:1) (down), PC (29:0) (down), and PC (32:2) (down). DMs related to the tricarboxylic acid cycle, such as D-glucose, cis-aconitic acid and citric acid were abundant in the tail fat of Large-tailed Han sheep. The DEGs, DLs and DMs were enriched mainly in the ferroptosis, the extracellular matrix (ECM)-receptor interaction, cGMP-PKG, calcium signaling and pathways related to cardiomyopathy and the tricarboxylic acid cycle.

CONCLUSION

This study obtained profiles of the transcriptome, lipidome and metabolome of the tail fat tissues of sheep with long and short fat tails. The findings suggested that ELOVL7, UCP3 and ferroptosis, ECM-receptor interaction pathways contributed to the difference in fat deposition, and phosphatidylcholines biosynthesis and tricarboxylic acid cycle may affect lipid metabolism in sheep tails. The results enhance our understanding of the differences in fat deposition in sheep tail.

Comprehensive multi-tissue epigenome atlas in sheep: A resource for complex traits, domestication, and breeding

Comprehensive functional genome annotation is crucial to elucidate the molecular mechanisms of agronomic traits in livestock, yet systematic functional annotation of the sheep genome is lacking. Here, we generated 92 transcriptomic and epigenomic data sets from nine major tissues, along with whole-genome data from 2357 individuals across 29 breeds worldwide, and 4006 phenotypic data related to tail fat weight. We constructed the first multi-tissue epigenome atlas in terms of functional elements, chromatin states, and their functions and explored the utility of the functional elements in interpreting phenotypic variation during sheep domestication and improvement. Particularly, we identified a total of 753,723 nonredundant functional elements, with over 60% being novel. We found tissue-specific promoters and enhancers related to sensory abilities and immune response that were highly enriched in genomic regions influenced by domestication, while longissimus dorsi tissue-specific active enhancers and tail fat tissue-specific active promoters were highly enriched in genomic regions influenced by breeding and improvement. Notably, a variant, Chr13:51760995A>C, located in an enhancer region, was identified as a causal variant for tail fat deposition based on multi-layered data sets. Overall, this research provides foundational resources and a successful case for future investigations of complex traits in sheep through the integration of multi-omics data sets.

Genetic analysis of key agronomic traits of local sheep breeds in Xinjiang, China

The formation of sheep (Ovis aries) breeds is influenced by different ecological environments and populations with different living habits, resulting in the development of germplasm resources with stable genetic key agronomic traits. Thus, investigating the genetic mechanisms behind various agronomic traits can enhance the conservation and utilization of diverse sheep breeds. Here, we explored the sheep variome and selection signatures using the Ovine Infinium HD SNP BeadChip (600 K SNPs) from 23 sheep breeds, comprising a total of 1215 sheep. The genetic mechanisms of wool quality and tail morphology were analyzed by selective sweep and genome-wide association study. Based on the results of within-population selective sweep analysis, we performed gene network analysis and divided them into 6 gene communities. We identified genetic regions containing genes linked to sheep wool and tail, which have been and may continue to be important targets for breeding and selection. Furthermore, our results revealed the expression profiles of genes in these regions across different biological systems. Our study provides insights into categorizing sheep breeds into distinct gene communities, as well as references for constructing genetic network pathways related to key agronomic traits in sheep and other domestic animals.

Identification of key modules and hub genes involved in regulating the feather follicle development of Wannan chickens using WGCNA

Carcass appearance is important economic trait, which affects customers in making purchase decisions. Both density and diameter of feather follicles are two important indicators of carcass appearance. However, the regulatory network and key genes be involved in feather follicle development remain poorly understood. To identify key genes and modules that involved in feather follicle development in chickens, 16 transcriptome datasets of Wannan chickens skin tissue (3 birds at the E9, E11, and E14, respectively, and 7 birds at the 12W) were used for weighted gene co-expression network analysis (WGCNA) analysis, and 12 skin tissue samples (3 birds for each stage) were selected for DEGs analysis. A total of 5,025, 2,337, and 10,623 DEGs were identified in 3 comparison groups, including the E9 vs. E11, the E11 vs. E14, and the E14 vs. 12W. Additionally, 31 co-expression gene modules were identified by WGCNA and the dark-orange, cyan, and blue module were found to be significantly associated with feather follicle development (p < 0.01). In total, 92,898 and 8,448 hub genes were obtained in the dark-orange, cyan, and blue modules, respectively. We focused on the cyan and blue modules, as 6 and 336 hub genes of these modules were identified to overlap with the DEGs of the three comparison groups, respectively. The 6 overlapped genes such as LAMC2, COL6A3, and COL6A2 etc., were over-represented in 12 categories such as focal adhesion and ECM-receptor interaction signaling pathway. Among the 336 genes that overlapped between the blue module and different DEGs comparison groups several genes including WNT7A and WNT9B were enriched in Wnt and ECM-receptor interaction signaling pathway. These results suggested that the LAMC2, COL6A3, COL6A2, WNT7A, and WNT9B genes may play a crucial role in the regulation of feather follicle development in Wannan chickens. Our results provided a reference for the molecular regulatory network and key genes in the development of feather follicles and contribute to molecular breeding for carcass appearance traits in chickens.

Full-Length Transcriptome and Gene Expression Analysis of Different Ovis aries Adipose Tissues Reveals Transcript Variants Involved in Lipid Biosynthesis

Sheep have historically been bred globally as a vital food source. To explore the transcriptome of adipose tissue and investigate key genes regulating adipose metabolism in sheep, adipose tissue samples were obtained from F1 Dorper × Hu sheep. High-throughput sequencing libraries for second- and third-generation sequencing were constructed using extracted total RNA. Functional annotation of differentially expressed genes and isoforms facilitated the identification of key regulatory genes and isoforms associated with sheep fat metabolism. SMRT-seq generated 919,259 high-accuracy cDNA sequences after filtering. Full-length sequences were corrected using RNA-seq sequences, and 699,680 high-quality full-length non-chimeric (FLNC) reads were obtained. Upon evaluating the ratio of total lengths based on FLNC sequencing, it was determined that 36,909 out of 56,316 multiple-exon isoforms met the criteria for full-length status. This indicates the identification of 330,375 full-length FLNC transcripts among the 370,114 multiple-exon FLNC transcripts. By comparing the reference genomes, 60,276 loci and 111,302 isoforms were identified. In addition, 43,423 new genes and 44,563 new isoforms were identified. The results identified 185 (3198), 394 (3592), and 83 (3286) differentially expressed genes (transcripts) between tail and subcutaneous, tail and visceral, and subcutaneous and visceral adipose tissues, respectively. Functional annotation and pathway analysis revealed the following observations. (1) Among the differentially expressed genes (DEGs) of TF and SF tissues, the downregulation of ACADL, ACSL6, and NC_056060.1.2536 was observed in SF, while FFAR4 exhibited upregulation. (2) Among the DEGs of TF and VF tissues, expressions of ACADL, ACSL6, COL1A1, COL1A2, and SCD were downregulated in VF, with upregulation of FFAR4. (3) Among SF and VF expressions of COL1A1, COL1A2, and NC_056060.1.2536 were downregulated in VF. Specific differentially expressed genes (ACADL, ACSL6, COL1A1, COL1A2, FFAR4, NC_056060.1.2536, and SCD) and transcripts (NC_056066.1.1866.16 and NC_056066.1.1866.22) were identified as relevant to fat metabolism. These results provide a dataset for further verification of the regulatory pathway associated with fat metabolism in sheep.

Dynamic alterations in metabolomics and transcriptomics associated with intestinal fibrosis in a 2,4,6-trinitrobenzene sulfonic acid-induced murine model

BACKGROUND & AIMS

Intestinal fibrosis is a common and severe complication of inflammatory bowel disease without clear pathogenesis. Abnormal expression of host genes and metabolic perturbations might associate with the onset of intestinal fibrosis. In this study, we aimed to investigate the relationship between the development of intestinal fibrosis and the dynamic alterations in both fecal metabolites and host gene expression.

METHODS

We induced intestinal fibrosis in a murine model using 2,4,6-trinitrobenzene sulfonic acid (TNBS). TNBS-treated or control mice were sacrificed after 4 and 6 weeks of intervention; alterations in colonic genes and fecal metabolites were determined by transcriptomics and metabolomics, respectively. Differential, tendency, enrichment, and correlation analyses were performed to assess the relationship between host genes and fecal metabolites.

RESULTS

RNA-sequencing analysis revealed that 679 differential genes with enduring changes were mainly enriched in immune response-related signaling pathways and metabolism-related biological processes. Among them, 15 lipid metabolism-related genes were closely related to the development of intestinal fibrosis. Moreover, the fecal metabolic profile was significantly altered during intestinal fibrosis development, especially the lipid metabolites. Particularly, dynamic perturbations in lipids were strongly associated with alterations in lipid metabolism-related genes expression. Additionally, six dynamically altered metabolites might serve as biomarkers to identify colitis-related intestinal fibrosis in the murine model.

CONCLUSIONS

Intestinal fibrosis in colitis mice might be related to dynamic changes in gene expression and metabolites. These findings could provide new insights into the pathogenesis of intestinal fibrosis.

Transcriptome analysis of perirenal fat from Spanish Assaf suckling lamb carcasses showing different levels of kidney knob and channel fat

Introduction

Suckling lamb meat is highly appreciated in European Mediterranean countries because of its mild flavor and soft texture. In suckling lamb carcasses, perirenal and pelvic fat depots account for a large fraction of carcass fat accumulation, and their proportions are used as an indicator of carcass quality.

Material and Methods

This study aimed to characterize the genetic mechanisms that regulate fat deposition in suckling lambs by evaluating the transcriptomic differences between Spanish Assaf lambs with significantly different proportions of kidney knob and channel fat (KKCF) depots in their carcasses (4 High-KKCF lambs vs. 4 Low-KKCF lambs).

Results

The analyzed fat tissue showed overall dominant expression of white adipose tissue gene markers, although due to the young age of the animals (17-36 days), the expression of some brown adipose tissue gene markers (e.g., UCP1, CIDEA) was still identified. The transcriptomic comparison between the High-KKCF and Low-KKCF groups revealed a total of 80 differentially expressed genes (DEGs). The enrichment analysis of the 49 DEGs with increased expression levels in the Low-KKCF lambs identified significant terms linked to the biosynthesis of lipids and thermogenesis, which may be related to the higher expression of the UCP1 gene in this group. In contrast, the enrichment analysis of the 31 DEGs with increased expression in the High-KKCF lambs highlighted angiogenesis as a key biological process supported by the higher expression of some genes, such as VEGF-A and THBS1, which encode a major angiogenic factor and a large adhesive extracellular matrix glycoprotein, respectively.

Discussion

The increased expression of sestrins, which are negative regulators of the mTOR complex, suggests that the preadipocyte differentiation stage is being inhibited in the High-KKCF group in favor of adipose tissue expansion, in which vasculogenesis is an essential process. All of these results suggest that the fat depots of the High-KKCF animals are in a later stage of development than those of the Low-KKCF lambs. Further genomic studies based on larger sample sizes and complementary analyses, such as the identification of polymorphisms in the DEGs, should be designed to confirm these results and achieve a deeper understanding of the genetic mechanisms underlying fat deposition in suckling lambs.

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.

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.