(2011) [16]). Sequence data generated in this study were submitted to the Sequence Read Archive with the study accession
number ERP001705. The dataset is available at http://www.ebi.ac.uk/ena/data/view/ERP001705. Taxonomical analysis For taxonomic grouping of the sequence reads, MEGAN V3.4 http://www-ab.informatik.Cell Cycle inhibitor uni-tuebingen.de/software/megan/welcome.html[23, 24] was used. First, the sequence reads were compared to a curated version of the SSUrdp database [25] using blastn with a maximum expectation value (E) of 10-5. To reflect the actual abundance behind every denoised sequence cluster, each entry in the blast result file was replicated as many times as the total number of reads that mapped to that query see more Chk inhibitor sequence (for detailed procedure and parameters see Siddiqui et al. (2011) [16]). When comparing the individual datasets using MEGAN, numbers of reads were normalized up to 100,000 for every dataset. Metastats, statistical methods ( http://metastats.cbcb.umd.edu/, [26, 27]) for detecting differentially abundant taxa, was used to reveal significant differences between IC urine microbiota and HF urine microbiota (taxonomy assessed in Siddiqui et al. 2011 [16]). This method employs a false discovery rate to improve specificity in high-complexity environments, and in addition handles sparsely sampled features
using Fisher’s exact test. The Metastats p – values at different taxon levels, which were assigned using MEGAN, are listed in Additional file 1: Table S1. A p – value ≤ 0.05 was considered significant. Comparative OTU based clustering analysis of IC and HF urine Numbers of operational taxonomical units (OTUs), rarefaction curves and diversity indices were calculated using MOTHUR v1.22.2 [28, 29] (see Table 1). To enable comparisons, the HF sequences generated in Siddiqui et al. (2011) [16] were reanalyzed along with the IC dataset from this study. Briefly, the sequences were aligned to the Silva 16S alignment as recommended by MOTHUR [29] – sequences not aligned or aligned outside of
where 95% of all of the sequences aligned were removed from the datasets. For an improved OTU clustering single linkage preclustering [30] was performed, allowing two nucleotides to differ between sequences, before clustering using average linkage. The processing was done both on each separate Pyruvate dehydrogenase sample and on pooled V1V2 and V6 sequences for both IC and HF samples. We also calculated the OTUs and Shannon index for normalized numbers of sequences for each separate sample [31]. A random number of reads, corresponding to the lowest number of sequences in a sample group, i.e. 2,720 for V1V2 and 2,988 for V6, was picked 100 times from each sequence set. These new sequence sets were processed through MOTHUR in the same fashion as the full sequence sets and the average of the resulting OTUs and Shannon values are shown in Additional file 2: Table S2.