Next Article in Journal
Impact of Feeding Pattern on the Structure and the Economic Performance of Dairy Cow Sector
Next Article in Special Issue
Quality Control in Fiore Sardo PDO Cheese: Detection of Heat Treatment Application and Production Chain by MRI Relaxometry and Image Analysis
Previous Article in Journal
Validation and Use of the RumiWatch Noseband Sensor for Monitoring Grazing Behaviours of Lactating Dairy Cows
Previous Article in Special Issue
Functional Odd- and Branched-Chain Fatty Acid in Sheep and Goat Milk and Cheeses
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

LC-QTOF/MS Untargeted Metabolomics of Sheep Milk under Cocoa Husks Enriched Diet

1
Dipartimento di Scienze della Vita e dell’Ambiente, Campus Universitario, 09042 Monserrato, Italy
2
Dipartimento di Agraria, Sezione di Scienze Zootecniche, Università degli Studi di Sassari, Viale Italia 39, 07100 Sassari, Italy
*
Author to whom correspondence should be addressed.
Submission received: 15 December 2020 / Revised: 18 January 2021 / Accepted: 26 January 2021 / Published: 22 February 2021
(This article belongs to the Special Issue Innovation Meets Tradition in the Sheep and Goat Dairy Industry)

Abstract

:
The aim of this work was to evaluate, by an untargeted metabolomics approach, changes of milk metabolites induced by the replacement of soybean hulls with cocoa husks in the ewes’ diet. Animals were fed with a soybean diet integrated with 50 or 100 g/d of cacao husks. Milk samples were analyzed by an ultra high performance liquid chromatograph coupled to a time of flight mass spectrometer (UHPLC-QTOF-MS) platform. Multivariate statistical analysis showed that the time of sampling profoundly affected metabolite levels, while differences between treatments were evident at the fourth week of sampling. Cocoa husks seem to induce level changes of milk metabolites implicated in the thyroid hormone metabolism and ubiquinol-10 biosynthesis.

1. Introduction

In the Mediterranean basin, where dairy sheep and goats breeding is widespread, the use of agricultural by-products in the diet of small ruminants is an ancient practice. Today the great availability of agro-industrial by-products produced worldwide opens the door to other products to be tested. Noteworthy, the use of co-feeds can contribute to reducing the impact associated with ruminant’s management. Several positive aspects were observed when agro-industrial by-products (e.g., grape, olive, tomato, citrus pulp and myrtle residues) were included in the diet of small dairy ruminants. In particular, it has been reported that the diet implementation with these by-products has beneficial effects on ruminal metabolism, animal health, and quality of derived products [1,2,3].
Cocoa husks represent the part of cocoa pod left over and are the principal by-products derived from Theobroma cacao L., representing an important crop for many tropical developing countries. Cocoa husks are obtained after the removal of the cocoa beans, the principal commercial product of cocoa processing [3]. Husks represent 70–75% of the fruit weight [4] and, besides limited local use, are considered as waste [5]. Cocoa pod husks along with other industrial by-products such as cocoa bean shells, cocoa bean meal and cocoa germ can be used as feed especially in the developing countries.
Purine alkaloids, including theobromine, theophylline and caffeine, and other methyl-xanthines, play a substantial role in pharmacology and food chemistry. A limited number of plant species accumulates purine alkaloids, such as caffeine and theobromine, which are synthesized from xanthosine, a catabolite of purine nucleotides. The most widely distributed methylxanthine in the plant kingdom is the antioxidant caffeine, which accumulates in leaves and seeds of tea (Camellia sinensis), coffee (Coffea arabica) and a limited number of other species. Considerable amounts of theobromine are stored in the seeds of cacao (Theobroma cacao) [6]. Moreover, from a pharmacokinetic point of view, theobromine can be easily adsorbed, with a low protein binding affinity and distributed into body tissues with no specific tissue accumulation [7]. Despite the positive health benefits of consuming foods containing purine alkaloids [8] one of the limits in the use of cocoa husks as supplementary feed for animals is the fact that in certain amount, theobromine, can be toxic [9]. The EU Scientific Opinion of the Panel on Contaminants in the Food Chain reported that when exposed to theobromine dairy cows showed reduced milk yield, increase in fat levels, and adverse effects such as hyperexcitability, sweating, increase in respiration and heart rates [10]. Reported levels of theobromine in cocoa husk meal are in the range of 1.5–4.0 g per kg of feed material and UE regulations set for ruminants the maximum levels of theobromine in feed materials at 300 mg/kg. On the other side, cocoa husks on a dry matter basis (DM) are a good source of neutral detergent fiber (NFC 38–44%) and high variable concentrations of non-fiber carbohydrate (NFC 17.5–47% on DM), crude proteins (CP 2.1–9.1% on DM), and lipids (ethyl esters 0.6–4.7% on DM); phenolic compounds account for 4.6–6.9% of gallic acid/100 g DM. However, the high content of acid detergent lignin (19.6% on DM) derived from cocoa pericarp affects the digestibility of this by-product, limiting its further use in animal feeding [11].
Metabolomics is a largely used branch of omics sciences used to explore the composition, relative levels, dynamics, and interactions amongst metabolites in an organism or biological system in response to various stimuli, diet or treatments. In the field of dairy, metabolomics has been used for example in the authentication of beef production systems [12], in the milk to investigate metabolite differences of healthy, subclinical, and clinical mastitis cows [13] or sheep milk grazed on different grazing systems [14], while in cheese samples used for discriminate cheese produced from raw or thermized milk [15].
To our knowledge, experimental data on theobromine carry over in milk are not available, as well as the metabolic effects of cocoa husks diet supplementation. In a previous work, milk and blood parameters of ewes supplemented with different doses of cocoa husks were evaluated [5]. In this experiment DMI and milk yield were not affected; however, milk protein content increased with cocoa husk supplementation. Furthermore, the cocoa husks supplementation lowered milk somatic cell count (SCC), suggesting a beneficial role of these by-products on mammary health status. Both systemic and local health conditions were good as evidenced by hematological parameters and electrophoretic profile of serum protein fractions [5]. Both systemic and local health conditions were good as evidenced by hematological parameters and electrophoretic profile of serum protein fractions [5]. Taking into consideration these results, we decided to investigate the potential beneficial nutritional outcomes by exploring, for the first time, milk metabolite profile changes of healthy ewes fed with a diet supplemented with 100 g/d per head soybean hulls that were replaced with 50 and 100 g/d per head of cocoa husks, respectively. In order to pursue these objectives, we used a liquid chromatography mass spectrometry metabolomics platform and multivariate statistical analysis.

