Metabolomics meets machine learning: Longitudinal metabolite profiling in serum of normal versus overconditioned cows and pathway analysis

Bile acid profiles and mRNA expression of bile acid-related genes in the liver of dairy cows with high versus normal body condition

Bile acids (BA) play a crucial role not only in lipid digestion but also in the regulation of overall energy homeostasis, including glucose and lipid metabolism. The aim of this study was to investigate BA profiles and mRNA expression of BA-related genes in the liver of high versus normal body condition in dairy cows. We hypothesized that body condition and the transition from gestation to lactation affect hepatic BA concentrations as well as the mRNA abundance of BA-related receptors, regulatory enzymes, and transporters. Therefore, we analyzed BA in the liver as well as the mRNA abundance of BA-related synthesizing enzymes, transporters, and receptors in the liver during the transition period in cows with different body conditions around calving. In a previously established animal model, 38 German Holstein cows were divided into groups with high body condition score (BCS) (HBCS; n = 19) or normal BCS (NBCS; n = 19) based on BCS and backfat thickness (BFT). Cows were fed diets aimed at achieving the targeted differences in BCS and BFT (NBCS: BCS <3.5, BFT <1.2 cm; HBCS: BCS >3.75, BFT >1.4 cm) until they were dried off at wk 7 before parturition. Both groups were fed identical diets during the dry period and subsequent lactation. Liver biopsies were taken at wk -7, 1, 3, and 12 relative to parturition. For BA measurement, a targeted metabolomics approach with LC-ESI-MS/MS was used to analyze BA in the liver. The mRNA abundance of targeted genes related to BA-synthesizing enzymes, transporters, and receptors in the liver was analyzed using microfluidic quantitative PCR. In total, we could detect 14 BA in the liver: 6 primary and 8 secondary BA, with glycocholic acid (GCA) being the most abundant one. The increase of glycine-conjugated BA after parturition, in parallel to increasing serum glycine concentrations may originate from an enhanced mobilization of muscle protein to meet the high nutritional requirements in early lactating cows. Higher DMI in NBCS cows compared with HBCS cows was associated with higher liver BA concentrations such as GCA, deoxycholic acid (DCA), and cholic acid (CA). The mRNA abundance of BA-related enzymes measured herein suggests the dominance of the alternative signaling pathway in the liver of HBCS cows. Overall, BA profiles and BA metabolism in the liver depend on both, the body condition and lactation-induced effects in periparturient dairy cows.

Longitudinal characterization of the metabolome of dairy cows transitioning from one lactation to the next: Investigations in the liver

This study aimed to investigate the metabolic changes in the livers of dairy cows from 1 wk before dry off to 1 wk after calving. Twelve high-yielding Holstein cows were included in a longitudinal study and housed in a tiestall barn. The cows were dried off at 6 wk before the expected calving date (dry period length = 42 d). During the entire lactation, the cows were milked twice daily at 0600 and 1700 h. Liver biopsies were taken from each cow at 4 different times: wk -7 (before drying off), -5 (after drying off), -1 and +1 relative to calving. A targeted metabolomics approach was performed by liquid chromatography and flow injection with electrospray ionization triple quadrupole mass spectrometry using the MxP Quant 500 kit (Biocrates Life Sciences AG). A total of 185 metabolites in the liver were used for the final data analysis. Principal component analysis revealed a clear separation by days of sampling, indicating a notable shift in metabolic phenotype from late lactation to the dry period and further changes after calving. Changes were observed in several classes of compounds, including AA and biogenic amines. In particular, the changes in acylcarnitines (AcylCN), phosphatidylcholines (PC), sphingomyelins (SM), and bile acids (BA) indicated extensive remodeling of the hepatic lipidome. The changes in AcylCN concentrations in early lactation suggest incomplete fatty acid oxidation in the liver, possibly indicating mitochondrial dysfunction or enzymatic imbalance. In addition, the changes in PC and SM species in early lactation indicate altered cell membrane composition, which may affect cell signaling and functionality. In addition, changes in BA concentrations and profiles indicate dynamic adaptations in BA synthesis, as well as lipid digestion and absorption during the observation period. In particular, principal component analysis showed an overlapping distribution of liver metabolites in primiparous and multiparous cows, indicating no significant difference between these groups. In addition, Volcano plots showed similar liver metabolism between primiparous and multiparous cows, with no significant fold changes (>1.5) in any metabolite at significant P-values (false discovery rate <0.05). These results provide valuable insight into the physiological ranges of liver metabolites during dry period and calving in healthy dairy cows and should contribute to the design and interpretation of future metabolite-based studies of the transition dairy cow.

Longitudinal characterization of serum metabolome and lipidome reveals that the ceramide profile is associated with metabolic health in early postpartum cows experiencing different lipolysis

Reduced feed intake in early lactation prompts increased fat mobilization to meet dairy cows' energy needs for milk production. The increased lipolysis in cows presents significant health risks with unclear mechanisms. The objectives of our study were to compare the longitudinal profiles of metabolites and lipids of serum from high and low-lipolysis cows. Forty multiparous Holstein dairy cows were enrolled in the retrospective study. Serum samples were collected on d 7 before expected calving, as well as on d 5, d 7, d 14, and d 21 postpartum. Dairy cows were grouped according to mean serum nonesterified fatty acids on d 5 and 7 after parturition as low (<0.600 mmol/L; n = 8; LFM) and high (>0.750 mmol/L; n = 8; HFM), indicating fat mobilization during early lactation. Lactational performance and serum metabolic parameters related to glucose and lipid metabolism, liver functions, oxidative status, and inflammatory responses were determined. Serum samples were subjected to LC-MS-based metabolomics and lipidomics. Despite differences in postpartum BW change, there were no observed variations in milk yield and composition between 2 groups. Serum BHBA, glucose, leptin, aspartate aminotransferase, IL-6, and TNF-α were greater in cows with HFM than in LFM. Serum adiponectin, revised quantitative insulin sensitivity check index and albumin were lower in cows with HFM than LFM. Intensified fat mobilization in the HFM cows came along with reduced estimated insulin sensitivity, impaired liver functions, and increased oxidative stress and inflammatory responses. Differences in metabolic patterns were observed across the transition period when comparing serum blood matrices (e.g., in different amino acids, acylcarnitines, and sphingolipids). The serum metabolome of the HFM cows was characterized by higher concentrations of glycine, acylcarnitines, carnosine, Cer(d20:0/18:0), Cer(d18:1/16:0), and Cer(t18:0/24:0) compared with LFM. The differential serum metabolites and lipids at different sampling times during the peripartum period were enriched in the sphingolipid metabolism. Differences in serum metabolic status parameters suggest that cows adopt varied metabolic adaptation strategies to cope with energy deficits postpartum. Our investigation found a comprehensive remodeling of the serum metabolic profiles in transition dairy cattle, highlighting the significance of alterations in sphingolipid species, as they play a crucial role in insulin resistance and metabolic disorders.

Prostaglandin D2 synthase controls Schwann cells metabolism

We previously reported that in the absence of Prostaglandin D2 synthase (L-PGDS) peripheral nerves are hypomyelinated in development and that with aging they present aberrant myelin sheaths. We now demonstrate that L-PGDS expressed in Schwann cells is part of a coordinated program aiming at preserving myelin integrity. In vivo and in vitro lipidomic, metabolomic and transcriptomic analyses confirmed that myelin lipids composition, Schwann cells energetic metabolism and key enzymes controlling these processes are altered in the absence of L-PGDS. Moreover, Schwann cells undergo a metabolic rewiring and turn to acetate as the main energetic source. Further, they produce ketone bodies to ensure glial cell and neuronal survival. Importantly, we demonstrate that all these changes correlate with morphological myelin alterations and describe the first physiological pathway implicated in preserving PNS myelin. Collectively, we posit that myelin lipids serve as a reservoir to provide ketone bodies, which together with acetate represent the adaptive substrates Schwann cells can rely on to sustain the axo-glial unit and preserve the integrity of the PNS.

Liver transcriptome profiles of dairy cows with different serum metabotypes

In a previously established animal model, 38 multiparous Holstein cows were assigned to 2 groups fed different diets to achieve either a normal (NBCS) or high (HBCS) body condition score (BCS) and backfat thickness (BFT) until dry-off at -49 d before calving (NBCS: BCS <3.5 [3.02 ± 0.24) and BFT <1.2 cm [0.92 ± 0.21]; HBCS: BCS >3.75 [3.82 ± 0.33] and BFT >1.4 cm [2.36 ± 0.35], mean ± SD). The groups were also stratified for comparable milk yields (NBCS: 10,361 ± 302 kg; HBCS: 10,315 ± 437 kg; mean ± SD). The cows were then fed the same diet during the dry period and subsequent lactation, maintaining the differences in BFT and BCS throughout the study. Using the serum metabolomics data, we created a classification model that identified different metabotypes. Machine learning classifiers revealed a distinct cluster labeled HBCS-PN (HBCS predicted normal BCS) among over-conditioned cows. These cows showed higher feed intake and better energy balance than the HBCS-PH (high BCS predicted high BCS) group, while milk yield was similar. The aim of this study was to investigate the changes in the hepatic transcriptome of cows differing in serum-metabotype postpartum. We performed hepatic transcriptome analysis in cows from 3 metabolic clusters: HBCS-PH (n = 8), HBCS-PN (n = 6), and normal BCS predicted normal BCS (NBCS-PN, n = 8) on d 21 (±2) postpartum. Liver tissue from cows expressed a total of 13,118 genes aligned with the bovine genome. A total of 48 differentially expressed genes (DEG; false discovery rate ≤0.1 and fold-change >1.5) were found between NBCS-PN and HBCS-PH cows, whereas 24 DEG (14 downregulated and 10 upregulated) were found between HBCS-PN and HBCS-PH cows. The downregulated DEG (n = 31) in NBCS-PN cows compared with HBCS-PH cows are involved in biosynthetic processes such as lipid, lipoprotein, and cholesterol synthesis (e.g., APOA1, MKX, RPL3L, CANT1, CHPF, FUT1, ZNF696), cell organization, biogenesis, and localization (e.g., SLC12A8, APOA1, BRME1, RPL3L, STAG3, FBXW5, TMEM120A, SLC16A5, FGF21), catabolic processes (e.g., BREH1, MIOX, APOBEC2, FBXW5, NUDT16), and response to external stimuli (e.g., APOA1, FGF21, TMEM120A, FNDC4), whereas upregulated DEG (n = 17) are related to signal transduction and cell motility (e.g., RASSF2, ASPN, SGK1, KIF7, ZEB2, MAOA, ACKR4, TCAF1), suggesting altered metabolic adaptations during lactation. Our results showed 24 DEG between HBCS-PN and HBCS-PH in the liver. The expression of SLC12A8, SLC16A5, FBXW5, OSGIN1, LAMA3, KDELR3, OR4X17, and INHBE, which are responsible for regulating cellular processes was downregulated in HBCS-PN cows compared with HBCS-PH cows. In particular, the downregulation of SLC12A8 and SLC16A5 expression in HBCS-PN cows indicates lower metabolic load and reduced need for NAD+ biosynthesis to support mitochondrial respiratory processes. The upregulation of MAOA, ACKR4, KIF27, SFRP1, and CAV2 in the liver of HBCS-PN cows may indicate adaptive mechanisms to maintain normal liver function in response to increased metabolic demands from over-conditioning. These molecular differences underscore the existence of distinct metabolic types in cows and provide evidence for the role of the liver in shaping different metabolic patterns.

