Comparative genomic and transcriptome analyses of pathotypes of Xanthomonas citri subsp. citri provide insights into mechanisms of bacterial virulence and host range

Abstract Background Citrus bacterial canker is a disease that has severe economic impact on citrus industries worldwide and is caused by a few species and pathotypes of Xanthomonas. X. citri subsp. citri strain 306 (XccA306) is a type A (Asiatic) strain with a wide host range, whereas its variant X. citri subsp. citri strain Aw12879 (Xcaw12879, Wellington strain) is restricted to Mexican lime. Results To characterize the mechanism for the differences in host range of XccA and Xcaw, the genome of Xcaw12879 that was completed recently was compared with XccA306 genome. Effectors xopAF and avrGf1 are present in Xcaw12879, but were absent in XccA306. AvrGf1 was shown previously for Xcaw to cause hypersensitive response in Duncan grapefruit. Mutation analysis of xopAF indicates that the gene contributes to Xcaw growth in Mexican lime but does not contribute to the limited host range of Xcaw. RNA-Seq analysis was conducted to compare the expression profiles of Xcaw12879 and XccA306 in Nutrient Broth (NB) medium and XVM2 medium, which induces hrp gene expression. Two hundred ninety two and 281 genes showed differential expression in XVM2 compared to in NB for XccA306 and Xcaw12879, respectively. Twenty-five type 3 secretion system genes were up-regulated in XVM2 for both XccA and Xcaw. Among the 4,370 common genes of Xcaw12879 compared to XccA306, 603 genes in NB and 450 genes in XVM2 conditions were differentially regulated. Xcaw12879 showed higher protease activity than XccA306 whereas Xcaw12879 showed lower pectate lyase activity in comparison to XccA306. Conclusions Comparative genomic analysis of XccA306 and Xcaw12879 identified strain specific genes. Our study indicated that AvrGf1 contributes to the host range limitation of Xcaw12879 whereas XopAF contributes to virulence. Transcriptome analyses of XccA306 and Xcaw12879 presented insights into the expression of the two closely related strains of X. citri subsp. citri. Virulence genes including genes encoding T3SS components and effectors are induced in XVM2 medium. Numerous genes with differential expression in Xcaw12879 and XccA306 were identified. This study provided the foundation to further characterize the mechanisms for virulence and host range of pathotypes of X. citri subsp. citri.

Antigen-presenting cells transfected with Hsp65 messenger RNA fail to treat experimental tuberculosis

In the last several years, the use of dendritic cells has been studied as a therapeutic strategy against tumors. Dendritic cells can be pulsed with peptides or full-length protein, or they can be transfected with DNA or RNA. However, comparative studies suggest that transfecting dendritic cells with messenger RNA (mRNA) is superior to other antigen-loading techniques in generating immunocompetent dendritic cells. In the present study, we evaluated a new therapeutic strategy to fight tuberculosis using dendritic cells and macrophages transfected with Hsp65 mRNA. First, we demonstrated that antigen-presenting cells transfected with Hsp65 mRNA exhibit a higher level of expression of co-stimulatory molecules, suggesting that Hsp65 mRNA has immunostimulatory properties. We also demonstrated that spleen cells obtained from animals immunized with mock and Hsp65 mRNA-transfected dendritic cells were able to generate a mixed Th1/Th2 response with production not only of IFN-γ but also of IL-5 and IL-10. In contrast, cells recovered from mice immunized with Hsp65 mRNA-transfected macrophages were able to produce only IL-5. When mice were infected with Mycobacterium tuberculosis and treated with antigen-presenting cells transfected with Hsp65 mRNA (therapeutic immunization), we did not detect any decrease in the lung bacterial load or any preservation of the lung parenchyma, indicating the inability of transfected cells to confer curative effects against tuberculosis. In spite of the lack of therapeutic efficacy, this study reports for the first time the use of antigen-presenting cells transfected with mRNA in experimental tuberculosis.

RNA Interference and cyclooxygenase-2 (COX-2) regulation in colon cancer cells

