Skip to main content
Help us understand how you use our websites. Take part in our 30 minute study now.

Assessment of the impact of different fecal storage protocols on the microbiota diversity and composition: a pilot study



Fecal samples are currently the most commonly studied proxy for gut microbiota. The gold standard of sample handling and storage for microbiota analysis is maintaining the cold chain during sample transfer and immediate storage at − 80 °C. Gut microbiota studies in large-scale, population-based cohorts require a feasible sample collection protocol. We compared the effect of three different storage methods and mock shipment: immediate freezing at − 80 °C, in 95% ethanol stored at room temperature (RT) for 48 h, and on blood collection card stored at RT for 48 h, on the measured composition of fecal microbiota of eight healthy, female volunteers by sequencing the V4 region of the 16S rRNA gene on an Illumina MiSeq.


Shared operational taxonomic units (OTUs) between different methods were 68 and 3% for OTUs > 0.01 and < 0.01% mean relative abundance within each group, respectively. α and β-diversity measures were not significantly impacted by different storage methods. With the exception of Actinobacteria, fecal microbiota profiles at the phylum level were not significantly affected by the storage method. Actinobacteria was significantly higher in samples collected on card compared to immediate freezing (1.6 ± 1.1% vs. 0.4 ± 0.2%, p = 0.005) mainly driven by expansion of Actinobacteria relative abundance in fecal samples stored on card in two individuals. There was no statistically significant difference at lower taxonomic levels tested.


Consistent results of the microbiota composition and structure for different storage methods were observed. Fecal collection on card could be a suitable alternative to immediate freezing for fecal microbiota analysis using 16S rRNA gene amplicon sequencing.


Gut microbiota bridges the interaction of host and environment and is centrally located in the mechanistic framework of multiple complex intestinal and extraintestinal diseases with multifactorial aetiologies such as inflammatory bowel disease and cardiovascular disease [1]. Accordingly, gut microbiota characterization is an important pillar of precision medicine [2] and is investigated in population-based cohorts [3,4,5]. Analysis of fecal samples is frequently used as a surrogate for studying the gut microbiota with immediate freezing being the gold standard method for preserving the microbial community profile [6]. However, this method introduces logistical difficulties and may not be feasible for all large-scale population-based cohort studies especially in rural and remote settings [7].

Appropriate fecal sampling methods and storage conditions are essential to avoid distorting the original microbial community profile and to minimize the levels of introduced biases in microbiota assessment [6, 8]. Although immediate freezing is commonly considered the gold standard, this approach is not always possible for home sample collection, collection from remote areas, or in large-scale, population-based studies [7]. A variety of alternative methods of sample collection have been proposed in cohort studies, such as the use of preservatives (e.g., ethanol or other solutions) [9,10,11] or collection on fecal occult blood test (FOBT) cards, which are commonly used for fecal collection for colorectal cancer screening [12]. Furthermore, the impact of delayed freezing has been assessed by simulating mock shipments at ambient temperature [13, 14]. These methods inhibit the continued microbial growth after defecation and, thereby preserve the in vivo microbial community composition [7]. Previous studies have detected an overall high agreement between alternative fecal collection/storage methods compared to immediate freezing [13,14,15,16,17]. Sample storage in ethanol or card could be ideal methods for population-based cohorts given their relatively low cost and/or storage volume. Therefore, our objective was to assess different storage methods and mock shipment of fecal sample in terms of the microbial community structure in a pilot study in healthy female volunteers to identify a feasible method for large scale sample collection in our population-based cohorts.


After merging, trimming and quality control, the sequencing yielded total counts of 1,089,909 sequences with a mean of 45,412 (SD = 12,728) per sample. In all samples from healthy female participants, sequences derived from the phyla Firmicutes and Bacteroidetes dominated in all individuals and treatments. At finer taxonomic resolution, 488 operational taxonomic units (OTUs) were identified to have mean relative abundances of at least 0.01% in one of the methods (immediate freezing n = 408, ethanol n = 424, and card n = 397) with only 68% of them being shared in all methods. The shared proportion of OTUs was 3% for rare OTUs (< 0.01% within each group) (Fig. 1a). Approximately 5% of abundant OTUs were unique to each storage condition (Fig. 1a). Similar results were observed for abundant species and only 37% of all species were common to different methods (Additional file 1: Figure S1).

