The concentration of bovine milk fat changes regularly with lactation stages. In particular, milk fat percentage is higher in late lactation than mid lactation. Furthermore, milk fat composition is highly subject to a few genes. Thus, transcriptome sequencing was performed to explore the expression patterns of differentially-expressed genes (DEGs) in the parenchymal mammary gland of Holstein dairy cows between mid and late lactation. The 725 DEGs were screened (fold change > 2 and p-value < 0.05), and the peroxisome proliferator-activated receptor (PPAR) signaling pathway associated with lipid synthesis had a significant variation between the two periods (p-value < 0.05). The activation of the PPAR signal pathway may a key factor in the increasing of milk fat content in late lactation compared to mid lactation. Acyl-CoA synthetase long-chain family member 4 (ACSL4), a member of the PPAR signaling pathway, was upregulated in late lactation compared to mid lactation (p < 0.05). ACSL4 catalyzes the activation of long-chain fatty acids for cellular lipid synthesis. However, it remains uncertain that the molecular mechanism of milk fat synthesis is regulated by ACSL4 in dairy cows. Subsequently, the function verification of ACSL4 was performed in bovine mammary epithelial cells (BMECs). The upregulated expression of ACSL4 was accompanied by the increase of the concentration of intracellular triglycerides, whereas knockdown of ACSL4 decreased the concentration of intracellular triglycerides, which demonstrated that ACSL4 plays an important role in modulating milk fat synthesis. In conclusion, the results displayed that ACSL4 expression regulates triglyceride metabolism in ruminant mammary cells.

The beef aging process is essential for compliance with certain major requisites, such as sensory characteristics for cooking and meat processing. Meat quality analysis of Yunling cattle, a new hybrid beef cattle bred by Chinese researchers, during the aging process, represents a major research gap. To explore Yunling beef initially, indicators associated with meat quality during the aging process of Yunling, Simmental, and Wenshan cattle were measured. In addition, some important economic traits were detected in the three breeds, including growth performance and carcass characteristics. The results showed that the growth performance, carcass traits, pH, and water holding capacity of Yunling and Simmental cattle were basically the same and better, respectively, than those of Wenshan cattle. The proportions of individual fatty acids in Yunling beef were healthier than in the other two breeds. Aging time did not affect the fatty acid profiles of the beef (p > 0.05). The contents of certain fatty acids in the three beef types displayed some differences in terms of days of aging (p < 0.05). The tenderness and meat color were better in the Yunling beef as the aging time increased, indicating that Yunling beef aged for 7 days was more suitable for cooking, exhibiting better sensory characteristics. Thus, a 7-day short-term aging process is very effective in improving the quality of Yunling beef. Our study attempted to fill a gap in the Yunling beef quality analysis during aging, providing further evidence for Yunling beef improvement.

The existing research on dairy cow mammary gland genes is extensive, but there have been few reports about dynamic changes in dairy cow mammary gland genes as milk yield decrease. For the first time, transcriptome analysis based on short time-series expression miner (STEM) and histological observations were performed using the Holstein dairy cow mammary gland to explore gene expression patterns in this process of decrease (at peak, mid-, and late lactation). Histological observations suggested that the number of mammary acinous cells at peak/mid-lactation was significantly higher than that at mid-/late lactation, and the lipid droplets area secreted by dairy cows was almost unaltered across the three stages of lactation (p > 0.05). Totals of 882 and 1439 genes were differentially expressed at mid- and late lactation, respectively, compared to peak lactation. Function analysis showed that differentially expressed genes (DEGs) were mainly related to apoptosis and energy metabolism (fold change ≥ 2 or fold change ≤ 0.5, p-value ≤ 0.05). Transcriptome analysis based on STEM identified 16 profiles of differential gene expression patterns, including 5 significant profiles (false discovery rate, FDR ≤ 0.05). Function analysis revealed DEGs involved in milk fat synthesis were downregulated in Profile 0 and DEGs in Profile 12 associated with protein synthesis. These findings provide a foundation for future studies on the molecular mechanisms underlying mammary gland development in dairy cows.