Despite new methods and combined strategies, conventional cancer chemotherapy still lacks specificity and induces drug resistance. Gene therapy can offer the potential to obtain the success in the clinical treatment of cancer and this can be achieved by replacing mutated tumour suppressor genes, inhibiting gene transcription, introducing new genes encoding for therapeutic products, or specifically silencing any given target gene. Concerning gene silencing, attention has recently shifted onto the RNA interference (RNAi) phenomenon. Gene silencing mediated by RNAi machinery is based on short RNA molecules, small interfering RNAs (siRNAs) and microRNAs (miRNAs), that are fully o partially homologous to the mRNA of the genes being silenced, respectively. On one hand, synthetic siRNAs appear as an important research tool to understand the function of a gene and the prospect of using siRNAs as potent and specific inhibitors of any target gene provides a new therapeutical approach for many untreatable diseases, particularly cancer. On the other hand, the discovery of the gene regulatory pathways mediated by miRNAs, offered to the research community new important perspectives for the comprehension of the physiological and, above all, the pathological mechanisms underlying the gene regulation. Indeed, changes in miRNAs expression have been identified in several types of neoplasia and it has also been proposed that the overexpression of genes in cancer cells may be due to the disruption of a control network in which relevant miRNA are implicated. For these reasons, I focused my research on a possible link between RNAi and the enzyme cyclooxygenase-2 (COX-2) in the field of colorectal cancer (CRC), since it has been established that the transition adenoma-adenocarcinoma and the progression of CRC depend on aberrant constitutive expression of COX-2 gene. In fact, overexpressed COX-2 is involved in the block of apoptosis, the stimulation of tumor-angiogenesis and promotes cell invasion, tumour growth and metastatization. On the basis of data reported in the literature, the first aim of my research was to develop an innovative and effective tool, based on the RNAi mechanism, able to silence strongly and specifically COX-2 expression in human colorectal cancer cell lines. In this study, I firstly show that an siRNA sequence directed against COX-2 mRNA (siCOX-2), potently downregulated COX-2 gene expression in human umbilical vein endothelial cells (HUVEC) and inhibited PMA-induced angiogenesis in vitro in a specific, non-toxic manner. Moreover, I found that the insertion of a specific cassette carrying anti-COX-2 shRNA sequence (shCOX-2, the precursor of siCOX-2 previously tested) into a viral vector (pSUPER.retro) greatly increased silencing potency in a colon cancer cell line (HT-29) without activating any interferon response. Phenotypically, COX-2 deficient HT-29 cells showed a significant impairment of their in vitro malignant behaviour. Thus, results reported here indicate an easy-to-use, powerful and high selective virus-based method to knockdown COX-2 gene in a stable and long-lasting manner, in colon cancer cells. Furthermore, they open up the possibility of an in vivo application of this anti-COX-2 retroviral vector, as therapeutic agent for human cancers overexpressing COX-2. In order to improve the tumour selectivity, pSUPER.retro vector was modified for the shCOX-2 expression cassette. The aim was to obtain a strong, specific transcription of shCOX-2 followed by COX-2 silencing mediated by siCOX-2 only in cancer cells. For this reason, H1 promoter in basic pSUPER.retro vector [pS(H1)] was substituted with the human Cox-2 promoter [pS(COX2)] and with a promoter containing repeated copies of the TCF binding element (TBE) [pS(TBE)]. These promoters were choosen because they are partculary activated in colon cancer cells. COX-2 was effectively silenced in HT-29 and HCA-7 colon cancer cells by using enhanced pS(COX2) and pS(TBE) vectors. In particular, an higher siCOX-2 production followed by a stronger inhibition of Cox-2 gene were achieved by using pS(TBE) vector, that represents not only the most effective, but also the most specific system to downregulate COX-2 in colon cancer cells. Because of the many limits that a retroviral therapy could have in a possible in vivo treatment of CRC, the next goal was to render the enhanced RNAi-mediate COX-2 silencing more suitable for this kind of application. Xiang and et al. (2006) demonstrated that it is possible to induce RNAi in mammalian cells after infection with engineered E. Coli strains expressing Inv and HlyA genes, which encode for two bacterial factors needed for successful transfer of shRNA in mammalian cells. This system, called “trans-kingdom” RNAi (tkRNAi) could represent an optimal approach for the treatment of colorectal cancer, since E. Coli in normally resident in human intestinal flora and could easily vehicled to the tumor tissue. For this reason, I tested the improved COX-2 silencing mediated by pS(COX2) and pS(TBE) vectors by using tkRNAi system. Results obtained in HT-29 and HCA-7 cell lines were in high agreement with data previously collected after the transfection of pS(COX2) and pS(TBE) vectors in the same cell lines. These findings suggest that tkRNAi system for COX-2 silencing, in particular mediated by pS(TBE) vector, could represent a promising tool for the treatment of colorectal cancer. Flanking the studies addressed to the setting-up of a RNAi-mediated therapeutical strategy, I proposed to get ahead with the comprehension of new molecular basis of human colorectal cancer. In particular, it is known that components of the miRNA/RNAi pathway may be altered during the progressive development of colorectal cancer (CRC), and it has been already demonstrated that some miRNAs work as tumor suppressors or oncomiRs in colon cancer. Thus, my hypothesis was that overexpressed COX-2 protein in colon cancer could be the result of decreased levels of one or more tumor suppressor miRNAs. In this thesis, I clearly show an inverse correlation between COX-2 expression and the human miR- 101(1) levels in colon cancer cell lines, tissues and metastases. I also demonstrate that the in vitro modulating of miR-101(1) expression in colon cancer cell lines leads to significant variations in COX-2 expression, and this phenomenon is based on a direct interaction between miR-101(1) and COX-2 mRNA. Moreover, I started to investigate miR-101(1) regulation in the hypoxic environment since adaptation to hypoxia is critical for tumor cell growth and survival and it is known that COX-2 can be induced directly by hypoxia-inducible factor 1 (HIF-1). Surprisingly, I observed that COX-2 overexpression induced by hypoxia is always coupled to a significant decrease of miR-101(1) levels in colon cancer cell lines, suggesting that miR-101(1) regulation could be involved in the adaption of cancer cells to the hypoxic environment that strongly characterize CRC tissues.

