Comparison of MST and UPGMA The geographic dependency found in UP

Comparison of MST and UPGMA The geographic dependency found in UPGMAs but not in MSTs could be explained by the different approaches of sequence-based versus allelic profile-based comparison. Sequences with fewer differences are grouped close together in the UPGMA whereas in MSTs all sequences which differ in at least one nucleotide have the same distance to each

selleck chemical other. Thus the UPGMA seems to be more suitable for showing geographical relationships between strains of highly diverse populations. The CCs identified by goeBURST were grouped together also in UPGMA analysis. Similarly Yan et al. observed the grouping of CCs identified by eBURST in high monophyletic clades of UPGMA analysis [15]. Conclusions The generated data reveal a high genetic diversity

for all V. parahaemolyticus strain subsets analyzed, with a high proportion of new alleles and STs discovered, typical for environmental strain collections. Clusters of strains on nucleotide level contained mainly strains originating from one continent, but no exclusive clusters for the distinct continents were identified. STs and pSTs were either supra-regionally distributed or exclusively present in one region. Using AA-MLST instead of MLST in the goeBURST analysis allowed reliable identification of closely related strains (pSTs were SLVs), independent of their geographic origin. In contrast the application of MLST is more useful to recognize relationships in an epidemiological context by creating distinct CCs. In general AZD1480 pubMLST database reflects only the diversity of so far S63845 purchase analyzed strains, and may not represent the natural diversity of the V. parahaemolyticus population as also indicated by our rarefaction analysis. Further analysis of strains of diverse origins may help to complete the database and to keep pace with new evolving genotypes. Availability of supporting data The data sets and additional figures supporting the results of this article are included in Additional files 1, 2, 3, 4 and 5. Acknowledgements We acknowledge Kathrin Oeleker for assistance in performing PCR and strain cultivation. The Montelukast Sodium project was funded by

the German Ministry of Education and Research (BMBF) within the VibrioNet project. Electronic supplementary material Additional file 1: Table S1: Characteristics and allelic profiles of V. parahaemolyticus strains included within this study. (PDF 151 KB) Additional file 2: Tables S2: AA-MLST profiles and properties of each allele on peptide level (numbers, sequences and frequencies). (XLSX 42 KB) Additional file 3: Figure S1: Population snapshot based on MLST profiles of pubMLST dataset. Coloring depends on geographical origin of isolates: Asia (red), South America (light green), North America (dark green), Africa (yellow) and Europe (blue). Size of circles represents number of isolates with the corresponding ST. STs that differ in one allele are connected via black lines.

Cells with annexin V (+) and PI (−) were deemed

Cells with annexin V (+) and PI (−) were deemed GDC-0068 cell line early apoptotic cells. Cells with both annexin V (+) and PI (+) were deemed late apoptotic cells. TUNEL assay To identify apoptosis in the transfected cells, we utilized the dead-end colorimetric TUNEL system kit (Promega, Madison, USA) to measure DNA fragmentation and caspase-3 activation in the GKN1 transfected cells, according to the manufacturer’s instructions. Briefly, cells were fixed in 4% paraformaldehyde solution for 25 min at room temperature, rinsed in PBS, and permeabilized by incubating the slides in 0.2%

Triton X-100 solution. Cells were then incubated with a terminal deoxynucleotidyl transferase (TdT) reaction mixture containing biotinylated nucleotides and TdT at 37°C for 60 min, and rinsed with 1 × SSC (sodium chloride-sodium citrate buffer) and PBS. Next, streptavidin HRP was added to the cells, and the cell

