Search

Abstract Background Superoxide dismutase (SOD) can greatly scavenge reactive oxygen species (ROS) in plants. SOD activity is highly related to plant stress tolerance that can be improved by overexpression of SOD genes. Identification of SOD activity-related loci and potential candidate genes is essential for improvement of grain quality in wheat breeding. However, the loci and candidate genes for relating SOD in wheat grains are largely unknown. In the present study, grain SOD activities of 309 recombinant inbred lines (RILs) derived from the ‘Berkut’ × ‘Worrakatta’ cross were assayed by photoreduction method with nitro-blue tetrazolium (NBT) in four environments. Quantitative trait loci (QTL) of SOD activity were identified using inclusive composite interval mapping (ICIM) with the genotypic data of 50 K single nucleotide polymorphism (SNP) array. Results Six QTL for SOD activity were mapped on chromosomes 1BL, 4DS, 5AL (2), and 5DL (2), respectively, explaining 2.2 ~ 7.4% of the phenotypic variances. Moreover, QSOD.xjau-1BL, QSOD.xjau-4DS, QSOD.xjau-5 A.1, QSOD.xjau-5 A.2, and QSOD.xjau-5DL.2 identified are likely to be new loci for SOD activity. Four candidate genes TraesCS4D01G059500, TraesCS5A01G371600, TraesCS5D01G299900, TraesCS5D01G343100LC, were identified for QSOD.xjau-4DS, QSOD.xjau-5AL.1, and QSOD.xjau-5DL.1 (2), respectively, including three SOD genes and a gene associated with SOD activity. Based on genetic effect analysis, this can be used to identify desirable alleles and excellent allele variations in wheat cultivars. Conclusion These candidate genes are annotated for promoting SOD production and inhibiting the accumulation of ROS during plant growth. Therefore, lines with high SOD activity identified in this study may be preferred for future wheat breeding.

The rumen microbiota—a symbiont to its host and consists of critical functional substances—plays a vital role in the animal body and represents a new perspective in the study of adaptive evolution in animals. This study used Slide Viewer slicing analysis system, gas chromatography, RT-qPCR and other technologies, as well as 16S and metabolomics determination methods, to measure and analyze the microstructure of rumen epithelium, rumen fermentation parameters, rumen transport genes, rumen microbiota and metabolites in Tibetan sheep and Hu sheep. The results indicate that the rumen nipple height and cuticle thickness of Tibetan sheep are significantly greater than those of Hu sheep (p

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.

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.

IntroductionHigh-yield and high-quality production of silage maize in cold regions is crucial for ensuring the sustainable development of livestock industry.MethodsThis study first conducted an experiment to select the optimized silage maize varieties and densities using a split-plot design. The tested maize varieties were Xuntian 3171, Xuntian 16, Xunqing 858, and Fengtian 12, with each variety planted at densities of 67,500, 79,500, and 90,000 plants ha-1. Following the variety and density selection, another experiment on optimizing nitrogen management for silage maize was carried out using a completely randomized design: no nitrogen fertilizer (T1), applying urea-N 320 kg ha-1 (T2), applying urea-N 240 kg ha-1 (T3), applying polymer-coated urea-N 240 kg N ha-1 (T4), and ratios of polymer-coated urea-N to urea-N at 9:1 (T5), 8:2 (T6), 7:3 (T7), and 6:4 (T8). T5-T8 all applied 240 kg N ha-1. The yield and quality of silage maize, nitrogen use efficiency and balance, and economic benefits were evaluated.ResultsResults showed that Xunqing 858 had significantly higher plant height (8.7%-22.6% taller than the other three varieties) and leaf area (30.9% larger than Xuntian 3171), resulting in yield 11.5%-51.6% higher than the other three varieties. All varieties achieved maximum yields at a planting density of 79,500 plants ha-1. Integrated management strategy 7 (T7: Xunqing 858, 79,500 plants ha-1, polymer-coated urea-N to urea-N ratio of 7:3) achieved the highest yield of 73.1 t ha-1, a 6.1%-58.1% increase over other treatments. This strategy also produced the highest crude protein (11.1%) and starch (19.1%) contents, and the lowest neutral detergent fiber content (50.6%), with economic benefits improved by 10.3%-97.8% compared to other strategies. Additionally, T7 improved nitrogen use efficiency by 15.4%-94.5%, reduced soil nitrate leaching by 4.4%-36.5%, and decreased nitrogen surplus by 7.0%-46.6%.Conclusion and discussionComprehensive analysis revealed that the integrated management strategy 7 significantly improved silage maize yield and quality in cold regions while enhancing nitrogen use efficiency and reducing the risk of nitrate leaching, aligning with green agriculture development requirements. These findings will provide vital theoretical insights and practical guidance for high-yield and high-quality silage maize production in cold regions worldwide.

Microbiota density plays an important role in maintaining host metabolism, immune function, and health, and age has a specific effect on the composition of intestinal microbiota. Therefore, the age-specific effects of age differences on the structure and function of the ileum microbiota in Tibetan sheep were investigated by determining the density of the ileum microbiota, the content of VFAs, and the expression levels of their transporter-related genes at different ages. The results showed that the contents of acetic acid and propionic acid in the ileum of Tibetan sheep in the 1.5-year-old group were significantly higher (p < 0.05) than those in other age groups, and that the contents of total VFAs were also significantly higher (p < 0.05) than those in other age groups. The relative densities of ileum Rf, Ra, and Fs were significantly higher in the 1.5-year-old group than in the other age groups (p < 0.05). The ileum epithelial VFAs transport-related genes AE2, MCT-4, and NHE1 had the highest expression in the 1.5-year-old group, and the expression of DRA was significantly lower in the 1.5-year-old group than in the 6-year-old group (p < 0.05). Correlation analysis showed that Cb, Sr, and Tb were significantly positively correlated with butyric acid concentration (p < 0.05) and negatively correlated with acetic acid, but the difference was not significant (p > 0.05); MCT-1, MCT-4, and AE2 were significantly positively correlated (p < 0.05) with acetic, propionic, and isobutyric acid concentrations; NHE1, NHE2, and MCT-4 were highly significantly positively correlated (p < 0.01) with Romboutsia and unclassified_Peptostreptococcaceae, while acetic acid was significantly positively correlated (p < 0.05) with NK4A214_group; Romboutsia, and unclassified_Peptostreptococcaceae were significantly positively correlated (p < 0.05). Therefore, compared with other ages, the 1.5-year-old Tibetan sheep had a stronger fermentation and metabolic capacity in the ileum under traditional grazing conditions on the plateau, which could provide more energy for Tibetan sheep during plateau acclimatization.

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.

IntroductionImprovement of root architecture is crucial to increasing nutrient acquisition.MethodsTwo pot experiments were conducted to investigate the effects of different concentrations of urea ammonium nitrate solution (UAN) and ammonium polyphosphate (APP) on lettuce root architecture and the relationship between roots and nitrogen (N) and phosphorus (P) absorption.ResultsThe results showed that lettuce yield, quality, and root architecture were superior in the APP4 treatment compared to other P fertilizer treatments. The N480 treatment (480 mg N kg-1 UAN) significantly outperformed other N treatments in terms of root length, root surface area, and root volume. There were significant quantitative relationships between root architecture indices and crop uptake of N and P. The relationships between P uptake and root length and root surface area followed power functions. Crop N uptake was significantly linearly related to the length of fine roots with a diameter of

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.

The Gannan yak, a superior livestock breed found on the Tibetan Plateau, exhibits significantly enhanced body size, weight, and growth performance in comparison to the Tianzhu white yak. MiRNAs play a pivotal role in regulating muscle growth by negatively modulating target genes. In this study, we found the average diameter, area, and length of myofibers in Gannan yaks were significantly higher than those of Tianzhu white yaks. Further, we focused on analyzing the longissimus dorsi muscle from both Gannan yaks and Tianzhu white yaks through transcriptome sequencing to identify differentially expressed (DE)miRNAs that influence skeletal muscle development. A total of 254 DE miRNAs were identified, of which 126 miRNAs were up-regulated and 128 miRNAs were down-regulated. GO and KEGG enrichment analysis showed that the target genes of these DE miRNAs were significantly enriched in signaling pathways associated with muscle growth and development. By constructing a DE miRNA- DE mRNA interaction network, we screened 18 key miRNAs, and notably, four of the candidates (novel-m0143-3p, novel-m0024-3p, novel-m0128-5p, and novel-m0026-3p) targeted six genes associated with muscle growth and development (DDIT4, ADAMTS1, CRY2, AKIRIN2, SIX1, and FOXO1). These findings may provide theoretical references for further studies on the role of miRNAs in muscle growth and development in Gannan yaks.

During the adaptive evolution of animals, the host and its gut microbiota co-adapt to different elevations. Currently, there are few reports on the rumen microbiota–hepato-intestinal axis of Tibetan sheep at different altitudes. Therefore, the purpose of this study was to explore the regulatory effect of rumen microorganism–volatile fatty acids (VFAs)–VFAs transporter gene interactions on the key enzymes and genes related to gluconeogenesis in Tibetan sheep. The rumen fermentation parameters, rumen microbial densities, liver gluconeogenesis activity and related genes were determined and analyzed using gas chromatography, RT-qPCR and other research methods. Correlation analysis revealed a reciprocal relationship among rumen microflora–VFAs-hepatic gluconeogenesis in Tibetan sheep at different altitudes. Among the microbiota, Ruminococcus flavefaciens (R. flavefaciens), Ruminococcus albus (R. albus), Fibrobactersuccinogenes and Ruminobacter amylophilus (R. amylophilus) were significantly correlated with propionic acid (p < 0.05), while propionic acid was significantly correlated with the transport genes monocarboxylate transporter 4 (MCT4) and anion exchanger 2 (AE2) (p < 0.05). Propionic acid was significantly correlated with key enzymes such as pyruvate carboxylase, phosphoenolpyruvic acid carboxylase and glucose (Glu) in the gluconeogenesis pathway (p < 0.05). Additionally, the expressions of these genes were significantly correlated with those of the related genes, namely, forkhead box protein O1 (FOXO1) and mitochondrial phosphoenolpyruvate carboxykinase 2 (PCK2) (p < 0.05). The results showed that rumen microbiota densities differed at different altitudes, and the metabolically produced VFA contents differed, which led to adaptive changes in the key enzyme activities of gluconeogenesis and the expressions of related genes.

The rumen microbiota and metabolites play an important role in energy metabolism and immune regulation of the host. However, the regulatory mechanism of rumen microbiota and metabolite interactions with host on Tibetan sheep’s plateau adaptability is still unclear. We analyzed the ruminal microbiome and metabolome, host transcriptome and serum metabolome characteristics of Tibetan sheep at different ages. Biomarkers Butyrivibrio, Lachnospiraceae_XPB1014_group, Prevotella, and Rikenellaceae_RC9_gut_group were found in 4 months, 1.5 years, 3.5 years, and 6 years Tibetan sheep, respectively. The rumen microbial metabolites were mainly enriched in galactose metabolism, unsaturated fatty acid biosynthesis and fatty acid degradation pathways, and had significant correlation with microbiota. These metabolites further interact with mRNA, and are co-enriched in arginine and proline metabolism, metabolism of xenobiotics by cytochrome P450, propanoate metabolism, starch and sucrose metabolism, gap junction pathway. Meanwhile, serum metabolites also have a similar function, such as chemical carcinogenesis − reactive oxygen species, limonene and pinene degradation, and cutin, suberine and wax biosynthesis, thus participating in the regulation of the body’s immune and energy-related metabolic processes. This study systematically revealed that rumen microbiota, metabolites, mRNA and serum metabolites of Tibetan sheep were involved in the regulation of fermentation metabolic function and immune level of Tibetan sheep at different ages, which provided a new perspective for plateau adaptability research of Tibetan sheep at different ages.

Background and aims: Defective enzymes, cofactors, or transporters of metabolic pathways cause inherited metabolic disorders (IMDs), a group of genetic disorders. Several IMDs have serious consequences for the affected neonates. Newborn screening for IMDs is conducted by measuring specific metabolites between 3 and 7 days of life. Herein, we analyzed the incidence, spectrum, and genetic characteristics of IMDs in newborns in the Zhuzhou area.Methods: Tandem mass spectrometry was conducted on 90,829 newborns who were admitted to the Women and Children Healthcare Hospital of Zhuzhou and requested for screening for IMDs. These newborns were subsequently subjected to next-generation sequencing and further validated using Sanger sequencing.Results: 30 IMDs cases were found in 90,829 cases of newborns screened for IMDs, and the overall incidence was 1/3,027. The incidence of amino acid, organic acid, fatty acid oxidation and urea cycle disorders were 1/8,257, 1/18,165, 1/7,569, and 1/45,414, respectively. Additionally, 9 cases of maternal IMDs were found in our study, and unreported gene mutations of 3 cases IMDs were identified.Conclusion: Our data indicated that IMDs are never uncommon in zhuzhou, meanwhile, we also found that primary carnitine deficiency was the only disorder of fatty acid oxidation in Zhuzhou, and the incidence (1/7,569) was higher than the national level, organic acid metabolic diseases are mostly inherited. Therefore, our study has clarified the disease spectrum and genetic backgrounds, contributing to the treatment and prenatal genetic counseling of these disorders in this region.

ABSTRACT Different reproductive stages of mammals involve complex biological processes, and the intestinal microbiota, as an endocrine organ or an “invisible organ,” is involved in the regulation of hormone levels, immune function, and metabolism. However, the effects of the rumen microbiota, its metabolites, and blood metabolites on the reproductive performance of ruminants remain unclear. This study revealed that the Prevotella abundance increased significantly during pregnancy (P < 0.01); the Fibrobacter abundance increased significantly during lactation (P < 0.05); and rumen microbial carbohydrate metabolism, glucose biosynthesis, and metabolic functions were significantly enriched during pregnancy (P < 0.05). Microbial metabolic profile analysis showed that the differentially abundant microbial metabolites during pregnancy and lactation were mainly enriched in the biosynthesis of ubiquinone and other terpenoid quinones, and there was a certain correlation with the microbiota. Among them, sapindoside A was increased during pregnancy, nicotinamide riboside and β-cryptoxanthin were reduced during pregnancy, and L-tryptophan was significantly increased during lactation. In addition, the volatile fatty acid levels in lactation were significantly higher than those in non-pregnancy and pregnancy (P < 0.05), and the NH3-N content during pregnancy was significantly higher than that during lactation and non-pregnancy (P < 0.05). Moreover, there were differences in the serum metabolite levels at different reproductive stages, and similar metabolites existed when comparing the rumen metabolites, which were mainly enriched in arachidonic acid metabolism, vitamin B6 metabolism, and ABC transporter protein, resulting in significantly higher serum IgA and IgM levels during lactation than during non-pregnancy and pregnancy (P < 0.05). IMPORTANCE Our study illustrates the succession of the rumen microbiota and its metabolites in Small-tailed Han sheep at different reproductive stages. Among them, Firmicutes and Prevotella, which are related to energy metabolism, increased in abundance during pregnancy, while Fibrobacter, a fiber-degrading bacterium, increased in abundance during lactation. At the same time, the microbial metabolic profile and serum metabolic profile characteristics of different reproductive stages were revealed, and some functional pathways and metabolites related to energy and immunity were found. This study provides a reference for the health management of ruminants during non-pregnancy, pregnancy, and lactation.

Unsaturated fatty acids (UFAs) in beef play a vital role in promoting human health. Long-chain fatty acyl-CoA synthase 1 (ACSL1) is a crucial gene for UFA synthesis in bovine adipocytes. To investigate the protein expression profile during UFA synthesis, we performed a proteomic analysis of bovine adipocytes by RNA interference and non-interference with ACSL1 using label-free techniques. A total of 3558 proteins were identified in both the NC and si-treated groups, of which 1428 were differentially expressed proteins (DEPs; fold change ≥ 1.2 or ≤ 0.83 and p-value < 0.05). The enrichment analysis of the DEPs revealed signaling pathways related to UFA synthesis or metabolism, including cAMP, oxytocin, fatty acid degradation, glycerol metabolism, insulin, and the regulation of lipolysis in adipocytes (p-value < 0.05). Furthermore, based on the enrichment analysis of the DEPs, we screened 50 DEPs that potentially influence the synthesis of UFAs and constructed an interaction network. Moreover, by integrating our previously published transcriptome data, this study established a regulatory network involving differentially expressed long non-coding RNAs (DELs), highlighting 21 DEPs and 13 DELs as key genes involved in UFA synthesis. These findings present potential candidate genes for further investigation into the molecular mechanisms underlying UFA synthesis in bovines, thereby offering insights to enhance the quality of beef and contribute to consumer health in future studies.

Jeryak is the F1 generation of the cross between Gannan yak and Jersey cattle, which has the advantages of fast growth and high adaptability. The growth and development of skeletal muscle is closely linked to meat production and the quality of meat. However, the molecular regulatory mechanisms of muscle growth differences between Gannan yak and Jeryak analyzed from the perspective of chromatin opening have not been reported. In this study, ATAC-seq was used to analyze the difference of chromatin openness in longissimus muscle of Gannan yak and Jeryak. It was found that chromatin accessibility was more enriched in Jeryak compared to Gannan yak, especially in the range of the transcription start site (TSS) ± 2 kb. GO and KEGG enrichment analysis indicate that differential peak-associated genes are involved in the negative regulation of muscle adaptation and the Hippo signaling pathway. Integration analysis of ATAC-seq and RNA-seq revealed overlapping genes were significantly enriched during skeletal muscle cell differentiation and muscle organ morphogenesis. At the same time, we screened FOXO1, ZBED6, CRY2 and CFL2 for possible involvement in skeletal muscle development, constructed a genes and transcription factors network map, and found that some transcription factors (TFs), including YY1, KLF4, KLF5 and Bach1, were involved in skeletal muscle development. Overall, we have gained a comprehensive understanding of the key factors that impact skeletal muscle development in various breeds of cattle, providing new insights for future analysis of the molecular regulatory mechanisms involved in muscle growth and development.

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.

The gut microbiota is highly symbiotic with the host, and the microbiota and its metabolites are essential for regulating host health and physiological functions. Astragalus, as a feed additive, can improve animal immunity. However, the effects of Astragalus root powder on the rumen microbiota and their metabolites in lambs are not apparent. In this study, thirty healthy Hu sheep lambs with similar body weights (17.42 ± 2.02 kg) were randomly selected for the feeding experiment. Lambs were fed diets supplemented with 0.3% Astragalus root powder, and the rumen microbiota density and metabolome were measured to determine the effects of Astragalus on the health of lambs in the rumen. The results showed that the relative abundance of Butyrivibrio fibrisolvens (Bf), Ruminococcus flavefaciens (Rf), Succiniclasticum (Su), and Prevotella (Pr) in the rumen was increased in the Astragalus group (p < 0.01), and metabolic profiling showed that the metabolites, such as L-lyrosine and L-leucine, were upregulated in the Astragalus group (p < 0.01). KEGG functional annotation revealed that upregulated metabolites were mainly enriched in the pathways of amino acid metabolism, lipid metabolism, fatty acid biosynthesis, and bile secretion in the Astragalus group, and downregulated metabolites were enriched in the pathways of methane metabolism and other pathways. Correlation analysis revealed that butyric acid was positively correlated with Roseburia and Blautia (p < 0.05) and negatively correlated with Desulfovibrio (p < 0.05). Thus, by analyzing the interactions of Astragalus root powder with the density of rumen microorganisms and their metabolites in lambs, it was shown that Astragalus root powder could improve the structure of rumen microbiota and their metabolites and then participate in the regulation of amino acid metabolism, lipid metabolism, immune metabolism, and other pathways to improve the efficiency of energy absorption of the lambs.

The intestinal microbiota of ruminants is an important factor affecting animal production and health. Research on the association mechanism between the intestinal microbiota and meat quality of ruminants will play a positive role in understanding the formation mechanism of meat quality in ruminants and improving production efficiency. In this study, the fatty acid composition and content, expression of related genes, and structural characteristics of the ileum microbiota of ewes of Tibetan sheep at different ages (4 months, 1.5 years, 3.5 years, and 6 years) were detected and analyzed. The results revealed significant differences in fatty acid composition and content in the muscle of Tibetan sheep at different ages (p < 0.05); in addition, the content of MUFAs in the longissimus dorsi muscle and leg muscle was higher. Similarly, the expressions of muscle-related genes differed among the different age groups, and the expression of the LPL, SCD, and FABP4 genes was higher in the 1.5-year-old group. The ileum microbiota diversity was higher in the 1.5-year-old group, the Romboutsia abundance ratio was significantly higher in the 1.5-year-old group (p < 0.05), and there was a significant positive correlation with oleic acid (C18:1n9c) (p < 0.05). In conclusion, the content of beneficial fatty acids in the longissimus dorsi muscle and leg muscle of Tibetan sheep was higher at 1.5 years of age, and the best slaughter age was 1.5 years. This study provides a reference for in-depth research on the mechanism of the influence of the gut microbiota on meat quality and related regulation.

Muscle formation directly determines meat production and quality. The non-SMC condensin I complex subunit G (NCAPG) is strongly linked to the growth features of domestic animals because it is essential in controlling muscle growth and development. This study aims to elucidate the tissue expression level of the bovine NCAPG gene, and determine the key transcription factors for regulating the bovine NCAPG gene. In this study, we observed that the bovine NCAPG gene exhibited high expression levels in longissimus dorsi and spleen tissues. Subsequently, we cloned and characterized the promoter region of the bovine NCAPG gene, consisting of a 2039 bp sequence, through constructing the deletion fragment double-luciferase reporter vector and site-directed mutation-identifying core promoter region with its key transcription factor binding site. In addition, the key transcription factors of the core promoter sequence of the bovine NCAPG gene were analyzed and predicted using online software. Furthermore, by integrating overexpression experiments and the electrophoretic mobility shift assay (EMSA), we have shown that cAMP response element binding protein 1 (CREB1) and myogenic differentiation 1 (MYOD1) bind to the core promoter region (−598/+87), activating transcription activity in the bovine NCAPG gene. In conclusion, these findings shed important light on the regulatory network mechanism that underlies the expression of the NCAPG gene throughout the development of the muscles in beef cattle.

The FOXO3 gene, a prominent member of the FOXO family, has been identified as a potential quantitative trait locus for muscle atrophy and lipid metabolism in livestock. It is also considered a promising candidate gene for meat quality traits such as Warner–Bratzler shear force (WBSF) and water holding capacity (WHC). The aim of this study was to identify sequence mutations in the FOXO3 gene of yaks and to analyze the association of genotypes and haplotypes with meat traits such as WBSF and WHC. Quantitative reverse-transcriptase PCR (RT-qPCR) was applied to determine the expression levels of FOXO3 in yak tissues, with the results revealing a high expression in the yak longissimus dorsi muscle. Exons of the FOXO3 gene were then sequenced in 572 yaks using hybrid pool sequencing. Five single nucleotide polymorphisms were identified. Additionally, four effective haplotypes and four combined haplotypes were constructed. Two mutations of the FOXO3 gene, namely C>G at exon g.636 and A>G at exon g.1296, were associated with cooked meat percentage (CMP) (p < 0.05) and WBSF (p < 0.05), respectively. Furthermore, the WBSF of the H2H3 haplotype combination was significantly lower than that of other combinations (p < 0.05). The findings of this study suggest that genetic variations in FOXO3 could be a promising biomarker for improving yak meat traits.

Homobox C13 (Hoxc13) is an important transcription factor in hair follicle cycle development, and its deletion had been found in a variety of animals leading to abnormal hair growth and disruption of the hair follicle system. In this study, we used immunofluorescence, immunohistochemistry, real-time fluorescence quantitative PCR (RT-qPCR), and Kompetitive Allele-Specific PCR (KASP) genotyping to investigate molecular genetic characteristics of the Hoxc13 gene in Gansu alpine fine-wool sheep. The results revealed that Hoxc13 was significantly expressed during both the anagen and catagen phases (p < 0.05). It was found to be highly expressed predominantly in the dermal papillae and the inner and outer root sheaths, showing a distinct spatiotemporal expression pattern. Two single nucleotide polymorphisms (SNPs) in the exon 1 of Hoxc13, both the individual locus genotypes and the combined haplotypes were found to be correlated with wool length (p < 0.05). It was determined the mutations led to changes in mRNA expression, in which higher expression of this gene was related with longer wool length. In summary, this unique spatiotemporal expression pattern of the Hoxc13 gene may regulate the wool length of Gansu alpine fine-wool sheep, which can be used as a molecular genetic marker for wool traits and thus improve the breed.

During the postnatal stages, skeletal muscle development undergoes a series of meticulously regulated alterations in gene expression. However, limited studies have employed chromatin accessibility to unravel the underlying molecular mechanisms governing muscle development in yak species. Therefore, we conducted an analysis of both gene expression levels and chromatin accessibility to comprehensively characterize the dynamic genome-wide chromatin accessibility during muscle growth and development in the Tianzhu white yak, thereby elucidating the features of accessible chromatin regions throughout this process. Initially, we compared the differences in chromatin accessibility between two groups and observed that calves exhibited higher levels of chromatin accessibility compared to adult cattle, particularly within ±2 kb of the transcription start site (TSS). In order to investigate the correlation between alterations in chromatin accessible regions and variations in gene expression levels, we employed a combination of ATAC-seq and RNA-seq techniques, leading to the identification of 18 central transcriptional factors (TFs) and 110 key genes with significant effects. Through further analysis, we successfully identified several TFs, including Sp1, YY1, MyoG, MEF2A and MEF2C, as well as a number of candidate genes (ANKRD2, ANKRD1, BTG2 and LMOD3) which may be closely associated with muscle growth and development. Moreover, we constructed an interactive network program encompassing hub TFs and key genes related to muscle growth and development. This innovative approach provided valuable insights into the molecular mechanism underlying skeletal muscle development in the postnatal stages of Tianzhu white yaks while also establishing a solid theoretical foundation for future research on yak muscle development.

The Tibetan sheep has an intricate mechanism of adaptation to low oxygen levels, which is influenced by both genetic and environmental factors. The heart plays a crucial role in the adaptation of Tibetan sheep to hypoxia. In the present study, we utilized transcriptomic and proteomic technologies to comprehensively analyze and identify the long non-coding RNAs (lncRNAs), genes, proteins, pathways, and gene ontology (GO) terms associated with hypoxic adaptation in Tibetan sheep at three different altitudes (2500 m, 3500 m, and 4500 m). By integrating the differentially expressed (DE) lncRNA target genes, differentially expressed proteins (DEPs), and differentially expressed genes (DEGs), we were able to identify and characterize the mechanisms underlying hypoxic adaptation in Tibetan sheep. Through this integration, we identified 41 shared genes/proteins, and functional enrichment analyses revealed their close association with lipid metabolism, glycolysis/gluconeogenesis, and angiogenesis. Additionally, significant enrichment was observed in important pathways such as the PPAR signaling pathway, glycolysis/gluconeogenesis, the oxoacid metabolic process, and angiogenesis. Furthermore, the co-expression network of lncRNAs and mRNAs demonstrated that lncRNAs (MSTRG.4748.1, ENSOART00020025894, and ENSOART00020036371) may play a pivotal role in the adaptation of Tibetan sheep to the hypoxic conditions of the plateau. In conclusion, this study expands the existing database of lncRNAs and proteins in Tibetan sheep, and these findings may serve as a reference for the prevention of altitude sickness in humans.

Hormone-sensitive lipase (HSL) is involved in the breakdown of triacylglycerols in adipose tissue, which influences muscle tenderness and juiciness by affecting the intramuscular fat content (IMF). This study analyzed the association between different genotypes and haplotypes of the yak HSL gene and carcass and meat quality traits. We used hybridization pool sequencing to detect exon 2, exon 8, and intron 3 variants of the yak HSL gene and genotyped 525 Gannan yaks via KASP to analyze the effects of the HSL gene variants on the carcass and meat quality traits in yaks. According to the results, the HSL gene is highly expressed in yak adipose tissue. Three single nucleotide polymorphisms (SNPs) were identified, with 2 of them located in the coding region and one in the intron region. Variants in the 2 coding regions resulted in amino acid changes. The population had 3 genotypes of GG, AG, and AA, and individuals with the AA genotype had lower WBSF values (p < 0.05). The H3H3 haplotype combinations could improve meat tenderness by reducing the WBSF values and the cooking loss rate (CLR) (p < 0.05). H1H1 haplotype combinations were associated with the increased drip loss rate (DLR) (p < 0.05). The presence of the H1 haplotype was associated the increased CLR in yaks, while that of the H2 haplotype was associated with the decreased DLR in yaks (p < 0.05). These results demonstrated that the HSL gene may influence the meat quality traits in yaks by affecting the IMF content in muscle tissues. Consequently, the HSL gene can possibly be used as a biomarker for improving the meat quality traits in yaks in the future.

ObjectiveThis study aims to observe the changes of the umbilical venous–arterial index (VAI) and investigate its predictive power for fetal outcome during the second half of pregnancy.MethodsFetuses with gestational age (GA) at 24–39 weeks were collected. According to the outcome score, neonates with outcome scores of 0, 1, or 2 were assigned to the control group, whereas those with scores of 3–12 were assigned to the compromised group. VAI was calculated as the ratio of normalized umbilical vein blood flow volume and umbilical artery pulsatility index. Regression analysis was performed to obtain the best-fitting curves between VAI and GA in the controls. Doppler parameters and perinatal outcomes were compared in both groups. Receiver operating characteristic analysis was used to assess the diagnostic performance of the VAI.ResultsA total of 833 (95%) fetuses had Doppler parameters and pregnancy outcomes documented. Compared with the controls, the VAI was significantly lower in the compromised group (83.2 vs. 184.8 ml/min/kg, p

Abstract Background The gut microbiota play an important role in maintaining host metabolism, the immune system and health, while sex, genotype, diet and health have specific effects on the composition of the gut microbiota. Therefore, to explore the sex differences in the structure and function of rumen microbiota in Tibetan goats, herein we analyzed sex differences in rumen fermentation parameters, rumen microbiota and the expression of genes related to VFA transport in Tibetan goats. Results The results showed that the contents of acetic acid and propionic acid in the rumen of TGM (Tibetan goat male) were significantly higher than those in TGFm (Tibetan goat female) (P

The production performance of Jeryak, resulting from the F1 generation of the cross between Gannan yak and Jersey cattle, exhibits a significantly superior outcome compared with that of Gannan yak. Therefore, we used an RNA-seq approach to identify differentially expressed mRNAs (DEMs) and differentially expressed lncRNAs (DELs) influencing muscle growth and development in Gannan yaks and Jeryaks. A total of 304 differentially expressed lncRNAs and 1819 differentially expressed mRNAs were identified based on the screening criteria of |log 2 FC| > 1 and FDR < 0.05. Among these, 132 lncRNAs and 1081 mRNAs were found to be down-regulated, while 172 lncRNAs and 738 mRNAs were up-regulated. GO and KEGG analyses showed that the identified DELs and DEMs were enriched in the entries of pathways associated with muscle growth and development. On this basis, we constructed an lncRNA–mRNA interaction network. Interestingly, two candidate DELs (MSTRG.16260.9 and MSTRG.22127.1) had targeting relationships with 16 (MYC, IGFBP5, IGFBP2, MYH4, FGF6, etc.) genes related to muscle growth and development. These results could provide a basis for further studies on the roles of lncRNAs and mRNAs in muscle growth in Gannan yaks and Jeryak breeds.

Based on the serious phenomenon of Tibetan sheep “growing strong in warm seasons and losing weight in cold seasons”, this study explores the regulation of lncRNAs, microbiota, and metabolites in the cold season adaptation of Tibetan sheep from the perspective of the co-evolution of the host genome (first genome) and microbiome (second genome). RNA-seq results showed that 172 DE lncRNAs were identified in the rumen epithelium of Tibetan sheep in warm and cold seasons, of which 87 DE lncRNAs were significantly up-regulated in cold seasons. KEGG enrichment showed that target genes of up-regulated lncRNAs were significantly enriched in TNF signaling and oxidative phosphorylation pathways. LncRNA-mRNA regulatory network indicated that DE lncRNAs were involved in nutrient stress in the cold season by targeting ATP1B2, CADPS, TLR5, and UGT1A6. Correlation analysis showed some lncRNAs were significantly correlated with acetic acid, propionic acid, butyric acid, and rumen epithelial histomorphology and had a negative correlation with Butyrivibrio-2 and Succiniclasticum (p < 0.05). In addition, differential metabolites bilirubin and lncRNAs were co-enriched in the bile secretion pathway. lncRNAs played an important role in the adaptation process of Tibetan sheep in the cold season, and mediate the host to participate in nutrient absorption, energy utilization, and immune response, indicating that the host genome and microbial genome promote Tibetan sheep to adapt to nutrient stress in the cold season through co-evolution.

Acetyl-CoA carboxylase beta (ACACB) is a functional candidate gene that impacts fat deposition. In the present study, we sequenced exon 37–intron 37, exon 46–intron 46, and intron 47 of yak ACACB using hybrid pool sequencing to search for variants and genotyped the gene in 593 Gannan yaks via Kompetitive allele-specific polymerase chain (KASP) reaction to determine the effect of ACACB variants on carcass and meat quality traits. Seven single nucleotide polymorphisms were detected in three regions. Eight effective haplotypes and ten diplotypes were constructed. Among them, a missense variation g.50421 A > G was identified in exon 37 of ACACB, resulting in an amino acid shift from serine to glycine. Correlation analysis revealed that this variation was associated with the cooking loss rate and yak carcass weight (p = 0.024 and 0.012, respectively). The presence of haplotypes H5 and H6 decreased Warner–Bratzler shear force (p = 0.049 and 0.006, respectively), whereas that of haplotypes H3 and H4 increased cooking loss rate and eye muscle area (p = 0.004 and 0.034, respectively). Moreover, the presence of haplotype H8 decreased the drip loss rate (p = 0.019). The presence of one and two copies of haplotypes H1 and H8 decreased the drip loss rate (p = 0.028 and 0.004, respectively). However, haplotype H1 did not decrease hot carcass weight (p = 0.011), whereas H3 increased the cooking loss rate (p = 0.007). The presence of one and two copies of haplotype H6 decreased Warner–Bratzler shear force (p = 0.014). The findings of the present study suggest that genetic variations in ACACB can be a preferable biomarker for improving yak meat quality.

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.

Tibetan sheep are already well adapted to cold season nutrient stress on the Tibetan Plateau. Rumen, an important nutrient for metabolism and as an absorption organ in ruminants, plays a vital role in the cold stress adaptations of Tibetan sheep. Ruminal microbiota also plays an indispensable role in rumen function. In this study, combined multiomics data were utilized to comprehensively analyze the interaction mechanism between rumen epithelial miRNAs and microbiota and their metabolites in Tibetan sheep under nutrient stress in the cold season. A total of 949 miRNAs were identified in the rumen epithelium of both cold and warm seasons. A total of 62 differentially expressed (DE) miRNAs were screened using FC > 1.5 and p value < 0.01, and a total of 20,206 targeted genes were predicted by DE miRNAs. KEGG enrichment analysis revealed that DE miRNA-targeted genes were mainly enriched in axon guidance(ko04360), tight junction(ko04530), inflammatory mediator regulation of TRP channels(ko04750) and metabolism-related pathways. Correlation analysis revealed that rumen microbiota, rumen VFAs and DE miRNAs were all correlated. Further study revealed that the targeted genes of cold and warm season rumen epithelial DE miRNAs were coenriched with differential metabolites of microbiota in glycerophospholipid metabolism (ko00564), apoptosis (ko04210), inflammatory mediator regulation of TRP channels (ko04750), small cell lung cancer (ko05222), and choline metabolism in cancer (ko05231) pathways. There are several interactions between Tibetan sheep rumen epithelial miRNAs, rumen microbiota, and microbial metabolites, mainly through maintaining rumen epithelial barrier function and host homeostasis of choline and cholesterol, improving host immunity, and promoting energy metabolism pathways, thus enabling Tibetan sheep to effectively respond to cold season nutrient stress. The results also suggest that rumen microbiota have coevolved with their hosts to improve the adaptive capacity of Tibetan sheep to cold season nutrient stress, providing a new perspective for the study of cold season nutritional stress adaptation in Tibetan sheep.

The Cashmere goat is an excellent local goat breed in Gansu Province of China, and it is expected to improve cashmere production and cashmere quality through selection and breeding to enhance its commercial value. Keratin-associated proteins (KAPs) play an important role in maintaining wool structure. The gene encoding the keratin-associated protein 22-2 (KAP22-2) gene has been identified in selected species other than goats, such as humans, mice, and sheep. In this study, the sequence of the sheep KAP22-2 gene (KRTAP22-2) was aligned into the goat genome, and the sequence with the highest homology was assumed to be the goat KRTAP22-2 sequence and used to design primers to amplify the goat gene sequence. A total of 356 Longdong Cashmere goats (Gansu Province, China) were used for screening of genetic variants. Four specific bands were detected by polymerase chain reaction-single-stranded conformational polymorphism (PCR-SSCP) analysis, and they formed a total of six band types individually or in combination. Four alleles were identified by DNA sequencing of PCR amplification products. A total of four single nucleotide polymorphic sites (SNPs) were detected in the four sequenced KRTAP22-2 alleles. Two of them are in the 5’UTR region and the other two are in the coding region, and the variants in the coding region are all non-synonymous mutations. In addition, there was a 6 bp length variation in allele C. The gene was expressed in the cortical layer of primary and secondary hair follicles, the inner root sheath, as well as hair papillae and hair maternal cells in goats. The results of the correlation analysis between genotypes and cashmere traits showed that after excluding genotypes with a gene frequency of less than 5%, the mean fiber diameter (MFD) of cashmere was significantly higher in the AB genotype than in the AA and AC genotypes. That is, the KRTAP22-2 gene variants are associated with mean fiber diameter in cashmere. The above results suggest that the goat KRTAP22-2 variant can be utilized as a molecular marker candidate gene for cashmere traits.

Long-chain fatty acyl-CoA synthase 1 (ACSL1) plays a vital role in the synthesis and metabolism of fatty acids. The proportion of highly unsaturated fatty acids in beef not only affects the flavor and improves the meat’s nutritional value. In this study, si-ACSL1 and NC-ACSL1 were transfected in bovine preadipocytes, respectively, collected cells were isolated on the fourth day of induction, and then RNA-Seq technology was used to screen miRNAs related to unsaturated fatty acid synthesis. A total of 1,075 miRNAs were characterized as differentially expressed miRNAs (DE-miRNAs), of which the expressions of 16 miRNAs were upregulated, and that of 12 were downregulated. Gene ontology analysis indicated that the target genes of DE-miRNAs were mainly involved in biological regulation and metabolic processes. Additionally, KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis identified that the target genes of DE-miRNAs were mainly enriched in metabolic pathways, fatty acid metabolism, PI3K-Akt signaling pathway, glycerophospholipid metabolism, fatty acid elongation, and glucagon signaling pathway. Combined with the previous mRNA sequencing results, several key miRNA-mRNA targeting relationship pairs, i.e., novel-m0035-5p—ACSL1, novel-m0035-5p—ELOVL4, miR-9-X—ACSL1, bta-miR-677—ACSL1, miR-129-X—ELOVL4, and bta-miR-485—FADS2 were screened via the miRNA-mRNA interaction network. Thus, the results of this study provide a theoretical basis for further research on miRNA regulation of unsaturated fatty acid synthesis in bovine adipocytes.

In our previous study, microRNA (miR)-199a-3p was found to be the most upregulated miRNA in mammary gland tissue during the non-lactation period compared with the peak-lactation period. However, there have been no reports describing the function of miR-199a-3p in ovine mammary epithelial cells (OMECs) and the biological mechanisms by which the miRNA affects cell proliferation and milk fat synthesis in sheep. In this study, the effect of miR-199a-3p on viability, proliferation, and milk fat synthesis of OMECs was investigated, and the target relationship of the miRNA with very low-density lipoprotein receptor (VLDLR) was also verified. Transfection with a miR-199a-3p mimic increased the viability of OMECs and the number of Edu-labeled positive OMECs. In contrast, a miR-199-3p inhibitor had the opposite effect with the miR-199a-3p mimic. The expression levels of three marker genes were also regulated by both the miR-199a-3p mimic and miR-199-3p inhibitor in OMECs. Together, these results suggest that miR-199a-3p promotes the viability and proliferation of OMECs. A dual luciferase assay confirmed that miR-199a-3p can target VLDLR by binding to the 3′-untranslated regions (3'UTR) of the gene. Further studies found a negative correlation in the expression of miR-199a-3p with VLDLR. The miR-199a-3p mimic decreased the content of triglycerides, as well as the expression levels of six milk fat synthesis marker genes in OMECs, namely, lipoprotein lipase gene (LPL), acetyl-CoA carboxylase alpha gene (ACACA), fatty acid binding protein 3 gene (FABP3), CD36, stearoyl-CoA desaturase gene (SCD), and fatty acid synthase gene (FASN). The inhibition of miR-199a-3p increased the level of triglycerides and the expression of LPL, ACACA, FABP3, SCD, and FASN in OMECs. These findings suggest that miR-199a-3p inhibited milk fat synthesis of OMECs. This is the first study to reveal the molecular mechanisms by which miR-199a-3p regulates the proliferation and milk fat synthesis of OMECs in sheep.

Tibetan sheep are mainly distributed in the Qinghai–Tibet Plateau. Its meat is not only essential for the local people but also preferred by the non-inhabitant of this plateau also. To investigate the salient development features and molecular mechanism of the meat difference of LT muscle caused by different growth stages in Tibetan sheep, the carcass performance, meat quality, and comparative transcriptome analysis were performed for investigating the potential molecular mechanism of the meat quality difference of the LT muscle caused by four growth stages [4-months old (4 months), 1.5-years old (1.5 years), 3.5-years old (3.5 years), and 6-years old (6 years)] in the Tibetan sheep. The shear force increased with the increase of age (p < 0.05) while the intramuscular fat (IMF) was the highest at 1.5 y. The AMPK signaling pathway was significantly enriched in the four comparative groups. The weighted gene co-expression network analysis (WGCNA) results showed that the hub genes P4HA2, FBXL4, and PPARA were identified to regulate the meat quality. In summary, 1.5 years was found to be the most suitable slaughter age of the Tibetan sheep which ensured better meat tenderness and higher IMF content. Moreover, the genes LIPE, LEP, ADIPOQ, SCD, and FASN may regulate the transformation of the muscle fiber types through the AMPK signaling pathway, further affecting the meat quality.

Circular RNAs (circRNAs) are a kind of non-coding RNA that have an important molecular function in mammary gland development and lactation of mammals. In our previous study, circ_015343 was found to be highly expressed in the ovine mammary gland tissue at the peak-lactation period by using RNA sequencing (RNA-seq). In the present study, the authenticity of circ_015343 was confirmed by using reverse transcriptase-polymerase chain reaction (RT-PCR) analysis and Sanger sequencing. The circ_015343 was derived from the complete 10 exons of aminoadipic semialdehyde synthase (AASS), ranging from exon 2 to exon 11 and mainly located in cytoplasm of ovine mammary epithelial cells. The circRNA was found to be expressed in eight ovine tissues, with the highest expression level in the mammary gland and the least expression in Longissimus dorsi muscle. The circ_015343 had a lower level of expression in a sheep breed with higher milk yield and milk fat content. The disturbed circ_015343 increased the viability and proliferation of the ovine mammary epithelial cells. The inhibition of circ_015343 also increased the expression levels of three milk fat synthesis marker genes: acetyl-coenzyme A carboxylase alpha (ACACA), fatty acid-binding protein 4 (FABP4), and sterol regulatory element-binding protein 1 (SREBP1), as well as three proliferation-related genes: cyclin dependent kinase 2 (CDK2), cyclin dependent kinase 4 (CDK4) and proliferating cell nuclear antigen (PCNA), but decreased the expression level of its parent gene AASS. A circRNA-miRNA-mRNA interaction network showed that circ_015343 would bind some microRNAs (miRNAs) to regulate the expression of functional genes related to the development of mammary gland and lactation. This study contributes to a better understanding of the roles of circ_015343 in the mammary gland of sheep.

MicroRNAs (miRNAs) are a class of small non-coding RNAs that have been shown to play important post-transcriptional regulatory roles in the growth and development of skeletal muscle tissues. However, limited research into the effect of miRNAs on muscle development in goats has been reported. In this study, Liaoning cashmere (LC) goats and Ziwuling black (ZB) goats with significant phenotype difference in meat production performance were selected and the difference in Longissimus dorsi muscle tissue expression profile of miRNAs between the two goat breeds was then compared using small RNA sequencing. A total of 1,623 miRNAs were identified in Longissimus dorsi muscle tissues of the two goat breeds, including 410 known caprine miRNAs, 928 known species-conserved miRNAs and 285 novel miRNAs. Of these, 1,142 were co-expressed in both breeds, while 230 and 251 miRNAs were only expressed in LC and ZB goats, respectively. Compared with ZB goats, 24 up-regulated miRNAs and 135 miRNAs down-regulated were screened in LC goats. A miRNA-mRNA interaction network showed that the differentially expressed miRNAs would target important functional genes associated with muscle development and intramuscular fat deposition. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that the target genes of differentially expressed miRNAs were significantly enriched in Ras, Rap 1, FoxO, and Hippo signaling pathways. This study suggested that these differentially expressed miRNAs may be responsible for the phenotype differences in meat production performance between the two goat breeds, thereby providing an improved understanding of the roles of miRNAs in muscle tissue of goats.

The nuclear factor Y (NF-Y) transcription factor contains three subfamilies: NF-YA, NF-YB, and NF-YC. The NF-Y family have been reported to be key regulators in plant growth and stress responses. However, little attention has been given to these genes in melon (Cucumis melo L.). In this study, twenty-five NF-Ys were identified in the melon genome, including six CmNF-YAs, eleven CmNF-YBs, and eight CmNF-YCs. Their basic information (gene location, protein characteristics, and subcellular localization), conserved domains and motifs, and phylogeny and gene structure were subsequently analyzed. Results showed highly conserved motifs exist in each subfamily, which are distinct between subfamilies. Most CmNF-Ys were expressed in five tissues and exhibited distinct expression patterns. However, CmNF-YA6, CmNF-YB1/B2/B3/B8, and CmNF-YC6 were not expressed and might be pseudogenes. Twelve CmNF-Ys were induced by cold stress, indicating the NF-Y family plays a key role in melon cold tolerance. Taken together, our findings provide a comprehensive understanding of CmNF-Y genes in the development and stress response of melon and provide genetic resources for solving the practical problems of melon production.

Diacylglycerol acyltransferase-1 (DGAT1) has been recognised as one of the functional genes during fat deposition and closely related to meat quality. It is well known that yak meat has low intramuscular fat (IMF) and poor tenderness, influencing the purchase decision for consumers. In this study, the variations of DGAT1 gene were detected using polymerase chain reaction single-stranded conformational polymorphism (PCR-SSCP) analysis so as to ascertain the association of DGAT1 variation with yak carcase and meat quality traits. Two variants (A1 and B1) in intron 1 – exon 2 and another three variants (A2 to C2) in intron 15 – exon 17 were identified in the yaks investigated. Three sequence variations (c.192-117 G > A, c.1252-23 C > T and c.1339 C > T) were detected among these variants, with c.192-117 G > A located in intron 1, c.1252-23 C > T located in intron 15 and c.1339 C > T located in exon 17. The mutation c.1339 C > T resulted in an amino acid change (p. Arg 447 Cys). Variants A1 and B2 were associated with a decrease and increase in Warner–Bratzler shear force (WBSF) (p = 0.027 and p

The Tibetan sheep is an indigenous species on the Tibetan plateau with excellent adaptability to high-altitude hypoxia and is distributed at altitudes of 2500–5000 m. The high-altitude hypoxia adaptation of Tibetan sheep requires adaptive reshaping of multiple tissues and organs, especially the lungs. To reveal the mechanisms of adaptation at the tissue and molecular levels in the lungs of Tibetan sheep under hypoxic conditions at different altitudes, we performed light and electron microscopic observations, transcriptomic sequencing, and enzyme-linked immunosorbent assay studies on the lungs of Tibetan sheep from three altitudes (2500, 3500, and 4500 m). The results showed that in addition to continuous increase in pulmonary artery volume, thickness, and elastic fiber content with altitude, Tibetan sheep increase the hemoglobin concentration at an altitude of 3500 m, while they decrease the Hb concentration and increase the surface area of gas exchange and capacity of the blood at an altitude of 4500 m. Other than that, some important differentially expressed genes related to angiogenesis (FNDC1, HPSE, and E2F8), vasomotion and fibrogenesis (GJA4, FAP, COL1A1, COL1A2, COL3A1, and COL14A1), and gas transport (HBB, HBA1, APOLD1, and CHL1) were also identified; these discoveries at the molecular level explain to some extent the physiological findings. In conclusion, the lungs of Tibetan sheep adopt different strategies when adapting to different altitudes, and these findings are valuable for understanding the basis of survival of indigenous species on the Tibetan plateau.

Circular RNAs (circRNAs) have a regulatory role in animal skeletal muscle development. In this study, RNA sequencing was performed to reveal the temporal regularity of circRNA expression and the effect of the circRNA–miRNA–mRNA ceRNA regulatory network on the meat quality of longissimus thoracis (LT) muscle in Tibetan sheep at different growth stages (4 months old, 4 m; 1.5 years old, 1.5 y; 3.5 years old, 3.5 y; 6 years old, 6 y). There were differences in the carcass performance and meat quality of Tibetan sheep at different ages. Especially, the meat tenderness significantly decreased (p < 0.05) with the increase of age. GO functional enrichment indicated that the source genes of the DE circRNAs were mainly involved in the protein binding, and myofibril and organelle assembly. Moreover, there was a significant KEGG enrichment in the adenosine 5′-monophosphate (AMP)-activated protein kinase (AMPK) signaling pathway, as well as the calcium signaling pathway, regulating the pluripotency of the stem cells. The circRNA–miRNA–mRNA ceRNA interaction network analysis indicated that circRNAs such as circ_000631, circ_000281, and circ_003400 combined with miR-29-3p and miR-185-5p regulate the expression of LEP, SCD, and FASN related to the transformation of muscle fiber types in the AMPK signaling pathway. The oxidized muscle fibers were transformed into the glycolytic muscle fibers with the increase of age, the content of intramuscular fat (IMF) was lowered, and the diameter of the muscle fiber was larger in the glycolytic muscle fibers, ultimately increasing the meat tenderness. The study revealed the role of the circRNAs in the transformation of skeletal muscle fiber types in Tibetan sheep and its influence on meat quality. It improves our understanding of the role of circRNAs in Tibetan sheep muscle development.

Tibetan sheep can maintain a normal life and reproduce in harsh environments under extreme cold and lack of nutrition. However, the molecular and metabolic mechanisms underlying the adaptability of Tibetan sheep during the cold season are still unclear. Hence, we conducted a comprehensive analysis of rumen epithelial morphology, epithelial transcriptomics, microbiology and metabolomics in a Tibetan sheep model. The results showed that morphological structure of rumen epithelium of Tibetan sheep in cold season had adaptive changes. Transcriptomics analysis showed that the differential genes were primarily enriched in the PPAR signaling pathway (ko03320), legionellosis (ko05134), phagosome (ko04145), arginine and proline metabolism (ko00330), and metabolism of xenobiotics by cytochrome P450 (ko00980). Unique differential metabolites were identified in cold season, such as cynaroside A, sanguisorbin B and tryptophyl-valine, which were mainly enriched in arachidonic acid metabolism, arachidonic acid metabolism and linolenic acid metabolism pathways, and had certain correlation with microorganisms. Integrated transcriptome-metabolome-microbiome analysis showed that epithelial gene-GSTM3 expression was upregulated in the metabolism of xenobiotics by the cytochrome P450 pathway during the cold season, leading to the downregulation of some harmful metabolites; TLR5 gene expression was upregulated and CD14 gene expression was downregulated in the legionellosis pathway during the cold season. This study comprehensively described the interaction mechanism between the rumen host and microbes and their metabolites in grazing Tibetan sheep during the cold season. Rumen epithelial genes, microbiota and metabolites act together in some key pathways related to cold season adaptation.

Long noncoding RNA (lncRNA) plays an important regulatory role in mammalian adipogenesis and lipid metabolism. However, their function in the longissimus thoracis (LT) muscle of fatty acid metabolism of Tibetan sheep remains undefined. In this study, fatty acid and fat content in LT muscle of Tibetan sheep were determined, and RNA sequencing was performed to reveal the temporal regularity of lncRNA expression and the effect of lncRNA-miRNA-mRNA ceRNA regulatory network on lipid metabolism of LT muscle in Tibetan sheep at four growth stages (4-month-old, 4 m; 1.5-year-old, 1.5 y; 3.5-year-old, 3.5 y; 6-year-old, 6 y). The results indicated that the intramuscular fat (IMF) content was highest at 1.5 y. Moreover, the monounsaturated fatty acid (MUFA) content in 1.5 y of Tibetan sheep is significantly higher than those of the other groups (P < 0.05), and it was also rich in a variety of polyunsaturated fatty acids (PUFA). A total of 360 differentially expressed lncRNAs (DE lncRNAs) were identified from contiguous period transcriptome comparative groups of 4 m vs. 1.5 y, 1.5 y vs. 3.5 y, 3.5 y vs. 6 y, and 4 m vs. 6 y, respectively. Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis found that the target genes in lncRNA trans-mRNA were significantly related to the protein digestion, absorption, and fatty acid biosynthesis pathways (P < 0.05), which demonstrated that DE lncRNA trans-regulated the target genes, and further regulated the growth and development of the LT muscle and intramuscular fatty acid metabolism in Tibetan sheep. We further analyzed the role of the lncRNA-miRNA-mRNA regulatory network in the lipid metabolism of Tibetan sheep. Additionally, GPD2, LIPE (lipase E hormone-sensitive enzyme), TFDP2, CPT1A, ACACB, ADIPOQ, and other mRNA related to fatty acid and lipid metabolism and the corresponding lncRNA-miRNA regulatory pairs were identified. The enrichment analysis of mRNA in the regulatory network found that the AMPK signaling pathway was the most significantly enriched (P = 0.0000112361). Comprehensive transcriptome analysis found that the LIPE, ADIPOQ, ACACB, and CPT1A that were regulated by lncRNA might change the formation of energy metabolism in Tibetan sheep muscle through the AMPK signaling pathway, and oxidized muscle fibers are transformed into glycolytic muscle fibers, reduced IMF content, and the fatty acid profile also changed.

The Tianzhu white yaks (Bos grunniens) live on the Qinghai–Tibet Plateau. During winter, a lack of resources and low nutritional levels seriously affect their growth performance. In this study, we aimed to explore the effect of supplementation straw pellets on the growth performance and ruminal microbiota of yaks. Overall, at 6 (6M, n = 24), 18 (18M, n = 26), 30 (30M, n = 20), 42 (42M, n = 24), and 54 (54M, n = 22) month old Tianzhu white yaks were selected (total n = 116) and divided into the mixed straw + grazing (MSG), corn straw + grazing (CSG), and the grazing control (G) groups according to age and gender. Their growth performance was measured as per different dietary treatments. The rumen microbial community structure and levels of VFAs were analyzed from the 6M, 30M, and 54M male yaks from each group. The supplementary diets led to an increase in the ADG, which was the highest in the MSG group. The MSG group exhibited the highest level of acetate and total VFAs (TVFAs) among the three groups (p < 0.05). In addition, the 16S rRNA sequencing results proved that the microbial composition was dominated by the members of Firmicutes and Bacteroidetes. Christensenellaceae R-7 group was significantly abundant in the CSG and MSG groups compared to the G group (p < 0.05). Principal coordinate analysis (PCoA) revealed that the bacterial community structure of rumen in the MSG and CSG groups was considerably different from that in the G group; 6M samples exhibited different rumen microbial diversity compared with the other samples. Correlation analysis revealed that Christensenellaceae_R-7_group was positively correlated with the levels of acetate, TVFAs, and ADG. These results demonstrated that mixed straw pellets improved the growth performance of yaks, increased the abundance of Christensenellaceae R-7_group involved in cellulose degradation in the rumen, and produced large amounts of VFAs, which were absorbed by yaks, thus increasing their ADG. This study provides new insights into the effects of straw pellet supplementation on the changes in the rumen microbiota and growth performance of yaks.

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.

Milk fat is the foremost nutrient of milk and a vital indicator in evaluating milk quality. Accumulating evidence suggests that microRNAs (miRNAs) are involved in the synthesis of milk fat. The miR-200c is closely related to lipid metabolism, but little is known about its effect on the synthesis of milk fat in MECs of ewes. Herein, the effect of miR-200c on the proliferation of ovine mammary epithelial cells (MECs) and its target relationship with a predicted target gene were investigated. The regulatory effects of miR-200c on the expression of the target genes and the content of triglycerides in ovine MECs were further analyzed. The results revealed that the expression level of miR-200c was differentially expressed in both eight tissues selected during lactation and in mammary gland tissues at different physiological periods. Overexpression of miR-200c inhibited the viability and proliferation of ovine MECs, while inhibition of miR-200c increased cell viability and promoted the proliferation of ovine MECs. Target gene prediction results indicated that miR-200c would bind the 3′UTR region of pantothenate kinase 3 (PANK3). Overexpression of miR-200c reduced the luciferase activity of PANK3, while inhibition of miR-200c increased its luciferase activity. These findings illustrated that miR-200c could directly interact with the target site of the PANK3. It was further found that overexpression of miR-200c reduced the expression levels of PANK3 and, thus, accelerated the synthesis of triglycerides. In contrary, the inhibitor of miR-200c promoted the expression of PANK3 that, thus, inhibited the synthesis of triglycerides in ovine MECs. Together, these findings revealed that miR-200c promotes the triglycerides synthesis in ovine MECs via increasing the lipid synthesis related genes expression by targeting PANK3.

The rumen is an important hallmark organ of ruminants and plays an important role in the metabolism and immune barrier of Tibetan sheep on the Plateau. However, there are few studies on rumen development and metabolism regulation in Tibetan sheep at different ages. Here, we comprehensively analyzed the immune function, fermentation function, rumen epithelial micromorphology and transcriptome profile of Tibetan sheep at different ages. The results showed that the concentration of IgG decreased and the concentration of IgM increased with age (p < 0.05), and the highest concentration of IgA was observed at 1.5 and 3.5 years of age. In terms of rumen fermentation characteristics, VFAs of 4-month-old lambs were the highest, followed by VFAs and NH3-N of Tibetan sheep at 3.5 years of age. Hematoxylin-eosin staining and transmission electron microscopy section examination of rumen epithelial tissue showed that the rumen papilla width increased with age (p < 0.001), the thickness of the stratum corneum decreased, the cells in the stratum corneum showed accelerated migration and the thickness of the rumen muscle layer increased (p < 0.001). Desmosomal junctions between the layers of rumen epithelium increased at 1.5 and 3.5 years old, forming a compact barrier structure, and the basal layer had more mitochondria involved in the regulation of energy metabolism. RNA-seq analysis revealed that a total of 1006 differentially expressed genes (DEGs) were identified at four ages. The DEGs of Tibetan sheep aged 4 months and 6 years were mainly enriched in the oxidation–reduction process and ISG15-protein conjugation pathway. The 1.5 and 3.5-year-olds were mainly enriched in skeletal muscle thin filament assembly, mesenchyme migration and the tight junction pathway. WGCNA showed that DEGs related to rumen microbiota metabolite VFAs and epithelial morphology were enriched in “Metabolism of xenobiotics by cytochrome P450, PPAR signaling pathway, Butanoate metabolism pathways” and participated in the regulation of rumen epithelial immune and fermentation metabolism functions of Tibetan sheep at different ages. This study systematically revealed the regulatory mechanism of rumen epithelial development and metabolism in the plateau adaptation of Tibetan sheep, providing a new approach for the study of plateau adaptation.

The aim of this study was to investigate the impacts of inulin and Chinese gallotannin on the meat fatty acids and urinary metabolites in sheep. Twenty-four healthy (25.80 ± 3.85 kg) weaned Hu lambs of approximately 4.5 months old were equally divided into four groups: control group (basal diet), treatment group I (basal diet + 0.1% inulin), treatment group II (basal diet + 0.1% inulin + 2% Chinese gallotannin), and treatment group III (basal diet + 0.1% inulin + 2% Chinese gallotannin + 4% PEG). The contents of myristic acid (C14:0) and palmitic acid (C16:0) were found to be lower in treatment group II than in the control group (p < 0.05). Moreover, the palmitoleic acid (C16:1) content in treatment group II was notably higher than that in the control group (p < 0.05), while the elaidic acid (C18:1n9t) content in treatment group II was higher than that in other groups (p < 0.05). Besides, the linoleic acid (C18:2n6c) content was higher in the treatment II and control groups than in the treatment I and III groups. Furthermore, compared with the control group, both 4-pyridoxic acid and creatinine in treatment groups I and II were upregulated (p < 0.05), while other metabolites, such as nicotinuric acid, l-threonine, palmitic acid, and oleic acid, were drastically downregulated (p < 0.05). These differential metabolites were found to be mainly involved in nicotinate and nicotinamide metabolism (ko00760), vitamin B6 metabolism (ko00750), and the fatty acid biosynthesis pathway (ko00061). It is concluded that the combination of inulin and Chinese gallotannin in the diet could improve the energy and lipid metabolism of sheep, which may improve both mutton quality and production performance.