2. Materials and Methods

The experiment was approved by the Ethics Committee of the University of Sassari (no. 54584/2018).
Animals and diets. As previously described [5], twenty-four Sarda dairy ewes in mid lactation stage (120 days in milking), at third parity, with an average body weight of 42.5 kg were divided into three homogenous groups. Each experimental group received a total mixed ration (2.58 kg/d per head) as the basal diet and received (a) a supplement of 100 g/d per head of soybean hulls (CH0 group); (b) half of the soybean hulls supplement was replaced with 50 g/d per head of cocoa husks (CH50 group); and (c) soybean hulls was totally replaced with 100 g/d per head of cocoa husks (CH100 group). Theobromine concentration was 130 ± 12 and 253 ± 27 mg/kg of DM for the CH50 and CH100 diets, respectively. The CH was administered individually during milking, mixed with beet pulp to increase palatability. The experiment lasted 8 weeks (from May to June 2018), with the last 4 weeks of sampling collection. Individual milk samples were collected weekly at morning milking.
Chemicals. Sigma Aldrich (Milan, Italy) reagents were used for the analysis. Bi-distilled water was obtained with a MilliQ purification system (Millipore, Milan, Italy).
Sample preparation for UHPLC-QTOF/MS analysis. Individual milk samples from the morning milking were stored into sterile plastic Falcon tubes at −80 °C before analysis. Prior to UHPLC-QTOF/MS analysis, the 107 individual sheep milk samples were thawed and thoroughly vortex mixed. A total of 150 µL of milk samples were transferred to Eppendorf tube containing 10 uL of the internal mixture of standards (Splash, Lipidomics, Sigma Aldrich, Milan, Italy) and added with 525 µL of methanol and 525 µL of MTBE. Samples were then mixed by vortexing for 1 min. Next, samples were centrifuged at 4000 rpm for 15 min. Subsequently, a second round of centrifuge was performed at 12000 rpm for 5 min. Before transferring to autosampler vials, the supernatant was filtered through a 0.22 µm MS nylon syringe filter.
UHPLC-QTOF/MS analysis. After extraction with MTBE/methanol, the supernatant of milk samples was analyzed with a 6560 Q-TOF (Agilent Technologies, Palo Alto, CA, USA) coupled with an Agilent 1290 Infinity II LC system. An aliquot of 1.0 μL from each sample was injected in a Kinetex 5 µm EVO C18 100 A, 150 mm × 2.1 μm column (Agilent Technologies, Palo Alto, CA, USA). The column was maintained at 50 °C at a flow rate of 0.5 mL/min. The mobile phases consisted of (A) methanol:water (90:10, v/v) with ammonium acetate (10 mM) and (B) acetonitrile:methanol:2-propanol (20:30:50, v/v) with ammonium acetate (10 mM). The chromatographic separation was obtained using the following gradient—0 min 70% B kept for 1 min; 1–3.5 min 86% B; 3.5–10 min 86% B; 10.1–17 min 100% B; 17.1–19 min 70% B. The analytical setup used was equipped with an Agilent jet stream technology source which was operated in both positive and negative ion modes with the following parameters: gas temperature,200 °C; gas flow (nitrogen) 10 L/min; nebulizer gas (nitrogen), 50 psig; sheath gas temperature, 300 °C; sheath gas flow, 12 L/min; capillary voltage 3500 V for positive and 3000 V for negative; nozzle voltage 0 V; fragmentor 150 V; skimmer 65 V, octapole RF 7550 V; mass range, 40–1700 m/z; capillary voltage, 3.5 kV; collision energy 20 eV in positive and 25 eV in negative mode. An Agilent MassHunter software was used for instrument control (revision B.09.00).
Multivariate statistical data analysis (MVA). The QTOF-MS data were submitted to the web platform XCMS [16]. The workflow allows feature detection, retention time correction, alignment, annotation, statistical analysis, and data visualization. Retention time, m/z values, and intensities of each feature obtained from XCMS were then submitted to multivariate statistical analysis (MVA), as implemented in SIMCA-P+ software (version 14.1, Umetrics, Umeå, Sweden). Prior to MVA, QTOF-MS features were mean centered and scaled to unit variance column-wise. Principal component analysis (PCA) was performed to investigate sample distributions, deviating features and prevailing trends. This technique reduced the dimensionality of data set retaining most of the variance. Results are shown as scatter plots of scores and loadings in the first principal components, reporting sample and variable displacement in the hyperplane, respectively. The partial least squares-discriminant analysis (PLS-DA) and its orthogonal variant (OPLS-DA) were performed for classification of samples and identification of the most discriminant variables. The quality of the models was evaluated based on the cumulative parameters R2Y and Q2Y, the latter estimated by cross validation, and tested for overfitting using a y-table permutation test (n = 400). The variable importance in projection (VIP) score summarizes the contribution of each variable to the model. The VIP scores in the predictive component were analyzed and only those metabolites having VIP values > 1 were considered as discriminant between the classes [17].

