Search

Abstract Studies of preadipocyte differentiation and fat deposition in sheep have mainly focused on functional genes, and with no emphasis placed on the role that long non-coding RNAs (lncRNAs) may have on the activity of those genes. Here, the expression profile of lncRNAs in ovine preadipocyte differentiation was investigated and the differentially expressed lncRNAs were screened on day 0 (D0), day 2(D2) and day 8(D8) of ovine preadipocyte differentiation, with their target genes being predicted. The competing endogenous RNA (ceRNA) regulatory network was constructed by GO and KEGG enrichment analysis for functional annotation, and some differentially expressed lncRNAs were randomly selected to verify the RNA-Seq results by RT-qPCR. In the study, a total of 2517 novel lncRNAs and 3943 known lncRNAs were identified from ovine preadipocytes at the three stages of differentiation, with the highest proportion being intergenic lncRNAs. A total of 3455 lncRNAs were expressed at all three stages of preadipocyte differentiation, while 214, 226 and 228 lncRNAs were uniquely expressed at day 0, day 2 and day 8, respectively. By comparing the expression of the lncRNAs between the three stages of differentiation stages, a total of 405, 272 and 359 differentially expressed lncRNAs were found in D0-vs-D2, D0-vs-D8, and D2-vs-D8, respectively. Functional analysis revealed that the differentially expressed lncRNAs were enriched in signaling pathways related to ovine preadipocyte differentiation, such as mitogen-activated protein kinase (MAPK) pathway, the phosphoinositide 3-kinase protein kinase B (PI3K-Akt) pathway, and the transforming growth factor beta (TGF-β) pathway. In summary, lncRNAs from preadipocytes at different stages of differentiation in sheep were identified and screened using RNA-Seq technology, and the regulatory mechanisms of lncRNAs in preadipocyte differentiation and lipid deposition were explored. This study provides a theoretical reference for revealing the roles of lncRNAs in ovine preadipocyte differentiation and also offers a theoretical basis for further understanding the regulatory mechanisms of ovine preadipocyte differentiation.

Yaks are tough animals living in Tibet’s hypoxic stress environment. However, the metabolite composition of yak milk and its role in hypoxic stress tolerance remains largely unexplored. The similarities and differences between yak and human milk in hypoxic stress tolerance are also unclear. This study explored yak colostrum (YC) and yak mature milk (YMM) using GC–MS, and 354 metabolites were identified in yak milk. A comparative metabolomic analysis of yak and human milk metabolites showed that over 70% of metabolites were species-specific. Yak milk relies mainly on essential amino acids- arginine and essential branched-chain amino acids (BCAAs): L-isoleucine, L-leucine, and L-valine tolerate hypoxic stress. To slow hypoxic stress, human breast milk relies primarily on the neuroprotective effects of non-essential amino acids or derivates, such as citrulline, sarcosine, and creatine. In addition, metabolites related to hypoxic stress were significantly enriched in YC than in YMM. These results reveal the unique metabolite composition of yak and human milk and provide practical information for applying yak and human milk to hypoxic stress tolerance.

High-altitude stress threatens the survival rate of Tibetan sheep and reduces their fertility. However, the molecular basis of this phenomenon remains elusive. Here, we used RNA-seq to elucidate the transcriptome dynamics of high-altitude stress in Tibetan sheep ovaries. In total, 104 genes were characterized as high-altitude stress-related differentially expressed genes (DEGs). In addition, 36 DEGs contributed to ovarian follicle development, and 28 of them were downregulated under high-altitude stress. In particular, high-altitude stress significantly suppressed the expression of two ovarian lymphatic system marker genes: LYVE1 and ADAMTS-1. Network analysis revealed that luteinizing hormone (LH)/follicle-stimulating hormone (FSH) signaling-related genes, such as EGR1, FKBP5, DUSP1, and FOS, were central regulators in the DEG network, and these genes were also suppressed under high-altitude stress. As a post-transcriptional regulation mechanism, alternative splicing (AS) is ubiquitous in Tibetan sheep. High-altitude stress induced 917 differentially alternative splicing (DAS) events. High-altitude stress modulated DAS in an AS-type-specific manner: suppressing skipped exon events but increasing retained intron events. C2H2-type zinc finger transcription factors and RNA processing factors were mainly enriched in DAS. These findings revealed high-altitude stress repressed ovarian development by suppressing the gene expression of LH/FSH hormone signaling genes and inducing intron retention of C2H2-type zinc finger transcription factors.

MicroRNAs (miRNAs) are crucial regulatory molecules in lipid deposition and metabolism. However, the effect of miR-200b on the regulation of proliferation and adipogenesis of ovine preadipocytes is unknown in the sheep (Ovis aries). In this study, the expression profiles of miR-200b were investigated in the seven tissues of Tibetan ewes and differentiated preadipocytes. The effect of miR-200b, as well as its target genes p27 and KLF9, on the proliferation of ovine preadipocytes and adipogenesis was also investigated, using cell viability analysis, EdU staining, Oil Red O staining and reverse transcription-quantitative PCR (RT-qRCR). The miR-200b was expressed in all the tissues investigated, and it was highly expressed in lung, liver, subcutaneous adipose and spleen tissues. The expression of miR-200b continuously decreased when the differentiation of ovine preadipocytes initiated. The miR-200b mimic dramatically accelerated the proliferation but inhibited differentiation of ovine preadipocytes. The miR-200b inhibitor resulted in an opposite effect on the proliferation and differentiation of ovine preadipocytes. The dual luciferase reporter assay results showed that miR-200b mimic significantly decreased the luciferase activity of p27 and KLF9 in HEK293 cells transfected with wild-type dual luciferase reporter vectors. This suggests that p27 and KLF9 are the target genes of miR-200b. In over-expressed-p27 preadipocytes, the number of EdU-labeled preadipocytes and the expression levels of proliferation marker genes CDK2, CDK4, CCND1 and PCNA significantly decreased. In addition, the transfection of over-expressed-KLF9 vector into adipocytes remarkably increased the accumulation of lipid droplets and the expression levels of differentiation marker genes aP2, PPARγ, LPL and GLUT4. These results suggest that miR-200b accelerated the proliferation but inhibited the adipogenic differentiation of ovine preadipocytes by targeting p27 and KLF9, respectively.

