Skip to main content
  • Research article
  • Open access
  • Published:

Evaluating gut microbiota profiles from archived fecal samples

Abstract

Background

Associations between colorectal cancer and microbiota have been identified. Archived fecal samples might be valuable sample sources for investigating causality in carcinogenesis and biomarkers discovery due to the potential of performing longitudinal studies. However, the quality, quantity and stability of the gut microbiota in these fecal samples must be assessed prior to such studies. We evaluated i) cross-contamination during analysis for fecal blood and ii) evaporation in stored perforated fecal immunochemical tests (iFOBT) samples, iii) temperature stability as well as iv) comparison of the gut microbiota diversity and composition in archived, iFOBT and fresh fecal samples in order to assess feasibility of large scale microbiota studies.

Methods

The microbiota profiles were obtained by sequencing the V3-V4 region of 16S rDNA gene.

Results

The iFOBT does not introduce any cross-sample contamination detectable by qPCR. Neither could we detect evaporation during freeze-thaw cycle of perforated iFOBT samples. Our results confirm room temperature stability of the gut microbiome. Diverse microbial profiles were achieved in 100% of fresh, 81% of long-term archived and 96% of iFOBT samples. Microbial diversity and composition were comparable between fresh and iFOBT samples, however, diversity differed significantly between long-term archived, fresh and iFOBT samples.

Conclusion

Our data showed that it is feasible to exploit archived fecal sample sets originally collected for testing of fecal blood. The advantages of using these sample sets for microbial biomarker discovery and longitudinal observational studies are the availability of high-quality diagnostic and follow-up data. However, care must be taken when microbiota are profiled in long-term archived fecal samples.

Peer Review reports

Background