slides were stained with 3,3′-diaminobenzidine color solution. Finally, cells were examined under a light microscope and the number of positive cells was counted and summarized from a total of 10 high power fields. CP673451 in vivo Cell cycle analysis To analyze cell cycle distribution, transfected cells were grown and treated with 25 M olomoucine (Santa Cruz Biotechnologies, Santa Cruz, USA) for 1 h, and then incubated with regular culture medium for an additional 1 h [13]. Cells were then collected and subjected to cell cycle analysis by flow cytometry as described in the previous section. Captisol concentration Sensitivity to 5-FU treatment To detect the role of GKN1 in mediating sensitivity of gastric cancer cells to 5-FU Amisulpride treatment, we grew and treated GKN1 transfected tumor cells with 5-FU (Sigma) or DMSO for 24 h and 48 h. Concentrations of

5-FU ranged from 0.25 to 1.0 mmol/L. The apoptosis rate from these cells was detected by flow cytometry as previously described. cDNA microarray analysis To perform cDNA microarray analysis, total cellular RNA from GKN1-transfected and vector-control tumor cells were isolated with the Trizol® Reagent (Invitrogen). RNA was then reversely transcribed into cDNA using the TrueLabeling-AMP Linear RNA amplification kit (Superarray, Frederick, MD, USA), and then converted into biotin-labeled cRNA using biotin-16-UTP and an in vitro transcription kit (Roche, Basel, Switzerland). The newly synthesized cRNA probes were then purified with the ArrayGrade cRNA cleanup kit (Superarray) and then added to the pretreated Oligo GEArrays Human Apoptosis Microarray (OHS-012 from Superarray) that contains 112 apoptosis-related genes. The microarray was then hybridized overnight at 42oC. The next day, the hybridized arrays were washed and detected by chemiluminescence according to the manufacturer’s instructions (Pierce). The data were analyzed using GEArray Expression Analysis software (Superarray). If spot intensity increased by more than two fold, this gene was deemed upregulated.

4 [16] Hydrate Ridge ANME-2a/2c and SRB consortia 1 4 MPa Fed-bat

4 [16] Hydrate Ridge ANME-2a/2c and SRB consortia 1.4 MPa Fed-batch 7.5 [9] Conclusions After 286 days incubation in a simulated cold seep environment under

high methane pressure, ANME-2 and SRB in the selleck kinase inhibitor sediment from Captain Arutyunov Mud Volcano were enriched. Based on biovolume calculation, the populations of ANME-2 and SRB increased for 12.5 times and 8.4 times. Within total biomass volume, 99.7% was accounted from aggregates. Therefore the incubation condition apparently favoured the cells to form aggregates, especially in small size (2<Ø≤5 μm), rather than to selleck chemicals live as single cells. No aggregate bigger than 15 μm in diameter was observed; they apparently divided after reaching a critical size. Based on the 16S rRNA gene clone library, the archaeal diversity was low, and contained only ANME-2 (88%) and MBG-D (12%). In contrast, the bacterial community was highly diverse. Methods Incubation condition In a previous Nepicastat study, the sediment sample originally from Captain Arutyunov Mud Volcano (Gulf of Cadiz, North East Atlantic) was diluted 12 times with artificial sea water medium and incubated in a continuous high-pressure bioreactor at 15°C [11]. This bioreactor system was a simulator for cold seep ecosystems, where sulphate and high-pressure methane were supplied. Because the high apparent affinity for methane (37 mM) in SR-AOM reaction

and low dissolubility of methane in seawater (1.3 mM at 15°C at ambient pressure), it is necessary to supply high pressure methane to obtain high concentration of dissolved methane which can be directly used by microorganisms for high in vitro SR-AOM activity [11]. During this research, the reactor was operated in a fed-batch mode or a continuous mode. When it was in fed-batch mode, the methane pressures were switched between 1, 4.5 and 8 MPa. When it was in continuous mode, the methane pressure was either 1 or 8 MPa and the flow rate was 0.1 ml/min (HRT 100 hours). The SR-AOM activities under different operational conditions have been described previously [11]. To take a slurry sample, the

