Effects of Air Travel Stress on the Canine Microbiome: A Pilot Study
Venable EB1*, Bland SD2, Holscher HD2, Swanson KS3
1 Animal Science Food and Nutrition; Southern Illinois University, Lincoln Drive, Carbondale, IL, USA.
2 Department of Food Science and Human Nutrition, University of Illinois, Urbana, IL. USA.
3 Department of Animal Sciences, University of Illinois, Urbana, IL, USA.
Erin B. Venable
Assistant Professor, Animal Science Food and Nutrition;
Southern Illinois University, 1205 Lincoln Drive, Carbondale, IL 62901, USA.
Tel: (618) 453-1358
Fax: (618) 453-2090
Received: March 22, 2016; Accepted: June 28, 2016; Published: July 01, 2016
Citation: Venable EB, Bland SD, Holscher HD, Swanson KS (2016) Effects of Air Travel Stress on the Canine Microbiome: A Pilot Study. Int J Vet Health Sci Res. 4(6), 132-139.DOI : dx.doi.org/10.19070/2332-2748-1600028
Copyright: Venable EB© 2016. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution and reproduction in any medium, provided the original author and source are credited.
Behavior assessment is commonly used to identify canines that may have a higher threshold for environmental stress, but no work has established a connection between behavior (as indicated by search performance) and travel stress (as indicated by fecal scores and microbial stability). Six canines (aged 18 months to 8 years), trained according to the standard established by the Federal Emergency Management Agency (FEMA), were utilized to test the effects of airline travel stress on working canines. Our objectives were to test the hypotheses that: 1) working canines can overcome air travel stress with little or no impact on their performance; and 2) fecal score and microbial composition is impacted by airline travel stress. Two groups of dogs, (n=3 per group), were randomly selected from FEMA canine teams rostered in New York City, NY (CONTROL) and in Miami, FL (TRAVEL). TRAVEL dogs were flown in the cabin on a commercial airline to New York. Blood and fecal samples were collected each morning prior to travel (d0) and search work (d1-d3). Fecal bacterial DNA was extracted and 16S rRNA amplicon sequencing was completed using an Illumina MiSeq followed by analysis with QIIME 1.8.0. Fecal scores from TRAVEL were significantly higher (P = 0.01) than CONTROL indicating softer stool in the group that travelled. Pre- and post-travel blood samples for the TRAVEL group were compared and demonstrated significant decreases in lactate, bicarbonate, total carbon dioxide, and base excess (P < 0.05) following travel. However, these decreases were still within normal range, therefore, may not be biologically significant. In addition, blood samples from TRAVEL and CONTROL were compared on search days and increases were observed in the TRAVEL group for ionized calcium, bicarbonate, total carbon dioxide, base excess, creatinine, and blood urea nitrogen (P < 0.05). In contrast, blood glucose was decreased in the TRAVEL group (P < 0.05). Search behavior scores were not significant in TRAVEL compared to CONTROL. Principal coordinates analysis (PCoA) of UniFrac distances between samples based on their 97% OTU composition indicated that TRAVEL bacterial communities (P = 0.01) and bacterial community abundances (P = 0.02) were significantly different from CONTROL. These data demonstrate that airline travel of 2.5h impacts the working canine gut microbiota and some blood metabolites, but has no observed effect on working canine performance.
3.Materials and Methods
3.1.Animals and Diets
3.3.Blood and Fecal Assays
Microbiome; Travel Stress; Working Canine.
Stress is defined as a state during which the organism reacts to endogenous and exogenous threats and focuses its energies on coping with the situation of danger . Stress can cause a wide variety of physiological and psychological health issues in canines, including increased activity of the sympathetic nervous system, increased release of catecholamine, and increased blood pressure . The physiological effects of stress on the dog are poorly understood and therefore most likely under-reported. It is important to understand the level of stress working dogs can handle before experiencing physiological consequences in order to prevent negative impacts on job performance.
Working canines are often called upon to travel by air as part of their mission. However, no studies have been conducted identifying the effects of travel on working canines’ physiological state and performance ability. Stress can cause a wide variety of physiological and psychological health issues in canines, including increased activity of the sympathetic nervous system,catecholamine, blood pressure, and permeability of the intestinal epithelial lining to microorganisms [7, 28]. In addition, a past study demonstrated that gastrointestinal permeability is increased as part of the stress response which may be correlated to acute stress symptoms such as diarrhea . This may result in a change in fecal score. The only data available to assess how working canines handle the stress of air travel is extrapolated from companion animal work. According to a study conducted by Bergeron et al. (2002) , both ground and air transportation caused significantly increased levels of stress in dogs. The second most common cause of animal deaths associated with air transportation are secondary illnesses that are thought to be triggered by the environmental stress associated with travel. This travel was followed by acute onset of diseases caused by stress and high altitude, and also may have been due to mishandling of the animals in their carriers .
In contrast, search and rescue (SAR) canines have been conditioned to work around and through stressful situations . However, more investigation is warranted regarding the physiological impacts of stress on working canines’ performance. Prior research has shown a variety of gastrointestinal issues associated with working canines during deployment [3, 24]. However, no information is available that has connected issues associated with gastrointestinal upset and performance ability. The objectives of the current investigation were to test two hypotheses: 1) fecal scores and fecal microbial composition of working canines are impacted by airline travel stress; 2) working canines can overcome this stress with little or no impact on their performance in the field.
Materials and Methods
Southern Illinois University Institutional Animal Care and Use approval was obtained prior to the initiation of this pilot study (#14-047). Standards for animal care were adopted from previously published recommendations . Six canines (two German Shepherds, one Belgian Malinois, two Labrador Retrievers, one mixed breed), aged 5 years ± 3 with a body condition score (BCS) of 4.5 ± 0.5, were randomly selected from two Federal Emergency Management Agency (FEMA) teams and fed commercially available complete and balanced dry kibble diets. Subjects used in this study and for FEMA work were screened according to similar criteria and were trained according to the established federal standard . Dogs were divided into two groups (TRAVEL and CONTROL). Canines in the CONTROL group were housed individually in New York, NY, in temperature-and light-controlled kennels with ad libitum access to water. Canines in the TRAVEL group were located in Miami, FL and flew to New York as part of the experimental procedure. Canines in the TRAVEL group were offered water prior to departure and upon landing and were housed individually overnight in temperature and light-controlled hotel rooms. All canines were fed twice daily at 0600 and 1700h. Canines in the TRAVEL group were screened through standard commercial airline screening procedures (www.tsa.gov/traveler-information) and were flown in the cabin (2.5h flight time). Dogs were transported to the search site and were housed in temperature-controlled kennels (0.962m x 0.635m x 0.685m) at the New York Police Department training site for the duration of the search scenarios. All TRAVEL dogs returned to their hotels each night. All CONTROL dogs returned to their individual kennels each night.
The experiment was conducted in October of 2013, with ambient air temperature of 23.8 ± 2.1°C and relative humidity of 55.4 ± 9.2%, with search behavior scoring taking place on the second and third day of the four day study. The search site was unfamiliar to the TRAVEL dogs, but was in compliance with rubble pile standards established and published by FEMA . The principal investigator and FEMA instructors assigned ranked behavioral scores during the search performance for each subject using an adaptation of a previously published behavior scoring system [17, 18] and by utilizing focal point sampling.
Samples were collected from subject seach day prior to flight (d0) and prior to search (d1 to d3). Parameters collected included blood, rectal temperature, heart rate, and morning fecal samples (Purina Fecal Scoring System, Nestle-Purina, St. Louis, MO). Whole blood was used to monitor blood metabolites, including lactate, pH, partial pressure carbon dioxide, partial pressure oxygen, bicarbonate, sodium, potassium, ionized calcium, hematocrit, base excess, hemoglobin, and oxygen saturation with an Abaxis iSTAT (Union City, CA). Remaining blood samples were collected in EDTA tubes and stored on ice until centrifuged(1,300×g, 15min, 21°C) for serum analysis (PowerSpin LX™, Dayton, NJ). Serum samples were analyzed for liver function, alanine aminotransferase (ALT), kidney function, blood urea nitrogen (BUN), total protein, and creatinine, and blood glucose concentrations with an AbaxisVetScan V2 (Union City, California).
Fecal samples were collected in sterile 15mL vials, stored on ice within 15min of defecation, and frozen overnight prior to overnight shipment for laboratory analysis. Upon arriving at the laboratory, samples were stored at -80°C. Fecal bacterial DNA was extracted using MoBio PowerSoil DNA extraction kits (MoBio Laboratories, Carlsbad, CA) and quantified using a NanoDrop (ND-1000 spectrophotometer, (Nanodrop Technologies) and diluted to 10 ng/mL. Electrophoresis using 1% agarose gel, with a 1xTAE buffer at 170V for 60min, evaluated the quality of the extracted DNA (BioRad Cell-Sub, Hercules, CA). A Fluidigm Access Array was used to generate barcoded amplicons of the V4 region of bacterial 16S rRNA from the extracted genomic DNA. Equimolar concentrations of amplicons were then combined and quality was assessed using a 2100 BioAnalyzer (Agilent Tecnologies). High-throughput sequencing was performed at the W.M. Keck Center for Biotechnology at University of Illinois using a MiSeq with Version 3 chemistry (Illumina, San Diego CA). Sequence analysis was completed using QIIME 1.8.0 and data were evaluated using previously published methods [6, 15]. Sequences were clustered into operational taxonomic units (OTUs) using closed reference OTU against the Greengenes 13-8 reference database with a 97% sequence similarity threshold. Following an additional quality filtering , samples were rarefied to a depth of 4,088 for subsequent diversity and taxonomic analysis.
Fecal and search behavioral scores were analyzed as a chi-square using the PROC FREQ procedure. Blood work was run as an independent t-test using PROC TTEST procedure, for preand post-travel (TRAVEL dogs only, d0). Additional TTEST procedures were utilized to compare differences in blood parameters of TRAVEL dogs vs CONTROL dogs on search days (d1-d3).
Microbial data were analyzed in three different ways. First, microbial profiles were developed using principal coordinate analysis (PCoA) of UniFrac distances based on a 97% OTU composition. Data were evaluated using both weighted and unweighted models. This technique provides a visual representation of the contrast between both relative abundance and presence or absence of the phylogenetic composition for each treatment group. Second, microbial data were analyzed as a series of t-tests comparing day 0 (baseline) for both groups against each successive day. In this way, all dogs served as controls for themselves and potential confounding factors associated with breed, gender, or diet were neutralized. Lastly, microbial taxa sequence percentages were analyzed using the MIXED MODELS procedure with a Tukeypost-hoc adjustment. Statistical analysis was completed using SAS version 9.4 (SAS Institute, Inc., Cary, NC) with significance established at (P<0.05).
Perfect behavior scores (12 of 12 points, Table 1) were assessed at a rate of 44.4% and 33.3% for CONTROL and TRAVEL, respectively. No difference was observed in search behavior between the groups (Table 2). In addition, ranked search behavior scores were not significantly different (P> 0.05).
Table 2. Effects of air travel on search behavior scores for working Federal Emergency Management Agency (FEMA) canine teams.
Although all dogs for both groups started with fecal scores between 2-2.5, scores between groups were significantly different (P=0.0133) as shown in Figure 1. TRAVEL dogs had higher fecal scores indicating softer stool samples.
Blood work values were different for the TRAVEL group when pre-travel values were compared against post-travel (Table 3). Additionally, TRAVEL and CONTROL blood work demonstrated significant differences between the two groups when measured on days that included search work (Table 4). However, it is important to note that all blood work values were within the normal range and although statistically significant, it is unlikely that there is a biological significance attached.
The phylogenic composition of the treatment groups were compared using weighted UniFrac PCoA  of distances between samples based on their 97% OTU composition. Betadiversity measures revealed that TRAVEL bacterial communities and bacterial community abundances were significantly different (P=0.02) from CONTROL (Figure 2A). In addition, when unweighted UniFrac PCoA was applied to test for the presence or absence of the identified microbes, TRAVEL was different (P=0.01) from CONTROL (Figure 2B).
Overall, five phyla, 34 families, and 61 genera were identified in these adult canine fecal samples. However, > 70% of all sequences were accounted for by the Firmicutes phyla. Most abundant families for CONTROL and TRAVEL, respectively,included Clostridiaceae (27.06% and 28.32%), Lachnospiraceae (16.07% and 28.58%), Streptococcaceae (9.71% and 0.14%), and Fusobacteriaceae (9.33% and 11.29%). This is in agreement with work previously published identifying similar composition structure  in adult dogs. In addition, the Clostridia class constituted > 50% of the Firmicutes phyla and was significantly higher (P < 0.005) for dogs in the TRAVEL group (67.78%) as compared to the CONTROL group (50.99%). Also, the Bacteroidaceae family was greater (P < 0.02) in TRAVEL (6.43%) compared to CONTROL (1.09%).
When baseline (d0) sample data was compared against day 1 (24 hours post travel) and each successive search day within each group, there were several changes observed across taxonomic levels for TRAVEL. Relative abundance of Bifidobacteriaceae increased (P = 0.03) from d0 to d3, from 0% to 0.23%, respectively, in TRAVEL. At the genus level, Blautiaincreased (P = 0.05) from 8.99% at d0 to 15.02% at d3. Bifidobacterium increased (P = 0.001) from 0% at d0 to 0.17% at d3. In addition, Clostridium tended to decrease (P = 0.09) when d0 was compared to d1 (0.08% and 0.02%, respectively). No changes were observed in the CONTROL group.
Figure 2A. Weighted UniFrac principal coordinates analysis (PCoA) of TRAVEL and CONTROL working canines.
Figure 2B. Unweighted UniFrac principal coordinates analysis (PCoA) of TRAVEL and CONTROL working canines.
The effects of travel on companion animals have not been extensively studied and even fewer data are available on working canines [16, 23, 26, 29, 30]. Past research using laboratory canines revealed that during air transportation , loading and unloading seemed to be the most stressful periods . However, these dogs were transported in kennels and flown in cargo, which provides different traveling conditions than those experienced by canines flown in-cabin. In addition, FEMA canines are vigorously screened prior to acceptance into the training program and are subject to stringent certification requirements . The screening process includes testing for noise sensitivity, surface sensitivity, height sensitivity, as well as a socialization test. Therefore, it seems unlikely that a working canine would have a similar response to the airport and environmental stresses when compared to companion animals. Our results highlight the importance of further investigation of travel stress relative to the working canines’ gastrointestinal microbial composition.
Although pre-and post-travel blood parameters were affected, it is unlikely that there is any biological significance attached as all values remained within identified normal reference ranges. However, these data indicate that there is an effect on the blood metabolites for canines that travel in cabin. Therefore it is feasible that more stressful travel conditions, e.g. larger flight durations or transport in cargo, could result in biologically significant changes. More work is needed to identify how this change may impact the working canine and what, if any, steps may be taken to prevent these shifts from occurring.
Analysis of data for fecal scores revealed that the TRAVEL group was affected by the travel experience, which resulted in softer stool samples. Softer stool samples may have been due to acute stress during search performance, from the air transportation, or the significant change in the gut microbiome. Fecal scores have been widely used to identify potential fluctuations in the gastrointestinal state of dogs [11, 12]. Changes in fecal consistency are also visible indicators of environmental stress and hydration status [11, 20]. In addition, special consideration should be given to maintaining proper hydration status in working canines. It has been reported that water absorption in the colon is the primary factor associated with fecal quality in the dog, the negative impacts of which may be seen with higher fecal scores . In addition, hydration concerns are often reported during deployment and have been associated with required veterinary intervention in the past .
The effects of travel on companion animals have not been well documented. Some research in beagles  demonstrated that transportation is stressful, particularly during loading and unloading. However, these dogs were described as “laboratory animals” and were transported in kennels and flown in cargo. Additionally, it seems unlikely that the socialization available for laboratory animals would compare with the socialization experience of a pet or a service animal. Other work has shown a significant (and permanent) effect when animals were exposed to simulated air transport involving hyperthermia . Again, this simulation was typical of dogs flown in cargo, not in cabin. Kennel size was also investigated using greyhounds and showed little effect on the stress response .
The stress of travel has been well documented in human travelers. Potential stressors of travel on the ground may include commotion in the airport, baggage transport, long walks through terminals, timely arrival at the gate and potential flight delays. Inflight stressors were identified as turbulence, noise, barometric pressure changes, temperature and humidity fluctuations and general fatigue. Other factors investigated included direction of travel, altitude, and meal and sleep as well as circadian rhythm disruption. Although it is not known how canines perceive ground stressors, it seems logical that the environmental changes identified would likely be experienced similarly in the dog as in the human. This is another area that warrants further investigation.
In addition to work in humans, some work in horses has shown a sensitivity to travel that was linked to a disruption in the circadian rhythm and an interruption in the typical photoperiod. Although the canines utilized in this study traveled during the day with no perceived disruption to their circadian rhythm, it is unclear whether or not the time of travel would have an impact on the working canine. In fact, the working canines used in this study are often unable to plan their travel in advance. Frequently they are called in response to a disaster that has happened elsewhere and travel must occur immediately. More work needs to be done on time of travel in order to more fully understand this variable in the working canine.
The bacterial diversity demonstrated in our novel data set are in agreement with prior work published characterizing a healthy adult canine microbiome . Although the pilot study was limited in its scope, these data reveal an interesting perspective associated with the microbial shift that takes place in the gut of working canines following air travel. Bacterial abundance as well as communities present were different between groups as indicated by the PCoA imaging for this four day study. However, since our data were limited by a four day collection period, it is unknown whether or not there would have been additional differences in the microbial profile for subsequent days. Due to the limitations associated with the sample size in this pilot study, future studies should definitely incorporate methods to test for effects of breed, gender, and diet.
The aim of this pilot study was to examine the effect of travel stress on physical status, search performance, and composition of the gastrointestinal microbiome of the working canine. While search performance scores, and physical exam parameters were not significantly different, the blood parameters, composition of the gastrointestinal microbiome and fecal scores were impacted by travel.Further assessment of working canine fecal microbial composition and subsequent stress response may provide more information to better understand how to effectively manage working dogs in the field after air transportation. In addition, future research needs to be conducted assessing travel stress on a larger scale and utilizing a controlled diet. Lastly, future studies should be directed at identifying factors such as breed, age, dietary changes, gender and others that may result in variability of the canine gastrointestinal microbiome.
All of the research conducted in this original manuscript was funded by Southern Illinois University. No external sponsors participated in its funding. Institutional Animal Care and Use Committee (IACUC) approval was obtained prior to the initiation of this study. Special thanks to Detective Scott Mateyaschuk, Chief Alan Perry, and the canine handlers of the New York Police Department and Florida Task Force One.
- Bland A (2013) Is taking your pet on an airplane worth the risk? www.smithsonianmag.com/travel/is-taking-your-pet-on-an-airplane-worththe- risk-6241533/?no-ist
- Bergeron R, Scott SL, Emond JP, Mercier F, Cook NJ, Schaefer AL (2002) Physiology and behavior of dogs during air transport. Can J Vet Res 66(3): 211-216.
- Bokulich NA, Subramanian S, Faith JJ, Gevers D, Gordon JI, Knight R, et al., (2013) Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nat Methods10(1): 57–59.
- Bravo, J., P. Forsythe, M. Chew, E. Escaravage, H. Savignnac (2011) Ingestion of Lactobacillus strain regulates emotional behavior and central GABA receptor expression in a mouse via the vagus nerve. In PNAS 108(38):16050-16055.
- Burghardt W (2003) Behavioral considerations in the management of working dogs. Veterinary Clinics of North America: Small Animal Practice 33(2):417-446.
- Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al., (2010) QIIME allows analysis of High-throughput community sequencing data. Nat Methods 7(5): 335–336.
- Dinan TG, Cryan JF (2012) Regulation of the stress response by the gut microbiota: Implications for psychoneuroendocrinology. Psychoneuroendo 37(9): 1369-1378.
- FEMA (2008) Federal Emergency Management Agency. Canine Search Specialists Certification Process. National Urban Search and Rescue Response System (US& R).
- Gagné JW, Wakshlag JJ, Simpson KW, Dowd SE, Latchman S, Brown DA, et al., (2013) Effects of a symbiotic on fecal quality, short-chain fatty acidconcentrationsand the microbiome of healthy sled dogs. BMC Vet Res 9. doi:10.1186/1746-6148-9-246.
- Gordon LE (2012) Injuries and illness among urban search-and-rescue dogs deployed to Haiti following the January 12, 2010 earthquake. JAVMA 240(4): 396-403.
- Grellet A, Feugier A, Chastant-Maillard S, Carrez B, Boucraut-Baralon C, Casseleux G, et al., (2012) Validation of a fecal scoring scale in puppies during the weaningperiod. Prev Vet Med 106(3-4): 315-323.
- Guevara MA, Bauer LL, Garleb KA, Fahey GC, de Godoy MR (2015) Serum lipid profiles, total tract nutrient digestibility, and gastrointestinal tolerance by dogs of an α-cyclodextrin. J Anim Sci 93(5): 2201-2207.
- Hanneman GD, Higgins EA, Price GT, Funkhouser GE, Grape PM, et al., (1977) Transient and permanent effects of hyperthermia in dogs: a study of a simulated air transported environmental stress. Am J Vet Res should be journal 38(7): 955-958.
- Hernot DC, Weber MP, Biourge VC, Martin LJ, Dumon HJ, et al., (2004) Relationship between electrolyte apparent absorption and fecal quality in adult dogs differing in body size. J Nutr 134(8): 20315-20345.
- Holscher HD, Bauer LL, Gourineni V, Pelkman CL, Fahey GC, et al.,(2015) Agave inulin supplementation affects the fecal microbiota of healthy adults participating in a randomized, double-blind, placebo-controlled, crossover trial. J Nutr. 145(9): 2025-2032. doi:10.3945/jn.115.217331.
- Hooda S, Minamoto Y, Suchodolski JS, Swanson KS (2012) Current state of knowledge: The canine gastrointestinal micro biome. Anim Health Rev 13(1): 78-88.
- Horváth Z, Igyártó BZ, Magyar A, Miklósi A (2007) Three different coping styles in Police dogs exposed to short-term challenge. Hormones Behav 52(5): 621-630.
- Horváth Z, Dóka A, Miklósi A (2008) Affiliative and disciplinary behavior of human handlers during play with their dog affect cortisol concentrations in opposite directions. Hormones Behav 54(1): 107-114.
- Jacobs, DS, Kasten Jr., BL, De Mott, WR, & Wolfson, WL. (1990) Laboratory Test Handbook, 2nd ed. Baltimore: Williams & Wilkins; p.302-303.
- Kondo T, Naruse S, Hayakawa T, Shibata T (1994) Effect of exercise on gastroduodenal functions in untrained dogs. Int J Sports Med 15(4): 186-191.
- Leadon DP, E Mullins (1991) Relationship between kennel size and stress in greyhounds transported short distances by air. Vet Record 129(4): 70-73.
- Lozupone C, Lladser ME, Knights D, Stombaugh J, Knight R (2011) Uni- Frac: An effective distance metric for microbial community comparison. ISME J 5(2): 169-172.
- Middelbos, I, Fastinger, N, & Fahey, G (2007) Evaluation of fermentable oligosaccharides in diets fed to dogs in comparison to fiber standards. J Anim Sci 85(11): 3033-3044.
- Otto C, Franz MA, Kellogg B, Lewis R, Murphy L, Lauber G (2002) Field treatment of search dogs: Lessons learned from the World Trade Center disaster. Vet Emer Crit Care Soc 12(1): 33-42.
- Prescott, M, Morton, D, Anderson, D, Buckwell, A, Heath, S, Hubrecht, R, Jennings, M, Robb, D, Ruane, B, Swallow, J, & Thompson, P (2004) Refining dog husbandry and Care. Eighth report of the BAAWF/FRAME/RSPCA/UFAW joint working group on refinement. Lab Anim 1: 1-94.
- Rastall RA (2004) Bacteria in the gut: Friends and foes and how to alter the balance. J Nutr 134(8): 2022s-2026s.
- Schneider, M, & Slotta-Bachmayr, L (2009) The science of working dogs. Canine Ergonomics.CRC Press. Boca Raton, FL.
- Scholz, M, & Reinhardt, C (2007) Stress in dogs. Wenatchee, Wash.: Dogwise Pub.
- Swanson KS, Grieshop CM, Flickinger EA, Bauer LL, Healy HP, et al., (2002) Supplemental fructooligosaccharides and mannanoligosaccharidesinfluence immune function, ileal and total tract nutrient digestibilities, microbial populations and concentrations of protein catabolites in the large bowel of dogs. J Nutr 132(5): 980-989.
- Swanson K, Grieshop C, Flickinger E, Bauer L, Chow J, Wolf B, et al., (2002) Fructooligosaccharides and Lactobacillus acidophilus modify gut microbial populations, total tract nutrient digestibilities and fecal protein catabolite concentrations in healthy adult dogs. Nutr Immun 132(12):3721-3731.
- Yaguiyan-Colliard L, Grandjean D (2013) Digestive issues of working and athletic dogs. Vet Focus 23: 1-2.