The human gastrointestinal tract harbours a highly diverse population of commensal microbes, estimated as 1014 bacterial cells, participating in the breakdown of indigested food(Reference Hooper, Midtvedt and Gordon1, Reference Mueller, Saunier and Hanisch2). The composition of the microbiota develops during the first years of life and remains relatively stable throughout adulthood(Reference Ouwehand, Isolauri and Salminen3). Interestingly, after birth the maturation of immune responses reflects the microbial colonisation in the intestine(Reference Bauer, Williams and Smidt4, Reference Rautava and Walker5). It has been suggested that towards seniority specific changes in the microbial community(Reference Benno, Endo and Mizutani6) and immune defence occur again(Reference Gill, Darragh and Cross7). In the elderly, immune responses are weakened, predisposing seniors to a higher incidence of morbidity and mortality by infectious diseases(Reference Kumar and Burns8, Reference van Duin and Shaw9). Furthermore, the common use of non-steroidal anti-inflammatory drugs (NSAID) among the elderly is known to increase gastrointestinal ailments(Reference Bjarnason, Hayllar and MacPherson10–Reference Laine, Connors and Reicin12). However, the elderly immune responses can be enhanced by supplementing the diet with probiotics, implying that not only in children are the immune responses modified by intestinal bacteria(Reference Gill, Rutherfurd and Cross13). With improved life expectancy the share of the population over the age of 60 years is estimated to reach 2 billion by 2050(Reference Cohen14). In order to improve the quality of life of the growing elderly population and decrease health-care costs, the role of intestinal microbes in health necessitates deeper understanding.
Age- and country-specific differences in intestinal microbiota have been reported by several studies. Most commonly, decreasing numbers of bifidobacteria and increasing numbers of enterobacteria have been suggested(Reference Mueller, Saunier and Hanisch2, Reference Benno, Endo and Mizutani6, Reference Hayashi, Sakamoto and Kitahara15, Reference Bartosch, Fite and Macfarlane16) in the elderly. Some reports indicate changes also in the Lactobacillus and Bacteroides groups(Reference Tiihonen, Tynkkynen and Ouwehand17). Recently, also a reduced microbial activity has been described in elderly subjects(Reference Tiihonen, Tynkkynen and Ouwehand17, Reference Woodmansey, McMurdo and Macfarlane18). It is noteworthy that changes in the microbial composition have not only been linked to age but also to a condition with multiple health risk factors, namely obesity(Reference Turnbaugh, Ley and Mahowald19), and several gastrointestinal diseases such as Crohn's disease, irritable bowel disease and irritable bowel syndrome (IBS)(Reference Malinen, Rinttilä and Kajander20, Reference Kassinen, Krogius-Kurikka and Makivuokko21). This underlies the importance of increasing knowledge of the intestinal microbial composition.
The availability of microbial DNA sequences has revealed that a great proportion of the gastrointestinal microbes belong to unknown and uncultivable genera and species(Reference Hayashi, Sakamoto and Benno22). Thus it is important to include DNA-based methods that are independent of prior knowledge of bacteria when evaluating the overall changes in the microbial communities. The percentage guanine+cytosine (%G+C) profiling technique has been previously used to reveal changes in the microbial community structure(Reference Malinen, Rinttilä and Kajander20, Reference Kassinen, Krogius-Kurikka and Makivuokko21, Reference Apajalahti, Kettunen and Kettunen23–Reference Mäkivuokko, Saarinen and Ouwehand25). The technique exploits the diverging weight of microbial genomes due to different proportions of guanidine and cytosine bases in their genomes, and thus needs no prior information of specific sequences, allowing analysis of previously uncultivated microbes. The 16S rDNA sequencing technique, however, requires some prior information of the microbes present in the sample, since the broad-range primers used by the technique are common for known cultivated bacteria(Reference Lane, Pace and Olsen26). The combination of these two methods can provide useful information from samples with mixed populations containing hundreds of different species. In fact, the two methods have been used successfully previously in describing differences in the microbiota between IBS patients and healthy subjects(Reference Kassinen, Krogius-Kurikka and Makivuokko21).
In the present study we have characterised the intestinal microbiota of a subgroup of previously well-defined elderly subjects, with or without NSAID, and young adults(Reference Tiihonen, Tynkkynen and Ouwehand17) using %G+C profiling and 16S rDNA sequencing techniques in order to expose differences caused by age and NSAID use.
Materials and methods
Study subjects
A total of nine elderly NSAID users (mean age 80 (range 77–85) years) and nine elderly non-users (mean age 78 (range 70–83) years) living in a nursing home and fourteen young adults (mean age 28 (range 21–39) years) were recruited to the study in the city of Kuopio, Finland. The number of study subjects per group was originally planned to be fifteen but due to an unexpected high number of drop-outs, such as use of antibiotics, or diarrhoea, the number of especially elderly study subjects decreased. Significant differences in microbial composition and metabolites were found for the two elderly groups(Reference Tiihonen, Tynkkynen and Ouwehand17), prompting the separate analysis of the sequence data for the elderly NSAID users and non-users. The elderly were defined as NSAID users when they used the medication three or more times per week. The NSAID used included indomethacin, ibuprofen, ketoprofen, diclofenac, piroxicam, celecoxib, nimesulid and aspirin. The exclusion criteria were critical illnesses, inflammatory bowel disease, coeliac diseases and major malignancies in the gastrointestinal tract. Subjects using antibiotics during the study or within 2 months before the screening were excluded. The use of probiotic products was not allowed during the study. The inclusion and exclusion criteria were monitored during the run-in period before randomisation.
Throughout the sampling both the elderly and the young study subjects followed their habitual diet and most of them used the uniform meal service provided at least once per d. The calculated dietary intakes of the elderly living in the nursing home, based on the diets provided, were: energy, 6·17 MJ/d (6·5 MJ/d fibre included); protein, 66 g/d (18 % energy); fat, 52 g/d (32 % energy); cholesterol, 146 mg/d; carbohydrates, 181 g/d (50 % energy); dietary fibres, 20 g/d.
The present study was conducted according to the guidelines laid down in the Declaration of Helsinki and all procedures involving human subjects were approved by the Research Ethics Committee of the Hospital Districts of Northern Savo, Finland. The purpose of the study was explained to the participants, and written consent to the study was given by all of the subjects. Information on the age, sex, bowel function, diseases, allergies, current medication and pre-trial use of probiotics of the study subjects was obtained at the beginning of the study. The pre-trial characteristics of the study groups are presented in Table 1.
NSAID, non-steroidal anti-inflammatory drugs.
* Several drugs can have gastrointestinal side effects. No side effects by the medications used by study subjects were observed.
† Macrocol, sodium bicarbonate, sodium chloride, potassium chloride, sodium sulfate, lactulose, sodium picosulfate (common (>1/100) gastrointestinal side effects: diarrhoea; others – flatus).
‡ Enalapril, bisoprolol, spironolactone, losartan, aspirin, warfarin, metoprolol, dipyridamole+aspirin, bisoprolol, clopidogrel, dipyridamole, dipyridamole, furosemide, isosorbide dinitrate, indapapamide, propranolol (common (>1/100) gastrointestinal side effects: diarrhoea; others – flatus, vomiting, gastrointestinal bleeding).
§ Magnesium hydroxide, finasteride, ranitidine, macrocol, sodium bicarbonate, sodium chloride, potassium chloride, sodium sulfate, lactulose, sodium picosulfate (common (>1/100) side effects: diarrhoea, constipation; others – flatus).
∥ Olanzapine, mirtazapine, oxazepam, risperidone, melperone, thioridazine, chlorprothixene, temazepam, zopiclone, quetiapine, citalopram (very common (>1/10) gastrointestinal side effects: constipation; others – flatus, diarrhoea, vomiting).
¶ Galantamine, oxcarbazepine, levodopa+benserazide, entacapone, pramipexole, donepezil, memantine, sodium valporate, memantine, rivastigmine (very common (>1/10) gastrointestinal side effects: vomiting, diarrhoea; others – constipation).
Faecal samples
Spot samples of faeces were obtained three times from each subject over 3 consecutive weeks. Subjects were personally instructed to take the faecal samples according to the written protocol. All the equipment and documentation papers were provided to the subjects before sampling. All faecal samples were initially stored at − 20°C, and, transferred to the laboratory within 12 h of defecation and stored at − 70°C until analysed. The three samples obtained from each study subject were pooled to result in a representative individual sample.
Percentage guanine+cytosine profiling
Total genomic DNA of faecal microbes was extracted as previously described(Reference Tiihonen, Tynkkynen and Ouwehand17). In short, 1 g of faecal material was washed, microbes were collected and freeze-thawed to loosen their cell membranes, after which the microbes were lysed, the cell debris was removed and genomic DNA was precipitated.
The composition of microbiota in the faecal samples was monitored by %G+C profiling(Reference Apajalahti, Sarkilahti and Maki24). The method enables fractionation of total microbial DNA into bacterial subpopulations by their genomic G+C content. For fractioning, three individual DNA samples were pooled for each elderly donor pool (five pools in total; three pools of elderly without NSAID and two with NSAID) and four or five individual DNA samples in each of the young adult donor pools (three pools in total). One pool of three subjects for the elderly with NSAID group had to be removed from the %G+C-profile analysis due to breakdown of genomic DNA into smaller fragments during the extraction phase, making %G+C-profiling of that particular pool inaccurate compared with pools with intact DNA. However, this fragmented pool was included in the sequencing procedure, in which only intact 16S rDNA was needed. Each of the pools contained 300 μg genomic DNA.
Microbial enumeration
The total bacteria cell counts in faecal samples were determined by flow cytometry (FACSCalibur; Becton Dickinson, Franklin Lakes, NJ, USA) as described previously(Reference Apajalahti, Sarkilahti and Maki24). The bacterial fractions were recovered by suspending faecal samples in a buffer, followed by centrifugation and washing. The cell samples were diluted, fixed and stained with the fluorescent nucleic acid-binding dye SYTO 24 (Molecular Probes, Leiden, The Netherlands). The results are expressed as the quantity of bacteria/g faeces.
Sequencing procedure
For sequencing, six fractions containing %G+C areas of 20–30, 30–40, 40–50, 50–60, 60–70 and 70–80 %G+C were collected from the pools. All corresponding %G+C fractions in one study group (i.e. young adults, elderly and elderly with NSAID) were pooled, resulting in six %G+C fractions for each study group, and totally eighteen fractions. The DNA of all the eighteen fractions were used as templates to amplify 900 bp 16S rDNA regions for sequencing. The PCR amplification was performed using a Dynazyme II polymerase kit (Finnzymes, Espoo, Finland) according to the manufacturer's instructions and broad-range primers pA 5′-AGAGTTTGATCCTGGCTCAG-3′ and pE′ 5′-CCGTCAATTCCTTTGAGTTT-3′(Reference Lane, Pace and Olsen26, Reference Edwards, Rogall and Blocker27) corresponding to the 8–28 and 928–908 positions in the Escherichia coli 16S rDNA sequence. Amplification products were purified with Sephacryl S-400 colon (Amersham Biosciences, Little Chalfont, Bucks, UK) and cloned in pDrive cloning vector (Qiagen, Hilden, Germany) in E. coli. Colonies were picked and sequenced from both directions with T7 and SP6 primers using BigDye Chemistry and analysed on an ABI 3730 Capillary DNA Sequencer (Applied Biosystems, Foster City, CA, USA). The obtained sequences were assembled and edited using Gap4 from the Staden package(Reference Staden, Beal and Bonfield28).
Spermidine analysis
The spermidine concentrations in faeces were analysed using the method described by Saarinen(Reference Saarinen29).
Phylogenetic and statistical analysis of sequence data
%G+C fractions were analysed with ANOVA applying linear mixed models(Reference Pinheiro and Bates30) inside each fraction, and then computing post hoc t tests whenever a significant difference between the groups was found. Spermidine data were analysed with ANOVA applying a linear model and post hoc t tests. The sequences were phylogenetically aligned with the Ribosomal Database Project (RDP) II My RDP on-line-tool (https://rdp.cme.msu.edu/login/myrdp) and statistically compared using the RDP II Library Compare Tool Analysis on-line-tool (http://rdp.cme.msu.edu/seqmatch/) with a 95 % confidence threshold for genus determination, which uses a RDP-naive Bayesian classifier to provide rapid classification of library sequences into the bacterial taxonomy. The equality of the sequence proportions between the groups in Fig. 1 was tested with a standard multi-population proportions test whenever all the frequencies were larger than five, and if statistically significant non-equality was found, the proportions were analysed further with a pair-wise proportions test. If the proportions were smaller than five, a special proportions test for small counts was applied as explained in Wang et al. (Reference Wang, Garrity and Tiedje31), and first introduced for digital gene expression profiles in Audic & Claverie(Reference Audic and Claverie32). The comparison analysis enabled the calculation of the most abundant strain of each genus using a similarity score or a sequence match score (S_ab) comparison reporting the number of (unique) seven-base oligomers shared between an input sequence and a given RDP sequence divided by the lowest number of unique oligos in either of the two sequences. In addition to RDP software (Center for Microbial Ecology at Michigan State University, East Lansing, MI, USA), the statistical programming language R was used to perform analyses (http://www.r-project.org/).
Results
The total microbial numbers and the corresponding %G+C-profiles were compared to reveal the proportional differences in the total microbiota profiles of the three study groups (Fig. 2). The young adult group had statistically higher numbers of microbes in all cumulated 10 %G+C-fractions between the 30–69 %G+C area when compared with either of the elderly groups (P < 0·001 and P < 0·0003, respectively, in each 5 % fraction) and in the 20–29 %G+C area when compared with the elderly group without NSAID (P < 0·006). The two elderly groups differed from each other in the 20–59 %G+C area (P < 0·0002, respectively, in each 5 % fraction), the microbial numbers being higher in the elderly group with NSAID.
A total of 2112 high-quality sequences were obtained by sequencing of the six fractioned DNA samples of each study group. The total sums of sequences obtained were 803 sequences in the young adults, 702 sequences in the elderly without NSAID and 610 sequences in the elderly with NSAID. In Fig. 1 these aligned sequences are divided into six phyla and fifty-one presently known microbial genera; the bacterial genera differing significantly (P < 0·001) between the study groups are marked with arrows. The majority of the microbes in each study group were found to be associated to the Firmicutes phylum followed by the Bacteroidetes and Actinobacteria phyla. Comparison of young adults with either of the elderly groups revealed more differences in the number of genus-specific sequences than the comparison of the two elderly groups. The young adults exhibited higher faecal numbers of Rumicococcus, Roseburia, Coprobacillus and Dialister spp., all belonging to Firmicutes, than the elderly without NSAID, suggesting an age-dependent shift in genera. Lactobacillus, Streptococcus and Bacteroides were present in higher numbers in the elderly. Compared with the elderly NSAID users, the samples obtained from the young adults exhibited more Dialister, Coprobacillus and Collinsella. The significant differences between the elderly groups were the higher amounts of Lactobacillus and Collinsella in the elderly without NSAID and Roseburia in the NSAID users.
In order to indicate the degree of match of study sequences to each named phylotype in the RDP-II database (http://rdp.cme.msu.edu/), S_ab-scores were calculated on single species in each genus with statistical significance between study groups present in Fig. 1; species with S_ab-scores higher than 0·97 were considered significant. The S_ab-scores for the Ruminococcus and Coprobacillus genera had no significant matches, while the Roseburia, Lactobacillus, Streptococcus, Dialister, Bacteroides and Collinsella genera each had one or two matching strains with the RDP-II database; the most abundant bacterial species were, respectively, Roseburia intestinalis, Lactobacillus ruminis, L. sobrius, Streptococcus lutiensis, Dialister invisus, Bacteroides coprocola, B. dorei and Collinsella aerofaciens.
The concentrations of spermidine in the faeces of the elderly (median 76·7 and 67·6 nmol/g faeces, respectively, in the elderly without NSAID and with NSAID) were monitored to be significantly lower when compared with the concentrations in samples obtained from the younger adults (median 160·4 nmol/g faeces; P = 0·004 and P = 0·007, respectively).
Discussion
In the present study we characterised the intestinal microbiota of a subgroup of previously well-defined elderly subjects, with or without NSAID drugs, and young adults(Reference Tiihonen, Tynkkynen and Ouwehand17) using %G+C profiling and 16S rDNA sequencing techniques in order to expose differences caused by age and NSAID use. Our previous study(Reference Tiihonen, Tynkkynen and Ouwehand17) focused mainly on the microbial metabolism and intestinal immunological status but analysed also differences in faecal microbiota. The study showed that specific immune responses as well as generally the microbial metabolism were decreased, and furthermore that, for example, the proportion of facultative anaerobic microbes was increased in the elderly.
The %G+C-profiling technique has recently been used as a tool for describing the total microbial community in human(Reference Kassinen, Krogius-Kurikka and Makivuokko21, Reference Mäkivuokko, Saarinen and Ouwehand25), mouse(Reference Apajalahti, Kettunen and Kettunen23) and piglet(Reference Fava, Makivuokko and Siljander-Rasi33) intestines. The %G+C-profiling enables microbiota analysis that includes previously unknown or presently unculturable microbial genera since it exploits only the different proportion of guanidine and cytosine bases in their genomes. Moreover, %G+C fractioning of microbiota DNA before 16S rRNA gene sequencing enhances the possibility of less abundant microbial species to be revealed in specific fractions.
The obtained %G+C-profiles support the sequencing data well by depicting that young adults have the highest numbers of microbes especially in the Firmicutes-rich %G+C area (low %G+C bacteria). In the young adults 87 % of sequences were classified to Firmicutes. Similarly, in the elderly groups, 78–82 % of the sequences grouped to this phylum; however, the significantly higher number of bacteria in the NSAID users shown in %G+C profiling did not result in higher numbers of sequences. Also in the high %G+C areas (>50 %G+C), a difference between the %G+C and sequencing data can be noted. In the elderly NSAID users the 50–59 %G+C area was significantly higher compared with that of the non-users, but our sequencing data do not reveal an increase in bacterial species. As %G+C profiling is a totally culture- and probe-independent analysis method portraying only the concentration differences of the microbial genomic DNA, this discrepancy is most probably due to higher numbers of bacteria belonging to the same genera or problems in amplifying high-G+C sequences, which affected especially the Bifidobacterium genus. One explanation for the lack of bifidobacteria in the present study is due to the primer used in PCR. The pE′ primer 5′-CCGTCAATTCCTTTGAGTTT-3′ has a mismatch C/T (highlighted in bold) in the middle. This mismatch can abolish the amplification of bifidobacteria 16S rRNA genes. Others have also reported the lack of bifidobacteria in 16S rRNA libraries, possibly due to mismatches in primers or imperfect DNA denaturation during PCR(Reference Hayashi, Sakamoto and Kitahara15, Reference Suau, Bonnet and Sutren34).
The %G+C profiling showed that the young adults had statistically higher numbers of microbes in the %G+C area of Firmicutes, Bacteroidetes and Actinobacteria, the major phyla reported in the human intestine, compared with either elderly group. The most prominent difference dependent on age, however, was a clear reduction in numbers of Ruminococcus and Roseburia spp., the known butyrate producers in Clostridium cluster XIVa, in the elderly. This is consistent with our earlier studies showing a decrease in the Clostridium coccoides–Eubacterium rectale group in the elderly(Reference Tiihonen, Tynkkynen and Ouwehand17). Similarly, the elderly with NSAID had higher numbers of bacteria in most %G+C fractions compared with the non-users. Therefore, both ageing and the use of specific medication can deeply alter the abundance and composition of the microbial community. The sequencing revealed higher amounts of lactobacilli in the elderly groups compared with the young adults, which is in accordance with our earlier findings(Reference Tiihonen, Tynkkynen and Ouwehand17).
The major anaerobic microbes in the gastrointestinal tract are Bacteroides species. These organisms can function as symbiotic partners in the host physiology or they can cause serious diseases(Reference Wexler35). Recent studies have indicated the link between levels of the Bacteroides and Firmicutes and obesity(Reference Ley, Backhed and Turnbaugh36–Reference Turnbaugh, Backhed and Fulton38). Ley et al. (Reference Ley, Backhed and Turnbaugh36) showed that obesity alters gut microbiota by increasing the relative abundance of Firmicutes over Bacteroidetes. Evidently in mice with diet-induced obesity the Western diet nutrient milieu provides a competitive advantage to the Firmicutes(Reference Turnbaugh, Backhed and Fulton38).
Interestingly, also in the present study, major shifts in the microbiota occurred in the levels of Bacteroidetes and Firmicutes, suggesting perhaps differences in the diet of the elderly and young adults and also in the physiology of the intestinal tract. This is consistent with our earlier finding of a higher proportion of anaerobic bacteria in faecal samples from elderly subjects(Reference Tiihonen, Tynkkynen and Ouwehand17, Reference Woodmansey, McMurdo and Macfarlane18). Furthermore, some of the major differences between the microbiota of healthy elderly and young adults reported previously in German(Reference Mueller, Saunier and Hanisch2) and Japanese populations(Reference Hayashi, Sakamoto and Kitahara15) correlate with those presented in the present study, that is, a higher proportion of Bacteroides in the elderly and a higher proportion of Ruminococcus in the young adults. However, Italian elderly were reported to have lower levels of Bacteroides when compared with young adults(Reference Mueller, Saunier and Hanisch2). These results indicate strong country specificity in diversity of the elderly microbiota.
The metabolic activity of the intestinal microbes has been reported by us to be decreased in elderly subjects(Reference Tiihonen, Tynkkynen and Ouwehand17, Reference Woodmansey, McMurdo and Macfarlane18). Changes in the microbial community structure described in the present study may offer some explanations or hypotheses to the detected differences. R. intestinalis and Eubacterium rectale are known butyrate-producing species belonging to Firmicutes(Reference Barcenilla, Pryde and Martin39, Reference Duncan, Hold and Barcenilla40). Roseburia sp. can convert acetate to butyrate, which is the major energy source for colonocytes and has been reported to reduce risk for developing colon cancer and colitis(Reference Louis, Scott and Duncan41). During a low-carbohydrate diet, the levels of Roseburia and Eubacterium rectale decrease, causing a drop in butyrate concentrations(Reference Duncan, Belenguer and Holtrop42), which emphasises the importance of the diet. In the present study, reduction of Roseburia was described in the elderly, and it may partially explain the described reduced concentrations of butyrate in elderly subjects compared with young adults as reported by Tiihonen et al. (Reference Tiihonen, Tynkkynen and Ouwehand17). C. aerofaciens is abundantly found in the human intestine(Reference Kageyama, Sakamoto and Benno43). Kassinen et al. (Reference Kassinen, Krogius-Kurikka and Makivuokko21) reported that the amount of Collinsella in IBS patients was significantly reduced. In the present study the number of C. aerofaciens was reduced in the elderly subjects with NSAID, while elderly subjects without NSAID and young adults had high numbers of C. aerofaciens. The use of NSAID in the study with IBS patients was not recorded but due to common gastrointestinal discomfort and pain in subjects with IBS, NSAID may have been used by the subjects.
Ruminococci are abundant in the human gastrointestinal tract, comprising 2·5 % of total bacteria(Reference Zoetendal, Ben-Amor and Harmsen44). These Clostridium coccoides group XIVa-belonging bacteria are cellulolytic, found on solid particles in faeces(Reference Walker, Duncan and Harmsen45) and have been linked to acetate and butyrate production(Reference Duncan, Holtrop and Lobley46). Furthermore, R. intestinalis has been suggested to be a polyamine producer(Reference Matsumoto and Benno47). Putrescine, spermidine and spermine are important for cell proliferation and differentiation and mucosal barrier function (for reviews, see Thomas & Thomas(Reference Thomas and Thomas48) and Larque et al. (Reference Larque, Sabater-Molina and Zamora49)). In the present study we found that the concentration of spermidine in the faecal samples was decreased in the elderly individuals when compared with the young adults. The reduced production of this polyamine, which is essential for the maintenance and function of rapidly renewing gastrointestinal epithelia(Reference Wang, Feith and Welsh50, Reference Seiler and Raul51), can most probably alter the structure of the epithelial layer, thus affecting the adherence of the commensal intestinal microbiota, and therefore reducing the overall microbial numbers in the gut. Lowered microbial numbers, together with an altered epithelial layer surface, can have a significant effect on the digestibility and breakdown of food components entering the colon, enhancing development of constipation and reducing nutrient intake.
In conclusion, further research using a larger set of faecal donors of different ages could confirm the relationship between changes in shares of specific species and the consequences in their activity, especially regarding metabolism of important compounds for colon health, such as butyrate or spermidine. We propose that in future research, measurements of, for example, specifically two prominent intestinal strains, R. intestinalis and C. aerofaciens, should be included as potentially useful biomarkers for monitoring changes in the intestinal conditions.
Acknowledgements
We especially appreciate the help of the organisers of the sampling in the Pohjois-Savo Vocational Institute in Siilinjärvi and the Harjula Hospital in the District of Northern Savo, and Essi Sarkkinen and Henna Karvonen from Foodfiles Ltd for organising the sampling procedure. Tanja Horko is acknowledged for technical assistance and Janne Nikkilä from Genolyze Ltd for statistical analyses. The present study was supported by the National Technology Agency of Finland (TEKES), grant no. 129/31/03.
H. M. was responsible for the design of the study, performing %G+C-analyses, sequence comparisons, interpretation of the results and writing the manuscript. K. T. was responsible for the design of the study, managing the clinical trial, interpretation of the results and writing the manuscript. S. T. was responsible for the design of the study, interpretation of the results and writing the manuscript. L. P. was responsible for the sequencing procedure, interpretation of results and writing the manuscript. N. R. was responsible for the design of the study, interpretation of the results and writing the manuscript.
There are no conflicts of interests.