The Erythrocyte Membrane Lipidome of Healthy Dogs: Creating a Benchmark of Fatty Acid Distribution and Interval Values

Introduction

Lipidomics is the discipline that gathers lipids, not only considering lipid structures and their transformations, but also providing a “dynamic” interpretation of lipid diversity and functions in view of cellular, metabolic, and environmental conditions influencing living organisms (1). Investigations in lipidomics encompass a broad range of topics, going from physiological to oxidative, and pathological processes (2, 3), and large sets of data can be acquired, obtaining the map of the entire lipidome and its variations as a fingerprint of health, age, or sex status (4, 5). The development of shotgun lipidomics accelerated the process of gathering information on thousands of lipid molecules in different organisms and conditions, going from yeast, virus, and bacteria to murine models and humans (6–9). In particular, membrane lipidome has a central place in a multidisciplinary context of biophysical, biological, pharmacological, and clinical studies, which testify to its importance in fundamental processes such as cell formation, regeneration, and metabolic regulation (10–12). Focusing on membrane glycerophospholipids, the analysis of the two fatty acid chains esterified to the glycerol moiety led to the discovery of positional and geometrical isomers and of the homeostatic balance created by these hydrophobic molecules, based on the fact that the fatty acid structures influence bilayer thickness, permeability, fluidity, and membrane protein functions (10–14). In Figure 1, the main fatty acid types are shown, that can be saturated (SFA), monounsaturated (MUFA), and polyunsaturated fatty acids (PUFA). Application of fatty acid-based membrane lipidomics to human health is a hot topic of research, comparing healthy controls with disease subjects (15–18). Fatty acid analysis can be carried out starting from a non-invasive blood withdrawal, separating blood fractions that contain molecular information at different metabolic levels (18, 19). In fact, while plasma informs on the daily fatty acid dietary intakes, cell membrane fatty acids result from more complex contributions of physiological and metabolic processes together with nutritional intakes. Erythrocytes are the most numerous cell type in blood (20), and their membrane homeostasis derives from the distribution of the various fatty acid types, playing roles in several functions, such as for example oxygen distribution. Having lost DNA very early, red blood cell (RBC) fatty acids cannot come from biosynthesis but derive from exchange with lipoproteins and tissues. In this scenario, it is worth mentioning that eukaryotic cells need essential fatty acids (EFA), such as omega-3 and omega-6 PUFA, that must be taken daily from the diet in order to be processed enzymatically in the body and provide long-chain PUFAs. The PUFA pathways, together with SFA and MUFA, contribute to the formation of the cellular lipid pool and subsequently to the specific fatty acid biodistribution in RBC membranes as well as in the various organs. The RBC fatty acid composition is informative of levels reached in tissues, like muscle, and liver, and also in tissues not withdrawable during life, like retina, and brain (21–24). In human studies, we contributed to research in membrane fatty-acid-based lipidomics, making the arbitrary choice of analyzing a specific fatty acid cohort of erythrocyte membrane glycerophospholipids, made of 10 fatty acids representative of SFA, MUFA, and PUFA families. Their interval values in healthy cohorts were reported in membrane fatty acids by several studies, also by us (25–27), and in one of the most complete meta-analysis appeared in the literature (28). Taking into account that novel fatty acid pathways and transformations are continuously studied (29), it is important to make the strategic choice of a specific number of fatty acids considering their fundamental biochemical/biological roles and the consensus reached about their levels in cell membranes.

FIGURE 1

Figure 1. (A) The structure of fatty acids considered in the present work. Together with the common names, the abbreviations describe the position and geometry of the double bonds (e.g., 9c for palmitoleic acid), as well as the notation of the carbon chain length and total number of double bonds (e.g., C18:1). For PUFA in the parentheses, their acronyms are indicated. (B) Molecular structures of glycerophospholipids with the two fatty acid residues at the C-1 and C-2 positions of the glycerol moiety (in red), whereas in the C-3 position the polar head group is connected (examples of polar heads: ethanolamine, serine, and choline).

By these premises, we considered whether analogous data are reported in domestic animals, finding only a few studies on plasma or other tissues (brain, sperm) and on correlations between dietary uptake and lipidomic profiles of dogs (30–34). Limited studies presented a small number of samples (35, 36), focusing for example on PUFA (omega-6 and omega-3), and not deepening the different fatty acid families, which are instead important for molecular contribution to membrane homeostasis (37). Although fatty acids are recognized to be crucial for all animals (38), an effort to build up the methodology for fatty-acid-based membrane lipidomics in the veterinary field is missing. In the present work, we evaluated the composition of fatty acids in erythrocyte membrane phospholipids in a group of clinically healthy dogs (n = 68), focusing on the interval ranges of the 10 fatty acid cohort previously used for humans. We considered also correlations with dog characteristics such as age and bodyweight. The general aim of this work is to create a benchmark of fatty acid interval values in the membrane lipidome of healthy animals, useful to start a systematic approach for the examination of metabolic and nutritional status in healthy and diseased dogs.

Materials and Methods

For the list of materials and providers see Supplementary Information.

Inclusion Criteria and Samples Collection

Inclusion criteria for the study were: healthy dogs without any clinical or pathological evidence of disease accordingly to unremarkable history, physical examination and results of CBC, serum biochemistry, and voided urine analysis (Supplementary Table 1) within the reference range for each dog. Furthermore, recruited dogs were on commercial diet and did not receive supplements and medications in the previous 4 months (except regular preventive treatments for ecto- and endoparasites and prophylactic vaccination), as ascertained by an interview with the owners. Supplementary Table 2 presents the main dog characteristics. Blood samples (1 mL each) from healthy dogs recruited from Italian owners were collected in Ethylenediaminetetraacetic acid as tripotassium salt (EDTA) tubes from the medical staff of the Veterinary Teaching Hospital (VTH) of Teramo.