MicroRNAs (miRNAs) have been found to be involved in lipid deposition and metabolism. However, there have been no reports on the roles of miR-148a in the proliferation and adipogenesis of preadipocytes in sheep. In this study, the expression of miR-148a was profiled in the eight tissues of Tibetan ewes and differentiated preadipocytes, and the role of miR-148a in differentiation and proliferation of ovine preadipocytes was investigated using Oil Red O staining, CCK-8, EdU staining, cell cycle detection, and RT-qPCR. The effect of PTEN on the differentiation of ovine preadipocytes was also investigated. The miR-148a was widely expressed in the eight tissues investigated and had significantly increased expression in liver, spleen and subcutaneous adipose tissues, and the heart. The expression of miR-148a continued to increase with the differentiation of ovine preadipocytes. The over-expression of miR-148a significantly promoted differentiation but inhibited the proliferation of ovine preadipocytes. The inhibition of miR-148a had the opposite effect on the differentiation and proliferation of ovine preadipocytes with over-expressed miR-148a. The results from the dual luciferase reporter assays showed that miR-148a mimic significantly decreased the luciferase activity of PTEN-3′UTR dual luciferase reporter vector, suggesting that PTEN is a target gene of miR-148a. In over-expressed-PTEN preadipocytes, the number of lipid droplets remarkably decreased, and the expression levels of adipogenesis marker genes PPARγ, FASN, FATP4, GLUT4, C/EBPβ and LPL were also significantly down-regulated. These results suggest that miR-148a accelerated the adipogenic differentiation of ovine preadipocytes by inhibiting PTEN expression, and also inhibited the proliferation of ovine preadipocytes.

Long non-coding RNAs (lncRNAs) are a kind of non-coding RNA with >200 nucleotides in length. Some lncRNAs have been proven to have clear regulatory functions in many biological processes of mammals. However, there have been no reports on the roles of lncRNAs in ovine mammary gland tissues. In the study, the expression profiles of lncRNAs were studied using RNA-Seq in mammary gland tissues from lactating Small-Tailed Han (STH) ewes and Gansu Alpine Merino (GAM) ewes with different milk yield and ingredients. A total of 1894 lncRNAs were found to be expressed. Compared with the GAM ewes, the expression levels of 31 lncRNAs were significantly up-regulated in the mammary gland tissues of STH ewes, while 37 lncRNAs were remarkably down-regulated. Gene Ontogeny (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis found that the target genes of differentially expressed lncRNAs were enriched in the development and proliferation of mammary epithelial cells, morphogenesis of mammary gland, ErbB signaling pathway, and Wnt signaling pathway. Some miRNA sponges of differentially expressed lncRNAs, reported to be associated with lactation and mammary gland morphogenesis, were found in a lncRNA-miRNA network. This study reveals comprehensive lncRNAs expression profiles in ovine mammary gland tissues, thereby providing a further understanding of the functions of lncRNAs in the lactation and mammary gland development of sheep.

Variation in some caprine keratin-associated protein (KAP) genes has been associated with cashmere fiber traits, but many KAP genes remain unidentified in goats. In this study, we confirm the identification of a KAP27-1 gene (KRTAP27-1) and describe its effect on cashmere traits in 248 Longdong cashmere goats. A polymerase chain reaction–single strand conformation polymorphism (PCR-SSCP) analysis was used to screen for sequence variation in this gene, and three sequence variants (named A to C) were found. These sequences have the highest similarity (77% identity) to a human KRTAP27-1 sequence, while sharing some homology with a predicted caprine KRTAP27-1 sequence ENSCHIG00000023347 in the goat genome construct (ARS1:CM004562.1) at chromosome 1 position 3,966,193–3,973,677 in the forward strand. There were two single nucleotide polymorphisms (SNPs) detected in the coding sequence, including one nonsynonymous SNP (c.413C/T; p.Ala138Val) and one synonymous SNP (c.495C/T). The C variant differed from A and B at c.413C/T, having cytosine in its nucleotide sequence, while the B variant differed from A and C at c.495C/T, having thymine in its nucleotide sequence. Goats of the genotypes AB and BB produced cashmere fibers of higher mean fiber diameter (MFD) than goats of genotype AA, but no difference in MFD was detected between the AB and BB goats. These results suggest that B is associated with increased MFD. Expression of the caprine KRTAP27-1 sequence was predominantly detected in the skin tissue of goats but not or only weakly detected in other tissues, including longissimus dorsi muscle, heart, kidney, liver, lung and spleen.