Identification of Candidate Genes and Functional Pathways Associated with Body Size Traits in Chinese Holstein Cattle Based on GWAS Analysis.

Body size is one of the most economically important traits of dairy cattle, as it is significantly associated with cow longevity, production, health, fertility, and environmental adaptation. The identification and application of genetic variants using a novel genetic approach, such as genome-wide association studies (GWASs), may give more insights into the genetic architecture of complex traits. The identification of genes, single nucleotide polymorphisms (SNPs), and pathways associated with the body size traits may offer a contribution to genomic selection and long-term planning for selection in dairy cows. In this study, we performed GWAS analysis to identify the genetic markers and genes associated with four body size traits (body height, body depth, chest width, and angularity) in 1000 Chinese Holstein cows. We performed SNPs genotyping in 1000 individuals, based on the GeneSeek Genomic Profiler Bovine 100 K. In total, we identified 11 significant SNPs in association with body size traits at the threshold of Bonferroni correction (5.90 × 10-7) using the fixed and random model circulating probability unification (FarmCPU) model. Several genes within 200 kb distances (upstream or downstream) of the significant SNPs were identified as candidate genes, including MYH15, KHDRBS3, AIP, DCC, SQOR, and UBAP1L. Moreover, genes within 200 kb of the identified SNPs were significantly enriched (p ≤ 0.05) in 25 Gene Ontology terms and five Kyoto Encyclopedia of Genes and Genomes pathways. We anticipate that these results provide a foundation for understanding the genetic architecture of body size traits. They will also contribute to breeding programs and genomic selection work on Chinese Holstein cattle.

Analysis of Non-Genetic Factors Affecting Wood's Model of Daily Milk Fat Percentage of Holstein Cattle.

This research paper aimed to explore the characteristics of Holstein cattle's milk fat percentage lactation curve and its influencing factors. The Wood model was used for fitting the lactation curve of 398,449 DHI test-day milk fat percentage records of Holstein cows from 2018 to 2020 in 12 dairy farms in Jiangsu province, and the influencing factors­including farm size, parity, calving season, calving interval, and 305-days milk production­on the parameters of the lactation curve were analyzed. The results showed that the non-genetic factors such as dairy farm size, calving season, parity, calving interval, and 305-days milk yield have a significant impact on milk fat percentage (p < 0.01); the average R2 of the daily milk fat percentage curve was 0.9699; the lowest milk fat percentage was 3.54%; the time to reach the lowest milk fat percentage was 126 days; and the persistence of milk fat percentage was 3.59%. All of these factors explored in this study fit at different levels above 0.96. The Wood model performed well in the fitting and analysis of the milk fat percentage curve of Holstein cattle in Jiangsu Province. This study provides a reference for improving the milk fat percentage of Holstein cattle.

Effects of Seasonal Heat Stress during Late Gestation on Growth Performance, Metabolic and Immuno-Endocrine Parameters of Calves.

Heat stress during late gestation could affect subsequent lactation performance, resulting in damage to the immune function, health, and growth performance of calves. This study aimed to compare the effects of 33 days of summer stress (Summer group, 70.15 < THI < 74.28) with 33 days of winter during late gestation (Winter group, 57.55 < THI < 67.25) on the growth, hormones, oxidative stress, and immune function of calves. Calves (Summer, n = 28; Winter, n = 23) were separated from cows immediately after birth and fed with 2 L colostrum within 2 h and 8−10 h after birth, respectively, and weaned at 60 days of age. Bodyweight (BW) was measured at birth and weaning. Withers height (WH), body length, and chest girth were measured at birth, 30 days, and 60 days of age. The health of calves ranging in age from 1 to 7 days was recorded. Plasma interferon-γ (IFN-γ), superoxide dismutase (SOD), adrenocorticotropin (ACTH), gonadotropin-releasing hormone (GnRH), IgG, cortisol, heat shock protein 70 (Hsp70), growth hormone (GH), insulin, lipid peroxide (LPO), and tumor necrosis factor-α (TNF-α) levels were measured in calves at 0 (before colostrum feeding), 3, 7, 14, 28, and 56 days of age. The pregnancy period of the Summer group was shortened by 1.44 days. The Winter and Summer groups had the same birth weight. One week after birth, the incidence of diarrhea was 57.14% and 21.74% in Summer and Winter groups, respectively. Compared with the Winter group, TNF-α in the Summer group increased significantly before colostrum feeding. ACTH and LPO decreased significantly at 3 days of age, ACTH and TNF-α decreased significantly at 7 days of age, Hsp70 increased significantly, ACTH was significantly reduced at 14 days of age, and Hsp70 increased dramatically at 7 days of age. SOD and TNF-α increased statistically at 28 days of age, LPO decreased significantly, and IFN-γ decreased significantly at 56 days of age, while IgG and GH increased significantly. We conclude that maternal heat stress during late gestation can damage the oxidative stress and immune plasma indexes of offspring before weaning.