Ethical Statement

The project has been approved by the Health Ministry and the Committee on Animal Research and Ethics of the Universities of Chieti-Pescara, Teramo and Experimental Zooprophylactic Institute of AeM (CEISA), Protocol UNICHD12 n. 1168.

Isolation of Phospholipids From Erythrocyte Membranes and Transesterification Procedure

The steps of the protocol procedure are depicted in Scheme 1. Step 1: separation of red blood cells from plasma, effected on 1 mL of fresh EDTA-treated whole blood sample, by two consecutive centrifugations (3,000 g × 5 min, each) followed by plasma removal. Step 2: cell washings (two times) with phosphate buffer (0.5 mL) followed by centrifugation (3,000 g × 5 min, each) and elimination of supernatant. Step 3: lysis of erythrocytes by twice mixing cells with distilled water (1 mL), followed by centrifugation (15,000 g × 15 min), in order to eliminate the aqueous layers and to obtain the erythrocyte membrane pellet. Step 4: the pellet was added by a mix of 2:1 chloroform:methanol (2 mL) and partitioned with pure water (1 mL) according to the Folch's procedure (39), for the lipid extraction. The organic layer was separated and evaporated under vacuum to dryness. Step 5: thin layer chromatography (TLC) using chloroform/methanol/water 65:25:4 to determine the purity of the phospholipid fraction (40). Step 6: the phospholipid extract was transesterified at room temperature for 10 min with 0.5 M KOH/MeOH to obtain the fatty acid methyl esters (FAMEs), derived from the fatty acid residues present in membrane glycerophospholipids. This chemical procedure avoids oxidative and degradation reactions and gives the fatty acid composition, which represents that of the membrane, as ascertained using appropriate MUFA and PUFA internal standards. FAMEs were extracted by partition between n-hexane (2 mL) and water (0.5 mL), followed by evaporation of the organic phase under vacuum to dryness. Steps 7 and 8: analysis using gas chromatography (GC) as described below.

SCHEME 1

Scheme 1. The eight steps of the fatty acid-based erythrocyte membrane lipidomic analysis. Details in the Experimental section.

GC Analysis of FAME

First, GC analysis of the commercially available reference standard materials for the 10 fatty acids of the cluster was performed as described in Supplementary Information (Step 7 of Figure 1). Calibration curves were obtained for the quantitative analysis of each peak of the chromatogram and are shown in Supplementary Figure 1. Step 8: the FAME mixture obtained from the erythrocyte membrane pellet was dissolved in 20 μL of n-hexane and 1 μL was directly injected to the Agilent 7890B GC (Agilent, Milan) system equipped with a flame ionization detector and a (50%-cyanopropyl)-methylpolysiloxane capillary column (60 m, 0.25 mm i.d., 0.25 μm film thickness) (DB-23, Agilent; Milan). The initial oven temperature was 165°C, held for 3 min, followed by an increase of 1°C/min up to 195°C, held for 40 min, followed by a second increase of 10°C/min up to 240°C, held for 10 min. The carrier gas was hydrogen, held at a constant pressure of 16.482 psi. FAMEs were identified by comparison with the retention times of standard references (see Supplementary Figure 2 for a representative chromatogram) in agreement with previously reported procedures, satisfactorily separating all the 10 fatty acids, without superimposition of other peaks (15, 16, 26, 29). The group of 10 fatty acids corresponds to chromatographic peak areas >97%, and consisted of: two SFAs (palmitic and stearic acids); three MUFAs (palmitoleic, oleic, and cis-vaccenic acids); three PUFA omega-6 (linoleic, dihomogamma linolenic, arachidonic acids); and two PUFA omega-3 (eicosapentaenoic and docosahexaenoic acids). Using the calibration curves, the quantitative values for each peak were obtained as μg/mL, allowing to calculate also the total fatty acid contents (total SFA, total MUFA, total PUFA), ratios between the families (SFA/MUFA, omega-6/omega-3), and three indexes [unsaturation index (UI), peroxidation index (PI), PUFA balance]. Table 1 reports the interval values and the median of each value taking into account that each fatty acid is expressed as % of its quantity (% quant. rel.) respect to the sum of the 10 fatty acid quantities taken as 100%. Table 2 reports such data divided according to the sex of the dogs as median value ± standard deviation (sd).

TABLE 1

Table 1. Interval values of the 10 fatty acids cohort expressed as relative percentages of the μg/mL quantities detected by gas chromatographic analyses^a after isolation and workup of erythrocyte membrane glycerophospholipids of healthy dogs (n = 68).

TABLE 2

Table 2. Cohort of 10 fatty acids expressed as percentages of the found μg/mL quantities, detected by gas chromatographic analyses of the fatty acid methyl esters (FAME)^a after isolation and workup of erythrocyte membrane glycerophospholipids, of male dogs (M = 30) and female dogs (F = 38).

Statistical Methods

Statistical analysis was performed using GraphPad Prism 6.01 software (GraphPad Software, Inc., San Diego, CA). All data were evaluated using a standard descriptive statistic and reported as mean ± sd or as median and range (minimum-maximum), based on their distribution. Normality was checked graphically or using the D'Agostino Pearson test. A comparison between male and female subgroups were done using the unpaired t-test or the Mann-Whitney test. According to the distribution of data, Pearson or Spearman correlations were used to evaluate the relationship between fatty acids percentages or fatty acids indexes with bodyweight and age. The threshold for the statistical significance (p-value) was set up at 0.05. The distribution graphs were produced using Past 3.14 software (free download; Øyvind Hammer, Natural History Museum, University of Oslo).

Results