3. Results and Discussion

This untargeted metabolomics approach on dairy ewes was conceived to provide a molecular and biological insights on milk metabolite changes induced by a diet integrated with 50 or 100 g/d of cacao husks. Moreover, we designed this cost-efficient metabolomics study to measure predictive metabolite biomarkers to be used in exploring feed suitability and efficiency.
For these purposes, a total of 107 sheep milk samples were analyzed with a UHPLC-QTOF/MS platform, (Figure S1 reported the LC-MS chromatograms). After analysis, raw mass spectrometry data acquired were then uploaded to the XCMS platform, which generated an average of 7700 and 3300 features for the positive and negative modes, respectively. The XCMS outputs were submitted to multivariate statistical analysis using the software SIMCA.
A PCA was performed to explore similarities of milk sample metabolic profiles of the three diets over the experimental sampling period. Results indicated that metabolic similarities of milk samples were mainly due to the time of sampling, as clearly visible in the score plots shown in Figure 1 for the positive and negative ionization modes.
Information collected for the two ionization modes yielded similar results (i.e., samples collected in the two first weeks of treatment clustered together as a result of similar metabolite profile, while milk samples from the third and the last week of experiment can be easily differentiated along the second principal component (PC2) in the positive and negative ionization modes). Interestingly, the same pattern was observed for milk yield, as shown in Figure S2, whereby the first two weeks gave comparable average values for all the group samples, followed by a drastic milk yield drop after the second week of sampling (from the 13 to 26 of June). A decrease of milk yield and lactose content and an increase in fat and protein content, due to a milk concentration effect, during the late-spring/summer seasons were already reported as a natural evolution of the lactation curve in dairy Sarda sheep [18]. Probably these milk qualitative and quantitative changes due to lactation seasonality have caused an overall modification of the metabolic profiles captured by the PCA.
Attempts to differentiate samples based on the diets (CH0, CH50, and CH100) were carried out by performing three-classes PLS-DA for each time point for positive and negative modes. The multivariate statistical analysis failed to discriminate samples for the first three weeks of sampling (data not shown) whereas at the fourth week of treatment, samples were properly classified, with high statistical confidence (Figure S3), indicating that the effects of the cocoa husks supplementation on milk metabolic profile were evident only in the last time point. Subsequently, with the purpose of finding discriminant metabolites, pairwise OPLS-DA were carried out comparing the three groups of milk samples collected in the last week of treatment (Figure 2). In Table 1 and Table 2 we reported the list of discriminant metabolites, their biological significance, and regulation in sheep milk samples.
From the analysis of discriminant metabolites, as summarized in Table 1 and Table 2, we found that the addition of cocoa husks in the dairy lactating sheep diet affected levels of different metabolites which are involved in different metabolic pathways, including ubiquinol synthesis and thyroid hormone metabolism. In particular, we found that in milk samples of lactating ewes the thyroid hormone L-thyroxine (T4) and tetraiodothyroacetate, or tetrac (T4A), were down regulated for both CH50 and CH100 groups, thus suggesting a modification of the thyroid hormone metabolism. Levels of T4 and T4A are shown in Figure 3, indicating that, in the first two weeks of treatment, group samples showed similar average values. Levels where higher in the third week and in the last week CH50 and CH100 groups showed lower values when compared to CH0.
The thyroid hormones T4 and 3,5,3′-triiodo-L-thyronine (T3) are regulated at many metabolic levels. In addition to deiodination (thyroid hormone metabolism I pathway), other routes of thyroid hormone metabolism have been described and metabolically active metabolites have been identified. In the thyroid hormone metabolism pathway (via degradation), small amounts of T4 may also be oxidatively deaminated at the alanine side chain to yield T4A. The oxidative deamination due the kidney L-amino acid oxidase seems the most probable explanation [19,20]. The acetic acid analogues of thyroid hormones are formed in the liver and other peripheral tissue by a two-stage process involving transamination followed by oxidative decarboxylation [21]. T4A may be conjugated to form ether glucuronides or ester glucuronides of these metabolites may form by conjugation of the phenolic hydroxyl group, or carboxyl group respectively [20].
Thyroid hormones in milk should derive from circulating blood, then undergoing both desiodative and non-desiodative metabolism. A reduction in milk T4 could; therefore, be derived from a reduction in circulating T4 levels or from the reduction in T4 level available at the mammary gland, due to the effect of interferents on the transporters of thyroid hormones, or, finally, by increased desiodation of T4. The reduction of milk T4A level seems the logical direct consequence of the reduction of T4, deriving the T4A directly from T4. Thyroid hormones play a relatively important role in lactational processes. In small ruminants and, in particular, in free grazing animals which are subjected to seasonality, endogenous (breed, age, gender, physiological state) and environmental factors (climate, season, with a primary role of nutrition) are able to impact on thyroid activity and hormone blood levels [22]. During heat stress, an increase of blood T4 concentration and a decrease of milk production were reported in sheep [23,24]. This trend was evident also in our milk samples, which were collected in June, where, during the experimental period, in the CH0 group the milk T4 levels inversely follow the milk production trend (Figures S2 and S4). Differently from CH0, in the last week of sampling, milk T4 levels dropped in the CH50 and CH100 groups (Figure S4), thus indicating a goitrogenic effect of cocoa husks on milk T4. In the rat, Wolff and Varrone (1969) reported that methyl-xanthines (i.e., theobromine, theophylline, and caffeine) showed a goitrogenic effect and are; therefore, to be considered weak anti-thyroid agents if present in the diet at concentrations varying between 5–20 microM [25]. Further studies will be carried out to better evaluate this effect in sheep, measuring the volume of the thyroid glands and the concentration of T4 and TSH in ewes’ blood after consumption of cocoa husks. At the milk level, which only contains traces of iodothyronines, which; therefore, do not represent a significant nutrient of the milk itself, as opposed to iodine, any reduction in concentration should not be particularly relevant for its nutritional value.
Notably, phenylalanine was found up-regulated in CH50 and CH100 milk samples. Considering that phenylalanine is the precursor of tyrosine and the thyroid hormone T4 is a tetraiodinated derivative of tyrosine [22], this fact confirms a dysregulation of thyroid hormones metabolism.
The discriminant metabolites: 2-methoxy-6-(all-trans-octaprenyl)phenol, 6-methoxy-3-methyl-2-all-trans-decaprenyl-1,4-benzoquinol, and 3,4-dihydroxy-5-all-trans-decaprenylbenzoate, involved in the ubiquinol-10 biosynthesis, were found down regulated in milk samples of the CH50 and CH100 groups with respect to CH0 (Table 1). Coenzyme Q (CoQ10), also known as ubiquinone, is an essential component of the mitochondrial respiratory chain and participates in aerobic cell respiration. The reduced ubiquinol form of CoQ10 is a chain-breaking antioxidant, decreasing oxidative damage caused by lipid peroxidation within mitochondria. The cocoa husks feeding regime, showed markedly lower levels of intermediate metabolites of the CoQ10 pathway compared to the controls. To the best of our knowledge, no studies of theobromine effects on CoQ10 pathway have been widely reported, but interference of phenolic compounds has been found in liver of rats [26]. Specifically, the authors found that pentagalloylglucose, which is a component of tannic acid, was able to interfere with the CoQ10 pathway: At low concentration inhibiting the electron transport system, while at high concentration impairing the structural integrity of the mitochondrial membrane. Consistently, CoQ10 was found to affect the pharmacokinetic parameters of theophylline [27], an alkaloid of the xanthine family.
Remarkably, in the negative ionization mode we found different discriminant metabolites involved in the purine metabolism (Table 1). Inosine is a purine nucleoside formed by deamination of adenosine and, in the context of purine metabolism, is further transformed to hypoxanthine by the enzyme phosphate alpha-D-ribosyltransferase and then to xanthine and to uric acid. The fact that inosine was found upregulated in milk samples of sheep consuming cocoa husks may be associated to the inhibition of the enzyme phosphate alpha-D-ribosyltransferase, in order to reduce the increase of xanthines and uric acid derived from the metabolism of theobromine.
Our experimental findings are in accordance with scientific opinion of EFSA [10] on theobromine as an undesired substance in animal feed. In the EFSA report, theobromine was reported as showing moderate acute toxicity on dog and rodents, while in horses, which are particularly susceptible to theobromine, the liver and thyroid were affected, and pigs showed growth retardation, diarrhea and lethargy. Particular attention should be paid on the levels of theobromine and other alkaloids present in the industrial by-products in the preparing of the dairy sheep’s diet.

