ORIGINAL RESEARCH article

Front. Vet. Sci., 17 November 2022

Sec. Veterinary Infectious Diseases

Volume 9 - 2022 | https://doi.org/10.3389/fvets.2022.923792

Randomized placebo-controlled trial of feline-origin Enterococcus hirae probiotic effects on preventative health and fecal microbiota composition of fostered shelter kittens

  • 1. Department of Clinical Sciences, College of Veterinary Medicine and Comparative Medicine Institute, North Carolina State University, Raleigh, NC, United States

  • 2. Department of Environmental Services, Wake County Animal Center, Raleigh, NC, United States

  • 3. Orange County Animal Services, Chapel Hill, NC, United States

  • 4. Department of Plant and Microbial Biology, North Carolina State University, Raleigh, NC, United States

  • 5. Veterinary Hospital Pharmacy, College of Veterinary Medicine and Comparative Medicine Institute, North Carolina State University, Raleigh, NC, United States

  • 6. University of Wisconsin Veterinary Care, Madison, WI, United States

  • 7. Department of Clinical Sciences, College of Veterinary Medicine, Mississippi State University, Starkville, MS, United States

  • 8. Division of Gastroenterology and Hepatology, and UNC Microbiome Core, Department of Medicine, Center for Gastrointestinal Biology and Disease, School of Medicine, University of North Carolina, Chapel Hill, NC, United States

  • 9. IDEXX Laboratories, Inc., West Sacramento, CA, United States

  • 10. Vet Med Labor GmbH Division, IDEXX Laboratories, Inc., Kornwestheim, Germany

Abstract

Introduction:

Diarrhea is the second most common cause of mortality in shelter kittens. Studies examining prevention strategies in this population are lacking. Probiotics are of particular interest but studies in cats are largely limited to healthy adults or those with induced disease. Only one study in domestic cats describes the use of host-derived bacteria as a probiotic. We previously identified Enterococcus hirae as a dominant species colonizing the small intestinal mucosa in healthy shelter kittens. Oral administration of a probiotic formulation of kitten-origin E. hirae (strain 1002-2) mitigated the increase in intestinal permeability and fecal water loss resulting from experimental enteropathogenic E. coli infection in purpose-bred kittens. Based on these findings, we hypothesized that administration of kitten-origin E. hirae to weaned fostered shelter kittens could provide a measurable preventative health benefit.

Methods:

We conducted a randomized, placebo-controlled, blinded clinical trial to determine the impact of a freeze-dried E. hirae probiotic on body weight gain, incidence of diarrhea, carriage of potential diarrheal pathogens, and composition of the intestinal microbiota in weaned fostered shelter kittens.

Results:

One-hundred thirty kittens completed the study. Fifty-eight kittens received the probiotic and 72 received the placebo. There were no significant differences in age, weight upon initiation of the study, number of days in the study, average daily gain in body weight, or weight at completion of the study. Kittens treated with E. hirae were 3.4 times less likely to develop diarrhea compared to kittens treated with placebo (odds ratio = 0.294, 95% CI 0.109–0.792, p = 0.022). A significant impact of E. hirae was not observed on the presence or abundance of 30 different bacterial, viral, protozoal, fungal, algal, and parasitic agents in feces examined by qPCR. With exception to a decrease in Megamonas, administration of the E. hirae probiotic did not alter the predominant bacterial phyla present in feces based on 16S rRNA gene amplicon sequencing.

Discussion:

Decreased incidence of diarrhea associated with preventative administration of E. hirae to foster kittens supports a rationale for use of E. hirae for disease prevention in this young population at high risk for intestinal disease though additional studies are warranted.

Introduction

It is estimated that 180 million kittens are born each year in the United States and hundreds of thousands of orphaned and abandoned kittens are fostered by U.S. animal shelters (–). Diarrheal disease is the second most common illness to afflict shelter kittens and many will die from diarrhea before they reach an age of 8 weeks (, ). The causes of diarrhea are likely multifactorial including an immature immune system, weaning stress, independent eating, change in food type, waning maternal antibodies, and increased exposure to diarrheal pathogens (, , , ). Efforts aimed at diagnosis and treatment of infectious causes of diarrhea in this population is challenging due to the large number and diversity of potential pathogens (viral, bacterial, parasitic and protozoal), frequent carriage of these pathogens by kittens without diarrhea, many agents for which there is no specific treatment, and limited resources. Consequently, orphan kittens that develop diarrhea are often treated empirically for intestinal bacterial, protozoal, or parasitic infections and are commonly provided supportive care in the form of antibiotics, probiotics, vitamin supplements, subcutaneous fluids, or assisted feeding ().

Among the many knowledge gaps including causes and treatment outcomes of foster kitten diarrhea is a lack of studies examining prevention strategies for diarrhea in this population. In particular, probiotics () are of particular interest due to the multiple proposed beneficial effects such as promoting resistance to colonization by pathogens, stabilization of the intestinal microbiota, immunomodulatory effects on gut function, and metabolic impact on nutrient digestion. To date however, research examining the impact of probiotics in cats are mostly limited to investigations of healthy, adult-aged, purpose-bred cats, or cats with experimentally-induced disease (e.g., FIV and FHV-1) or administered antibiotics (–). Only 2 studies have examined the effect of probiotics on cats with diarrhea. In a study of older stray and feral cats housed short-term in a shelter setting, administration of a probiotic containing Enterococcus faecium SF68 did not protect against diarrhea but decreased the number of cats having diarrhea of ≥2 days duration (). In a second study of adult cats with naturally-occurring chronic diarrhea, treatment with a synbiotic containing Bifidobacterium, Enterococcus, and Lactobacillus was reported to significantly improve fecal score ().

In particular, the enterococci are components of the normal microbiota of the gastrointestinal tract and are ubiquitous in ecosystems ranging from soil, water, plants and food products (). We have previously identified that a specific species of enterococcus, Enterococcus hirae, is the principal species inhabiting the mucosa-associated microbiota of the small intestine in healthy fostered shelter kittens () and is frequently observed to extensively colonize the surface of the intestinal epithelium (). In kittens succumbing to severe illness, E. hirae are displaced by opportunistic and virulent members of enterococci such as Enterococcus faecalis and such kittens are more likely to be diagnosed with enteroadherent enteropathogenic Escherichia coli infection (EPEC) (). Using a strain of E. hirae, cultured from the mucosa-associated microbiota of a healthy kitten and manufactured into a freeze-dried probiotic, we demonstrated that oral administration of E. hirae could mitigate the increase in intestinal permeability and fecal water loss resulting from experimental EPEC infection in purpose-bred kittens ().

Based on these findings, we hypothesized that administration of feline-origin E. hirae to weaned fostered shelter kittens could provide a measurable preventative health benefit. To test this hypothesis we conducted a randomized, placebo-controlled, blinded clinical trial to determine the impact of a freeze-dried E. hirae probiotic on parameters of body weight gain, incidence of diarrhea requiring examination by a veterinarian, carriage of potential diarrheal pathogens, and composition of the intestinal microbiota in weaned fostered shelter kittens.

Materials and methods

Enterococcus hirae probiotic

Formulation

