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.

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

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

Identification of quantitative trait loci and candidate genes for grain superoxide dismutase activity in wheat.

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.

Effects of the Interaction between Rumen Microbiota Density-VFAs-Hepatic Gluconeogenesis on the Adaptability of Tibetan Sheep to Plateau.

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.

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.

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

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

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.

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.

Functional differentiation of the ovine preadipocytes -insights from gene expression profiling.

The preadipocytes differentiation is a vital process of lipogenesis; exploring the molecular mechanisms of lipogenesis contributes to improve the meat quality and final commercial income. Lipogenesis has been widely reported in other livestock, but little is known about the gene expression profiles at different stages during preadipocytes differentiation in sheep. In this study, ovine preadipocytes were cultured in vitro and then induced to begin differentiation. Then, the gene expression profiles of preadipocytes collected on day 0 (D0), day 2 (D2), and day 8 (D8) of differentiation were analyzed by RNA-seq technology. According to the findings, 2254 differentially expressed genes (DEGs) were found in D2 vs D0; 1817 DEGs and 1902 DEGs were found in D8 vs D0 and D8 vs D2, respectively. The DEGs were found to be enriched in several biological processes, including focal adhesion, ECM-receptor interaction, PI3K-Akt signaling pathway, steroid biosynthesis, and MAPK signaling pathway, according to Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. The regulatory network of the DEGs related to ovine preadipocytes differentiation was systematically constructed, which showed that hub genes might modulate ovine preadipocytes differentiation. In summary, preadipocyte differentiation is regulated by several key genes, including ACACB, CXCL6, SREBF1, INSIG1, APOE, GJA1, CDH11, SYNE1, PCSK1, S100A4, FN1, PLIN2, CXCL6, FN1, PTX3, and FABP3. This study provides a deeper knowledge of the roles of genes in sheep lipogenesis by revealing global gene expression profiles during preadipocyte differentiation.

Integrated transcriptome analysis reveals roles of long non-coding RNAs (lncRNAs) in caprine skeletal muscle mass and meat quality.

Long non-coding RNAs (lncRNAs) play important roles in the growth and development of skeletal muscle. However, there is limited information on goats. In this study, expression profiles of lncRNAs in Longissimus dorsi muscle from Liaoning cashmere (LC) goats and Ziwuling black (ZB) goats with divergent meat yield and meat quality were compared using RNA-sequencing. Based on our previous microRNA (miRNA) and mRNA profiles obtained from the same tissues, the target genes and binding miRNAs of differentially expressed lncRNAs were obtained. Subsequently, lncRNA-mRNA interaction networks and a ceRNA network of lncRNA-miRNA-mRNA were constructed. A total of 136 differentially expressed lncRNAs were identified between the two breeds. Fifteen cis target genes and 143 trans target genes were found for differentially expressed lncRNAs, and they were enriched in muscle contraction, muscle system process, muscle cell differentiation, and p53 signaling pathway. A total of 69 lncRNA-trans target gene pairs were constructed, with close relationship with muscle development, intramuscular fat deposition, and meat tenderness. A total of 16 lncRNA-miRNA-mRNA ceRNA pairs were identified, of which some reportedly associated with skeletal muscle development and fat deposition were found. The study will provide an improved understanding of the roles of lncRNAs in caprine meat yield and meat quality.

Interactions of rumen microbiota and metabolites with meat quality-related genes to regulate meat quality and flavor of Tibetan sheep under nutrient stress in the cold season.

AIM: The meat of Tibetan sheep has a unique flavor, delicious taste, and superior nutritional value. However, the change of grass will lead to a change in meat quality. This study aimed to explore the potential regulatory mechanisms of microbial metabolites with respect to meat quality traits of Tibetan sheep under nutrient stress in the cold season. METHODS AND RESULTS: We determined and analyzed the longissimus dorsi quality, fatty acid composition, expression of genes, and rumen microbial metabolites of Tibetan sheep in cold and warm seasons. The shear force was decreased (P < .05), the meat color a*24 h value was increased (P < .05), and the contents of crude fat (EE) and protein (CP) were decreased in the cold season. Polyunsaturated fatty acids (PUFAs)-linoleic acid and docosahexaenoic acid increased significantly in the cold season (P < .05). The expressions of meat quality genes MC4R, CAPN1, H-FABP, and LPL were significantly higher in the warm season (P < .05), and the CAST gene was significantly expressed in the cold season (P < .01). The different microbial metabolites of Tibetan sheep in the cold and warm seasons were mainly involved in amino acid metabolism, lipid metabolism, and digestive system pathway, and there was some correlation between microbiota and meat quality traits. There are similarities between microbial metabolites enriched in the lipid metabolism pathway and muscle metabolites. CONCLUSION: Under nutritional stress in the cold season, the muscle tenderness of Tibetan sheep was improved, and the fat deposition capacity was weakened, but the levels of beneficial fatty acids were higher than those in the warm season, which was more conducive to healthy eating.

Astragalus additive in feed improved serum immune function, rumen fermentation and the microbiota structure of early-weaned lambs.

AIM: The purpose of this study was to determine the mechanism of Astragalus activity on the immune function, rumen microbiota structure, and rumen fermentation of early-weaned lambs. METHODS AND RESULTS: Thirty healthy early-weaned lambs with similar body weights (17.42 ± 2.02 kg) were selected for the feeding experiment. The control group (KB) was fed a basal diet, and the Astragalus group (HQ) was fed 0.3% Astragalus additive on the basis of a basic diet. The formal trial period was 60 days. The results showed that the concentrations of blood immunoglobulin A (IgA) and immunoglobulin M (IgM) in the HQ group were significantly higher than those in the KB group (P < 0.05). Compared with the KB group, the concentrations of acetic acid, butyric acid, and total volatile fatty acids (VFAs) in the HQ group were higher (P < 0.01). The expression levels of the rumen epithelial-related genes MCT1, MCT4, NHE2, and ZO1 in the Astragalus group were significantly higher than those in the KB group (P < 0.05). 16S rRNA analysis showed that at the phylum level, Bacteroidetes in the HQ group significantly increased (P < 0.01); at the genus level, Prevotella (P < 0.01) and Succiniclasticum (P < 0.01) in the HQ group were found at significantly higher abundances than those in the KB group, and the results of microbiota gene and function prediction showed that "energy metabolism," "glycan biosynthesis and metabolic" pathways were significantly enriched in the HQ group (P < 0.05). CONCLUSION: As a feed additive, Astragalus can improve the immunity of early-weaned lambs, the structure of the rumen microbiota of lambs, and the fermentation capacity of the rumen.

Genome-Wide Identification and Chilling Stress Analysis of the NF-Y Gene Family in Melon.

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.

Interaction between Rumen Epithelial miRNAs-Microbiota-Metabolites in Response to Cold-Season Nutritional Stress in Tibetan Sheep.

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.

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

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

Variation in Acetyl-CoA Carboxylase Beta Gene and Its Effect on Carcass and Meat Traits in Gannan Yaks.

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.

LncRNA and Protein Expression Profiles Reveal Heart Adaptation to High-Altitude Hypoxia in Tibetan Sheep.

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.

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.

Integration of ATAC-Seq and RNA-Seq Analysis to Identify Key Genes in the Longissimus Dorsi Muscle Development of the Tianzhu White Yak.

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.