The Erythrocyte Membrane Lipidome in Healthy Dogs

Blood samples were collected from 68 clinically healthy dogs, including 30 males (6 neutered) and 38 females (12 sterilized), weighting from 2.6 to 43 kg, aged from 2 to 156 months (median 41). Supplementary Table 2 shows age, sex, bodyweight, and breed of the cohort. The work up of the erythrocyte membrane glycerophospholipids to obtain the corresponding FAMEs followed the procedure described in Scheme 1, and each step is detailed in the Experimental section. The FAMEs were separated, recognized, and quantified by GC following known procedures (15, 16, 26, 29) (see Materials and Methods, and Supplementary Figure 2 for a representative GC chromatogram). The cluster of 10 fatty acids contains the following molecules (see Figure 1 for their structures): palmitic (C16:0) and stearic (C18:0) acids as SFA; palmitoleic (C16:1), oleic (9c,C18:1), and vaccenic (11c,C18:1) acids as MUFA; linoleic (LA, C18:2), dihomogammalinolenic (DGLA;C20:3), and arachidonic (AA, C20:4) acids as PUFA omega-6; eicosapentaenoic (EPA, C20:5) and docosahexaenoic (DHA, C22:6) acids as PUFA omega-3. Evaluation of the 10 peaks of the GC chromatogram was performed using calibration curves (see Supplementary Figure 2), obtaining μg/mL values and then expressing each fatty acid as relative quantitative percentage (% rel. quant.) of the total quantities of the 10 fatty acids recognized and calibrated in the GC analysis. The intervals and the median values obtained for our cohort of healthy dogs are reported in Table 1, together with lipid indexes, derived from the 10 fatty acid values, as follows: omega-6/omega-3 ratio, PUFA balance (omega-3/omega-3 + omega-6), SFA/MUFA ratio, UI (UI = MUFA tot × 1 + C18:2 × 2 + C20:3 × 3 + C20:4 × 4 + C20:5 × 5 + C22:6 × 6), and PI (PI = MUFA tot × 0.025 + C18:2 × 1 + C20:3 × 2 + C20:4 × 4 + C20:5 × 6 + C22:6 × 8) (41). A few fatty acids of dog erythrocyte membranes were previously reported (35–37), and their values are listed in Supplementary Table 3.

In Figure 2, the distribution graphics for each fatty acid (% rel. quant.) are shown together with their families (SFA, MUFA, and PUFA omega-6 and omega-3), whereas Figure 3 shows the distribution graphics of the total PUFA and lipid indexes obtained from the 10 fatty acid values. Figures 2, 3 also show minimum and maximum values, as well as the median value of each measured parameter, as listed in Table 1.

FIGURE 2

Figure 2. Distribution of the values in the population of Healthy Dogs (n = 68) for the individual fatty acids and the corresponding families using the data in Tables 1, 2 with 95% confidence interval. Each member of the fatty acid family is given in a row, the last column being the sum of the corresponding fatty acid family. Black: the minimum and the maximum value obtained. Green: the median value.

FIGURE 3

Figure 3. Distribution of the lipid indexes of the population of Healthy Dogs (n = 68) with 95% confidence interval as reported in Tables 1, 2. 1st row: total PUFA, omega-6/omega-3 and SFA/PUFA ratios; 2nd row: SFA/MUFA ratio, unsaturation and peroxidation indexes. Black: the minimum and the maximum value obtained. Green: the median value.

The Correlations of Fatty Acids With Dog Characteristics

We used the above described interval values to evaluate correlations between FAME types, families, and lipid indexes with the dog characteristics. In Supplementary Information, the correlation graphs are shown (Supplementary Figures 3–6). For the age the only correlation was found with eicosapentaenoic acid (EPA, omega-3; p < 0.001, r = 0.396). With bodyweight, several fatty acid parameters correlated significantly, as follows: positive correlation with palmitic acid (C16:0; p = 0.001, r = 0.385) and the total amount of SFA (p = 0.001, r = 0.402). The MUFA palmitoleic acid (C16:1) significantly correlated in a positive manner to bodyweight (p = 0.007; r = 0.326). The bodyweight negatively correlated with arachidonic acid (p = 0.036, r = 0.257), the total content of the omega-6 PUFA (p = 0.003, r = 0.361), and the total amount of PUFA (p = 0.004, r = 0.347). Moreover, both UI and PI showed negative correlation with bodyweight (p =0.031, r = 0.301 and p = 0.037, r = 0.256, respectively). A statistical analysis based on the different breeds of the healthy dogs was not possible because of the limited number of samples for each breed and in addition 24 out of the total 68 (35%) of the healthy dogs recruited for this study were mixed-breed.

Subsequently, two groups of male and female healthy dogs were examined to envisage differences of the membrane RBC fatty acid composition between the two genders. Because of the small number of sterilized females (n = 12) and neutered males (n = 6), we did not consider this additional characteristic to further discriminate groups. The values for the individual FAME and the lipid indexes for male and female healthy dogs are reported as mean ± sd values in Table 2, expressed as quantitative relative percentages (% rel. quant.) of each value over the total 10 fatty acids values, as explained above. Significant differences were found for the lower levels of stearic and linoleic acids (p = 0.0126 and p = 0.0266, respectively) and higher levels of palmitoleic acid (p = 0.0344) in female than in male dogs. The corresponding graphical representations are shown Supplementary Figures 7, 8.

Discussion

