Mammary Gland Transcriptome and Proteome Modifications by Nutrient Restriction in Early Lactation Holstein Cows Challenged with Intra-Mammary Lipopolysaccharide

DNA methylation haplotype block signatures responding to Staphylococcus aureus subclinical mastitis and association with production and health traits

BACKGROUND

DNA methylation has been documented to play vital roles in diseases and biological processes. In bovine, little is known about the regulatory roles of DNA methylation alterations on production and health traits, including mastitis.

RESULTS

Here, we employed whole-genome DNA methylation sequencing to profile the DNA methylation patterns of milk somatic cells from sixteen cows with naturally occurring Staphylococcus aureus (S. aureus) subclinical mastitis and ten healthy control cows. We observed abundant DNA methylation alterations, including 3,356,456 differentially methylated cytosines and 153,783 differential methylation haplotype blocks (dMHBs). The DNA methylation in regulatory regions, including promoters, first exons and first introns, showed global significant negative correlations with gene expression status. We identified 6435 dMHBs located in the regulatory regions of differentially expressed genes and significantly correlated with their corresponding genes, revealing their potential effects on transcriptional activities. Genes harboring DNA methylation alterations were significantly enriched in multiple immune- and disease-related pathways, suggesting the involvement of DNA methylation in regulating host responses to S. aureus subclinical mastitis. In addition, we found nine discriminant signatures (differentiates cows with S. aureus subclinical mastitis from healthy cows) representing the majority of the DNA methylation variations related to S. aureus subclinical mastitis. Validation of seven dMHBs in 200 cows indicated significant associations with mammary gland health (SCC and SCS) and milk production performance (milk yield).

CONCLUSIONS

In conclusion, our findings revealed abundant DNA methylation alterations in milk somatic cells that may be involved in regulating mammary gland defense against S. aureus infection. Particularly noteworthy is the identification of seven dMHBs showing significant associations with mammary gland health, underscoring their potential as promising epigenetic biomarkers. Overall, our findings on DNA methylation alterations offer novel insights into the regulatory mechanisms of bovine subclinical mastitis, providing further avenues for the development of effective control measures.

Milk fat miRNome changes in response to LPS challenge in Holstein cows

Mastitis is an inflammatory disease in dairy cows, causing economic losses and reducing animal welfare. In order to contribute for the discovery of early and noninvasive indicators, our objective was to determine the effects of a lipopolysaccharide (LPS) challenge on the microRNA profile (miRNome) of milk fat, using microarray analyses in cows. Cows were fed a lactation diet at ad libitum intake (n = 6). At 27 ± 3 days in milk, cows were injected with 50 µg of LPS Escherichia coli in one healthy rear mammary quarter. Milk samples were collected just before LPS challenge (LPS-) and 6.5 h after LPS challenge (LPS +) from the same cows. Microarray analysis was performed using customized 8 × 60 K ruminant miRNA microarrays to compare LPS- to LPS + miRNome. In silico functional analyses were performed using OmicsNet and Mienturnet software. MiRNome comparison between LPS- and LPS + identified 37 differentially abundant miRNAs (q-value ≤ 0.05). The predicted target genes of the 37 differentially abundant miRNAs are mostly involved in cell life including apoptosis, cell cycle, proliferation and differentiation and in gene expression processes. MiRNome analyses suggest that miRNAs profile is related to the inflammation response of the mammary gland. In conclusion, we demonstrated that milk fat might be an easy and rapid source of miRNAs that are potential indicators of early mastitis in cows.

Identification of Hub Genes and Target miRNAs Crucial for Milk Production in Holstein Friesian Dairy Cattle