Fig. 1

Effect of storage condition on the fecal microbiota profile. a Venn diagram of abundant and rare operational taxonomic units (OTUs) defined as having mean relative abundance of > 0.01% and < 0.01% within each method, respectively. b α diversity across methods, and c across individuals. d β diversity based on Bray-Curtis dissimilarity for all and rare OTUs (mean relative abundance < 0.01%)

Alpha diversity was assessed using observed OTUs and the inverse Simpson index on the rarefied data. α diversity metrics were not significantly different based on the method of storage (Fig. 1b). However, within each individual, the measured α diversity metrics differed between sample storage treatments. The difference was more pronounced for the number of observed OTUs compared to inverse Simpson index. While some samples had consistent α diversity across different storage methods (e.g., individuals 1, 2, and 8), other samples had lower (individuals 3 and 5) or higher (individual 4) α diversity relative to immediate freezing (Fig. 1c).

We assessed the impact of different storage methods on the overall fecal microbiota composition at the OTU level using Bray-Curtis dissimilarity index (takes into account the relative abundances) and Jaccard distance (examines the presence/absence of taxa). The overall microbial community assessed by Bray-Curtis index was not significantly associated with the storage method both for abundant and rare OTUs (PERMANOVA, p > 0.05). In both cases, between individual differences were statistically significant (PERMANOVA, p < 0.001, R2 = 0.87 and 0.44 for abundant and rare OTUs, respectively, Fig. 1d). Similar results were also observed with Jaccard distance (not shown).

With the exception of Actinobacteria, fecal microbiota profiles at the phylum level were not significantly associated with the storage method. Indeed, inter-individual variation in the relative abundances of different phyla persisted across the different storage methods (Fig. 2a; ANOVA by individual not shown). Actinobacteria was significantly higher in samples collected on card compared to immediate freezing (1.6 ± 1.1% vs. 0.4 ± 0.2%, ANOVA, p = 0.005, Fig. 2b and Additional file 1: Figure S2). This observation was mainly driven by higher relative abundance of Actinobacteria in fecal samples stored on card in individuals 3 and 4 (Fig. 2a and b). The Actinobacteria phylum was dominated by Bifidobacteriaceae and Coriobacteriaceae. While relative abundances of Bifidobacteriaceae remained consistent between different storage methods, Coriobacteriaceae was enriched in samples stored on card (ANOVA, p = 0.002, Additional file 1: Figure S2). In addition, there was a trend in the difference of Firmicutes relative abundances across fecal storage conditions (ANOVA, p = 0.051). Although, there was not a significant difference between storage on card and immediately frozen samples, samples stored in ethanol had lower levels of Firmicutes compared to the card storage (Fig. 2b). Individual-level relative abundances at the family level remained consistent between methods for Bacteroidetes, Firmicutes, and Proteobacteria (Additional file 1: Figure S3). At the genus level, the most abundant taxa were Bacteroides, Prevotella, unclassified Clostridiales, unclassified Lachnospiraceae, and unclassified Ruminococcaceae. Storage method was not significantly associated with relative abundances of the abundant genera (> 0.01% mean relative abundance; Fig. 3). Patterns of microbiota taxonomic structure differences of the most abundant OTUs (> 1% mean relative abundance) did not show prominent differences when comparing storage in 95% ethanol and on card to immediate freezing (Additional file 1: Table S1 and Additional file 1: Table S2). The largest difference was observed for the genus Succinivibrio in one individual (individual 4). The relative abundance of this genera was 14% (3 folds) less abundant in samples stored on card relative to the sample preserved immediately at − 80 °C (6.4% vs. 20.1%; Additional file 1: Table S2).

Fig. 2

Comparison of relative abundances of gut microbiota profile at phylum level a across individuals, and b across methods

Fig. 3

Comparison of relative abundances of abundant genera (> 0.01% mean relative abundance) across individuals and fecal storage methods