For the first time, the erythrocyte membrane lipidome characterization of the fatty acid content of glycerophospholipids is provided in a cohort of clinically healthy dogs (n = 68). The rationale for the choice of 10 fatty acids takes into account the most representative components in terms of structural and functional roles played by such hydrophobic moieties in membrane phospholipids, as well as the nutritionally important PUFAs (see Figure 1). We arbitrarily chose this cohort for the application to veterinary analyses, taking into consideration our and others' work done in humans with physiological and pathological conditions (15–18, 25–28). It is worth underlining that studies reported in literature on fatty acids of erythrocyte membranes in both humans and animals do not report the same cohort of fatty acids. For humans, the meta-analysis reported by Hodson et al. (28) highlighted 13 SFA, MUFA, and PUFA moieties in erythrocyte membrane phospholipids obtained from nine studies of 321 men and 70 women. In another study on infants, only the total SFA and MUFA values are given without specifying each fatty acid type, whereas for the omega-6 pathway together with the C18 and C20 PUFA long-chain fatty acids such as 22:4 and 22:5 are evaluated, which are not commonly reported in other studies (21). With the emerging field of lipidomics, including our own research, new fatty acid structures are evidenced, such as for example the presence of the n-10 MUFA family, with 6cis-C16:1 (sapienic acid) as positional isomer of palmitoleic acid, and the corresponding geometrical trans isomer (6trans-C16:1). These new markers indicate the metabolic partition of delta-6 desaturase enzymes between palmitic and linoleic acids and the free radical-based transformations (11), respectively, but it is not possible to define yet their participation and effects for the membrane composition. In model animals as well as in humans, increased or decreased fatty acid levels were found in several diseases such as for example, for kidney (42) or liver diseases (43) being such alterations correlated with other biochemical markers. Newly identified fatty acid markers need further research to establish their normal and altered levels. Considering such complexity, it was evaluated of strategic importance to concentrate on a specific number of fatty acids, that is, the 10 fatty acid cohort, and on a specific cell compartment, that is, the erythrocyte membrane, in order to get a simplified, yet representative and significative, approach for the examination of known metabolic and nutritional influences on the biophysically regulated homeostatic control of membrane formation and remodeling. The cohort includes fatty acids that cannot be missed in biological membranes, that are characteristic for each tissue type (8, 14, 24, 44) and their interval values (minimum and maximum values) in the specific membrane type, such as in erythrocytes, derive from population studies, as previously mentioned.

It is worth mentioning that, when fatty acid analysis is performed, it is very important to ascertain that the subjects did not follow any supplementation during the last 4 months. We asked this information to the owners during the recruitment, and no specific supplementation or regular consumption of specifically enriched food emerged. We remark that in canine diets, it is becoming more and more popular the use of omega-3 supplements or commercial “omega-3 enriched” diets, that lead to the uptake of this EFA family. It is known that such supplementations can change the results of lipid analyses (45).

The analytical methodology is also important to ensure the efficient separation of the 10 fatty acid cohort especially of unsaturated fatty acids with the same fatty acid chain length but different geometry and position of the double bonds. In this respect, the library of geometrical and positional isomers in our hands plays an important role to check peak superimposition (11, 15, 16, 29). The overall fatty acid distribution reported in Table 1 and Figures 2, 3 highlights that in the dog erythrocyte membrane lipidome the omega-6 PUFA content is prevalent followed by SFA and MUFA. It is worth underlining that omega-3 fatty acids are present in minimal concentrations, marking the difference with human lipidome (11, 18, 24, 30–32). Moreover, in healthy dogs the SFAs, stearic and palmitic acids, and the omega-6 PUFAs, arachidonic, and linoleic acids, sum up to nearly the 94% of the total membrane RBC fatty acids, as already observed in other mammalian species (24). In the distribution graphs of Figure 2, large interval ranges for the palmitic acid percentage (8.2–25.8) and arachidonic acid (17.5–43.7) can be appreciated. It is also interesting to note that in 80% of the dogs' narrower ranges are described (palmitic acid: 8.5–20%, arachidonic acid: 30–40%). On the other hand, the distribution of linoleic acid values follows a gaussian-like behavior, with a probability distribution symmetric about the mean, showing that data near the mean are more frequent in occurrence than data far from the mean. The different distributions of SFA, C18, and C20 PUFA omega-6 indicate that it is necessary to gather large sets of these data for mathematical treatment to model their roles in the membrane compositions (46). It is also worth pointing out again that the omega-6 are EFA, not prepared directly by biosynthesis, and also dogs need the intake of linoleic acid as precursor of other biosynthetically prepared omega-6 fatty acids. The importance of ascertaining the omega-6 dietary supply to animals represents an historical highlight in veterinary medicine, having been noted in canine dermatological problems by Burr and Burr (47). Therefore, it is of extreme importance to follow up the intakes of linoleic acid and its metabolic transformations to DGLA and arachidonic acid, together with their incorporation at the level of cell membranes, since they regulate fundamental structural and functional properties for health in the whole. This information can be integrated in a panel with other biochemical, biomolecular, clinical parameters for a thorough evaluation of health in dogs.

The strength of our cohort of fatty acids is to provide a valuable set of metabolic information to examine: (1) data on the PUFA pathways that start from precursors that must be necessarily taken from the diet thus reporting on the omega-6 and omega-3 fatty acid balance, that in dogs such balance is shifted toward the omega-6 family; and (2) data on the SFA-MUFA pathway that starts from the de-novo biosynthesis of the first fatty acid in the body, that is, the SFA palmitic acid, thus including the effects of insulin response and liver conditions (1, 24). From these crucial data sets, the influence of enzymatic activities, hormones, and nutrients can be studied, and it is possible to use the specific SFA-MUFA-PUFA cohort to examine metabolic profiles, decide nutritional-based strategies, and follow up on the changes during physiological and pathological conditions, including oxidative stress (1–3, 9, 11). The limitations are due to the small size of our data, compared with the large data set provided by other powerful lipidomic techniques, although the balance cost-effectiveness are important parameters to evaluate in wide applications to population. Another limitation of our study is in the non-homogeneous distribution of male and female dogs and in the diversity of breeds, but the prosecution of data collection will certainly compensate such weaknesses.

