Diversity of zooplankton assemblages at Station L4, Western English Channel, measured using next generation DNA sequencing

Background: Zooplankton play an important role in our oceans, in biogeochemical cycling and providing a food source for commercially important fish larvae. However, difficulties in correctly identifying zooplankton hinder our understanding of their roles in marine ecosystem functioning, and can prevent detection of long term changes in their community structure. The advent of massively parallel Next Generation Sequencing technology allows DNA sequence data to be recovered directly from whole community samples. Here we assess the ability of such sequencing to quantify the richness and diversity of a mixed zooplankton assemblage from a productive monitoring site in the Western English Channel. Methodology/Principle Findings: Plankton WP2 replicate net hauls (200 µm) were taken at the Western Channel Observatory long-term monitoring station L4 in September 2010 and January 2011. These samples were analysed by microscopy and metagenetic analysis of the 18S nuclear small subunit ribosomal RNA gene using the 454 pyrosequencing platform. Following quality control a total of 419,042 sequences were obtained for all samples. The sequences clustered in to 205 operational taxonomic units using a 97% similarity cut-off. Allocation of taxonomy by comparison with the National Centre for Biotechnology Information database identified 138 OTUs to species level, 11 to genus level and 1 to order, <2.5% of sequences were classified as unknowns. By comparison a skilled microscopic analyst was able to routinely enumerate only 75 taxonomic groups. Conclusions: The percentage of OTUs assigned to major eukaryotic taxonomic groups broadly aligns between the metagenetic and morphological analysis and are dominated by Copepoda. However, the metagenetics reveals a previously hidden taxonomic richness, especially for Copepoda and meroplankton such as Bivalvia, Gastropoda and Polychaeta. It also reveals rare species and parasites. We conclude that Next Generation Sequencing of 18S amplicons is a powerful tool for estimating diversity and species richness of zooplankton communities.

(Table S2) Cell concentrations of individual and combined Alexandrium spp. and phylogenetic groups as obtained by qPCR assay and Utermöhl counts, as well as particulate PSP toxin concentrations from discrete water depths

Molecular methods provide promising tools for routine detection and quantification of toxic microalgae in plankton samples. To this end, novel TaqMan minor groove binding probes and primers targeting the small (SSU) or large (LSU) ribosomal subunit (rRNA) were developed for two species of the marine dinoflagellate genus Alexandrium (A. minutum, A. tamutum) and for three groups/ribotypes of the A. tamarense species complex: Group I/North American (NA), Group II/Mediterranean (ME) and Group III/Western European (WE). Primers and probes for real-time quantitative PCR (qPCR) were species-specific and highly efficient when tested in qPCR assays for cross-validation with pure DNA from cultured Alexandrium strains. Suitability of the qPCR assays as molecular tools for the detection and estimation of relative cell abundances of Alexandrium species and groups was evaluated from samples of natural plankton assemblages along the Scottish east coast. The results were compared with inverted microscope cell counts (Utermöhl technique) of Alexandrium spp. and associated paralytic shellfish poisoning (PSP) toxin concentrations. The qPCR assays indicated that A. tamarense (Group I) and A. tamutum were the most abundant Alexandrium taxa and both were highly positively correlated with PSP toxin content of plankton samples. Cells of A. tamarense (Group III) were present at nearly all stations but in low abundance. Alexandrium minutum and A. tamarense (Group II) cells were not detected in any of the samples, thereby arguing for their absence from the specific North Sea region, at least at the time of the survey. The sympatric occurrence of A. tamarense Group I and Group III gives further support to the hypothesis that the groups/ribotypes of the A. tamarense species complex are cryptic species rather than variants belonging to the same species.

Effects of increased CO2 on fish gill and plasma proteome

Ocean acidification and warming are both primarily caused by increased levels of atmospheric CO2, and marine organisms are exposed to these two stressors simultaneously. Although the effects of temperature on fish have been investigated over the last century, the long-term effects of moderate CO2 exposure and the combination of both stressors are almost entirely unknown. A proteomics approach was used to assess the adverse physiological and biochemical changes that may occur from the exposure to these two environmental stressors. We analysed gills and blood plasma of Atlantic halibut (Hippoglossus hippoglossus) exposed to temperatures of 12°C (control) and 18°C (impaired growth) in combination with control (400 µatm) or high-CO2 water (1000 µatm) for 14 weeks. The proteomic analysis was performed using two-dimensional gel electrophoresis (2DE) followed by Nanoflow LC-MS/MS using a LTQ-Orbitrap. The high-CO2 treatment induced the up-regulation of immune system-related proteins, as indicated by the up-regulation of the plasma proteins complement component C3 and fibrinogen beta chain precursor in both temperature treatments. Changes in gill proteome in the high-CO2 (18°C) group were mostly related to increased energy metabolism proteins (ATP synthase, malate dehydrogenase, malate dehydrogenase thermostable, and fructose-1,6-bisphosphate aldolase), possibly coupled to a higher energy demand. Gills from fish exposed to high-CO2 at both temperature treatments showed changes in proteins associated with increased cellular turnover and apoptosis signalling (annexin 5, eukaryotic translation elongation factor 1 gamma, receptor for protein kinase C, and putative ribosomal protein S27). This study indicates that moderate CO2-driven acidification, alone and combined with high temperature, can elicit biochemical changes that may affect fish health.

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.

Environmental conditions and habitats of GOS sites and bottle data of HOT-ALOHA station

