Ultra high performance liquid chromatography–high resolution mass spectrometry plasma lipidomics can distinguish between canine breeds despite uncontrolled environmental variability and non-standardized diets

Introduction and objectives The purpose of this study was to use high accurate mass metabolomic profiling to investigate differences within a phenotypically diverse canine population, with breed-related morphological, physiological and behavioural differences. Previously, using a broad metabolite fingerprinting approach, lipids appear to dominate inter- and intra- breed discrimination. The purpose here was to use Ultra High Performance Liquid Chromatography–High Resolution Mass Spectrometry (UHPLC–HRMS) to identify in more detail, inter-breed signatures in plasma lipidomic profiles of home-based, client-owned dogs maintained on different diets and fed according to their owners’ feeding regimens. Methods Nine dog breeds were recruited in this study (Beagle, Chihuahua, Cocker Spaniel, Dachshund, Golden Retriever, Greyhound, German Shepherd, Labrador Retriever and Maltese: 7–12 dogs per breed). Metabolite profiling on a MTBE lipid extract of fasted plasma was performed using UHPLC-HRMS. Results Multivariate modelling and classification indicated that the main source of lipidome variance was between the three breeds Chihuahua, Dachshund and Greyhound and the other six breeds, however some intra-breed variance was evident in Labrador Retrievers. Metabolites associated with dietary intake impacted on breed-associated variance and following filtering of these signals out of the data-set unique inter-breed lipidome differences for Chihuahua, Golden Retriever and Greyhound were identified. Conclusion By using a phenotypically diverse home-based canine population, we were able to show that high accurate mass lipidomics can enable identification of metabolites in the first pass plasma profile, capturing distinct metabolomic variability associated with genetic differences, despite environmental and dietary variability. Electronic supplementary material The online version of this article (doi:10.1007/s11306-016-1152-0) contains supplementary material, which is available to authorized users.


3
15 Page 2 of 11 home-based dogs resulted in less distinct breed discrimination (Beckmann et al. 2010;Lloyd et al. 2016). In a more phenotypically diverse and 'free-living' canine population, despite the diet and environmental confounders, some breeds remained distinct from all others (Greyhound, Chihuahua and Dachshund). Furthermore in the two breeds that were present in both the controlled environmental and the client-owned environmental study cohorts (Labrador Retriever and Cocker Spaniel) phosphatidylcholines were identified common to dogs of the same breeds but of different genetic stock and from different environmental conditions (Lloyd et al. 2016). It was speculated that other differences in lipid metabolism and lipid functionality may exist between breeds.
Breed differences in lipids have been reported previously. In a targeted metabolite profiling study investigating individual variability when fed the same diet, metabolites that showed high inter-individual and low intra-individual variance were also related to lipid/sterol chemistry (Colyer et al. 2011). These confirm other non-metabolomics studies: for example total plasma cholesterol and triacylglyceride concentrations have been shown to vary between canine breeds and with diet (Downs et al. 1993;Usui et al. 2014) and Beagles and Greyhounds have been shown to have different pulmonary surfactants (Clercx et al. 1989). Primary hyperlipidemia is usually associated with certain breeds such as Miniature Schnauzers (Xenoulis et al. 2013) and significant differences have been shown between lipoprotein density profiles of different breeds (Downs et al. 1993;Xenoulis et al. 2013), with and without hyperlipidemia (Xenoulis et al. 2013). Additionally, within the same breed of dog, the lipoprotein profile and plasma lipids can be influenced by dietary fat, activity level and lifestyle (Downs et al. 1997). These data from the literature and the consistent finding in our previous studies show that lipids dominate the metabolites that discriminate between breeds and individuals within a phenotypically diverse dog population (Lloyd et al. 2016). In the current study we present more systematic and detailed characterization of the plasma lipidome in dog breeds using Ultra High Performance Liquid Chromatography-High Resolution Mass Spectrometry and suggest a strategy for working with owner-maintained animals in the context of uncontrolled environmental variability and non-standardized diets.