Differences of serum glucose and lipid metabolism and immune parameters and blood metabolomics regarding the transition cows in the antepartum and postpartum period

This study aims to investigate differences in metabolism regarding the transition cows. Eight cows were selected for the test. Serum was collected on antepartum days 14th (ap14) and 7th (ap7) and postpartum days 1st (pp1), 7th (pp7), and 14th (pp14) to detect biochemical parameters. The experiment screened out differential metabolites in the antepartum (ap) and postpartum (pp) periods and combined with metabolic pathway analysis to study the relationship and role between metabolites and metabolic abnormalities. Results: (1) The glucose (Glu) levels in ap7 were significantly higher than the other groups (p < 0.01). The insulin (Ins) levels of ap7 were significantly higher than pp7 (p = 0.028) and pp14 (p < 0.01), and pp1 was also significantly higher than pp14 (p = 0.016). The insulin resistance (HOMA-IR) levels of ap7 were significantly higher than ap14, pp7, and pp14 (p < 0.01). The cholestenone (CHO) levels of ap14 and pp14 were significantly higher than pp1 (p < 0.01). The CHO levels of pp14 were significantly higher than pp7 (p < 0.01). The high density lipoprotein cholesterol (DHDL) levels of pp1 were significantly lower than ap14 (p = 0.04), pp7 (p < 0.01), and pp14 (p < 0.01), and pp14 was also significantly higher than ap14 and ap7 (p < 0.01). (2) The interferon-gamma (IFN-γ) and tumor necrosis factor α (TNF-α) levels of ap7 were significantly higher than pp1 and pp7 (p < 0.01); the immunoglobulin A (IgA) levels of pp1 were significantly higher than ap7 and pp7 (p < 0.01); the interleukin-4 (IL-4) levels of pp7 were significantly higher than ap7 and pp1 (p < 0.01), the interleukin-6 (IL-6) levels of ap7 and pp1 were significantly higher than pp7 (p < 0.01). (3) Metabolomics identified differential metabolites mainly involved in metabolic pathways, such as tryptophan metabolism, alpha-linolenic acid metabolism, tyrosine metabolism, and lysine degradation. The main relevant metabolism was concentrated in lipid and lipid-like molecules, organic heterocyclic compounds, organic acids, and their derivatives. The results displayed the metabolic changes in the transition period, which laid a foundation for further exploring the mechanism of metabolic abnormalities in dairy cows in the transition period.

Examination of Serum Metabolome Altered by Dietary Carbohydrate, Milk Protein, and Soy Protein Interventions Identified Novel Metabolites Associated with Blood Pressure: The ProBP Trial

SCOPE

This study aims to discover metabolites of dietary carbohydrate, soy and milk protein supplements and evaluate their roles in blood pressure (BP) regulation in the protein and blood pressure (ProBP), a cross-over trial.

METHODS AND RESULTS

Plasma metabolites are profiled at pre-trial baseline and after 8 weeks of supplementation with carbohydrate, soy protein, and milk protein, respectively, among 80 ProBP participants. After Bonferroni correction (α = 6.49 × 10-4 ), dietary interventions significantly changed 40 metabolites. Changes of erucate (22:1n9), an omega-9 fatty acid, are positively associated with systolic BP changes (Beta = 1.90, p = 6·27 × 10-4 ). This metabolite is also associated with higher odds of hypertension among 1261 participants of an independent cohort (odds ratio per unit increase = 1.34; 95% confidence interval: 1.07-1.68). High levels of acylcholines dihomo-linolenoyl-choline (p = 4.71E-04) and oleoylcholine (p = 3.48E-04) at baseline predicted larger BP lowering effects of soy protein. Increasing cheese intake during the trial, as reflected by isobutyrylglycine and isovalerylglycine, reduces the BP lowering effect of soy protein.

CONCLUSIONS

The study identifies molecular signatures of dietary interventions. Erucate (22:1n9) increases systolic BP. Acylcholine enhances and cheese intake reduces the BP lowering effect of soy protein supplement.

Venlafaxine's effect on resilience to stress is associated with a shift in the balance between glucose and fatty acid utilization

Brain metabolism is a fundamental process involved in the proper development of the central nervous system and in the maintenance of the main higher functions in humans. As consequence, energy metabolism imbalance has been commonly associated to several mental disorders, including depression. Here, by employing a metabolomic approach, we aimed to establish if differences in energy metabolite concentration may underlie the vulnerability and resilience in an animal model of mood disorder named chronic mild stress (CMS) paradigm. In addition, we have investigated the possibility that modulation of metabolite concentration may represent a pharmacological target for depression by testing whether repeated treatment with the antidepressant venlafaxine may normalize the pathological phenotype by acting at metabolic level. The analyses were conducted in the ventral hippocampus (vHip) for its key role in the modulation of anhedonia, a core symptom of patients affected by depression. Interestingly, we showed that a shift from glycolysis to beta oxidation seems to be responsible for the vulnerability to chronic stress and that vHip metabolism contributes to the ability of the antidepressant venlafaxine to normalize the pathological phenotype, as shown by the reversal of the changes observed in specific metabolites. These findings may provide novel perspectives on metabolic changes that could serve as diagnostic markers and preventive strategies for the early detection and treatment of depression as well as for the identification of potential drug targets.

Untargeted metabolomics and lipidomics to assess plasma metabolite changes in dairy goats with subclinical hyperketonemia

Subclinical hyperketonemia (SCHK) is the major metabolic disease observed during the transition period in dairy goats, and is characterized by high plasma levels of nonesterified fatty acids (NEFA) and β-hydroxybutyrate (BHB). However, no prior study has comprehensively assessed metabolomic profiles of dairy goats with SCHK. Plasma samples were collected within 1 h after kidding from SCHK goats (BHB concentration >0.8 mM, n = 7) and clinically healthy goats (BHB concentration <0.8 mM, n = 7) with similar body condition score (2.75 ± 0.15, mean ± standard error of the mean) and parity (primiparous). A combination of targeted and untargeted mass spectrometric approaches was employed for analyzing the various changes in the plasma lipidome and metabolome. Statistical analyses were performed using the GraphPad Prism 8.0, SIMCA-P software (version 14.1), and R packages (version 4.1.3). Plasma aminotransferase, nonesterified fatty acids, and BHB concentrations were greater in the SCHK group, but plasma glucose concentrations were lower. A total of 156 metabolites and 466 lipids were identified. The analysis of untargeted metabolomics data by principal component analysis and orthogonal partial least squares discriminant analysis revealed a separation between SCHK and clinically healthy goats. According to the screening criteria (unpaired t-test, P < 0.05), 30 differentially altered metabolites and 115 differentially altered lipids were detected. Pathway enrichment analysis identified citrate cycle, alanine, aspartate and glutamate metabolism, glyoxylate and dicarboxylate metabolism, and phenylalanine metabolism as significantly altered pathways. A greater concentration of plasma isocitric acid and cis-aconitic acid levels was observed in SCHK goats. In addition, AA such as lysine and isoleucine were greater, whereas alanine and phenylacetylglycine were lower in SCHK dairy goats. Dairy goats with SCHK also exhibited greater oleic acid, acylcarnitine, and phosphatidylcholine and lower choline and sphingomyelins. Acylcarnitines, oleic acid, and tridecanoic acid displayed positive correlations with several lipid species. Alanine, hippuric acid, and histidinyl-phenylalanine were negatively correlated with several lipids. Overall, altered metabolites in SCHK dairy goats indicated a more severe degree of negative energy balance. Data also indicated an imbalance in the tricarboxylic acid (TCA) cycle, lipid metabolism, and AA metabolism. The findings provide a more comprehensive understanding of the pathogenesis of SCHK in dairy goats.

Impact of dry-off and lyophilized Aloe arborescens supplementation on plasma metabolome of dairy cows

Positive effects have been observed as a result of Aloe arborescens supplementation in the dry-off phase in dairy cows. Metabolomic approaches can provide additional information about animal physiology. Thus, we characterized plasma metabolome around dry-off in 12 cows supplemented (AL) or not (CTR) with 10 g/d of lyophilized A. arborescens with an untargeted metabolomic approach. Overall, 1658 mass features were annotated. Regardless of treatment, multivariate statistics discriminated samples taken before and after dry-off. Overall, 490 metabolites were different between late lactation and early dry period, of which 237 were shared between AL and CTR. The most discriminant compounds (pentosidine and luteolin 7-O-glucoside) were related to the more fibrous diet. Pathway analysis indicated that pyrimidine and glycerophospholipid metabolisms were down-accumulated, suggesting reduced rumen microbial activity and liver load. Samples from AL were discriminated from CTR either the day of dry-off or 7 days after. At dry-off, aloin and emodin were the most discriminant metabolites, indicating that Aloe's bioactive compounds were absorbed. Seven days later, 534 compounds were different between groups, and emodin was among the most impacted. Pathway analysis highlighted that glycerophospholipid, pyrimidine, and folate metabolisms were affected. These results might indicate that Aloe has positive effects on liver function and a modulatory effect on rumen fermentation.

Serum 1H nuclear magnetic resonance-based metabolomics of sole lesion development in Holstein cows

Sole hemorrhage and sole ulcers, referred to as sole lesions, are important causes of lameness in dairy cattle. We aimed to compare the serum metabolome of dairy cows that developed sole lesions in early lactation with that of cows that remained unaffected. We prospectively enrolled a cohort of 1,169 Holstein dairy cows from a single dairy herd and assessed animals at 4 time points: before calving, immediately after calving, early lactation, and late lactation. Sole lesions were recorded by veterinary surgeons at each time point, and serum samples were collected at the first 3 time points. Cases were defined by the presence of sole lesions in early lactation and further subdivided by whether sole lesions had been previously recorded; unaffected controls were randomly selected to match cases. Serum samples from a case-control subset of 228 animals were analyzed with proton nuclear magnetic resonance spectroscopy. Spectral signals, corresponding to 34 provisionally annotated metabolites and 51 unlabeled metabolites, were analyzed in subsets relating to time point, parity cohort, and sole lesion outcome. We used 3 analytic methods (partial least squares discriminant analysis, least absolute shrinkage and selection operator regression, and random forest) to determine the predictive capacity of the serum metabolome and identify informative metabolites. We applied bootstrapped selection stability, triangulation, and permutation to support the inference of variable selection. The average balanced accuracy of class prediction ranged from 50 to 62% depending on the subset. Across all 17 subsets, 20 variables had a high probability of being informative; those with the strongest evidence of being associated with sole lesions corresponded to phenylalanine and 4 unlabeled metabolites. We conclude that the serum metabolome, as characterized by proton nuclear magnetic resonance spectroscopy, does not appear able to predict sole lesion presence or future development of lesions. A small number of metabolites may be associated with sole lesions although, given the poor prediction accuracies, these metabolites are likely to explain only a small proportion of the differences between affected and unaffected animals. Future metabolomic studies may reveal underlying metabolic mechanisms of sole lesion etiopathogenesis in dairy cows; however, the experimental design and analysis need to effectively control for interanimal and extraneous sources of spectral variation.

Metabolic Profiling in Ruminant Diagnostics