Identification and Characterization of Circular RNAs in Mammary Tissue from Holstein Cows at Early Lactation and Non-Lactation.

In this study, circular RNAs (circRNAs) from Holstein cow mammary tissues were identified and compared between early lactation and non-lactation. After analysis, 10,684 circRNAs were identified, ranging from 48 to 99,406 bp, and the average size was 882 bp. The circRNAs were mainly distributed on chromosomes 1 to 11, and 89.89% of the circRNAs belonged to sense-overlapping circRNA. The exons contained with circRNAs ranged from 1 to 47 and were concentrated from 1 to 5. Compared with the non-lactating cows, 87 circRNAs were significantly differentially expressed in the peak lactation cows. There were 68 upregulated circRNAs and 19 downregulated circRNAs. Enrichment analysis of circRNAs showed that GO analysis mainly focused on immune response, triglyceride transport, T cell receptor signaling pathway, etc. Pathway analysis mainly focused on cytokine-cytokine receptor interaction, T helper 17 cell differentiation, fatty acid biosynthesis, the JAK-STAT signaling pathway, etc. Specific primers were designed for two proximal ends of the circRNA junction sites to allow for PCR validation of four randomly selected circRNAs and carry out circRNA-miRNA interaction research. This study revealed the expression profile and characteristics of circRNAs in mammary tissue from Holstein cows at early lactation and non-lactation, thus providing rich information for the study of circRNA functions and mechanisms, as well as potential candidate miRNA genes for studying lactation in Holstein cows.

Influence of Heat Stress on Intestinal Epithelial Barrier Function, Tight Junction Protein, and Immune and Reproductive Physiology.

The potential threat of global warming in the 21st century is on the ecosystem through many aspects, including the negative impact of rising global temperature on the health of humans and animals, especially domestic animals. The damage caused by heat stress to animals has been more and more significant as the worldwide climate continues to rise, along with the breeding industry's expanding scale and stocking density, and it has become the most important stress-causing factor in southern China. In this review, we described the effects of heat stress on animal immune organs and immune system. The much-debated topic is how hyperthermia affects the tight junction barrier. Heat stress also induces inflammation in the body of animals causing low body weight and loss of appetite. This review also discussed that heat stress leads to hepatic disorder, and it also damages the intestine. The small intestine experiences ischemia, and the permeability of the intestine increases. Furthermore, the oxidative stress and mitogen-activated protein kinase (MAPK) pathways have a significant role in stress-induced cellular and organ injury. The study has shown that MAPK activity in the small intestine was increased by heat stress. Heat stress caused extreme small intestine damage, enhanced oxidative stress, and activated MAPK signaling pathways.

MicroRNA-193a-5p Regulates the Synthesis of Polyunsaturated Fatty Acids by Targeting Fatty Acid Desaturase 1 (FADS1) in Bovine Mammary Epithelial Cells.

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.

Metformin Inhibits Lipoteichoic Acid-Induced Oxidative Stress and Inflammation Through AMPK/NRF2/NF-κB Signaling Pathway in Bovine Mammary Epithelial Cells.