The omega-6 EFA linoleic acid and its transformation to DGLA and ARA, which is predominant in dogs, affect not only from the point of view of the membrane structural properties but also affects the membrane response to stimuli. In fact, upon different kinds of stimuli, phospholipase A₂ is activated and liberates the fatty acid moieties from the cell membrane phospholipids (48); therefore, the membrane fatty acid composition can be considered a precious information of the pro- and anti- inflammatory predisposition of the organism (49, 50). Indeed, liberation of SFA, MUFA, and PUFA fatty acids generates potent mediators, such as prostaglandins, leukotrienes, endocannabinoids, and many others that influence dog metabolism in health and diseases (51–53). There is a growing attention to the lipidomic analysis in animal studies, as reported recently for variations of phospholipids of whole blood and plasma in canine inflammatory bowel disease (53). In this context we believe that the present work on erythrocyte membrane fatty acid cohort in healthy animals could represent the benchmark for assessment of healthy dog molecular profile.

A focus of our approach in dogs was to observe the omega-3 value and its distribution in healthy dogs, since omega-3 are EFAs, and EPA and DHA levels derive from the dietary intake of the precursor (alpha-linolenic acid). It is well-known that this pathway in animals is less prevalent than omega-6 (52). Our data confirmed low levels and a narrow range (0.5–4) (Table 1) together with a significant age-related increase of EPA (Supplementary Figure 3). It is worth noting that in humans EPA and DHA levels are strongly influenced by dietary intakes, but some studies have shown that age remains a significant determinant of EPA and DHA status regardless of dietary intake (54). Omega-3 are known to have specific anti-inflammatory activities and, being the omega-6/omega-3 ratio naturally high in dogs (12-7-83.7), it is evident that their increase in the cellular environment, influencing this ratio, can exert significant metabolic effects. In this context, the establishment of omega-3 levels under various dog health conditions can be interesting to precisely individuate levels of omega-3 to introduce in the regular diets of animals.

Considering the correlations found in our work, it is worth underlining that de novo lipogenesis with formation of palmitic acid and palmitoleic acid are known biomarkers of bodyweight increase in animals and humans (15, 55), and their increase over the normal levels are indicators of metabolic derangements in humans at older ages (56). We found positive correlation between bodyweight and the values of palmitoleic, stearic acids, and total SFA of our healthy dog cohort. At the same time, the negative correlation between bodyweight with all unsaturated indicators (arachidonic acid, PUFA omega-6, UI, PI) indicates a shift of membrane components from unsaturated to saturated moieties, that will be worth of further deepening, especially in the dog obesity conditions. An increase of SFA and MUFA components has been correlated with suppression of inflammatory response in animal models (55).

Correlations of the fatty acid data with different breeds need more data, whereas interesting observations could be done related to dog sex. The lower levels of stearic acid and the higher levels of palmitoleic acid detected in male dogs indicate that SFA to MUFA pathway functions more than in females. Palmitoleic acid is a lipokine that travels to muscles and liver, where it improves cell sensitivity to insulin, blocks fat accumulation in the liver (57), and seems to be involved in regulating muscle mass (57, 58). Although preliminarily, this different distribution of FA indirectly indicates differences in metabolism of dogs according to the sex, to be deepened also considering larger numbers of sterilized females and neutered males.

Conclusion

In this paper, we provided the first panel of erythrocyte membrane fatty acids in healthy dogs choosing a cohort representative for the main structural and functional roles of these hydrophobic molecules in the cell membrane compartment. The panel of 10 fatty acids represents a set of information, chosen for their biochemical, nutritional, and metabolic significance, and their reference ranges for healthy dogs can represent a benchmark to evaluate pathological conditions. It is worth to remark that the effects of fatty acids in cellular metabolism departing from the membrane release and remodeling processes are more and more studied also for immune functions (59); therefore, the knowledge of the membrane asset can be a useful indication for the examination of animal response at molecular level. Indeed, the protocol herein described, including calibration and quantitation procedures, can be applied to large animal populations, with molecular information that are ready to be combined with other medical data sets for mining techniques, which are nowadays under study for the exploitation of disease risks and prevention plans (60).

Clinically speaking, since imbalances in lipid metabolism contributes to diverse animal phenotypes and disease states, ranging from inflammation and cancer to metabolic diseases, the study of membrane lipids in normal and pathological states can successfully contribute to monitor health conditions as well as the result of treatments. Membrane lipidomic research in veterinary will provide, as done for humans, deep insights in mechanisms of several diseases, helping to identify molecular targets, patient selection for better treatment response, prediction and monitoring of treatment efficacy, and response to therapeutic measures. In the light of this, membrane lipidomics may become part of the arsenal adopted for clinical diagnostics in veterinary medicine, useful for the evaluation of a disease onset and its progression, representing a potential avenue to individualized treatment and monitoring in different pathological states.

Data Availability Statement

All datasets generated for this study are included in the article/Supplementary Material.

Ethics Statement

This animal study was reviewed and approved by Health Ministry and the Committee on Animal Research and Ethics of the Universities of Chieti-Pescara, Teramo and Experimental Zooprophylactic Institute of AeM (CEISA), Protocol UNICHD12 n. 1168. Written informed consent was obtained from the owners for the participation of their animals in this study.

Author Contributions

CF, CC, and AB: conceptualization. CF, CC, AB, and PC: methodology and supervision. PC, AG, MD, BB, and FD: recruitment and clinical data. AS and PP: lipidomic analyses. CF, PC, AB, CC, and AG: data curation. CF: writing—original draft preparation. CF, CC, AB, PC, AG, and PP: writing—review and editing. CF, AB, and AG: funding acquisition. All authors contributed to the article and approved the submitted and revised versions.