The Global Ocean Sampling (GOS) expedition is currently the largest and geographically most comprehensive metagenomic dataset, including samples from the Atlantic, Pacific, and Indian Oceans. This study makes use of the wide range of environmental conditions and habitats encompassed within the GOS sites in order to investigate the ecological structuring of bacterial and archaeal taxon ranks. Community structures based on taxonomically classified 16S ribosomal RNA (rRNA) gene fragments at phylum, class, order, family, and genus rank levels were examined using multivariate statistical analysis, and the results were inspected in the context of oceanographic environmental variables and structured habitat classifications. At all taxon rank levels, community structures of neritic, oceanic, estuarine biomes, as well as other exotic biomes (salt marsh, lake, mangrove), were readily distinguishable from each other. A strong structuring of the communities with chlorophyll a concentration and a weaker yet significant structuring with temperature and salinity were observed. Furthermore, there were significant correlations between community structures and habitat classification. These results were used for further investigation of one-to-one relationships between taxa and environment and provided indications for ecological preferences shaped by primary production for both cultured and uncultured bacterial and archaeal clades.

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.

Experiment: Marine fungi may benefit from ocean acidification

Recent studies have discussed the consequences of ocean acidification for bacterial processes and diversity. However, the decomposition of complex substrates in marine environments, a key part of the flow of energy in ecosystems, is largely mediated by marine fungi. Although marine fungi have frequently been reported to prefer low pH levels, this group has been neglected in ocean acidification research. We present the first investigation of direct pH effects on marine fungal abundance and community structure. In microcosm experiments repeated in 2 consecutive years, we incubated natural North Sea water for 4 wk at in situ seawater pH (8.10 and 8.26), pH 7.82 and pH 7.67. Fungal abundance was determined by colony forming unit (cfu) counts, and fungal community structure was investigated by the culture-independent fingerprint method Fungal Automated Ribosomal Intergenic Spacer Analysis (F-ARISA). Furthermore, pH at the study site was determined over a yearly cycle. Fungal cfu were on average 9 times higher at pH 7.82 and 34 times higher at pH 7.67 compared to in situ seawater pH, and we observed fungal community shifts predominantly at pH 7.67. Currently, surface seawater pH at Helgoland Roads remains >8.0 throughout the year; thus we cannot exclude that fungal responses may differ in regions regularly experiencing lower pH values. However, our results suggest that under realistic levels of ocean acidification, marine fungi will reach greater importance in marine biogeochemical cycles. The rise of this group of organisms will affect a variety of biotic interactions in the sea.

Small Changes in pH Have Direct Effects on Marine Bacterial Community Composition: A Microcosm Approach

As the atmospheric CO2 concentration rises, more CO2 will dissolve in the oceans, leading to a reduction in pH. Effects of ocean acidification on bacterial communities have mainly been studied in biologically complex systems, in which indirect effects, mediated through food web interactions, come into play. These approaches come close to nature but suffer from low replication and neglect seasonality. To comprehensively investigate direct pH effects, we conducted highly-replicated laboratory acidification experiments with the natural bacterial community from Helgoland Roads (North Sea). Seasonal variability was accounted for by repeating the experiment four times (spring, summer, autumn, winter). Three dilution approaches were used to select for different ecological strategies, i.e. fast-growing or low-nutrient adapted bacteria. The pH levels investigated were in situ seawater pH (8.15-8.22), pH 7.82 and pH 7.67, representing the present-day situation and two acidification scenarios projected for the North Sea for the year 2100. In all seasons, both automated ribosomal intergenic spacer analysis and 16S ribosomal amplicon pyrosequencing revealed pH-dependent community shifts for two of the dilution approaches. Bacteria susceptible to changes in pH were different members of Gammaproteobacteria, Flavobacteriaceae, Rhodobacteraceae, Campylobacteraceae and further less abundant groups. Their specific response to reduced pH was often context-dependent. Bacterial abundance was not influenced by pH. Our findings suggest that already moderate changes in pH have the potential to cause compositional shifts, depending on the community assembly and environmental factors. By identifying pH-susceptible groups, this study provides insights for more directed, in-depth community analyses in large-scale and long-term experiments.

(Table T1) Total procaryotes, and living bacteria and archaea in black shales of ODP Leg 207 sites

Subseafloor sediments harbor over half of all prokaryotic cells on Earth (Whitman et al., 1998). This immense number is calculated from numerous microscopic acridine orange direct counts (AODCs) conducted on sediment cores drilled during the Ocean Drilling Program (ODP) (Parkes et al., 1994, doi:10.1038/371410a0, 2000, doi:10.1007/PL00010971). Because these counts cannot differentiate between living and inactive or even dead cells (Kepner and Pratt, 1994; Morita, 1997), the population size of living microorganisms has recently been enumerated for ODP Leg 201 sediment samples from the equatorial Pacific and the Peru margin using ribosomal ribonucleic acid targeting catalyzed reporter deposition-fluorescence in situ hybridization (CARD-FISH) (Schippers et al., 2005, doi:10.1038/nature03302). A large fraction of the subseafloor prokaryotes were alive, even in very old (16 Ma) and deep (>400 m) sediments. In this study, black shale samples from the Demerara Rise (Erbacher, Mosher, Malone, et al., 2004, doi:10.2973/odp.proc.ir.207.2004) were analyzed using AODC and CARD-FISH to find out if black shales also harbor microorganisms.