The objective of this research was to explore the effect of metformin on the lipoteichoic acid (LTA)-induced mastitis model using isolated primary bovine mammary epithelial cells (PBMECs). The PBMECs were exposed to either 3 mM metformin for 12 h as a metformin group (MET) or 100 µg/mL LTA for 6 h as LTA group (LTA). Cells pretreated with 3 mM metformin for 12 h followed by washing and 100 µg/mL LTA exposure for 6 h served as the MET + LTA group. Phosphate-buffered saline was added to cells as the control group. PBMECs pretreated with different metformin doses were analyzed by a flow cytometry (annexin V-fluorescein isothiocyanate assay) to detect the cell apoptotic rate. We performed quantitative reverse transcriptase-polymerase chain reaction and Western blot analysis to evaluate the inflammatory and oxidative responses to metformin and LTA by measuring cellular cytotoxicity, mRNA expression, and protein expression. Immunofluorescence was used to evaluate nuclear localization. The results showed that the gene expression of COX2, IL-1ß, and IL-6 significantly increased in the cells challenged with LTA doses compared to control cells. In inflammatory PBMECs, metformin attenuated LTA-induced expression of inflammatory genes nuclear factor κB (NF-κB) p65, tumor necrosis factor α, cyclooxygenase 2, and interleukin 1ß, as well as the nuclear localization and phosphorylation of NF-κBp65 protein, but increased the transcription of nuclear factor erythroid 2-related factor 2 (Nrf2) and Nrf2-targeted antioxidative genes heme oxygenase-1 (HO-1) and Gpx1, as well as the nuclear localization of HO-1 protein. Importantly, metformin-induced activation of Nrf2 is AMP-activated protein kinase (AMPK)-dependent; as metformin-pretreated PBMECs activated AMPK signaling via the upregulation of phosphorylated AMPK levels, cell pretreatment with metformin also reversed the translocation of Nrf2 that was LTA inhibited. This convergence between AMPK and Nrf2 pathways is essential for the anti-inflammatory effect of metformin in LTA-stimulated PBMECs. Altogether, our results indicate that metformin exerts anti-inflammation and oxidative stress through regulation of AMPK/Nrf2/NF-κB signaling pathway, which highlights the role of AMPK as a potential therapeutic strategy for treatment of bovine mastitis.

Genome-Wide Association Study Candidate Genes on Mammary System-Related Teat-Shape Conformation Traits in Chinese Holstein Cattle.

In the dairy industry, mammary system traits are economically important for dairy animals, and it is important to explain their fundamental genetic architecture in Holstein cattle. Good and stable mammary system-related teat traits are essential for producer profitability in animal fitness and in the safety of dairy production. In this study, we conducted a genome-wide association study on three traits-anterior teat position (ATP), posterior teat position (PTP), and front teat length (FTL)-in which the FarmCPU method was used for association analyses. Phenotypic data were collected from 1000 Chinese Holstein cattle, and the GeneSeek Genomic Profiler Bovine 100K single-nucleotide polymorphisms (SNP) chip was used for cattle genotyping data. After the quality control process, 984 individual cattle and 84,406 SNPs remained for GWAS work analysis. Nine SNPs were detected significantly associated with mammary-system-related teat traits after a Bonferroni correction (p < 5.92 × 10-7), and genes within a region of 200 kb upstream or downstream of these SNPs were performed bioinformatics analysis. A total of 36 gene ontology (GO) terms and 3 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were significantly enriched (p < 0.05), and these terms and pathways are mainly related to metabolic processes, immune response, and cellular and amino acid catabolic processes. Eleven genes including MMS22L, E2F8, CSRP3, CDH11, PEX26, HAL, TAMM41, HIVEP3, SBF2, MYO16 and STXBP6 were selected as candidate genes that might play roles in the teat traits of cows. These results identify SNPs and candidate genes that give helpful biological information for the genetic architecture of these teat traits, thus contributing to the dairy production, health, and genetic selection of Chinese Holstein cattle.

Genetic Parameter Estimation and Genome-Wide Association Study-Based Loci Identification of Milk-Related Traits in Chinese Holstein.

Accurately estimating the genetic parameters and revealing more genetic variants underlying milk production and quality are conducive to the genetic improvement of dairy cows. In this study, we estimate the genetic parameters of five milk-related traits of cows-namely, milk yield (MY), milk fat percentage (MFP), milk fat yield (MFY), milk protein percentage (MPP), and milk protein yield (MPY)-based on a random regression test-day model. A total of 95,375 test-day records of 9,834 cows in the lower reaches of the Yangtze River were used for the estimation. In addition, genome-wide association studies (GWASs) for these traits were conducted, based on adjusted phenotypes. The heritability, as well as the standard errors, of MY, MFP, MFY, MPP, and MPY during lactation ranged from 0.22 ± 0.02 to 0.31 ± 0.04, 0.06 ± 0.02 to 0.15 ± 0.03, 0.09 ± 0.02 to 0.28 ± 0.04, 0.07 ± 0.01 to 0.16 ± 0.03, and 0.14 ± 0.02 to 0.27 ± 0.03, respectively, and the genetic correlations between different days in milk (DIM) within lactations decreased as the time interval increased. Two, six, four, six, and three single nucleotide polymorphisms (SNPs) were detected, which explained 5.44, 12.39, 8.89, 10.65, and 7.09% of the phenotypic variation in MY, MFP, MFY, MPP, and MPY, respectively. Ten Kyoto Encyclopedia of Genes and Genomes pathways and 25 Gene Ontology terms were enriched by analyzing the nearest genes and genes within 200 kb of the detected SNPs. Moreover, 17 genes in the enrichment results that may play roles in milk production and quality were selected as candidates, including CAMK2G, WNT3A, WNT9A, PLCB4, SMAD9, PLA2G4A, ARF1, OPLAH, MGST1, CLIP1, DGAT1, PRMT6, VPS28, HSF1, MAF1, TMEM98, and F7. We hope that this study will provide useful information for in-depth understanding of the genetic architecture of milk production and quality traits, as well as contribute to the genomic selection work of dairy cows in the lower reaches of the Yangtze River.