Dairy milk production is a quantitative trait that is controlled by many biological and environmental factors. This study employs a network-driven systems approach and clustering algorithm to uncover deeper insights into its genetic associations. We analyzed the GSE33680 dataset from the GEO database to understand the biological importance of milk production through gene expression and modules. In this study, we employed CytoNCA and ClusterONE plugins within Cytoscape for network analysis. Moreover, miRWalk software was utilized to detect miRNAs, and DAVID was employed to identify gene ontology and pathways. The results revealed 140 up-regulated genes and 312 down-regulated genes. In addition, we have identified 91 influential genes and 47 miRNAs that are closely associated with milk production. Through our examination of the network connecting these genes, we have found significant involvement in important biological processes such as calcium ion transit across cell membranes, the BMP signaling pathway, and the regulation of MAPK cascade. The conclusive network analysis further reveals that GAPDH, KDR, CSF1, PYGM, RET, PPP2CA, GUSB, and PRKCA are closely linked to key pathways essential for governing milk production. Various mechanisms can control these genes, making them valuable for breeding programs aiming to enhance selection indexes.

Goat mammary gland metabolism: An integrated Omics analysis to unravel seasonal weight loss tolerance

Seasonal weight loss (SWL), is a major limitation to animal production. In the Canary Islands, there are two dairy goat breeds with different levels of tolerance to SWL: Majorera (tolerant) and Palmera (susceptible). Our team has studied the response of these breeds to SWL using different Omics tools. The objective of this study was to integrate such results in a data driven approach and using dedicated tools, namely the DIABLO method. The outputs of our analysis mainly separate unrestricted from restricted goats. Metabolites behave as "hub" molecules, grouping interactions with several genes and proteins. Unrestricted goats upregulated protein synthesis, along with arginine catabolism and adipogenesis pathways, which are related with higher anabolic rates and a larger proportion of secretory tissue, in agreement with their higher milk production. Contrarily, restricted goats seemingly increased the synthesis of acetyl-CoA through serine and acetate conversion into pyruvate. This may have occurred to increase fatty acid synthesis and/or to use them as an energy source in detriment to glucose, which was more available in the diet of unrestricted goats. Lastly, restricted Palmera upregulated the expression of PEBP4 and GPD1 genes compared to all other groups, which might support their use as putative biomarkers for SWL susceptibility. SIGNIFICANCE: Seasonal weight loss (SWL) is a major issue influencing animal production in the tropics and Mediterranean. By studying its impact on the mammary gland of tolerant and susceptible dairy goat breeds, using Omics, we aim at surveying the tissue for possible biomarkers that reflect these traits. In this study, data integration of three Omics (transcriptomics, proteomics and metabolomics) was performed using bioinformatic tools, to relate putative biomarkers and evaluate all three levels of information; in a novel approach. This information can enhance selection programs, lowering the impact of SWL on food production systems.

Nutrition Genomics

This Special Issue is devoted to nutrition genomics, which is the characterization of the whole genome response to nutrients, in an effort to gather all the available pertinent information and to establish the foundation for a future encyclopedia of genomic responses driven by diets or nutrients [...]

Comparative transcriptome analysis of goat (Capra hircus) adipose tissue reveals physiological regulation of body reserve recovery after the peak of lactation

Adipose tissue is the energy storage organ providing energy to other tissues, including mammary gland, that supports the achievement of successive lactation cycles. Our objective was to investigate the ability of goats to restore body fat reserves by comparing lipogenic enzyme activities and by transcriptomic RNA-Seq data at two different physiological stages, mid- and post-lactation. Key lipogenic enzyme activities were higher in goat omental adipose tissue during mid-lactation (74 days in milk) than during the post-lactation period (300 days postpartum). RNA-Sequencing analysis revealed 19,271 expressed genes in the omental adipose tissue. The comparison between adipose transcriptome analysis from mid- and post-lactation goats highlighted 252 differentially expressed genes (p_adj < 0.05) between these two physiological stages. The differential expression of 11 genes was confirmed by RT-qPCR. Functional genomic analysis revealed that 31% were involved in metabolic processes among which 38% in lipid metabolism. Most of the genes involved in lipid synthesis and those in lipid transport and storage were upregulated in adipose tissue of mid- compared to post-lactation goats. In addition, adipose tissue plasticity was emphasized by genes involved in cellular signaling and tissue integrity. Network analyses also highlighted three key regulators of lipid metabolism (LEP, APOE and HNF4A) and a key target gene (VCAM1). The greatest lipogenic enzyme activities with the upregulation of genes involved in lipid metabolism highlighted a higher recovery of lipid reserves after the lactation peak than 4 months post-lactation. This study contributes to a better understanding of the molecular mechanisms controlling the body lipid reserves management during the successive lactations.