A herd-based approach and interpretative perspective is necessary in using metabolic profile testing in contrast to individual animal disease diagnostics. Metabolic profile testing requires formulating a question to be answered, followed by the appropriate selection of animals for testing. A range of blood analytes and nutrients can be determined with newer biomarkers being developed. Sample collection and handling and herd-based reference criteria adjusted to time relative to parturition are critical for interpretation. The objective of this article is to review the concepts and practical applications of metabolic profile testing in ruminants.

Brown goat yogurt: Metabolomics, peptidomics, and sensory changes during production

Brown goat milk products have gained popularity for their unique taste and flavor. The emergence of chain-reversal phenomenon makes the design and development of goat milk products gradually tend to a consumer-oriented model. However, the precise mechanism of how browning and fermentation process causes characteristics is not clear. In an effort to understand how the treatments potentially lead to certain metabolite profile changes in goat milk, comprehensive, quantitative metabolomics and peptidomics analysis of goat milk samples after browning and fermentation were undertaken. An intelligent hybrid z-score standardization-principal components algorithm-multimodal denoizing autoencoder was used for feature fusion and hidden layer fusion in high-dimensional variable space. The fermentation process significantly improved the flavor of brown goat yogurt through the tricarboxylic acid-urea-glycolysis composite pathway. Bitter peptides HPFLEWAR, PPGLPDKY, and PPPPPKK have strong interactions with both putative dipeptidyl peptidase IV and angiotensin-converting enzyme, proving that brown goat yogurt can be considered as effective provider of potential putative dipeptidyl peptidase IV and angiotensin-converting enzyme inhibitors. The level of health-promoting bioactive components and sensory contributed to consumer selection. The proposed multimodal data integrative analysis platform was applicable to explain the effect of the dynamic changes of metabolites and peptides on consumer preferences.

Filling gaps in animal welfare assessment through metabolomics

Sustainability has become a central issue in Italian livestock systems driving food business operators to adopt high standards of production concerning animal husbandry conditions. Meat sector is largely involved in this ecological transition with the introduction of new label claims concerning the defense of animal welfare (AW). These new guarantees referred to AW provision require new tools for the purpose of authenticity and traceability to assure meat supply chain integrity. Over the years, European Union (EU) Regulations, national, and international initiatives proposed provisions and guidelines for assuring AW introducing requirements to be complied with and providing tools based on scoring systems for a proper animal status assessment. However, the comprehensive and objective assessment of the AW status remains challenging. In this regard, phenotypic insights at molecular level may be investigated by metabolomics, one of the most recent high-throughput omics techniques. Recent advances in analytical and bioinformatic technologies have led to the identification of relevant biomarkers involved in complex clinical phenotypes of diverse biological systems suggesting that metabolomics is a key tool for biomarker discovery. In the present review, the Five Domains model has been employed as a vademecum describing AW. Starting from the individual Domains-nutrition (I), environment (II), health (III), behavior (IV), and mental state (V)-applications and advances of metabolomics related to AW setting aimed at investigating phenotypic outcomes on molecular scale and elucidating the biological routes most perturbed from external solicitations, are reviewed. Strengths and weaknesses of the current state-of-art are highlighted, and new frontiers to be explored for AW assessment throughout the metabolomics approach are argued. Moreover, a detailed description of metabolomics workflow is provided to understand dos and don'ts at experimental level to pursue effective results. Combining the demand for new assessment tools and meat market trends, a new cross-strategy is proposed as the promising combo for the future of AW assessment.

Invited review: Assessment of body condition score and body fat reserves in relation to insulin sensitivity and metabolic phenotyping in dairy cows

The purpose of this article is to review body condition scoring and the role of body fat reserves in relation to insulin sensitivity and metabolic phenotyping. This article summarizes body condition scoring assessment methods and the differences between subcutaneous and visceral fat depots in dairy cows. The mass of subcutaneous and visceral adipose tissue (AT) changes significantly during the transition period; however, metabolism and intensity of lipolysis differ between subcutaneous and visceral AT depots of dairy cows. The majority of studies on AT have focused on subcutaneous AT, and few have explored visceral AT using noninvasive methods. In this systematic review, we summarize the relationship between body fat reserves and insulin sensitivity and integrate omics research (e.g., metabolomics, proteomics, lipidomics) for metabolic phenotyping of cows, particularly overconditioned cows. Several studies have shown that AT insulin resistance develops during the prepartum period, especially in overconditioned cows. We discuss the role of AT lipolysis, fatty acid oxidation, mitochondrial function, acylcarnitines, and lipid insulin antagonists, including ceramide and glycerophospholipids, in cows with different body condition scoring. Nonoptimal body conditions (under- or overconditioned cows) exhibit marked abnormalities in metabolic and endocrine function. Overall, reducing the number of cows with nonoptimal body conditions in herds seems to be the most practical solution to improve profitability, and dairy farmers should adjust their management practices accordingly.

Associations of automated body condition scores at dry-off and through early lactation with milk yield of Holstein cows

The objective of this study was to analyze the associations of body condition score (BCS) and BCS change (∆BCS) during the dry period and the first 100 d of lactation with daily milk yield. Examining the involvement of health status in the associations between BCS and milk yield was a secondary objective of this research. Data included 12,042 lactations in 7,626 Holstein cows calving between April 2019 and January 2022 in a commercial dairy operation located in Colorado, USA. BCSs were generated daily by an automated BCS camera system located at the exit of the milking parlor. The assessment points selected for this study were dry-off (BCSdry), calving (BCS1), 7 DIM (BCS7), 14 DIM (BCS14), 21 DIM (BCS21), and nadir (nBCS; defined as the lowest daily BCS from calving to 100 DIM). Subsequently, these BCS were categorized considering quartiles (Q1 = 25% lowest BCS; Q4 = 25% greatest BCS), separately for primiparous and multiparous cows. Changes in BCS were calculated from dry-off to calving (multiparous); and from calving to 7 DIM, 14 DIM, 21 DIM, and nadir and assigned into quartile categories considering Q1 as the 25% of cows with the greatest decrease of BCS. Lactations were classified based on the number of health events before nadir as healthy, affected by one event, or having multiple events. Data were examined in primiparous and multiparous cows separately using ANOVA. The least square means for daily milk at 60 DIM and 305 DIM were calculated by category of BCS and ∆BCS at multiple time points and time periods. Subsequently, lactation curves were created by BCS and ∆BCS categories and by health status. Multivariable models included calving season and BCS1 as covariables. The largest differences in milk yield among categories of BCS and ∆BCS were identified for BCS originated at nadir and for the ∆BCS between calving and nadir. The differences in average daily milk yield between cows in the lowest and the greatest nBCS category (Q1 vs. Q4) were 3.3 kg/d (60 DIM) and 3.4 kg/d (305 DIM) for primiparous cows and 2.4 kg/d (60 DIM) and 2.1 kg/d (305 DIM) for multiparous cows. During the period from calving to nadir, primiparous cows in Q1 (greatest decrease of BCS) produced 4.3 kg/d (60 DIM) and 3.8 kg/d (305 DIM) more than cows in Q4. For multiparous cows, the differences were 3.0 kg/d (60 DIM) and 1.9 kg/d (305 DIM) in favor of Q1 cows. Overall, the associations between BCS and ∆BCS categories and milk yield were not consistent across time and they depended on the parity category. Nonetheless, as the assessment of BCS and ∆BCS approached the nadir, the association between greater milk yield and lower BCS or greater reduction in BCS became more evident.

Untargeted metabolomics confirms the association between plasma branched chain amino acids and residual feed intake in beef heifers

This study explored plasma biomarkers and metabolic pathways underlying feed efficiency measured as residual feed intake (RFI) in Charolais heifers. A total of 48 RFI extreme individuals (High-RFI, n = 24; Low-RFI, n = 24) were selected from a population of 142 heifers for classical plasma metabolite and hormone quantification and plasma metabolomic profiling through untargeted LC-MS. Most efficient heifers (Low-RFI) had greater (P = 0.03) plasma concentrations of IGF-1 and tended to have (P = 0.06) a lower back fat depth compared to least efficient heifers. However, no changes were noted (P ≥ 0.10) for plasma concentrations of glucose, insulin, non-esterified fatty acids, β-hydroxybutyrate and urea. The plasma metabolomic dataset comprised 3,457 ions with none significantly differing between RFI classes after false discovery rate correction (FDR > 0.10). Among the 101 ions having a raw P < 0.05 for the RFI effect, 13 were putatively annotated by using internal databases and 6 compounds were further confirmed with standards. Metabolic pathway analysis from these 6 confirmed compounds revealed that the branched chain amino acid metabolism was significantly (FDR < 0.05) impacted by the RFI classes. Our results confirmed for the first time in beef heifers previous findings obtained in male beef cattle and pointing to changes in branched-chain amino acids metabolism along with that of body composition as biological mechanisms related to RFI. Further studies are warranted to ascertain whether there is a cause-and-effect relationship between these mechanisms and RFI.

Metabolomic exploration of the effects of habituation to livestock trailer and extended transportation in goats

Goats raised for meat production are often transported long distances. Twelve-month-old male Spanish goats were used to determine the effects of habituation to trailers on plasma metabolomic profiles when transported for extended periods. In a split-plot design, 168 goats were separated into two treatment (TRT; whole plot) groups and maintained on two different paddocks. Concentrate supplement was fed to one group inside two livestock trailers (habituated group, H), while the other group received the same quantity of concentrate, but not inside the trailers (non-habituated, NH). Goats were subjected to a 10-h transportation stress in 4 replicates (n = 21 goats/replicate/TRT) after 4 weeks of habituation period. Blood samples were collected prior to loading, 20 min after loading (0 h), and at 2, 4, 6, 8, and 10 h of transportation (Time; subplot). A targeted quantitative metabolomics approach was employed to analyze the samples. The data were analyzed using R software and MIXED procedures in SAS. Several amino acids (alanine, serine, glycine, histidine, glutamate, trans-hydroxyproline, asparagine, threonine, methylhistidine, ornithine, proline, leucine, tryptophan) were higher (p < 0.05) in the H group compared to the NH group. Six long-chain acylcarnitines were higher (p < 0.05), while free (C0) and short-chain (C3, C5) carnitines were lower (p < 0.05) in the NH goats compared to the H goats. In general, amino acid concentrations decreased and long-chain acylcarnitine (>C10) levels increased with transportation time (p < 0.05). Butyric acid, α-ketoglutaric acid, and α-aminoadipic acid concentrations were lower (p < 0.05) and β-hydroxybutyric acid concentrations were higher in the NH goats compared to the H goats. Plasma glucose, non-esterified fatty acid (NEFA) and urea nitrogen concentrations were significantly influenced by Time (p < 0.01). Plasma NEFA concentrations were significantly lower (p < 0.01) in the H group than the NH group. Habituation to trailers can be beneficial in enhancing stress coping abilities in goats due to higher concentrations of metabolites such as butyrate and certain amino acids that support antioxidant activities and immune function. Plasma long-chain acylcarnitines may be good indicators of stress during long-distance transportation in goats.

Metabolomics and modelling approaches for systems metabolic engineering