4. Conclusions

In this paper we show evidence, with an untargeted metabolomics approach, that sheep milk metabolite profiles are profoundly altered on the basis of an eight-week diet containing 50 and 100 g/d of cocoa husks. In particular these effects were more pronounced in the eighth week, at the end of the treatment. Moreover, this approach allowed us to highlight changes in multiple metabolic pathways. Importantly, the thyroid hormone L-thyroxine (T4) and tetraiodothyroacetate (T4A) were found down regulated for both CH50 and CH100 groups, thus suggesting a modification of the thyroid hormone metabolism. Furthermore, different ubiquinol-10 biosynthesis metabolites were found down regulated in milk samples of the CH50 and CH100. Biochemical mechanisms related to these changes remain unclear, while thyroid hormones levels should be further investigated in the ewes’ blood. In conclusion, taking into account these results, particular attention should be paid when using cocoa by-products in small ruminant feed management.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2624-862X/2/1/11/s1: Figure S1. Milk sample overlaid chromatograms of CH0, CH50, and CH100 groups acquired in positive (a) and negative (b) ionization mode. Figure S2. Temporal evolution of milk yield in the control (CH0) and in the two cocoa husks supplemented groups (CH50 and CH100). Figure S3. Partial least square discriminant analysis (PLS-DA) score plots of milk samples of the experimental groups (CH0, CH50, and CH100) collected in the last week of treatment. Lipid data collected in the LCMS positive (top, R2Y = 0.62, Q2Y = 0.27) and negative (bottom, R2Y = 0.98, Q2Y = 0.63) ionization mode. Figure S4. Temporal evolution of thyroxine (T4) milk levels in the control (CH0) and in the two cocoa husks supplemented groups (CH50 and CH100).