incubation vessel was open under a nitrogen atmosphere and manually stirred to make the slurry sample homogeneous. The slurry samples before (S1) and Dimethyl sulfoxide after (S2) 286 days incubation were fixed in 4% formaldehyde and stored at 4°C for cell staining. Additional slurry from S2 was stored at -20°C for DNA extraction and clone library analysis. Cell and aggregates quantification To assess the number and the size of cells and aggregates, DAPI (4′, 6-diamidino-2-phenylindole) staining was performed on S1 (after 2000 times dilution) and S2 (after 5600 times dilution). Subsequently, the samples were filtrated onto a circular GTTP polycarbonate filter (0.2 μm, Millipore, Germany) with a diameter of 2.5 cm. The number of cells (or aggregates) was quantified under a microscope (Zeiss, Carl Zeiss Microimaging GmbH, Germany) at 1,000 times magnification. The diameter of a single cell was assumed as 0.

7 1 7 16 5 22 3 318 1 4 16 5 24 7 Serogroup C1                  

7 1.7 16.5 22.3 318 1.4 16.5 24.7 Serogroup C1                             Choleraesuis c 1,5 0 0 0.03 4.2 0 0 0.05 4.3 0.03 0 0.02 2.0 Grampian r l,w 0 0 0 0 0 0 0 0 0 0 0 0 Hissar c 1,2 0 0 0 0 0 0 0 0 0 0 0 0 Redba z10 z35 0 0 0 0 0 0 0 0 0 0 0 0 Serogroup C2-C3                             Blockley k 1,5 0 0 0.18 0 0 0 0.23 0 0.05 0 0.14 0 Albany Z4,z24 – 0 0 0.05 4.7 0.6 0 0.09 3.4 0.03 0 0.10 4.9 Serogroup D1                             Enteritidis [f],g.m. [p] [1, 7] 3.8 5.2 13.1 22.7 9.8 1.8 14.10 22.9 4.7 4.5 18.6 24.4 Serogroup E                             Anatum e,h 1,6: [z64] 0.5 0.6 0.47 1.0 0 0 0.7 1.1 0.64 0.6 0.54 0.7 Serogroup G                      

      Havana f,g, [s] – 0.2 1.2 0.08 0 0.6 0.7 0.089 0.1 0.27 0.8 0.07 0 Total Salmonellae find more   2038 924 37442 529 164 717 35661 2557 3743 665 36214 2228 adata from Salmonella Annual Summary for clinical Salmonella isolates from nonhuman and human sources reported to the Disease Control and Prevention (CDC) and the USDA National Veterinary Services Laboratory (NSVL), USA. bdata from Annual Report and Accounts 2008/2009 of Veterinary Laboratory Agency, Department check details of Environment, Food and Rural Affairs, United Kingdom. cdata from the Disease Control and Prevention (CDC), Taiwan. Discussion As one of main pathogen to cause foodborne diseases, Salmonella has been frequently reported

among different animal sources, especially more divergent Salmonella serovars found in chickens [34]. With the limited serovars in 164 chicken isolates, serogroups C2, D, E and G were restricted in one SYN-117 datasheet county and serogroup B and C1 were found in all three counties (Table 2), suggesting possibly that serogroup B and C1 isolates may be