In this study, we compared the combined effects of mock shipment and alternative collection and/or preservation methods (collection on card and preservation in 95% ethanol) to the gold standard of immediate freezing at − 80 °C on the observed microbial community structure in fecal samples. All methods delivered remarkably similar microbial community profiles, as indicated by the high similarity of replicate samples from the same individual. We observed that inter-individual variations were stronger than differences between different methods of fecal collection and storage in all microbiota features, including α and β diversity measures and relative abundances of dominant taxa at phylum and genus levels. In other words, the microbiota profiles remained consistent within individuals regardless of the storage method. This observation is in accordance with previous reports [13,14,15,16,17]. Strikingly, however, we observed that only 70% of abundant OTUs were shared among all different storage methods. This proportion was even less for rare OTUs, with 3% being shared, an indication of the stronger impact of the storage method on the observation of low abundance taxa within the fecal microbiota.

While the majority of studies focus on the correlation of results from different storage methods mainly using interclass correlation assessment [14,15,16,17], it is of interest that approximately two-thirds of the most abundant OTUs are shared across different methods in our study suggesting a less prominent impact of storage method on the most abundant taxa. However, only 3% of rare OTUs were detected in all three storage methods which could be attributed to the specifically higher inter-individual variability in less abundant taxa as well as higher sensitivity to the technical variation in the sequencing limit of detection. Our small sample size limits our power to examine the impact of storage conditions on the rare taxa in more details. Even though β diversity was not impacted by the storage method for the rare sub-community in our study, in cases where the research interest is in rare taxa, caution should be exercised in the choice of fecal storage.

Fecal collection methods are often in agreement with immediate freezing on the population level [13,14,15,16,17]. In accordance, we did not observe any differences in microbiota composition and structure on average. However, we observed individual differences specifically in microbial richness and relative abundances between different methods at individual level. For example, in one individual, richness was substantially low in the ethanol sample and higher in the card sample, relative to the frozen control. Conversely, in another individual, diversity was higher in ethanol and card compared to immediate freeze. Additionally, we observed higher relative abundance of Actinobacteria in samples stored on card in two individuals. Although, this observation could be attributed to different sample handling for these samples, the possibility of different behaviour of microbial communities in different storage conditions based on their original profile and internal dynamic cannot be excluded.

The main strength of population-based studies is a large sample size (usually on the scale of 50,000 to 200,000 samples), the availability of rich and comprehensive metadata, and longitudinal follow up [18]. Until recently, however, fecal samples were not routinely bio-banked in population-based cohorts. This trend is changing [19], leading to an interest in proper storage techniques. Samples collected on cards such as FOBT cards occupy considerably less freezer space and are preferable from a logistics perspective. However, other considerations need to be addressed, such as yield and quality of genomic DNA extracts. Under ideal conditions, where microbiota profiling from cards is consistent with that of immediate freezing, samples collected on these cards will still pose a technical challenge for the initial step of the DNA extraction by hindering automation. However, novel technical advances such as the use of specific punching tools may overcome such difficulty.

Although, we and others have observed consistent results of DNA-based microbiome analysis using alternative collection and storage methods, the performance of these approaches for other methods of gut microbiota interrogation such as metatrasncriptomics, metaproteomics, and metabolomics is not very commonly studied. One of the few available studies identified a subset of transcripts which were differentially abundant in alternative methods of sample storage [13]. Currently, DNA-based methods such as 16S rRNA gene sequencing is the method of choice for many projects as a first step of hypothesis generating which could then be followed up with more advanced techniques with higher taxonomic resolution. It could also be complemented with other readouts such as RNA and protein to be able to infer biological functions of the microbiota [20]. Although, fecal collection in ethanol was found to be suitable for metatrasncriptomics [13], preservation of RNA on card has not been investigated. However, card was shown to be suitable for fecal metabolomics [21]. Whether samples collected on card are also useful for bacterial culturomics and investigation of nonbacterial components of the gut microbiota (fungi, viruses, and protozoa) is another open question.

Similar to immediately frozen fecal sample, standardised protocols need to be developed for alternative methods of fecal collection and storage to minimize sources of lab-to-lab variation and to ensure repeatability and reproducibility [22,23,24]. Regardless of the collection method, attempts should be made to avoid repeated temperature fluctuation and to shorten the transportation time [25]. Results on the short-term storage of samples with alternative methods do not necessarily translate into unbiased microbiota results in long-term storage, which is often the case in cohort biobanks. In two small pilots, microbiome profile remained consistent as late as 10 years following storage on FOBT card [26, 27] suggesting the suitability of FOBT card (or similar products) for long-term cohort biobanking.


