- Methodology article
- Open Access
Comparison of stool versus rectal swab samples and storage conditions on bacterial community profiles
BMC Microbiology volume 17, Article number: 78 (2017)
Sample collection for gut microbiota analysis from in-patients can be challenging. Collection method and storage conditions are potential sources of variability. In this study, we compared the bacterial microbiota from stool stored under different conditions, as well as stool and swab samples, to assess differences due to sample storage conditions and collection method.
Using bacterial 16S rRNA gene sequence analysis, we compared the microbiota profiles of stool samples stored and collected under various conditions. Stool samples (2 liquid, 1 formed) from three different patients at two hospitals were each evaluated under the following conditions: immediately frozen at -80°C, stored at 4°C for 12-48 hours before freezing at -80°C and stored at -20°C with 1-2 thaw cycles before storage at -80°C. Additionally, 8 stool and 30 rectal swab samples were collected from 8 in-patients at one hospital. Microbiota differences were assessed using the Yue and Clayton dissimilarity index (θYC distance) and analysis of molecular variance (AMOVA).
Regardless of the storage conditions, the bacterial communities of aliquots from the same stool samples were very similar based on θYC distances (median intra-sample θYC distance: 0.035, IQR: 0.015-0.061) compared to aliquots from different stool samples (median inter-sample θYC distance: 0.93, IQR: 0.85-0.97) (Wilcoxon test p-value: <0.0001). For the stool and rectal swab comparison, samples from different patients, regardless of sample collection method, were significantly different (AMOVA p-values: <0.001-0.029) compared to no significant difference between all stool and swab samples (AMOVA p-value: 0.976). The θYC dissimilarity index between swab and stool samples was significantly lower within individuals (median 0.17, IQR: 0.10-0.27) than between individuals (median 0.93, IQR: 0.85-0.97) (Wilcoxon test p-value: <0.0001), indicating minimal differences between stool and swab samples collected from the same individual over the sampling period.
For gastrointestinal microbiota studies based on bacterial 16S rRNA gene sequence analysis, interim stool sample storage at 4 °C or -20 °C, rather than immediate storage at -80 °C, does not significantly alter results. Additionally, stool and rectal swab microbiotas from the same subject were highly similar, indicating that these sampling methods could be used interchangeably to assess the community structure of the distal GI tract.
The diverse communities of microorganisms that compose the human gut microbiota play key roles in health and disease. Advances in sequencing technology have facilitated the wide use of bacterial 16S rRNA-encoding gene sequence analysis for the identification of bacterial lineages as well as their relative abundances in microbial communities. Alterations in the gut microbiota are associated with numerous diseases including cardiovascular disease, inflammatory bowel diseases (IBD) and colorectal cancer as well as increased susceptibility to infections [1–5]. Carriage of multidrug-resistant opportunistic enteric pathogens, such as vancomycin-resistant enterococcus (VRE) and extended spectrum beta-lactamase (ESBL) producing Enterobacteriaceae, has also been associated with changes in intestinal bacterial communities among hospital patients . This observation has resulted in an interest in understanding the gastrointestinal microbiota features that may predispose patients to colonization with multidrug-resistant organisms (MDROs), which could lead to the development of interventions to prevent MDRO colonization and subsequent infection.
Much of the work in gastrointestinal microbiota analyses from human subjects has been done using stool samples (e.g. ). In hospitals, patient-level factors, such as fecal incontinence, and facility-level factors, such as heavy nursing workloads, can make collection of freshly passed stool challenging or impractical. In contrast, collection of rectal swab samples for surveillance cultures among hospitalized patients to determine colonization with MDROs is a routine infection control practice [8–12].
Fecal samples for routine culture are often preserved using chemicals, refrigeration, or freezing depending on the testing that will occur on the specimen. In microbiota analyses, it is important to utilize a procedure that will minimize DNA alteration in samples prior to the analysis. If DNA extraction is not done immediately after collection, the gold standard is to store specimens at -80 °C. In some clinical settings, however, there may not be immediate access to an ultralow-temperature freezer, and samples may need to be transported at higher temperatures before reaching the lab.
In this study, we compared the bacterial profiles of rectal swab samples to stool samples collected from patients at one hospital. We also assessed the effects of common storage conditions on the composition of the fecal microbiota. The objective of our study was to determine the effects of storage and sampling method on the gastrointestinal microbiota.
Evaluation of stool storage conditions
We first compared the effects of different storage conditions on the microbial composition of stool samples (collected into a sterile container without the use of a preservative) from three patients. Sample A was a diarrheal stool from a patient who had been in Hospital A for 11 days. The patient was currently receiving enteral nutrition via a gastrointestinal tube and had received intravenous colistin, daptomycin, and vancomycin. Sample B was a clear, watery stool from a patient who had been hospitalized for seven days at Hospital B. The patient had received 7 doses of oral levofloxacin and 1 dose of intravenous vancomycin. Sample C was a formed stool from an outpatient at Hospital B who had completed a two-week course of oral clarithromycin and metronidazole approximately 2 months prior to sample collection.
Each sample was split into 15 aliquots and tested in triplicate under the 5 storage conditions described in Table 1. For stool samples, dilution of DNA often improves PCR amplification. So, in addition to comparing storage conditions, we also tested the effect of dilution on community analysis by comparing DNA from the same samples that was undiluted and diluted 1:10 for the PCR. We then used 16S rRNA gene sequence analysis to assess the microbial community profile of each sample condition. Aliquots that did not amplify or were poorly sequenced (<1000 sequences per sample) were not included in the analysis. After sequence processing we obtained 1,882,278 sequences from the V4 region of the 16S rRNA gene from 74 sample aliquots with an average of 25,436 ± 10,999 (SD) sequences per sample.
The bacterial community composition was not strongly affected by storage condition or DNA dilution (Fig. 1). However, the gut microbiota of each patient was significantly different from that of the other patients based on θYC distances (AMOVA p-value: <0.001 for all comparisons) (Fig. 2a). Although each patient displayed a distinct microbiota, community structure was markedly similar within each patient for all storage conditions and dilutions (Figs. 1 and 2). Regardless of storage condition or dilution, θYC distances between microbiota of aliquots from the same samples (median: 0.035, IQR: 0.015-0.061) were significantly lower than θYC distances between microbiota of aliquots from different samples (median: 0.93, IQR: 0.85-0.97) (Wilcoxon test p-value: <0.0001) (Fig. 2b). Additionally, overall community richness (the number of OTUs per sample) did not differ significantly between different storage conditions or dilutions within each patient (Kruskal-Wallis test) (OTUs per sample: A median: 89, IQR: 83-99.8; B median: 62.5, IQR: 58-70.8; C median: 112.5, IQR: 107-118.8). Our results indicate high similarity between the bacterial communities of aliquots from the same sample even when different storage conditions or dilutions were used.
Rectal swab compared to stool specimens
To evaluate stool versus swab sample collection, we collected one stool sample and multiple rectal swabs during the next 24 to 27 h from 8 patients each: 6 women and 2 men, with a median age of 55 years (IQR: 47-58 years). To detect possible contamination, an unused swab and a reagents only/no sample control were processed through DNA isolation and PCR with the stool and swab samples. These control samples did not yield a PCR product that was visible on a gel (data not shown). After sequence processing we obtained 754,371 sequences from the V4 region of the 16S rRNA gene from 8 stool specimens and 30 swab samples with an average of 19,852 ± 8484 (SD) sequences per sample.
Overall, the bacterial community structure was similar between the freshly passed stool and the rectal swabs collected at various time points within each patient (Figs. 3 and 4). PCoA of θYC distances indicated that the microbiota of stool samples and rectal swabs clustered by subject (Fig. 4a). In addition, regardless of sample collection method, there were significant differences between all subjects (AMOVA p-value: <0.001-0.029) compared to no overall difference between all stool and swab samples from within subjects (AMOVA p-value: 0.976). The θYC distance between swab and stool samples was significantly lower within subjects (median: 0.17, IQR: 0.10-0.27) compared to between subjects (median: 0.93, IQR: 0.85-0.97) (Wilcoxon test p-value <0.0001) (Fig. 4b). These results indicate that there were minimal differences present between stool and swab samples collected from the same subject over the sampling period up to 27 h after the baseline bowel movement.
Sample s1_2_3 from subject 1 was found to be an outlier, as the bacterial community profile does not correlate with other specimens from subject 1 (Figs. 3 and 4). We attempted to confirm the bacterial community composition of s1_2_3 from the second swab head of the dual swab sample. DNA isolation from the second swab head of s1_2_3 yielded no detectable DNA and PCR of the V4 region of the bacterial 16S rRNA gene did not yield a detectable product (data not shown), suggesting that contamination from another sample may explain these results.
Understanding links between the gastrointestinal microbiota and health in the clinical setting has the potential to improve patient care in the context of infection prevention and beyond. Optimizing study feasibility without altering results is critical for research on the gut microbiota in hospitalized patients. In this study we investigated the effects of sample storage conditions and collection methods on analysis of the gut microbiota.
Accurate analysis of the gastrointestinal microbiota based on bacterial 16S rRNA gene sequences did not require immediate storage of samples at -80 °C. Interim storage (24-48 h) of stool aliquots at temperatures likely to be available in a hospital (4 °C or -20 °C), even with 1 or 2 freeze/thaw cycles, didn’t significantly alter the microbiota. This agrees with a recent study showing minimal changes in the microbiota of stool samples stored without a preservative at -20 °C or 4 °C up to 8 weeks, although fungal growth is a likely complication with extended periods at 4 °C . Storage at warmer temperatures have been shown to alter microbiota from stool [13, 14] and sputum .
The utility of rectal swab cultures for the surveillance of MDROs among hospitalized patients is widely recognized. Rectal swabs are relatively simple samples to collect, require no patient preparation, and can be transported easily from the bedside to the laboratory. In the hospital setting, which often includes medically complex patients and heavy nursing workloads, rectal swabs are more convenient to collect than stool samples. With healthier subjects, swabs can be self-collected with minimal instruction. In our study, we found bacterial communities in individual patients to be highly similar from stool and rectal swab samples. The overall composition of bacterial communities was comparable to a freshly passed stool specimen even in swabs collected up to 27-h after stool passage. A potential limitation of using swab samples compared to stool samples is that the amount of sample collected is smaller. In our study, successful bacterial community analysis was possible from most (29/30) of the rectal swab samples. One rectal swab sample yielded DNA levels too low for accurate bacterial community analysis, likely rendering that sample especially susceptible to contamination. Our findings suggest that rectal swabs are an acceptable and practical proxy for the collection of fecal specimens for stool microbiota analysis. Similarly, a previous study found that rectal swabs were a suitable alternative to stool for analyzing the intestinal microbiota using IS-pro, a method that differentiates bacteria based on internal transcribed spacer (ITS) length and phylum-specific fluorescent primers .
Gastrointestinal microbiota studies based on bacterial 16S rRNA gene sequencing have options for interim sample storage conditions (4 °C or -20 °C vs. -80 °C) and sample collection methods (stool vs. rectal swab) that may increase sampling feasibility in the hospital setting without altering results.
Specimen selection and collection
To assess the effects of different storage conditions on bacterial community profiles, salvaged stool samples submitted to the clinical microbiology laboratory of a 108-bed long-term acute care hospital (hospital A) and a 720-bed tertiary, short-stay acute care hospital (hospital B) in Chicago, IL were tested. To compare rectal swabs with freshly passed stool, a convenience sample of subjects from hospital A (2 women, 6 men) was selected from those patients who were present in the facility on the day of sample collection. Stool that would have been otherwise discarded was collected in a sterile container without preservative from each patient immediately after a bowel movement. Rectal swab samples were collected within 5 min after the bowel movement and 3, 6 and 12-27 h later, by inserting a dual Dacron swab moistened with sterile liquid Stuart medium (Becton Dickenson, Sparks, MD) 1-2 cm past the anal verge and rotating the swab gently 360°. Swab samples were stored in the original swab collection container with liquid Stuart medium. Stool and rectal swab samples were stored up to 27 h at 1-8 °C before being frozen at -80 °C.
Specimen processing for storage conditions analysis
For the storage conditions analysis, each specimen was divided into 15 aliquots to evaluate each of the different storage conditions being tested (Table 1). Aliquots of 0.2 g of stool were prepared in Sarstedt tubes in triplicate for each storage condition. After the samples were subjected to the storage conditions as described, the samples were transferred to an ultralow-temperature freezer.
DNA isolation, library preparation and sequencing
All specimens were shipped overnight on dry ice to the University of Michigan. Samples were transferred to a 96-well bead plate and then submitted to the University of Michigan Microbial Systems Laboratory for DNA isolation and sequencing. DNA was isolated with a PowerMag Soil DNA Isolation Kit (Mo Bio Laboratories, Inc.) using an epMotion 5075 liquid handling system (Eppendorf). The V4 region of the 16S rRNA gene was amplified and sequenced with a MiSeq (Illumina) as described previously . Fastq files were deposited in the SRA (Bioproject: PRJNA317493).
Analysis of 16S rRNA gene sequences
The 16S rRNA gene sequence data was processed and analyzed using the software package mothur (v.1.34.4) and MiSeq standard operating procedure described in Kozich et al. [17–19]. After sequence processing and alignment to the SILVA reference alignment (release 109) , sequences were binned into operational taxonomic units (OTUs) based on 97% sequence similarity using the average neighbor method. Sample aliquots were removed from the analysis if the number of sequences was below 1000. By calculating θYC distances (a metric that takes relative abundances of both shared and non-shared OTUs into account)  between communities and using analysis of molecular variance (AMOVA)  it was possible to determine if there were statistically significant differences between the microbiota of different groups. Principle coordinates analysis (PCoA) was used to visualize the θYC distances between samples. To assess the effect of storage conditions on community composition, θYC distances between aliquots of the same sample were compared to θYC distances between aliquots from different samples using a Wilcoxon test in Prism 6 for Mac OS X (GraphPad Software, Inc.). We also used a Wilcoxon test to compare θYC distances between swab and stool samples from the same subject and between swab and stool samples from different subjects. The taxonomic composition of the bacterial communities was investigated by classifying sequences within mothur using a modified version of the Ribosomal Database Project (RDP) training set [23, 24].
Operational taxonomic unit
Shreiner AB, Kao JY, Young VB. The gut microbiome in health and in disease. Curr Opin Gastroenterol. 2015;31:69–75.
Becattini S, Taur Y, Pamer EG. Antibiotic-induced changes in the intestinal microbiota and disease. Trends Mol Med. 2016;22:458–78.
Leslie JL, Young VB. The rest of the story: the microbiome and gastrointestinal infections. Curr Opin Microbiol. 2015;23:121–5.
Denny JE, Powell WL, Schmidt NW. Local and long-distance calling: conversations between the gut microbiota and intra- and extra-gastrointestinal tract infections. Front Cell Infect Microbiol. 2016;6:41.
Hold GL. Gastrointestinal microbiota and colon cancer. Dig Dis. 2016;34:244–50.
Taur Y, Xavier JB, Lipuma L, Ubeda C, Goldberg J, Gobourne A, Lee YJ, Dubin KA, Socci ND, Viale A, et al. Intestinal domination and the risk of bacteremia in patients undergoing allogeneic hematopoietic stem cell transplantation. Clin Infect Dis. 2012;55:905–14.
Human Microbiome Project Consortium. Structure, function and diversity of the healthy human microbiome. Nature. 2012;486:207–14.
Siegel JD, Rhinehart E, Jackson M, Chiarello L, Healthcare Infection Control Practices Advisory Committee. Management of Multidrug-Resistant Organisms In Healthcare Settings, 2006. http://www.cdc.gov/hicpac/mdro/mdro_0.html: Centers for Disease Control and Prevention; 2006.
Facility Guidance for Control of Carbapenem-resistant Enterobacteriaceae (CRE) – November 2015 Update CRE Toolkit. http://www.cdc.gov/hai/organisms/cre/cre-toolkit/index.html: Centers for Disease Control and Prevention; 2015.
Bonten MJ, Hayden MK, Nathan C, van Voorhis J, Matushek M, Slaughter S, Rice T, Weinstein RA. Epidemiology of colonisation of patients and environment with vancomycin-resistant enterococci. Lancet. 1996;348:1615–9.
Trick WE, Weinstein RA, DeMarais PL, Kuehnert MJ, Tomaska W, Nathan C, Rice TW, McAllister SK, Carson LA, Jarvis WR. Colonization of skilled-care facility residents with antimicrobial-resistant pathogens. J Am Geriatr Soc. 2001;49:270–6.
Hayden MK, Lin MY, Lolans K, Weiner S, Blom D, Moore NM, Fogg L, Henry D, Lyles R, Thurlow C, et al. Prevention of colonization and infection by Klebsiella pneumoniae carbapenemase-producing enterobacteriaceae in long-term acute-care hospitals. Clin Infect Dis. 2015;60:1153–61.
Song SJ, Amir A, Metcalf JL, Amato KR, Xu ZZ, Humphrey G, Knight R. Preservation methods differ in fecal microbiome stability, affecting suitability for field studies. mSystems. 2016;1(3):e00021–16. DOI: 10.1128/mSystems.00021-16.
Shaw AG, Sim K, Powell E, Cornwell E, Cramer T, McClure ZE, Li M-S, Kroll JS. Latitude in sample handling and storage for infant faecal microbiota studies: the elephant in the room? Microbiome. 2016;4:1–14.
Cuthbertson L, Rogers GB, Walker AW, Oliver A, Hafiz T, Hoffman LR, Carroll MP, Parkhill J, Bruce KD, van der Gast CJ. Time between collection and storage significantly influences bacterial sequence composition in sputum samples from cystic fibrosis respiratory infections. J Clin Microbiol. 2014;52:3011–6.
Budding AE, Grasman ME, Eck A, Bogaards JA, Vandenbroucke-Grauls CM, van Bodegraven AA, Savelkoul PH. Rectal swabs for analysis of the intestinal microbiota. PLoS One. 2014;9:e101344.
Seekatz AM, Theriot CM, Molloy CT, Wozniak KL, Bergin IL, Young VB. Fecal microbiota transplantation eliminates Clostridium difficile in a murine model of relapsing disease. Infect Immun. 2015;83:3838–46.
Kozich JJ, Westcott SL, Baxter NT, Highlander SK, Schloss PD. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq illumina sequencing platform. Appl Environ Microbiol. 2013;79:5112–20.
Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75:7537–41.
Schloss PD. A high-throughput DNA sequence aligner for microbial ecology studies. PLoS One. 2009;4:e8230.
Yue JC, Clayton MK. A similarity measure based on species proportions. Commun Stat-Theory Methods. 2005;34:2123–31.
Anderson MJ. A new method for non-parametric multivariate analysis of variance. Austral Ecol. 2001;26:32–46.
Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73:5261–7.
Cole JR, Wang Q, Fish JA, Chai B, McGarrell DM, Sun Y, Brown CT, Porras-Alfaro A, Kuske CR, Tiedje JM. Ribosomal Database Project: data and tools for high throughput rRNA analysis. Nucleic Acids Res. 2014;42:D633–642.
This study was presented in part at IDWeek 2015, San Diego, CA, October 7-11, 2015, abstract number 755.
We would like to thank the study subjects, administration, medical staff, nursing staff, and the infection preventionists at participating hospitals and the Host Microbiome Initiative at the University of Michigan. This research was supported in part through computational resources and services provided by Advanced Research Computing at the University of Michigan, Ann Arbor.
This project was funded through the CDC Prevention Epicenters Program under Cooperative Agreement #U54CK000161-03S1. AMS was supported by grant # UL1TR000433 from the National Center for Advancing Translational Sciences.
Availability of data and materials
Sequence data (fastq files) were deposited in the NCBI’s SRA (Bioproject: PRJNA317493).
CMB processed samples, analyzed sequence data and wrote significant portions of the manuscript. NMM collected and processed samples and was a major contributor to drafting the manuscript. KL collected and processed samples and contributed to study design. AMS processed samples and analyzed sequence data. RAW contributed to discussions of study conception. VBY and MKH played key roles in the conception, design and coordination of the study. All authors edited the manuscript and read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
This study was reviewed and approved by the Institutional Review Board at Rush University Medical Center. Written informed consent was waived because the study presented no more than minimal risk of harm to subjects and involved no procedures for which written consent is normally required outside of the research context.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.