microRNA Temporal-Specific Expression Profiles Reveal longissimus dorsi Muscle Development in Tianzhu White Yak.

As a class of regulatory factors, microRNAs (miRNAs) play an important role in regulating normal muscle development and fat deposition. Muscle and adipose tissues, as major components of the animal organism, are also economically important traits in livestock production. However, the effect of miRNA expression profiles on the development of muscle and adipose tissues in yak is currently unknown. In this study, we performed RNA sequencing (RNA-Seq) on Tianzhu white yak longissimus dorsi muscle tissue obtained from calves (6 months of age, M6, n = 6) and young (30 months of age, M30, n = 6) and adult yak (54 months of age, M54, n = 6) to identify which miRNAs are differentially expressed and to investigate their temporal expression profiles, establishing a regulatory network of miRNAs associated with the development of muscle and adipose. The results showed that 1191 miRNAs and 22061 mRNAs were screened across the three stages, of which the numbers of differentially expressed miRNAs (DE miRNAs) and differentially expressed mRNAs (DE mRNAs) were 225 and 450, respectively. The expression levels of the nine DE miRNAs were confirmed using a reverse transcription quantitative PCR (RT-qPCR) assay, and the trend of the assay results was generally consistent with the trend of the transcriptome profiles. Based on the expression trend, DE miRNAs were categorized into eight different expression patterns. Regarding the expression of DE miRNAs in sub-trends Profile 1 and Profile 2 (p < 0.05), the gene expression patterns were upregulated (87 DE miRNAs). Gene ontology (GO) and Kyoto Encyclopedia of Genes Genomes (KEGG) analyses showed that the identified DE miRNAs and DE mRNAs were enriched in pathway entries associated with muscle and intramuscular fat (IMF) growth and development. On this basis, we constructed a DE miRNA-mRNA interaction network. We found that some DE mRNAs of interest overlapped with miRNA target genes, such as ACSL3, FOXO3, FBXO30, FGFBP4, TSKU, MYH10 (muscle development), ACOX1, FADS2, EIF4E2, SCD1, EL0VL5, and ACACB (intramuscular fat deposition). These results provide a valuable resource for further studies on the molecular mechanisms of muscle tissue development in yak and also lay a foundation for investigating the interactions between genes and miRNAs.

MicroRNA-148a Targets DNMT1 and PPARGC1A to Regulate the Viability, Proliferation, and Milk Fat Synthesis of Ovine Mammary Epithelial Cells.

In this study, the expression profiles of miR-148a were constructed in eight different ovine tissues, including mammary gland tissue, during six different developmental periods. The effect of miR-148a on the viability, proliferation, and milk fat synthesis of ovine mammary epithelial cells (OMECs) was investigated, and the target relationship of miR-148a with two predicted target genes was verified. The expression of miR-148a exhibited obvious tissue-specific and temporal-specific patterns. miR-148a was expressed in all eight ovine tissues investigated, with the highest expression level in mammary gland tissue (p < 0.05). Additionally, miR-148a was expressed in ovine mammary gland tissue during each of the six developmental periods studied, with its highest level at peak lactation (p < 0.05). The overexpression of miR-148a increased the viability of OMECs, the number and percentage of Edu-labeled positive OMECs, and the expression levels of two cell-proliferation marker genes. miR-148a also increased the percentage of OMECs in the S phase. In contrast, transfection with an miR-148a inhibitor produced the opposite effect compared to the miR-148a mimic. These results indicate that miR-148a promotes the viability and proliferation of OMECs in Small-tailed Han sheep. The miR-148a mimic increased the triglyceride content by 37.78% (p < 0.01) and the expression levels of three milk fat synthesis marker genes in OMECs. However, the miR-148a inhibitor reduced the triglyceride level by 87.11% (p < 0.01). These results suggest that miR-148a promotes milk fat synthesis in OMECs. The dual-luciferase reporter assay showed that miR-148a reduced the luciferase activities of DNA methyltransferase 1 (DNMT1) and peroxisome proliferator-activated receptor gamma coactivator 1-A (PPARGC1A) in wild-type vectors, suggesting that they are target genes of miR-148a. The expression of miR-148a was highly negatively correlated with PPARGC1A (r = -0.789, p < 0.001) in ovine mammary gland tissue, while it had a moderate negative correlation with DNMT1 (r = -0.515, p = 0.029). This is the first study to reveal the molecular mechanisms of miR-148a underlying the viability, proliferation, and milk fat synthesis of OMECs in sheep.