Generation and analysis of a mouse intestinal metatranscriptome through Illumina based RNA-sequencing

With the advent of high through-put sequencing (HTS), the emerging science of metagenomics is transforming our understanding of the relationships of microbial communities with their environments. While metagenomics aims to catalogue the genes present in a sample through assessing which genes are actively expressed, metatranscriptomics can provide a mechanistic understanding of community inter-relationships. To achieve these goals, several challenges need to be addressed from sample preparation to sequence processing, statistical analysis and functional annotation. Here we use an inbred non-obese diabetic (NOD) mouse model in which germ-free animals were colonized with a defined mixture of eight commensal bacteria, to explore methods of RNA extraction and to develop a pipeline for the generation and analysis of metatranscriptomic data. Applying the Illumina HTS platform, we sequenced 12 NOD cecal samples prepared using multiple RNA-extraction protocols. The absence of a complete set of reference genomes necessitated a peptide-based search strategy. Up to 16% of sequence reads could be matched to a known bacterial gene. Phylogenetic analysis of the mapped ORFs revealed a distribution consistent with ribosomal RNA, the majority from Bacteroides or Clostridium species. To place these HTS data within a systems context, we mapped the relative abundance of corresponding Escherichia coli homologs onto metabolic and protein-protein interaction networks. These maps identified bacterial processes with components that were well-represented in the datasets. In summary this study highlights the potential of exploiting the economy of HTS platforms for metatranscriptomics.

Intestinal bacterial translocation in rats with cirrhosis is related to compromised Paneth cell antimicrobial host defense

Liver cirrhosis is associated with bacterial translocation (BT) and endotoxemia. Most translocating bacteria belong to the common intestinal microbiota, suggesting a breakdown of intestinal barrier function. We hypothesized that diminished mucosal antimicrobial host defense could predispose to BT. Two rodent models of portal hypertension with increased BT were used, CCl(4)-induced ascitic cirrhosis and 2-day portal vein-ligated (PVL) animals. BT was assessed by standard microbiological techniques on mesenteric lymph nodes. Total RNA was isolated systematically throughout the intestinal tract, and expression of Paneth cell α-cryptdins and β-defensins was determined by real-time quantitative polymerase chain reaction (qPCR). To determine functional consequences, mucosal antimicrobial activity was assessed with a fluorescence-activated cell sorting assay. BT was detectable in 40% of rats with cirrhosis. Compared with the group without BT, these animals exhibited diminished intestinal Paneth cell α-cryptdin 5 and 7 expression. In contrast, PVL was associated with BT in all animals but did not affect antimicrobial peptides. The decrease in Paneth cell antimicrobials was most pronounced in the ileum and the coecum. Other antimicrobials showed no changes or even an induction in the case of BT at different sites. Antimicrobial activity toward different commensal strains was reduced, especially in the distal ileum and the cecum in experimental cirrhosis with BT (excluding PVL). Conclusion: Compromised Paneth cell antimicrobial host defense seems to predispose to BT in experimental cirrhosis. Understanding this liver-gut axis including the underlying mechanisms could help us to find new treatment avenues.