In conclusion, acknowledging the small sample size in this pilot study, we provide evidence that microbiota profile of fecal samples collected in ethanol or on card and delayed storage at − 80 °C for up to 2 days is comparable to immediately frozen samples. Therefore, sample collection on card could potentially provide a practical alternative method of fecal biobanking in large-scale population-based cohorts. Although, the suitability of samples collected on card for microbiome studies other than 16S rRNA gene sequencing has not been confirmed, the cost- and space-effectiveness of this approach would be of importance to large-scale cohorts significantly facilitating the translational aspect of the microbiome research.


Study design

We performed a pilot study to assess the impact of different storage methods on fecal microbiota composition. Eight healthy volunteers provided self-collected fecal samples. To minimise variability, we recruited female research staff from Digestive Disease Research Institute with an age range of 25–40 years. Exclusion criteria included pregnancy or lactation at the time of sample collection, history of major gastrointestinal diseases including colorectal cancer or inflammatory bowel disease, and major gastrointestinal surgery. The study was approved by the Institutional Review Board of the Digestive Disease Research institute. Written informed consents were obtained from participants.

Sample collection

Individuals were instructed to collect the fecal samples using sterile fecal collection device (Pole Ideal Pars co., Iran). The samples were immediately transferred to the lab, where they were aliquoted. One aliquot was immediately stored at − 80 °C without preservative. Using a sterile swab, fecal material was transferred to tubes containing 95% ethanol and onto blood collection cards (Kawsar Biotech Company, Iran). Fecal samples in ethanol and on card were kept at room temperature (RT) for 48 h prior to storing at − 80 °C. The samples were stored for two months before DNA extraction.

Library construction and 16S rRNA gene sequencing

DNA extraction was performed with the QIAamp DNA Stool Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instruction. The V4 variable region of the microbial 16S rRNA gene was amplified with modified 515F/806R primers [28], and prepared for next-generation sequencing using a modified two-step targeted amplicon sequencing approach, as described previously [29]. Sequencing was performed using an Illumina MiSeq (Illumina, San Diego, CA, USA), with a V2 kit and paired-end 250 base reads at the Sequencing Core at the University of Illinois at Chicago. Raw fastq files for each sample were processed using the software package PEAR (Paired-end read merger) (v0.9.8) [30]. The merged fastq files were imported into the software package CLCGenomics Workbench v10.0 (CLC Bio, Aarhus, Denmark, Qiagen, Venlo, The Netherlands). Sequences without both the forward and reverse primers were discarded and subsequently, the primer sequences were removed. Sequences with low quality base calling scores (< 20) and insufficient length (< 250 bp) were discarded. Chimeric reads were filtered using USEARCH (v6.1) [31] within the open-source software QIIME (v1.8) environment [32]. Reads were clustered into OTUs using denovo OTU picking based on 97% similarity using UCLUST [31]. Representative sequences for each OTU were selected, and these sequences were assigned taxonomy using UCLUST and aligned to the 2013 release of the Greengenes reference database at 97% sequence similarity [33]. Demultiplexed sequencing data was deposited into the Sequence Read Archive of NCBI and can be accessed via accession number SRR8979312.

Microbial data pre-processing

Data analysis was conducted in R [34]. Initial preprocessing of the OTU table was conducted using the Phyloseq package [35]. Data was rarefied to the minimum 16,304 sequencing reads per sample. The numbers of sequencing reads of taxa were then relativized to the total sum of 16,304 without excluding rare taxa. Overall, 30,244 unique OTUs were detected. OTUs belonging to Archea, phylum Cyanobacteria, family of mitochondria, and class of chloroplast were excluded from the analysis resulting in 27,263 remaining OTUs.

Statistical analysis

The shared taxa were visualized using VennDiagram [36]. Alpha (α) diversity was assessed on the rarefied OTU table by the observed OTUs (richness) and inverse Simpson index (diversity). Association of α diversity with the storage method was assessed by one-way analysis of variance (ANOVA). Between-sample dissimilarity (β diversity) was assessed on Bray-Curtis dissimilarity and Jaccard distance indices for all and rare OTUs (< 0.01% mean relative abundance) by permutational ANOVA (PERMANOVA) using the vegan package [37]. To control for the compositional nature of the data, OTU counts were centre log-ratio transformed following zero-replacement [38, 39]. After this transformation, taxa relative abundances were compared at phylum and genus levels by ANOVA. Significant overall associations were further evaluated by post hoc pairwise comparison. The P values were corrected with Benjamini-Hochberg’s false discovery rate (FDR) method [40]. Corrected FDR p values of less than 0.05 were considered significant.