Short Time-Series Expression Transcriptome Data Reveal the Gene Expression Patterns of Dairy Cow Mammary Gland as Milk Yield Decreased Process.

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.

Genome-Wide Association Study Identifies Candidate Genes Associated with Feet and Leg Conformation Traits in Chinese Holstein Cattle.

Feet and leg conformation traits are considered one of the most important economical traits in dairy cattle and have a great impact on the profitability of milk production. Therefore, identifying the single nucleotide polymorphisms (SNPs), genes and pathways analysis associated with these traits might contribute to the genomic selection and long-term plan selection for dairy cattle. We conducted genome-wide association studies (GWASs) using the fixed and random model circulating probability unification (FarmCPU) method to identify SNPs associated with bone quality, heel depth, rear leg side view and rear leg rear view of Chinese Holstein cows. Phenotypic measurements were collected from 1000 individuals of Chinese Holstein cattle and the GeneSeek Genomic Profiler Bovine 100 K SNP chip was utilized for individual genotyping. After quality control, 984 individual cows and 84,906 SNPs remained for GWAS work; as a result, we identified 20 significant SNPs after Bonferroni correction. Several candidate genes were identified within distances of 200 kb upstream or downstream to the significant SNPs, including ADIPOR2, INPP4A, DNMT3A, ALDH1A2, PCDH7, XKR4 and CADPS. Further bioinformatics analyses showed 34 gene ontology terms and two signaling pathways were significantly enriched (p ≤ 0.05). Many terms and pathways are related to biological quality, metabolism and development processes; these identified SNPs and genes could provide useful information about the genetic architecture of feet and leg traits, thus improving the longevity and productivity of Chinese Holstein dairy cattle.

Metformin acts to suppress ß-hydroxybutyric acid-mediated inflammatory responses through activation of AMPK signaling in bovine hepatocytes.

The occurrence of bovine ketosis involves the accumulation of ß-hydroxybutyric acid (BHBA), which contributes to the initiation and acceleration of hepatic metabolic stress and inflammation. Metformin has other beneficial effects apart from its medical intervention for diabetes, such as prevention of laminitis and hyper-triglyceridemic. AMPK maintains energy homeostasis and is the intracellular target of metformin action. This study aims to uncover the role of metformin in modulating BHBA-induced inflammatory responses through the activation of AMPK signaling. The hepatocytes were isolated from the liver tissue of mid-lactation multiparous Holstein cows (~160 d postpartum). Treatments were conducted as follows: treated with PBS for 18 h (control); pretreated with PBS for 12 h followed by treatment of 1.2 mM BHBA for 6 h (BHBA); pretreated with 1.5 mM or 3 mM metformin for 12 h followed by the BHBA treatment (1.2 mM) for 6 h (M(1.5)+B; M(3)+B). The inhibitor of AMPK, Compound C, at a concentration of 10 µM, was applied to substantiate the AMPK-dependent responses. RT-qPCR were applied for the mRNA expression while Western-blots and immunofluorescence were conducted for the target proteins expression. Among dose-dependent assays for BHBA, the concentration of BHBA at 1.2 mM activated NF-κB signaling by upregulating the expression of phosphorylated NF-κB and pro-inflammatory cytokines compared with the control cells (P < 0.05). Along with the upregulation of phosphorylated AMPKα and ACCα, metformin at 1.5 and 3 mM inactivated NF-κB signaling components (p65 and IκBα) and the inflammatory genes (TNFA, IL6, IL1B and COX-2) which were activated by BHBA. Additionally, BHBA inhibited cells staining intensity in EdU assay were increased by pretreatment with metformin. The activation of AMPK resulted in the increased gene and protein expression of SIRT1, along with the deacetylation of H3K9 and H3K14. However, the AMPK inhibitor compound C blocked this effect. Compared with BHBA treated cells, the protein expression of COX-2 and IL-1ß were decreased by the pretreatment with metformin, and the inhibitory effect of metformin was released by compound C. The bound of NF-κB onto IL1B promoter displayed higher in BHBA group and this was suppressed by pretreatment with metformin (P < 0.05). Altogether, metformin attenuates the BHBA-induced inflammation through the inactivation of NF-κB as a target for AMPK/SIRT1 signaling in bovine hepatocytes.