Comparison of amplification methods for transcriptomic analyses of low abundance prokaryotic RNA sources

Microarrays have established as instrumental for bacterial detection, identification, and genotyping as well as for transcriptomic studies. For gene expression analyses using limited numbers of bacteria (derived from in vivo or ex vivo origin, for example), RNA amplification is often required prior to labeling and hybridization onto microarrays. Evaluation of the fidelity of the amplification methods is crucial for the robustness and reproducibility of microarray results. We report here the first utilization of random primers and the highly processive Phi29 phage polymerase to amplify material for transcription profiling analyses. We compared two commercial amplification methods (GenomiPhi and MessageAmp kits) with direct reverse-transcription as the reference method, focusing on the robustness of mRNA quantification using either microarrays or quantitative RT-PCR. Both amplification methods using either poly-A tailing followed by in vitro transcription, or direct strand displacement polymerase, showed appreciable linearity. Strand displacement technique was particularly affordable compared to in vitro transcription-based (IVT) amplification methods and consisted in a single tube reaction leading to high amplification yields. Real-time measurements using low-, medium-, and highly expressed genes revealed that this simple method provided linear amplification with equivalent results in terms of relative messenger abundance as those obtained by conventional direct reverse-transcription.

Routine phenotypic identification of bacterial species of the family Pasteurellaceae isolated from animals

Pasteurellaceae are bacteria with an important role as primary or opportunistic, mainly respiratory, pathogens in domestic and wild animals. Some species of Pasteurellaceae cause severe diseases with high economic losses in commercial animal husbandry and are of great diagnostic concern. Because of new data on the phylogeny of Pasteurellaceae, their taxonomy has recently been revised profoundly, thus requiring an improved phenotypic differentiation procedure to identify the individual species of this family. A new and simplified procedure to identify species of Actinobacillus, Avibacterium, Gallibacterium, Haemophilus, Mannheimia, Nicoletella, and Pasteurella, which are most commonly isolated from clinical samples of diseased animals in veterinary diagnostic laboratories, is presented in the current study. The identification procedure was evaluated with 40 type and reference strains and with 267 strains from routine diagnostic analysis of various animal species, including 28 different bacterial species. Type, reference, and field strains were analyzed by 16S ribosomal RNA (rrs) and rpoB gene sequencing for unambiguous species determination as a basis to evaluate the phenotypic differentiation schema. Primary phenotypic differentiation is based on beta-nicotinamide adenine dinucleotide (beta-NAD) dependence and hemolysis, which are readily determined on the isolation medium. The procedure divides the 28 species into 4 groups for which particular biochemical reactions were chosen to identify the bacterial species. The phenotypic identification procedure allowed researchers to determine the species of 240 out of 267 field strains. The procedure is an easy and cost-effective system for the rapid identification of species of the Pasteurellaceae family isolated from clinical specimens of animals.

Gingival fluid cytokine expression and subgingival bacterial counts during pregnancy and postpartum: a case series