Funding

PP was supported in the frame of the Rep-Eat doctorate programme—funding: European Union's Horizon 2020 under the Marie Skłodowska-Curie Grant Agreement No. 713714. Co-funding institution: University of Teramo, Abruzzo Region. Partner institution: ISOF, CNR.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors acknowledge the Fondazione Di Bella for covering some laboratory expenses.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2020.00502/full#supplementary-material

References

1. Han W. Lipidomics for studying metabolism. Nat Rev Endocrinol. (2016) 12:668–79. doi: 10.1038/nrendo.2016.98

CrossRef Full Text | Google Scholar

2. Ni Z, Goracci L, Cruciani G, Fedorova M. Computational solutions in redox lipidomics – current strategies and future perspectives. Free Radic Biol Med. (2019) 144:110–23. doi: 10.1016/j.freeradbiomed.2019.04.027

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Zhao Y-T, Cheng X-I, Lin E-C. Lipidomics applications for discovering biomarkers of diseases in clinical chemistry. Int Rev Cell Mol Biol. (2014) 313:1–26. doi: 10.1016/B978-0-12-800177-6.00001-3

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Quehenberger O, Armando AM, Brown AH, Milne SB, Myers DS, Merrill AH, et al. Lipidomics reveals a remarkable diversity of lipids in human plasma. J Lipid Res. (2010) 51:3299–305. doi: 10.1194/jlr.M009449

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Wong MWK, Braidy N, Pickford R, Vafaee F, Crawford J, Muenchhoff J, et al. Plasma lipidome variation during the second half of the human lifespan is associated with age and sex but minimally with BMI. PLoS ONE. (2019) 14:e0214141. doi: 10.1371/journal.pone.0214141

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Jain M, Ngoy S, Sheth SA, Swanson RA, Rhee EP, Liao R, et al. A systematic survey of lipids across mouse tissues. Am J Physiol Endocrinol Metab. (2014) 306:E854–68. doi: 10.1152/ajpendo.00371.2013

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Klose C, Surma MA, Gerl MJ, Meyenhofer F, Shevchenko A, Simons K. Flexibility of a eukaryotic lipidome – insights from yeast lipidomics. PLoS ONE. (2012) 7:e35063. doi: 10.1371/journal.pone.0035063

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Pradas I, Huynh K, Cabré R, Ayala V, Meikle PJ, Jové M, et al. Lipidomics reveals a tissue-specific fingerprint. Front Physiol. (2018) 9:1165. doi: 10.3389/fphys.2018.01165

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Zhao YY, Cheng XL, Lin RC, Wei F. Lipidomics applications for disease biomarker discovery in mammal models. Biomark Med. (2015) 9:153–68. doi: 10.2217/bmm.14.81

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Casares D, Escribá PV, Rosselló CA. Membrane lipid composition: effect on membrane and organelle structure, function and compartmentalization and therapeutic avenues. Int J Mol Sci. (2019) 20:2167. doi: 10.3390/ijms20092167

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Ferreri C, Masi A, Sansone A, Giacometti G, Larocca AV, Menounou G, et al. Fatty acids in membranes as homeostatic, metabolic and nutritional biomarkers: recent advancements in analytics and diagnostics. Diagn Basel Switz. (2016) 7:1. doi: 10.3390/diagnostics7010001

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Zhu D, Bungart BL, Yang X, Zhumadilov Z, Lee JC-M, Askarova S. Role of membrane biophysics in Alzheimer's-related cell pathways. Front Neurosci. (2015) 9:186. doi: 10.3389/fnins.2015.00186

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Diaz M, Retamal MA. Modulation of ion channels and ionic pumps by fatty acids: implications in physiology and pathology. Front Physiol. (2019) 9:1625. doi: 10.3389/fphys.2018.01625

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Ferreri C, Chatgilialoglu C. Role of fatty acid-based functional lipidomics in the development of molecular diagnostic tools. Expert Rev Mol Diagn. (2012) 12:767–80. doi: 10.1586/erm.12.73

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Sansone A, Tolika E, Louka M, Sunda V, Deplano S, Melchiorre M, et al. Hexadecenoic fatty acid isomers in human blood lipids and their relevance for the interpretation of lipidomic profiles. PLoS ONE. (2016) 11:e0152378. doi: 10.1371/journal.pone.0152378

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Amézaga J, Arranz S, Urruticoechea A, Ugartemendia G, Larraioz A, Louka M, et al. Altered red blood cell membrane fatty acid profile in cancer patients. Nutrients. (2018) 10:1853. doi: 10.3390/nu10121853

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Gai Z, Wang T, Visentin M, Kullak-Ublick GA, Fu X, Wang Z. Lipid accumulation and chronic kidney disease. Nutrients. (2019) 11:722. doi: 10.3390/nu11040722

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Harris WS, Pottala JV, Varvel SA, Borowski JJ, Ward JN, McConnell JP. Erythrocyte omega-3 fatty acids increase and linoleic acid decreases with age: observations from 160,000 patients. Prostaglandins Leukot Essent Fatty Acids. (2013) 88:257–63. doi: 10.1016/j.plefa.2012.12.004

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Furtado JD, Beqari J, Campos H. Comparison of the utility of total plasma fatty acids versus those in cholesteryl ester, phospholipid, and triglyceride as biomarkers of fatty acid intake. Nutrients. (2019) 11:2081. doi: 10.3390/nu11092081

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Sender R, Fuchs S, Milo R. Revised estimates for the number of human and bacteria cells in the body. PLoS Biol. (2016) 14:e1002533. doi: 10.1371/journal.pbio.1002533

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Makrides M, Neumann MA, Byard RW, Simmer K, Gibson RA. Fatty acid composition of brain, retina, and erythrocytes in breast- and formula-fed infants. Am J Clin Nutr. (1994) 60:189–94. doi: 10.1093/ajcn/60.2.189

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Baur LA, O'Connor J, Pan DA, Wu BJ, O'Connor MJ, Storlien LH. Relationships between the fatty acid composition of muscle and erythrocyte membrane phospholipid in young children and the effect of type of infant feeding. Lipids. (2000) 35:77–82. doi: 10.1007/s11745-000-0497-y

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Acar N, Berdeaux O, Grégoire S, Cabaret S, Martine L, Gain P, et al. Lipid composition of the human eye: are red blood cells a good mirror of retinal and optic nerve fatty acids? PLoS ONE. (2012) 7:e35102. doi: 10.1371/journal.pone.0035102

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Abbott SK, Else PL, Atkins TA, Hulbert AJ. Fatty acid composition of membrane bilayers: importance of diet polyunsaturated fat balance. Biochim Biophys Acta. (2012) 1818:1309–17. doi: 10.1016/j.bbamem.2012.01.011

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Heude B, Ducimetière P, Berr C, EVA Study. Cognitive decline and fatty acid composition of erythrocyte membranes–the EVA study. Am J Clin Nutr. (2003) 77:803–8. doi: 10.1093/ajcn/77.4.803

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Giacometti G, Ferreri C, Sansone A, Chatgilialoglu C, Marzetti C, Spyratou E, et al. High predictive values of RBC membrane-based diagnostics by biophotonics in an integrated approach for autism spectrum disorders. Sci Rep. (2017) 7:9854. doi: 10.1038/s41598-017-10361-7

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Pallot C, Mazzocco J, Meillon C, Semama DS, Chantegret C, Ternoy N, et al. Alteration of erythrocyte membrane polyunsaturated fatty acids in preterm newborns with retinopathy of prematurity. Sci Rep. (2019) 9:7930. doi: 10.1038/s41598-019-44476-w

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Hodson L, Skeaff CM, Fielding BA. Fatty acid composition of adipose tissue and blood in humans and its use as a biomarker of dietary intake. Prog Lipid Res. (2008) 47:348–80. doi: 10.1016/j.plipres.2008.03.003

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Sansone A, Melchiorre M, Chatgilialoglu C, Ferreri C. Hexadecenoic fatty acid isomers: a chemical biology approach for human plasma biomarker development. Chem Res Toxicol. (2013) 26:1703–9. doi: 10.1021/tx400287u

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Sieber-Ruckstuhl NS, Burla B, Spoerel S, Schmid F, Venzin C, Cazenave-Gassiot A, et al. Changes in the canine plasma lipidome after short- and long-term excess glucocorticoid exposure. Sci Rep. (2019) 9:6015. doi: 10.1038/s41598-019-42190-1

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Hall JA, Jewell DE. Feeding healthy beagles medium-chain triglycerides, fish oil, and carnitine offsets age-related changes in serum fatty acids and carnitine metabolites. PLoS ONE. (2012) 7:e49510. doi: 10.1371/journal.pone.0049510

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Boretti FS, Burla B, Deuel J, Gao L, Wenk MR, Liesegang A, et al. Serum lipidome analysis of healthy beagle dogs receiving different diets. Metabolomics. (2020) 16:1. doi: 10.1007/s11306-019-1621-3

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Snigdha S, Astarita G, Piomelli D, Cotman CW. Effects of diet and behavioral enrichment on free fatty acids in the aged canine brain. Neuroscience. (2012) 202:326–33. doi: 10.1016/j.neuroscience.2011.12.002

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Risso A, Pellegrino FJ, Relling AE, Corrada Y. Effect of long-term fish oil supplementation on semen quality and serum testosterone concentrations in male dogs. Int J Fertil Steril. (2016) 10:223–31. doi: 10.22074/ijfs.2016.4913

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Mehta JR, Braund KG, Hegreberg GA, Thukral V. Lipid fluidity and composition of the erythrocyte membrane from healthy dogs and Labrador retrievers with hereditary muscular dystrophy. Neurochem Res. (1991) 16:129–35. doi: 10.1007/BF00965700

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Fuhrmann H, Zimmermann A, Gück T, Oechtering G. Erythrocyte and plasma fatty acid patterns in dogs with atopic dermatitis and healthy dogs in the same household. Can J Vet Res. (2006) 70:191–6.

