[Spatio-temporal Evolution and Prediction of Carbon Storage in Huaibei City Based on InVEST-PLUS Model].

This study aimed to investigate the impact of spatiotemporal changes in land use on ecosystem carbon storage. The study analyzed the spatiotemporal changes in carbon storage in the study area based on land use data from five periods (1985, 1995, 2005, 2015, and 2020) using the InVEST model. The PLUS model was used to predict land use changes in the study area under four different scenarios (natural development, farmland protection, ecological protection, and double protection of farmland and ecology) in 2035, and the ecosystem carbon storage under different scenarios was estimated. The results of the study indicated that the farmland in the area under investigation had been decreasing consistently from 1985 to 2020, with a more rapid rate of change observed between 2015 and 2020. During this period, the overall dynamic attitude towards land use reached 34.62 %. Additionally, the carbon storage in the area showed a decreasing trend over the years, with a decrease of 1.55×105 t from 1985 to 2020. Between 2005 and 2015, the carbon storage showed a decrease of 1.22×105 t, with an average annual decrease of 1.22×104 t. The areas with higher carbon storage were located in the eastern part of the study area, whereas areas with lower carbon storage were found in the central and northwestern parts. Although the proportion of carbon storage in farmland decreased from 66.89 % to 57.73 %, farmland remained the most important carbon pool in the study area. The conversion of other land use types to grassland and forestland was advantageous for increasing ecosystem carbon storage. Finally, the study projected that by 2035, the carbon storage in the natural development scenario, the farmland protection scenario, the ecological protection scenario, and the dual protection scenario would be 81.77×105, 82.45×105, 82.82×105, and 82.51×105 t, respectively.

Proteome Analysis Related to Unsaturated Fatty Acid Synthesis by Interfering with Bovine Adipocyte ACSL1 Gene.

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.

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.

Integrative ATAC-seq and RNA-seq Analysis of the Longissimus Dorsi Muscle of Gannan Yak and Jeryak.

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.

A conjoint analysis of renal structure and omics characteristics reveal new insight to yak high-altitude hypoxia adaptation.

BACKGROUND: Yaks have unique adaptive mechanisms to the hypoxic environment, in which the kidney plays an important role. The aim of this study was to explore the histological changes of yak kidney at different altitudes and the metabolites and genes associated with adaptation to the hypoxic environment. METHODS: We analyzed the tissue structure and transcriptomic metabolomic data of yak kidney tissue at two altitudes, 2600 and 4400 m. We compared and identified the morphological adaptations of the kidney and the metabolites and genes associated with hypoxia adaptation in yaks. Changes in renal morphological adaptations, differential metabolites and genes were compared and identified, combining the two in a joint analysis. RESULTS: High-altitude yak kidneys showed significant adaptive changes: increased mitochondria, increased glomerular thylakoid area, and decreased localized ribosomes. Transcriptomics and metabolomics identified 69 DAMs (Differential metabolites) and 594 DEGs (differential genes). Functional enrichment analysis showed that the DAMs were associated with protein digestion and absorption, ABC transporter, and MTOR signaling pathway; the DEGs were significantly enriched in Cholesterol metabolism and P53 signaling pathway. The joint analysis indicated that metabolites such as lysine and arginine, as well as key genes such as ABCB5 and COL1A2, were particularly affected under hypoxic conditions, whereas changes in mitochondria in the tissue structure may be related to the expression of MFN1 and OPA1, and changes in glomerular thylakoid membranes are related to VEGFA and TGFB3. CONCLUSION: The kidney regulates metabolites and gene expression related to hormone synthesis, protein metabolism, and angiogenesis by adjusting the mitochondrial and glomerular thylakoid membrane structure to support the survival of yaks in high-altitude environments.

Unlocking the Transcriptional Control of NCAPG in Bovine Myoblasts: CREB1 and MYOD1 as Key Players.

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.

Contrastive representation learning on dynamic networks.

Representation learning for dynamic networks is designed to learn the low-dimensional embeddings of nodes that can well preserve the snapshot structure, properties and temporal evolution of dynamic networks. However, current dynamic network representation learning methods tend to focus on estimating or generating observed snapshot structures, paying excessive attention to network details, and disregarding distinctions between snapshots with larger time intervals, resulting in less robustness for sparse or noisy networks. To alleviate these challenges, this paper proposes a contrastive mechanism for temporal representation learning on dynamic networks, inspired by the success of contrastive learning in visual and static network representation learning. This paper proposes a novel Dynamic Network Contrastive representation Learning (DNCL) model. Specifically, contrast objective functions are constructed using intra-snapshot and inter-snapshot contrasts to capture the network topology, node feature information, and network evolution information, respectively. Rather than estimating or generating ground-truth network features, the proposed approach maximizes mutual information between nodes from different time steps and views generated. The experimental results of link prediction, node classification, and clustering on several real-world and synthetic networks demonstrate the superiority of DNCL over state-of-the-art methods, indicating the effectiveness of the proposed approach for dynamic network representation learning.