Metabolic engineering involves the manipulation of microbes to produce desirable compounds through genetic engineering or synthetic biology approaches. Metabolomics involves the quantitation of intracellular and extracellular metabolites, where mass spectrometry and nuclear magnetic resonance based analytical instrumentation are often used. Here, the experimental designs, sample preparations, metabolite quenching and extraction are essential to the quantitative metabolomics workflow. The resultant metabolomics data can then be used with computational modelling approaches, such as kinetic and constraint-based modelling, to better understand underlying mechanisms and bottlenecks in the synthesis of desired compounds, thereby accelerating research through systems metabolic engineering. Constraint-based models, such as genome scale models, have been used successfully to enhance the yield of desired compounds from engineered microbes, however, unlike kinetic or dynamic models, constraint-based models do not incorporate regulatory effects. Nevertheless, the lack of time-series metabolomic data generation has hindered the usefulness of dynamic models till today. In this review, we show that improvements in automation, dynamic real-time analysis and high throughput workflows can drive the generation of more quality data for dynamic models through time-series metabolomics data generation. Spatial metabolomics also has the potential to be used as a complementary approach to conventional metabolomics, as it provides information on the localization of metabolites. However, more effort must be undertaken to identify metabolites from spatial metabolomics data derived through imaging mass spectrometry, where machine learning approaches could prove useful. On the other hand, single-cell metabolomics has also seen rapid growth, where understanding cell-cell heterogeneity can provide more insights into efficient metabolic engineering of microbes. Moving forward, with potential improvements in automation, dynamic real-time analysis, high throughput workflows, and spatial metabolomics, more data can be produced and studied using machine learning algorithms, in conjunction with dynamic models, to generate qualitative and quantitative predictions to advance metabolic engineering efforts.

Transient rapamycin treatment during developmental stage extends lifespan in Mus musculus and Drosophila melanogaster

Lifespan is determined by complex and tangled mechanisms that are largely unknown. The early postnatal stage has been proposed to play a role in lifespan, but its contribution is still controversial. Here, we show that a short rapamycin treatment during early life can prolong lifespan in Mus musculus and Drosophila melanogaster. Notably, the same treatment at later time points has no effect on lifespan, suggesting that a specific time window is involved in lifespan regulation. We also find that sulfotransferases are upregulated during early rapamycin treatment both in newborn mice and in Drosophila larvae, and transient dST1 overexpression in Drosophila larvae extends lifespan. Our findings unveil a novel link between early‐life treatments and long‐term effects on lifespan.

Effects of dietary branched-chain amino acid supplementation on serum and milk metabolome profiles in dairy cows during early lactation

The 3 branched-chain AA (BCAA), Val, Leu, and Ile, are essential AA used by tissues as substrates for protein synthesis and energy generation. In addition, BCAA are also involved in modulating cell signaling pathways, such as nutrient sensing and insulin signaling. In our previous study, dietary BCAA supplementation was shown to improve protein synthesis and glucose homeostasis in transition cows. However, a more detailed understanding of the changes in metabolic pathways associated with an increased BCAA availability is desired to fine-tune nutritional supplementation strategies. Multiparous Holstein cows (n = 20) were enrolled 28 d before expected calving and assigned to either the BCAA treatment (n = 10) or the control group (n = 10). Cows assigned to BCAA were fed 550 g/d of rumen-protected BCAA mixed with 200 g/d of dry molasses from calving until 35 DIM, whereas the cows assigned to the control were fed only 200 g/d of dry molasses. Serum samples were collected on d 10 before expected calving, as well as on d 4 and d 21 postpartum. Milk samples were collected on d 14 postpartum. From a larger cohort, we selected 20 BCAA-supplemented cows with the greatest plasma urea nitrogen concentration, as an indicator for greater BCAA availability, for the metabolomics analysis herein. Serum and milk samples were subjected to a liquid chromatography-mass spectrometry-based assay, detecting and measuring the abundance of 241 serum and 211 milk metabolic features, respectively. Multivariable statistical analyses revealed that BCAA supplementation altered the metabolome profiles of both serum and milk samples. Increased abundance of serum phosphocholine and glutathione and of milk Val, Ile, and Leu, and decreased abundance of milk acyl-carnitines were associated with BCAA supplementation. Altered phosphocholine and glutathione abundances point to altered hepatic choline metabolism and antioxidant balance, respectively. Altered milk acyl-carnitine abundances suggest changes in mammary fatty acid metabolism. Dietary BCAA supplementation was associated with a range of alterations in serum and milk metabolome profiles, adding to our understanding of the role of BCAA availability in modulating dairy cow protein, lipid, and energy metabolism on a whole-body level and how it affects milk composition.

Rewiring of Glucose and Lipid Metabolism Induced by G Protein-Coupled Receptor 17 Silencing Enables the Transition of Oligodendrocyte Progenitors to Myelinating Cells

In the mature central nervous system (CNS), oligodendrocytes (OLs) provide support and insulation to axons thanks to the production of a myelin sheath. During their maturation to myelinating cells, OLs require energy and building blocks for lipids, which implies a great investment of energy fuels and molecular sources of carbon. The oligodendroglial G protein-coupled receptor 17 (GPR17) has emerged as a key player in OL maturation; it reaches maximal expression in pre-OLs, but then it has to be internalized to allow terminal maturation. In this study, we aim at elucidating the role of physiological GPR17 downregulation in OL metabolism by applying transcriptomics, metabolomics and lipidomics on differentiating OLs. After GPR17 silencing, we found a significant increase in mature OL markers and alteration of several genes involved in glucose metabolism and lipid biosynthesis. We also observed an increased release of lactate, which is partially responsible for the maturation boost induced by GPR17 downregulation. Concomitantly, GPR17 depletion also changed the kinetics of specific myelin lipid classes. Globally, this study unveils a functional link between GPR17 expression, lactate release and myelin composition, and suggests that innovative interventions targeting GPR17 may help to foster endogenous myelination in demyelinating diseases.

Changes in the lipidome of water buffalo milk during intramammary infection by non-aureus Staphylococci

This study aimed to determine the lipidome of water buffalo milk with intramammary infection (IMI) by non-aureus staphylococci (NAS), also defined as coagulase-negative staphylococci, using an untargeted lipidomic approach. Non-aureus Staphylococci are the most frequently isolated pathogens from dairy water buffalo milk during mastitis. A total of 17 milk samples from quarters affected by NAS-IMI were collected, and the lipidome was determined by liquid chromatography-quadrupole time-of-flight mass spectrometry. The results were compared with the lipidome determined on samples collected from 16 healthy quarters. The study identified 1934 different lipids, which were classified into 15 classes. The abundance of 72 lipids changed in NAS-IMI milk compared to healthy quarters. Significant changes occurred primarily in the class of free fatty acids. The results of this study provided first-time insight into the lipidome of dairy water buffalo milk. Moreover, the present findings provide evidence that NAS-IMI induces changes in water buffalo milk's lipidome.

Identification and characterization of dairy cows with different backfat thickness antepartum in relation to postpartum loss of backfat thickness: A cluster analytic approach

The objectives of this study were (1) to characterize the interindividual variation in the relationship between antepartum (ap) backfat thickness (BFT) and subsequent BFT loss during early lactation in a large dairy herd using cluster analysis; (2) to compare the serum concentrations of metabolites (nonesterified fatty acids, β-hydroxybutyrate), metabolic hormones (leptin and adiponectin), and an inflammatory marker (haptoglobin) among the respective clusters; and (3) to compare lactation performance and uterine health status in the different clusters. An additional objective was (4) to investigate differences in these serum variables and in milk yield of overconditioned (OC) cows that differed in the extent of BFT loss. Using data from a large study of 1,709 multiparous Holstein cows, we first selected those animals from which serum samples and BFT results (mm) were available at d 25 (±10) ap and d 31 (±3 d) postpartum (pp). The remaining 713 cows (parity of 2 to 7) were then subjected to cluster analysis: different approaches based on the BFT of the cows were performed. K-means (unsupervised machine learning algorithm) clustering based on BFT-ap alone identified 5 clusters: lean (5-8 mm BFT, n = 50), normal (9-12 mm, n = 206), slightly fat (SF; 13-16 mm, n = 203), just fat (JF; 16-22 mm, n = 193), and very fat (VF; 23-43 mm, n = 61). Clustering by difference between BFT-ap and BFT-pp (ΔBFT) also revealed 5 clusters: extreme loss (17-23 mm ΔBFT, n = 16), moderate loss (9-15 mm, n = 119), little loss (4-8 mm, n = 326), no loss (0-3 mm, n = 203), and gain (-8 to -1 mm, n = 51). Based on the blood variables measured, our results confirm that cows with greater BFT losses had higher lipid mobilization and ketogenesis than cows with less BFT loss. The serum variables of cows that gained BFT did not differ from normal cows. Milk yield was affected by the BFT-ap cluster, but not by the ΔBFT cluster. Cows categorized as VF had lesser milk yield than other clusters. We further compared the OC cows that had little or no BFT loss (i.e., 2% of VF, 12% of JF, and 31% of SF, OC-no loss, n = 85) with the OC cows that lost BFT (OC-loss, n = 135). Both NEFA and BHB pp concentrations and milk yield were greater in OC-loss cows compared with the OC-no loss cows. The serum concentration of leptin ap was greater in OC-loss than in the OC-no loss cows. Overall, OC cows lost more BFT than normal or lean cows. However, those OC cows with a smaller loss of BFT produced less milk than OC cows with greater losses.

Metabolomic signature and mitochondrial dynamics outline the difference between vulnerability and resilience to chronic stress

Stress is the foremost environmental factor involved in the pathophysiology of major depressive disorder (MDD). However, individual differences among people are critical as some people exhibit vulnerability while other are resilient to repeated exposure to stress. Among the others, a recent theory postulates that alterations of energy metabolism might contribute to the development of psychopathologies. Here we show that the bioenergetic status in the ventral hippocampus (vHip), a brain subregion tightly involved in the regulation of MDD, defined the development of vulnerability or resilience following two weeks of chronic mild stress. Among the different metabolomic signatures observed, the glycolysis and tricarboxylic acid cycle may be specifically involved in defining vulnerability, revealing a previously unappreciated mechanism of sensitivity to stress. These findings point to mitochondrial morphology and recycling as critical in the ability to cope with stress. We show that vulnerable rats favor mitochondrial fusion to counteract the overproduction of reactive oxidative species whereas resilient rats activate fission to guarantee metabolic efficiency. Our results indicate that the modulation of the energetic metabolite profile in vHip under chronic stress exposure may represent a mechanism to explain the difference between vulnerable and resilient rats, unraveling novel and promising targets for specific therapeutic interventions.

Oxidative pentose phosphate pathway controls vascular mural cell coverage by regulating extracellular matrix composition

Vascular mural cells (vMCs) play an essential role in the development and maturation of the vasculature by promoting vessel stabilization through their interactions with endothelial cells. Whether endothelial metabolism influences mural cell recruitment and differentiation is unknown. Here, we show that the oxidative pentose phosphate pathway (oxPPP) in endothelial cells is required for establishing vMC coverage of the dorsal aorta during early vertebrate development in zebrafish and mice. We demonstrate that laminar shear stress and blood flow maintain oxPPP activity, which in turn, promotes elastin expression in blood vessels through production of ribose-5-phosphate. Elastin is both necessary and sufficient to drive vMC recruitment and maintenance when the oxPPP is active. In summary, our work demonstrates that endothelial cell metabolism regulates blood vessel maturation by controlling vascular matrix composition and vMC recruitment.