OBJECTIVES The aim of this study was to assess gingival fluid (GCF) cytokine messenger RNA (mRNA) levels, subgingival bacteria, and clinical periodontal conditions during a normal pregnancy to postpartum. MATERIALS AND METHODS Subgingival bacterial samples were analyzed with the checkerboard DNA-DNA hybridization method. GCF samples were assessed with real-time PCR including five proinflammatory cytokines and secretory leukocyte protease inhibitor. RESULTS Nineteen pregnant women with a mean age of 32 years (S.D. ± 4 years, range 26-42) participated in the study. Full-mouth bleeding scores (BOP) decreased from an average of 41.2% (S.D. ± 18.6%) at the 12th week of pregnancy to 26.6% (S.D. ± 14.4%) at the 4-6 weeks postpartum (p < 0.001). Between week 12 and 4-6 weeks postpartum, the mean probing pocket depth changed from 2.4 mm (S.D. ± 0.4) to 2.3 mm (S.D. ± 0.3) (p = 0.34). Higher counts of Eubacterium saburreum, Parvimonas micra, Selenomonas noxia, and Staphylococcus aureus were found at week 12 of pregnancy than at the 4-6 weeks postpartum examinations (p < 0.001). During and after pregnancy, statistically significant correlations between BOP scores and bacterial counts were observed. BOP scores and GCF levels of selected cytokines were not related to each other and no differences in GCF levels of the cytokines were observed between samples from the 12th week of pregnancy to 4-6 weeks postpartum. Decreasing postpartum counts of Porphyromonas endodontalis and Pseudomonas aeruginosa were associated with decreasing levels of Il-8 and Il-1β. CONCLUSIONS BOP decreased after pregnancy without any active periodontal therapy. Associations between bacterial counts and cytokine levels varied greatly in pregnant women with gingivitis and a normal pregnancy outcome. Postpartum associations between GCF cytokines and bacterial counts were more consistent. CLINICAL RELEVANCE Combined assessments of gingival fluid cytokines and subgingival bacteria may provide important information on host response.

A bacterial cytotoxin identifies the RhoA exchange factor Net1 as a key effector in the response to DNA damage.

BACKGROUND: Exposure of adherent cells to DNA damaging agents, such as the bacterial cytolethal distending toxin (CDT) or ionizing radiations (IR), activates the small GTPase RhoA, which promotes the formation of actin stress fibers and delays cell death. The signalling intermediates that regulate RhoA activation and promote cell survival are unknown. PRINCIPAL FINDINGS: We demonstrate that the nuclear RhoA-specific Guanine nucleotide Exchange Factor (GEF) Net1 becomes dephosphorylated at a critical inhibitory site in cells exposed to CDT or IR. Expression of a dominant negative Net1 or Net1 knock down by iRNA prevented RhoA activation, inhibited the formation of stress fibers, and enhanced cell death, indicating that Net1 activation is required for this RhoA-mediated responses to genotoxic stress. The Net1 and RhoA-dependent signals involved activation of the Mitogen-Activated Protein Kinase p38 and its downstream target MAPK-activated protein kinase 2. SIGNIFICANCE: Our data highlight the importance of Net1 in controlling RhoA and p38 MAPK mediated cell survival in cells exposed to DNA damaging agents and illustrate a molecular pathway whereby chronic exposure to a bacterial toxin may promote genomic instability.

Role of a ribosomal RNA phosphate oxygen during the EF-G–triggered GTP hydrolysis

Elongation factor-catalyzed GTP hydrolysis is a key reaction during the ribosomal elongation cycle. Recent crystal structures of G proteins, such as elongation factor G (EF-G) bound to the ribosome, as well as many biochemical studies, provide evidence that the direct interaction of translational GTPases (trGTPases) with the sarcin-ricin loop (SRL) of ribosomal RNA (rRNA) is pivotal for hydrolysis. However, the precise mechanism remains elusive and is intensively debated. Based on the close proximity of the phosphate oxygen of A2662 of the SRL to the supposedly catalytic histidine of EF-G (His87), we probed this interaction by an atomic mutagenesis approach. We individually replaced either of the two nonbridging phosphate oxygens at A2662 with a methyl group by the introduction of a methylphosphonate instead of the natural phosphate in fully functional, reconstituted bacterial ribosomes. Our major finding was that only one of the two resulting diastereomers, the SP methylphosphonate, was compatible with efficient GTPase activation on EF-G. The same trend was observed for a second trGTPase, namely EF4 (LepA). In addition, we provide evidence that the negative charge of the A2662 phosphate group must be retained for uncompromised activity in GTP hydrolysis. In summary, our data strongly corroborate that the nonbridging proSP phosphate oxygen at the A2662 of the SRL is critically involved in the activation of GTP hydrolysis. A mechanistic scenario is supported in which positioning of the catalytically active, protonated His87 through electrostatic interactions with the A2662 phosphate group and H-bond networks are key features of ribosome-triggered activation of trGTPases.