Milk metabolites and fatty acids as noninvasive biomarkers of metabolic status and energy balance in early-lactation cows

The objective was to study the effects of week of lactation (WOL) and experimental nutrient restriction on concentrations of selected milk metabolites and fatty acids (FA), and assess their potential as biomarkers of energy status in early-lactation cows. To study WOL effects, 17 multiparous Holstein cows were phenotyped from calving until 7 WOL while allowed ad libitum intake of a lactation diet. Further, to study the effects of nutrient restriction, 8 of these cows received a diet containing 48% straw (high-straw) for 4 d starting at 24 ± 3 days in milk (mean ± SD), and 8 cows maintained on the lactation diet were sampled to serve as controls. Blood and milk samples were collected weekly for the WOL data set, and daily from d -1 to 3 of nutrient restriction (or control) for the nutritional challenge data set. Milk β-hydroxybutyrate (BHB), isocitrate, glucose, glucose-6-phosphate (glucose-6P), galactose, glutamate, creatinine, uric acid, and N-acetyl-β-d-glucosaminidase activity (NAGase) were analyzed in p.m. and a.m. samples, and milk FA were analyzed in pooled p.m. and a.m. samples. Average energy balance (EB) per day ranged from -27 MJ/d to neutral when cows received the lactation total mixed ration, and from -109 to -87 ± 7 MJ/d for high-straw (least squares means ± standard error of the mean). Plasma nonesterified FA concentration was 1.67 ± 0.13 mM and BHB was 2.96 ± 0.39 mM on the d 3 of high-straw (least squares means ± standard error of the mean). Milk concentrations of BHB, glucose, glucose-6P, glutamate, and uric acid differed significantly between p.m. and a.m. milkings. Milk isocitrate, glucose-6P, creatinine, and NAGase decreased, whereas milk glucose and galactose increased with WOL. Changes in milk BHB, isocitrate, glucose, glucose-6P, and creatinine were concordant during early lactation and in response to nutrient restriction. Milk galactose and NAGase were modulated by WOL only, whereas glutamate and uric acid concentrations responded to nutrient restriction only. The high-straw increased milk concentrations of FA potentially mobilized from adipose tissue (e.g., C18:0 and cis-9 C18:1 and sum of odd- and branched-chain FA (OBCFA) with carbon chain greater than 16; ∑ OBCFA >C16), and decreased concentrations of FA synthesized de novo by the mammary gland (e.g., sum of FA with 6 to 15 carbons; ∑ C6:0 to C15:0). Similar observations were made during early lactation. Plasma nonesterified FA concentrations had the best single linear regression with EB (R² = 0.62). Milk isocitrate, Σ C6:0 to C15:0. and cis-9 C18:1 had the best single linear regressions with EB (R² ≥ 0.44). Milk BHB, isocitrate, galactose, glutamate, and creatinine explained up to 64% of the EB variation observed in the current study using multiple linear regression. Milk concentrations of ∑ C6:0 to C15:0, C18:0, cis-9 C18:1, and ∑ OBCFA >C16 presented some of the best correlations and regressions with other indicators of metabolic status, lipomobilization, and EB, and their responses were concordant during early lactation and during experimental nutrient restriction. Metabolites and FA secreted in milk may serve as noninvasive indicators of metabolic status and EB of early-lactation cows.

Nutritional Regulation of Mammary Gland Development and Milk Synthesis in Animal Models and Dairy Species