Over 20 Years of Machine Learning Applications on Dairy Farms: A Comprehensive Mapping Study

Machine learning applications are becoming more ubiquitous in dairy farming decision support applications in areas such as feeding, animal husbandry, healthcare, animal behavior, milking and resource management. Thus, the objective of this mapping study was to collate and assess studies published in journals and conference proceedings between 1999 and 2021, which applied machine learning algorithms to dairy farming-related problems to identify trends in the geographical origins of data, as well as the algorithms, features and evaluation metrics and methods used. This mapping study was carried out in line with PRISMA guidelines, with six pre-defined research questions (RQ) and a broad and unbiased search strategy that explored five databases. In total, 129 publications passed the pre-defined selection criteria, from which relevant data required to answer each RQ were extracted and analyzed. This study found that Europe (43% of studies) produced the largest number of publications (RQ1), while the largest number of articles were published in the Computers and Electronics in Agriculture journal (21%) (RQ2). The largest number of studies addressed problems related to the physiology and health of dairy cows (32%) (RQ3), while the most frequently employed feature data were derived from sensors (48%) (RQ4). The largest number of studies employed tree-based algorithms (54%) (RQ5), while RMSE (56%) (regression) and accuracy (77%) (classification) were the most frequently employed metrics used, and hold-out cross-validation (39%) was the most frequently employed evaluation method (RQ6). Since 2018, there has been more than a sevenfold increase in the number of studies that focused on the physiology and health of dairy cows, compared to almost a threefold increase in the overall number of publications, suggesting an increased focus on this subdomain. In addition, a fivefold increase in the number of publications that employed neural network algorithms was identified since 2018, in comparison to a threefold increase in the use of both tree-based algorithms and statistical regression algorithms, suggesting an increasing utilization of neural network-based algorithms.

Dietary Concentrate Supplementation Alters Serum Metabolic Profiles Related to Energy and Amino Acid Metabolism in Grazing Simmental Heifers

Supplementation plays a vital role in the growth performance of grazing heifers. We investigated the effects of maize-based concentrate supplementation on the serum metabolome in grazing heifers. Twenty-four 7-month-old heifers (211.65 ± 4.25 kg BW) were randomly divided into a supplement (SUP) group and a control (CON) group. The results indicated that concentrate supplementation increased the final body weight (BW) of grazing heifers, and the average daily gain (ADG) was 61.5% (P = 0.011) higher in the SUP group than in the CON group. Serum concentrations of total protein (TP), triglyceride (TG), and leptin were higher in the SUP group than in the CON group (p < 0.05). Supplementation increased serum metabolites and amino acids and markedly altered glucose, lipid, and protein metabolism, which contributed to the heifer growth. Furthermore, by multivariate analysis, 45 serum metabolites were identified as significantly different between the two groups. Enrichment analysis revealed that arginine biosynthesis and tryptophan metabolism as well as glycerophospholipid metabolism were significantly enriched between the two groups. We concluded that the growth potential of heifers could be improved by maize-based concentrate supplementation, and the main biological pathways affected were those related to energy and amino acid metabolism.

Overview of Metabolomic Analysis and the Integration with Multi-Omics for Economic Traits in Cattle

Metabolomics has been applied to measure the dynamic metabolic responses, to understand the systematic biological networks, to reveal the potential genetic architecture, etc., for human diseases and livestock traits. For example, the current published results include the detected relevant candidate metabolites, identified metabolic pathways, potential systematic networks, etc., for different cattle traits that can be applied for further metabolomic and integrated omics studies. Therefore, summarizing the applications of metabolomics for economic traits is required in cattle. We here provide a comprehensive review about metabolomic analysis and its integration with other omics in five aspects: (1) characterization of the metabolomic profile of cattle; (2) metabolomic applications in cattle; (3) integrated metabolomic analysis with other omics; (4) methods and tools in metabolomic analysis; and (5) further potentialities. The review aims to investigate the existing metabolomic studies by highlighting the results in cattle, integrated with other omics studies, to understand the metabolic mechanisms underlying the economic traits and to provide useful information for further research and practical breeding programs in cattle.

Branched-Chain Amino Acid Supplementation Alters the Abundance of Mechanistic Target of Rapamycin and Insulin Signaling Proteins in Subcutaneous Adipose Explants from Lactating Holstein Cows

Simple Summary

Branched-chain amino acids (BCAAs) are import regulators of mechanistic target of rapamycin (mTOR). In humans and rodents, increased circulating BCAA levels are positively associated with changes in protein abundance of insulin and amino acid (AA) signaling pathways in organs such as skeletal muscle and adipose. Unlike aspects of fatty acid metabolism (e.g., lipolysis, lipogenesis), it is unknown if BCAA directly affect subcutaneous adipose tissue (SAT) AA metabolism and insulin signaling. We propose that BCAA availability within SAT could enhance aspects of AA and insulin function by promoting increases in the abundance of key proteins.

Abstract

The objective of this study was to investigate changes in protein abundance of mTOR and insulin signaling pathway components along with amino acid (AA) transporters in bovine s.c. adipose (SAT) explants in response to increased supply of Leu, Ile, or Val. Explants of SAT from four lactating Holstein cows were incubated with high-glucose serum-free DMEM, to which the 10 essential AAs were added to create the following treatments: ideal mix of essential AA (IPAA; Lys:Met 2.9:1; Lys:Thr 1.8:1; Lys:His 2.38:1; Lys:Val 1.23:1; Lys:Ile 1.45:1; Lys:Leu 0.85:1; Lys:Arg 2.08:1) or IPAA supplemented with Ile, Val, or Leu to achieve a Lys:Ile of 1.29:1 (incIle), Lys:Val 1.12:1 (incVal), or Lys:Leu (incLeu) 0.78:1 for 4 h. Compared with IPAA, incLeu or incIle led to greater activation of protein kinase B (AKT; p-AKT/total AKT) and mTOR (p-mTOR/total mTOR). Total EAA in media averaged 7.8 ± 0.06 mmol/L across treatments. Incubation with incLeu, incIle, or incVal led to greater protein abundance of solute carrier family 38 member 1 (SLC38A1), a Gln transporter, and the BCAA catabolism enzyme branched-chain α-keto acid dehydrogenase kinase (BCKDK) compared with IPAA. Activation of eukaryotic elongation factor 2 (eEF2; p-eEF2/total eEF2) was also greater in response to incLeu, incIle, or incVal. Furthermore, compared with incLeu or incIle, incVal supplementation led to greater abundance of SLC38A1 and BCKDK. BCKDK is a rate-limiting enzyme regulating BCAA catabolism via inactivation and phosphorylation of the BCKD complex. Overall, data suggested that enhanced individual supplementation of BCAA activates mTOR and insulin signaling in SAT. Increased AA transport into tissue and lower BCAA catabolism could be part of the mechanism driving these responses. The potential practical applications for enhancing post-ruminal supply of BCAA via feeding in rumen-protected form support in vivo studies to ascertain the role of these AAs on adipose tissue biology.

Association between alterations in plasma metabolome profiles and laminitis in intensively finished Holstein bulls in a randomized controlled study

Metabolic consequences of an energy and protein rich diet can compromise metabolic health of cattle by promoting a pro-inflammatory phenotype. Laminitis is a common clinical sign, but affected metabolic pathways, underlying pathophysiology and causative relationships of a systemic pro-inflammatory phenotype are unclear. Therefore, the aim of this study was to elucidate changes in metabolome profiles of 20 months old Holstein bulls fed a high energy and protein diet and to identify novel metabolites and affected pathways, associated with diet-related laminitis. In a randomized controlled feeding trial using bulls fed a high energy and protein diet (HEP; metabolizable energy [ME] intake 169.0 ± 1.4 MJ/day; crude protein [CP] intake 2.3 ± 0.02 kg/day; calculated means ± SEM; n = 15) versus a low energy and protein diet (LEP; ME intake 92.9 ± 1.3 MJ/day; CP intake 1.0 ± 0.01 kg/day; n = 15), wide ranging effects of HEP diet on metabolism were demonstrated with a targeted metabolomics approach using the AbsoluteIDQ p180 kit (Biocrates Life Sciences). Multivariate statistics revealed that lower concentrations of phosphatidylcholines and sphingomyelins and higher concentrations of lyso-phosphatidylcholines, branched chain amino acids and aromatic amino acids were associated with an inflammatory state of diet-related laminitis in Holstein bulls fed a HEP diet. The latter two metabolites share similarities with changes in metabolism of obese humans, indicating a conserved pathophysiological role. The observed alterations in the metabolome provide further explanation on the underlying metabolic consequences of excessive dietary nutrient intake.

Longitudinal changes in fatty acid metabolism and in the mitochondrial protein import system in overconditioned and normal conditioned cows: A transcriptional study using microfluidic quantitative PCR

This study investigated the effect of body condition around calving on the hepatic mRNA expression of genes involved in fatty acid (FA) metabolism and mitochondrial protein import system of dairy cows during the transition period. Fifteen weeks before their anticipated calving date, 38 multiparous Holstein cows were selected based on their current and previous body condition scores (BCS) and allocated to either a high or a normal BCS group (19 cows each). They received different diets to reach targeted differences in BCS and backfat thickness (BFT) until dry-off. At dry-off, normal BCS (NBCS) cows had a BCS <3.5 and BFT <1.2 cm, and the high BCS (HBCS) cows had a BCS >3.75 and BFT >1.4 cm. The expression of targeted genes in the liver was assayed by reverse-transcription quantitative real-time PCR using microfluidics integrated fluidic circuit chips on a subset of 5 cows from each group. Liver biopsies were collected at d -49, +3, +21, and +84 relative to parturition. The mRNA abundance of 47 genes related to lipid metabolism including carnitine metabolism, FA uptake and transport, lipoprotein export, carnitine metabolism, mitochondrial and proximal FA oxidation, ketogenesis, AMP-activated protein kinase/mammalian target of rapamycin pathway, and mitochondrial protein import system was assessed in liver tissue. The mRNA abundances of FA binding protein (FABP)6 (in both groups), and FABP1 and solute carrier family 22 member 5 (SLC22A5) in HBCS were upregulated (>1.5-fold change, FC) in early lactation (at d +3 and +21 postpartum) compared with antepartum (d -49), indicating promoted FA uptake and intracellular transport in the liver due to the metabolic adaptations of elevated lipo-mobilization after parturition. The upregulation of SLC22A5 and SLC25A20 after parturition was more pronounced in HBCS than in NBCS cows, suggesting a need for increasing the capacity of FA uptake, and FA transport into the hepatocyte. The increased mRNA abundance of carnitine palmitoyltransferase 1A, after parturition and to a greater extent in HBCS (FC = 4.1) versus NBCS (FC = 2.1) indicates a physiological increase in the capacity of long-chain fatty acyl-CoA entry into the liver mitochondria compared with antepartum (ap; d -49 relative to calving). The greater hepatic mRNA abundance of genes encoding enzymes involved in mitochondrial FA oxidation in HBCS than in NBCS points to an increased rate of mitochondrial β-oxidation. The hepatic mRNA abundance of 3-hydroxy-3-methylglutaryl-CoA synthase 2 and 3-hydroxy-3-methylglutaryl-CoA were upregulated after parturition (d +21/d +3 pp) to a greater extent in HBCS than in NBCS cows, indicating that excess acetyl-CoA generated via β-oxidation was increasingly used for ketogenesis. We observed for the first time that the mRNA abundance of genes involved in the translocase of the inner membrane (TIM) complex (TIM22 and TIM23) in the hepatic mitochondrial protein import system were undergoing distinct changes during the transition from late pregnancy to early lactation in dairy cows. Even though sample size in this study was relatively small, the results support that overconditioning around calving may contribute to mitochondrial FA overload and greater ketogenesis at the level of transcription in the liver of early lactation cows.

