The immune system undergoes marked changes during aging characterized by a state of chronic, low-grade inflammation termed ‘inflammaging’. We explore this phenomenon in domestic dogs, which are the most morphologically and physiologically diverse group of mammals, with the widest range in body sizes for a single species. Additionally, smaller dogs tend to live significantly longer than larger dogs across all breeds. Body size is intricately linked to mass-specific metabolism and aging rates, which suggests that dogs are exemplary for studies in inflammaging. Dermal fibroblast cells play an important role in skin inflammation, making them a good model for inflammatory patterns across dog breed, body sizes and ages. Here, we examined aerobic and glycolytic cellular metabolism, and IL-6 concentrations in primary fibroblast cells isolated from small and large dog breeds, that were either recently born puppies or old dogs after death. We found no differences in cellular metabolism when isolated fibroblasts were treated with lipopolysaccharide (LPS) from Escherichia coli to stimulate an inflammatory phenotype. Unlike responses observed in mice and humans, there was a less drastic amplification of IL-6 concentration after LPS treatment in the geriatric population of dogs compared with recently born dogs. In young dogs, we also found evidence that untreated fibroblasts from large breeds had significantly lower IL-6 concentrations than observed for smaller breeds. This implies that the patterns of inflammaging in dogs may be distinct and different from other mammals commonly studied.

Physiological aging, which is difficult to define biologically, describes the decrease in physiological function and fitness with age (Ricklefs, 2010). It is progressive, endogenously derived, irreversible and expands across multiple processes and systems within an organism (Cohen, 2018). Owing to the accumulation of deleterious traits in its forward progress (Sanz et al., 2006), every aspect of an organism's phenotype is modified during aging (Kirkwood 2005); thus, lifespan is limited by the rate of an organism's aging. Various processes, including accumulation of oxidative damage, increases in inflammation, telomere shortening, mitochondrial dysfunction and misfolded proteins may affect the rate of aging (Cohen, 2018). Age-related physiological changes to the immune system are considered a main culprit of co-morbidities in old age and are often associated with life expectancy and survival across species (Cossarizza et al., 1997; Strasser et al., 2000). A decline in mitochondrial quality has been linked to specific mechanisms of aging including cellular dysfunction and inflammation (Sun et al., 2016). Mitochondria dysfunction with age can lead to activation of NIrp3 inflammasome, a multiprotein complex that leads to the activation of proinflammatory cytokines IL-1β and IL-18, and the induction of chronic reactive oxygen species (ROS) production (Franceschi and Campisi, 2014). Similarly, the immune system undergoes marked changes during aging characterized by a state of chronic, low-grade inflammation: a process called ‘inflammaging’ (Franceschi and Campisi, 2014). Specifically, the aging process activates inflammatory molecules, including cytokines, inevitably resulting in systemic inflammation (Baylis et al., 2013). Pro-inflammatory cytokines such as interleukin (IL)-1β and IL-6 in particular are termed the ‘cytokines of gerontology’ (Franceschi et al., 2000). These inflammatory molecules contribute to inducible, innate immune responses and their release is triggered during immune challenges (Lee, 2006; Zimmerman et al., 2014). For example, IL-6 induces B cell differentiation, activates T cells and thymocytes, and activates acute phase proteins, among other immunological roles (Yamashita et al., 1994). Mechanistically, inflammaging seems to be the consequence of cumulative lifetime exposure to antigenic load (Baylis et al., 2013). It can be caused by accumulation of pro-inflammatory tissue damage, accumulation of pathogens due to an ill-functioning immune system, dysfunctional host cells, to name a few (López-Otín et al., 2013).

Body size is closely related to lifespan and aging rates so that, typically, larger mammals are longer lived than smaller mammals (Promislow, 1993). Because the maintenance of the immune system is a metabolic process, we would expect body mass to have an effect on immunity, particularly in innate responses (Downs et al., 2020; Dingli and Pacheco, 2006). For example, constitutive concentrations of neutrophils, a white blood cell recruited to sites of infection early during an immune response that help eliminate bacterial challenges and the propagation of subsequent immune defense (Kolaczkowska and Kubes, 2013), increases disproportionately with body mass in mammals (Downs et al., 2020). The relationship between a trait and body size is also connected with life history evolution of that species because an individual can only grow large by having an extended maturation period to grow for longer, which is often associated with increased longevity (Harrison, 2017). Much research in ecoimmunology has focused on exploring the prediction that larger animals that tend to be longer-lived would have a higher investment in safety because they live longer (Lee, 2006).

Domestic dogs are a morphologically and physiologically diverse group of mammals. Underlying observed morphological differences there must be cellular-level disparities that could lead to elucidating aging rates and lifespan discrepancies across dog breeds. Smaller dogs tend to live significantly longer than larger dogs across all breeds (Jimenez, 2016). Smaller dogs also have lower cancer risks (Michell, 1999), and demonstrate lower prevalence of age-related diseases such as cataracts (Urfer et al., 2011), implying slower rates of aging. Often, the aging literature refers to the GH–IGF-1 (growth hormone–insulin-like growth factor 1) pathway as the most conserved aging pathway across species (López-Otín et al., 2013; Cohen, 2018). The somatotrophic axis in mammals is composed of growth hormone, produced by the anterior pituitary and IGF-1, which is produced in response to GH in many cell types. This pathway informs cells on the presence of glucose. Insulin is involved in cell growth, and glucose homeostasis (Partridge and Gems, 2002) whereas IGF-1 is thought to influence cell growth, differentiation and survival (Dantzer and Swanson, 2012). Insulin and IGF-1 concentration seem to be negatively correlated with lifespan in several mouse strains selected for disparate lifespans (Fontana et al., 2010; Tian et al., 2017), where mice selected for dwarfism due to lack of GH production and a reduction of IGF-1 in plasma had longer lives than their control counterparts (Dantzer and Swanson, 2012). A single IGF1 haplotype seems to substantially contribute to size variation in dogs (Sutter et al., 2007), and serum IGF-1 concentration in small dogs is reduced relative to that found in large breeds, providing a potential link to their longer lives (Eigenmann et al., 1988; Greer et al., 2007).

