Serum Metabolomic Profiling in Healthy Dogs Supplemented with Increasing Levels of Purified Beta-1,3/1,6-Glucans
Pedro Henrique Marchi, Leonardo de Andrade Príncipe, Felipe Sesti Trindade, Luana Dias dos Santos, Gabriela Luiza Fagundes Finardi, Eduarda Lorena Fernandes, Thaila Cristina Putarov, Gabriel Henrique Ribeiro, Luiz Alberto Colnago, Júlio Cesar de Carvalho Balieiro

TL;DR
This study shows how beta-glucan supplements affect the metabolism of healthy dogs, particularly in amino acid and energy pathways.
Contribution
The study provides new insights into how beta-glucans modulate canine metabolism through detailed serum metabolomic profiling.
Findings
Beta-glucan levels of 0.07% and 0.14% significantly influenced amino acid and energy metabolism in dogs.
Key metabolic pathways affected include phenylalanine, tyrosine, tryptophan, and glyoxylate metabolism.
Higher beta-glucan levels modulated lipid metabolism and gut microbiota-related pathways.
Abstract
Diet plays a crucial role in the health of dogs, and researchers are constantly exploring ways to improve nutrition through dietary supplements. This study investigated how beta-1,3/1,6-glucans, a nutraceutical known for its potential immune and gut health benefits, affect the metabolism of healthy dogs. Eight dogs participated in a controlled experiment using a double Latin square design, ensuring that each animal received all diet treatments in a balanced way. This study lasted a total of 140 days, divided into four experimental periods of 35 days each. The dogs were fed diets with different amounts of beta-glucans (0.0%, 0.07%, 0.14%, and 0.28%), and blood samples were collected at the end of each period to analyze metabolic changes. The results showed that certain levels of beta-glucans (0.07% and 0.14%) had a greater impact on metabolism, particularly in pathways related to amino…
Genes, proteins, chemicals, diseases, species, mutations and cell lines named across the full text — each resolved to its canonical identifier and authoritative record.
Click any figure to enlarge with its caption.
Figure 1
Figure 2
Figure 3
Figure 4
Figure 5
Figure 6- —Biorigin (Açucareira Quatá S.A.)
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsMetabolomics and Mass Spectrometry Studies · Gut microbiota and health · Human-Animal Interaction Studies
1. Introduction
Omics sciences aim to identify, characterize, and quantify the components involved in cellular mechanisms of information and gene expression, from DNA sequences to synthesize proteins and circulating metabolites. This field establishes the connection between an organism’s genotype and its phenotypic manifestation through the application of techniques such as genomics, transcriptomics, proteomics, lipidomics, and metabolomics [1,2,3]. Metabolomics makes it possible to characterize the metabolome of an organ, tissue, or specific sample, which can identify and quantify all of the low-molecular-weight compounds, collectively known as metabolites [4,5].
The metabolome corresponds to the set of metabolites present in a biological system, encompassing well-known molecules such as amino acids, carbohydrates, lipids, and nucleotides, as well as their respective metabolic pathways [6]. In this way, metabolomic analysis provides detailed information on aspects of an animal’s digestive physiology, such as the influence of different ingredients and gut microbiota on metabolism [7]. In this context, metabolomics has proven to be an effective tool for elucidating mechanisms and evaluating the effectiveness of dietary interventions in small animal nutrition [8].
Among the most relevant nutraceuticals for promoting health in dogs are beta-glucans, homopolysaccharides of glucose linked by beta 1→3 glycosidic bonds [9,10]. These bioactive compounds are present in the cellular structure of cereals, fungi, yeast, algae, and bacteria [11,12] and exhibit distinct molecular characteristics and biological functions depending on their source [13]. Among them, beta-glucans extracted from the cell wall of Saccharomyces cerevisiae are considered immunostimulant compounds [14]. Their molecular structure is characterized by the predominance of beta 1→3 linkages, with a lower proportion of beta 1→6 linkages and branches [15].
The beta-1,3/1,6-glucans are relevant immunomodulator compounds, with substantial evidence demonstrating their role in enhancing the immune status of dogs and other animals [16]. Due to their structure, they are recognized by specific receptors on myeloid-derived cells, such as dectin-1 [17]. This interaction stimulates innate and adaptive immune responses, activating various immune cells, including neutrophils, B cells, T cells, NK cells, and especially macrophages, while also regulating cytokine production [14,18].
In addition to the immunomodulatory effects, studies conducted by the same research group on the supplementation of beta-1,3/1,6-glucans purified from the cell wall of Saccharomyces cerevisiae have shown additional benefits for canine health. These include a reduction in plasma levels of basal glucose, serum cholesterol, and triglycerides in obese dogs, as well as the modulation of the intestinal microbiota in both healthy dogs and those with mild inflammatory bowel disease [19,20,21]. However, the mechanisms by which these compounds exert such metabolic and intestinal effects have yet to be fully elucidated.
In view of this, metabolomic approaches represent a promising strategy for deepening knowledge about the effects of beta-glucans in dogs. Therefore, this study aimed to conduct a metabolomic analysis of the serum of dogs subjected to supplementation with increasing levels of beta-1,3/1,6-glucans to generate evidence and gain a deeper understanding of the metabolic responses associated with this supplementation.
2. Materials and Methods
2.1. Ethical Statement
This study adhered to the Ethical Principles in Animal Research established by the Ethics Committee on Animal Use of the School of Veterinary Medicine and Animal Science, University of Sao Paulo (CEUA/FMVZ), under protocol number 2866090223.
2.2. Animals, Diets, and Experimental Design
The experiment was conducted at the Pet Nutrology Research Center (CEPEN Pet), part of the School of Veterinary Medicine and Animal Science of the University of Sao Paulo, Brazil. Eight dogs from the research facility were used, consisting of four Border Collies and four English Cocker Spaniels, both male and female, aged 3.5 ± 0.5 years, with an ideal body condition score [22]. They were assigned to two contemporary balanced Latin squares with four treatments, according to their breed. Their health status was confirmed through a comprehensive physical examination, nutritional history assessment, complete blood count, and biochemical profile tests, including albumin, glucose, total protein, urea, creatinine, alkaline phosphatase, cholesterol, triglycerides, aspartate aminotransferase (AST), and alanine aminotransferase (ALT).
All dogs used in this study were housed in the same facility, in collective kennels with 3.42 m^2^ of covered area and 7.21 m^2^ of uncovered area, with concrete floors and tiled walls. Dogs were released into grass-covered play areas twice daily to allow for social interaction. Five days prior to blood sampling, the animals were housed individually for fecal sample collection intended for other evaluations and were not released into the play areas during this period. To minimize stress, supervised walks were conducted twice daily in a clean, concrete area. All dogs had ad libitum access to fresh water.
Four nutritionally similar diets were formulated and extruded, differing only in beta-glucan content (0.0%, 0.07%, 0.14%, and 0.28%; Table 1). These levels of beta-glucans were chosen to encompass the doses reported in previous studies [16,19,20,21]. The beta-1,3/1,6-glucan used was a purified extract from the cell wall of Saccharomyces cerevisiae with a purity level of at least 60% (MacroGard^®^, Biorigin, Lençois Paulistas, Brazil). The dogs’ maintenance energy requirements were initially calculated using the equation 95 kcal × body weight^0.75^ [23]. This study lasted a total of 140 days, and each experimental period took 35 days, comprising 34 days of diet acclimation and one day for blood sampling. On day 35, 5 mL of blood was collected via jugular venipuncture to obtain serum for metabolomic analysis. Throughout the four periods, blood samples were consistently collected after 12 h of fasting to standardize the metabolic state of all animals.
2.3. Metabolomic Analysis
For each sample, an aliquot of 200 µL of serum was diluted in 400 µL of phosphate buffer (D2O-based PBS; 100 mM; pD = 7.4) containing 4.5 mM of sodium 3-(Trimethylsilyl)-1-propanesulfonic acid-d6 (DSS-d6, from Cambridge Isotopes, Leicestershire, UK), used as the nuclear magnetic resonance (NMR) internal standard. Posteriorly, a volume of 600 µL of each sample was transferred to a standard 5 mm NMR tube for NMR measurements. The 1D NMR experiments were performed on a Bruker Avance III 14.1 T spectrometer (Bruker Corporation, Karlsruhe, Baden-Württemberg, Germany) fitted with a 5 mm Broadband Observe probe with ATMA^®^ (Automatic Tunning Matching Adjustment), a z-field gradient, BCU-I variable temperature unit, field gradient generator unit, and a Sample-Xpress™ automatic sample changer.
Standard 1H spectra were acquired at 298.15 K using a NOESY-1D pulse sequence (named noesypr1d in the TopSpin Bruker software version 3.6.1) with water signal suppression by irradiation at 2822.79 Hz (O1). The acquisition parameters were as follows: number scans (ns) of 128, spectral width (sw) of 12,019 Hz (20.0276 ppm), 90° pulse (P1) of 13.53 µs, acquisition time (aq) of 4.52 s, relaxation delay (d1) of 4 s, data points (TD) of 108,544 (106 K), mixing time (d8) of 5 ms, and dummy scans (ds) of 4. The spectra were processed utilizing TopSpin™ 3.6.1 software (Bruker Biospin, Ettlingen, Germany). All the chemical shifts of 1H and 13C were observed in ppm in relation to the sign of DSS in 0.00 ppm as an internal reference, and an exponential enlargement of line of 0.3Hz was applied. The metabolites were assigned based on chemical displacements and multiplied by the signal using the software Chenomx 12.0 [25].
2.4. Statistical Analysis
The experimental design used was a contemporary double Latin square with four treatments. The statistical model included the fixed effect of treatment and the random effects of the animal within square, period, and residue. For pre-treatment, simple rectangular buckets of 0.04 ppm were calculated using the Mvova V.15 software (www.mestrelab.com, accessed on 4 February 2024) with the auto-integration mode and scaling to total intensity. The analysis covered regions between 0.2 and 10 ppm, excluding the residual water signal (4.7–5.1 ppm). Metabolome data were uploaded to MetaboAnalyst 5.0 [26], median-normalized, and Pareto-scaled before statistical and bioinformatics analyses (Figures S1 and S2).
The first principal component analysis (PCA, 32 spectra × 235 columns) was performed using the intensity of all obtained metabolites and PERMANOVA validation. Subsequently, key metabolites from the loading analysis of each PC were identified and used for another PCA and hierarchical cluster analysis to compare metabolite intensity between treatments. Additionally, pathway analysis was conducted to identify the most relevant pathways associated with the identified metabolites. Compound names were standardized according to HMDB ID. The pathway analysis parameters included the global test and relative-betweenness centrality, with the Homo sapiens library selected. Metabolic pathways related to the identified metabolites were mapped based on p-values (p ≤ 0.10) and pathway impacts (impact ≥ 0.001).
3. Results
Figure 1 represents the spectrum generated by the pooled samples from all experimental units and demonstrates the regions and peaks of the following main serum metabolites: formate, histidine, phenylalanine, tyrosine, glucose, lactate, carnitine, creatine, creatine phosphate, glutamine, citrate, alanine, isoleucine, leucine, threonine, valine, low-density lipoprotein, very-low-density lipoprotein, and high-density lipoprotein. The PCA (Figure 2A) revealed some overlaps between the treatments, in which PC1 explained 24.4% and PC2 explained 9.5% of the total variance of the dataset. Moreover, the PCA loading plot analysis indicated some possible associations between key metabolites and each PC (Figure 2B; Table 1). The top 15 key metabolites observed in PC1 were related to lipids, while those related to PC2 were identified and listed in Table 2.
The complete dataset analysis revealed 15 different lipids associated with PC1 and 12 main metabolites linked to PC2. Based on the 12 identified metabolites, a new PCA was performed (Figure 3A, p < 0.05), revealing additional regions of overlap between the treatments. PC1 accounted for 40.8% of the total variance, while PC2 captured 16.5% (Figure 3B, p < 0.05). Furthermore, analysis of variance showed no significant differences in metabolite abundance between groups (Figure 4, p > 0.05). However, hierarchical heatmap analysis indicated differences in the intensity of the 12 metabolites across the treatments (Figure 5, p < 0.05), with the control group displaying greater variability and lower intensity for most compounds compared to the treatment groups.
Based on these metabolite intensity differences, an enrichment analysis was conducted to identify the most relevant metabolic pathways (Figure 6, p < 0.05). The main pathways are represented by darker nodes (higher p-values) and larger sizes (greater impact analysis). All identified pathways had p-values ≤ 0.10 and pathway impacts ≥ 0.001 (Table 3, p < 0.05). These pathways were associated with protein metabolism—including phenylalanine, histidine, and alanine, aspartate, and glutamate metabolism—as well as phenylalanine, tyrosine, and tryptophan biosynthesis and energy metabolism, which included glyoxylate and dicarboxylate metabolism and starch and sucrose metabolism.
4. Discussion
The application of metabolomics to animal nutrition poses challenges in data interpretation due to the large volume of information generated [1]. To address this complexity, different analysis strategies can be employed depending on the study’s aim. The untargeted analysis approach used in this study involves the global profiling of all metabolites present in a sample, without prior selection [27,28]. To organize and facilitate data interpretation, principal component analysis (PCA) is the most commonly used technique for metabolomic data analysis [29]. In this study, the metabolites identified through nuclear magnetic resonance were plotted in loading graphs and evaluated for potential discriminatory patterns between the treatments. In the initial PCA performed using all metabolites from the dataset, a significant overlap was observed. The metabolites with the highest loadings on PC1 were all classified as lipids but could not be individually identified due to technical limitations and the physicochemical characteristics of fats (Han & Gross, 2022) [30]. However, the analysis of PC2 enabled the identification of twelve metabolites with the highest scores, which were subsequently examined in greater detail. The absence of a clear and complete separation in the metabolome profile may be attributed to the subtle inclusion of beta-glucan as the only modified ingredient among the experimental groups. This study was specifically designed to isolate, as much as possible, the effects of beta-glucan inclusion on canine metabolism. With this in mind, it is understandable that a more pronounced separation in the serum metabolome would be unlikely, especially considering that the primary site of beta-glucan action is the intestine. A similar pattern was observed by Ephraim et al. [31], who reported that variations in dietary protein levels did not lead to a clear separation in the metabolite profile.
As a subsequent step, a more detailed analysis was conducted by selecting only the main metabolites. The comparison of mean values before and after the normalization of the twelve most relevant metabolites, along with the overlaps observed in the principal component analysis, indicated minimal differences between the groups. However, hierarchical clustering with a heatmap provided a clearer distinction, highlighting the importance of specific metabolites in differentiating treatments. The heatmap revealed evident differences between the control and beta-glucan groups, where red indicated positive correlations and darker shades of blue represented negative correlations. Overall, all levels of beta-glucan inclusion showed differences in intensity when compared to the control group. Nonetheless, the 0.07% and 0.14% beta-glucan inclusions exhibited the most pronounced differences. This finding corroborates previous studies by our research group, which, although employing different evaluation methods, also concluded that supplementation of approximately 10 mg/kg of beta-glucan promotes beneficial effects in dogs [19,20]. Furthermore, the first results of this experiment indicate that the 0.14% inclusion level (approximately 10.6 mg/kg) led to the best immune response and a favorable fecal microbiota profile [21].
Among the most relevant metabolites, cholesterol and L-acetylcarnitine stood out as being strongly associated with changes in the metabolic profile of dogs in the 0.14% beta-glucan group. Furthermore, the 0.07% and 0.28% groups showed lower signal intensities, characterized by an inverted U-shaped curve. This suggests that the optimal range for beta-glucan inclusion and metabolism lies around 0.14% (or 10 mg/kg), particularly in relation to lipid metabolism. Cholesterol is an essential component of cell membrane structure and lipid metabolism. L-acetylcarnitine plays a crucial role in beta-oxidation, facilitating the transport of acetyl-CoA across mitochondrial membranes for fatty acid oxidation [32]. The interaction between these metabolites suggests a potential influence on lipid metabolism regulation and the physiological response of dogs in the 0.14% beta-glucan group [33,34]. In view of this, studies have shown that beta-glucan inclusion can reduce basal and mean plasma concentrations of glucose, basal insulin, triglycerides, and cholesterol in obese dogs to levels similar to those of healthy dogs [10,19,35]. Incorporating beta-glucans into the diet may be a promising strategy to mitigate dyslipidemia, a condition often associated with hyperinsulinemia and obesity [36].
In addition, this study provided insights into alanine, aspartate, and glutamate metabolism (most significant p-value and high pathway impact), highlighted by the production of alanine via the transamination of pyruvate in the presence of glutamate [37]. Alanine plays a crucial role in metabolic regulation, protein synthesis, and ATP supply and can be synthesized from other amino acids, such as valine [38]. In the present study, the heatmap revealed a positive correlation between valine and the 0.07% beta-glucan group. Moreover, these metabolites play a central role in nitrogen metabolism, as compounds such as aspartate, glutamate, and glutamine are synthesized from ammonia during the cycle. Another relevant aspect is the utilization of carbon atoms under extreme conditions for glucose production [37]. Glutamine is involved in virtually all metabolic pathways of proliferating cells [39]. In the present study, glutamine concentration was significantly reduced in the 0.14% group, which may be related to a more pronounced use of this amino acid as an energy source by intestinal cells. Glutamine is the only amino acid extracted from the bloodstream for this purpose and, when metabolized in the intestine, produces alanine [37], which aligns with the observed behavior of both metabolites in the 0.14% group.
In this study, the biosynthesis of phenylalanine, tyrosine, and tryptophan exhibited the highest pathway impact value. These aromatic amino acids are essential for dogs, except for tyrosine, which is considered semi-essential, as it can be synthesized endogenously from phenylalanine when available [40,41]. Their biosynthesis occurs exclusively in microorganisms and plants via the shikimate pathway [42]. Given that beta-glucans can modulate gut microbiota composition [21], it is likely that these changes enhanced the microbial production of aromatic amino acids, which were subsequently absorbed and became metabolically available in the bloodstream of the dogs. This is a significant finding considering the crucial role that these essential amino acids play in canine health. Both phenylalanine and tyrosine are extensively studied for their neurological functions, serving as precursors to neurotransmitters such as norepinephrine, epinephrine, and dopamine [43]. Highlighting its beneficial effects on brain function, Kano et al. [44] observed that oral tyrosine supplementation positively influenced dogs’ cognitive performance, particularly in tasks related to learning and attention. Beyond its neurological role, tyrosine also contributes to various physiological processes, including coat quality and pigmentation [45,46]. A deficiency of tyrosine, or even its precursor phenylalanine, can lead to coat fading, as metabolites such as eumelanin and pheomelanin, the primary melanin pigments, are synthesized through its metabolic pathway [47]. Additionally, L-phenylalanine, a bioactive form of phenylalanine, has been shown to support intestinal barrier function and act as an inhibitor of epithelial alkaline phosphatase, reducing endotoxemia and inflammation [48,49].
Tryptophan metabolites perform various functions in the body, especially in the central nervous system, where they act as precursors of serotonin and melatonin [50,51,52]. However, the percentage of tryptophan converted into serotonin is very low, whereas approximately 90% of this amino acid is metabolized through the kynurenine pathway in mammals [51,53]. Kynurenine is closely linked to immune and inflammatory responses [52]. During inflammatory conditions, such as colitis, metabolites from this pathway can modulate the immune response by activating the aryl hydrocarbon receptor (AhR). The AhR plays a key role in mediating almost all toxicological effects of aromatic hydrocarbons [54]. Additionally, indoles produced through this pathway have been shown to reduce TNF-alpha activity while increasing the expression of the anti-inflammatory cytokine IL-10 [55]. A study on oat-derived beta-glucans suggests that beta-glucan consumption may alter tryptophan metabolism and could be responsible for reducing the incidence of colorectal cancer in rats [56]. In dogs, kynurenic acid inhibits intestinal motility and may exert an anti-inflammatory effect by inhibiting xanthine oxidase [57]. Additionally, Guard et al. [58] identified that dogs with acute diarrhea exhibited a reduction in urinary tryptophan metabolites. In this sense, an increase in intestinal tryptophan biosynthesis may represent an additional mechanism through which beta-glucans modulate immune activity. In contrast, supplementation with oat-derived beta-1,3/1,4-glucans was shown to reduce plasma kynurenine concentrations in dairy calves [59]. Although this represents a divergent finding, it reinforces the complexity of the microbiome–immune-system interaction. Emerging evidence suggests that tryptophan and kynurenine are key intermediates in these pathways, which have been increasingly recognized for their roles in modulating immune and inflammatory responses.
Additionally, beta-glucan consumption appears to have influenced glyoxylate metabolism. Although its pathway impact value was below 0.1, increased activity in this pathway may also be interconnected with the gut microbiota. This cycle involves the utilization of acetate for glucose production and is generally absent in dog metabolism, occurring only in plants and microorganisms [60]. In the study published by Marchi et al. [21], no differences were observed in fecal short-chain fatty acid levels after the consumption of beta-glucans at any concentration. However, there was a significant increase in Erysipelotrichaceae, Ruminococcaceae, Faecalibacterium, and Prevotella, which are known producers of short-chain fatty acids [61]. Furthermore, in a study conducted with older hens that evaluated the effects of beta-1,3/1,6-glucan supplementation on the gut microbiome and metabolome, the modulation of several beneficial bacterial families, as well as positive modulation in glyoxylate metabolism, were also observed [62]. Nevertheless, due to limited supporting evidence, it is not possible to establish concrete hypotheses regarding this finding. The application of metabolomics in animal nutrition provides a comprehensive understanding of the effects of specific ingredients and enables the identification of risk factors for various diseases [63,64]. By offering a broader perspective on metabolism, this approach supports the discovery of biomarkers [65]. To better explore the therapeutic potential of the hypotheses discussed in this paper, further studies involving animals with conditions such as obesity, diabetes, or dyslipidemia are needed. Additionally, integrating metabolomic analysis with the assessment of other parameters will allow for a deeper correlation of metabolic alterations, enhancing our understanding of the underlying mechanisms of these diseases and how beta-1,3/1,6-glucans interact with them.
5. Conclusions
This study highlighted metabolic changes in dogs that consumed increasing levels of purified beta-1,3/1,6-glucans from the cell wall of Saccharomyces cerevisiae. The main modulated pathways included the metabolism of alanine, aspartate, and glutamate; glyoxylate and dicarboxylate; and phenylalanine, tyrosine, and tryptophan. These alterations in lipid and amino acid metabolism, as well as their relationship with gut microbiota, may contribute to a better understanding of the mechanisms of action of beta-glucans on canine health. Moreover, the results indicate a better metabolic response following the consumption of the 0.14% diet (approximately 10 mg/kg), which is consistent with previous findings in the literature.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Roessner U. Bowne J. What is metabolomics all about?Bio Techniques 20094636336510.2144/00011313319480633 · doi ↗ · pubmed ↗
- 2Aleksandrova K. Egea Rodrigues C. Floegel A. Ahrens W. Omics biomarkers in obesity: Novel etiological insights and targets for precision prevention Curr. Obes. Rep.2020921923010.1007/s 13679-020-00393-y 32594318 PMC 7447658 · doi ↗ · pubmed ↗
- 3Avela H.F. Sirén H. Advances in lipidomics Clin. Chim. Acta 202051012314110.1016/j.cca.2020.06.04932622966 · doi ↗ · pubmed ↗
- 4Zhou Z. Luo M. Chen X. Yin Y. Xiong X. Wang R. Zhu Z.-J. Ion mobility collision cross-section atlas for known and unknown metabolite annotation in untargeted metabolomics Nat. Commun.202011433410.1038/s 41467-020-18171-832859911 PMC 7455731 · doi ↗ · pubmed ↗
- 5Qiu S. Cai Y. Yao H. Lin C. Xie Y. Tang S. Zhang A. Small molecule metabolites: Discovery of biomarkers and therapeutic targets Signal Transduct. Target. Ther.2023813210.1038/s 41392-023-01399-336941259 PMC 10026263 · doi ↗ · pubmed ↗
- 6Bar N. Korem T. Weissbrod O. Zeevi D. Rothschild D. Leviatan S. Kosower N. Pompan-Lotan M. Weinberger A. Le Roy C.I. A reference map of potential determinants for the human serum metabolome Nature 202058813514010.1038/s 41586-020-2896-233177712 · doi ↗ · pubmed ↗
- 7Patti G.J. Yanes O. Siuzdak G. Innovation: Metabolomics: The apogee of the omics trilogy Nat. Rev. Mol. Cell Biol.20121326326910.1038/nrm 331422436749 PMC 3682684 · doi ↗ · pubmed ↗
- 8Carlos G. dos Santos F.P. Fröehlich P.E. Canine metabolomics advances Metabolomics 2020161610.1007/s 11306-020-1638-731955274 · doi ↗ · pubmed ↗