The untargeted lipidomic profile of quarter milk from dairy cows with subclinical intramammary infection by non-aureus staphylococci

This observational study determined the lipidome of cow milk during subclinical intramammary infection (IMI) by non-aureus staphylococci (NAS), also defined as coagulase-negative staphylococci, using an untargeted approach. Among the pathogens causing bovine IMI, NAS have become the most frequently isolated bacteria from milk samples. Although the application of system biology approaches to mastitis has provided pivotal information by investigating the transcriptome, proteome, peptidome, and metabolome, the milk lipidome during mammary gland inflammation remains undisclosed. To cover this gap, we determined the milk lipidome of 17 dairy cows with IMI caused by NAS (NAS-IMI), and we compared the results with those of healthy quarter milk from 11 cows. The lipidome was determined following a liquid chromatography-quadrupole time-of-flight mass spectrometry approach. Sixteen subclasses of lipids were identified in both groups of animals. From 2,556 measured lipids, the abundance of 597 changed more than 10-fold in quarter milk with NAS-IMI compared with healthy quarters. The results demonstrate the influence of NAS-IMI on the milk lipidome, implying significant changes in lipid species belonging to the family of triacylglycerols and sphingomyelins, and contribute to the understanding of inflammatory processes in the bovine udder, highlighting potential novel biomarkers for improving mastitis diagnostics.

Aspects of transition cow metabolomics-Part I: Effects of a metaphylactic butaphosphan and cyanocobalamin treatment on the metabolome in liver, blood, and urine in cows with different liver metabotypes

Dairy cows in modern production systems are at risk to develop metabolic disorders during the transition period. Reasons for individual differences in susceptibility, as well as the underlying pathomechanisms, are still only partially understood. The development of metaphylactic treatment protocols is needed. In this context, an on-farm prospective 3-fold blinded randomized study involving 80 German Holstein cows was performed throughout 1 yr. The trial involved a thorough recording of the production and clinical traits, clinical chemistry, and liver biopsies and blood and urine sampling at d 14 (mean: 12 d, range: 1-26 d) antepartum (AP), and d 7 (7, 4-13) and 28 (28, 23-34) postpartum (PP) for metabolomics analyses. Two groups received a treatment with butaphosphan and cyanocobalamin (BCC) at either the dosage recommended by the manufacturer or the double dosage (5 or 10 mL/100 kg of body weight 10% butaphosphan and 0.005% cyanocobalamin (Catosal, Bayer Animal Health), n = 20 in each group, parity: 4.2 ± 2.0 and 3.4 ± 1.3, respectively (mean ± SD)] and one group a placebo treatment (NaCl 0.9%, n = 40, parity: 4.0 ± 1.9). The animals were treated at 6 time points (7, 6, and 5 d AP, and 1, 2, and 3 d PP) via intravenous injection. Mass spectroscopy-based targeted metabolomics analysis of blood plasma and liver samples were performed using the AbsoluteIDQ p180 kit (Biocrates Life Sciences), whereas the urine samples were analyzed by nuclear magnetic resonance spectroscopy. Statistical analysis was performed using multivariate [partial least squares discriminant analysis (PLS-DA)] and univariate methods (linear mixed model). Multivariate data analysis (PLS-DA plots) of the liver metabolome revealed 3 different metabotypes (A = medium, B = minor, C = large alterations in liver metabolome profile between AP and PP status). Metabotype B animals were characterized by higher PP lipomobilization (stronger PP body condition decrease and higher blood bilirubin, fatty acids, gamma-glutamyltransferase, and triglyceride levels) and a higher occurrence of transition cow diseases, compared with the animals in metabotype C. Analysis of the feeding data showed that the period of metabotype B animals (calving in a distinct time frame) was characterized by a decreased grass silage quality. The PP liver metabolome of the metabotype C animals was characterized by higher concentrations of AA, acylcarnitines, lysoPC and sphingomyelins compared with metabotype B. For the metaphylactic treatment with BCC a dose-dependent effect was confirmed, differing between the metabotypes. In all matrices and metabotypes at various time points significant treatment effects were observed, with different profiles in clinical chemistry and as well in metabolomics data. The most clear-cut treatment effect was observed in metabotype B in the liver at 7 d PP, characterized by an increase in several acylcarnitines and phosphatidylcholines, indicating a more efficient influx and oxidation of fatty acids in mitochondria and thereby an increase in energy supply and more efficient triglyceride export in the liver. The results from the liver metabolomics analysis support the application of an indication-based metaphylactic treatment with BCC.

Aspects of transition cow metabolomics-Part III: Alterations in the metabolome of liver and blood throughout the transition period in cows with different liver metabotypes

The liver plays a central role in the postpartum (PP) energy metabolism of the transition dairy cow; however, studies describing the liver metabolome during this period were lacking. The aim of the presented study was therefore to compare the alterations in the liver and blood metabolome of transition dairy cows. For this purpose, an on-farm trial with 80 German Holstein cows (mean lactation number: 3.9; range: 2-9) was performed, with thorough documentation of clinical traits and clinical chemistry, as well as production data. Liver biopsies and blood samples were collected at d 14 (mean: 12 d, range: 1-26 d) antepartum (AP), d 7 (7, 4-13) and 28 (28, 23-34; mean, earliest-latest) PP for targeted mass spectroscopy-based metabolomics analysis using the AbsoluteIDQ p180 kit (Biocrates Life Sciences). Statistical analysis was performed using multivariate (partial least squares discriminant analysis) as well as univariate methods (linear mixed model). Multivariate data analysis of the liver metabolome revealed 3 different metabotypes (A = medium, B = minor, C = large alterations in the liver metabolome profile between AP and PP). In metabotype C, an increase of almost all acylcarnitines, lysophosphatidylcholines (lysoPC), sphingomyelins, and some phosphatidylcholines (PC, mainly at 7 d PP) was observed after calving. In contrast to metabotype C, the clinical data of the metabotype B animals indicated a higher PP lipomobilization and occurrence of transition cow diseases. The liver metabolome profile of these animals most likely mirrors a failure of adaptation to the PP state. This strong occurrence of metabotypes was much less pronounced in the blood metabolome. Additionally, differences in metabolic patterns were observed across the transition period when comparing liver and blood matrices (e.g., in different biogenic amines, acylcarnitines and sphingolipids). In summary, the blood samples at 7 d PP showed lower acylcarnitines and PC, with minor alterations and a heterogeneous pattern in AA, biogenic amines, and sphingomyelins compared with 14 d AP. In contrast to 7 d PP, the blood samples at 28 PP revealed an increase in several AA, lysoPC, PC, and sphingomyelins in comparison to the AP state, irrespective of the metabotype. In the liver biopsies metabotype B differed from metabotype C animals ante partum by following metabolites: higher α aminoadipic acid, lower AA, serotonin, taurine, and symmetric dimethylarginine levels, lower or higher concentrations of certain acylcarnitines (higher: C2, C3, C5, C4:1; lower: C12:1, C14:1-OH, C16:2), and lower lysoPC (a C16:0, C18:0, C20:3, C20:4) and hexose levels. In blood samples, fewer differences were observed, with lower serotonin, acylcarnitine C16:2, lysoPC (a C16:0, C17:0, C18:0 and C18:1), PC aa C38:0, and PC ae C42:2. The results show that the use of only the blood metabolome to assess liver metabolism may be hampered by the fact that blood profiles are influenced by the metabolism of many organs, and metabolomics analysis from liver biopsies is a more suitable method to identify distinct metabotypes. Future studies should investigate the stability and reproducibility of the metabotype and phenotypes observed, and the possible predictive value of the metabolites already differing AP between metabotype B and C.

Body condition loss during the dry period: Insights from feeding behavior studies

Body condition change during the dry period (ΔBCS) has been associated with a myriad of transition cow diseases. We used data from 3 studies to assess the relationship between ΔBCS, feeding behavior, and body condition score (BCS) at dry-off. We also studied the mediation effect that dry matter intake (DMI) has on ΔBCS and the association between dry period feeding time and DMI. A total of 100 parous cows were enrolled in 3 studies to investigate differences in dry period diet on behavior, health, and performance pre- and postcalving. Body weight was measured and BCS was assessed by the same trained observer after dry-off and 1 wk from calving date. The ΔBCS was calculated by subtracting the BCS at calving minus the BCS at dry-off. The BCS at dry-off was categorized as overconditioned (≥3.5) or not overconditioned (<3.5); no cows had a BCS <2.75. Feeding behavior data were collected using electronic feed bins. Parity at dry-off (median = 2; min = 1, max = 6) and 305-d milk production (mean = 10,235 kg, SD = 1,625 kg) from the previous lactation were considered. Data sets were split into 2 time periods: d -56 to -22 (early) and -21 to 0 (late) in relation to calving. Selected feeding behaviors (DMI, DMI as a percentage of body weight, and feeding time) were used to evaluate the associations between each feeding behavior and BCS at dry-off in each period using mixed linear regression models. Each model included the following covariates: parity, previous 305-d milk yield, and trial treatment. Experimental day was included as random slope, and cow was included as random intercept. A mediation analysis was used to evaluate the potential causal direct effect of BCS at dry-off on ΔBCS and the potential indirect effect mediated by differences in DMI. The BCS at dry-off was associated with changes in feeding behavior, such that overconditioned cows had lesser daily DMI and feeding time during the early and late dry periods compared with not overconditioned animals. We also noted an effect of previous 305-d milk yield on DMI; cows that produced more milk had greater DMI throughout the dry period. The ΔBCS was only partially mediated by DMI, and BCS at dry-off still had a direct effect on ΔBCS. This result indicated that mechanisms other than DMI were associated with BCS loss during the dry period. Feeding time correlated weakly and moderately with DMI during the early and late dry periods, respectively. To conclude, strategies to improve intake during the dry period should take dry-off BCS into account or, preferably, efforts should be made to minimize the number of overconditioned cows at the end of lactation.

Systems biology approaches integrated with artificial intelligence for optimized metabolic engineering

Metabolic engineering aims to maximize the production of bio-economically important substances (compounds, enzymes, or other proteins) through the optimization of the genetics, cellular processes and growth conditions of microorganisms. This requires detailed understanding of underlying metabolic pathways involved in the production of the targeted substances, and how the cellular processes or growth conditions are regulated by the engineering. To achieve this goal, a large system of experimental techniques, compound libraries, computational methods and data resources, including multi-omics data, are used. The recent advent of multi-omics systems biology approaches significantly impacted the field by opening new avenues to perform dynamic and large-scale analyses that deepen our knowledge on the manipulations. However, with the enormous transcriptomics, proteomics and metabolomics available, it is a daunting task to integrate the data for a more holistic understanding. Novel data mining and analytics approaches, including Artificial Intelligence (AI), can provide breakthroughs where traditional low-throughput experiment-alone methods cannot easily achieve. Here, we review the latest attempts of combining systems biology and AI in metabolic engineering research, and highlight how this alliance can help overcome the current challenges facing industrial biotechnology, especially for food-related substances and compounds using microorganisms.