Lactation Associated Genes Revealed in Holstein Dairy Cows by Weighted Gene Co-Expression Network Analysis (WGCNA).

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.

Effect of Aging Time on Meat Quality of Longissimus Dorsi from Yunling Cattle: A New Hybrid Beef Cattle.

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.

Comparative Transcriptomic Analysis of the Pituitary Gland between Cattle Breeds Differing in Growth: Yunling Cattle and Leiqiong Cattle.

The hypothalamic-pituitary-thyroid (HPT) axis hormones regulate the growth and development of ruminants, and the pituitary gland plays a decisive role in this process. In order to identify pivotal genes in the pituitary gland that could affect the growth of cattle by regulating the secretion of hormones, we detected the content of six HPT hormones related to growth in the plasma of two cattle breeds (Yunling and Leiqiong cattle, both also known as the zebu cattle) with great differences in growth and compared the transcriptome data of their pituitary glands. Our study found that the contents of GH, IGF, TSH, thyroxine, triiodothyronine, and insulin were significantly different between the two breeds, which was the main cause of the difference in growth; 175 genes were identified as differentially expressed genes (DEGs). Functional association analyses revealed that DEGs were mainly involved in the process of transcription and signal transduction. Combining the enrichment analysis and protein interaction analysis, eight DEGs were predicted to control the growth of cattle by affecting the expression of growth-related hormones in the pituitary gland. In summary, our results suggested that SLC38A1, SLC38A3, DGKH, GNB4, GNAQ, ESR1, NPY, and GAL are candidates in the pituitary gland for regulating the growth of Yunling and Leiqiong cattle by regulating the secretion of growth-related hormones. This study may help researchers further understand the growth mechanisms and improve the artificial selection of zebu cattle.

Identification of Milk Fat Metabolism-Related Pathways of the Bovine Mammary Gland during Mid and Late Lactation and Functional Verification of the ACSL4 Gene.

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.

Polymorphisms of the ACSL1 Gene Influence Milk Production Traits and Somatic Cell Score in Chinese Holstein Cows.