Animal maintenance
The samples are derived from the previously reported study (Lloyd et al. 2016). The client-owned dogs were recruited in Australia (University of New England, UNE) and maintained on different diets and fed according to their owners' feeding regimen and amounts (see electronic supplementary Table S1 for full diet information). Exclusion criteria were dogs with obvious signs of ill health, pregnant or lactating dogs, bitches in season and dogs considered to be puppy or senior for the breed (veterinarian's discretion). The breeds represented in the present study were Beagle (Be, n = 12), Chihuahua (Ch, n = 7), Cocker Spaniel (CS, n = 12), Dachshund (Da, n = 8), Golden Retriever (GR, n = 12), Greyhound (Gh, n = 12), German Shepherd (GS, n = 12), Labrador Retriever (LR, n = 12) and Maltese (Ma, n = 9).

Plasma collection and extraction
A fasted (>12 h fast) blood sample was taken from the jugular vein by trained veterinarians from each individual dog. Blood was collected for both RNA and metabolite fingerprinting analysis. In keeping with reduction and refinement principles of research with animals, the sample volume required was minimised by filtering blood to collect leukocytes for RNA extraction prior to centrifugation of the leukocyte-depleted blood filtrate to obtain plasma. Fasting blood samples (up to 9 ml), collected in EDTA tubes (BD Vacutainer, 10 ml 367,525), were mixed by inversion and passed through a Leukolock filter using a 25G needle into an evacuated plain tube (BD Vacutainer, 10 ml 366,636) on ice. The leukocyte-depleted blood sample was kept on ice until centrifuged (2000 g for 15 min). Plasma samples were collected and stored, either on dry ice for transport or for long-term storage at −80 °C until analysis.

Lipid extraction of plasma
Plasma was extracted for metabolomic analysis using the following method. Aliquots of chilled plasma (20 μl) on ice were diluted with 10 µl of pre-chilled water (on ice) and 150 µl of pre-chilled methanol (−20 °C), vortexed and then 500 µl MTBE (methyl-tert-butyl ether) was added (Matyash et al. 2008). Samples were then shaken for 60 min at room temperature and 125 µl of water was added before they were centrifuged for 1 min at 14,000 g (room temperature). The organic phase (400 µl) was transferred into a new Eppendorf tube and stored at −80 °C until required. Immediately prior to use, the samples were centrifuged for 3 min at 14, 000 g (room temperature).

Ultra high performance liquid chromatographyhigh resolution mass spectrometry (UHPLC-HRMS) analysis
Extracted plasma samples were analysed on an Exactive Orbitrap (Thermo Fisher Scientific) mass spectrometer, which was coupled to an Accela Ultra High Performance Liquid Chromatography (UHPLC) system (Thermo Fisher Scientific). Chromatographic separation was performed on a reverse phase (RP) Acquity UPLC BEH C 18 1.7 µm, 2.1 × 100 mm column (Waters, Milford, USA) using ACN: H 2 O (40:60) with 10 mM ammonium acetate as mobile phase solvent A and ACN: isopropanol (10:90) with 10 mM ammonium acetate as mobile phase solvent B. Each sample (20 μl) was analysed using 30-100% gradient of B in 10 min after washing the column with 100% of B for 2 min and re-equilibrating the column with 70: 30 of A: B for 3.5 min. A flow rate of 400 µl/min was used for running the samples. Column oven temperature was set to 60 °C and the data were acquired in both positive and negative ESI, using a heated electospray ionisation source (HESI). Mass spectra were acquired from 70 to 1400 mass-to-charge ratio (m/z) using a mass resolution of 100,000. The spray voltage was 4 kV for both ionization modes. The temperature of the ion transfer capillary was 370 °C and sheath and auxiliary gas was 30 and 15 arbitrary units, respectively. The data were recorded using the Xcalibur 2.0.0 software package (Thermo Fisher Scientific). Mass calibration was performed for both ESI polarities before the analysis using a mixture of caffeine, MRFA (L-methionyl-arginyl-phenylalanyl-alanine), Ultramark 1621, sodium dodecyl sulfate (SDS), and sodium taurocholate dissolved in acetonitrilemethanol-water solution with 1% acetic acid.