Colorectal cancer (CRC) is the third most frequently diagnosed cancer with more than 1.4 million new cases diagnosed annually worldwide [1, 2]. The World Health Organization and the European Union recommend CRC screening to detect CRC and its precursors at a curable stage. However, the most frequently used non-invasive fecal tests (guaiac-based fecal occult blood test (gFOBT) or an immunochemical FOBT (iFOBT) have inadequate sensitivity [1, 3]. A multi target stool DNA test combined with iFOBT has shown high sensitivity, but poor specificity [3]. Highly specific and sensitive non-invasive screening tests are urgently needed. Thus, identifying fecal biomarkers for CRC risk assessment, early detection and prognosis is a priority.

Cross-sectional human studies have shown associations between bacteria and colorectal cancer [4, 5] and the microbiota differs between neoplastic lesions and healthy mucosa [6, 7]. The gut microbiome is less studied in the early stages of carcinogenesis. Hale et al. [8], showed that the bile-tolerant microbes Bilophila, Desulfovibrio, proinflammatory bacteria in the genus Mogibacterium, and multiple Bacteroidetes species are more abundant in adenoma cases. A microbiome biomarker based on 16S sequencing combined with FOBT screening was shown to increase CRC detection compared to FOBT screening [9]. However, the studies lack co-variables and could only marginally improve detection of adenomas [9, 10]. It is probable that, with further development, the gut microbiome can be used to detect the presence of precancerous and cancerous lesions [11]. However, human studies have been limited to cross-sectional case-control series showing bacterial associations to CRC rather than a causative link to future risk which requires a prospective, observational study design [12, 13].

Fecal samples from CRC screening can provide large number of samples with diagnostic data needed for biomarker discovery in longitudinal studies [14]. Long-term stored biobanks, collected before preservatives were common, may have follow-up data from long-term projects or national registries, such as the NORCCAP cohort [15, 16]. Biobanks of FOBT samples are often large and collected at all stages of CRC development, such as the The Bowel Cancer Screening in Norway (BCSN trial) [17]. These biobanks are invaluable for assessing the temporal dynamics of the microbiota-CRC relationship from years before diagnosis to times after treatment.

Nonetheless, using fecal samples not originally intended for microbial profiling may introduce technical challenges due to inadequate material and varying sample handling and storage. Technical aspects regarding storage conditions, freeze-thaw effects and use of preservative media have previously been studied in fresh fecal samples [14, 18,19,20,21]. Good microbiota concordance has been shown between fresh stool and FOBT cards [22]. However, all aspects of how sample processing and storage may potentially influence the microbiota need to be investigated.

The aims of this study were therefore to investigate the feasibility of using iFOBT and long-term archived fecal samples frozen for up to 16 years for 16S rDNA sequencing. Specifically, cross-contamination during the test for fecal blood, reduction of sample quality due to storage and handling of perforated iFOBT samples, temperature stability and comparison of the gut microbiota diversity and composition was evaluated.

Methods

Samples and design

The BCSN trial is a comparative effectiveness research trial randomizing 140.000 persons aged 50–74 years to be screened for CRC either with flexible sigmoidoscopy (FS) or iFOBT [17]. About 70.000 individuals have been invited for iFOBT screening and will deliver iFOBT samples every other year for 10 years. The sampling kits for iFOBT are mailed to the participants and collected by return mail with a shipping time of 2–6 days. The stool samples are collected on plastic sticks designed to catch about 10 μg stool and then stored in 2 ml buffer containing HEPES buffer (4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid), BSA (Bovine serum albumin) and sodium azide. The participants were asked to mail the sample to the lab as soon as possible. On the day of arrival to the laboratory, the samples were analyzed with an immunochemical test for human blood (globin) (Eiken Chemicals Ltd., Tokyo, Japan) and the perforated iFOBT tube was frozen at − 80 °C immediately thereafter. We randomly selected 50 anonymized iFOBT samples among the participants in the BCSN trial for this study. The samples have not previously been thawed.

The NORCCAP study was a randomized flexible sigmoidoscopy screening clinical trial of 100,210 individuals aged 50 to 64 years from the populations in two Norwegian counties (Oslo and Telemark). The trial was performed during 1999–2000 (55–64–year age group) and in 2001 (50–54–year age group) [15, 16]. The intervention group received flexible sigmoidoscopy and delivered stool samples for iFOBT. To-date the trial has about 4800 stool samples stored at − 30 °C. 50 feces samples randomly selected from NORCCAP participants with normal sigmoidoscopy were included in this study. Stools were sampled in 20 ml vials by the screening participants at home and kept for 1–7 days in their home freeze (− 20 °C) before delivery to the screening centre when attending FS screening. Samples were transported in the 20 ml vials with an outer transport tube. Most participants had less than 30 min travel to the screening centre, but partial thawing during this transport is likely. At the screening centers, stools were stored at − 20 °C for up to 1 week before transport for up to 2 h in insulated transport boxes to a central − 30 °C storage room with alarm systems. Some samples may have been partially thawed to provide material for non-microbial research purposes.

Fresh samples from eight anonymous, presumably healthy individuals, four men and four women, aged 50 to 64 years were collected in 2015. Each individual froze one part of their stool sample directly and one part was stored at room temperature for 48 h before it was frozen at home. Partial thawing may have occurred during transport to the lab. All samples were frozen at − 20 °C in the lab until DNA extraction.

The collection of iFOBT and archived samples was approved by the Regional Committees for Medical and Health Research Ethics in South-Eastern Norway (2011/1272 and 2010/3087 A, respectively). An additional ethical approval for the feasibility study was not required since all samples were anonymized and the purpose of the study was method related, as stated by the Norwegian regional ethics committee (Regional Committees for Medical and Health Research Ethics in South-Eastern Norway, ref. 2015/9).

Cross-contamination during iFOBT analysis

In the automatized analysis of iFOBT, the same instrument needle extracts liquid from all samples in a consecutive order with intermediary washing. To investigate if this may result in cross-contamination of the samples, a test sample series alternating between 200 ng/μl human DNA (Roche #11691112001) and water was prepared. This sample series was run on the desktop instrument OC-Sensor Diana (Eiken Chemicals,Tokyo, Japan). A qPCR measurement of human DNA beta-globin in the 40 water samples and two positive controls were analyzed in two replicates in 25 μl total TaqMan universal PCR master mix (Thermofisher Scientific, Waltham, USA) reactions using 5 μl template [23].

Sample handling and storage effects on perforated iFOBT samples

The iFOBT tubes are perforated during iFOBT analysis. Evaporation during the freezing process, sublimation during storage or condensation during thawing of these “open” tubes might increase or decrease sample volume indicating that contamination of samples due to these processes might be possible. To asses if samples quality could be diminished by such processes we tested for changes in sample volume by freezing 50 perforated iFOBT tubes with and without parafilm and parafilm with rubber bands. We compared the weight of the 50 tubes before and after 6 months in the freezer (− 80 °C), 3 h after thawing.

DNA isolation, PCR and sequencing

All samples were processed in Lysing Matrix E tubes with silica beads. Sample input was about 10 μl, approximately 10 mg, for the solid archived fecal samples and 500 μl for the iFOBT samples. Dry weight of feces varies in the iFOBT samples, however, we estimate that about 2.5–5 μg feces were used in the protocol. PBS (phosphate-buffered saline) buffer was added to a total liquid volume of 1.1 ml prior to mechanical lysis on a FastPrep 24 and centrifugation at 400×g for 1 min. Nucleic acids were extracted on an automated platform, QIAsymphony, using a standard nucleic acid extraction protocol that included pretreatment with ProteinaseK and a final elution in 85 μl. To account for the level of microbiological contamination from the environment and reagents an extraction-negative control (ENC) containing all the reagents except for sample material was processed in each batch of samples. Two additional negative control samples were included in the PCR, one with water and one with 10 ng/μl human DNA. Due to experimentally encountered inhibition of PCR in the archived samples, these were diluted 1:5 before used as template for PCR. All other samples were used directly in PCR without dilution.

A two-step PCR protocol was used. The primers S-D-Bact-0341-b-S-17 and S-D-Bact-0785-a-A-21 [24] were first used to amplify a 16S rDNA in the V3-V4 region, to produce a 464 nucleotide (nt) long amplicon. PCR was performed in 20 μl reaction volumes with 1× Phusion Master Mix, 0.25 μM each 16S primer and 2 μl template under the following conditions; initial denaturation 98 °C for 30 s and then 27 cycles of 98 °C for 10 s, 62 °C for 30 s and 72 °C for 15 s, before a final extension at 72 °C for 10 min and storage at 10 °C. PCR cleanup was done using partly modified PERFORMA DTR V3 96-well Short Plates and QuickStep™2, 96 Well PCR Purification kit protocols. All samples PCR products were tested on 1% agarose gels using 5 μl sample. Then in the second round of PCR, indexes were added. The index PCR was also performed in 20 μl volumes using 1× Phusion Master Mix, 0.5 μM each index primer, and 1 μl of template DNA from the previous PCR, diluted 1:100.

Sequencing of the amplicons was performed on an Illumina MiSeq (Illumina) desktop sequencer using V3 chemistry 2 × 300 cycle kits. This protocol has shown high reproducibility and repeatability and very little contamination [25].

Bioinformatics and statistics

Sequencing quality filtering and processing were carried out using the MiSeq SOP and the Mothur pipeline (v.1.36.1) [26, 27]. To ensure high-quality data for analysis, sequence reads containing ambiguous bases, homopolymers > 8 bp, more than one mismatch in the primer sequence, or sequences under the default per base quality score were removed. Assembled reads > 460 bp in length, singletons and chimeric sequences identified with the UCHIME algorithm [27] were excluded from the analyses. The high-quality assembled sequences (contigs) were aligned to the Silva 16S rRNA database (v.119) [28] and clustered into operational taxonomic units (OTUs) at a cut-off value > 98% in a closed reference OTU approach. OTUs were calculated at distance 0.02, and alpha diversity (Shannon index and inverse Simpson index) was calculated per sample. The number of assembled sequences for each sample was rarefied to 2000 to minimize the impacts of uneven sampling. Inverse Simpson’s index is weighted on dominant species, whereas the Shannon index assumes all species are represented [29]. To identify differences in species diversity between the groups, a two-way ANOVA and post hoc comparisons (Tukey Honest test of Significant Differences - TukeyHSD) of the observed number of OTUs, Inverse Simpson diversity and the Shannon index were used. Beta diversity (i.e., the variation in community composition between microbiota samples) was calculated between the sample types using the Bray–Curtis dissimilarity index [30]. Dissimilarities between sample groups were tested using permutational analysis of variance (PERMANOVA) implemented in the vegan (v. 2.5–2) R package and the “adonis” function. Differentially abundant OTUs were identified using the zero-inflated Gaussian model (fitZig), implemented in the metagenomeSeq package [31, 32]. The community composition data was transformed using the “decostand” function and divided by the margin total. Principal components analyses were done using the prcomp function in R. Heatmaps of the log transformed OTUs (for OTUs with more than 10 in counts) were produced using hierarchical clustering and Euclidean distance. To improve visualization only OTUs with a log sum > 20 for all samples, and > 2 for the fresh samples were shown. The clustering of these heatmaps were indistinguishable from heatmap with all OTUs included.

All of the ENC and negative controls, but also four iFOBTs and 14 archived samples did not produce sufficient aligned reads to calculate diversity.

Results

No cross-contamination during iFOBT or change in volume in perforated samples during storage

Potential transfer of DNA from sample to sample during iFOBT analysis was tested with qPCR. 79 of the 80 qPCR measurements were negative. i.e. no contamination. One sample had a cycle threshold (Ct) value of 40.7, also considered a negative result. The positive controls were positive on the qPCR (Ct values < 30).

The weight difference of 50 perforated iFOBT tubes before freezing (383.6 g) and after 6 months in the freezer (383.7 g) was 0.1 g or 0.03%. The weight difference of 50 perforated iFOBT tubes packed in parafilm to prevent change in volume before freezing (384,5 g) and after freezing (384,9 g) was 0.4 g or 0.1%. Some condensation on the tube surface was observed.

Microbiota differs in iFOBT, fresh and archived samples

In total, 4.2 million 16S rDNA contigs were assembled from the sequencing data. The fresh samples were sequenced separately, with a larger fraction of the sequencing pool per sample than the other sample groups, thus producing more reads, contigs and OTUs (Fig. 1a, Additional files 1, 2 and 3). Archived samples produced less OTUs than fresh and iFOBT samples (Table 1). Gel electrophoresis of the 16S rDNA PCR product from the archived fecal samples showed increased PCR efficiency when the DNA solutions were diluted 5-fold compared to undiluted DNA, signifying PCR inhibition. There were no signs of PCR inhibition in iFOBT samples.

Fig. 1
figure 1

Alpha and beta diversity in iFOBT, fresh and archived samples. a Shows a boxplot of the number of observed OTUs in each sample group. b and c Shows boxplots of the Inverse Simpson index (b) and Shannon (c) index in fecal immunochemical tests (iFOBT samples, fresh fecal samples and fecal samples archived for approximately 16 years. The indexes are based on rarefied OTU data to minimize the impacts of uneven sampling. The Bray-Curtis dissimilarity index for comparisons of groups are shown (d)

Table 1 Overview of operational taxonomical unit (OTU) counts in each sample type

Diversity

The mean, maximum and minimum numbers of contigs (assembled reads) assigned to OTUs per sample group are shown in Table 1. The ENC and negative control samples produced very few OTUs. 96% of the iFOBT, 81% of the archived and 100% of the fresh samples passed a minimum relative abundance of 200 OTUs. The minimum relative OTU abundance chosen as a threshold was based on the sequencing depth and distribution of OTU counts (Additional file 1).

We compared the diversity in iFOBT, fresh and archived samples using observed OTUs, and the Inverse Simpson and Shannon indexes from rarefied OTU tables (Fig. 1a, b, c and Additional file 4). Fresh and iFOBT samples were similar in mean diversity while the archived samples had statistically significant less diversity compared to iFOBT (P-value = 0.027 for Inverse Simpson and P-value = 0.06 for Shannon) and fresh (P-value < 0.001 for both Inverse Simpson and Shannon) in the ANOVA and TukeyHSD test. The Bray-Curtis dissimilarity index showed that the microbiota in iFOBT, archived and fresh samples differ (Fig. 1d and Additional file 5) and the PERMANOVA analysis showed that this difference is significant (P-value < 0.001, F = 1.4587). The first and second principal components showed overlapping community composition (Fig. 2).

Fig. 2
figure 2

Principal component plot of the community composition of iFOBT, archived and fresh sample. The first and second principal components of the community composition of iFOBT samples from a screening trial in Norway (BCSN), archived samples from the NORCCAP cohort stored for about 16 years and fresh samples

Relative abundance

The microbiota profiles differed between fresh and iFOBT vs. archived fecal samples (Fig. 3). More than 62% of the archived fecal samples clustered together, while the fresh and iFOBT samples were intermixed (Fig. 3). In particular, the OTUs 1 to 6 (rightmost vertical clade), including Lachnospiraceae and Ruminococcaceae family OTUs, were of high abundance in fresh and iFOBT samples, whereas OTU 7 and 11, representing the Peptostreptococcaceae and Ruminococcaceae family, were abundant in archived samples.

Fig. 3
figure 3

Clustering of archived, fresh and fecal immunochemical tests (iFOBT) fecal samples. Heatmap of the log transformed OTU table for all samples produced by hierarchical clustering and Euclidean distance. Only OTUs with a log sum > 20 were illustrated. iFOBT samples are marked in green, archived fecal samples stored for 14 to 16 years at − 30 °C are marked in blue and fresh samples are marked in pink. All samples are from presumably healthy individuals

The zero-inflated log-normal mixture model identified 11 OTUs that significantly (P-value < 0.05) differed between archived and iFOBT samples (Table 2 and Additional file 6). The OTUs of the Ruminococcaceae family (OTU4, 14 and 9) and OTU5 Lachnospiraceae were significantly less abundant in archived samples. Peptostreptococcaceae (OTU11) and Clostridiaceae (OTU18) were more abundant in archived fecal samples.

Table 2 Differences in OTU abundance in archived and iFOBT samples

Stability of fecal microbial composition in room temperature

The mean alpha-diversity did not differ significantly between samples that had been directly frozen and samples that were frozen after 48 h in room temperature (Observed OTU count P-value = 0.18, Inverse Simpson P-value = 0.65 and Shannon P-value = 0.94) (Fig. 4a and Additional file 4). Dissimilarities (Bray-Curtis) between sample pairs from the same individual with different room temperature storage time, was lower than between samples from different individuals (Fig. 4b). Unsupervised clustering of the log-transformed OTU relative abundance supports similarities between sample pairs (Fig. 4c). The log-normal mixture model identified no significant differences in relative abundance (Additional file 6) between samples stored in room temperature for 48 h and samples directly frozen (Additional file 7).

Fig. 4
figure 4

Microbiota profiles in fresh samples frozen directly and frozen after 48 h in room temperature. a Boxplot of Inverse Simpson alpha-diversity index in samples frozen directly and frozen after 48 h. b Bray-Curtis dissimilarities between samples from the same individuals with different storage conditions, and between different individuals regardless of room temperature storage. c Heatmap of the log transformed OTU table from paired fecal samples from 8 presumably healthy individuals (serial number 01 to 10). One part of the samples was directly frozen (0 h, marked with green text) and the other part of the samples was frozen after 48 h in room temperature (48 h, marked with blue text). Hierarchical clustering and Euclidean distance produced the clustering showing higher inter-person variability than intra-person variability. Only OTUs with a log sum > 2 were illustrated

Discussion

Fecal samples collected as part of CRC screening are usually intended for analyses of occult blood and not for studying the microbiota. Thus, it is necessary to assess the microbial DNA quantity, quality and biases introduced by sample collection and storage before these samples can be used in microbial epidemiological studies.

Our study has shown that it is feasible to profile the microbiota directly from iFOBT samples, also shown by others [9, 33]. Fecal samples archived for up to 16 years, collected as part of the NORCCAP study [17, 22], also produced microbial profiles sufficient for identification of inter-individual differences and association studies, although further improvement of the protocol such as homogenization of the samples [34], is required.

Microbial profiling with highly sensitive methods such as Next-Generation Sequencing will likely detect potential cross-sample contamination that may have occurred during analysis for blood in the iFOBT samples. This type of cross-contamination would distort the microbial profiles and make the samples unsuitable for epidemiological studies. However, our experiment using qPCR tests on runs with alternating blanks and concentrated DNA added to the iFOBT samples could not detect any cross-sample contamination. Our simple freezing tests indicated that storage of perforated iFOBT tubes will not change the volume of the samples and therefore likely not reduce the microbiota quality. We therefore believe that the iFOBT samples are eligible for analysis using sensitive methodologies. These results are also relevant for other tests, such as DNA methylation detection combined with testing for occult blood from the same collection device [35].

Proficient microbial profiling depends on complete lysis, here ensured by the use of both mechanical and enzymatic steps during the DNA isolation, primers that amplify as wide a target as possible (evaluated by Klindworth et al. [24]) and effective amplification. Bile salts and complex polysaccharides in feces [36, 37] and heme in blood [38] are substances that inhibit the PCR and reduce the amplification efficiency. Reducing these components in the PCRs with fecal DNA extracts will likely increase OTU counts and potentially the diversity in the archived fecal samples.

Our data show differences in microbial composition between the fresh, iFOBT and archived fecal samples. Species diversity in fresh samples is more similar to the iFOBT samples than to the long-term archived samples. OTU count differences are strongly related to sequencing depth in a non-saturated setting. Taking this into account, the observed OTUs are comparable between iFOBT and archived samples. The community compositions of these sample groups are overlapping, although there are significant differences in beta-diversity. RNAlater® has been shown to preserve the microbial components in fecal samples [19]. It is likely that the iFOBT buffer containing HEPES, BSA and sodium azide may also preserve the microbiota. The sodium azide should stop any bacterial growth as it is a strong anti-microbial agent.

To our knowledge, there are no studies of the stability of the microbiota in long-term stored archived samples. We showed significant reduction of species diversity in samples archived for up to 16 years. However, long-term stored samples have successfully been used to identify associations in fecal microbiota between colorectal cancer cases and controls [33, 39]. Samples used in these studies were freeze-dried and the lyophilates were pooled, mixed, and stored at − 40 °C. Large cohorts with follow up data and fecal samples that are processed in this manner are very rare to come by. Reduced alpha and beta diversity estimates have previously been shown for samples collected in no buffer compared to iFOBT samples and attributed to the thawing process [19]. Storage in domestic freezers, fluctuating in temperatures, has been shown to shift the abundances of the major taxa [34]. Thus, freeze-thaw cycles or difference in storage time in frost-free domestic freezers regardless of the long-term storage time might explain the difference in species diversity.

Differences in bacterial relative abundance were also observed between iFOBT and archived samples in our study. Particularly the common gut Ruminococcaceae family within the Clostridia order was less abundant in archived samples. The relative abundance of Ruminococcus species has also been shown to decline with storage time in samples from different persons [21]. Therefore, a study with mixed sample types may detect mostly sample type differences and mask true associations. Although, on their own and with an improved protocol, archived samples are of high value, particularly due to the possibility of doing longitudinal studies. These observed differences emphasize the need to analyze the iFOBT and the archived NORCCAP materials separately in future studies on prevalent and incidents CRC during years of follow-up.

Several studies have investigated the stability of microbiota in fecal samples stored for different length of time in room temperature using different study design, sample conservation, limitations and endpoint, with conflicting results. Most studies have reported small or insignificant changes of the microbial composition or diversity for feces samples stored in preservative solutions [18, 20, 21, 33, 34, 38, 40,41,42]. A comprehensive test of preservation methods showed that variability in microbiota was only minimally explained by preservation medium and that storage effects are small compared to individual differences [42]. Fecal samples without preservative, comparable to our long-term archived samples, stored in freezers for 8 weeks, showed lower diversity than between individuals [42].

Flores et al. [19] have showed significant changes between replicates for rare OTUs (< 1% abundance), but otherwise stable microbiota for feces stored in RNAlater®. Beta diversity was stable for samples both with and without preservative. Contrary to our results, an increased variance for freezing delays was shown [19]. Two other studies have reported effect of storage on microbial composition and diversity [43, 44]. There is a need for longitudinal studies subjecting aliquots of the same samples to different storage for different length of time over years.

Ahn and co-workers [39] identified large differences in relative microbial abundance with reported odd ratios in the range of 4 to 5 between colorectal cancer cases and controls, reproduced by Vogtmann et al. [33]. More studies are needed in order to determine if the effects of delayed freezing and preservation media may distort the associative signals from for example colorectal cancer - control studies.

Although limited by samples stored in a buffer designed for iFOBT analyses, potential exposure to thawing during transport, use of home freezers and storage without preservative, our results indicate that samples collected during screening, not initially intended for the purpose, may be used for microbial profiling. The disadvantages of using screening samples are variable quality of material and insufficient control of sampling, storage and freeze-thaw cycles. However, we suggest that biases introduced from provenance of sample, sequencing protocol, primers and bioinformatics analyses may potentially pose a larger effect on the outcome than sample storage, media and freeze-thaw cycles, also proposed by Lauber et al. [20], and future meta-analyses may confirm this suggestion.

The advantages of using screening samples are large number of sample, population wide selections, well-characterized participants often with colonoscopy data and potentially long follow up times enabling longitudinal design. Some of the available samples and data sources are screening pilots and trials such as NORCCAP and the BCSN trial.

The established association between infections with certain microbes and colorectal cancer [1] can now be studied in depth with new technology and unprecedented accuracy and resolution. There is a need for longitudinal observational studies to define modes of microbial action and confirm a causative relationship for one of the world’s leading causes of cancer death. The results reported here may contribute to expand the use of archived fecal material for this purpose.

Conclusion

We have shown that it is feasible to use 1) fecal samples used for iFOBT screening and 2) fecal samples frozen for multiple years and potentially thawed several times, in large-scale epidemiological studies on microbiota. The test for fecal occult blood in the iFOBT sample does not introduce cross sample contamination and substantial storage effects, therefore the sample from the collection device may be used in highly sensitive methods such as next generation sequencing. The iFOBT samples are comparable in alpha diversity and composition to fresh fecal samples. However, storage and sampling methods have been shown to affect the microbiome composition and diversity, therefore in the case-control design, both cases and controls should be selected from the same sample group and processed in the same way. Care must be taken when using archived fecal samples and optimization of the protocol is critical.

Abbreviations

BCSN:

The Bowel Cancer Screening in Norway

BSA:

Bovine serum albumin

CRC:

colorectal cancer

ENC:

extraction-negative control

FIT:

fecal immunochemical tests (=iFOBT)

FS:

flexible sigmoidoscopy

gFOBT:

guaiac-based fecal occult blood test

HEPES:

4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid

iFOBT:

immunochemical-based fecal occult blood test

NORCCAP:

NORwegian Colorectal Cancer Prevention

OTU:

Operational taxonomic unit

PBS:

phosphate-buffered saline

References

  1. Schreuders EH, Ruco A, Rabeneck L, Schoen RE, Sung JJY, Young GP, Kuipers EJ. Colorectal cancer screening: a global overview of existing programmes. Gut. 2015;64(10):1637–49. https://doi.org/10.1136/gutjnl-2014-309086.

    Article  Google Scholar 

  2. Colorectal cancer estimated incidence, mortality and prevalence worldwide in 2012. [http://globocan.iarc.fr/old/FactSheets/cancers/colorectal-new.asp].

  3. Imperiale TF, Ransohoff DF, Itzkowitz SH, Levin TR, Lavin P, Lidgard GP, Ahlquist DA, Berger BM. Multitarget stool DNA testing for colorectal-cancer screening. N Engl J Med. 2014;370(14):1287–97.

    Article  CAS  Google Scholar 

  4. Bultman SJ. Emerging roles of the microbiome in cancer. Carcinogenesis. 2014;35(2):249–55.

    Article  CAS  Google Scholar 

  5. Garrett WS. Cancer and the microbiota. Science (New York, NY). 2015;348(6230):80–6.

    Article  CAS  Google Scholar 

  6. Kostic AD, Gevers D, Pedamallu CS, Michaud M, Duke F, Earl AM, Ojesina AI, Jung J, Bass AJ, Tabernero J, et al. Genomic analysis identifies association of fusobacterium with colorectal carcinoma. Genome Res. 2012;22(2):292–98. https://doi.org/10.1101/gr.126573.111.

    Article  Google Scholar 

  7. Castellarin M, Warren RL, Freeman JD, Dreolini L, Krzywinski M, Strauss J, Barnes R, Watson P, Allen-Vercoe E, Moore RA, et al. Fusobacterium nucleatum infection is prevalent in human colorectal carcinoma. Genome Res. 2012;22(2):299–306.

    Article  CAS  Google Scholar 

  8. Hale VL, Chen J, Johnson S, Harrington SC, Yab TC, Smyrk TC, Nelson H, Boardman LA, Druliner BR, Levin TR. Shifts in the fecal microbiota associated with adenomatous polyps. Cancer Epidemiol Biomark Prev. 2017;26(1):85–94.

    Article  CAS  Google Scholar 

  9. Baxter NT, Ruffin MT, Rogers MA, Schloss PD. Microbiota-based model improves the sensitivity of fecal immunochemical test for detecting colonic lesions. Genome Med. 2016;8(1):37.

    Article  Google Scholar 

  10. Zeller G, Tap J, Voigt AY, Sunagawa S, Kultima JR, Costea PI, Amiot A, Böhm J, Brunetti F, Habermann N. Potential of fecal microbiota for early-stage detection of colorectal cancer. Mol Syst Biol. 2014;10(11):766.

    Article  Google Scholar 

  11. Zackular JP, Rogers MA, MTt R, Schloss PD. The human gut microbiome as a screening tool for colorectal cancer. Cancer Prev Res (Philadelphia, Pa). 2014;7(11):1112–21.

    Article  CAS  Google Scholar 

  12. Gagniere J, Raisch J, Veziant J, Barnich N, Bonnet R, Buc E, Bringer MA, Pezet D, Bonnet M. Gut microbiota imbalance and colorectal cancer. World J Gastroenterol. 2016;22(2):501–18.

    Article  CAS  Google Scholar 

  13. Yu J, Feng Q, Wong SH, Zhang D, Liang QY, Qin Y, Tang L, Zhao H, Stenvang J, Li Y, et al. Metagenomic analysis of faecal microbiome as a tool towards targeted non-invasive biomarkers for colorectal cancer. Gut. 2017;66(1):70–78. https://doi.org/10.1136/gutjnl-2015-309800.

    Article  Google Scholar 

  14. 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 Biomarkers Prev. 2016;25(2):407–16.

    Article  Google Scholar 

  15. Gondal G, Grotmol T, Hofstad B, Bretthauer M, Eide TJ, Hoff G. The Norwegian colorectal Cancer prevention (NORCCAP) screening study: baseline findings and implementations for clinical work-up in age groups 50-64 years. Scand J Gastroenterol. 2003;38(6):635–42.

    Article  CAS  Google Scholar 

  16. Holme O, Loberg M, Kalager M, Bretthauer M, Hernan MA, Aas E, Eide TJ, Skovlund E, Schneede J, Tveit KM, et al. Effect of flexible sigmoidoscopy screening on colorectal cancer incidence and mortality: a randomized clinical trial. Jama. 2014;312(6):606–15.

    Article  CAS  Google Scholar 

  17. Lange T, Randel KR, Schult AL, Knudsen MD, Kirkoen B, Botteri E, Berstad P, Jorgensen A, Ursin G, Bretthauer M, et al. Sigmoidoscopy and faecal occult blood test - a comparative screening trial. Tidsskr Nor Laegeforen. 2017;137(10):727–30.

    Article  Google Scholar 

  18. Cardona S, Eck A, Cassellas M, Gallart M, Alastrue C, Dore J, Azpiroz F, Roca J, Guarner F, Manichanh C. Storage conditions of intestinal microbiota matter in metagenomic analysis. BMC Microbiol. 2012;12:158.

    Article  CAS  Google Scholar 

  19. Flores R, Shi J, Yu G, Ma B, Ravel J, Goedert JJ, Sinha R. Collection media and delayed freezing effects on microbial composition of human stool. Microbiome. 2015;12;3:33. https://doi.org/10.1186/s40168-015-0092-7.

  20. 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–6.

    Article  CAS  Google Scholar 

  21. Roesch LF, Casella G, Simell O, Krischer J, Wasserfall CH, Schatz D, Atkinson MA, Neu J, Triplett EW. Influence of fecal sample storage on bacterial community diversity. Open Microbiol J. 2009;3:40–6.

    Article  Google Scholar 

  22. 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–606. https://doi.org/10.1136/jclinpath-2016-204165.

    Article  Google Scholar 

  23. Saiki RK, Scharf S, Faloona F, Mullis KB, Horn GT, Erlich HA, Arnheim N. Enzymatic amplification of beta-globin genomic sequences and restriction site analysis for diagnosis of sickle cell anemia. Science (New York, NY). 1985;230(4732):1350–4.

    Article  CAS  Google Scholar 

  24. Klindworth A, Pruesse E, Schweer T, Peplies J, Quast C, Horn M, Glockner FO. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 2013;41(1):e1.

    Article  CAS  Google Scholar 

  25. Raju SC, Lagstrom S, Ellonen P, de Vos WM, Eriksson JG, Weiderpass E, Rounge TB. Reproducibility and repeatability of six high-throughput 16S rDNA sequencing protocols for microbiota profiling. J Microbiol Methods. 2018;147:76–86.

    Article  CAS  Google Scholar 

  26. 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(23):7537–41.

    Article  CAS  Google Scholar 

  27. Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics (Oxford, England). 2011;27(16):2194–200.

    Article  CAS  Google Scholar 

  28. Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, Peplies J, Glockner FO. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41(Database issue):D590–6.

    CAS  Google Scholar 

  29. Haegeman B, Hamelin J, Moriarty J, Neal P, Dushoff J, Weitz JS. Robust estimation of microbial diversity in theory and in practice. ISME J. 2013;7(6):1092–101.

    Article  Google Scholar 

  30. Bashan A, Gibson TE, Friedman J, Carey VJ, Weiss ST, Hohmann EL, Liu YY. Universality of human microbial dynamics. Nature. 2016;534(7606):259–62.

    Article  CAS  Google Scholar 

  31. Paulson JN, Stine OC, Bravo HC, Pop M. Differential abundance analysis for microbial marker-gene surveys. Nat Methods. 2013;10(12):1200–2.

    Article  CAS  Google Scholar 

  32. Paulson JN, Pop M, Bravo H. metagenomeSeq: Statistical analysis for sparse high-throughput sequencing. Bioconductor Package. 2013;1(0).

  33. Vogtmann E, Hua X, Zeller G, Sunagawa S, Voigt AY, Hercog R, Goedert JJ, Shi J, Bork P, Sinha R. Colorectal Cancer and the human gut microbiome: reproducibility with whole-genome shotgun sequencing. PLoS One. 2016;11(5):e0155362.

    Article  Google Scholar 

  34. Gorzelak MA, Gill SK, Tasnim N, Ahmadi-Vand Z, Jay M, Gibson DL. Methods for improving human gut microbiome data by reducing variability through sample processing and storage of stool. PLoS One. 2015;10(8):e0134802.

    Article  Google Scholar 

  35. Bosch LJ, Mongera S, Terhaar Sive Droste JS, Oort FA, van Turenhout ST, Penning MT, Louwagie J, Mulder CJ, van Engeland M, Carvalho B, et al. Analytical sensitivity and stability of DNA methylation testing in stool samples for colorectal cancer detection. Cell Oncol (Dordrecht). 2012;35(4):309–15.

    Article  CAS  Google Scholar 

  36. Matsson M, Wadstrom T. Methods removal of PCR inhibitors from human faecal samples through the use of an aqueous two-phase system for sample preparation prior to PCR. J Microbiol Methods. 1997;3(28):159–67.

    Google Scholar 

  37. Monteiro L, Bonnemaison D, Vekris A, Petry KG, Bonnet J, Vidal R, Cabrita J, Megraud F. Complex polysaccharides as PCR inhibitors in feces: helicobacter pylori model. J Clin Microbiol. 1997;35(4):995–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  38. Akane A, Matsubara K, Nakamura H, Takahashi S, Kimura K. Identification of the heme compound copurified with deoxyribonucleic acid (DNA) from bloodstains, a major inhibitor of polymerase chain reaction (PCR) amplification. J Forensic Sci. 1994;39(2):362–72.

    Article  CAS  Google Scholar 

  39. Ahn J, Sinha R, Pei Z, Dominianni C, Wu J, Shi J, Goedert JJ, Hayes RB, Yang L. Human gut microbiome and risk for colorectal cancer. J Natl Cancer Inst. 2013;105(24):1907–11.

    Article  CAS  Google Scholar 

  40. Dolfing J, Vos A, Bloem J, Ehlert PA, Naumova NB, Kuikman PJ. Microbial diversity in archived soils. Science (New York, NY). 2004;306(5697):813.

    Article  CAS  Google Scholar 

  41. Goberna M, Insam H, Klammer S, Pascual JA, Sanchez J. Microbial community structure at different depths in disturbed and undisturbed semiarid Mediterranean forest soils. Microb Ecol. 2005;50(3):315–26.

    Article  CAS  Google Scholar 

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

  43. Ott SJ, Musfeldt M, Timmis KN, Hampe J, Wenderoth DF, Schreiber S. In vitro alterations of intestinal bacterial microbiota in fecal samples during storage. Diagn Microbiol Infect Dis. 2004;50(4):237–45.

    Article  CAS  Google Scholar 

  44. Tzeneva VA, Salles JF, Naumova N, de Vos WM, Kuikman PJ, Dolfing J, Smidt H. Effect of soil sample preservation, compared to the effect of other environmental variables, on bacterial and eukaryotic diversity. Res Microbiol. 2009;160(2):89–98.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We would like to thank Catherine Stenstad, Elina Vinberg, Sofia Strand and Jorun Bratlie for coordination, biobanking and providing samples, Sonja Lagström and Pekka Ellonen for design of amplicon protocol and advise on sequencing, Hotan Esnaashari for DNA extraction, Alexander Kirkeby Eieland and Ellen Myrvang for amplicon generation and sequencing and Steinar Auensen for contributing with the high-performance computing knowledge.

Funding

The study was funded by The Cancer Registry of Norway funds. The funding body had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Availability of data and materials

The 16S rDNA sequences have been submitted to The European Bioinformatics Institute (EBI) metagenomics database (study accession number is ERP016903).

Author information

Authors and Affiliations

Authors

Contributions

TBR, OHA, TdL and GH designed the research. TBR, RM and JIN performed the experiments. TBR analysed the data and drafted the manuscript. All authors contributed to the writing and approved the final manuscript.

Corresponding author

Correspondence to Trine B. Rounge.

Ethics declarations

Ethics approval and consent to participate

The collection of iFOBT and archived samples was approved by the Regional Committees for Medical and Health Research Ethics in South-Eastern Norway (2011/1272 and 2010/3087 A, respectively). An additional ethical approval for the feasibility study was not required since all samples were anonymized and the purpose of the study was method related, as stated by the Norwegian regional ethics committee (Regional Committees for Medical and Health Research Ethics in South-Eastern Norway, ref. 2015/9).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests, with the exception of GH received payment from Amgen Norway for giving a lecture at a medical conference in 2017.

Publisher’s Note

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

Additional files

Additional file 1:

Distribution of OTU counts in all samples analysed. (PDF 513 kb)

Additional file 2:

Table with OTU counts for all samples. (TXT 13318 kb)

Additional file 3:

Taxonomical information for all OTUs. (TXT 9091 kb)

Additional file 4:

Alpha diversity measurements. (TXT 22 kb)

Additional file 5:

Beta diversity measurements. (XLSX 126 kb)

Additional file 6:

Differences in microbial composition between iFOBT and archived samples. (PDF 115 kb)

Additional file 7:

Statistical differences in microbiota composition between fresh frozen directly and fresh frozen after 48 h in room temperature. (PDF 92 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), 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 (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Rounge, T.B., Meisal, R., Nordby, J.I. et al. Evaluating gut microbiota profiles from archived fecal samples. BMC Gastroenterol 18, 171 (2018). https://doi.org/10.1186/s12876-018-0896-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12876-018-0896-6

Keywords