Availability of data and materials

The datasets generated and analysed during the current study are available in Sequence Read Archive of NCBI and can be accessed via accession number SRR8979312.



One-way analysis of variance


False discovery rate


Fecal occult blood test


Operational taxonomic unit


Permutational ANOVA


Room temperature


  1. 1.

    Marchesi JR, Adams DH, Fava F, Hermes GD, Hirschfield GM, Hold G, Quraishi MN, Kinross J, Smidt H, Tuohy KM, et al. The gut microbiota and host health: a new clinical frontier. Gut. 2016;65(2):330–9.

    PubMed  Article  Google Scholar 

  2. 2.

    Zeevi D, Korem T, Zmora N, Israeli D, Rothschild D, Weinberger A, Ben-Yacov O, Lador D, Avnit-Sagi T, Lotan-Pompan M, et al. Personalized nutrition by prediction of glycemic responses. Cell. 2015;163(5):1079–94.

    CAS  PubMed  Article  Google Scholar 

  3. 3.

    Lloyd-Price J, Abu-Ali G, Huttenhower C. The healthy human microbiome. Genome Med. 2016;8(1):51.

    PubMed  PubMed Central  Article  Google Scholar 

  4. 4.

    Backhed F, Fraser CM, Ringel Y, Sanders ME, Sartor RB, Sherman PM, Versalovic J, Young V, Finlay BB. Defining a healthy human gut microbiome: current concepts, future directions, and clinical applications. Cell Host Microbe. 2012;12(5):611–22.

    PubMed  Article  CAS  Google Scholar 

  5. 5.

    Mallick H, Ma S, Franzosa EA, Vatanen T, Morgan XC, Huttenhower C. Experimental design and quantitative analysis of microbial community multiomics. Genome Biol. 2017;18(1):228.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  6. 6.

    Fouhy F, Deane J, Rea MC, O'Sullivan O, Ross RP, O'Callaghan G, Plant BJ, Stanton C. The effects of freezing on faecal microbiota as determined using MiSeq sequencing and culture-based investigations. PLoS One. 2015;10(3):e0119355.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  7. 7.

    Vandeputte D, Tito RY, Vanleeuwen R, Falony G, Raes J. Practical considerations for large-scale gut microbiome studies. FEMS Microbiol Rev. 2017;41(Supp_1):S154–67.

    PubMed  Google Scholar 

  8. 8.

    Wu GD, Lewis JD, Hoffmann C, Chen Y-Y, Knight R, Bittinger K, Hwang J, Chen J, Berkowsky R, Nessel L. Sampling and pyrosequencing methods for characterizing bacterial communities in the human gut using 16S sequence tags. BMC Microbiol. 2010;10(1):1.

    CAS  Article  Google Scholar 

  9. 9.

    Han M, Hao L, Lin Y, Li F, Wang J, Yang H, Xiao L, Kristiansen K, Jia H, Li J. A novel affordable reagent for room temperature storage and transport of fecal samples for metagenomic analyses. Microbiome. 2018;6(1):43.

    PubMed  PubMed Central  Article  Google Scholar 

  10. 10.

    Adler CJ, White A, Bockmann M, Browne GV, Townsend G, Hughes T. VMG II transport medium stabilises oral microbiome samples for next-generation sequencing. J Microbiol Methods. 2018;144:91–8.

    CAS  PubMed  Article  Google Scholar 

  11. 11.

    Hao L, Xia Z, Yang H, Wang J, Han M. Ionic liquid-based reagents improve the stability of midterm fecal sample storage. J Microbiol Methods. 2017;139:68–73.

    PubMed  Article  Google Scholar 

  12. 12.

    Wong WSW, Clemency N, Klein E, Provenzano M, Iyer R, Niederhuber JE, Hourigan SK. Collection of non-meconium stool on fecal occult blood cards is an effective method for fecal microbiota studies in infants. Microbiome. 2017;5(1):114.

    PubMed  PubMed Central  Article  Google Scholar 

  13. 13.

    Franzosa EA, Morgan XC, Segata N, Waldron L, Reyes J, Earl AM, Giannoukos G, Boylan MR, Ciulla D, Gevers D, et al. Relating the metatranscriptome and metagenome of the human gut. Proc Natl Acad Sci U S A. 2014;111(22):E2329–38.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  14. 14.

    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).

  15. 15.

    Penington JS, Penno MAS, Ngui KM, Ajami NJ, Roth-Schulze AJ, Wilcox SA, Bandala-Sanchez E, Wentworth JM, Barry SC, Brown CY, et al. Influence of fecal collection conditions and 16S rRNA gene sequencing at two centers on human gut microbiota analysis. Sci Rep. 2018;8(1):4386.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  16. 16.

    Vogtmann E, Chen J, Kibriya MG, Chen Y, Islam T, Eunes M, Ahmed A, Naher J, Rahman A, Amir A, et al. Comparison of fecal collection methods for microbiota studies in Bangladesh. Appl Environ Microbiol. 2017;83(10):pii: e00361–17.

    Article  Google Scholar 

  17. 17.

    Sinha R, Chen J, Amir A, Vogtmann E, Shi J, Inman KS, Flores R, Sampson J, Knight R, Chia N. Collecting fecal samples for microbiome analyses in epidemiology studies. Cancer Epidemiol Biomark Prev. 2016;25(2):407–16.

    Article  Google Scholar 

  18. 18.

    Sorlie P, Wei GS. Population-based cohort studies: still relevant? J Am Coll Cardiol. 2011;58(19):2010–3.

    PubMed  Article  Google Scholar 

  19. 19.

    Poustchi H, Eghtesad S, Kamangar F, Etemadi A, Keshtkar AA, Hekmatdoost A, Mohammadi Z, Mahmoudi Z, Shayanrad A, Roozafzai F, et al. Prospective epidemiological research studies in IrAN (the PERSIAN cohort): rationale, objectives and design. Am J Epidemiol. 2018;187(4):647–55.

    PubMed  Article  Google Scholar 

  20. 20.

    Morgan XC, Huttenhower C. Chapter 12: human microbiome analysis. PLoS Comput Biol. 2012;8(12):e1002808.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  21. 21.

    Loftfield E, Vogtmann E, Sampson JN, Moore SC, Nelson H, Knight R, Chia N, Sinha R. Comparison of collection methods for fecal samples for discovery metabolomics in epidemiologic studies. Cancer Epidemiol Biomark Prev. 2016;25(11):1483–90.

    Article  Google Scholar 

  22. 22.

    Tabb DL, Vega-Montoto L, Rudnick PA, Variyath AM, Ham AJ, Bunk DM, Kilpatrick LE, Billheimer DD, Blackman RK, Cardasis HL, et al. Repeatability and reproducibility in proteomic identifications by liquid chromatography-tandem mass spectrometry. J Proteome Res. 2010;9(2):761–76.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  23. 23.

    Sinha R, Abnet CC, White O, Knight R, Huttenhower C. The microbiome quality control project: baseline study design and future directions. Genome Biol. 2015;16:276.

    PubMed  PubMed Central  Article  Google Scholar 

  24. 24.

    Sinha R, Abu-Ali G, Vogtmann E, Fodor AA, Ren B, Amir A, Schwager E, Crabtree J, Ma S. Microbiome quality control project C et al: assessment of variation in microbial community amplicon sequencing by the microbiome quality control (MBQC) project consortium. Nat Biotechnol. 2017;35(11):1077–86.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  25. 25.

    Wu WK, Chen CC, Panyod S, Chen RA, Wu MS, Sheen LY, Chang SC. Optimization of fecal sample processing for microbiome study - the journey from bathroom to bench. J Formos Med Assoc. 2019;118(2):545-55.

    PubMed  Article  Google Scholar 

  26. 26.

    Taylor M, Wood HM, Halloran SP, Quirke P. Examining the potential use and long-term stability of guaiac faecal occult blood test cards for microbial DNA 16S rRNA sequencing. J Clin Pathol. 2017;70(7):600–6.

    CAS  PubMed  Article  Google Scholar 

  27. 27.

    Couto Furtado Albuquerque M, van Herwaarden Y, Kortman GAM, Dutilh BE, Bisseling T, Boleij A. Preservation of bacterial DNA in 10-year-old guaiac FOBT cards and FIT tubes. J Clin Pathol. 2017;70(11):994–6.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  28. 28.

    Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Huntley J, Fierer N, Owens SM, Betley J, Fraser L, Bauer M, et al. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012;6(8):1621–4.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  29. 29.

    Naqib A, Poggi S, Wang W, Hyde M, Kunstman K, Green SJ. Making and sequencing heavily multiplexed, high-throughput 16S ribosomal RNA gene amplicon libraries using a flexible, two-stage PCR protocol. In: Gene expression analysis. New York, NY: Humana Press; 2018. p. 149–69.

    Google Scholar 

  30. 30.

    Zhang J, Kobert K, Flouri T, Stamatakis A. PEAR: a fast and accurate Illumina paired-end reAd mergeR. Bioinformatics. 2014;30(5):614–20.

    CAS  PubMed  Article  Google Scholar 

  31. 31.

    Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26(19):2460–1.

    CAS  Article  Google Scholar 

  32. 32.

    Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Pena AG, Goodrich JK, Gordon JI, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7(5):335–6.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  33. 33.

    DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, Huber T, Dalevi D, Hu P, Andersen GL. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006;72(7):5069–72.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  34. 34.

    R Core Team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2017.

    Google Scholar 

  35. 35.

    McMurdie PJ, Holmes S. Phyloseq: an R package for reproducible intekractive analysis and graphics of microbiome census data. PLoS One. 2013;8(4):e61217.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  36. 36.

    Chen H. VennDiagram: Generate High-Resolution Venn and Euler Plots. R package version 1.6.17. 2016.

  37. 37.

    Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, Minchin PR, O'Hara RB, Simpson GL, Solymos P, et al. Vegan: community ecology package. R package version 2.4–3. 2017

    Google Scholar 

  38. 38.

    Palarea-Albaladejo J, Martin-Fernandez JA. zCompositions -- R package for multivariate imputation of left-censored data under a compositional approach. Chemom Intell Lab Syst. 2015;143:85–96.

    CAS  Article  Google Scholar 

  39. 39.

    Gloor GB, Reid G. Compositional analysis: a valid approach to analyze microbiome high-throughput sequencing data. Can J Microbiol. 2016;62(8):692–703.

    CAS  PubMed  Article  Google Scholar 

  40. 40.

    Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Statist Soc B. 1995;57(1):289–300.

    Google Scholar 