The administered freeze-dried probiotic was manufactured from a pure, well-characterized isolate of Enterococcus hirae. The strain (1002-2) was isolated by swab of the ileum mucosa performed during autopsy of a healthy kitten in which histological evidence of E. hirae attachment to the intestinal epithelium was concurrently documented (). The isolate lacked all virulence characteristics, the results of which were previously reported (). The selected strain was cultured for use as a probiotic in Biostat B-plus reactors continuously fed with MB medium () containing 20 g/L glucose. Chemostat maintained a dilution rate (D = 0.17 h−1), temperature (37°C), and pH (5.5). The viable biomass was collected from the reactor overflow, concentrated, and suspended in bacteriological peptone prior to lyophilization. Glucose consumption, accumulation of lactate and absence of alternative fermentation products were monitored with high-performance liquid chromatography (Shimadzu Corporation, Kyoto, Japan) performed under isocratic conditions at 65°C, a mobile phase of 5 mM sulfuric acid (H2SO4) at a 0.4 ml/min flow rate using an Alltech IOA-1000 organic acids column (300 mm x 7.8 mm, Alltech, IL, USA), and coupled to a refractive-index detector. Numbers of viable cells were confirmed by plating on M17 agar () and incubating at 37°C for 72 h. For each sample, three plates were prepared from each dilution and recorded as the average of three independent counts. Finalized freeze-dried peptone containing probiotic and vehicle (sterile ultrapure bacteriological peptone, Thermo Scientific) were stored at −80°C.

Compounding

Capsules containing either the probiotic or placebo (bacteriological peptone) were compounded by the NC State Veterinary Hospital Pharmacy into blue or yellow opaque capsules whose contents remained blinded to study investigators.

Pre-compounding E. hirae stock containing 2.24 × 107 colony forming units (CFU) per milligram was stored at −80°C and transferred to dry ice when in use. A desired dose of 1 × 108 CFU was achieved by compounding 4.5 mg of E. hirae stock with 110.5 mg of bacteriological peptone (ThermoFisher Scientific, Waltham, MA) to achieve an average pack weight of 115 mg ± 5 mg. Placebo consisted of 115 mg of bacteriological peptone only. Compounding was performed in a 4°C walk-in environmental chamber. Batch sizes of 300 capsules were compounded using a pre-chilled mortar and pestle, ProFiller 3600 capsule machine, and Letco #4 size gelatin capsules. Capsules were transferred to 30–count vials labeled with capsule color and date of manufacture and stored at −20°C.

Viability

Sustained numbers of E. hirae in the compounded probiotic was confirmed by culture of the contents of a subset of capsules immediately after formulation, after 5 months of storage at −20°C on location at the shelter, and following return of leftover (unused) capsules to the shelter after storage in the home freezer of a foster care provider. Briefly, probiotic and placebo capsule contents were dissolved separately in 1 ml of sterile phosphate-buffered saline (PBS) to generate a working stock. Serial dilutions were performed to generate concentrations ranging from 1 × 107 to 1 × 102 CFUs per 1 ml and the equivalent for the placebo. For each of the 104, 103, and 102 dilutions and sterile PBS, 100 μl was spread onto three separate plates of Difco m-Enterococcus Agar (Becton Dickinson, Franklin Lakes, NJ). Plates were incubated at 37°C for 24 h. Following incubation, colonies were counted on the plates of the dilution(s) that resulted in 30 and 300 CFUs. The average colony count was calculated, allowing for the current CFU concentration of the working stock to be extrapolated. Viability of the probiotic capsules were then determined by dividing the extrapolated concentration by the initial presumed concentration of 1 × 108 CFUs/ml.

Clinical trial design

A prospective, randomized, placebo-controlled, blinded clinical study was conducted at a large municipal animal shelter [Wake County Animal Center (WCAC)] over the time interval from May 1 to October 1 in the year 2018 and 2019.

Eligibility criteria included any weaned and apparently healthy kitten aged <12 weeks and admitted to theshelter with a body weight < 0.9 kg (2 lbs) and therefore requiring temporary care prior to fitness for surgical sterilization and adoption. Kittens destined for fostering by a WCAC-trained care provider were eligible for participation in the study. Kittens that were alternatively transferred to the care of rescue partner organizations were not eligible to participate.

Preventative care

All kittens were housed for an average of 2–3 days in a dedicated room within the shelter prior to transfer to the home of a volunteer foster care provider. During this time, each kitten was assigned an individual identification/medical record number for tracking using a web-based shelter software program (ShelterBuddy, Wacol Brisbane QLD AUSTRALIA). Kittens were administered a single oral dose of ponazuril paste (Marquis®, Boehringer Ingelheim; 30mg/kg) and pyrantel pamoate (Columbia Laboratories; 11 mg/kg) and topical administration of selamectin (Revolution®, Zoetis) for treatment of endo- and ectoparasites. Kittens that were ≥4 weeks of age and ≥0.45 kg were vaccinated by subcutaneous injection with feline viral rhinotracheitis, calicivirus, and panleukopenia antigen (NobivacR:Feline 1-HCP, Intervet INC Merck Animal Health division US Omaha NE 68103) and revaccinated every 2–3 weeks until 16–20 weeks ().

Randomization

In advance of the study, each day of each month from May 1 to October 1 was randomly designated as blue or yellow using an online program (Random.org). The designation of blue or yellow indicated the treatment group to which any kittens entering foster care on that day would be assigned and was unconcealed. Participation in the study was voluntary and each foster volunteer signed an informed consent. To prevent cross-contamination of treatments, participating foster-care providers were only allowed to simultaneously foster kittens designated to the same treatment group at any given time.

Blinding

The placebo and probiotic were formulated into opaque size number 4 gelatin capsules colored blue and yellow to correspond with the treatment group allocation color. Study investigators, shelter staff, and foster care providers were blinded to identity of the capsule contents. Only the pharmacists (SW, EW, and GD) responsible for compounding the capsules and technician (SHS) responsible for documenting bacterial viability of capsule contents were aware of the capsule contents. All other investigators remained blinded to group allocation until all data were analyzed and results of the study were written in manuscript form.

Treatment and monitoring

Participating foster care providers were provided with capsules containing probiotic or placebo based on prior random allocation, two zip-lock bags for collection of a pre- and end-of-study fecal sample with labels for providing the collection date and kitten shelter identification number, observation worksheet for recording pre- and end-of-study body weight, and instructions for treatment and sample collections. Foster care providers were instructed to store the dispensed capsules in their home freezer throughout the study.

During the study, the type of diet provided to kittens was chosen at the discretion of the foster care provider. Supplements such as vitamins were permitted but not administration of non-study probiotics. Prior to beginning administration of the probiotic or placebo, foster care providers were instructed to obtain an initial body weight for which a scale (AccuWeight AW-KS001BB) was provided, and to collect a voided fecal sample with instructions to store the sample in their home freezer. After recording of an initial body weight and collection of a fecal sample, foster care providers were instructed to administer a single capsule of probiotic/placebo daily by opening the contents of the capsule onto a small portion of canned kitten food and monitoring consumption by the kitten until complete. Each foster-care provider was provided with an observation sheet for documenting daily probiotic administration and encouraged, but not required, to take periodic body weight measurements over the course of the study. Upon reaching a weight of ≥2 pounds, foster-care providers were instructed to obtain a final body weight and fecal sample while the kitten was still receiving the probiotic and store as per the initial sample. Upon completion of the study, fecal samples were transferred by the foster care provider back to Wake County Animal Center and stored at −20°C prior to transport to the laboratory on dry ice where they were subsequently stored at −80°C until analysis.