more adapted to chicken. In human isolates, we found that the serovar number in each serogroup were not associated positively with the serogroup prevalence, such as highest serovar number in low prevalent serogroup C1 vs lower serovar number in high prevalent serogroup B and serogroup D (Table 4). These results imply that serogroup C1 may occasionally infect human isolates. Further, serovars are determined by flagellins: H1 and H2 antigens encoded by fliC and fljB. As one of the most important immunogens, flagellin interacts with the toll-like receptor 5 (TLR5) to activate NFκB pathway and proinflammatory PtdIns(3,4)P2 genes to regulate innate and adaptive immune system [35–38]. However, aflagellar serovars S. Pullorum and S. Gallinarum cause more severe infection than flagellar serovars in chicken because of aflagellar S. Typhimurium could avoid the TLR5 regulation of IL-1β expression and polymorphonuclear cell infiltration in gut [39]. Such evasion of TLR5 is critical for survival of flagellar bacteria at muscos [40]. [In the present study, we found that i of H1 antigen and lack of H2 antigen were the common antigens for all serogroups in human isolates (Table 4).

Whatever results Stuart et al achieved between picosecond and fe

Whatever results Stuart et al. achieved between picosecond and femtosecond pulses, we acquired it within the femtosecond pulse regime. www.selleckchem.com/products/elacridar-gf120918.html For example, they discovered that the damage area generated by the 500-fs pulse in fused silica glass was twice as much smaller than that produced by the 900-ps pulse. Figure 2 Interaction of femtosecond laser pulses of different

pulse-width sizes with glass surface. Schematic representation of glass irradiation with femtosecond laser pulses with pulse widths of (a) 214, (b) 428, and (c) 714 fs (schematic not to the scale). Figure 3 Microholes drilled via different pulse-width sizes. Microholes drilled by femtosecond laser pulses with pulse widths of (a) 214 and (b) 714 fs at 16-W average laser power and 0.5-ms dwell time, 13-MHz repetition rate. Even though we did not work in the picosecond pulse duration regime,

we obtained similar result as we increased the pulse width in the femtosecond click here regime. Figure 3 shows the SEM images of the microholes drilled by femtosecond laser pulses at 13-MHz repetition rate for 0.5-ms dwell time with pulse widths of 214 and 714 fs, respectively. The diameters of these microholes are approximately on average 12 and 21 μm, respectively. The size of microhole represent the GSK2245840 ic50 amount of material removed from the target; larger diameter means larger amount of material removal compare to smaller hole diameter. The life span of the plasma is also an important factor. In the current investigation, the turbulence created in the plasma due to the interactions between nitrogen gas and plasma species lengthens the plasma life. Since the longer pulses spend a significant portion of their (-)-p-Bromotetramisole Oxalate duration traveling through

previously formed plasma, as depicted in Figure 2, the energy transmitted via longer pulse is not enough to ablate the material upon contact with the target material. Rather, this transmitted energy gets stored in the top part of the lattice and gets transferred into the bulk in all directions, making the target temperature rise in the area surrounding the irradiated spot. This makes molecules to become loose to form a larger pool of molten material. As a result, the subsequent longer pulses expel large particles and droplets into the plasma upon contacting the molten pool. On the contrary, the interaction of the short pulses with the target surface does not rise as much high temperature which creates shallow molten pool. Hence, the material removed from the target is composed of smaller particles and droplets. The size of the plasma species and the temperature rise of the target surface greatly affect the type of nanotips that grow on the target surface. Figure 4 shows SEM images of the randomly selected spots from the irradiated target surface with 214-fs laser pulses.

Interestingly, tumor lysates from TaxMTD–treated mice contained h

Interestingly, tumor lysates from TaxMTD–treated mice contained higher levels of cathepsin activity and mRNA. As infiltrating immune cells are the primary source of cathepsins in these tumors, we reasoned that tumors may mobilize cathepsin-positive cells from the bone marrow after TaxMTD treatment to promote recovery from the cytotoxic assault, potentially explaining why cathepsin inhibition in the context of TaxMTD treatment is more

effective than treating with either drug alone. Indeed, increased cathepsin activity-positive cells were found in the blood 48 hours after TaxMTD treatment. Our current data also suggests selleckchem that cathepsin inhibition specifically impairs the development of lung metastases. These analyses clearly support a therapeutic IWP-2 purchase benefit from adding cathepsin inhibition to chemotherapeutics in the treatment of breast SAR302503 research buy Cancer and the prevention of metastases. O180 The Effect of the PAX2 Oncogene on the Tumor Microinvironment, Tumor Progression and its Potential as a Therapeutic

Target for Prostate Cancer Carlton Donald 1 1 Phigenix, Inc., Atlanta, GA, USA Inhibition of cell death is a critical pathophysiological factor that contributes to the initiation and progression of cancer. Recently, much attention has focused on developing therapeutic agents aimed at cancer cell survival pathways involving factors such as MEK kinase Astemizole and AKT. Unattenuated, tumour-associated expression of PAX2, a transcriptional regulator implicated in oncogenesis and cancer development,

has been observed to play a direct role in these pathways. PAX2 expression is aberrantly turned on in a number of cancers such as Wilm’s Tumor, breast, ovarian, bladder and prostate. We have discovered a novel mechanism by which PAX2 promotes cancer cell survival through the suppression of the host defense peptide and putative tumor suppressor Human Beta Defensin-1 (hBD-1). Our current findings provide the first indication of the cellular factors responsible for deregulated PAX2 expression in prostate cancer and how targeting these factors promote cancer cell death. Collectively, these data offers substantial evidence of the therapeutic potential of inhibiting PAX2 for the treating prostate cancer. O181 Targeting the Tumor Stroma – a Novel Therapeutic Strategy Based on Separate Analysis of the Malignant and Stromal Cell Compartments in Brain Tumors Jian Wang 1 , Anne M.

On the other hand, biotinylated purified Bt 18 toxin

On the other hand, biotinylated purified Bt 18 toxin tagged with FITC-conjugated streptavidin appeared green or greenish yellow (as a result of see more overlap with Alexa Fluor® 594). For CEM-SS (Figure 5), biotinylated purified Bt 18 toxin appeared at all test intervals (1, 2, 12 and 24 hours). The intensity and extent of staining increased relatively with increased incubation period. The biotinylated toxin was seen at the periphery of

the cell where the cell membrane was located. For human T lymphocytes (Figure 6), biotinylated purified Bt 18 toxin did not appear at all test intervals except at 24 hours. Compared to CEM-SS at the same Adavosertib interval, the intensity and extent of staining was very much less remarkable. Figure 5 Confocal microscopic appearance of biotinylated purified Bt 18 toxin-treated CEM-SS. CEM-SS cells were treated with fixed concentration of biotinylated purified Bt 18 toxin at various time intervals. The biotinylated toxin was then tagged

with FITC-conjugated streptavidin (green). The biotinylated toxin was found at the periphery of treated CEM-SS cells. Increased binding of the biotinylated toxin was observed with increased incubation period. (A) Untreated negative control at 24 hours. (B) to (E) Treated cells at 1, 2, 12 and 24 hours respectively. Magnification: 630×. B = Biotinylated purified Bt 18 toxin, M = Cell membrane, N = Nucleus. Selleck GDC0068 ID-8 Bar = 10 μm. Figure 6 Confocal microscopic appearance of biotinylated purified Bt 18 toxin-treated human T lymphocytes. Human T lymphocytes were treated with fixed concentration of biotinylated purified Bt 18 toxin at various time intervals. The biotinylated toxin was then tagged with FITC-conjugated streptavidin (green). Biotinylated toxin was not found at all intervals except 24 hours, where the binding was very minimal. (A) Untreated negative control at 24 hours. (B) to (E) Treated cells at 1, 2, 12 and 24 hours respectively. Magnification:

630 ×. B = Biotinylated purified Bt 18 toxin, M = Cell membrane, N = Nucleus. Bar = 10 μm. Discussion Binding studies related to Bt have largely been carried out in insects mainly for characterisation of Bt toxins and the determination of Bt toxin resistance in insects [14]. According to data from cell viability assays in this study, biotinylation did not affect the biological activity of the toxin. This phenomenon was observed in other studies [15, 16]. Results from the homologous competitive assays suggested that purified Bt 18 toxin had a higher affinity for CEM-SS, followed by CCRF-SB and CCRF-HSB-2 as the dissociation constant is inversely proportional to the binding affinity. For MCF-7, the dissociation constant could not be obtained because the inhibitory concentration (IC50) was not achieved.

There is a large component of ecological restoration that still p

There is a large component of ecological restoration that still places considerable value on past ecosystems and seeks to restore the system’s characteristics to its past state. Valuing the past when the past is not an accurate indicator for the future may fulfill a nostalgic need but may ultimately be VX-680 cost counterproductive in terms of achieving realistic and lasting restoration outcomes. Our results indicate a significant gap

between theory and practice—understandable for the early stages of climate adaptation. We hypothesize that climate adaptation in reality may require a greater preponderance of transformative strategies, and that scientists and institutions should accelerate exploring such approaches to define and develop the next generation of conservation strategies. Acknowledgements We would like to express Selleckchem TGFbeta inhibitor appreciation

to everyone involved in the buy Erismodegib climate adaptation process and workshop, especially the 20 conservation project teams and class facilitators and knowledge managers who contributed their ideas and experiences to the collective wisdom presented in this paper. Appreciation goes to Kristin Richards Betz and Anne Wallach Thomas for building and maintaining TNC’s climate adaptation website (http://​conserveonline.​org/​workspaces/​climateadaptatio​n) before, during, and after the workshop. This paper also benefitted from the input of our colleagues Peter Kareiva, Stacey Solie, Karen Lombard, and Dan Majka, and all the participants of the January 2010 TNC writing workshop in Tucson, ADP ribosylation factor Arizona. Open Access This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium,

provided the original author(s) and source are credited. Electronic supplementary material Below is the link to the electronic supplementary material. Supplementary material 1 (PDF 324 kb) Supplementary material 2 (PDF 182 kb) References Araujo MB, Rahbek C (2006) How does climate change affect biodiversity? Science 313:1396–1397PubMedCrossRef Bierwagen BG, Thomas R, Kane A (2008) Capacity of management plans for aquatic invasive species to integrate climate change. Conserv Biol 22:568–574PubMedCrossRef CMP (2007) Open standards for the practice of conservation. Version 2.0. http://​www.​conservationmeas​ures.​org/​CMP/​Site_​Docs/​CMP_​Open_​Standards_​Version_​2.​0.​pdf. Cited 22 Apr 2010 Dunwiddie PW, Hall SA, Ingraham MW, Bakker JD, Nelson KS, Fuller R, Gray E (2009) Rethinking conservation practice in light of climate change. Ecol Restor 27:320–329CrossRef Galatowitsch S, Frelich L, Phillips-Mao L (2009) Regional climate change adaptation strategies for biodiversity conservation in a midcontinental region of North America. Biol Conserv 142:2012–2022CrossRef Girvetz EH, Zganjar C, Raber GT, Maurer EP, Kareiva P, Lawler JJ (2009) Applied climate-change analysis: the Climate Wizard tool.

There were only five association rules that involved epitopes fro

There were only five association rules that involved epitopes from the Env gene. Four of these five were from Gag-Env

and one from Pol-Env associations. Notably, associations with antibody epitopes were limited to these five Env association rules, which can partially be attributed to the high degree of sequence divergence among the Env FK228 mouse sequences that can differ by as much as 30% I-BET151 at the amino acid level [76]. Figure 2 Relative composition of unique association rules involving multiple genes ( Gag , Pol and Nef ) and epitope types (Cytotoxic T Lymphocyte (CTL), T-Helper (Th) and antibody (Ab) epitopes). The 6142 unique association rules are classified according to the genes that harbor these epitopes. The pie-chart inside each segment represents the division according to the epitope region types involved. The single association rule in Nef-only category involved CTL and Th epitopes, while that in Pol-Env category involved CTL and Ab epitopes. Out of four association rules involving epitopes from Gag and Env, three belonged to CTL-Ab and one belonged to Th-Ab epitope regions types. No association rules included all three types of epitopes (CTL, Th and Ab) and

four genes (Gag, Pol, Env and Nef). MAPK inhibitor However, several “”multi-type”" association rules comprised of two different epitope types (CTL and Th) and three genes (Gag, Pol and Nef) were discovered (Figure 1, Additional file 5). For example, in the association rule: GHQAAMQML (CTL, Gag) – PKEPFRDYV (Th, Gag) – KLNWASQIY (CTL, Pol) – FLKEKGGL (CTL, Nef) (Figure 1), GHQAAMQML, KLNWASQIY and FLKEKGGL are CTL epitopes from the Gag, Pol and Nef genes, respectively, while PKEPFRDYV is a Th epitope from the Gag gene. Overall, there were 137 “”multi-type”" associations involving

epitopes from two types and three genes (2T-3G) among a total of 21 CTL and Th epitopes from the Gag, Pol and Nef genes (Additional Abiraterone file 5). These 21 epitopes can be mapped to 14 different non-overlapping genomic regions (Table 3) and a single association rule is generally spread across 3 to 5 of such regions. Interestingly, even though the association rule with the maximum number of epitopes in a single rule (9 epitopes) involved four non-overlapping genomic regions, it included epitopes from only two genes, Gag and Pol. Epitope-associations in the reference genome are representative of the global HIV-1 population Presence of association rules discovered in the reference genome set was verified by analyzing a larger worldwide set of 978 HIV-1 genomes (including 888 sequences from the 2008 web alignment and 90 reference sequences from the HIV Sequence database). The Gag, Pol and Nef genes in each sequence were concatenated for the purpose of the analysis, and presence of each association rule (as a complete match of all epitope regions involved) was noted. The results showed that most of the epitope-associations were present in the majority of genomes from the global HIV-1 population.

, 2009 [34] 3The number of reads

, 2009 [34]. 3The number of reads Selleck Autophagy inhibitor per dataset after removal of sequences that could be from the same source as those in the contamination control dataset. 4OTUs: Operational Taxonomic Units at 3% or 6% nucleotide difference. 5Number of phyla

and genera are based on taxonomic classification by MEGAN V3.4 [36, 37], with the total number of phyla and genera detected in parenthesis. 6Chao1 is an estimator of the minimum richness and is based on the number of rare OTUs (singletons and doublets) within a sample. 7The Shannon index combines estimates of richness (total number of OTUs) and evenness (relative abundance). 8The Shannon index after normalization of the number of sequences (as described in Methods). The 454 pyrosequencing method has a characteristic error rate in the form of insertion/deletion errors at homopolymer runs. To correct for this phenomenon, the raw reads were processed with PyroNoise [34] with a minimum length cutoff of 218 and 235 nt for the V1V2 and V6 regions, respectively. The PyroNoise program clusters all reads find more whose flowgrams indicate that they could stem from the same sequence, while also considering read abundance. After denoising, one sequence per cluster

together with the number of reads mapping to that cluster is reported. Next, the sequences (at this stage one sequence per denoised cluster) that did not have

an exact match to the primer were removed, and the forward primer sequence itself was also trimmed. Finally, the urine sample sequence sets were stripped for sequences that could be from the same source as those in the contamination control dataset. This was done by using Oxymatrine the program selleck chemicals ESPRIT http://​www.​biotech.​ufl.​edu/​people/​sun/​esprit.​html[35] to do a complete linkage clustering at 1% genetic difference of each sample together with its respective control. Before clustering, the control sequences were weighed so that there were the same number of reads stemming from both the sample and the control going into the process. Within each cluster the frequency of sample vs control sequence was calculated, and any sample sequences found in clusters where 50% or more of the sequences belonged to the control were removed. For taxonomic grouping we used MEGAN V3.4 http://​www-ab.​informatik.​uni-tuebingen.​de/​software/​megan/​welcome.​html[36, 37], which uses blast hits to place reads onto a taxonomy by assigning each read to a taxonomic group at a level in the NCBI taxonomy. The sequence reads (one read per denoised cluster from the pyronoise step) that passed the filtering steps were compared to a curated version of the SSUrdp database [38] using blastn with parameters set to a maximum expectation value (E) of 10-5. The 25 best hits were kept.