Role of a ribosomal RNA phosphate oxygen during the EF-G-triggered hydrolysis

Elongation factor-catalyzed GTP hydrolysis is a key reaction during the ribosomal elongation cycle. Recent crystal structures of G proteins, such as elongation factor G (EF-G) bound to the ribosome, as well as many biochemical studies, provide evidence that the direct interaction of translational GTPases (trGTPases) with the sarcin-ricin loop (SRL) of ribosomal RNA (rRNA) is pivotal for hydrolysis. However, the precise mechanism remains elusive and is intensively debated. Based on the close proximity of the phosphate oxygen of A2662 of the SRL to the supposedly catalytic histidine of EF-G (His87), we probed this interaction by an atomic mutagenesis approach. We individually replaced either of the two nonbridging phosphate oxygens at A2662 with a methyl group by the introduction of a methylphosphonate instead of the natural phosphate in fully functional, reconstituted bacterial ribosomes. Our major finding was that only one of the two resulting diastereomers, the SP methylphosphonate, was compatible with efficient GTPase activation on EF-G. The same trend was observed for a second trGTPase, namely EF4 (LepA). In addition, we provide evidence that the negative charge of the A2662 phosphate group must be retained for uncompromised activity in GTP hydrolysis. (1) In summary, our data strongly corroborate that the nonbridging proSP phosphate oxygen at the A2662 of the SRL is critically involved in the activation of GTP hydrolysis. A mechanistic scenario is supported in which positioning of the catalytically active, protonated His87 through electrostatic interactions with the A2662 phosphate group and H-bond networks are key features of ribosome-triggered activation of trGTPases.

Role of a ribosomal RNA phosphate oxygen during the EF-G-triggered hydrolysis

Elongation factor-catalyzed GTP hydrolysis is a key reaction during the ribosomal elongation cycle. Recent crystal structures of G proteins, such as elongation factor G (EF-G) bound to the ribosome, as well as many biochemical studies, provide evidence that the direct interaction of translational GTPases (trGTPases) with the sarcin-ricin loop (SRL) of ribosomal RNA (rRNA) is pivotal for hydrolysis. However, the precise mechanism remains elusive and is intensively debated. Based on the close proximity of the phosphate oxygen of A2662 of the SRL to the supposedly catalytic histidine of EF-G (His87), we probed this interaction by an atomic mutagenesis approach. We individually replaced either of the two nonbridging phosphate oxygens at A2662 with a methyl group by the introduction of a methylphosphonate instead of the natural phosphate in fully functional, reconstituted bacterial ribosomes. Our major finding was that only one of the two resulting diastereomers, the SP methylphosphonate, was compatible with efficient GTPase activation on EF-G. The same trend was observed for a second trGTPase, namely EF4 (LepA). In addition, we provide evidence that the negative charge of the A2662 phosphate group must be retained for uncompromised activity in GTP hydrolysis. (1) In summary, our data strongly corroborate that the nonbridging proSP phosphate oxygen at the A2662 of the SRL is critically involved in the activation of GTP hydrolysis. A mechanistic scenario is supported in which positioning of the catalytically active, protonated His87 through electrostatic interactions with the A2662 phosphate group and H-bond networks are key features of ribosome-triggered activation of trGTPases.

Histone-specific RNA 3' processing in nuclear extracts from mammalian cells.

The mature 3' ends of histone mRNAs are formed by endonucleolytic cleavage of longer precursor transcripts. This process occurs in the nucleus and can be regarded as the equivalent of the polyadenylation reaction involved in 3′-end-generation of all other mRNAs. A sea urchin H3 gene that failed to be properly processed in the Xenopus oocyte system proved particularly useful, because it allowed the identification of a processing component from sea urchins by a complementation assay. Nuclear extracts prepared from cells under various growth conditions have helped to reveal proliferation-dependent changes in the efficiency of histone RNA 3′ processing. RNA substrates for in vitro processing are best prepared by runoff transcription of specific DNA templates with bacterial or phage RNA polymerases. For this purpose, a restriction fragment containing the 3′-terminal region of a histone gene and including the conserved palindrome and spacer motifs is cloned into a polylinker sequence downstream of a strong promoter.