Improving the quality of milk is a challenge for zootechnicians and dairy farms across the globe. Long-chain acyl-CoA synthetase 1 (ACSL1) is a significant member of the long-chain acyl-CoA synthetase gene family. It is widely found in various organisms and influences the lactation performance of cows, including fat percentage, milk protein percentage etc. Our study was aimed to investigate the genetic effects of single nucleotide polymorphisms (SNPs) in ACSL1 on milk production traits. Twenty Chinese Holstein cows were randomly selected to extract DNA from their blood samples for PCR amplification and sequencing to identify SNPs of the bovine ACSL1 gene, and six SNPs (5'UTR-g.20523C>G, g.35446C>T, g.35651G>A, g.35827C>T, g.35941G>A and g.51472C>T) were discovered. Then, Holstein cow genotyping (n = 992) was performed by Sequenom MassARRAY based on former SNP information. Associations between SNPs and milk production traits and somatic cell score (SCS) were analyzed by the least-squares method. The results showed that SNP g.35827C>T was in high linkage disequilibrium with g.35941G>A. Significant associations were found between SNPs and test-day milk yield (TDMY), fat content (FC), protein content (PC) and SCS (p < 0.05). Among these SNPs, SNP 5'UTR-g.20523C>G showed an extremely significant effect on PC and SCS (p < 0.01). The SNP g.35446C>T showed a statistically significant effect on FC, PC, and SCS (p < 0.01), and also TDMY (p < 0.05). The SNP g.35651G>A had a statistically significant effect on PC (p < 0.01). The SNP g.35827C>T showed a highly significant effect on TDMY, FC, and SCS (p < 0.01) and significantly influenced PC (p < 0.05). Lastly, SNP g.51472C>T was significantly associated with TDMY, FC, and SCS (p < 0.05). In summary, the pleiotropic effects of bovine ACSL1 for milk production traits were found in this paper, but further investigation will be required on the intrinsic correlation to provide a theoretical basis for the research on molecular genetics of milk quality traits of Holstein cows.

A Functional 3' UTR Polymorphism of FADS2 Affects Cow Milk Composition through Modifying Mir-744 Binding.

This study determined the associations of FADS2 c.1571G>A with milk FAs content and revealed that cows with the GG genotype had improved levels of delta-6 desaturase substrates (linoleic acid, C18:2n-6; p < 0.001) and decreased levels of desaturase products (gamma-linolenic acid, C18:3n-6; p < 0.001), indicating a reduction in FADS2 expression or delta-6 desaturase activity caused by this polymorphism. Computer alignment demonstrated that c.1571G>A occurred within a potential miR-744 binding site. When the c.1571G allele was present, the luciferase activity of reporter constructs was significantly suppressed by miR-744, while no such effect was observed with the A allele. Overexpression of miR-744 in bovine mammary epithelial cells (with the 1571GG genotype) downregulated FADS2 expression at both mRNA and protein levels. In contrast, inhibition of endogenous miR-744 with a specific inhibitor dramatically upregulated FADS2 expression. Taken together, these lines of evidence indicated that the c.1571A minor allele abolished the ability of miR-744 to bind FADS2, with a consequent increase in FADS2 expression levels and synthesis of omega-6 LC-PUFAs.

MicroRNA-106b Regulates Milk Fat Metabolism via ATP Binding Cassette Subfamily A Member 1 ( ABCA1) in Bovine Mammary Epithelial Cells.

Research on the mechanisms that regulate milk fat synthesis in dairy cows is essential to identify potential molecular targets that in the long term can help develop appropriate molecular breeding programs. Although some studies have revealed that microRNA (miRNA) affect lipid metabolism by targeting specific genes, joint analysis of miRNA and target mRNA data from bovine mammary tissue has revealed few clues regarding the underlying mechanisms controlling milk fat synthesis. The objective of the present study was to use high-throughput sequencing and bioinformatics analysis to identify miRNA and mRNA pairs and explore further their potential roles in regulating milk fat synthesis. A total of 233 pairs of negatively associated miRNA and mRNA pairs were detected. Among those, there were 162 pairs in which the miRNAs were down-regulated and the target mRNAs were up-regulated. Among the identified miRNA, miR-106b can bind the 3'-UTR of the ATP binding cassette subfamily A member 1 ( ABCA1), a gene previously identified as having a positive association with bovine milk fat synthesis. The overexpression of miR-106b in bovine mammary epithelial cells caused a decrease in triglyceride and cholesterol content while the inhibition of miR-106b increased triglyceride and cholesterol content, confirming its role in lipid metabolism. The present study allowed for the construction of a miR-106b- ABCA1 regulatory network map, thus providing a theoretical basis to target this gene in the molecular breeding of dairy cows.

Pathological Features of Staphylococcus aureus Induced Mastitis in Dairy Cows and Isobaric-Tags-for-Relative-and-Absolute-Quantitation Proteomic Analyses.

In part as a result of the production of an enterotoxin, Staphylococcus aureus is a highly infectious pathogen and is a considerable threat to food hygiene and safety. Clinical mastitis models were established by S. aureus nipple-tube perfusion. The influence of mastitis on the mammary-gland-tissue proteomic profile was investigated using isobaric tags for relative and absolute quantitation. In this study, healthy and mastitic tissues from different mammary-gland areas of the same dairy cows were screened, and differentially expressed proteins were identified. Bioinformatic analysis identified proteins related to the inflammation and immunization of dairy cows. Histology, immunoblotting, and immunohistochemical-staining analyses were used to determine the expression of PGLYRP1 and PTX3 proteins in the acquired mammary-gland-tissue samples. PGLYRP1 and PTX3 in mastitic mammary glands may be associated with tissue damage and immune responses to late stages of infection. This further contributes to the understanding of the molecular theory of the treatment of mastitis caused by S. aureus.