Machine learning approach to integrated endometrial transcriptomic datasets reveals biomarkers predicting uterine receptivity in cattle at seven days after estrous

The main goal was to apply machine learning (ML) methods on integrated multi-transcriptomic data, to identify endometrial genes capable of predicting uterine receptivity according to their expression patterns in the cow. Public data from five studies were re-analyzed. In all of them, endometrial samples were obtained at day 6–7 of the estrous cycle, from cows or heifers of four different European breeds, classified as pregnant (n = 26) or not (n = 26). First, gene selection was performed through supervised and unsupervised ML algorithms. Then, the predictive ability of potential key genes was evaluated through support vector machine as classifier, using the expression levels of the samples from all the breeds but one, to train the model, and the samples from that one breed, to test it. Finally, the biological meaning of the key genes was explored. Fifty genes were identified, and they could predict uterine receptivity with an overall 96.1% accuracy, despite the animal’s breed and category. Genes with higher expression in the pregnant cows were related to circadian rhythm, Wnt receptor signaling pathway, and embryonic development. This novel and robust combination of computational tools allowed the identification of a group of biologically relevant endometrial genes that could support pregnancy in the cattle.

Metabolic Profiling of Plasma in Different Calving Body Condition Score Cows Using an Untargeted Liquid Chromatography-Mass Spectrometry Metabolomics Approach

This study was undertaken to identify metabolite differences in plasma of dairy cows with a normal or high calving body condition score (CBCS), using untargeted liquid chromatography-mass spectrometry (LC-MS) metabolomics. Sixteen multiparous dairy cows were assigned to one of two groups based on CBCS (0 to 5 scale): Normal group (NBCS, 3.25 ≤ BCS ≤ 3.5, n = 8), and high BCS group (HBCS, BCS ≥ 4, n = 8). Plasma samples were collected for metabolomics analysis and evaluation of biomarkers of lipid metabolism (nonesterified fatty acid (NEFA) and β-hydroxybutyrate (BHB)), and cytokines (leptin, adiponectin, tumor necrosis factor-α (TNF-α) and interleukin 6 (IL-6)). A total of 23 differential metabolites were identified, and functional analyses were performed using the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Among these metabolites, the concentrations of six lysophosphatidylcholines and one phosphatidylethanolamine, were lower in the HBCS group than in the NBCS group (p < 0.01). Furthermore, these metabolites were involved in these four pathways, among others: glycerophospholipid metabolism, retrograde endocannabinoid signaling, autophagy, and glycosylphosphatidylinositol (GPI)-anchor biosynthesis (p < 0.05). In addition, plasma concentrations of leptin (p = 0.06) and TNF-α (p = 0.08) tended to be greater while adiponectin (p = 0.09) lower in HBCS cows than in NBCS cows. The concentrations of NEFA, BHB, or IL-6 did not differ between NBCS and HBCS groups. More importantly, based on the results of the Spearman's correlation analysis, the seven important metabolites were negatively correlated with indices of lipid metabolisms, proinflammatory cytokines, and leptin, but positively correlated with adiponectin. These results demonstrate that CBCS has a measurable impact on the plasma metabolic profile, even when NEFA and BHB are not different. In addition, the identified differential metabolites were significantly correlated to lipid metabolism and inflammation in the over-conditioned fresh cows, which are expected to render a metabolic basis for the diseases associated with over-conditioned dry cows.

Deep learning meets metabolomics: a methodological perspective

Deep learning (DL), an emerging area of investigation in the fields of machine learning and artificial intelligence, has markedly advanced over the past years. DL techniques are being applied to assist medical professionals and researchers in improving clinical diagnosis, disease prediction and drug discovery. It is expected that DL will help to provide actionable knowledge from a variety of 'big data', including metabolomics data. In this review, we discuss the applicability of DL to metabolomics, while presenting and discussing several examples from recent research. We emphasize the use of DL in tackling bottlenecks in metabolomics data acquisition, processing, metabolite identification, as well as in metabolic phenotyping and biomarker discovery. Finally, we discuss how DL is used in genome-scale metabolic modelling and in interpretation of metabolomics data. The DL-based approaches discussed here may assist computational biologists with the integration, prediction and drawing of statistical inference about biological outcomes, based on metabolomics data.

Molecular networks of insulin signaling and amino acid metabolism in subcutaneous adipose tissue are altered by body condition in periparturient Holstein cows

Peripartal cows mobilize not only body fat but also body protein to satisfy their energy requirements. The objective of this study was to determine the effect of prepartum BCS on blood biomarkers related to energy and nitrogen metabolism, and mRNA and protein abundance associated with AA metabolism and insulin signaling in subcutaneous adipose tissue (SAT) in peripartal cows. Twenty-two multiparous Holstein cows were retrospectively classified into a high BCS (HBCS; n = 11, BCS ≥ 3.5) or normal BCS (NBCS; n = 11, BCS ≤ 3.17) group at d 28 before expected parturition. Cows were fed the same diet as a total mixed ration before parturition and were fed the same lactation diet postpartum. Blood samples collected at -10, 7, 15, and 30 d relative to parturition were used for analyses of biomarkers associated with energy and nitrogen metabolism. Biopsies of SAT harvested at -15, 7, and 30 d relative to parturition were used for mRNA (real time-PCR) and protein abundance (Western blotting) assays. Data were subjected to ANOVA using the MIXED procedure of SAS (v. 9.4; SAS Institute Inc., Cary, NC), with P ≤ 0.05 being the threshold for significance. Cows in HBCS had greater overall plasma nonesterified fatty acid concentrations, due to marked increases at 7 and 15 d postpartum. This response was similar (BCS × Day effect) to protein abundance of phosphorylated (p) protein kinase B (p-AKT), the insulin-induced glucose transporter (SLC2A4), and the sodium-coupled neutral AA transporter (SLC38A1). Abundance of these proteins was lower at -15 d compared with NBCS cows, and either increased (SLC2A4, SLC38A1) or did not change (p-AKT) at 7 d postpartum in HBCS. Unlike protein abundance, however, overall mRNA abundances of the high-affinity cationic (SLC7A1), proton-coupled (SLC36A1), and sodium-coupled amino acid transporters (SLC38A2) were greater in HBCS than NBCS cows, due to upregulation in the postpartum phase. Those responses were similar to protein abundance of p-mTOR, which increased (BCS × Day effect) at 7 d in HBCS compared with NBCS cows. mRNA abundance of argininosuccinate lyase (ASL) and arginase 1 (ARG1) also was greater overall in HBCS cows. Together, these responses suggested impaired insulin signaling, coupled with greater postpartum AA transport rate and urea cycle activity in SAT of HBCS cows. An in vitro study using adipocyte and macrophage cocultures stimulated with various concentrations of fatty acids could provide some insights into the role of immune cells in modulating adipose tissue immunometabolic status, including insulin resistance and AA metabolism.

Distinct serum metabolomic signatures of multiparous and primiparous dairy cows switched from a moderate to high-grain diet during early lactation

INTRODUCTION

Feeding of high-grain diets is common in cows during early lactation, but increases the odds of metabolic derailments, which can likely be detected as undesirable shifts in the serum metabolome signature.

OBJECTIVES

The present study aimed to identify the metabolic signatures of the serum metabolome of early lactation dairy cows switched from a moderate to a high-grain diet.

METHODS

Targeted ESI-LC-MS/MS-based metabolomics was used to characterize metabolic alterations in the serum of early lactation multiparous (MP, n = 16) and primiparous (PP, n = 8) Simmental cows, according to parity and feeding phase. Data were analysed using different data mining approaches.

RESULTS

Carnitine, acetylcarnitine, propionoylcarnitine, amino acid related compounds cis-4-hydroxyproline, trans-4-hydroxyproline, proline betaine, lysophosphatidylcholine PC a C16:1 and phosphatidylcholine PC ae C36:0 were identified as the key metabolites distinguishing MP from PP cows. A different serum metabolite composition during moderate and high-grain diet was also evident. Notably, cows fed high grain diet had higher serum concentrations of primary bile acids and triglycerides, but lower levels of conjugated bile acids and carboxylic acids during the first week in grain. Amino acids valine, cystine and taurine together with lysophosphatidylcholine PC a C26:0 and several phosphatidylcholines were classified as important features for cluster separation.

CONCLUSIONS

Our study greatly expands earlier observations on dietary effects on serum metabolome composition of cows. The altered metabolomic fingerprints clearly distinguishable by diet and cow parity hold potential to be used as early diagnostic tools for cows experiencing grain-induced metabolic disturbances.

Short communication: Colostrum versus formula: Effects on mRNA expression of genes related to branched-chain amino acid metabolism in neonatal dairy calves

The objective of the current study was to elucidate the effect of feeding colostrum or milk-based formula on the tissue mRNA abundance of the most relevant branched-chain amino acids (BCAA) transporters and catabolizing enzymes in newborn calves. German Holstein calves were fed either colostrum (COL; n = 7) or milk-based formula (FOR; n = 7) with comparable nutrient composition but lower contents of free BCAA, insulin, and insulin-like growth factor-I in the formula than in the respective colostrum for up to 4 d of life. Tissue samples from liver, kidney fat, 3 different muscles [M. longissimus dorsi (MLD), M. semitendinosus (MST), and M. masseter (MM)], as well as duodenum, jejunum, and ileum were collected following euthanasia on d 4 at 2 h after feeding. The plasma-free BCAA were analyzed, and the tissue abundance of solute carrier family 1 member 5 (SLC1A5), SLC7A5, and SLC38A2 as well as mitochondrial isoform of branched-chain aminotransferase (BCATm), branched-chain α-keto acid dehydrogenase E1α (BCKDHA), and branched-chain α-keto acid dehydrogenase E1β (BCKDHB) were assessed. The preprandial plasma concentrations of free BCAA were affected by time but did not differ between groups. The plasma concentrations of free BCAA decreased in COL, whereas they increased in FOR after feeding, resulting in higher postprandial plasma total BCAA concentrations in FOR than in COL. The mRNA abundances of BCATm, BCKDHA, BCKDHB, as well as BCAA transporters in the liver, were not affected by the diet. In kidney fat, the mRNA abundance of BCAA catabolizing enzymes did not differ between groups, but that of SLC1A5 was lower in FOR than in COL. The mRNA abundance of BCAA catabolizing enzymes in different sections of the small intestine was not affected by the diet, whereas that of SLC7A5 was or tended to be lower in the duodenum, proximal jejunum, and mid jejunum of the COL calves compared with the FOR calves. The mRNA abundance of BCKDHA was lower in MLD and MM but greater in MS for the FOR calves compared with the COL calves. The mRNA abundance of SLC7A5 in MST was lower in FOR than in COL, whereas it was unaffected by the diet in MLD and MM. The differential effect of feeding colostrum on the mRNA abundance of BCKDHA in 3 different muscle tissues might point to a muscle type-specific response. The results also indicate that the colostral BCAA might be favorably used for anabolic metabolism in the small intestine of neonatal calves. Such effects are speculated to be due to the stimulatory effects of growth factors and hormones present in colostrum.