UHPLC-HRMS data processing
All data were acquired in profile mode, and converted to centroid mzXML using msconvert (ProteoWizard). All processing was performed using XCMS (Smith et al. 2006) (version 1.3.8) and in R (version 3.0.2). XCMS was used to process raw LC-MS data into a matrix of ~2500 features with corresponding raw intensity values.

Sample classification and selection of potentially explanatory signals
FIEmspro was used for all multivariate modelling, classification and feature selection Enot et al. 2008) (URL http://users.aber.ac.uk/jhd/). High accurate mass data was binned to 0.01 amu after raw data processing using bespoke routines, allowing direct comparison of plasma profiles, prior to signal annotation. Principal Component Analysis (PCA) was followed by PC-Linear Discriminant Analysis (PC-LDA). Plots of the first two PC-Discriminant Functions (PC-DFs) allowed visualization of the goodness of class separation as quantified by Tw values (Eigenvalues). Random Forest (RF) was employed in the analysis of the multivariate data and the RF classification 'margin', along with the area under the ROC (receiver operating characteristic) curve (AUC) and accuracy (ACC) was used to assess classification performance . Models were deemed adequate overall if RF margins > 0.2 and AUC and/or ACC values > 0.8, thresholds which we have implemented in previous publications . A combination of RF, AUC and student's t-test, were used to highlight potentially explanatory signals responsible for discriminating between sample classes in a full feature rank list . RF feature selection was performed by calculating Importance Scores (RFIS), being the mean decrease in accuracy over all classes when a feature is omitted from the data. AUC used the area under curve of the sensitivity (true-positive rate) against the specificity (false-positive rate) and student's t-test ranked the features by the absolute value of the P-values.
Randomized re-sampling strategies using bootstrapping were applied in the process of classification and feature selection to counteract the effect of any unknown, structured variance in the data. We used 100 bootstraps in pairwise comparisons for each of the applied statistical operations with 2/3 of data as training and 1/3 as test set. RF was set to ntree = 1000 for each bootstrap which is adequate considering the dimensionality of data.
Pearson correlation coefficients between selected variables were calculated using the function cor in R version 2.5.1 . Variables with correlation coefficient absolute values (|r|) > 0.8 were considered to belong to a cluster indicative of different ionization or potential biotransformation/breakdown products of a single metabolite.

Targeted flow infusion electrospray-ionization tandem mass spectrometry and annotation of signals
For metabolite signal annotation, accurate m/z values were extracted from the un-binned matrix to enable direct identification of metabolites at 1-5 ppm directly in the first pass profile. These were queried using MZedDB, an interactive accurate mass annotation tool which can be used directly to annotate signals by means of neutral loss and/or adduct formation rules (Draper et al. 2009). Selected ion signals were investigated further using tandem Mass Spectrometry (MS n ) on a Nano-Flow (TriVer-saNanoMate, AdvionBioSciences Ltd, UK) LTQ-Fourier Transform-Ion Cyclotron Resonance Ultra-Mass-Spectrometry (FT-ICR-MS; where Ultra refers to the high-sensitivity ICR-cell) as reported previously (Lloyd et al. 2011(Lloyd et al. , 2016. Samples were prepared as for UHPLC-HRMS and reconstituted in methanol/water (80/20, v/v). For each targeted m/z value a scan window was set for 20 scans, an isolation width of 1 m/z and using normalized collision energy of 30 V. An activation coefficient 'Q' of 0.250 was chosen and an activation time of 30 ms, with wideband activation turned on and a source fragmentation of 20-30 V. Mass range settings were dependent upon the molecular weight of the target ion.
Metabolites were putatively annotated to MSI level 2 without chemical reference standards due to the lack of standard availability, based upon physicochemical properties, retention times and spectral similarity with public/ commercial spectral libraries [Lipid Maps, HMDB, Metlin and Massbank (Horai et al. 2010;Sana et al. 2008;Sud et al. 2007;Wishart et al. 2009)]. Choline-containing phospholipid species, Phosphatidylcholine (PCs) and Sphingomyelin (SMs), both show a characteristic m/z 184 phosphocholine head group peak, as well as an [M + H-59] + peak corresponding to the neutral loss of (CH 3 ) 3 N. Whereas fragmentation of phosphatidylethanolamines exclusively yields one peak, an [M + H-141] + ion from the neutral loss of the phosphoethanolamine head group and fragmentation of phosphatidylserines yields an [M + H-185] + ion from the neutral loss of the phosphatidylserine head group (Pi et al. 2016).

Dog demographics
The full demographics of the study dogs are shown in electronic supplementary Table S1. The mean (±SEM) age of the 96 dogs was 56 months (±3) and had a range of 8-120 months. Using a 9 point body condition score (Laflamme 1997), the mean was 3.0 with a range of 2.0-4.5. The weight varied considerably between breeds (1.8-44.2 kg) with a mean (±SEM) of 20.6 (±0.9 kg). Due to logistical constraints and the focus on the study being breed, recruitment did not result in balanced age and gender within each breed. More female dogs (n = 58) participated in the investigation than males (n = 37), and there was one unrecorded gender; however, there was representation of both females and males in each dog breed class.

Canine breed clusters can be differentiated in the plasma lipidome
Previously, we used non-targeted nominal mass fingerprinting to characterize the main drivers of variance in the plasma metabolomes of dogs (Lloyd et al. 2016). The primary fingerprint data consisted of approximately 800 signals at 1 amu in each ionisation mode. In this current work, after pre-processing, data resulted in approximately 2500 signals per profile (with high sensitivity of 1-5 ppm identification level). These data were  Fig. 1a). LR samples were distributed across the first dimension showing intra-breed variation. Separation of breeds into two clusters was also observed in PC-DF2 in which Ch, Da, Ma and GR were distinct from the other five breeds (Tw 2.55, Fig. 1a).
The Tw values calculated for the first two dimensions of variance indicated clearly that the plasma lipidome profile of dogs in a fasting state was strongly influenced by breed. As with previous observations using non-targeted metabolome fingerprinting (Lloyd et al. 2016) both breed and breed size appeared significant factors in PC-DF1.

Identification of the main signals associated with the plasma lipidome of canine breed meta-classes
Prior to investigating possible unique inter-breed plasma lipidome differences, we elected to further characterise the signals causing breed cluster discrimination on PC-DF1 and PC-DF2 seen in Fig. 1a. PCA and PC-LDA models were determined with new meta-classes explaining PC-DF1 and PC-DF2. The first meta-class structure was to explain the variance on PC-DF1 [Group 1 (Ch, Da and Gh) vs Group 2 (Be, CS, GR, GS and Ma)], with LR data omitted due to the lack of distinctive clustering with either of the two groups. The second meta-class designations were structured to elucidate the differences on PC-DF2: Ch, Da, Ma and GR were placed into one meta-class and the other breeds including LR were placed in the comparison metaclass. RF margins, AUC values and ACC modelling scores for these models are shown in electronic supplementary Table S2. Previously, to characterize the main drivers of variance in the plasma metabolomes of dogs, we used metabolite fingerprinting with targeted high accurate mass analysis, with information acquired on an LTQ-FT-Ultra for metabolite identification, followed by FIE-MS n for structual confirmation. In this current work, m/z data were generated at high resolution and sufficient accuracy to separate signals contained previously in the same 1 amu mass bin, as well as enabling direct identification of metabolites directly in the first pass profile without the need for further targetted accurate mass analysis, only FIE-MS n for structual confirmation. Coupling of MS to LC additionally allowed retention times to be associated with every accurate m/z enabling structural isomers, adducts and isobaric compounds to be distinguished and ionisation adducts to be confirmed.
The explanatory signals associated with the discrimination between the meta-classes in Fig. 1a (along PC-DF1 and PC-DF2) were investigated further as the lipidome differences were above adequate thresholds [i.e. RFIS > 0.003 and false discovery rate (FDR) corrected P-values < 0.05; ]. The top ranked cluster of correlated (|r| > 0.8) signals driving separation on PC-DF1 between Ch, Da and Gh and the rest of the breeds is displayed in Table 1. The remainder of features above the implemented thresholds primarily responsible for the discrimination on PC-DF1 are shown in electronic supplementary Table S2, where correlated signals (|r| >0.8) are highlighted in different colours to show association. These signals were identified as a range of triglycerides, diglycerides and phosphatidylcholines which appeared much reduced in intensity in the Ch, Da and Gh breeds compared with the other breeds (Fig. 2a). Triacylglycerides are the most common and efficient form of stored energy in mammals and higher plasma levels can reflect differences in dietary intake but also hepatic health (Ginsberg 1998;Rifai et al. 1999).
Metabolite identifications driving the separation on PC-DF2 were also assigned using the same strategy (Table 2 and electronic supplementary Table S4). These metabolites appeared reduced in intensity in the Ch, Da, Ma and GR compared with the other breeds (Fig. 2b) and represented different chemical classes to those causing the   Table S4).

Diet is one of the main drivers associated with the intra-and inter-variation of plasma lipidomes of different canine breeds
By analysing the meta-data and the habitual eating patterns (see electronic supplementary Table S1) associated with each dog we identified that the discrimination seen in PC-DF1 in Fig. 1a also was impacted by diet types. This is supported by the inter-individual variation seen in LR, as those were distributed in line with whether the dogs followed a diet containing chicken meat (group1; on the negative side of PC-DF1) or red meat (group 2; on the positive side of PC-DF1) as their main protein source. We suggest that the metabolites seen at higher concentrations in group 2 (shown in Table 1 and electronic supplementary  Table S3) compared with group 1 correlate with higher red meat consumption by group 2. Raw red meat including beef, lamb and mutton has been shown to contain increased total saturated and monounsaturated fatty acids when compared with white meat such as chicken (Droulez et al. 2006;Williams 2007). Additionally concentrations of total glycerophosphocholines and phosphatidylcholines have also been shown to be much higher in pork/ beef muscle meat and sausages/burgers when compared with chicken muscle meat and processed meat (Zeisel et al. 2003). In human plasma total saturated and monounsaturated fatty acid concentrations increased after a diet high in red meat when compared with a diet high in fatty fish (Wolmarans et al. 1991). Triacylglycerols, derived from glycerol and three fatty acids (both saturated and monounsaturated), have also been shown to increase with high habitual consumption of red meat (Wolmarans et al. 1991) in human plasma. It is important to note that all plasma samples were derived from fasted animals and thus as domesticated dogs may be maintained for many years on the same daily diet it is suggested that any studies on endogenous metabolism should ensure that confounding diet-derived signals are dealt with appropriately.

Breed cluster differences are enhanced when confounding diet signals are removed from the plasma lipidome of dogs consuming varied diets
We have previously used a data filtering approach to allow the identification of metabolite signals associated specifically with exposure to test foods within a mix of 'standardized food items' in a human intervention study (Lloyd et al. 2011). The filtered list included all those explanatory signals with a RFIS > 0.002 after the consumption of the standardised food items when compared with fasting urines. A similar approach was employed here. From the breed model on PC-DF1 (Fig. 1a) the top RFIS ranked signals causing the differences between the binary comparison of the two classes (Ch, Da and Gh vs the other breeds in Table 1), were filtered and removed as they were assumed to be responsible for the confounding diet effect. We incrementally increased the filtering cut-off until the PC-LDA model (Fig. 1a) collapsed and the Ch, Da, Gh vs Be, CS, GR, GS, Ma clustering along PC-DF1 was no longer predominant. This heuristic RFIS threshold was >0.0003; additionally, features which were below the critical threshold, yet showed high correlation with the filtered top ranked signals with a |r| > 0.8, were also omitted from subsequent analyses as these have been shown to be isotopes and adducts in previous analyses. The resultant matrix contained 21.7% fewer features (539 signals out of 2496 were removed) than the complete matrix and when used to model the breed classes the main source of variance (PC-DF1) was now between four breeds (Ch, Da, Ma and GR) and the other five breeds (Tw 2.65, Fig. 1b), as observed previously in PC-DF2 in the original un-filtered model (see Fig. 1a). Notably the LR samples now formed a single distinctive cluster, suggesting the differences in diet were responsible for much of the variance in fasting plasma lipidome. Furthermore, with the exception of GR, the main axis of discrimination (PC-DF1) may be related to breed average weight (Table 1 and electronic supplementary Table S1) as observed previously (Beckmann et al. 2010) in the urine metabolome analysis. The top ranked signals discriminating the Ch, Da, Ma and GR from the other breeds of dogs were identical to those identified previously (Table 2), suggesting shared metabolome similarities in these four breeds. The sphingolipid C 17 sphinganine and the long chain hydrocarbons with amine groups were lower in the Ch, Da, Ma and GR breeds. Labrador was one of the breeds that had increased concentrations of these  Table 2); c chihuahua (m/z 429.3715 in Table 4); d golden retriever (m/z 705.5889 in Table 4); e greyhound (m/z 776.5786 in Table 4). Where the box indicates the interquartile range; the horizontal bar, the median; error bars represent the range of 7-12 dogs. Be beagle, Ch chihuahua, CS cocker spaniel, Da dachshund, GR golden retriever, Gh greyhound, GS german shepherd, LR labrador retriever, Ma maltese ◂ discriminatory sphingolipids and related amines when compared with Ch, Da, Ma and GR breeds. Differences in sphingolipid concentrations evident between the breeds in this current study may reflect how the dog breeds differentially transport or mobilise plasma sphingolipids, or may indicate different responses in the insulin signalling pathway between canine breeds (reviewed in Ilan 2016). Also, it has been suggested that endogenous sphinganine may inhibit cholesterol transport in Niemann-Pick Type C (NPC) disease (Roff et al. 1991) providing a suggestion that these dog breeds may have differential cholesterol esterification abilities.