Incidence of diarrhea

Kittens were diagnosed as having significant diarrhea if they met criteria for unscheduled examination by the WCAC veterinary team. These criteria included diarrhea lasting more than 24 h and accompanied by straining to defecate, discomfort during defecation, or decreased appetite. Foster care providers were not trained to provide fecal scoring and fecal scores were not requested or recorded as part of the study.

Criteria for inclusion of gathered data

Minimum criteria for inclusion in data analysis were (1) uninterrupted daily administration of the placebo/probiotic and (2) recording of a pre- and end-of-study body weight. An end-of-study body weight was not required for kittens that died or were euthanized while on the study.

Molecular assays

Fecal PCR testing

Paired fecal samples, collected from kittens prior to beginning and immediately prior to completion of the study, were batched and shipped on dry ice to a collaborating laboratory (IDEXX Diagnostics Laboratory, Sacramento CA) where they were stored at −80°C prior to assay.

Fecal samples were suspended in 3 ml guanidinium thiocyanate based lysis solution (), vortexed to facilitate organism detachment and rapid protein denaturation and incubated at room temperature for 15 min. The lysate was then used to extract total nucleic acid on a MagMax 96 Flex (Life Technologies, Valencia, CA) with magnetic beads (Roche, Indianapolis, IN) using manufacturer's guidelines. Total nucleic acid was eluted in 200 μl of PCR-grade nuclease-free water and 5 μl amplified in subsequent single plex real-time PCR reactions.