Profiling of circulating microRNA and pathway analysis in normal- versus over-conditioned dairy cows during the dry period and early lactation

The objective of this study was to determine the circulating microRNA (miRNA) profile in over-conditioned (HBCS) versus normal-conditioned (NBCS) dairy cows in combination with pathway enrichment analyses during the transition period. Thirty-eight multiparous Holstein cows were selected 15 wk before anticipated calving date based on their current and previous body condition scores (BCS) for forming either a HBCS group (n = 19) or a NBCS group (n = 19). They were fed different diets during late lactation to reach the targeted differences in BCS and backfat thickness until dry-off. A subset of 15 animals per group was selected based on their circulating concentrations of nonesterified fatty acids (on d 14 postpartum) and β-hydroxybutyrate (on d 21 postpartum), representing the greater or the lower extreme values within their BCS group. Blood serum obtained at d -49 and 21 relative to parturition (3 pools with 5 cows per each group and time point) were used to identify miRNA that were differentially expressed (DE) between groups or time points using miRNA sequencing. No DE-miRNA were discovered between NBCS versus HBCS. Comparing pooled samples from d -49 and d 21 resulted in 7 DE-miRNA in the NBCS group, of which 5 miRNA were downregulated and 2 miRNA were overexpressed on d 21 versus -49. The abundance of 5 of these DE-miRNA was validated in all individual samples via quantitative PCR and extended to additional time points (d -7, 3, 84). Group differences were observed for miR-148a, miR-122 as well as miR-455-5p, and most DE-miRNA (miR-148a, miR-122, miR-30a, miR-450b, miR-455-5p) were downregulated directly after calving. Subsequently, the DE-miRNA was used for bioinformatics analysis to identify putative target genes and the most enriched biological pathways. The most significantly enriched pathways of DE-miRNA were associated with cell cycle and insulin signaling as well as glucose and lipid metabolism. Overall, we found little differences in circulating miRNA in HBCS versus NBCS cows around calving.

Machine Learning Model Comparison in the Screening of Cholangiocarcinoma Using Plasma Bile Acids Profiles

Bile acids (BAs) assessments are garnering increasing interest for their potential involvement in development and progression of cholangiocarcinoma (CCA). Since machine learning (ML) algorithms are increasingly used for exploring metabolomic profiles, we evaluated performance of some ML models for dissecting patients with CCA or benign biliary diseases according to their plasma BAs profiles. We used ultra-performance liquid chromatography tandem mass spectrometry (UHPLC-MS/MS) for assessing plasma BAs profile in 112 patients (70 CCA, 42 benign biliary diseases). Twelve normalisation procedures were applied, and performance of six ML algorithms were evaluated (logistic regression, k-nearest neighbors, naïve bayes, RBF SVM, random forest, extreme gradient boosting). Naïve bayes, using direct bilirubin concentration for normalisation of BAs, was the ML model displaying better performance in the holdout set, with an Area Under Curve (AUC) of 0.95, 0.79 sensitivity, 1.00 specificity. This model, also characterised by 1.00 positive predictive value and 0.73 negative predictive value, displayed a globally excellent accuracy (86.4%). The accuracy of the other five models was lower, and AUCs ranged 0.75–0.95. Preliminary results of this study show that application of ML to BAs profile analysis can provide a valuable contribution for characterising bile duct diseases and identifying patients with higher likelihood of having malignant pathologies.

Discovery of different metabotypes in overconditioned dairy cows by means of machine learning

Using data from targeted metabolomics in serum in combination with machine learning (ML) approaches, we aimed at (1) identifying divergent metabotypes in overconditioned cows and at (2) exploring how metabotypes are associated with lactation performance, blood metabolites, and hormones. In a previously established animal model, 38 pregnant multiparous Holstein cows were assigned to 2 groups that were fed differently to reach either high (HBCS) or normal (NBCS) body condition score (BCS) and backfat thickness (BFT) until dryoff at -49 d before calving [NBCS: BCS < 3.5 (3.02 ± 0.24) and BFT < 1.2 cm (0.92 ± 0.21), mean ± SD; HBCS: BCS > 3.75 (3.82 ± 0.33) and BFT > 1.4 cm (2.36 ± 0.35)]. Cows were then fed the same diets during the dry period and the subsequent lactation, and maintained the differences in BFT and BCS throughout the study. Blood samples were collected weekly from 7 wk antepartum (ap) to 12 wk postpartum (pp) to assess serum concentrations of metabolites (by targeted metabolomics and by classical analyses) and metabolic hormones. Metabolic clustering by applying 4 supervised ML-based classifiers [sequential minimal optimization (SMO), random forest (RF), alternating decision tree (ADTree), and naïve Bayes-updatable (NB)] on the changes (d 21 pp minus d 49 ap) in concentrations of 170 serum metabolites resulted in 4 distinct metabolic clusters: HBCS predicted HBCS (HBCS-PH, n = 13), HBCS predicted NBCS (HBCS-PN, n = 6), NBCS predicted NBCS (NBCS-PN, n = 15), and NBCS predicted HBCS (NBCS-PH, n = 4). The accuracies of SMO, RF, ADTree, and NB classifiers were >70%. Because the number of NBCS-PH cows was low, we did not consider this group for further comparisons. Dry matter intake (kg/d and percentage of body weight) and energy intake were greater in HBCS-PN than in HBCS-PH in early lactation, and HBCS-PN also reached a positive energy balance earlier than did HBCS-PH. Milk yield was not different between groups, but milk protein percentage was greater in HBCS-PN than in HBCS-PH cows. The circulating concentrations of fatty acids (FA) increased during early lactation in both groups, but HBCS-PN cows had lower concentrations of β-hydroxybutyrate, indicating lower ketogenesis compared with HBCS-PH cows. The concentrations of insulin, insulin-like growth factor 1, leptin, adiponectin, haptoglobin, glucose, and revised quantitative insulin sensitivity check index did not differ between the groups, whereas serum concentrations of glycerophospholipids were lower before calving in HBCS-PH than in HBCS-PN cows. Glycine was the only amino acid that had higher concentration after calving in HBCS-PH than in HBCS-PN cows. The circulating concentrations of some short- (C2, C3, and C4) and long-chain (C12, C16:0, C18:0, and C18:1) acylcarnitines on d 21 pp were greater in HBCS-PH than in HBCS-PN cows, indicating incomplete FA oxidation. In conclusion, the use of ML approaches involving data from targeted metabolomics in serum is a promising method for differentiating divergent metabotypes from apparently similar BCS phenotypes. Further investigations, using larger numbers of cows and farms, are warranted for confirmation of this finding.

Plasma proteomic profiling and pathway analysis of normal and overconditioned dairy cows during the transition from late pregnancy to early lactation

This study applied a quantitative proteomics approach along with bioinformatics analyses to investigate changes in the plasma proteome of normal and overconditioned dairy cows during the transition period. Fifteen weeks before their anticipated calving date, 38 multiparous Holstein cows were selected based on their current and previous body condition scores (BCS) and allocated to either a high or a normal BCS group (19 cows each). They received different diets until dry-off to reach targeted differences in BCS and back fat thickness (BFT) until dry-off. At dry-off, normal BCS cows had a BCS <3.5 (minimum, 2.75) and BFT <1.2 cm (minimum, 0.58), and the high BCS cows had a BCS >3.75 (maximum, 4.50) and BFT >1.4 cm (maximum, 2.90). The proteomics study used a subset of 5 animals from each group. These cows were selected based on their circulating concentrations of fatty acids (FA) on d 14 postpartum and β-hydroxybutyrate (BHB) on d 21 postpartum, representing the greater or the lower extreme values within their BCS group, respectively. The high BCS subset (HE-HBCS) had 4.50 < BCS > 3.75, FA = 1.17 ± 0.46 mmol/L, and BHB = 2.15 ± 0.42 mmol/L (means ± SD), and the low BCS subset (LE-NBCS) had 3.50 < BCS > 2.75, FA = 0.51 ± 0.28 mmol/L, and BHB = 0.84 ± 0.17 mmol/L. Plasma samples from d -49, +7, and +21 relative to parturition were used for proteome profiling by applying the quantitative tandem mass tags (TMT) approach. Nondepleted plasma samples were subjected to reduction and digestion and then labeled with TMT 10plex reagents. High-resolution liquid chromatography-tandem mass spectrometry analysis of TMT-labeled peptides was carried out, and the acquired spectra were analyzed for protein identification and quantification. In total, 254 quantifiable proteins (criteria: 2 unique peptides and 5% false discovery rate) were identified in the plasma samples. From these, 24 differentially abundant proteins (14 more abundant, 10 less abundant) were observed in the LE-NBCS cows compared with the HE-HBCS cows during the transition period. Plasma α-2-macroglobulins were more abundant in HE-HBCS versus LE-NBCS cows at d +7 and +21. Gene Ontology enrichment analyses of differentially abundant proteins revealed that the acute inflammatory response, regulation of complement activation, protein activation cascade, and regulation of humoral immune response were the most enriched terms in the LE-NBCS group compared with the HE-HBCS group. In addition, we identified 24 differentially abundant proteins (16 in the LE-NBCS group, and 8 in the HE-HBCS group) during the transition period. The complement components C1q and C5 were less abundant, while C3 and C3d were more abundant in LE-NBCS compared with HE-HBCS cows. Overall, overconditioning around calving was associated with alterations in protein pathways related to acute inflammatory response and regulation of complement and coagulation cascades in transition cows.

Alterations of the acylcarnitine profiles in blood serum and in muscle from periparturient cows with normal or elevated body condition

The objective of the current study was to characterize muscle and blood serum acylcarnitine (AcylCN) profiles and to determine the mRNA abundance of muscle carnitine acyltransferases in periparturient dairy cows with high (HBCS) and normal body condition (NBCS). Fifteen weeks antepartum, 38 pregnant multiparous Holstein cows were assigned to 2 groups that were fed differently to reach the targeted BCS and backfat thickness (BFT) until dry-off at -49 d before calving (HBCS: BCS >3.75 and BFT >1.4 cm; NBCS: <3.5 and <1.2 cm). Thereafter, both groups were fed identical diets. Blood samples and biopsies from the semitendinosus muscle were collected on d -49, 3, 21, and 84 relative to calving. Actual BCS at d -49 were 3.02 ± 0.24 and 3.82 ± 0.33 (mean ± SD) for NBCS and HBCS, respectively. In both groups, serum profiles showed marked changes during the periparturient period, with decreasing concentrations of free carnitine and increasing concentrations of long-chain AcylCN. Compared with NBCS, HBCS had greater serum long-chain AcylCN in early lactation, which may point to an insufficient adaptation of their metabolism in response to the metabolic load of fatty acids around parturition. The muscle concentrations of C5-, C9-, C18:1-, and C18:2-AcylCN were lower and those of C14:2-AcylCN were greater in HBCS than in NBCS cows. The mRNA abundance of carnitine palmitoyltransferase (CPT)1, muscle isoform (CPT1b) and CPT2 increased from d -49 to early lactation (d 3, d 21), followed by a decline to nearly antepartum values by d 84; this change was not affected by group. In conclusion, over-conditioning around calving seems to be associated with mitochondrial overload, which can result in incomplete fatty acid oxidation in dairy cows.