Particularly in human studies, aging has been demonstrated to increase levels of IL-6 and tumor necrosis factor alpha (TNF-α) and decrease IGF-1 levels (Xia et al., 2016), a pattern unexplored in dogs. IL-6 concentrations significantly predict mortality and predicted 2-year survival in very old humans, and seem to be positively associated with other age-related diseases such as type 2 diabetes, cardiovascular disease, Parkinson's disease, cancer and arthritis (Wikby et al., 2006; Baylis et al., 2013; Xia et al., 2016). Thus, the increases of IL-6 associated with an aging immune system may contribute to increases in incidences of disease and cancer (Kearns et al., 1999; Wikby et al., 2006).

To what degree donor body mass and/or age influence the inflammatory response in vitro is unknown for any mammal, much less for dogs, but there is no question that inflammatory dysfunction is a hallmark of aging (Franceschi and Campisi, 2014). In aging dogs, there seems to be marked changes in T-cell immunity, including changes in the number, and subsets of T-cells, which may lead to a reduced ability to respond to non-specific mitogens (Day, 2010). Older dogs also seem to have more immunoglobins than puppies (Day, 2010). Older female dogs seem to demonstrate increases in IL-1 activity compared with younger dogs (Day, 2010). There have been studies in dogs demonstrating decreases in mitogen stimulation (Greeley et al., 1996), number of white blood cells, and immature neutrophils (Strasser et al., 1993) with increasing age. However, none of these studies have addressed whether these changes differ between the different breed sizes as they age.

Here, we examine aerobic and glycolytic cellular metabolism, and IL-6 concentration in primary fibroblast cells isolated from small and large breed, recently born and old deceased dogs. Isolated fibroblasts were then treated with lipopolysaccharide (LPS) from Escherichia coli. LPS molecules from the cell wall of Gram-negative bacteria are an immunostimulatory agent that initiates synthesis and release of proinflammatory cytokines, including IL-6 (Agarwal et al., 1995). We predicted that (1) aerobic and anaerobic metabolism would increase with LPS treatment; (2) that small breeds (longer lived) would have a better capacity to deal with LPS treatment, either through smaller increases in aerobic and anaerobic metabolism or no change to spare respiratory capacity, compared with large breed dogs; (3) that recently born dogs would be better suited to deal with an LPS treatment; and (4) IL-6 concentrations would be greater in deceased dogs compared with recently born dogs. Expression of IL-6 can be induced in different cell types, including fibroblasts, monocytes, T cells, B cells, epidermal cells, and even in tumor cells (Yamashita et al., 1994). Because fibroblasts are the main cell type responsible for wound healing (including inflammation, re-forming the epithelial matrix and remodeling), these cells play an important role in skin inflammation, and as such, are a good study model to determine inflammatory patterns in dogs (Tsuchiya et al., 2015).

Isolation of dog primary fibroblasts