PubMed Abstract | Google Scholar

37. Stoeckel K, Nielsen LH, Fuhrmann H, Bachmann L. Fatty acid patterns of dog erythrocyte membranes after feeding of a fish-oil based DHA-rich supplement with a base diet low in n-3 fatty acids versus a diet containing added n-3 fatty acids. Acta Vet Scand. (2011) 53:57. doi: 10.1186/1751-0147-53-57

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Siobhan S, Mostyn A, Rutland CS. Fatty acids in veterinary medicine and research. In: Catala A, editor. Fatty Acids Techopen. London: IntechOpen (2017). p. 177–97. doi: 10.5772/intechopen.68440

CrossRef Full Text | Google Scholar

39. Folch J, Lees M, Sloane Stanley GH. A simple method for the isolation and purification of total lipides from animal tissues. J Biol Chem. (1957) 226:497–509.

PubMed Abstract | Google Scholar

40. Fuchs B, Süss R, Teuber K, Eibisch M, Schiller J. Lipid analysis by thin-layer chromatography–a review of the current state. J Chromatogr A. (2011) 1218:2754–74. doi: 10.1016/j.chroma.2010.11.066

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Hulbert AJ, Pamplona R, Buffenstein R, Buttemer WA. Life and death: metabolic rate, membrane composition, and life span of animals. Physiol Rev. (2007) 87:1175–213. doi: 10.1152/physrev.00047.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Zhao Y-Y, Cheng X-L, Wei F, Bai X, Tan X-J, Lin R-C, et al. Intrarenal metabolomic investigation of chronic kidney disease and its TGF-β1 mechanism in induced-adenine rats using UPLC Q-TOF/HSMS/MS^E. J Proteome Res. (2012) 12:692–703 doi: 10.1021/pr3007792

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Svegliati-Baroni G, Pierantonelli I, Torquato P, Marinelli R, Ferreri C, Chatgilialoglu C, et al. Lipidomic biomarkers and mechanisms of lipotoxicity in non-alcoholic fatty liver disease. Free Radic Biol Med. (2019) 144:293–309. doi: 10.1016/j.freeradbiomed.2019.05.029

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Harayama T, Riezman H. Understanding the diversity of membrane lipid composition published correction appears. Nat Rev Mol Cell Biol. (2018) 19:281–96. doi: 10.1038/nrm.2017.138