MicroRNA expression profiles reveal wool development and fineness regulation in Gansu alpine fine-wool sheep.

The development of wool has a complex regulatory mechanism both influenced by genetic and environmental factors. MicroRNAs (miRNA) were involved in various biological processes of animals, and may play an important role in the regulation of wool development. In this study, we comprehensively analyzed and identified the histological parameters of hair follicles, as well as the miRNAs, target genes, pathways, and Gene Ontology terms related to wool fineness regulation and wool growth and development using HE staining and RNA-Seqs methods. Both coarse (group C, mean fiber diameter (MFD) = 22.26 ± 0.69 µm, n = 6) and fine (group F, MFD = 16.91 ± 0.29 µm, n = 6) of Gansu alpine fine-wool sheep with different wool fineness were used in this study. The results showed that the primary follicle diameter and secondary wool fiber diameter in group C were significantly higher than those in group F (P < 0.05). And the number of primary and secondary hair follicles in group C was significantly lower than that in group F (P < 0.05). Furthermore, a total of 67 DE miRNAs and 290 potential DE miRNAs target genes were screened in the skin tissues of sheep from groups F and C, and some potential target genes related to wool fineness regulation were screened, such as CDH2, KRT82, FOXN1, LOC101106296, KRT20, MCOLN3, KRT71, and TERT. These genes were closely related to Glutathione metabolism, epidermal cell differentiation, keratinization, and regulation of hair cycle. Moreover, the regulatory network of miRNAs-mRNAs suggested that miRNAs (miR-129-x, novel m0079-3p, miR-2484-z, novel m0025-5P, etc.) may play a key role in the wool development and wool fineness regulation of Gansu alpine fine-wool sheep. In summary, this study expands the existing miRNAs database and provides new information for studying the regulation of wool development in Gansu alpine fine wool sheep.

Circular RNA_015343 sponges microRNA-25 to regulate viability, proliferation, and milk fat synthesis of ovine mammary epithelial cells via INSIG1.

In our previous study, circ_015343 was found to inhibit the viability and proliferation of ovine mammary epithelial cells (OMECs) and the expression levels of milk fat synthesis marker genes, but the regulatory mechanism underlying the processes is still unclear. Accordingly in this study, the target relationships between circ_015343 with miR-25 and between miR-25 with insulin induced gene 1 (INSIG1) were verified, and the functions of miR-25 and INSIG1 were investigated in OMECs. The dual-luciferase reporter assay revealed that miR-25 mimic remarkably decreased the luciferase activity of circ_015343 in HEK293T cells cotransfected with a wild-type vector, while it did not change the activity of circ_015343 in HEK293T cells cotransfected with a mutant vector. These suggest that cic_015343 can adsorb and bind miR-25. The miR-25 increased the viability and proliferation of OMECs, and the content of triglycerides in OMECs. In addition, INSIG1 was found to be a target gene of miR-25 using a dual-luciferase reporter assay. Overexpression of INSIG1 decreased the viability, proliferation, and level of triglycerides of OMECs. In contrast, the inhibition of INSIG1 in expression had the opposite effect on activities and triglycerides of OMECs with overexpressed INSIG1. A rescue experiment revealed that circ_015343 alleviated the inhibitory effect of miR-25 on the mRNA and protein abundance of INSIG1. These results indicate that circ_015343 sponges miR-25 to inhibit the activities and content of triglycerides of OMECs by upregulating the expression of INSIG1 in OMECs. This study provided new insights for understanding the genetic molecular mechanism of lactation traits in sheep.

Molecular Characteristics of JAK2 and Its Effect on the Milk Fat and Casein Synthesis of Ovine Mammary Epithelial Cells.