Download references


Not applicable.


This study was funded by an internal grant from Digestive Disease Research Institute. The funding body did not have a role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information




S Moossavi, AK, RM conceived the study; S Moossavi, PAE, RG, SJG, AN performed the experiments; S Moossavi performed the statistical analysis; S Moossavi and RM wrote the initial draft of the manuscript; PAE, RG, SJG, AN, FB, S Merat, HP, AK revised and edited the manuscript; RM finalised the manuscript; all authors read and approved the final manuscript.

Corresponding author

Correspondence to Reza Malekzadeh.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Institutional Review Board of the Digestive Disease Research institute. Written informed consents were obtained from participants.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional file

Additional file 1:

Table S1. Percent difference in the relative abundances of fecal genera (> 1%) between immediate freezing at − 80 °C (A) and storage in 95% ethanol for 48 h at room temperature (B) in eight healthy volunteers. Table S2. Percent difference in the relative abundances of fecal genera (> 1%) between immediate freezing at − 80 °C (A) and storage on card stored for 48 h at room temperature (C) in eight healthy volunteers. Figure S1. Effect of storage condition on the fecal microbiota profile at species level for abundant species defined as having mean relative abundance of > 0.01% within each method. Figure S2. Comparison of relative abundances of members of Actinobacteria at family level. Figure S3. Comparison of relative abundances of gut microbiota profile at family level across individuals and methods. (PDF 1343 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Moossavi, S., Engen, P.A., Ghanbari, R. et al. Assessment of the impact of different fecal storage protocols on the microbiota diversity and composition: a pilot study. BMC Microbiol 19, 145 (2019).

Download citation


  • Gut microbiota
  • Cohort study
  • Fecal storage
  • Card
  • Mock shipment