Effect of increased pCO2 on bacterial assemblage shifts in response to glucose addition in Fram Strait seawater mesocosms

Ocean acidification may stimulate primary production through increased availability of inorganic carbon in the photic zone, which may in turn change the biogenic flux of dissolved organic carbon (DOC) and the growth potential of heterotrophic bacteria. To investigate the effects of ocean acidification on marine bacterial assemblages, a two-by-three factorial mescosom experiment was conducted using surface sea water from the East Greenland Current in Fram Strait. Pyrosequencing of the V1-V2 region of bacterial 16S ribosomal RNA genes was used to investigate differences in the endpoint (Day 9) composition of bacterial assemblages in mineral nutrient-replete mesocosms amended with glucose (0 µm, 5.3 µm and 15.9 µm) under ambient (250 µatm) or acidified (400 µatm) partial pressures of CO2 (pCO2). All mesocosms showed low richness and diversity by Chao1 estimator and Shannon index, respectively, with general dominance by Gammaproteobacteria and Flavobacteria. Nonmetric multidimensional scaling analysis and two-way analysis of variance of the Jaccard dissimilarity matrix (97% similarity cut-off) demonstrated that the significant community shift between 0 µm and 15.9 µm glucose addition at 250 µatm pCO2 was eliminated at 400 µatm pCO2. These results suggest that the response potential of marine bacteria to DOC input may be altered under acidified conditions.

Metadata und statistic analysis of archael and bacterial sequences originating from sediments of the Håkon Mosby mud volcano (Z1-Z3, all habitats)

DNA extraction was carried out as described on the MICROBIS project pages (http://icomm.mbl.edu/microbis ) using a commercially available extraction kit. We amplified the hypervariable regions V4-V6 of archaeal and bacterial 16S rRNA genes using PCR and several sets of forward and reverse primers (http://vamps.mbl.edu/resources/primers.php). Massively parallel tag sequencing of the PCR products was carried out on a 454 Life Sciences GS FLX sequencer at Marine Biological Laboratory, Woods Hole, MA, following the same experimental conditions for all samples. Sequence reads were submitted to a rigorous quality control procedure based on mothur v30 (doi:10.1128/AEM.01541-09) including denoising of the flow grams using an algorithm based on PyroNoise (doi:10.1038/nmeth.1361), removal of PCR errors and a chimera check using uchime (doi:10.1093/bioinformatics/btr381). The reads were taxonomically assigned according to the SILVA taxonomy (SSURef v119, 07-2014; doi:10.1093/nar/gks1219) implemented in mothur and clustered at 98% ribosomal RNA gene V4-V6 sequence identity. V4-V6 amplicon sequence abundance tables were standardized to account for unequal sampling effort using 1000 (Archaea) and 2300 (Bacteria) randomly chosen sequences without replacement using mothur and then used to calculate inverse Simpson diversity indices and Chao1 richness (doi:10.2307/4615964). Bray-Curtis dissimilarities (doi:10.2307/1942268) between all samples were calculated and used for 2-dimensional non metric multidimensional scaling (NMDS) ordinations with 20 random starts (doi:10.1007/BF02289694). Stress values below 0.2 indicated that the multidimensional dataset was well represented by the 2D ordination. NMDS ordinations were compared and tested using Procrustes correlation analysis (doi:10.1007/BF02291478). All analyses were carried out with the R statistical environment and the packages vegan (available at: http://cran.r-project.org/package=vegan), labdsv (available at: http://cran.r-project.org/package=labdsv), as well as with custom R scripts. Operational taxonomic units at 98% sequence identity (OTU0.03) that occurred only once in the whole dataset were termed absolute single sequence OTUs (SSOabs; doi:10.1038/ismej.2011.132). OTU0.03 sequences that occurred only once in at least one sample, but may occur more often in other samples were termed relative single sequence OTUs (SSOrel). SSOrel are particularly interesting for community ecology, since they comprise rare organisms that might become abundant when conditions change.16S rRNA amplicons and metagenomic reads have been stored in the sequence read archive under SRA project accession number SRP042162.