In addition to its association with milk protein synthesis via the Janus kinase-signal transducer and activator of transcription (JAK-STAT) pathway, JAK2 also affects milk fat synthesis. However, to date, there have been no reports on the effect of JAK2 on ovine mammary epithelial cells (OMECs), which directly determine milk yield and milk contents. In this study, the coding sequence (CDS) region of ovine JAK2 was cloned and identified and its tissue expression and localization in ovine mammary glands, as well as its effects on the viability, proliferation, and milk fat and casein levels of OMECs, were also investigated. The CDS region of ovine JAK2, 3399 bp in length, was cloned and its authenticity was validated by analyzing its sequence similarity with JAK2 sequences from other animal species using a phylogenetic tree. JAK2 was found to be expressed in six ovine tissues, with the highest expression being in the mammary gland. Over-expressed JAK2 and three groups of JAK2 interference sequences were successfully transfected into OMECs identified by immunofluorescence staining. When compared with the negative control (NC) group, the viability of OMECs was increased by 90.1% in the pcDNA3.1-JAK2 group. The over-expression of JAK2 also increased the number and ratio of EdU-labeled positive OMECs, as well as the expression levels of three cell proliferation marker genes. These findings show that JAK2 promotes the viability and proliferation of OMECs. Meanwhile, the triglyceride content in the over-expressed JAK2 group was 2.9-fold higher than the controls and the expression levels of four milk fat synthesis marker genes were also increased. These results indicate that JAK2 promotes milk fat synthesis. Over-expressed JAK2 significantly up-regulated the expression levels of casein alpha s2 (CSN1S2), casein beta (CSN2), and casein kappa (CSN3) but down-regulated casein alpha s1 (CSN1S1) expression. In contrast, small interfered JAK2 had the opposite effect to JAK2 over-expression on the viability, proliferation, and milk fat and milk protein synthesis of OMECs. In summary, these results demonstrate that JAK2 promotes the viability, proliferation, and milk fat synthesis of OMECs in addition to regulating casein expression in these cells. This study contributes to a better comprehension of the role of JAK2 in the lactation performance of sheep.

Expression Localization of the KRT32 Gene and Its Association of Genetic Variation with Wool Traits.

Changes in keratin gene expression and spatiotemporal regulation determine the compositional content and cellular localization of wool keratin, thereby affecting wool traits. Therefore, keratin gene family member 32 (KRT32) was selected for a study using RT-qPCR, immunofluorescence, and penta-primer amplification refractory mutation system (PARMS) techniques. The results showed that KRT32 mRNA was highly expressed in the skin and localized to the inner root sheath (IRS), outer root sheath (ORS) and dermal papilla (DP). Sequencing results identified eight SNPs in KRT32, and association analyses revealed that the variations were significantly associated with multiple traits in wool (p < 0.05), including MFD, CF and MFC. The constructed haplotype combination H2H3 has higher CF and smaller MFD than other haplotype combination (p < 0.05). In conclusion, KRT32 can be used as a candidate gene for molecular genetic improvement of wool in Gansu Alpine Fine-wool sheep.

Screening and identification of lncRNAs in preadipocyte differentiation in sheep.

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.

Spatiotemporal Expression and Haplotypes Identification of KRT84 Gene and Their Association with Wool Traits in Gansu Alpine Fine-Wool Sheep.

Keratin (K) is a major protein component of hair and is involved in hair growth and development. In this study, we analysed the expression, localization, and polymorphism of the K84 gene (KRT84) in Gansu Alpine Fine-wool sheep using immunofluorescence, RT-qPCR, and PARMS (penta-primer amplification refractory mutation system). Haplotypes of KRT84 were also constructed and their relationship with wool traits analysed. It was revealed that KRT84 was highly expressed in hair follicles, including the inner root sheath, outer root sheath, and hair medulla and at all six lamb ages investigated from 1 to 270 days of age. Three SNPs were detected in KRT84 exon 1, and they formed three haplotypes (named H1, H2, and H3) and six genotypes. Analyses revealed an association between haplotype combinations (diplotypes) and the mean fibre curvature, mean staple length, mean staple strength, mean fibre diameter, the coefficient of variation of fibre diameter, and comfort factor for these sheep. These results suggest that KRT84 is of importance in determining several key traits in Gansu Alpine Fine-wool sheep and that the gene could possibly be used as a genetic marker for wool trait selection in these sheep.