We isolated primary fibroblast cells from recently born and deceased dogs of two size classes. The small breed size class was composed of breeds with an adult body mass of 15 kg or less, and the large breed size class included breeds or mixes with an adult body mass of 20 kg or more These size classes are based on American Kennel club (AKC) standards of each breed, and described in Jimenez (2016). Small breed recently born puppies (N=45) were on average 3.4 days old at the time of sample collection. Small breed deceased dogs (N=12) were on average 13.7 years old at the time of sample collection. Large breed puppies (N=74) were on average 9.6 days old at the time of sample collection. Large breed deceased dogs (N=15) were on average 11.9 years old at the time of sample collection. Samples from recently born puppies were obtained from routine tail docks, ear clips and dew claw removals performed at veterinarian offices in Central New York and Michigan. Samples from deceased dogs were collected from ear clips immediately after euthanasia in Central New York. The samples were placed in cold transfer medium [Dulbecco's modified Eagle medium (DMEM), with 4.5 g l−1 glucose, sodium pyruvate, and 4 mmol l−1l-glutamine supplemented with 10% heat-inactivated fetal bovine serum and antibiotics (100 U ml−1 penicillin/streptomycin], containing 10 mmol l−1 HEPES) and transferred to Colgate University. To isolate primary fibroblast cells, skin samples were sterilized in 20% bleach and 70% ethanol. Once any fat and bone were removed, skin was minced and incubated in sterile 0.5% collagenase type 2 (Worthington Chemicals, cat. no. LS004176) overnight in an atmosphere of 37°C, 5% CO2 and 5% O2. After incubation, the collagenase mixture was filtered through a 20 μm sterile mesh, and centrifuged at 1000 rpm for 5 min. The resulting supernatant was removed, and the pellet was resuspended with 7 ml mammal medium [DMEM with 4.5 g l−1 glucose, sodium pyruvate, and 4 mmol l−1l-glutamine supplemented with 10% heat-inactivated fetal bovine serum and antibiotics (100 U ml−1 penicillin/streptomycin)]. Cells were grown in Corning T-25 culture flasks at 37°C in an atmosphere of 5% O2 and 5% CO2. When cells reached 90% confluence, they were trypsinized (0.25%) and cryopreserved at 106 cells ml−1 in DMEM supplemented with 40% fetal bovine serum and dimethylsulfoxide (DMSO) at a final concentration of 10%. We stored cells in liquid N2 prior to any experiments. All the procedures within this study were approved by the Colgate University Institutional Care and Use Committee under protocol number 1819-13.

Treatments

Each cell line was resuspended in 7 ml of complete medium (as above) and given 5 days to recover from freezing in an atmosphere of 37°C, 5% O2 and 5% CO2. Cells were plated at passage 2 (P2) in a concentration of 10,000 cells per well on a seahorse microculture plate. Cells from each dog were plated in duplicate per treatment (control or LPS). Treated cells were incubated with an LPS from E. coli (O111:B4, Sigma Aldrich, cat. no. L2630) at a concentration of 30 µg ml−1 for 24 h (Wendell and Stein, 2001; Tabeta et al., 2000), whereas control wells were not exposed to LPS. Immediately following the LPS incubation, oxygen consumption rate (OCR; aerobic metabolic rate) and extracellular acidification rates (ECAR; anaerobic/glycolytic metabolic rate) were performed as described in Patton et al. (2018), and below.

Oxygen consumption rate (OCR) in primary fibroblast cells

OCR was determined using XF-96 FluxPaks from Agilent Technologies. We measured OCRs after cells were equilibrated for 1 h to running medium, which contained 10 mmol l−1 glucose, 1 mmol l−1 sodium pyruvate and 2 mmol l−1 glutamine, pH=7.4. Baseline measurements of OCRs were made three times to establish basal OCR prior to injecting a final well concentration of 2 μmol l−1 oligomycin, which inhibits ATP synthesis by blocking the proton channel of the F0 portion of the ATP synthase, causing OCR to fall. Thus, the decrease in OCR from basal levels represents ATP coupled respiration, whereas the remaining OCR is attributed to O2 consumption required to overcome the natural proton leak across the inner mitochondrial membrane plus any non-mitochondrial O2 consumption. We then injected a well concentration of 0.125 μmol l−1 carbonyl cyanide-4-(trifluoromethoxy)phenylhydrazone (FCCP), an uncoupling agent that disrupts ATP synthesis by collapsing the proton gradient across the mitochondrial membrane leading to uncoupled consumption of energy and O2 without generating ATP, providing a theoretical maximal OCR. Finally, we injected a final well concentration of 0.5 μmol l−1 antimycin A, a complex III inhibitor and rotenone, a complex I inhibitor. This combination stops mitochondrial respiration and enables non-mitochondrial respiration to be evaluated (Gerencser et al., 2009; Brand and Nicholls, 2011; Hill et al., 2012). In addition to the OCR parameters we measured directly (as listed above), we also calculated spare respiratory capacity as the difference between maximal OCR and basal OCR, and ATP coupling efficiency as the percentage of basal OCR accounted for by ATP-coupled respiration (Divakaruni et al., 2014). After measurements were completed, we used a 1:200 concentration of CyQUANT dye to quantify final counts of cells in each well (Jones et al., 2001) and normalized all rates to 20,000 cells.

Extracellular acidification rate (ECAR) in primary fibroblast cells

ECAR values were measured in units of mpH, which is the change in pH in the medium surrounding the cells due to proton flux in glycolysis. Measurements of ECAR were performed after the cells were equilibrated to running medium for 1 h. Running medium contained no glucose and 2 mmol l−1l-glutamine in all experiments, pH=7.4. Baseline rates were measured three times prior to any injections. We injected a final well concentration of 10 mmol l−1 glucose into medium surrounding cells, which provides a measure of glycolytic rate, and then we injected a final well concentration of 2 µmol l−1 oligomycin, giving us an estimate of glycolytic capacity in cells. Finally, we injected a final well concentration of 50 mmol l−1 2-DG, a glucose analog that inhibits glycolysis, which provided an estimate of non-glycolytic acidification (Hill et al., 2012). After measurements were completed, we used a 1:200 concentration of CyQUANT dye to quantify final counts of cells in each well (Jones et al., 2001) and normalized all rates to 20,000 cells.

IL-6 assay

Immediately after 24 h of treatment, supernatants from control and LPS-conditioned cells were collected from each individual well (duplicated per dog as stated above) to assess the secretory profile of the pro-inflammatory cytokine IL-6 and were frozen at −20°C until further use. We used a commercially available kit (Canine IL-6 ELISA kit RayBiotech) to quantify IL-6 in 100 µl of supernatants following the manufacturer's protocol. We used a 1:200 concentration of CyQUANT dye to quantify final counts of cells in each well (Jones et al., 2001) and normalized all IL-6 concentrations to 20,000 cells. Intra-assay variability was 11.2–12.5%, whereas inter-assay variability was 13.1%

Statistics

Data from every assay were first tested for normality using a Shapiro–Wilks test. The data were transformed in a manner informed by a Box–Cox power transformation (Box and Cox, 1964). This meant a log transform for all variables except IL-6, which got an inverse-square transformation, i.e. (1/Y2). First, we analyzed all data in a single repeated measures MANOVA to test for differences in the collection of dependent variables across body size (small or large, as stated above), sample type (recently born or deceased) and treatment (control or LPS), their two-way interactions and a three-way interaction. Follow-up analyses were conducted using three-way repeated measures ANOVAs on each dependent variable. For dependent variables that suggested a treatment effect, we fitted mixed effects models using the lmerTest package (Kuznetsova et al., 2017) for R (https://www.r-project.org/) to quantify the effect. When conducting the analyses, we accounted for the repeated measures design by including a random intercept for subjects. Best-subset model selection according to the Bayesian information criterion (BIC) (Schwarz, 1978) was completed via complete enumeration (Morgan and Tatar, 1972) for each model. BIC penalizes false positives more than false negatives when there are eight or more observations, and is asymptotically consistent. We adjusted the P-values in the MANOVA analyses and within each model reported using the approach in Benjamini and Hochberg (1995) rather than using the overly conservative Bonferroni adjustment (Dunn, 1961). This approach allows for the control of the expected proportion of false discoveries instead of the probability of making at least one false discovery, thus preserving power. In the statistical models (see full details at https://figshare.com/s/ee6a5f273a542b703857), we included categorical size, which allowed us to consider the fact that recently born dogs would have a disproportionally small body mass, but could be part of the large size class, noting that this also better matches the goal of the experiment. Post hoc testing involved estimating marginal means for specified values of independent variables and comparing these means via statistical contrasts. We performed these post hoc tests via the emmeans package for R (https://cran.R-project.org/package=emmeans). Results were considered significant if the adjusted P-value was less than 0.05.

Limitations

Owing to the nature of working with pet dogs, our sample collection was limited in several ways. First, the number of recently born dogs per breed or group was limited to those breeds that are altered as a breed requirement (tail docks and dew claw removal). The number of euthanized dogs was limited owing to the owners' willingness to provide us ear clips at the time of euthanasia. This collection design also yielded a gap in observations between 0.20 and 6 years of age, hence using sample type as a surrogate for age; we note that this diminishes statistical power to detect the effect of age but we feel it better represents the data available. Secondly, we were not given complete medical charts for any dogs included in this study, thus, we have no records of diet or exercise, which may be variants in OCR, ECAR and background inflammation parameters. At times, veterinarians did not provide sex or weight for recently born dogs either. We do know that none of the dogs included in this study were taking any metabolic, endocrine or nervous system medications, and none were obese. Furthermore, young puppies may not have a fully developed immune system, and older dogs may have a declining immune system and underlying systemic diseases. However, we used primary fibroblast cells to elicit and measure the capabilities of mounting a response at the cellular level.

Previously published work (Jimenez et al., 2018) describes differences associated with size class and age class (sample type). Here, we focused on the differences between control and LPS-treated fibroblasts, as well as the interactive effects of size and sample type with treatment. Full details of statistical models and all MANOVA/ANOVA results can be found on Figshare at: https://figshare.com/s/ee6a5f273a542b703857.

The repeated measures MANOVA showed significant differences for the multivariate analysis in sample type (F=3.4030, P=0.0019), size class (F=4.7953, P=0.0001) and treatment (F=3.4842, P=0.0019) (see table 1 in https://figshare.com/s/ee6a5f273a542b703857). Additionally, there was a statistically significant Sample type×Treatment interaction (F=2.6885, P=0.0109). Thus, we probed each dependent variable separately using repeated measures ANOVA.

Results showed significant unique effects of treatment on non-mitochondrial respiration (P=0.0371; Table 1) and IL-6 concentrations (P<0.0001; Fig. 1 and Table 1). Additionally, the treatment effect did not reach traditional levels of statistical significance for non-glycolytic acidification, total glycolysis, total glycolytic capacity, basal OCR, proton leak, maximal respiration, or spare respiratory capacity through the unique or interaction effects (Table 1).

Full transformed model

Full models were fitted for non-mitochondrial respiration and IL-6 to follow up the statistically significant effect of treatment in the ANOVA analyses for these dependent variables. The non-mitochondrial respiration model required a log transformation and the IL-6 model required an inverse square transformation to ameliorate concerns with assumptions of the linear mixed effects model. Despite significant results in the three-way ANOVA, the treatment effect did not reach traditional levels of statistical significance for non-mitochondrial respiration through the unique effect of treatment (P=0.1298), the Treatment×Sample type interaction (P=0.9657), the Treatment×Size interaction (P=0.8892) or through the Sample type×Size×Treatment three-way interaction (P=0.6822) (Table 2). There was, however, a significant unique treatment effect on IL-6 (P<0.0001), and a significant interaction is the Sample type×Size (P=0.0271), whereas Treatment×Sample type (P=0.0880), Treatment×Size (P=0.7374) and Sample type×Size×Treatment (P=0.1090) interactions did not reach traditional levels of statistical significance (α=0.05). See ‘Post hoc Analysis of the Transformed Model’ in https://figshare.com/s/ee6a5f273a542b703857 for full details.

Post hoc analysis of the full transformed model for IL-6

Control (background) values of IL-6 showed that fibroblasts from recently born small breed dogs had a significantly higher concentration of IL-6 compared with those from young large breed dogs (P=0.0065). After LPS treatment, we found differences between cells from young and old large breed dogs (P=0.0141), young and old small breed dogs (P=0.0112) and large and small breed recently born dogs (P=0.0021). When treated with LPS, fibroblasts from young small breed dogs (P<0.0001), young large breed dogs (P<0.0001) and old small breed dogs (P=0.0002), but not old large breed dogs (P=0.6028) had higher IL-6 concentrations compared with untreated cells.

Furthermore, cells from old large breed dogs showed a less pronounced treatment effect compared with those from young large breed (P<0.0001), old small breed (P=0.0288) and young small breed (P<0.0001) dogs. There was also a larger treatment effect for cells from old small breed versus young large breed (P<0.0001) dogs. For full details of all these comparisons, see ‘Post-hoc Analysis of the Transformed Model’ section in: https://figshare.com/s/ee6a5f273a542b703857.

Best subsets model

When selecting independent variables using best subsets regression according to BIC, we identified the best model following hierarchy without including redundant or unnecessary terms. These models supported a 27.83% increase in non-mitochondrial respiration (P=0.0002) and increased IL-6 (P<0.0001) for LPS-treated compared with untreated fibroblasts.

Immunosenescence is well documented in dogs (Day, 2010; Alexander et al., 2018); however, only a handful of studies have investigated pro-inflammatory pathways in older dogs (Day, 2010). To our knowledge, this is the first study that examines the effects of body size, as well as age class in terms of cellular metabolism and response to a pro-inflammatory challenge in the domestic dog system. Constitutive components of the innate immune system are responsible for producing inflammatory cytokines resulting from local inflammation. At the organismal level, this upregulation of the inflammatory response is marked by increases in acute phase proteins in the liver, and costly changes in energy metabolism (Lee, 2006). LPS also stimulates an increase in synthesis of growth factors and cytokines at the cellular level (Pang et al., 1994). Our results do not support our predictions listed in the Introduction, as we found no differences in cellular metabolism with LPS treatment, across body size and age classes (sample type), with the exception of non-mitochondrial respiration. We also found the opposite pattern to that we predicted for IL-6 concentration. Others have found that decreased mitochondrial membrane uncoupling found in large breed dog primary fibroblast cells, coupled with lower potential of β-oxidation and accumulation of acylcarnitines, promotes inflammation (Nicholatos et al., 2019), which follows the pattern in our current study that young large breed dogs seem to have a higher IL-6 concentration.

Life history theory predicts that fast- and slow-paced species should invest differently in inflammatory responses because their cost–benefit structure differs (Lee, 2006). Slow-paced species live longer than fast-paced species and thus benefit more from responding to infections in a manner that limits immunopathology, whereas fast-paced species benefit from clearing infections and maximizing current reproductive output regardless of the long-term costs of the immune response. It follows that species on the ‘fast’ end of the life-history spectrum should have proinflammatory responses that cause more immunopathology, whereas species on the ‘slow’ end should invest in anti-inflammatory responses that reduce immunopathology (Sears et al., 2011). Determining life-history differences in domestic dogs is not clear-cut, as small dog breeds live longer, and have fewer pups per litter – both determinants of a slower pace of life – but also seem to have faster growth rates compared with large breed dogs (Jimenez, 2016) like many faster pace of life animals. IL-6 background concentration showed that recently born small breed dogs had significantly higher IL-6 background concentration compared with recently born large breed dogs, and IL-6 background (control) concentration is similar in deceased small breed dogs, and deceased large breed dogs. However, our sample size for large breed old dogs may have limited these interpretations owing to a lack of statistical power. This pattern of higher IL-6 background concentrations may point to an anti-inflammatory investment strategy in small dogs at the ‘slow’ end of the life history spectrum.

Similar to findings included in this study, other studies have found differences between small and large size classes with respect to glycolytic parameters in dogs, where larger breed puppies have significantly higher glycolytic capacity compared with smaller breeds (Jimenez et al., 2018), although it seems this pattern is related to increased glucose oxidation (Brookes and Jimenez, 2021). LPS treatment has been shown to increase lactate concentrations in the acute monocytic leukaemia cell line compared with controls (Ubanako et al., 2019). Intracellular inflammatory mediators, such as TNF-α directly inhibit insulin signaling (Odegaard and Chawla, 2013), which would alter aerobic metabolism. Neither of these were applicable in our study, as we did not see any changes in cellular metabolic rates with LPS. However, non-mitochondrial oxygen consumption includes oxygen-consuming oxidases functioning throughout metabolism, such as peroxisomal fatty acid oxidation (Rolfe and Brown, 1997), and we saw an increase in non-mitochondrial respiration in LPS treated cells compared with controls. This could be due to increases in NADH oxidases that seem to be linked to increases in IL-6 concentration (Behrens et al., 2008; Yu et al., 2005).

In dogs, there is little information about the role of IL-6 during the acute phase response due to inflammatory stimuli (Yamashita et al., 1994). In beagle dogs injected with turpentine oil to induce an inflammatory response, IL-6 concentration in blood increased by 4 h after the injection until reaching a maximum at 12 h post injection, and remained higher than levels in the controls up to 6 days after the injection (Yamashita et al., 1994). Stimulated LPS responses in whole blood of aging dogs demonstrated a significant increase in IL-6 concentrations compared with control values, regardless of age, and although their study included dogs of different breeds, there seemed to be no correction of the data for body size (Deitschel et al., 2010). Different dog tissues have demonstrated different abilities to upregulate IL-6 concentrations with peritoneal cells having a more exaggerated response compared with peripheral blood mononuclear cells (PBMCs) isolated from fox terriers and Labrador retrievers, although age did not seem to affect IL-6 concentrations in either tissue or breed (Kearns et al., 1999). Our results are similar to those found in other studies of dogs (Deitschel et al., 2010). Unlike geriatric mice and humans (Tateda et al., 1996; Gabriel et al., 2002), the adult population of dogs exhibited a less drastic amplification of IL-6 concentration after LPS treatment. Rather, cells from younger dogs of both small and large size classes demonstrated a drastic upregulation of IL-6 when treated with LPS. In contrast, background IL-6 data from our study does follow the patterns of geriatric humans and mice when considering just the large dog size class, such that fibroblasts from geriatric large dogs have an increase in background IL-6 compared with cells from their younger counterparts. In humans, IL-6 concentrations are undetectable in young people, but start to increase at age 50–60 and continue to increase so that even healthy centenarians have high levels of IL-6; however, increases in IL-6 in humans have also been demonstrated as predictors of morbidities and mortality in elderly people and are linked to chronic inflammatory diseases (Franceschi et al., 2000; Maggio et al., 2006). Although the main source of inflammation associated with chronic disease and aging is leukocytes, particularly macrophages and monocytes (Sarkar and Fisher, 2006). In Labrador retrievers, C-reactive protein (CRP), which is produced by the liver in response to IL-6 and TNF-α, increased with age (Alexander et al., 2018), similar to our results in old large dogs.

Previous work suggests that DNA oxidative damage (8-OHdg) negatively correlates with mean breed lifespan, so that shorter lived breeds may accumulate more oxidative DNA damage compared with longer lived breeds (Jimenez et al., 2018). Others have also found that plasma 8-OHdg levels increase with age in Labrador retrievers (Alexander et al., 2018). Damaged DNA in the nucleus including DNA double strand breaks (DSBs) commonly accumulated as oxidative damage (Hasty et al., 2003), can stimulate proinflammatory responses (Brzostek-Racine et al., 2011). Acutely, these responses are initiated through communication between the nucleus and the cytoplasm to include activation of nuclear factors, cell surface ligands or intercellular adhesion molecules, leading to an upregulation of pro-inflammatory cytokines (Chatzinikolaou et al., 2014). Chronically, DNA damage due to age-accumulated cellular senescence leads to a phenotype that entails the secretion of cytokines, such as IL-6 and IL-8, promoting chronic inflammation (Chatzinikolaou et al., 2014). Cell culture work with human cells suggests that an increase in oxidative stress causes the activation of nuclear factor-kappa B (NF-κB), a transcription factor that upregulates the synthesis of proinflammatory cytokines, including IL-6, thus linking higher IL-6 levels with chronic inflammation and aging mechanistically (Sarkar and Fisher, 2006). The advantage of this phenotype is that DNA-damage driven inflammation enables angiogenesis, which promotes the invasion and metastasis of tumor cells (Chatzinikolaou et al., 2014). In contrast, a continued accumulation of cellular damage, leading to increases in proinflammatory pathways over time, leads to chronic inflammation, tissue malfunction and degeneration in old age (Chatzinikolaou et al., 2014). Considering that small breed dogs are usually longer lived compared with larger breeds, and that we expect more of an increase in IL-6 for small dogs compared with large dogs, it seems that the mechanism or pathway for increases in IL-6 is not related to stimulation of proinflammatory responses due to DNA oxidative damage accumulation. Alternatively, if inflammation is a contributor to somatic damage associated with aging, the lower IL-6 in small dogs might facilitate a longer lifespan.

Cytokines produced by fibroblasts support the growth, differentiation and activation of inflammatory cells (Vancheri et al., 1991), suggesting that they are important in immune defenses. The lack of metabolic dysregulation or stimulation in conjunction with the observed increase in IL-6 is intriguing. When fibroblasts from clinically normal duodenum tissues of humans were stimulated with LPS, those cells synthesized new IL-6, IL-8 and granulocyte-macrophage colony-stimulating factor (Pang et al., 1994), rather than releasing stored cytokines. If the fibroblasts used in this study also synthesize new cytokines, it follows that the cells must be downregulating other processes to support the production of IL-6 without also increasing ATP production.

Here, we have demonstrated that most of cellular metabolism in primary fibroblasts does not change significantly with age, size or LPS treatment. However, we did see differences in IL-6 concentration in young animals whereby cells from small breed dogs had higher IL-6 concentrations than those from large breed dogs. We also found that IL-6 concentration increased more dramatically with aging after LPS treatment in small breed compared with large breed dogs. This implies that the patterns of inflammaging in dogs may be distinct and different from other mammals commonly studied.

We are grateful to the following veterinarians and veterinary practices for providing us with samples: Dr Kerri Hudson, Dr James Gilchrist, Dr Heather Culbertson and Morgan Peppenelli at Waterville Veterinary Clinic (New York); Dr Frank Capella from Village Vet in Wampsville, NY. Pet Street Station Animal Hospital (New York); Dr Jim Bader at Mapleview Animal Hospital (Michigan). We are also grateful to the following breeders for participating in our study: Rhonda Poe, Bob Stauffer, Allison Mitchell, Nancy Secrist, Valeria Rickard, Joanne Manning, Lita Long, Betsy Geertson, Susan Banovic, Lisa Uhrich, Sheryl Beitch, Al Farrier, Barbara Hoopes and Rachel Sann.

Author contributions

Conceptualization: A.G.J.; Methodology: A.G.J., C.J.D., W.C.; Validation: A.G.J., W.C.; Formal analysis: S.L., W.C.; Investigation: A.G.J.; Resources: A.G.J., C.J.D., S.L., W.C.; Data curation: S.L., W.C.; Writing - original draft: A.G.J.; Writing - review & editing: C.J.D., S.L., W.C.; Supervision: W.C.; Project administration: A.G.J., W.C.; Funding acquisition: A.G.J., C.J.D.

Funding

The Seahorse XF96e oxygen flux analyzer was purchased via a National Science Foundation Major Research Instrument grant (NSF MRI 1725841 to A.G.J.). A Research Council grant from Colgate University to A.G.J. partly funded this work. Additionally, this research was supported by a National Science Foundation grant to C.J.D. (NSF IOS 1656551) and Hamilton College’s Dean of Faculty.

Data availability

Full statistical models and data analyses are available on Figshare at: https://figshare.com/s/ee6a5f273a542b703857.

Agarwal
,
S.
,
Baran
,
C.
,
Piesco
,
N. P.
,
Quintero
,
J. C.
,
Langkamp
,
H. H.
,
Johns
,
L. P.
and
Chandra
,
C. S.
(
1995
).
Synthesis of proinflammatory cytokines by human gingival fibroblasts in response to lipopolysaccharides and interleukin–1β
.
J. Periodontal Res.
30
,
382
-
389
.
Alexander
,
J. E.
,
Colyer
,
A.
,
Haydock
,
R. M.
,
Hayek
,
M. G.
and
Park
,
J.
(
2018
).
Understanding how dogs age: longitudinal analysis of markers of inflammation, immune function, and oxidative stress
.
J. Gerontol. A
73
,
720
-
728
.
Baylis
,
D.
,
Bartlett
,
D. B.
,
Patel
,
H. P.
and
Roberts
,
H. C.
(
2013
).
Understanding how we age: insights into inflammaging
.
Longev. Healthspan
2
,
8
.
Behrens
,
M. M.
,
Ali
,
S. S.
and
Dugan
,
L. L.
(
2008
).
Interleukin-6 mediates the increase in NADPH-oxidase in the ketamine model of schizophrenia
.
J. Neurosci.
28
,
13957
-
13966
.
Benjamini
,
Y.
and
Hochberg
,
Y.
(
1995
).
Controlling the false discovery rate: a practical and powerful approach to multiple testing
.
J. R. Stat. Soc. B
57
,
289
-
300
.
Brand
,
M. D.
and
Nicholls
,
D. G.
(
2011
).
Assessing mitochondrial dysfunction in cells
.
Biochem. J.
435
,
297
-
312
.
Brookes
,
P. S.
and
Jimenez
,
A. G.
(
2021
).
Metabolomics of aging in primary fibroblasts from small and large breed dogs
.
bioRxiv
.
Cohen
,
A. A.
(
2018
).
Aging across the tree of life: the importance of a comparative perspective for the use of animal models in aging
.
Biochim. Biophys. Acta Mol. Basis Dis.
1864
,
2680
-
2689
.
Cossarizza
,
A.
,
Ortolani
,
C.
,
Monti
,
D.
and
Franceschi
,
C.
(
1997
).
Cytometric analysis of immunosenescence
.
Cytometry
27
,
297
-
313
.
Dantzer
,
B.
and
Swanson
,
E. M.
(
2012
).
Mediation of vertebrate life histories via insulin–like growth factor–1
.
Biol. Rev.
87
,
414
-
429
.
Day
,
M. J.
(
2010
).
Ageing, immunosenescence and inflammageing in the dog and cat
.
J. Comp. Pathol.
142
,
S60
-
S69
.
Deitschel
,
S. J.
,
Kerl
,
M. E.
,
Chang
,
C. H.
and
DeClue
,
A. E.
(
2010
).
Age–associated changes to pathogen–associated molecular pattern–induced inflammatory mediator production in dogs
.
J. Vet. Emerg. Crit. Care
20
,
494
-
502
.
Dingli
,
D.
and
Pacheco
,
J. M.
(
2006
).
Allometric scaling of the active hematopoietic stem cell pool across mammals
.
PLoS ONE
1
,
e2
.
Divakaruni
,
A. S.
,
Paradyse
,
A.
,
Ferrick
,
D. A.
,
Murphy
,
A. N.
and
Jastroch
,
M.
(
2014
).
Analysis and interpretation of microplate-based oxygen consumption and pH data
. In
Methods in Enzymology
, vol.
547
, pp.
309
-
354
.
Academic Press
.
Downs
,
C. J.
,
Dochtermann
,
N. A.
,
Ball
,
R.
,
Klasing
,
K. C.
and
Martin
,
L. B.
(
2020
).
The effects of body mass on immune cell concentrations of mammals
.
Am. Nat.
195
,
107
-
114
.
Dunn
,
O. J.
(
1961
).
Multiple comparisons among means
.
J. Am. Stat. Assoc.
56
,
52
-
64
.
Eigenmann
,
J. E.
,
Amador
,
A.
and
Patterson
,
D. F.
(
1988
).
Insulin-like growth factor I levels in proportionate dogs, chondrodystrophic dogs and in giant dogs
.
Eur. J. Endocrinol.
118
,
105
-
108
.
Fontana
,
L.
,
Partridge
,
L.
and
Longo
,
V. D.
(
2010
).
Extending healthy life span—from yeast to humans
.
Science
328
,
321
-
326
.
Franceschi
,
C.
,
Bonafè
,
M.
,
Valensin
,
S.
,
Olivieri
,
F.
,
De Luca
,
M.
,
Ottaviani
,
E.
and
De Benedictis
,
G.
(
2000
).
Inflamm–aging: an evolutionary perspective on immunosenescence
.
Ann. N. Y. Acad. Sci.
908
,
244
-
254
.
Franceschi
,
C.
and
Campisi
,
J.
(
2014
).
Chronic inflammation (inflammaging) and its potential contribution to age-associated diseases
.
J. Gerontol. A Biomed. Sci. Med. Sci.
69
(Suppl_1)
,
S4
-
S9
.
Gabriel
,
P.
,
Cakman
,
I.
and
Rink
,
L.
(
2002
).
Overproduction of monokines by leukocytes after stimulation with lipopolysaccharide in the elderly
.
Exp. Gerontol.
37
,
235
-
247
.
Greeley
,
E. H.
,
Kealy
,
R. D.
,
Ballam
,
J. M.
,
Lawler
,
D. F.
and
Segre
,
M.
(
1996
).
The influence of age on the canine immune system
.
Vet. Immunol. Immunopathol.
55
,
1
-
10
.
Gerencser
,
A. A.
,
Neilson
,
A.
,
Choi
,
S. W.
,
Edman
,
U.
,
Yadava
,
N.
,
Oh
,
R. J.
,
Ferrick
,
D. A.
,
Nicholls
,
D. G.
and
Brand
,
M. D.
(
2009
).
Quantitative microplate-based respirometry with correction for oxygen diffusion
.
Anal. Chem.
81
,
6868
-
6878
.
Greer
,
K. A.
,
Canterberry
,
S. C.
and
Murphy
,
K. E.
(
2007
).
Statistical analysis regarding the effects of height and weight on life span of the domestic dog
.
Res. Vet. Sci.
82
,
208
-
214
.
Harrison
,
J. F.
(
2017
).
Do performance–safety tradeoffs cause hypometric metabolic scaling in animals?
Trends Ecol. Evol.
32
,
653
-
664
.
Hill
,
B. G.
,
Benavides
,
G. A.
,
Lancaster
,
J. R.
,
Ballinger
,
S.
,
Dell'Italia
,
L.
,
Zhang
,
J.
and
Darley-Usmar
,
V. M.
(
2012
).
Integration of cellular bioenergetics with mitochondrial quality control and autophagy
.
Biol. Chem.
393
,
1485
-
1512
.
Jimenez
,
A. G.
(
2016
).
Physiological underpinnings in life-history trade-offs in man's most popular selection experiment: the dog
.
J. Comp. Physiol. B
186
,
813
-
827
.
Jimenez
,
A. G.
,
Winward
,
J.
,
Beattie
,
U.
and
Cipolli
,
W.
(
2018
).
Cellular metabolism and oxidative stress as a possible determinant for longevity in small breed and large breed dogs
.
PLoS ONE
13
,
1
-
20
.
Jones
,
L. J.
,
Gray
,
M.
,
Yue
,
S. T.
,
Haugland
,
R. P.
and
Singer
,
V. L.
(
2001
).
Sensitive determination of cell number using the CyQUANT® cell proliferation assay
.
J. Immunol. Methods
254
,
85
-
98
.
Kearns
,
R. J.
,
Hayek
,
M. G.
,
Turek
,
J. J.
,
Meydani
,
M.
,
Burr
,
J. R.
,
Greene
,
R. J.
,
Marshall
,
C. A.
,
Adams
,
S. M.
,
Borgert
,
R. C.
and
Reinhart
,
G. A.
(
1999
).
Effect of age, breed and dietary omega-6 (n-6): omega-3 (n-3) fatty acid ratio on immune function, eicosanoid production, and lipid peroxidation in young and aged dogs
.
Vet. Immunol. Immunopathol.
69
,
165
-
183
.
Kirkwood
,
T. B.
(
2005
).
Understanding the odd science of aging
.
Cell
120
,
437
-
447
.
Kolaczkowska
,
E.
and
Kubes
,
P.
(
2013
).
Neutrophil recruitment and function in health and inflammation
.
Nat. Rev. Immunol.
13
,
159
-
175
.
Kuznetsova
,
A.
,
Brockhoff
,
P. B.
and
Christensen
,
R. H.
(
2017
).
lmerTest package: tests in linear mixed effects models
.
J. Stat. Softw.
82
,
1
-
26
.
Lee
,
K. A.
(
2006
).
Linking immune defenses and life history at the levels of the individual and the species
.
Integr. Comp. Biol.
46
,
1000
-
1015
.
López-Otín
,
C.
,
Blasco
,
M. A.
,
Partridge
,
L.
,
Serrano
,
M.
and
Kroemer
,
G.
(
2013
).
The hallmarks of aging
.
Cell
153
,
1194
-
1217
.
Maggio
,
M.
,
Guralnik
,
J. M.
,
Longo
,
D. L.
and
Ferrucci
,
L.
(
2006
).
Interleukin-6 in aging and chronic disease: a magnificent pathway
.
J. Gerontol. A Biol. Sci. Med. Sci.
61
,
575
-
584
.
Michell
,
A. R.
(
1999
).
Longevit of British breeds of dog and its relationships with-sex, size, cardiovascular variables and disease
.
Veterinary Record
145
,
625
-
629
.
Morgan
,
J. A.
and
Tatar
,
J. F.
(
1972
).
Calculation of the residual sum of squares for all possible regressions
.
Technometrics
14
,
317
-
325
.
Nicholatos
,
J. W.
,
Robinette
,
T. M.
,
Tata
,
S. V.
,
Yordy
,
J. D.
,
Francisco
,
A. B.
,
Platov
,
M.
,
Yeh
,
T. K.
,
Ilkayeva
,
O. R.
,
Huynh
,
F. K.
and
Dokukin
,
M.
(
2019
).
Cellular energetics and mitochondrial uncoupling in canine aging
.
GeroScience
41
,
229
-
242
.
Odegaard
,
J. I.
and
Chawla
,
A.
(
2013
).
The immune system as a sensor of the metabolic state
.
Immunity
38
,
644
-
654
.
Pang
,
G.
,
Couch
,
L.
,
Batey
,
R.
,
Clancy
,
R.
and
Cripps
,
A.
(
1994
).
GM–CSF, IL–1α, IL–β, IL–6, IL–8, IL–10, ICAM–1 and VCAM–1 gene expression and cytokine production in human duodenal fibroblasts stimulated with lipopolysaccharide, IL–1α and TNF–α
.
Clin. Exp. Immunol.
96
,
437
-
443
.
Partridge
,
L.
and
Gems
,
D.
(
2002
).
Mechanisms of aging: public or private?
Nat. Rev. Genet.
3
,
165
.
Patton
,
M. G.
,
Gillum
,
T. L.
,
Szymanski
,
M. C.
,
Gould
,
L. M.
,
Lauterbach
,
C. J.
,
Vaughan
,
R. A.
and
Kuennen
,
M. R.
(
2018
).
Heat acclimation increases mitochondrial respiration capacity of C2C12 myotubes and protects against LPS-mediated energy deficit
.
Cell Stress Chaperones
23
,
871
-
883
.
Promislow
,
D. E.
(
1993
).
On size and survival: progress and pitfalls in the allometry of life span
.
J. Gerontol.
48
,
B115
-
B123
.
Ricklefs
,
R. E.
(
2010
).
Life-history connections to rates of aging in terrestrial vertebrates
.
Proc. Natl Acad. Sci. USA
107
,
10314
-
10319
.
Sanz
,
A.
,
Pamplona
,
R.
and
Barja
,
G.
(
2006
).
Is the mitochondrial free radical theory of aging intact?
Antioxid Redox Signal.
8
,
582
-
599
.
Sarkar
,
D.
and
Fisher
,
P. B.
(
2006
).
Molecular mechanisms of aging-associated inflammation
.
Cancer Lett.
236
,
13
-
23
.
Schwarz
,
G.
(
1978
).
Estimating the dimension of a model
.
Ann. Stat.
6
,
461
-
464
.
Sears
,
B. F.
,
Rohr
,
J. R.
,
Allen
,
J. E.
and
Martin
,
L. B.
(
2011
).
The economy of inflammation: when is less more?
Trends Parasitol.
27
,
382
-
387
.
Strasser
,
A.
,
Teltscher
,
A.
,
May
,
B.
,
Sanders
,
C.
and
Niedermüller
,
H.
(
2000
).
Age–associated changes in the immune system of German shepherd dogs
.
J. Vet. Med. Ser. A
47
,
181
-
192
.
Sun
,
N.
,
Youle
,
R. J.
and
Finkel
,
T.
(
2016
).
The mitochondrial basis of aging
.
Mol. Cell
61
,
654
-
666
.
Sutter
,
N. B.
,
Bustamante
,
C. D.
,
Chase
,
K.
,
Gray
,
M. M.
,
Zhao
,
K.
,
Zhu
,
L.
,
B.
Padhukasahasram
,
E.
Karlins
,
S.
Davis
,
P. G.
Jones
et al. 
(
2007
).
A single IGF1 allele is a major determinant of small size in dogs
.
Science
316
,
112
-
115
.
Tabeta
,
K.
,
Yamazaki
,
K.
,
Akashi
,
S.
,
Miyake
,
K.
,
Kumada
,
H.
,
Umemoto
,
T.
Yoshie
,
H.
(
2000
).
Toll-like receptors confer responsiveness to lipopolysaccharide from Porphyromonas gingivalis in human gingival fibroblasts
.
Infect. Immun.
68
,
3731
-
3735
.
Tateda
,
K.
,
Matsumoto
,
T.
,
Miyazaki
,
S.
and
Yamaguchi
,
K.
(
1996
).
Lipopolysaccharide-induced lethality and cytokine production in aged mice
.
Infect. Immun.
64
,
769
-
774
.
Tian
,
X.
,
Seluanov
,
A.
and
Gorbunova
,
V.
(
2017
).
Molecular mechanisms determining lifespan in short-and long-lived species
.
Trend. Endocrinol. Metabol.
28
,
722
-
734
.
Tsuchiya
,
H.
,
Nakano
,
R.
,
Konno
,
T.
,
Okabayashi
,
K.
,
Narita
,
T.
and
Sugiya
,
H.
(
2015
).
Activation of MEK/ERK pathways through NF-κB activation is involved in interleukin-1β-induced cyclooxygenease-2 expression in canine dermal fibroblasts
.
Vet. Immunol. Immunopathol.
168
,
223
-
232
.
Ubanako
,
P.
,
Xelwa
,
N.
and
Ntwasa
,
M.
(
2019
).
LPS induces inflammatory chemokines via TLR-4 signalling and enhances the Warburg Effect in THP-1 cells
.
PLoS ONE
14
,
e0222614
.
Urfer
,
S. R.
,
Greer
,
K.
and
Wolf
,
N. S.
(
2011
).
Age-related cataract in dogs: a biomarker for life span and its relation to body size
.
Age
33
,
451
-
460
.
Vancheri
,
C.
,
Ohtoshi
,
T.
,
Cox
,
G.
,
Xaubet
,
A.
,
Abrams
,
J. S.
,
Gauldie
,
J.
,
Dolovich
,
J.
,
Denburg
,
J.
and
Jordana
,
M.
(
1991
).
Neutrophilic differentiation induced by human upper airway fibroblast-derived granulocyte/macrophage colony-stimulating factor (GM-CSF)
.
Am. J. Respir. Cell Mol. Biol.
4
,
11
-
17
.
Warburg
,
O.
(
1956
).
On respiratory impairment in cancer cells
.
Science (New York, NY)
124
,
269
-
270
.
Wendell
,
K. J.
and
Stein
,
S. H.
(
2001
).
Regulation of cytokine production in human gingival fibroblasts following treatment with nicotine and lipopolysaccharide
.
J. Periodontol.
72
,
1038
-
1044
.
Xia
,
S.
,
Zhang
,
X.
,
Zheng
,
S.
,
Khanabdali
,
R.
,
Kalionis
,
B.
,
Wu
,
J.
,
Wan
,
W.
and
Tai
,
X.
(
2016
).
An update on inflamm-aging: mechanisms, prevention, and treatment
.
J. Immunol. Res.
2016
,
8426874
.
Yamashita
,
K.
,
Fujinaga
,
T.
,
Miyamoto
,
T.
,
Hagio
,
M.
,
Izumisawa
,
Y.
and
Kotani
,
T.
(
1994
).
Canine acute phase response: relationship between serum cytokine activity and acute phase protein in dogs
.
J. Vet. Med. Sci.
56
,
487
-
492
.
Yu
,
J. H.
,
Lim
,
J. W.
,
Kim
,
H.
and
Kim
,
K. H.
(
2005
).
NADPH oxidase mediates interleukin-6 expression in cerulein-stimulated pancreatic acinar cells
.
Int. J. Biochem. Cell Biol.
37
,
1458
-
1469
.
Zimmerman
,
L. M.
,
Bowden
,
R. M.
and
Vogel
,
L. A.
(
2014
).
A vertebrate cytokine primer for eco–immunologists
.
Funct. Ecol.
28
,
1061
-
1073
.

Competing interests

The authors declare no competing or financial interests.