Mutations in the FOXO3 Gene and Their Effects on Meat Traits in Gannan Yaks.

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.

Molecular Genetic Characteristics of the Hoxc13 Gene and Association Analysis of Wool 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.

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.

Oxygen vacancy-rich Fe2(MoO4)3 combined with MWCNTs for electrochemical sensors of fentanyl and its analogs.

Hundreds of thousands of people dying from the abuse of fentanyl and its analogs. Hence, the development of an efficient and highly accurate detection method is extremely relevant and challenging. Therefore, we proposed the introduction of oxygen defects into Fe2(MoO4)3 nanoparticles for improving the catalyst performance and combining it with multi-walled carbon nanotubes (MWCNTs) for electrochemical detection of fentanyl and its analogs. Oxygen vacancy-rich Fe2(MoO4)3 (called r-Fe2(MoO4)3) nanoparticles were successfully synthesized and characterized in detail by scanning electron microscopy (SEM), transmission electron microscopy (TEM), energy-dispersive spectrometry (EDS), X-ray diffraction (XRD), Fourier transform infrared (FT-IR), Raman spectra, BET, X-ray photoelectron spectroscopy (XPS), and electron paramagnetic resonance (EPR) and investigated by comparison with oxygen vacancy-poor Fe2(MoO4)3 (called p-Fe2(MoO4)3). The obtained oxygen vacancy-rich Fe2(MoO4)3 was ultrasonically composited with MWCNTs for modification of glassy carbon electrodes (GCEs) used for the electrochemical detection of fentanyl and its analogs. The modified MWCNT-GCE showed ultrasensitivity to fentanyl, sufentanil, alfentanil, and acetylfentanyl with limits of detection (LOD) of 0.006 µmol·L-1, 0.008 µmol·L-1, 0.018 µmol·L-1, and 0.024 µmol·L-1, respectively, and could distinguish among the four drugs based on their peak voltages. Besides, the obtained r-Fe2(MoO4)3/MWCNT composite also exhibited high repeatability, selectivity, and stability. It showed satisfactory detection performance on real samples, with recoveries of 70.53 ~ 94.85% and 50.98 ~ 82.54% in serum and urine for the four drugs in a concentration range 0.2 ~ 1 µM, respectively. The experimental results confirm that the introduction of oxygen vacancies effectively improves the sensitivity of fentanyl electrochemical detection, and this work provides some inspiration for the development of catalytic materials for electrochemical sensors with higher sensitivity.

Microgel-based carriers enhance skeletal stem cell reprogramming towards immunomodulatory phenotype in osteoarthritic therapy.

Skeletal stem cells (SSC) have gained attentions as candidates for the treatment of osteoarthritis due to their osteochondrogenic capacity. However, the immunomodulatory properties of SSC, especially under delivery operations, have been largely ignored. In the study, we found that Pdpn+ and Grem1+ SSC subpopulations owned immunoregulatory potential, and the single-cell RNA sequencing (scRNA-seq) data suggested that the mechanical activation of microgel carriers on SSC induced the generation of Pdpn+Grem1+Ptgs2+ SSC subpopulation, which was potent at suppressing macrophage inflammation. The microgel carriers promoted the YAP nuclear translocation, and the activated YAP protein was necessary for the increased expression of Ptgs2 and PGE2 in microgels-delivered SSC, which further suppressed the expression of TNF-ɑ, IL-1ß and promoted the expression of IL-10 in macrophages. SSC delivered with microgels yielded better preventive effects on articular lesions and macrophage activation in osteoarthritic rats than SSC without microgels. Chemically blocking the YAP and Ptgs2 in microgels-delivered SSC partially abolished the enhanced protection on articular tissues and suppression on osteoarthritic macrophages. Moreover, microgel carriers significantly prolonged SSC retention time in vivo without increasing SSC implanting into osteoarthritic joints. Together, our study demonstrated that microgel carriers enhanced SSC reprogramming towards immunomodulatory phenotype to regulate macrophage phenotype transformation for effectively osteoarthritic therapy by promoting YAP protein translocation into nucleus. The study not only complement and perfect the immunological mechanisms of SSC-based therapy at the single-cell level, but also provide new insight for microgel carriers in stem cell-based therapy.