Cardiovascular diseases (CVDs) are seriously threatening to human life and health. Polyunsaturated fatty acids (PUFAs) are known for their role in preventing CVDs. It is beneficial to population health to promote the content of PUFAs in bovine milk. In recent years, limited research based on molecular mechanisms has focused on this field. The biological roles of numerous microRNAs (miRNAs) remain unknown. In this study, a promising and negatively correlated pair of the miRNA (miRNA-193a-5p) and a fatty acid desaturase 1 (FADS1) gene are identified and screened to explore whether they are potential factors of PUFAs’ synthesis in bovine milk. The targeted relationship between miRNA-193a-5p and FADS1 in bovine mammary epithelial cells (BMECs) is demonstrated by dual luciferase reporter assays. qRT-PCR and western blot assays indicate that both the expression of mRNA and the protein FADS1 show a negative correlation with miRNA-193a-5p expression in BMECs. Also, miR-193a-5p expression is positively correlated with the expression of genes associated with milk fatty acid metabolism, including ELOVL fatty acid elongase 6 (ELOVL6) and diacylglycerol O-acyltransferase 2 (DGAT2). The expression of fatty acid desaturase 2 (FADS2) is negatively correlated with miR-193a-5p expression in BMECs. The contents of triglycerides (TAG), eicosapentaenoic acid (EPA), docosapentaenoic acid (DPA) and docosahexaenoic acid (DHA) have a significant positive correlation with the expression of FADS1 and a significant negative correlation with the expression of miR-193a-5p in BMECs. For the first time, this study confirms that miRNA-193a-5p regulates PUFAs metabolism in BMECs by targeting FADS1, indicating that miRNA-193a-5p and FADS1 are underlying factors that improve PUFAs content in bovine milk.

Weighted gene coexpression network analysis (WGCNA) is a novel approach that can quickly analyze the relationships between genes and traits. In this study, the milk yield, lactose, fat, and protein of Holstein dairy cows were detected in a lactation cycle. Meanwhile, a total of 18 gene expression profiles were detected using mammary glands from six lactation stages (day 7 to calving, −7 d; day 30 post-calving, 30 d; day 90 post-calving, 90 d; day 180 post-calving, 180 d; day 270 post-calving, 270 d; day 315 post-calving, 315 d). On the basis of the 18 profiles, WGCNA identified for the first time 10 significant modules that may be related to lactation stage, milk yield, and the main milk composition content. Genes in the 10 significant modules were examined with gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. The results revealed that the galactose metabolism pathway was a potential candidate for milk yield and milk lactose synthesis. In −7 d, ion transportation was more frequent and cell proliferation related terms became active. In late lactation, the suppressor of cytokine signaling 3 (SOCS3) might play a role in apoptosis. The sphingolipid signaling pathway was a potential candidate for milk fat synthesis. Dairy cows at 315 d were in a period of cell proliferation. Another notable phenomenon was that nonlactating dairy cows had a more regular circadian rhythm after a cycle of lactation. The results provide an important theoretical basis for the further molecular breeding of dairy cows.

Understanding how hybrids integrate lineage-specific regulatory variants at the haplotype level is crucial for elucidating the genetic basis of heterosis in livestock. In this study, we established three crossbred pig families derived from distant genetic lineages and systematically identified variants from different lineages, including single nucleotide polymorphisms (SNPs) and structural variations (SVs). At the phase level, we quantitatively analyzed gene expression, four histone modifications (H3K4me3, H3K27ac, H3K4me1, and H3K27me3), and the binding strength of transcription factor (CTCF) in backfat (BF) and longissimus dorsi (LD) muscle. By colocalization analysis of phased genetic variants with phased gene expression levels and with phased epigenetic modifications, we identified 18,670 expression quantitative trait loci (eQTL) (FDR < 0.05) and 8,652 epigenetic modification quantitative trait loci (epiQTL) (FDR < 0.05). The integration of eQTL and epiQTL allowed us to explore the potential regulatory mechanisms by which lineage-specific genetic variants simultaneously influence gene expression and epigenetic modifications. For example, we identified a Large White lineage-specific duplication (DUP) encompassing the KIT gene that was significantly associated with its promoter activity (FDR = 7.83 × 10⁻⁴) and expression levels (FDR = 9.03 × 10⁻⁴). Additionally, we found that a Duroc lineage-specific SNP located upstream of AMIGO2 was significantly associated with a Duroc-specific H3K27ac peak (FDR = 0.035) and also showed a significant association with AMIGO2 expression levels (FDR = 5.12 × 10⁻⁴). These findings underscore the importance of phased regulatory variants in shaping lineage-specific transcriptional programs and highlight how the haplotype-resolved integration of eQTL and epigenetic signals can reveal the mechanistic underpinnings of hybrid regulatory architecture. Our results offer insights for molecular marker development in precision pig breeding.