Search

The composition and metabolic profile of the ruminal microbiome have an impact on milk composition. To unravel the ruminal microbiome and metabolome affecting milk fat synthesis in dairy cows, 16S rRNA and internal transcribed spacer (ITS) gene sequencing, as well as ultra-high performance liquid chromatography-mass spectrometry (UHPLC-MS/MS) methods were used to investigate the significant differences in ruminal bacterial and fungal communities as well as metabolome among Chinese Holstein cows with contrasting milk fat contents under the same diet (H-MF 5.82 ± 0.41% vs. L-MF 3.60 ± 0.12%). Another objective was to culture bovine mammary epithelial cells (BMECs) to assess the effect of metabolites on lipid metabolism. Results showed that the acetate-to-propionate ratio and xylanase activity in ruminal fluid were both higher in H-MF. Microbiome sequencing identified 10 types of bacteria and four types of fungi differently abundant at the genus level. Metabolomics analysis indicated 11 different ruminal metabolites between the two groups, the majority of which were lipids and organic acids. Among these, lauric acid (LA) was enriched in fatty acid biosynthesis with its concentration in milk fat of H-MF cows being greater (217 vs. 156 mg per 100 g milk), thus, it was selected for an in vitro study with BMECs. Exogenous LA led to a marked increase in intracellular triglyceride (TG) content and lipid droplet formation, and it upregulated the mRNA abundance of fatty acid uptake and activation ( CD36 and ACSL1 ), TG synthesis ( DGAT1 , DGAT2 and GPAM ), and transcriptional regulation ( SREBP1 ) genes. Taken together, the greater relative abundance of xylan-fermenting bacteria and fungi, and lower abundance of bacteria suppressing short-chain fatty acid-producing bacteria or participating in fatty acid hydrogenation altered lipids and organic acids in the rumen of dairy cows. In BMECs, LA altered the expression of genes involved in lipid metabolism in mammary cells, ultimately promoting milk fat synthesis. Thus, it appears that this fatty acid plays a key role in milk fat synthesis.

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.

Udder conformation traits are one of the most economic traits in dairy cows, greatly affecting animal health, milk production, and producer profitability in the dairy industry. Genetic analysis of udder structure and scores have been developed in Holstein cattle. In our research, we conducted a genome-wide association study for five udder traits, including anterior udder attachment (AUA), central suspensory ligament (CSL), posterior udder attachment height (PUAH), posterior udder attachment width (PUAW), and udder depth (UD), in which the fixed and random model circulating probability unification (FarmCPU) model was applied for the association analysis. The heritability and the standard errors of these five udder traits ranged from 0.04 ± 0.00 to 0.49 ± 0.03. Phenotype data were measured from 1000 Holstein cows, and the GeneSeek Genomic Profiler (GGP) Bovine 100 K SNP chip was used to analyze genotypic data in Holstein cattle. For GWAS analysis, 984 individual cows and 84,407 single-nucleotide polymorphisms (SNPs) remained after quality control; a total of 18 SNPs were found at the GW significant threshold (p < 5.90 × 10−7). Many candidate genes were identified within 200kb upstream or downstream of the significant SNPs, which include MGST1, MGST2, MTUS1, PRKN, STXBP6, GRID2, E2F8, CDH11, FOXP1, SLF1, TMEM117, SBF2, GC, ADGRB3, and GCLC. Pathway analysis revealed that 58 Gene Ontology (GO) terms and 18 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were enriched with adjusted p values, and these GO terms and the KEGG pathway analysis were associated with biological information, metabolism, hormonal growth, and development processes. These results could give valuable biological information for the genetic architecture of udder conformation traits in dairy Holstein cattle.

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.

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.

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.

This study investigated the single nucleotide polymorphisms (SNPs) of Fatty acid desaturase 2 ( FADS2 ) gene and further explored their genetic effects on conventionally collected milk production traits in Chinese Holstein cows using 18,264 test-day records of 841 cows. One missense mutation c. 908 C > T (SNP site in the complementary DNA sequence), which caused an amino acid change from alanine to valine (294Ala > Val), and two 3' untranslated region (UTR) SNPs, c.1571 G > A and c.2776 A > G were finally identified. The SNP c.908 C > T was significantly associated with test-day milk yield, fat percentage and 305-day milk, fat and protein yield. In particular, the T allele of the SNP c.908 C > T showed a significant association with decreased somatic cell score (SCS) in the investigated population. Significant relationship between the SNP c.1571 G > A and 305-day milk yield showed that genotype GG was linked to the highest milk yield. Substituting the allele G for A at the c.2776 A > G locus resulted in a decrease of protein percentage. Our results demonstrated that FADS2 was an interesting candidate for selection to increase milk production and improve resistance against mastitis., Competing Interests: The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

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.