In mammals, milk is essential for the growth, development, and health. Milk quantity and quality are dependent on mammary development, strongly influenced by nutrition. This review provides an overview of the data on nutritional regulations of mammary development and gene expression involved in milk component synthesis. Mammary development is described related to rodents, rabbits, and pigs, common models in mammary biology. Molecular mechanisms of the nutritional regulation of milk synthesis are reported in ruminants regarding the importance of ruminant milk in human health. The effects of dietary quantitative and qualitative alterations are described considering the dietary composition and in regard to the periods of nutritional susceptibly. During lactation, the effects of lipid supplementation and feed restriction or deprivation are discussed regarding gene expression involved in milk biosynthesis, in ruminants. Moreover, nutrigenomic studies underline the role of the mammary structure and the potential influence of microRNAs. Knowledge from three lactating and three dairy livestock species contribute to understanding the variety of phenotypes reported in this review and highlight (1) the importance of critical physiological stages, such as puberty gestation and early lactation and (2) the relative importance of the various nutrients besides the total energetic value and their interaction.

Combined transcriptome and proteome analyses reveal differences in the longissimus dorsi muscle between Kazakh cattle and Xinjiang brown cattle

Objective

With the rapid development of proteomics sequencing and RNA sequencing technology, multi-omics analysis has become a current research hotspot. Our previous study indicated that Xinjiang brown cattle have better meat quality than Kazakh cattle. In this study, Xinjiang brown cattle and Kazakh cattle were used as the research objects.

Methods

Proteome sequencing and RNA sequencing technology were used to analyze the proteome and transcriptome of the longissimus dorsi muscle of the two breeds of adult steers (n = 3).

Results

In this project, 22,677 transcripts and 1,874 proteins were identified through quantitative analysis of the transcriptome and proteome. By comparing the identified transcriptome and proteome, we found that 1,737 genes were identified at both the transcriptome and proteome levels. The results of the study revealed 12 differentially expressed genes and proteins: troponin I1, crystallin alpha B, cysteine, and glycine rich protein 3, phosphotriesterase-related, myosin-binding protein H, glutathione s-transferase mu 3, myosin light chain 3, nidogen 2, dihydropyrimidinase like 2, glutamate-oxaloacetic transaminase 1, receptor accessory protein 5, and aspartoacylase. We performed functional enrichment of these differentially expressed genes and proteins. The Kyoto encyclopedia of genes and genomes results showed that these differentially expressed genes and proteins are enriched in the fatty acid degradation and histidine metabolism signaling pathways. We performed parallel reaction monitoring (PRM) verification of the differentially expressed proteins, and the PRM results were consistent with the sequencing results.

Conclusion

Our study provided and identified the differentially expressed genes and proteins. In addition, identifying functional genes and proteins with important breeding value will provide genetic resources and technical support for the breeding and industrialization of new genetically modified beef cattle breeds.

Population demographic history and population structure for Pakistani Nili-Ravi breeding bulls based on SNP genotyping to identify genomic regions associated with male effects for milk yield and body weight

The domestic Nili-Ravi water buffalo (Bubalus bubalis) is the best dairy animal contributing 68% to total milk production in Pakistan. In this study, we identified genome-wide single nucleotide polymorphisms (SNPs) to estimate various population genetic parameters such as diversity, pairwise population differentiation, linkage disequilibrium (LD) distribution and for genome-wide association study for milk yield and body weight traits in the Nili-Ravi dairy bulls that they may pass on to their daughters who are retained for milking purposes. The genotyping by sequencing approach revealed 13,039 reference genome-anchored SNPs with minor allele frequency of 0.05 among 167 buffalos. Population structure analysis revealed that the bulls were grouped into two clusters (K = 2), which indicates the presence of two different lineages in the Pakistani Nili-Ravi water buffalo population, and we showed the extent of admixture of these two lineages in our bull collection. LD analysis revealed 4169 significant SNP associations, with an average LD decay of 90 kb for these buffalo genome. Genome-wide association study involved a multi-locus mixed linear model for milk yield and body weight to identify genome-wide male effects. Our study further illustrates the utility of the genotyping by sequencing approach for identifying genomic regions to uncover additional demographic complexity and to improve the complex dairy traits of the Pakistani Nili-Ravi water buffalo population that would provide the lot of economic benefits to dairy industry.