Author Contributions

Sample analysis, statistical analysis, wrote original draft, C.M. and P.C.; statistical analysis, wrote original draft, P.S.; conceptualization, wrote original draft, A.N. and G.P.; animal experiment and milk sample collection, S.C. All authors have read and agreed to the published version of the manuscript.

Funding

The research was partly funded by Mignini-Petrini company.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Ethics Committee of the University of Sassari (no. 54584/2018).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors are grateful to Stefano Mariotti and Elio Acquas for helpful discussions.

Conflicts of Interest

The authors declare that the sponsor had no role in the study design and interpretation of the results.

References

  1. Mirzaei-Aghsaghali, A.; Maheri-Sis, N. Nutritive value of some agro-industrial by-products for ruminants—A review. World J. Zool. 2008, 3, 40–46. [Google Scholar]
  2. Bampidis, V.A.; Robinson, P.H. Citrus by-products as ruminant feeds: A review. Anim. Feed Sci. Technol. 2006, 128, 175–217. [Google Scholar] [CrossRef]
  3. Yusof, F.; Khanahmadi, S.; Amid, A.; Mahmod, S.S. Cocoa pod husk, a new source of hydrolase enzymes for preparation of cross-linked enzyme aggregate. Springerplus 2016, 5, 1–18. [Google Scholar] [CrossRef] [Green Version]
  4. Cruz, G. Production of Activated Carbon from Cocoa (Theobroma cacao) Pod Husk. J. Civ. Environ. Eng. 2012, 02. [Google Scholar] [CrossRef]
  5. Carta, S.; Nudda, A.; Cappai, G. Cocoa husks replace effectively soybean hulls in dairy sheep diet: Effects on milk production traits and metabolic profile. J. Dairy Sci. 2019, 3–29. [Google Scholar]
  6. Sato, F. Plant Alkaloid Engineering, 3rd ed.; Elsevier Inc.: Amsterdam, The Netherlands, 2019; ISBN 9780124095472. [Google Scholar]
  7. Lelo, A.; Birkett, D.; Robson, R.; Miners, J. Comparative pharmacokinetics of caffeine and its primary demethylated metabolites paraxanthine, theobromine and theophylline in man. Br. J. Clin. Pharmacol. 1986, 22, 177–182. [Google Scholar] [CrossRef] [PubMed]
  8. Martínez-Pinilla, E.; Oñatibia-Astibia, A.; Franco, R. The relevance of theobromine for the beneficial effects of cocoa consumption. Front. Pharmacol. 2015, 6, 1–5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Adamafio, N.A. Theobromine Toxicity and Remediation of Cocoa By-product: An Overview. J. Biol. Sci. 2013, 13, 570–576. [Google Scholar] [CrossRef] [Green Version]
  10. Alexander, J.; Benford, D.; Cockburn, A.; Cravedi, J.; Dogliotti, E.; Domenico, A.D.; Férnandez-cruz, M.L.; Fürst, P.; Fink-gremmels, J.; Galli, C.L.; et al. Theobromine as undesirable substances in animal feed—Scientific opinion of the Panel on Contaminants in the Food Chain. EFSA J. 2008, 6, 1–66. [Google Scholar]
  11. Campos-Vega, R.; Nieto-Figueroa, K.H.; Oomah, B.D. Cocoa (Theobroma cacao L.) pod husk: Renewable source of bioactive compounds. Trends Food Sci. Technol. 2018, 81, 172–184. [Google Scholar] [CrossRef]
  12. Osorio, M.T.; Moloney, A.P.; Brennan, L.; Monahan, F.J. Authentication of beef production systems using a metabolomic-based approach. Animal 2012, 6, 167–172. [Google Scholar] [CrossRef]
  13. Xi, X.; Kwok, L.Y.; Wang, Y.; Ma, C.; Mi, Z.; Zhang, H. Ultra-performance liquid chromatography-quadrupole-time of flight mass spectrometry MSE-based untargeted milk metabolomics in dairy cows with subclinical or clinical mastitis. J. Dairy Sci. 2017, 100, 4884–4896. [Google Scholar] [CrossRef]
  14. Scano, P.; Carta, P.; Ibba, I.; Manis, C.; Caboni, P. An Untargeted Metabolomic Comparison of Milk Composition from Sheep Kept Under Different Grazing Systems. Dairy 2020, 1, 30–41. [Google Scholar] [CrossRef] [Green Version]
  15. Caboni, P.; Maxia, D.; Scano, P.; Addis, M.; Dedola, A.; Pes, M.; Murgia, A.; Casula, M.; Profumo, A.; Pirisi, A. A gas chromatography-mass spectrometry untargeted metabolomics approach to discriminate Fiore Sardo cheese produced from raw or thermized ovine milk. J. Dairy Sci. 2019, 102, 5005–5018. [Google Scholar] [CrossRef] [PubMed]
  16. Tautenhahn, R.; Patti, G.J.; Rinehart, D.; Siuzdak, G. XCMS online: A web-based platform to process untargeted metabolomic data. Anal. Chem. 2012, 84, 5035–5039. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Eriksson, L.; Johansson, E.; Antti, H.; Holmes, E. Multi- and Megavariate Data Analysis. In Metabonomics in Toxicity Assessment; CRC Press: Boca Raton, FL, USA, 2005; pp. 263–336. [Google Scholar]
  18. Macciotta, N.P.P.; Cappio-Borlino, A.; Pulina, G. Analysis of environmental effects on test day milk yields of Sarda dairy ewes. J. Dairy Sci. 1999, 82, 2212–2217. [Google Scholar] [CrossRef]
  19. Wu, S.Y.; Green, W.L.; Huang, W.S.; Hays, M.T.; Chopra, I.J. Alternate pathways of thyroid hormone metabolism. Thyroid 2005, 15, 943–958. [Google Scholar] [CrossRef] [PubMed]
  20. Moreno, M.; Kaptein, E.; Goglia, F.; Theo, T.J. Vesser Rapid Glucuronidation of Tri- and Tetraiodothyroacetic Acid to Ester Glucuronides in Human Liver and to Ether Glucuronides in Rat Liver. Endocrinology 1994, 135, 1004–1009. [Google Scholar] [CrossRef]
  21. Crossley, D.N.; Ramsden, D.B. Serum tetraiodothyroacetate (t,a) levels in normal healthy euthyroid individuals determined by gas chromatography-mass fragmentography (gc-mf). Clin. Chim. Acta 1978, 94, 267–272. [Google Scholar] [CrossRef]
  22. Todini, L. Thyroid hormones in small ruminants: Effects of endogenous, environmental and nutritional factors. Animal 2007, 1, 997–1008. [Google Scholar] [CrossRef] [Green Version]
  23. Valtorta, S.; Hahn, L.; Johnson, H.D. Effect of High Ambient Temperature (35°), and Feed Intake on Plasma T4 Levels in Sheep. Proc. Soc. Exp. Biology Med. 1982, 169, 260–265. [Google Scholar] [CrossRef] [PubMed]
  24. Silanikove, N. Effects of heat stress on the welfare of extensively managed domestic ruminants. Livest. Prod. Sci. 2000, 67, 1–18. [Google Scholar] [CrossRef]
  25. Wolff, J.; Varrone, S. The methyl xanthines—A new class of goitrogens. Endocrinology 1969, 85, 410–414. [Google Scholar] [CrossRef] [PubMed]
  26. Adachi, H.; Konishi, K.; Horikoschi, I. The effects of 1,2,3,4,6-penta-O-galloyl-beta-D-glucose on rat liver mitochondrial respiration. Chem. Pharm. Bull. 2002, 57, 364–370. [Google Scholar]
  27. Baskaran, R.; Shanmugam, S.; Nagayya-Sriraman, S.; Kim, J.H.; Jeong, T.C.; Yong, C.S.; Choi, H.G.; Yoo, B.K. The effect of coenzyme Q10 on the pharmacokinetic parameters of theophylline. Arch. Pharm. Res. 2008, 31, 938–944. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Principal component analysis (PCA) score plots of milk samples. Samples are colored by week of treatment (first, blue; second, red; third, green; and fourth, light blue circles). Metabolite data collected in the liquid chromatography-mass spectrometry (LCMS) (a) positive and (b) negative ionization modes.