CrossRef Full Text | Google Scholar

45. Walker CG, West AL, Browning LM, Madden J, Gambell JM, Jebb SA, et al. The pattern of fatty acids displaced by EPA and DHA following 12 months supplementation varies between blood cell and plasma fractions. Nutrients. (2015) 7:6281–93. doi: 10.3390/nu7085285

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Mar JC. The rise of the distributions: why non-normality is important for understanding the transcriptome and beyond. Biophys Rev. (2019) 11:89–94. doi: 10.1007/s12551-018-0494-4

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Burr GO, Burr MM. Nutrition classics from the journal of biological chemistry 82:345-67, 1929. A new deficiency disease produced by the rigid exclusion of fat from the diet. Nutr Rev. (1973) 31:248–9.

PubMed Abstract | Google Scholar

48. Murakami M, Taketomi Y, Miki Y, Sato H, Hirabayashi T, Yamamoto K. Recent progress in phospholipase A2 research: from cells to animals to humans. Prog Lipid Res. (2011) 50:152–92. doi: 10.1016/j.plipres.2010.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Adler N, Schoeniger A, Fuhrmann H. Polyunsaturated fatty acids influence inflammatory markers in a cellular model for canine osteoarthritis. J Anim Physiol Anim Nutr. (2018) 102:e623–32. doi: 10.1111/jpn.12804

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Loef M, Schoones JW, Kloppenburg M, Ioan-Facsinay A. Fatty acids and osteoarthritis: different types, different effects. Joint Bone Spine. (2019) 86:451–8. doi: 10.1016/j.jbspin.2018.07.005

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Bauer JJE. Essential fatty acid metabolism in dogs and cats. Rev Bras Zootec. (2008) 37:20–27. doi: 10.1590/S1516-35982008001300004

CrossRef Full Text | Google Scholar

52. Abramo F, Campora L, Albanese F, della Valle MF, Cristino L, Petrosino S, et al. Increased levels of palmitoylethanolamide and other bioactive lipid mediators and enhanced local mast cell proliferation in canine atopic dermatitis. BMC Vet Res. (2014) 10:21. doi: 10.1186/1746-6148-10-21

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Kalenyak K, Heilmann RM, van de Lest CHA, Brouwers JF, Burgener IA. Comparison of the systemic phospholipid profile in dogs diagnosed with idiopathic inflammatory bowel disease or food-responsive diarrhea before and after treatment. PLoS ONE. (2019) 14:e0215435. doi: 10.1371/journal.pone.0215435

CrossRef Full Text | Google Scholar

54. Johnson ML, Lalia AZ, Dasari S, Pallauf M, Fitch M, Hellerstein MK, et al. Eicosapentaenoic acid but not docosahexaenoic acid restores skeletal muscle mitochondrial oxidative capacity in old mice. Aging Cell. (2015) 14:734–43. doi: 10.1111/acel.12352

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Guo X, Li H, Xu H, Halim V, Zhang W, Wang H, et al. Palmitoleate induces hepatic steatosis but suppresses liver inflammatory response in mice. PLoS ONE. (2012) 7:e39286. doi: 10.1371/journal.pone.0039286

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Zong G, Ye X, Sun L, Li H, Yu Z, Hu FB, et al. Associations of erythrocyte palmitoleic acid with adipokines, inflammatory markers, and the metabolic syndrome in middle-aged and older Chinese. Am J Clin Nutr. (2012) 96:970–6. doi: 10.3945/ajcn.112.040204

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Cao H, Gerhold K, Mayers JR, Wiest MM, Watkins SM, Hotamisligil GS. Identification of a lipokine, a lipid hormone linking adipose tissue to systemic metabolism. Cell. (2008) 134:933–44. doi: 10.1016/j.cell.2008.07.048

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Lipina C, Hundal HS. Lipid modulation of skeletal muscle mass and function. J Cachexia Sarcopenia Muscle. (2017) 8:190–201. doi: 10.1002/jcsm.12144

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Alarcón PA, Manosalva C, Quiroga J, Belmar I, Álvarez K, Díaz G, Taubert A, et al. Oleic and linoleic acids induce the release of neutrophil extracellular traps via pannexin 1-dependent ATP release and P2X1 receptor activation. Front Vet Sci. (2020) 7:260. doi: 10.3389/fvets.2020.00260

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Frigolet ME, Gutiérrez-Aguilar R. The role of the novel lipokine palmitoleic acid in health and disease. Adv Nutr Bethesda Md. (2017) 8:173–81S. doi: 10.3945/an.115.011130

PubMed Abstract | CrossRef Full Text