The Synergism of PGN, LTA and LPS in Inducing Transcriptome Changes, Inflammatory Responses and a Decrease in Lactation as Well as the Associated Epigenetic Mechanisms in Bovine Mammary Epithelial Cells

Mastitis is usually caused by a variety of pathogenic bacteria that include both Gram-positive and Gram-negative bacteria. Lipopolysaccharide (LPS) is the pathogen-associated molecular pattern (PAMP) of Gram-negative bacteria, and peptidoglycan (PGN) and lipoteichoic acid (LTA) are those of Gram-positive bacteria. The effects of LPS, PGN and/or LTA on inflammatory response and lactation in bovine mammary epithelial cells (BMECs) are well studied, but the epigenetic mechanisms of their effects received less attention. Furthermore, since the three PAMPs are often simultaneously present in the udder of cows with mastitis, it has implications in practice to study their additive effects. The results show that co-stimulation of bovine mammary epithelial cells with PGN, LTA, and LPS induced a higher number of differentially expressed genes (DEGs) and greater expressions of inflammatory factors including interleukin (IL)-1β, IL-6, IL-8, tumor necrosis factor-α (TNF-α), chemokine (C-X-C motif) ligand (CXCL)1, and CXCL6. In addition, co-stimulation further increased DNA hypomethylation compared with sole LPS stimulation. Co-stimulation greatly decreased casein expression but did not further decrease histone acetylation levels and affect the activity of histone acetyltransferase (HAT) and histone deacetylase (HDAC), compared with sole LPS stimulation. Collectively, this study demonstrated that PGN, LTA, and LPS had an additive effect on inducing transcriptome changes and inflammatory responses in BMECs, probably through inducing a greater decrease in DNA methylation. Co-stimulation with PGN, LTA, and LPS decreased casein expression to a greater degree, but it might not be linked to histone acetylation and HAT and HDAC activity.

Milk metabolites as noninvasive indicators of nutritional status of mid-lactation Holstein and Montbéliarde cows

The objective was to investigate the effects of feed restriction on concentrations of selected milk metabolites in mid-lactation Holstein and Montbéliarde cows and to explore their correlations with energy balance and classic plasma and milk indicators of nutritional status. Eight Holstein and 10 Montbéliarde cows (165 ± 21 d in milk) underwent 6 d of feed restriction during which feed allowance was reduced to meet 50% of their net energy for lactation (NE_L) requirements. The experiment was divided in 4 periods: control (CON; d -3 to -1), restriction (RES; d 1 to 6), wk 1 (W1; d 7 to 13), and wk 2 (W2; d 14 to 18) after refeeding at ad libitum intake. Intake, milk production, energy balance and plasma metabolites were used to validate the feed restriction model. Concentrations of 7 milk metabolites: β-hydroxybutyrate (BHB), glucose, glucose-6-phosphate, isocitrate, glutamate, uric acid, and free amino groups were measured in morning milk samples, and fatty acids were measured in pooled p.m. and a.m. samples. Feed restriction induced a negative energy balance (-42.5 ± 4.4 MJ/d), increased plasma nonesterified fatty acids and BHB, and decreased plasma glucose concentrations. Feed restriction increased milk glucose-6-phosphate and isocitrate (+38% and +39%, respectively) and decreased milk BHB, glucose, glutamate, uric acid and free amino group concentrations (-20%, -57%, -65%, -42%, and -14%, respectively), compared with pre- restriction. Milk concentrations of medium-chain fatty acids (e.g., sum of C10 to C15) decreased and those of long chain (e.g., 18:0, cis-9 18:1) increased during restriction. Breed differences were not detected for the majority of variables. All studied milk metabolites were significantly correlated with energy balance (Spearman correlation = 0.48, 0.63, -0.31, -0.45, and 0.61 for BHB, glucose, glucose-6-phosphate, isocitrate, and glutamate, respectively). Milk glucose and glutamate showed the strongest correlations with plasma metabolites and milk FA associated with lipomobilization. These results suggest that milk metabolites may be used as noninvasive indicators of negative energy balance and metabolic status of dairy cows.