Figure 1. Principal component analysis (PCA) score plots of milk samples. Samples are colored by week of treatment (first, blue; second, red; third, green; and fourth, light blue circles). Metabolite data collected in the liquid chromatography-mass spectrometry (LCMS) (a) positive and (b) negative ionization modes.
Dairy 02 00011 g001
Figure 2. Score plots of pairwise orthogonal partial least squares-discriminant analysis (OPLS-DA) of milk samples data in liquid chromatography-mass spectrometry positive ((a) CH0 vs. CH100 and (b) CH50 vs. CH0), and negative ((c) CH0 vs. CH100 and (d) CH50 vs. CH0) ionization modes. CH0 = green circles; CH50 = red boxes; CH100 = light blue triangles.
Figure 2. Score plots of pairwise orthogonal partial least squares-discriminant analysis (OPLS-DA) of milk samples data in liquid chromatography-mass spectrometry positive ((a) CH0 vs. CH100 and (b) CH50 vs. CH0), and negative ((c) CH0 vs. CH100 and (d) CH50 vs. CH0) ionization modes. CH0 = green circles; CH50 = red boxes; CH100 = light blue triangles.
Dairy 02 00011 g002
Figure 3. Box plots of (a) thyroxine (T4) and (b) tetraiodothyroacetate (T4A) levels in milk samples for the CH0, CH50, and CH100 groups in the last four weeks of treatments (week 1 = blue, 2 = red, 3 = green, 4 = cyan).
Figure 3. Box plots of (a) thyroxine (T4) and (b) tetraiodothyroacetate (T4A) levels in milk samples for the CH0, CH50, and CH100 groups in the last four weeks of treatments (week 1 = blue, 2 = red, 3 = green, 4 = cyan).
Dairy 02 00011 g003
Table 1. Orthogonal partial least squares-discriminant analysis (OPLS-DA) main discriminant metabolites as variable importance for prediction (VIP > 1) detected in positive ionization modes.
Table 1. Orthogonal partial least squares-discriminant analysis (OPLS-DA) main discriminant metabolites as variable importance for prediction (VIP > 1) detected in positive ionization modes.
Positive Ionization Mode
RT (min)m/zAdductMetaboliteVIPRegulationPathway Involved
1.14279.1603[M+H]+α-carboxyethyl-hydroxychroman2.21downVitamin E metabolism
6.57707.5197[M+K]+2-methoxy-6-(all-trans-octaprenyl)phenol1.93downUbiquinol-10 biosynthesis
5.77717.5652[M-H2O+H]+1,2-dipalmitoyl-phosphatidylcholine1.81down
1.26631.4714[M-H2O+H]+dipalmitoyl phosphatidate1.80down
0.64336.0713[M+H]+β-nicotinamide D-ribonucleotide1.55upNADH synthesis
0.61145.0504[M-H2O+H]+1,5-anhydro-D-fructose1.43up
0.64158.0415[M+2H]2+geranyl diphosphate1.40up
0.63118.0876[M-NH3+H]+L-ornithine1.26upUrea cycle
0.66321.0601[M-H2O+H]+5-amino-1-(5-phospho-D-ribosyl)imidazole-4-carboxamide1.15upPurine metabolism
8.69817.6885[M-H2O+H]+6-methoxy-3-methyl-2-all-trans-decaprenyl-1,4-benzoquinol1.07downUbiquinol-10 biosynthesis
Table 2. Orthogonal partial least squares-discriminant analysis (OPLS-DA) main discriminant metabolites (VIP>1) detected in negative ionization modes.
Table 2. Orthogonal partial least squares-discriminant analysis (OPLS-DA) main discriminant metabolites (VIP>1) detected in negative ionization modes.
Negative Ionization Mode
RT (min)m/zAdductMetaboliteVIPRegulation
0.62179.0572[M-H]theobromine2.43up
0.64211.0842[M-H]phenylalanine2.34up
6.17710.6315[M-2H2O-H]tetraiodothyroacetate1.75downThyroid hormone synthesis
6.48833.6484[M-H]3,4-dihydroxy-5-all-trans-decaprenylbenzoate1.61downUbiquinol-10 biosynthesis
6.15812.6624[M+Cl]L-thyroxine1.60downThyroid hormone synthesis
11.38369.3014[M+CH3COO]phytenate1.58downPhytol degradation
0.72815.0952[M+CH3COO]5′,5′′′-diadenosine triphosphate1.49upPurine metabolism
0.65259.0235[M-H]hexose 6-phosphate1.40up
0.64478.9864[M+Na-2H]8-oxo-GDP1.36up
0.61267.0740[M-H]inosine1.20upPurine metabolism
7.33311.2968[M-H]arachidic acid1.08up
11.41489.3955[M+CH3COO]α-tocopherol1.03down
0.66142.0214[M-2H2O-H]L-cysteinyl-glycine1.00up
1.58435.2918[M+CH3COO]all-trans/10′-apo-β-carotenal1.00upCarotenoid cleavage
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Manis, C.; Scano, P.; Nudda, A.; Carta, S.; Pulina, G.; Caboni, P. LC-QTOF/MS Untargeted Metabolomics of Sheep Milk under Cocoa Husks Enriched Diet. Dairy 2021, 2, 112-121. https://0-doi-org.brum.beds.ac.uk/10.3390/dairy2010011

AMA Style

Manis C, Scano P, Nudda A, Carta S, Pulina G, Caboni P. LC-QTOF/MS Untargeted Metabolomics of Sheep Milk under Cocoa Husks Enriched Diet. Dairy. 2021; 2(1):112-121. https://0-doi-org.brum.beds.ac.uk/10.3390/dairy2010011

Chicago/Turabian Style

Manis, Cristina, Paola Scano, Anna Nudda, Silvia Carta, Giuseppe Pulina, and Pierluigi Caboni. 2021. "LC-QTOF/MS Untargeted Metabolomics of Sheep Milk under Cocoa Husks Enriched Diet" Dairy 2, no. 1: 112-121. https://0-doi-org.brum.beds.ac.uk/10.3390/dairy2010011

Article Metrics

Back to TopTop