- Research article
- Open Access
Evaluation of the impact of refrigeration on next generation sequencing-based assessment of the canine and feline fecal microbiota
BMC Veterinary Research volume 10, Article number: 230 (2014)
Evaluation of factors that might impact microbiota assessment is important to avoid spurious results, particularly in field and multicenter studies where sample collection may occur distant from the laboratory. This study evaluated the impact of refrigeration on next generation sequence-based assessment of the canine and feline fecal microbiota. Fecal samples were collected from seven dogs and ten cats, and analysed at baseline and after 3, 7, 10 and 14 days of storage at 4°C.
There were no differences in community membership or population structure between timepoints for either dogs or cats, nor were there any differences in richness, diversity and evenness. There were few differences in relative abundance of phyla or predominant genera, with the only differences being significant increases in Actinobacteria between Days 0-14 (P = 0.0184) and 1-14 (P = 0.0182) for canine samples, and a decrease in Erysipelotrichaceae incertae sedis between Day 0 and Day 7 (median 4.9 vs 2.2%, P = 0.046) in feline samples.
Linear discriminant analysis effect size and indicator analysis identified a small number of genera that were over-represented in, or defining characteristics of, Day 14 samples. These were predominantly Proteobacteria and Actinobacteria, with Psychrobacter and Arthrobacter enriched in both canine and feline Day 14 samples.
Storage for at least 14 days at 4°C has limited impact on culture-independent assessment of the canine and feline fecal microbiota, although changes in some individual groups may occur.
Advances in next generation sequencing and bioinformatics have revolutionized the study of complex microbial populations. Recent studies in dogs and cats have characterized the fecal, intestinal, oral and skin microbiotas -, and provided important insights into both the composition of the microbiota and its relationship with various diseases. However, microbiota assessment is not infallible and there are many steps in the process that could potentially impact results. One is the time from sample collection to processing, an important factor for field studies where samples may be collected remotely from the laboratory. An understanding of the impact of storage conditions and time is important for study design and interpretation.
There has been limited study of the impact of storage on the fecal microbiota. A study of human fecal samples reported limited impact of storage at room temperature for 24 h or at -80°C for six months . Similarly little impact was identified in human vaginal samples stored at -20°C and -80°C . A study of human and soil samples detected no influence of storage, including refrigeration, on microbial population structure or diversity, but some shifts in relative abundances of different taxa , results that were similar to a later study of soil bacteria .
Since field and multicenter studies often involve delays from sample collection to processing, it is critical to understand the potential impacts of storage. Refrigeration (4°C) is the most convenient temperature for temporary storage and shipping, and it is important to understand if short-term refrigeration will significantly impact downstream analysis. Species- and sample (i.e. body site)-specific study is required because of the differences in the microbiota between species and body sites, and the potential that there could be different impacts of storage on populations within these samples. The objective of this study was to evaluate the impact of refrigeration on assessment of the canine and feline fecal microbiota.
Fecal samples were collected from seven dogs and ten cats from a local animal shelter immediately after defecation. Dogs were clinically normal with no history of antimicrobial exposure or gastrointestinal disease, although medical histories were limited based on the nature of the population. Samples were stored in plastic fecal containers at 4°C for up to 2 hours prior to arrival at the laboratory. Immediately after arrival, samples were manually homogenized and separated into five aliquots. One aliquot was processed immediately while the other four were stored in a refrigerator at 4°C. One of each aliquots was then tested after 3, 7, 10 and 14 days of refrigeration.
DNA was extracted using a commercial kita, and DNA quantity and quality were assessed by spectrophotometryb. The V4 region of the 16S rRNA gene was then amplified using the primers S-D-Bact-0564-a-S-15 (5'-AYTGGGYDTAAAGNG-3') and S-D-Bact-0785-b-A-18 (5'-TACNVGGGTATCTAATCC-3') . The amplicon library was purified with Agencourt AMPure XP beadsc with slight modification to the manufacturer’s protocol. Briefly, 72 μL of AMPure beads was added to 20 μL of library and incubated for 10 min at room temperature. Samples were washed twice with 80% ethanol, and eluted with 20 μL of PCR-grade H2O. Purified samples were quantified by spectrophotometry, evaluated by electrophoresis on a 1% agarose gel, and diluted to 5 ng/μL. Sequencing was performed using an Illumina MiSeq with 2X250 chemistry.d
Mothur v33.3 was used for analysis . After paired end reads were assembled, sequences were aligned with the Silva 16S rRNA reference database  and any sequences not consistent with the target amplicon size (240 bp), containing any ambiguous base calls or long runs (>8 bp) of holopolymers, or that did not align with the correct 16S rRNA gene region were removed. Chimeras were detected using uchime  and removed. Taxonomy was assigned using the RDP taxonomy database . Sequences were then binned into operational taxon units (OTUs) at a 3% dissimilarity level.
Subsampling was performed to normalize sequence number for analyses . This involved random selection of a number of sequences from each sample that corresponded to the smallest number of sequences from an individual sample. Population diversity (inverse Simpson’s index), evenness (Shannon’s evenness index) and richness (Chao1) were calculated and compared between groups using Wilcoxon and Steel-Dwass tests. Linear discriminatory analysis effective size (LefSe)  analysis and indicator analysis  were performed.
Community membership was described using the classical Jaccard index, while population structure (evaluation of membership and relative abundance of members) was evaluated using the Yue & Clayton index of dissimilarity and Bray-Curtis index. Unifrac was used to compare these indices between groups . Principal coordinate analysis (PCoA) and random forest analysis were also performed. The relative abundance of all phyla were compared between groups, along with the relative abundance of all genera accounting for ≥ 1% of sequences in the Day 0 or Day 14 samples using the Steel-Dwass test. A P <0.05 was considered significant for all comparisons.
A total of 3,722,656 sequences passed all quality control filters, with a median sequence count of 41970 per sample. Subsampling at 5177 sequences per sample was performed to normalize data for subsequent analysis.
There were no differences in community membership (Jaccard index, P = 0.084) or population structure (Yue & Clayton P = 0.41, Bray-Curtis P = 0.77). There were also no differences in richness, diversity and evenness between the different timepoints (all P > 0.05).
Median relative abundances of phyla and genera accounting for 1% of more of the Day 0 or Day 14 microbiota are presented in Figures 1 and 2. Significant increases in Actinobacteria between Days 0-14 (median 0.23 vs 1.14%, P = 0.0184) and 1-14 (0.16 vs 1.14%, P = 0.0182) were identified. No significant differences were identified in relative abundance of predominant genera.
LEfSe identified 14 genera that were enriched in the Day 14 group compared to Day 0 (Figure 3). Four indicator OTUs were identified for Day 14 vs Day 0 samples, Rhizobium, Psychrobacter, Serratia and Delftia, all of which are Proteobacteria. Random forest modeling did not provide any further indication of an impact of storage time on the microbiota (data not presented).
No differences in community membership (P = 0.24) or population structure (Yue & Clayton P = 1.0, Bray-Curtis P = 0.74) were identified by unifrac. There were also no differences in richness, evenness or diversity between groups.
No significant differences in phyla and predominant genera were identified, with the exception of a decrease in Erysipelotrichaceae incerate sedis between Day 0 and Day 7 (4.9 vs 2.2%, P = 0.046).
LEfSe results comparing Days 0 and 14 are presented in Figure 4. Indicator analysis showed a similar dominance of Proteobacteria and Actinobacteria as Day 14 indicators (Table 1). When Days 0 and 7 were compared, no significant indicators were identified and there were only two significant OTUs identified by LEfSe. These were both Firmicutes, with Eubacterium significantly associated with Day 0 and Oscillospira associated with Day 7. Principal coordinate analysis data are presented in Figure 5.
Results of this study indicate that there are limited changes in the fecal microbiota of dogs and cats with short-term refrigeration. Most changes that were evident were identified at Day 14, with very little apparent impact of storage of seven days duration. Even by Day 14, changes were limited, with no differences in diversity, evenness and richness.
Various other microbiota assessment tools are available and consideration of the potential for differential impacts of storage is important. Assessment of the microbiota often involves evaluation ecological indices that evaluate OTUs that are present (membership) and OTUs that are present along with their relative abundances (structure). No alterations of these indices were noted. LEfSe is another method that is useful for analysis of high dimension data such as were present here, and which identifies OTUs that are most likely to explain differences between groups . Canine samples after 14 days of storage were enriched in 14 genera, predominantly members of Actinobacteria and Proteobacteria. A similar pattern was noted in feline samples, with enrichment of certain Actinobacteria and Proteobacteria by Day 14. Two genera, Psychrobacter and Arthrobacter, were enriched in both canine and feline Day 14 samples.
Indicator analysis is an ecological tool that identifies members (in this case, OTUs) that `define’ a population, based on their presence and relative abundance. Indicator analysis yielded similar results to LEfSe, with three of the four genera identified as Day 14 indicators in dogs also identified by LEfSe. There was less agreement in of indicator analysis and LEfSe for feline samples, yet significant Day 14 results from both methods were dominated by Proteobacteria and Actinobacteria.
The increase in certain Proteobacteria and Actinobacteria was presumably a result of growth of these members during storage, as opposed to loss of other components. Psychrobacter can grow at a wide range of temperatures, including 4°C , something that likely accounts for its increased presence in both canine and feline Day 14 samples. A study of stored soil samples reported over-representation of Rhizobiales, Alphaproteobacteria and Thermomicrobia after refrigeration , some of which were also identified here.
Despite the changes that were noted, there was limited overall impact on the microbiota with fourteen days of refrigeration. However, these data suggest that consideration should be given to any lower level (e.g. genus) taxonomic differences that are noted in samples that have undergone storage, particularly those involving Proteobacteria and Actinobacteria.
This study evaluated dogs from an animal shelter. While historical data for these animals are limited, that should have limited impact on the objectives of this study, since it was designed to evaluate the effects of storage, not compare groups or define the canine microbiota. This study only evaluated samples for 14 days, so no assurance can be given that more profound changes do not occur with longer storage. Yet, it is reasonable to assume that samples could be processed or frozen within this timeframe, even with multicenter and field studies. The relatively small sample size that could have hampered identification of some differences. Another factor that should be considered is the potential for intra-sample variation. Another consideration is that if fecal samples are non-homogenous, some minor changes could potentially be the result of inherent intra-sample variation, an area that has not been carefully studied.
Nonetheless, there were few changes identified in the fecal microbiota in samples stored at 4C for 14 days, something that is consistent with previous studies of different sample types and storage conditions -.
Short term (up to 14 day) refrigeration should have limited impact on studies of the canine and feline fecal microbiota. However, the potential for changes in some members of the microbiota must be considered during study design and interpretation, particularly when there will be a focus on individual genera, as opposed to broader population-based analysis. Thus, it appears to be reasonable to use short term storage of canine and feline fecal samples, if required, but to strive for similar storage conditions between groups to remove any potential impact on subsequent analyses.
Availability of supporting data
The dataset supporting the results of this article is available at the MG-RAST metagenomics analysis server (project 9714, http://metagenomics.anl.gov).
aE.Z.N.A. Stool DNA Kit, Omega Bio-Tek Inc., Doraville, Georgia, USA.
bNanodrop, Roche, Mississauga, Canada.
cBeckman Coulter Inc., Mississauga, Canada.
dIllumina, San Diego, USA.
JSW designed the study. JSW coordinated fecal collection. MJ performed laboratory testing. JSW prepared the manuscript. Both authors read and approved the final manuscript.
Linear discriminant analysis effective size
Operational taxon unit
Principle coordinate analysis
Garcia-Mazcorro JF, Dowd SE, Poulsen J, Steiner JM, Suchodolski JS: Abundance and short-term temporal variability of fecal microbiota in healthy dogs. Microbiol Open. 2012, 1 (3): 340-347. 10.1002/mbo3.36.
Handl S, German A, Holden S, Dowd S, Steiner J, Heilmann R, Grant R, Swanson K, Suchodolski J: Faecal microbiota in lean and obese dogs. FEMS Microbiol Ecol. 2013, 84 (2): 332-343. 10.1111/1574-6941.12067.
Suchodolski JS, Dowd SE, Wilke V, Steiner JM, Jergens AE: 16S rRNA gene pyrosequencing reveals bacterial dysbiosis in the duodenum of dogs with idiopathic inflammatory bowel disease. PLoS One. 2012, 7 (6): e39333-10.1371/journal.pone.0039333.
Sturgeon A, Pinder SL, Costa MC, Weese JS: Characterization of the oral microbiota of healthy cats using next-generation sequencing.Vet J. doi:10.1016/j.tvjl.2014.01.024.
Rodrigues Hoffmann A, Patterson AP, Diesel A, Lawhon SD, Ly HJ, Elkins Stephenson C, Mansell J, Steiner JM, Dowd SE, Olivry T, Suchodolski JS: The skin microbiome in healthy and allergic dogs. PLoS One. 2014, 9 (1): e83197-10.1371/journal.pone.0083197.
Carroll IM, Ringel-Kulka T, Siddle JP, Klaenhammer TR, Ringel Y: Characterization of the fecal microbiota using high-throughput sequencing reveals a stable microbial community during storage. PloS One. 2012, 7 (10): e46953-10.1371/journal.pone.0046953.
Bai G, Gajer P, Nandy M, Ma B, Yang H, Sakamoto J, Blanchard MH, Ravel J, Brotman RM: Comparison of storage conditions for human vaginal microbiome studies. PloS One. 2012, 7 (5): e36934-10.1371/journal.pone.0036934.
Lauber CL, Zhou N, Gordon JI, Knight R, Fierer N: Effect of storage conditions on the assessment of bacterial community structure in soil and human-associated samples. FEMS Microbiol Lett. 2010, 307 (1): 80-86. 10.1111/j.1574-6968.2010.01965.x.
Rubin BER, Gibbons SM, Kennedy S, Hampton-Marcell J, Owens S, Gilbert JA: Investigating the impact of storage conditions on microbial community composition in soil samples. PLoS One. 2013, 8 (7): e70460-10.1371/journal.pone.0070460.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Peña AG, Goodrich JK, Gordon JI, Huttley GA, Kelley ST, Knights D, Koenig JE, Ley RE, Lozupone CA, McDonald D, Muegge BD, Pirrung M, Reeder J, Sevinsky JR, Turnbaugh PJ, Walters WA, Widmann J, Yatsunenko T, Zaneveld J, Knight R: QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010, 7 (5): 335-336. 10.1038/nmeth.f.303.
Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, Sahl JW, Stres B, Thallinger GG, Van Horn DJ, Weber CF: Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009, 75 (23): 7537-7541. 10.1128/AEM.01541-09.
Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, Peplies J, Gl°Ckner FO: The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013, 41: D590-D596. 10.1093/nar/gks1219.
Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R: UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011, 27 (16): 2194-2200. 10.1093/bioinformatics/btr381.
RDP. , [http://rdp.cme.msu.edu/index.jsp]
Gihring TM, Green SJ, Schadt CW: Massively parallel rRNA gene sequencing exacerbates the potential for biased community diversity comparisons due to variable library sizes. Env Microbiol. 2012, 14 (2): 285-290. 10.1111/j.1462-2920.2011.02550.x.
Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, Huttenhower C: Metagenomic biomarker discovery and explanation. Genome Biol. 2011, 12 (6): R60-10.1186/gb-2011-12-6-r60.
Dufrêne M, Legendre P: Species assemblages and indicator species: the need for a flexible asymmetrical approach. Ecol Monogr. 1997, 67: 345-366.
Lozupone C, Hamady M, Knight R: UniFrac-an online tool for comparing microbial community diversity in a phylogenetic context. BMC Bioinformatics. 2006, 7: 371-10.1186/1471-2105-7-371.
De Filippis F, La Storia A, Villani F, Ercolini D: Exploring the sources of bacterial spoilers in beefsteaks by culture-independent high-throughput sequencing. PLoS One. 2013, 8 (7): e70222-10.1371/journal.pone.0070222.
This study was supported by the Ontario Veterinary College Pet Trust.
The authors declare that they have no competing interests.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Weese, J.S., Jalali, M. Evaluation of the impact of refrigeration on next generation sequencing-based assessment of the canine and feline fecal microbiota. BMC Vet Res 10, 230 (2014). https://doi.org/10.1186/s12917-014-0230-7