Hybrid-FHR: a multi-modal AI approach for automated fetal acidosis diagnosis.

BACKGROUND: In clinical medicine, fetal heart rate (FHR) monitoring using cardiotocography (CTG) is one of the most commonly used methods for assessing fetal acidosis. However, as the visual interpretation of CTG depends on the subjective judgment of the clinician, this has led to high inter-observer and intra-observer variability, making it necessary to introduce automated diagnostic techniques. METHODS: In this study, we propose a computer-aided diagnostic algorithm (Hybrid-FHR) for fetal acidosis to assist physicians in making objective decisions and taking timely interventions. Hybrid-FHR uses multi-modal features, including one-dimensional FHR signals and three types of expert features designed based on prior knowledge (morphological time domain, frequency domain, and nonlinear). To extract the spatiotemporal feature representation of one-dimensional FHR signals, we designed a multi-scale squeeze and excitation temporal convolutional network (SE-TCN) backbone model based on dilated causal convolution, which can effectively capture the long-term dependence of FHR signals by expanding the receptive field of each layer's convolution kernel while maintaining a relatively small parameter size. In addition, we proposed a cross-modal feature fusion (CMFF) method that uses multi-head attention mechanisms to explore the relationships between different modalities, obtaining more informative feature representations and improving diagnostic accuracy. RESULTS: Our ablation experiments show that the Hybrid-FHR outperforms traditional previous methods, with average accuracy, specificity, sensitivity, precision, and F1 score of 96.8, 97.5, 96, 97.5, and 96.7%, respectively. CONCLUSIONS: Our algorithm enables automated CTG analysis, assisting healthcare professionals in the early identification of fetal acidosis and the prompt implementation of interventions.

Spatiotemporal Expression Characterization of KRTAP6 Family Genes and Its Effect on Wool Traits.

Keratin-related proteins (KAPs) are structural components of wool fibers and are thought to play a key role in regulating the physical and mechanical properties of fibers. Among all KAP genes (KRTAPs), KRTAP6 gene family (KRTAP6-1, KRTAP6-2, KRTAP6-3, KRTAP6-4, and KRTAP6-5) is a very important member with high polymorphism and notable association with some wool traits. In this study, we used real-time fluorescence quantitative PCR (RT-qPCR) and in situ hybridization to investigate spatiotemporal expression of KRTAP6s. The results revealed that KRTAP6 family genes were significantly expressed during anagen compared to other stages (p < 0.05). And it was found the five genes were expressed predominantly in the dermal papillae, inner and outer root sheaths, and showed a distinct spatiotemporal expression pattern. Also, it was found that KRTAP6-1 and KRTAP6-5 mRNA expression was negatively correlated with wool mean fiber diameter (MFD) and mean staple strength (MSS) (p < 0.05). In summary, the KRTAP6 family genes share a similar spatiotemporal expression pattern. And KRTAP6-1 and KRTAP6-5 may regulate the MFD and MSS of Gansu Alpine fine-wool sheep wool by changing the expression.

TNFAIP3 Derived from Skeletal Stem Cells Alleviated Rat Osteoarthritis by Inhibiting the Necroptosis of Subchondral Osteoblasts.

Recent investigations have shown that the necroptosis of tissue cells in joints is important in the development of osteoarthritis (OA). This study aimed to investigate the potential effects of exogenous skeletal stem cells (SSCs) on the necroptosis of subchondral osteoblasts in OA. Human SSCs and subchondral osteoblasts isolated from human tibia plateaus were used for Western blotting, real-time PCR, RNA sequencing, gene editing, and necroptosis detection assays. In addition, the rat anterior cruciate ligament transection OA model was used to evaluate the effects of SSCs on osteoblast necroptosis in vivo. The micro-CT and pathological data showed that intra-articular injections of SSCs significantly improved the microarchitecture of subchondral trabecular bones in OA rats. Additionally, SSCs inhibited the necroptosis of subchondral osteoblasts in OA rats and necroptotic cell models. The results of bulk RNA sequencing of SSCs stimulated or not by tumor necrosis factor α suggested a correlation of SSCs-derived tumor necrosis factor α-induced protein 3 (TNFAIP3) and cell necroptosis. Furthermore, TNFAIP3-derived from SSCs contributed to the inhibition of the subchondral osteoblast necroptosis in vivo and in vitro. Moreover, the intra-articular injections of TNFAIP3-overexpressing SSCs further improved the subchondral trabecular bone remodeling of OA rats. Thus, we report that TNFAIP3 from SSCs contributed to the suppression of the subchondral osteoblast necroptosis, which suggests that necroptotic subchondral osteoblasts in joints may be possible targets to treat OA by stem cell therapy.

Transcriptome Analysis of mRNA and lncRNA Related to Muscle Growth and Development in Gannan Yak and Jeryak.

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.

Variation in the HSL Gene and Its Association with Carcass and Meat Quality Traits in Yak.

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.

Comprehensive Analysis of miRNA and mRNA Expression Profiles during Muscle Development of the Longissimus Dorsi Muscle in Gannan Yaks and Jeryaks.

A hybrid offspring of Gannan yak and Jersey cattle, the Jeryak exhibits apparent hybrid advantages over the Gannan yak in terms of production performance and other factors. The small non-coding RNAs known as miRNAs post-transcriptionally exert a significant regulatory influence on gene expression. However, the regulatory mechanism of miRNA associated with muscle development in Jeryak remains elusive. To elucidate the regulatory role of miRNAs in orchestrating skeletal muscle development in Jeryak, we selected longissimus dorsi muscle tissues from Gannan yak and Jeryak for transcriptome sequencing analysis. A total of 230 (DE) miRNAs were identified in the longissimus dorsi muscle of Gannan yak and Jeryak. The functional enrichment analysis revealed a significant enrichment of target genes from differentially expressed (DE)miRNAs in signaling pathways associated with muscle growth, such as the Ras signaling pathway and the MAPK signaling pathway. The network of interactions between miRNA and mRNA suggest that some (DE)miRNAs, including miR-2478-z, miR-339-x, novel-m0036-3p, and novel-m0037-3p, played a pivotal role in facilitating muscle development. These findings help us to deepen our understanding of the hybrid dominance of Jeryaks and provide a theoretical basis for further research on the regulatory mechanisms of miRNAs associated with Jeryak muscle growth and development.

Long-stranded non-coding RNAs temporal-specific expression profiles reveal longissimus dorsi muscle development and intramuscular fat deposition in Tianzhu white yak.

The process of muscle development and intramuscular fat (IMF) deposition is quite complex and controlled by both mRNAs and ncRNAs. Long-stranded non-coding RNAs (LncRNAs) are involved in various biological processes in mammals while also playing a critical role in muscle development and fat deposition. In the present study, RNA-Seq was used to comprehensively study the expression of lncRNAs and mRNAs during muscle development and intramuscular fat deposition in postnatal Tianzhu white yaks at three stages, including 6 mo of age (calve, n = 6), 30 mo of age (young cattle, n = 6) and 54 mo of age (adult cattle, n = 6). The results indicated that a total of 2,101 lncRNAs and 20,855 mRNAs were screened across the three stages, of which the numbers of differential expression (DE) lncRNAs and DE mRNAs were 289 and 1,339, respectively, and DE lncRNAs were divided into eight different expression patterns based on expression trends. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that some DE mRNAs overlapped with target genes of lncRNAs, such as NEDD4L, SCN3B, AGT, HDAC4, DES, MYH14, KLF15 (muscle development), ACACB, PCK2, LIPE, PIK3R1, PNPLA2, and MGLL (intramuscular fat deposition). These DE mRNAs were significantly enriched in critical muscle development and IMF deposition-related pathways and GO terms, such as AMPK signaling pathway, PI3K-Akt signaling pathway, PPAR signaling pathway, etc. In addition, lncRNA-mRNA co-expression network analysis revealed that six lncRNAs (MSTRG.20152.2, MSTRG.20152.3, XR_001351700.1, MSTRG.8190.1, MSTRG.4827.1, and MSTRG.11486.1) may play a major role in Tianzhu white yak muscle development and lipidosis deposition. Therefore, this study enriches the database of yak lncRNAs and could help to further explore the functions and roles of lncRNAs in different stages of muscle development and intramuscular fat deposition in the Tianzhu white yak.

In-depth analysis of the molecular regulatory mechanisms involved in the development of yak skeletal muscle can aid in the identification of crucial functional genes that may enhance the quality of yak meat in the future. Our research is the first to investigate the expression of long-stranded non-coding RNAs (lncRNAs) during muscle development and intramuscular fat deposition in Tianzhu white yaks at different ages (6, 30, and 54 mo). We have discovered numerous novel lncRNAs and their associated target genes that are linked to muscle development and fat deposition, and have constructed relevant interaction network diagrams. This study is expected to promote a genetics-based approach to further enhance the Tianzhu white yak and other yak breeds.