Identification and screening of circular RNAs during adipogenic differentiation of ovine preadipocyte by RNA-seq.

Circular RNAs (circRNAs) are a class of non-coding RNAs that play important roles in preadipocyte differentiation and adipogenesis. However, little is known about genome-wide identification, expression profile, and function of circRNAs in sheep. To investigate the role of circRNAs during ovine adipogenic differentiation, the subcutaneous adipose tissue of Tibetan rams was collected in June 2022. Subsequently, the preadipocytes were immediately isolated from collected adipose tissue and then induced to begin differentiation. The adipocytes samples cultured on days 0, 2, and 8 of preadipocytes differentiation were used to perform RNA sequencing (RNA-seq) analysis to construct the expression profiles of circRNAs. Subsequently, the function of differentially expressed circRNAs was investigated by performing the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of their parent genes. Finally, a circRNAs-miRNAs-mRNAs network involved in adipogenic differentiation was been analyzed. As a result, a total of 6,449 candidate circRNAs were identified in ovine preadipocytes. Of these circRNAs identified, 63 candidate circRNAs were differentially expressed among the three differentiation stages and their parent genes were mainly enriched in acetyl-CoA metabolic process, positive regulation of lipid biosynthetic process, positive regulation of steroid biosynthetic process, and focal adhesion pathway (P < 0.05). Based on a circRNAs-miRNAs-mRNAs regulatory network constructed, circ_004977, circ_006132 and circ_003788 were found to function as competing endogenous RNAs (ceRNAs) to regulate ovine preadipocyte differentiation and lipid metabolism. The results provide an improved understanding of functions and molecular mechanisms of circRNAs underlying ovine adipogenesis in sheep.

The moderate fat deposition contributes to improve mutton quality, which is associated with the differentiation of preadipocytes. To investigate roles of circular RNAs (circRNAs) in preadipocyte differentiation, we identified circRNAs on days 0, 2, and 8 of preadipocytes differentiation and compared the expression profile of circRNAs at different adipogenic differentiation stages. A total of 6,449 candidate circRNAs were identified, among which 63 candidate circRNAs were differentially expressed among the three differentiation stages. The parent genes of differentially expressed circRNAs were enriched in several biological process and pathways related to lipid metabolism and synthesis. In addition, several circRNAs may regulate ovine preadipocyte differentiation by interacting with microRNAs (miRNAs). The results reveal the potential roles of circRNAs in adipogenic differentiation of sheep.

Synergistic Responses of Tibetan Sheep Rumen Microbiota, Metabolites, and the Host to the Plateau Environment.

Plateau adaptation in animals involves genetic mechanisms as well as coevolutionary mechanisms of the microbiota and metabolome of the animal. Therefore, the characteristics of the rumen microbiome and metabolome, transcriptome, and serum metabolome of Tibetan sheep at different altitudes (4500 m, 3500 m, and 2500 m) were analyzed. The results showed that the rumen differential metabolites at 3500 m and 4500 m were mainly enriched in amino acid metabolism, lipid metabolism, and carbohydrate metabolism, and there was a significant correlation with microbiota. The differentially expressed genes and metabolites at middle and high altitudes were coenriched in asthma, arachidonic acid metabolism, and butanoate and propanoate metabolism. In addition, the serum differential metabolites at 3500 m and 4500 m were mainly enriched in amino acid metabolism, lipid metabolism, and metabolism of xenobiotics by cytochrome P450, and they were also related to microbiota. Further analysis revealed that rumen metabolites accounted for 7.65% of serum metabolites. These common metabolites were mainly enriched in metabolic pathways and were significantly correlated with host genes (p < 0.05). This study found that microbiota, metabolites, and epithelial genes were coenriched in pathways related to lipid metabolism, energy metabolism, and immune metabolism, which may be involved in the regulation of Tibetan sheep adaptation to plateau environmental changes.