Real-time PCR was performed with proprietary forward and reverse primers and hydrolysis probes. Target genes for enteropathogen detection were as follows: Feline enteric coronavirus 7b gene (DQ010921.1), Feline panleukopenia virus VP2 gene (EU252145), Rotavirus VP7 gene (EU708950), Salmonella enterica invasion A gene (EU348366), Campylobacter (coli and jejuni) IpxA gene [AY531496 (C. coli) and AL111168 (C. jejuni)], Clostridium perfringens alpha toxin gene (L43545), Clostridium perfringens enterotoxin gene (X81849), Shigella spp. RFC gene (AE005674), E. coli virulence factor genes (F4 (K88), F5 (K99), F6 (987P) and F18ab fimbrial adhesin genes; STX2e gene; Intimin adhesion gene; LT1 gene, ST1a and ST1b genes) (EU570252; M35282; M35257; GQ325633; GU945540; ECU38618; S60731.1; M34916), Yersinia spp. invasion locus protein (ail) gene (AJ605740), Cryptosporidium spp. ssrRNA gene (A093489), Giardia small-subunit rRNA gene (DQ836339), Toxoplasma gondii internal transcribed spacer-1 gene (L49390, Tritrichomonas foetus 5.8S rRNA gene (AF339736), Blastocystis spp. SSU rRNA gene (AB023499; KP890050; AB070992; AB071000; AB070999; AB070990; AB070991; AB107971; AF408426), Toxocaris leoni ITS-2 gene (Y09490), Toxocara (malayasiensis and cati) ITS-1 gene, Internal transcribed spacer-2 (ITS-2) gene (AB110033; AM231609), Ancylostoma spp. ITS-1 gene (DQ438074; DQ438063; DQ780009), Trichuris spp. ITS-2 gene (AM234616), Uncinaria spp. ITS-1 gene (AF217890), Parascaris spp. ITS-1 gene (AJ007459), Echinococcus spp. rRNA between the Cox 1 and Cox 2 genes (NC_000928), Pythium spp. ITS-1 gene (EF016907.1), and Prototheca spp. ITS-1 gene (AJ245645).

Real-time PCR was run with six quality controls as previously described () including (1) PCR positive controls [synthetic DNA (Integrated DNA Technologies IDT, Coralville, IA), run quantitatively], (2) PCR negative controls (RNase-free PCR-grade water, Fisher Scientific, Waltham MA), (3) negative extraction controls (lysis solution only), (4) DNA pre-analytical quality control targeting mammalian ssr rRNA (18S rRNA) gene complex, (5) environmental contamination monitoring control (swab-based laboratory monitoring), and (6) spike-in internal positive control. These controls assessed the functionality of the PCR test protocols for the (1), functional assessment of the real-time PCR test performance (2) absence of contamination (both PCR product carry-over and sample cross-contamination), (3) absence of detectable cross-contamination during the extraction process, (4) quality and integrity of the DNA as a measure of sample validity (by quantitatively assessing 18S gene load), reverse transcription protocol (5) absence of aerosol-based contamination within the PCR laboratory space, and (6) absence of PCR inhibitory substances as a carryover from the sample matrix (internal sample control with spike-in DNA).

Analysis was performed on a Roche LightCycler 480 (Roche Applied Science, Indianapolis (IN) on which amplification data were analyzed using the 2nd derivative maximum method to generate crossing points (Cp values). Each Cp value was quantitatively validated by prior calibration to a standard curve of DNA for each target assay.

16S rRNA gene amplicon sequencing

Fecal samples that were collected immediately prior to completion of the study were selected for microbiota analysis. To isolate the impact of the probiotic intervention and minimize bias in the data, sibling kittens had only a single/representative fecal sample chosen for analysis.

DNA isolation from fecal samples was performed as previously described (). Briefly, 0.3 ml of Qiagen ATL buffer (Valencia, CA) supplemented with 60 mg/ml lysozyme (Thermo Fisher Scientific, Grand Island, NY) was added to the samples. The suspensions were transferred to 2 ml tubes containing 106/500 μm glass beads (Sigma, St. Louis, MO) and incubated for 1 h at 37°C with occasional agitation. The suspensions were agitated for 40 min on a Digital Vortex Mixer. Subsequently, the suspensions were supplemented with 20 μL of Qiagen proteinase K and 0.3 ml of Qiagen AL buffer and incubated at 55°C overnight. After brief centrifugation, supernatants were aspirated and transferred to a new tube containing 0.3 ml of ethanol. DNA was purified using a standard on-column purification method with Qiagen buffers AW1 and AW2 as washing agents and eluted in 27.5 μl of DNase free water and quantified using QuantIT® PicoGreen®.

Extracted DNA was quantified via PicoGreen analysis and used for bacterial 16S rRNA gene amplicon sequencing. Total DNA (12.5 ng) was amplified using universal primers targeting the V4 region of the bacterial 16S rRNA gene (). The sequences of the primers were: 515F−5′ TCGTCGGCAGCGTCAGATGTGTATAAGAGACAGGTGCCAGCMGCCGCGGTAA 3′ and 806R−5′GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAGG GACTACHVGGGTWTCTAAT 3′. Overhang adapters were appended to the 5′ end of each primer sequence for compatibility with the Illumina sequencing platform. Master mixes contained 12.5 ng of total DNA, 0.2 μM of each primer and 2× KAPA HiFi HotStart ReadyMix (KAPA Biosystems, Wilmington, MA). The thermal profile for the amplification of each sample had an initial denaturing step at 95°C for 3 min, followed by 25 cycles of denaturing of 95°C for 30 s, annealing at 55°C for 30 s for 16S rRNA and a 30 s extension at 72°C, a 5 min extension at 72°C and a final hold at 4°C. Each 16S amplicon was purified using the AMPure XP reagent (Beckman Coulter, Indianapolis, IN). In the next step each sample was amplified using a limited cycle PCR program, adding Illumina sequencing adapters and dual index barcodes [index 1(i7) and index 2(i5)] (Illumina, San Diego, CA) to the amplicon target. The thermal profile for the amplification of each sample had an initial denaturing step at 95°C for 3 min, followed by a denaturing cycle of 95°C for 30 s, annealing at 55°C for 30 s and a 30 second extension at 72°C (8 cycles), a 5 min extension at 72°C and a final hold at 4°C. The final libraries were again purified using the AMPure XP reagent (Beckman Coulter), quantified with Quant-iT™ PicoGreen® dsDNA Reagent (Molecular Probes, Thermo Fisher Scientific, Waltham, MA), and normalized prior to equimolar pooling. The DNA library pool was then denatured with NaOH, diluted with hybridization buffer and heat denatured before loading on the MiSeq reagent cartridge MiSeq instrument (Illumina). Automated cluster generation and paired–end sequencing with dual reads were performed according to the manufacturer's instructions.

Data analysis

Kitten populations were described in terms of age, sex, weight on intake to the shelter, number of kittens in litter, time spent on study, average daily gain in body weight, weight on exit from study, episodes of diarrhea requiring examination by shelter veterinary staff, and outcome (alive–vs.–died or euthanized). Categorical data (e.g., proportions) were tested for significant differences between treatment groups using a Chi-square with calculation of odds ratios and 95% confidence intervals. Continuous data were tested for normal distribution and equal variance followed by parametric (one-way ANOVA) or non-parametric (Kruskal-Wallis one-way ANOVA on ranks) analysis as appropriate.

The proportion of kittens testing PCR positive for detection of selected infectious agents in feces were compared for significant differences between pre- and post-study samples within each treatment group using a Chi Square test or Fisher Exact test. Infectious agents for which ≤5 kittens tested positive at both time points were not subjected to statistical analysis. The median burden of selected infectious agent RNA/DNA detected within feces of test-positive kittens, as represented by the PCR Cp values, was compared for significant differences between pre- and post-study samples within each treatment group using a Kruskal-Wallis one-way ANOVA on ranks. Testing was performed using Systat software (SigmaPlot 12).

Sequencing output from the Illumina MiSeq platform were converted to fastq format and demultiplexed using Illumina Bcl2Fastq 2.18.0.12. The resulting paired-end reads were processed using QIIME 2 () 2018.11. Index and linker primer sequences were trimmed using the QIIME 2 invocation of cutadapt. The resulting paired-end reads were processed with DADA2 through QIIME 2 including merging paired ends, quality filtering, error correction, and chimera detection. Amplicon sequencing units from DADA2 were assigned taxonomic identifiers with respect to Green Genes release 13_08.

Alpha diversity estimates were calculated within QIIME 2 using Evenness (Shannon) index and observed species number metrics at a rarefaction depth of 5,000 reads. Pairwise significance was tested using a Kruskal-Wallis ANOVA with Benjamini-Hochberg corrected q-values calculated as implemented in QIIME 2. Beta diversity estimates were calculated within QIIME 2 using weighted and unweighted Unifrac distances as well as Bray-Curtis dissimilarity between samples at a subsampling depth of 5,000 reads. Results were summarized, visualized through principal coordinate analysis in Emperor, and significance was estimated by PERMANOVA with Benjamini-Hochberg corrected q-values calculated as implemented in QIIME 2. Relative abundance data were normalized by sample library size and taxa were removed if they were present in <10% of all samples or if lower than 0.01% average abundance (). Differences in abundance were tested for significance using a Kruskal-Wallis One-Way ANOVA on Ranks. When significant, pairwise comparisons between groups were performed using a Dunn's test. Generated p-values were corrected for multiple testing using the Benjamini-Hochberg procedure () at a false discovery rate of 0.25. For selected taxa, odds ratio and 95% confidence intervals for association with treatment group were calculated using a Chi-square Test with Yates Continuity Correction. Statistical examination of abundance data was conducted using Systat Software (SigmaPlot 12.0, San Jose, CA). Results were represented graphically using GraphPad Software (Prism version 7.03, San Diego, CA).

Results

Population description

Over the time interval from May 1 to October 1 in the years 2018 and 2019, a total of 1,471 kittens were assessed for eligibility to participate in the study. Among these kittens, 220 (15%) were randomly allocated to treatment with probiotic or placebo. One-hundred thirty kittens completed the study. Fifty-eight kittens from 40 unrelated litters received the probiotic (blue group) and 72 kittens from 37 unrelated litters received the placebo (yellow group) (Figure 1). A total of 28 different foster care providers participated over the 2-year course of the study, each fostering a median of 3 kittens (range, 1–18 kittens).

Figure 1

Clinical outcome

Viability of the probiotic (target dose, 1 × 108 CFU) was not diminished when quantified after storage at −20°C at the shelter (2.2–8.9 × 108 CFU), nor followed by storage in the home freezer of a foster care provider (0.94 × 108 CFU). All kittens readily consumed the probiotic or placebo when offered. When comparing the probiotic and placebo groups, there were no significant differences in starting age, weight upon initiation of the study, number of days in the study, average daily gain in body weight, or weight at completion of the study (Table 1 and Figure 2).

Table 1

Parameters evaluatedProbiotic (n = 58) male = 32, female = 26 undetermined = 0Placebo (n = 72) male = 36, female = 35, undetermined = 1Kruskal-Wallis One-Way ANOVA on Ranks p-value
MedianIQRMedianIQR
Estimated age (weeks)65–865–70.429
Number of days in study1714–262214–310.123
Weight at start of study (grams)602455–704501414–6490.075
Weight at end of study (grams)933†883–998964908–1,0210.205
Average daily gain in body weight (grams/day)a17†15–211714–230.904
No.%No.%Chi-square test p- value
Kittens examined for diarrheal illness6b1020c280.02
Kittens euthanized for diarrheal illness0000–
Kittens euthanized for non-diarrheal illness35110.17
Kittens that survived55957199

Population description of 130 kittens treated with probiotic or placebo.

IQR, interquartile range.

a

Calculated as = (end weight–start weight)/number of days between weight measurements.

b

From 6 different litters.

c

From 8 different litters.

†

Excluding 2 kittens that were euthanized.

Figure 2

Kittens treated with the probiotic were 3.4 times less likely to develop diarrhea compared to kittens treated with the placebo (odds ratio = 0.294, 95% CI 0.109–0.792, p = 0.022). No kittens were euthanized due to diarrheal illness and overall survival was not significantly different between the 2 groups. Reasons for euthanasia of kittens in the study included failure to thrive (2 probiotic kittens), panleukopenia (1 probiotic kitten), and major medical (1 placebo kitten) (Table 1).

Prevalence and abundance of selected infectious agents

Quantitative PCR was performed on paired fecal samples, collected prior to beginning and immediately prior to completion of the study, from 37 kittens receiving the probiotic and 55 kittens receiving placebo. Reasons for exclusion of kittens were lack of either a pre or post-fecal sample (13 probiotic kittens, 15 placebo kittens), insufficient feces for both microbiome and PCR analysis (5 probiotic kittens, 1 placebo kitten), or euthanasia prior to final fecal sample collection (3 probiotic kittens, 1 placebo kitten).

Fecal samples were tested for the presence of 30 different bacterial, viral, protozoal, fungal, algal, and parasitic agents, among which 17 were detected in ≥1 kitten (Tables 2 and 3). Prior to initiation of the study, 92% (34/37) of kittens receiving probiotic and 87% (48/55) of kittens receiving placebo had a median of 2 infectious agents (range, 0–4) identified. Upon completion of the study, 86% (32/37) of kittens receiving probiotic and 78% (43/55) of kittens receiving placebo had a median of 2 (range, 0–4) and 1 (range, 0–4) respectively, infectious agents identified. The number of identified agents did not differ significantly from pre-study values in either group. Examination of individual kittens demonstrated that 22% (8/37) of kittens receiving the probiotic had a median gain or loss in number of infectious agents of 0 (range, −3 to +3) and 20% (11/55) of kittens receiving the placebo had a median gain or loss in number of infectious agents of 0 (range, −4 to +3) which did not differ significantly between the groups.

Table 2

Selected agentsPrevalence/presence of infectious agent(s)
Probiotic (n = 37)Placebo (n = 55)
Pre-studyPost-studyP-valuePre-studyPost-studyP-value
# kittens passing QC# kittens positive% kittens positive# kittens passing QC# kittens positive% kittens positive# kittens passing QC# kittens positive% kittens positive# kittens passing QC# kittens positive% kittens positive
Coronavirus371643.2371437.80.813552036.4551730.90.686
E.coli Intimin (eae)341029.433618.20.42953815.148714.60.835
Feline Panleukopenia34926.5332†6.10.054531018.94400.00.002††
Salmonella spp.3400.03326.1–5247.744613.60.505
Cryptosporidium spp.3425.93300.0–5335.74436.8–
Giardia spp.3400.03300.0–5200.04536.7–
Toxoplasma gondii3400.03313.0–5200.04400.0–
Campylobacter jejuni3400.03326.1–5323.84524.4–
Campylobacter campi3400.03326.1–5200.04524.4–
Toxocaris leoni3412.93300.0–5200.04400.0–
Toxocara malayasiensis3412.93300.0–5200.04412.3–
Toxocara cati3412.93313.0–5200.04412.3–
Ancylostoma spp.3400.03300.0–5211.94400.0–
Pythium spp.3525.73339.1–5223.84536.7–
Shigella spp.3400.03300.0–5247.74400.0–
Trichuris sp.3400.03300.0–5200.04512.2–
Clostridium perfringens
Alpha toxin
3528803528800.765554480483062.50.080
> 300 Thous/g feces35617.13525.70.25955712.748510.40.955
Clostridium perfringens
Enterotoxin
3438.83313.0–5235.84400–
> 300 Thous/g feces3412.93300–52004400–

Prevalence of infectious agents identified by RT-PCR or PCR in RNA/DNA extracted from voided feces collected from fostered kittens prior to (pre-study) and upon completion (post-study) of ongoing daily treatment with a probiotic or placebo.

QC, Quality Control. The proportion of kittens testing PCR positive for detection of selected infectious agents in feces were compared for significant differences between pre- and post-study samples within each treatment group using a Chi Square or Fischer Exact test.

– Infectious agents for which ≤ 5 kittens tested positive at both time points were not subjected to statistical analysis.

†

Both kittens negative for panleukopenia DNA at time of pre-study PCR testing.

††

Remains significant after correction for multiple testing using a Benjamini-Hochberg false discovery rate of 0.15.

Table 3

Selected agentsBurden of infectious agent(s)
Probiotic (n = 37)Placebo (n = 55)
Pre-studyPost-studyKruskal-Wallis ANOVAPre-studyPost-studyKruskal-Wallis ANOVA
Median
Cp
Min
Cp
Max
Cp
Median
Cp
Min
Cp
Max
Cp
Median
Cp
Min
Cp
Max
Cp
Median CpMin
Cp
Max
Cp
Coronavirus3823392923390.1513821392824390.035†
E.coli Intimin (eae)3429373526361.003220363526380.769
Feline Panleukopenia2017252015261.00221624––––
Salmonella spp.–––403940–3937403836380.352
Cryptosporidium spp.373737––––383638383638–
Giardia spp.––––––––––353440–
Toxoplasma gondii–––373737––––––––
Campylobacter jejuni–––383838–302832383739–
Campylobacter campi–––353437––––363636–
Toxocaris leoni373737–––––––––––
Toxocara malayasiensis363636–––––––343434–
Toxocara cati383838353535––––363636–
Ancylostoma spp.–––––––363636––––
Pythium spp.393840383739–393940393839–
Shigella spp.–––––––343037––––
Trichuris sp.––––––––––393939–
Clostridium perfringens Alpha toxin3122353127340.1643022353223350.186
gene copies/ g feces51,3103,33523,235,02635,0444,160465,7210.16489,7162,60120,236,83426,2223,2008,652,4450.186
Clostridium perfringens Enterotoxin292531292929–333233000–
gene copies/g feces109,17436,4012,547,614157,442157,442157,442–12,82710,21217,624000–

Abundance of infectious agents identified by RT-PCR or PCR in RNA/DNA extracted from voided feces collected from fostered kittens prior to (pre-study) and upon completion (post-study) of ongoing daily treatment with a probiotic or placebo.

Cp, PCR cycle threshold. – Infectious agent not detected.

†

Not significant after correction for multiple testing using a Benjamini-Hochberg false discovery rate of 0.15.

The most common infectious agents having RNA or DNA demonstrated in feces were feline enteric coronavirus, attaching and effacing E. coli (presence of eae), panleukopenia, and Clostridium perfringens (alpha-toxin). There were no significant differences in prevalence (Table 2) or abundance (Table 3) of infectious agents between the 2 groups of kittens for fecal samples collected prior to initiation of the study. Upon completion of the study there was a significant decrease in prevalence of detection of panleukopenia DNA in kittens treated with the placebo (Table 2). Decrease in abundance of coronavirus RNA in kittens treated with placebo was not significant after correction for multiple comparisons (Table 3).

Infectious agent PCR testing was performed on feces collected from 16/26 kittens that developed diarrhea over the course of the study. Identified infectious agents having a potential to contribute to diarrhea included panleukopenia (4 kittens from 2 litters), feline enteric coronavirus (5 kittens from 3 litters), Campylobacter spp. (3 kittens from 2 litters), attaching and effacing E. coli (2 kittens), and alpha-toxin positive Clostridium perfringens (alpha toxin gene copies >30,000/gram feces) (2 kittens).

Microbial community composition

Analysis of the fecal microbiota by 16S rRNA gene amplicon sequencing was performed for 32/58 (55%) unrelated kittens that received the probiotic (32/40 (80%) of litters represented) and 35/72 (49%) unrelated kittens that received the placebo (35/37 (94%) of litters represented). Included in the analysis were 6 kittens from 6 individual litters that received probiotic and were examined for diarrheal illness during the course of the study. Eight litters in the probiotic group did not have a kitten represented in the analysis due to lack of an end-of-study fecal sample (5 litters) or euthanasia (3 litters). Also included in the analysis were 7 kittens from 7 individual litters that received placebo and were examined for diarrheal illness during the course of the study. Two litters in the placebo group did not have a kitten represented in the analysis due to lack of an end-of-study fecal sample (one of which included a kitten that developed diarrhea during the course of the study).

The sequence analysis of all 67 kittens yielded 5,675,488 quality sequences (mean ± SD = 84,709 ± 34,177). The major phyla represented in the fecal microbiota of kittens were Firmicutes, Actinobacteria, Bacteroidetes, Proteobacteria, and Fusobacteria (Figure 3). There were no significant differences in alpha diversity measures of evenness (Shannon), diversity (Faith), or observed number of taxa between kittens receiving the probiotic compared to placebo. Both weighted and unweighted measures of beta diversity did not identify significant differences in phylogenetic composition of the microbiota between kittens receiving the probiotic vs. placebo (Figure 4).

Figure 3

Figure 4

Kittens receiving the probiotic had a significantly higher median relative abundance of the genus Enterococcus and lower abundance and prevalence of Megamonas (Figure 5). Abundance of specific taxa observed in the feces of kittens receiving probiotic vs. placebo are shown in Table 4. Kittens receiving the probiotic were 10.9 times less likely to have detectable Megamonas (odds ratio = 0.157; 95%CI = 0.054–0.453, p < 0.001, Chi-square test) present in the microbial community compared to kittens receiving the placebo.

Figure 5

Table 4

PhylumFamilyGenusRelative % abundancePrevalence among kittens
Probiotic (n=32 kittens)Placebo (n=35 kittens)Kruskal-Wallis ANOVA p-valueProbiotic (n=32 kittens)Placebo (n=35 kittens)Chi-square p-value
MedianRangeMedianRange#%#%
Unassigned0.00220–0.440.00–0.2300.1521856.31234.30.119
Bacteria0.0120–0.150.0160–0.520.2342371.92777.10.831
Actinobacteria12.90.043–76.010.91.08–86.30.8513210035100–
o__Actinomycetales0.00–0.0280.00–0.0520.75039.38411.41.00
ActinomycetaceaeActinomyces0.00–0.0870.00–0.0410.0391546.97200.038
CorynebacteriaceaeCorynebacterium0.00–0.0150.00–0.200.229412.5822.90.432
BifidobacteriaceaeBifidobacterium0.340–66.10.400–52.40.7022990.63497.10.342
Coriobacteriaceae0.400–63.10.160–5.590.1292887.53291.40.701
Adlercreutzia0.0220–1.260.0570–2.880.3811959.42365.70.777
Collinsella6.70–33.08.460.75–74.90.163196.9351000.478
Slackia0.170–1.310.270–1.620.4292990.63394.30.664
Bacteroidetes0.430–24.31.120.012–33.10.1323196.9351000.478
BacteroidaceaeBacteroides0.0310–23.30.160–18.50.0972681.33188.60.501
PorphyromonadaceaeParabacteroides0.00–1.00.00–4.080.021618.81645.70.037
PrevotellaceaePrevotella0.190–10.30.210–13.70.3722371.92674.30.957
[Prevotella]0.00–2.10.00–1.530.557412.5617.10.736
Cyanobacteria0.00–0.0870.00–2.060.173412.5925.70.291
o__Stramenopiles0.00–0.0870.00–0.00.865412.5514.31.00
Firmicutes81.423.2–99.878.713.5–96.60.533210035100–
o__Bacillales0.000–0.00.00–0.140.10513.13514.30.200
Bacillaceae0.0150–0.260.00860–0.760.7261959.421600.844
BacillaceaeAnaerobacillus0.00–0.120.00–0.200.466928.11131.40.978
Anoxybacillus0.00–0.190.00–0.610.525928.11131.40.978
Staphylococcus0.00–0.00.00–0.250.090412.51028.60.188
Enterococcaceae0.00–0.170.00–0.380.502825617.10.625
Enterococcus0.430–82.50.100–56.00.005†2990.62982.90.480
Lactobacillus0.0580–26.50.0140–69.60.1152371.92057.10.317
StreptococcaceaeLactococcus0.00–0.290.00–0.570.2831031.3617.10.286
Streptococcus0.0200–2.50.00–33.20.0992371.91748.60.090
TuricibacteraceaeTuricibacter0.00–0.130.00–4.860.372928.11234.30.780
o__Clostridiales0.00–1.30.00520–0.470.064928.11954.30.055
o__Clostridiales0.320–10.10.150–14.90.9002887.52982.90.736
Clostridiaceae0.00–0.00.00–0.200.02800.0514.30.054
Candidatus Arthromitus0.00–0.420.00–0.570.6188251028.60.957
Clostridium9.500.22–77.513.00.40–83.20.2803210035100–
Pseudoramibacter Eubacterium0.00–0.0940.00–0.0510.0511134.4514.30.101
Lachnospiraceae0.0310–1.060.0240–0.770.6172062.52468.60.791
Lachnospiraceae1.150.033–31.20.990.036–14.60.9303210035100–
Blautia7.570.29–35.07.820.20–30.80.9503210035100–
Clostridium1.330.39–50.10.290–59.90.023321003291.40.240
Coprococcus0.130–2.40.160–2.30.8503093.828800.153
Dorea0.870–15.11.020.015–15.10.7823196.9351000.478
Epulopiscium0.00–0.270.00–2.400.356928.11234.30.780
Hespellia0.00–0.0810.00–0.0850.148825411.40.259
Roseburia0.0340–1.440.00610–2.30.05024751851.40.082
[Ruminococcus]1.680.97–15.41.560.027–9.90.8213210035100–
PeptococcaceaePeptococcus0.0230–3.50.00190–6.070.8371753.11851.40.916
Peptostreptococcaceae0.00820–0.490.00–0.530.3571753.114400.406
Peptostreptococcaceae1.370.11–82.50.500–42.70.380321003497.11.00
Peptostreptococcus0.00–0.630.00–0.120.723928.1925.70.957
Ruminococcaceae0.130–6.460.200–6.50.5182990.63291.41.00
Ruminococcaceae0.0520–5.490.180–4.30.6002578.12674.30.935
Butyricicoccus0.000–6.680.00–0.390.219928.1617.10.433
Faecalibacterium0.00650–51.70.0150–3.600.42016502057.10.734
Oscillospira0.0990–8.30.0510–3.480.6502578.12777.10.844
Ruminococcus0.00950–0.980.0270–0.890.3501856.32674.30.195
Subdoligranulum0.0240–16.60.0180–64.30.6352165.62365.70.803
Acidaminococcus0.00–0.390.00–0.690.4021031.37200.438
Dialister0.0840–18.50.170–10.00.7282371.928800.622
Megamonas0.00–1.000.0220–2.01<0.001†928.12571.4<0.001†
Megasphaera0.310–31.90.340–18.50.9702887.528800.619
Phascolarctobacterium0.00–0.230.00–1.760.465721.91028.60.728
Veillonella0.00–0.140.00–0.470.80439.4411.41.00
[Mogibacteriaceae]0.00–3.330.00–0.0500.6531031.31028.60.978
Mogibacterium0.00–1.270.00–0.930.659618.8822.90.911
Anaerococcus0.00–0.0270.00–0.0110.34426.3514.30.431
Gallicola0.00–0.0350.00–0.0260.916412.5514.31.00
Parvimonas0.00–0.340.00–0.0610.0361340.6617.10.063
Erysipelotrichaceae0.00–0.350.00–0.270.42739.4617.10.48
Bulleidia0.00–3.320.00–4.220.805721.97200.911
Catenibacterium0.00–6.280.00–3.300.9221340.61542.90.950
Clostridium0.00–0.870.00–0.320.9891443.81645.70.933
[Eubacterium]0.430.0054–6.840.270–8.310.660321003291.40.240
p-75-a50.00–0.0700.00–0.0430.603412.538.60.701
Fusobacteria0.00–10.50.00–7.870.5601546.91748.60.916
Fusobacteriaceae0.00–0.190.00–7.370.156515.61028.60.329
Fusobacterium0.00–10.50.00–1.190.9541340.61337.10.967
Proteobacteria0.290.0088–40.80.550.019–0–24.70.3863210035100–
mitochondria0.00–0.00.00–0.310.25626.3514.30.431
Sutterella0.00–4.50.00310–2.00.2611443.81851.40.701
BurkholderiaceaeBurkholderia0.00–0.00.00–0.0710.911721.9822.90.844
Comamonadaceae0.00–0.00.00–0.190.982412.5411.41.00
Ralstonia0.00–0.180.00–0.160.3181443.81131.40.43
Desulfovibrio0.00–0.00.00–0.180.46226.3411.40.675
CampylobacteraceaeCampylobacter0.00–2.00.00320–2.370.2571237.51954.30.258
HelicobacteraceaeHelicobacter0.00–1.10.00–0.790.4221134.414400.824
Anaerobiospirillum0.00–1.30.00–2.510.101618.81234.30.247
Enterobacteriaceae0.0510–40.70.0470–22.60.7822887.53085.71.00
Morganella0.00–0.0880.00–0.300.56239.4514.30.711
HalomonadaceaeHalomonas0.00–0.00.00–0.130.01639.41131.40.055
MoraxellaceaeAcinetobacter0.00–0.0530.00–0.0490.306928.1617.10.433
Pseudomonas0.00220–0.230.00740–0.410.62516502057.10.734
VibrionaceaeVibrio0.00–0.00.00–0.100.02800514.30.054
TM70.00–0.760.00–0.0890.1351031.3617.10.286
c__TM7-30.00–0.0590.00–0.0080.017721.912.90.023
o__CW0400.00–0.570.00–0.0890.244618.838.60.292

Relative abundance and prevalence of taxa (percent of total sequences) observed by amplification and sequencing of 16S rDNA from feces of 67 kittens randomized to treatment with daily E. hirae probiotic (32 kittens) or placebo (35 kittens).

Taxa present at a prevalence of < 10% of kittens in both treatment groups are not shown in the table.

†

Benjamini-Hochberg corrected p-value < 0.25.

Discussion

Some unique features of this study are the high-risk target population (shelter kittens), realistic environmental conditions, preventative strategy, naturally-occurring disease outcome measurement, randomized controlled trial design, and use of a host-origin probiotic. Only 2 prior studies describe use of a host-origin probiotic in cats. The first involved administration of Bifidobacterium pseudocatenulatum previously isolated from the feces of adult healthy cats to other healthy adult client-owned cats (). In the second study, E. faecium and Lactobacillus isolated via endoscopic duodenal aspirate from a healthy cheetah, increased body weight and improved fecal quality when given to 8–13 month old cheetah cubs ().

The present study demonstrates that feline-origin E. hirae, administered as a probiotic, is voluntarily consumed by kittens when added to canned cat food, maintained a stable number of colony forming units of bacteria as kept frozen over the course of the study, and had no adverse effects on body weight gain during administration. Based on results of 16S rRNA gene amplicon sequencing, kittens receiving E. hirae had a significantly higher relative abundance of Enterococcus spp. identified in feces compared to those receiving the placebo. We have previously shown that administration of the same probiotic formulation of E. hirae to purpose-bred kittens results in significant increases in the specific identity of E. hirae in feces using both qPCR and quantitative culture (). Inadequate numbers of samples had sufficient DNA remaining after sequencing to perform a specific qPCR assay for E. hirae in this study.

The incidence of diarrhea among foster kittens over the course of the study was relatively low. Despite this, there was a significantly lower incidence of diarrhea among kittens receiving the probiotic compared to the placebo which supports a beneficial effect. It is interesting that diarrhea in the placebo group involved a number of litters with multiple affected kittens whereas diarrhea occurred only among singletons in the probiotic group. Individual kittens with diarrhea did not undergo additional diagnostic testing as this was outside the scope of the study. Another limitation of our study is that data on daily fecal consistency of kittens was not obtained. Instead, we chose to conservatively define diarrhea based on a severity prompting examination of the kitten by the shelter veterinary team. This is the precedent means by which the shelter recognizes and treats diarrhea in foster kittens and best represents the subset of kittens at risk for diarrheal morbidity and mortality in the population. An experimental design requiring foster care providers to first witness and then assign a daily fecal score to individual kittens was considered unfeasible and a probable deterrent to participation in the study. Given the lower incidence of diarrhea in kittens receiving the probiotic in this investigation, a targeted study to examine fecal quality of kittens receiving E. hirae is of interest.

Our interest in E. hirae as a probiotic was based on identification of this species as the predominant mucosa-associated flora in the small intestine of healthy kittens (). In healthy kittens, E. hirae are frequently observed to colonize the epithelial surface (, ). Unknown to us at the time of our observation was a fairly robust preexisting attraction to use of E. hirae as a probiotic. Numerous isolates of E. hirae, ranging in origin from the rumen (), intestinal tract of ocean () and freshwater fish () to goats milk () have been studied for probiotic effects both in-vivo and in-vitro. E. hirae is bile salt and acid tolerant () which promotes its survival through the upper gastrointestinal tract () and demonstrate cell wall hydrophobicity which enables it to colonize by interaction with host intestinal epithelial cells (, , ). A lipoteichoic acid of E. hirae has been demonstrated to ameliorate the loss of barrier function caused by the pro-inflammatory cytokine TNF-α on Caco-2 cells by modulating the expression of tight junction regulatory proteins (). We have previously shown that administration of E. hirae to purpose-bred kittens mitigates the increase in intestinal permeability resulting from experimental infection with feline-origin EPEC (). Other beneficial mechanisms attributed to E. hirae include free radical scavenging and lipase activity (, ).

We previously observed a decrease in prevalence of mucosa-associated E. hirae in the ileum of terminally-ill foster kittens. Loss of E. hirae was associated with colonization by virulent genotypes of E. faecalis and ability of EPEC to attach to the intestinal epithelium (). The enterococci, including E. hirae (, , ), produce peptides (i.e., bacteriocins) with a wide range of antibacterial and in some cases antiviral activity (). Because bacteriocins can promote niche competition with intestinal pathogens, their production is considered to be a favorable attribute of a probiotic (). For example, when administered to hybrid catfish as a probiotic, E. hirae protected against infection by Aeromonas hydrophila, enhanced disease protection, and stimulated immunity-related gene expression (). Similarly, we sought to explore whether administration of feline-origin E. hirae probiotic to foster kittens might augment host defense by promoting clearance of existing infectious agents or preventing acquisition of new infectious agents during the course of foster care. To examine this possibility we used a commercial multiplex qPCR-based assay to determine the presence and quantity of intestinal microbes in fecal samples collected prior to onset and upon conclusion of the study. A high prevalence for carriage of several infectious agents was demonstrated prior to onset of the study, which included feline enteric coronavirus, attaching and effacing E. coli (presence of eae), and panleukopenia virus. A similar prevalence for carriage of these agents in asymptomatic communally housed cats and kittens has been reported by others (–). Overall, our study did not provide strong evidence for an impact of E. hirae on promoting clearance or preventing acquisition of these selected agents during the course of foster care. The only difference observed between the placebo and E. hirae treatment groups was a significant decrease in prevalence of panleukopenia virus DNA among kittens that received the placebo (p = 0.002) and a non-significant decrease among kittens that received E. hirae (P = 0.054). Interpretation of this finding is confounded by concurrent vaccination of kittens with modified-live panleukopenia virus vaccine every 2–3 weeks over the course of the study and the fact that no kittens were clinically unwell at the time of fecal PCR testing. In effort to distinguish between positive PCR results due to modified live virus replication vs. those arising from clinical infection, the reference laboratory applies an a-priori cut-off for positive test results at a Cp value of ≤26. All panleukopenia PCR-positive samples in the study met these criteria. However, recent studies have shown that 21.6% of healthy cats and kittens will shed panleukopenia virus in the positive test result range within 7 days of vaccination (). Remarkably, the percentage of asymptomatic kittens testing positive for panleukopenia virus DNA in this study is nearly identical at 21.8%. Accordingly it is not unreasonable to conclude that vaccine-origin panleukopenia was largely responsible for PCR-positive results in these healthy kittens. In interpreting the results of the infectious disease testing in this study, it is important to recognize that the presence of microbial DNA does not prove agent viability or measure a potential impact of E. hirae on microbial behavior. For example we have previously documented that E. hirae can mitigate intestinal injury caused by EPEC infection in kittens while having no measureable impact on EPEC shedding ().

Based on 16S rRNA gene amplicon sequencing, administration of the E. hirae probiotic did not alter the predominant bacterial phyla present in feces which were similar to prior descriptions of the fecal microbiota in cats. These phyla were represented mainly by Firmicutes, Bacteroidetes, Proteobacteria, Actinobacteria, and Fusobacteria (, –). The presence of E. hirae also did not affect the overall number or diversity of different taxa present in the microbial community which is a desirable attribute of a probiotic. Kittens receiving E. hirae had a significant decrease in presence and relative abundance of the genus Megamonas (phylum Firmicutes, family Veillonellaceae). Because Megamonas was present at low abundance among the fecal microbiota from kittens, the clinical significance of decreased representation in kittens receiving E. hirae is unclear. Arguably the most pertinent location to determine an impact of the E. hirae probiotic would be on composition of the small intestinal microbiota where mucosa-associated E. hirae dominate the enterococci and from which the probiotic strain was isolated. Obtaining samples from the small intestine of kittens in this study was not feasible using the current study design. In prior reports, dietary supplementation with fermentable prebiotics increased the abundance of Megamonas (–) while decreased in abundance of Megamonas have been described in cats with diarrhea (). Interestingly, Megamonas was increased in feces of cats that were naturally infected with the diarrheal pathogen Tritrichomonas foetus, but decreased in feces of cats experimentally infected with the same pathogen (). These collective observations suggest a vulnerability of Megamonas to specific, albeit unidentified changes in the feline intestinal tract. Because the enterococci are lactic acid producers and Megamonas utilizes lactate to form propionate () we would not have predicted a decrease in Megamonas in kittens receiving the E. hirae probiotic. In fact, in dogs fed a raw diet, increases in lactate production were positively correlated with abundance of Megamonas (). Ultimately a deeper understanding of any interaction between E. hirae and Megamonas and the clinical relevance of our observation will require additional study.

In conclusion, a decreased incidence of diarrhea associated with preventative administration of feline-origin E. hirae to foster kittens in this study supports a rationale for use of E. hirae for disease prevention in this large, young, and vulnerable population.

Funding

This work was supported by grants from the Winn Feline Foundation and PetSmart Charities® (Grant W18-002) and a Research Innovation and Seed Grant from North Carolina State University. The UNC Microbiome Core is supported in part by P30 DK034987 Center for Gastrointestinal Biology and Disease (CGIBD) and the UNC Nutrition Obesity Research Center (NORC P30 DK056350).

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Statements

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/, PRJNA776011.

Ethics statement

The animal study was reviewed and approved by Institutional Animal Care and Use Committee of North Carolina State University (Protocol Number: 18-067-O).

Author contributions

JG, SJS, JB-B, MA-P, and GD designed the project and secured funding. JG, SHS, SW, EW, and ME performed the experiments. JG, SJS, MA-P, AS, and GD supervised the experiments. JG, SJS, SHS, MA-P, ME, and AS performed bioinformatics analysis. JG, SJS, JB-B, MA-P, ME, AS, and JB interpreted the data. JG, SJS, JB-B, SW, MA-P, ME, AS, and JB wrote sections of the manuscript. All authors have read and approved the final manuscript.

Acknowledgments

The authors thank Jennifer Federico, Sarah McCain, Joseph Threadcraft, Joanne Duda, Diadra Cooley, Cindy Lynch, Greg Ceseretti, Elliot Ser, Maddie Mesaros, Rebecca Williams, Melissa Reel, Abigail Hiley, Gabby Vasek, Lindsay Stone, and Taylor Raines for their help with this study.

Conflict of interest

Authors ME and AS were employed by IDEXX Laboratories, Inc. Author JB was employed by Vet Med Labor GmbH Division, IDEXX Laboratories, Inc. The remaining 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.

References

Summary

Keywords

16S rRNA gene amplicon sequencing, diarrheal pathogens, growth and survival, infection, polymerase chain reaction, shelter medicine

Citation

Gookin JL, Strong SJ, Bruno-Bárcena JM, Stauffer SH, Williams S, Wassack E, Azcarate-Peril MA, Estrada M, Seguin A, Balzer J and Davidson G (2022) Randomized placebo-controlled trial of feline-origin Enterococcus hirae probiotic effects on preventative health and fecal microbiota composition of fostered shelter kittens. Front. Vet. Sci. 9:923792. doi: 10.3389/fvets.2022.923792

Received

19 April 2022

Accepted

02 November 2022

Published

17 November 2022

Volume

9 - 2022

Edited by

Andrea Lauková, Slovak Academy of Sciences, Slovakia

Reviewed by

Jan S. Suchodolski, Texas A&M University, United States; Martin Fraga, National Institute for Agricultural Research (INIA), Uruguay

Updates

Copyright

*Correspondence: Jody L. Gookin

This article was submitted to Veterinary Infectious Diseases, a section of the journal Frontiers in Veterinary Science

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics