Development of bacterial and archaeal communities in erupted subsurface muds at the Håkon Mosby mud volcano

The microbial oxidation of methane controls the emission of the greenhouse gas methane from the ocean floor. However, some seabed structures such as mud volcanoes have leaky microbial methane filters and can be important sources of methane. We investigated the disturbance and recovery of a methanotrophic mud volcano microbiome (Håkon Mosby mud volcano, 1250 m water depth), to assess time scales of community succession and function in the natural deep-sea environment. We analyzed 10 surface and 5 subsurface sediment samples across HMMV mud flows from most recently discharged subsurface muds towards old consolidated muds as well as one reference site (REF) located approximately 0.5 km outside of the HMMV. Surface samples were obtained in 2003, 2009 and 2010. The surface of the new mud flows at the geographical center was sampled in 2009 and 2010. Around 100 m south of the center, we sampled more consolidated aged muds in 2003 and 2010. Old mud flows were sampled around 300 m southeast and 100 m north of the geographical center in 2003, 2009 and 2010. Surface sediment samples (0-20 cm) were recovered either by TV-guided Multicorer or by push cores using the remotely operated vehicle Quest (Marum, University Bremen). Subsurface sediments of all zones (>2 m below sea floor) were obtained in 2003 by gravity corer. After recovery, sediments were immediately subsampled in a refrigerated container (0°C) and further processed for biogeochemical analyses or preserved at -20°C for later DNA analyses. Our study show that freshly erupted muds hosted heterotrophic deep subsurface communities, which were replaced by surface communities within a few years of exposure. Aerobic methanotrophy was established at the top surface layer within less than a year, followed by anaerobic methanotrophy, sulfate reduction and finally thiotrophy. Our data indicate that it takes decades in cold environments before efficient methanotrophic communities establish to control methane emission. The observed succession provides insights to the response time of complex deep-sea communities to seafloor disturbances.

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.

Molecular fingerprinting and amplicon sequencing of the bacterial biofilm on settlement tiles deployed along natural pH gradients in Papua New Guinea

Bacterial biofilms provide cues for the settlement of marine invertebrates such as coral larvae, and are therefore important for the resilience and recovery of coral reefs. This study aimed to better understand how ocean acidification may affect the community composition and diversity of bacterial biofilms on surfaces under naturally reduced pH conditions. Settlement tiles were deployed at coral reefs in Papua New Guinea along pH gradients created by two CO2 seeps, and upper and lower tiles surfaces were sampled 5 and 13 months after deployment. Automated Ribosomal Intergenic Spacer Analysis were used to characterize more than 200 separate bacterial communities, complemented by amplicon sequencing of the bacterial 16S rRNA gene of 16 samples. The bacterial biofilm consisted predominantly of Alpha-, Gamma- and Deltaproteobacteria, as well as Cyanobacteria, Flavobacteriia and Cytophaga, whereas putative settlement-inducing taxa only accounted for a small fraction of the community. Bacterial biofilm composition was heterogeneous with approximately 25% shared operational taxonomic units between samples. Among the observed environmental parameters, pH only had a weak effect on community composition (R² ~ 1%) and did not affect community richness and evenness. In contrast, there were strong differences between upper and lower surfaces (contrasting in light exposure and grazing intensity). There also appeared to be a strong interaction between bacterial biofilm composition and the macroscopic components of the tile community. Our results suggest that on mature settlement surfaces in situ, pH does not have a strong impact on the composition of bacterial biofilms. Other abiotic and biotic factors such as light exposure and interactions with other organisms may be more important in shaping bacterial biofilms than changes in seawater pH.

Registry of all samples from the Tara Oceans Expedition (2009-2013)

The Tara Oceans Expedition (2009-2013) sampled the world oceans on board a 36 m long schooner, collecting environmental data and organisms from viruses to planktonic metazoans for later analyses using modern sequencing and state-of-the-art imaging technologies. Tara Oceans Data are particularly suited to study the genetic, morphological and functional diversity of plankton. Data sets in this collection provide methodological and environmental context to all samples collected during the Tara Oceans Expedition (2009-2013).

Table 2: Bacterial diversity in Ny-Ålesund, Spitzbergen

We describe the antibiotic resistance profiling of bacterial isolates collected from Ny-Alesund, Arctic, as part of the Indian Arctic Summer Expedition 2009. It was interesting to note that the bacterial isolates collected from the Arctic showed multidrug resistance. 32% of the isolates were found to be multi- drug resistant with several combinations of antibiotics. The 16S rRNA sequencing results shows a diverse group of bacteria belonging to Phyla Proteobacteria, Actinobacteria and Bacteriodetes and their relatedness was studied by phylogenetic analysis. While analysing the plasmid profiling, the most resistant two strains of Pseudomonas migulae showed multiple plasmids of varying sizes ~5.2-5.3 kb and ~9.5 kb. The extent and frequency of multidrug resistance in the polar bacteria deserves close monitoring and efforts to understand the various molecular mechanisms of drug resistance and control the spread of antibiotic resistance in polar environment is called for.