Unique breed-specific signals are evident when diet-associated signals are filtered from plasma lipidome dataset
One of the objectives of the study was to determine if there were any unique breed-specific signals (i.e. a signal that is not part of a common series of signals discriminating several breeds). Using the diet-filtered dataset, the modelling scores (RF margins, AUC and ACC values) for each breed pair-wise comparison are shown in Table 3. Adequate discrimination modelling scores were shown by CS and LR for only four pair-wise breed comparisons: therefore, we did not pursue these two breeds further (Table 3). For the remaining seven breeds (Ma, GS, Be, Ch, GR, Gh and Da) the modelling scores for ≥6 pair-wise breed comparisons were above the adequate thresholds criteria set, and we sought to identify the explanatory signals specific to each breed. By comparing top ranked signals (RFIS > 0.002 and P-values < 0.05) for each breed comparison, none of the Ma, GS, Da nor Be plasma lipidomes showed any unique breed-specific features consistent across all comparisons. In contrast, consistent unique breed specific features were discovered for the Ch (5 decreased in signal intensity compared to other breeds, example in Fig. 2c), Golden Retriever (7 elevated, example in Fig. 2d) and Greyhound (1 decreased, Fig. 2e). The full RF-IS and P-value data for all pair-wise comparisons are shown in electronic supplementary Table S5. It was interesting to note that despite the absence of any unique breed differences the Beagle plasma lipidome still discriminated from many other breeds suggesting a broad variance in the Beagle plasma lipidome, which is consistent with previous observations in relation to the urinary metabolome of this breed (Beckmann et al. 2010). Direct annotation of signals unique to each of the three breeds, Ch, GR and Gh using MZedDB in the first pass profile was made without the need for further accurate mass analysis (Table 4). We have recently shown that by using a non-targeted 'global' metabolomic approach, Ch and Gh were discriminated from other breeds in the study (Lloyd et al. 2016). In Ch the plasma metabolites that were lower than in other breeds appeared to be vitamin D3 sterol lipids and phosphatidylcholines, which is in agreement with a previous report that the size of a dog appears to influence vitamin D3 metabolism (Tryfonidou et al. 2003 (Usui et al. 2014). Cholesterol is the main sterol in animal tissues which plays a fundamental role in central metabolic pathways including steroid hormone and vitamin D synthesis (Ginsberg 1998;Rifai et al. 1999). Phosphatidylcholines are predominantly structural lipids important in plasma membranes and lung surfactants as well as being major components in lipoproteins. Interbreed differences in lipoproteins have been reported previously, with the Cairn Terrier having particularly low levels of Low Density Lipoproteins compared to the LR, Be, Dh and other breeds (Downs et al. 1993). The explanatory signal that was at a lower concentration in Greyhound plasma compared with other breeds could not be identified. As Gh generally are seen to have a lower percent body fat than other breeds (Jeusette et al. 2010) and a high percent muscle mass (Drost et al. 2006), it is speculated that this single explanatory signal may be related to the unusual body conformation of the Gh and a robust marker of the Gh breed. The signals elevated in intensity in the Golden Retriever plasma appeared to be sphingolipids, sphingomyelins and phosphatidylcholines. This observation may be consistent with previous evidence of the role of sphingolipids in skin disorders in Golden Retrievers, which are prone to atopic dermatitis (Shaw et al. 2004). Healthy dogs of mixed breeds have been shown to have erythrocytes that contain increased concentrations both of saturated sphingomyelin and also some phosphatidylcholines when compared to dogs suffering with atopic dermatitis (Fuhrmann et al. 2006). Additionally, there is evidence that sphingomyelin undergoes significant interactions with cholesterol (García-Arribas et al. 2016;Ridgway 2000) and Golden Retrievers have been shown to have the highest total cholesterol concentration out of 51 dog breeds studied and a top 14 concentration of total triacylglycerides (Usui et al. 2014).

Concluding remarks
High accurate mass lipidomic profiling generated data at sufficient accuracy to enable direct identification of metabolites in the first pass profile, without the need for further targeted accurate mass analysis. Within a phenotypicallydiverse, home-based canine population, plasma lipidomics data indicated unique breed-specific signals can be discovered and improved when diet-confounding signals were removed. Despite constraints of unbalanced populations for age and gender within the set, robust lipids associated with unique breeds were postulated. The present data, together with our observations in previous studies (Lloyd et al. 2016), provide further evidence that individuals of closed genetic breed groups can be consistently discriminated from others across different metabolomic platforms and also stress the importance of extensive metadata to evaluate metabolome models. Overall, these data confirm the suggestion that breed development in domesticated dogs has resulted in distinctive, stable metabolic biotypes which, when understood in more depth, may allow improved understanding of healthcare requirements in individual animals. ing, welfare and sampling of dogs used in the study. In particular, we would like to thank Wendy Brown who co-ordinated the sampling and data collection from those dogs. We would also like to thank all of the clients who allowed blood samples to be used in this study. The authors' contributions to the work were as follows: AJL, data analysis, produced figures, researched literature and wrote the manuscript; MB, developed plasma extraction procedures, supervised mass spectrometry support staff, pre-processed data for analysis and edited the manuscript; TW performed UHPLC-HRMS analysis, data processing and edited the manuscript; KT provided metabolomics technical support. DA, devised and planned the studies, co-ordinated plasma sampling and metadata collection and wrote the manuscript. JD, supervised research in Aberystwyth, designed figures and wrote the manuscript.
Funding Funding was provided by WALTHAM Centre for Pet Nutrition, a division of Mars Petcare. DA is employed by the WALTHAM Centre for Pet Nutrition.

Compliance with ethical standards
Conflict of interest The work was funded by WALTHAM Centre for Pet Nutrition, a division of Mars Petcare. DA is employed by the WALTHAM Centre for Pet Nutrition.

Compliance with ethical requirements
All studies were in accordance with the research ethical standards at UNE. Informed consent was obtained from all dog owners.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Table 4 Identification of signals explanatory of Chihuahua (Ch), Golden Retriever (GR) and Greyhound (Gh) UHPLC-HRMS ultra high performance liquid chromatography-high resolution mass spectrometry analysis, PC phosphatidylcholine, SM sphingomyelin *Flow Infusion Electrospray-ionization tandem Mass Spectrometry (FIE-MS n ) was used to confirm annotations due to the distinctive fragmentation patterns. Choline-containing phospholipid species, PCs and SMs, both show a characteristic m/z 184 phosphocholine head group peak, as well as an [M + H-59] + peak corresponding to the neutral loss of (CH 3 ) 3 N