Here we describe, the longest microbial time-series analyzed to date using high-resolution 16S rRNA tag pyrosequencing of samples taken monthly over 6 years at a temperate marine coastal site off Plymouth, UK. Data treatment effected the estimation of community richness over a 6-year period, whereby 8794 operational taxonomic units (OTUs) were identified using single-linkage preclustering and 21 130 OTUs were identified by denoising the data. The Alphaproteobacteria were the most abundant Class, and the most frequently recorded OTUs were members of the Rickettsiales (SAR 11) and Rhodobacteriales. This near-surface ocean bacterial community showed strong repeatable seasonal patterns, which were defined by winter peaks in diversity across all years. Environmental variables explained far more variation in seasonally predictable bacteria than did data on protists or metazoan biomass. Change in day length alone explains >65% of the variance in community diversity. The results suggested that seasonal changes in environmental variables are more important than trophic interactions. Interestingly, microbial association network analysis showed that correlations in abundance were stronger within bacterial taxa rather than between bacteria and eukaryotes, or between bacteria and environmental variables.

No abstract

Running title: Marine bacterial seasonal succession 2 SummaryVery few marine microbial communities are well characterized even with the weight of research effort presently devoted to it. Only a small proportion of this effort has been aimed at investigating temporal community structure. Here we present the first report of the application of high-throughput pyrosequencing to investigate intraannual bacterial community structure. Microbial diversity was determined for 12 time points at the surface of the L4 sampling site in the Western English Channel. This was were found at all time-points, yet these few lineages comprised 50% of the total reads sequenced. The most abundant phylum was Proteobacteria (50% of all sequenced reads), while the SAR11 clade comprised 21% of the ubiquitous reads and ~12 % of the total sequenced reads. In contrast, 78% of all OTUs were only found at one timepoint and 67% were only found once, evidence of a large and transient rare assemblage. This time-series shows evidence of seasonally structured community diversity. There is also evidence for seasonal succession, primarily reflecting changes among dominant taxa. These changes in structure were significantly correlated to a combination of temperature, phosphate and silicate concentrations.3

Microbial communities inhabiting a multipond solar saltern were analysed and compared using SSU rRNA polymerase chain reaction (PCR)-based fingerprintings carried out in parallel by four laboratories. A salinity gradient from seawater (3.7%) to NaCl precipitation (37%) was studied for Bacteria, Archaea and Eukarya, and laboratories applied their own techniques and protocols on the same set of samples. Members of all three domains were retrieved from all salt concentrations. Three fingerprinting techniques were used: denaturing gradient gel electrophoresis (DGGE), ribosomal internal spacer analysis (RISA), and terminal-restriction fragments length polymorphism (T-RFLP). In addition, each laboratory used its own biomass collection method and DNA extraction protocols. Prokaryotes were addressed using DGGE and RISA with different 'domain-specific' primers sets. Eukaryotes were analysed by one laboratory using DGGE and T-RFLP, but targeting the same 18S rDNA site. Fingerprints were compared through cluster analysis and non-metric multidimensional scaling plots. This exercise allowed fast comparison of microbial assemblages and determined to what extent the picture provided by each laboratory was similar to those of others. Formation of two main, salinity-based groups of samples in prokaryotes (4-15% and 22-37% salinity) was consistent for all the laboratories. When other clusters appeared, this was a result of the particular technique and the protocol used in each case, but more affected by the primers set used. Eukaryotic microorganisms changed more from pond to pond; 4-5% and 8-37% salinity were but the two main groups detected. Archaea showed the lowest number of bands whereas Eukarya showed the highest number of operational taxonomic units (OTUs) in the initial ponds. Artefacts appeared in the DGGE from ponds with extremely low microbial richness. On the other hand, different 16S rDNA fragments with the same restriction or internal transcribed spacer (ITS) length were the main limitations for T-RFLP and RISA analyses, respectively, in ponds with the highest OTUs richness. However, although the particular taxonomic composition could vary among protocols, the general structure of the microbial assemblages was maintained.

Bacterial and archaeal assemblages have been studied in a multipond solar saltern using a range of microbial ecology techniques by four laboratories simultaneously. These include 16S rDNA sequencing from both denaturing gradient gel electrophoresis (DGGE) and clone libraries, and culturing methods. Water samples from eight ponds were analysed, covering a salinity range from near sea water (4% salt) to saturated sodium chloride (37% salt; ponds called crystallizers). Clone libraries focused on ponds with salinity of 8%, 22% and 32%. Although different cloning strategies were able to retrieve the same type of dominant sequences, there were differing degrees of success with less abundant sequences. Thus, the use of two sets of primers recovered a higher number of phylotypes. Bacterial and archaeal isolates were, however, different from any of the retrieved environmental sequences. For Bacteria, most sequences in the 8% salt pond were related to organisms of marine origin. Thus, representatives of the alpha-, beta-, gamma- and epsilon-subdivisions of Proteobacteria, the Cytophaga-Flavobacterium-Bacteroides group (CFB), high-G+C Gram-positive bacteria and cyanobacteria were found. In the 22% salt pond, alpha- and gamma-Proteobacteria, cyanobacteria and CFB were the only groups found, and most of them were related to specialized halophilic bacteria. From the 32% salt pond, only CFB were found, and most of the sequences retrieved clustered with Salinibacter ruber, an extremely halophilic bacterium. A decrease in the richness of bacterial genera was therefore apparent along the gradient. Archaea behaved quite similarly. In the lowest salinity ponds, sequences were related to environmental clones of Marine Archaea Group III (Thermoplasmales relatives) and to unclassified branches of Euryarchaeaota. In the 8%, 22% and 32% ponds, most of the clones were related to different cultured strains of Halobacteriaceae. Finally, most sequences from the crystallizers clustered with the uncultured square archaeon SPhT. Crenarchaeaota were not detected. Despite the fact that higher prokaryotic richness was apparent in the lower salinity ponds than in the crystallizers, the diversity index from clone libraries calculated according to Shannon and Weaver did not show this trend. This was because diversity in the crystallizers can be considered as 'microdiversity', the co-existence of several closely related clones of Bacteria (the S. ruber cluster) and Archaea (the SPhT cluster). Regardless of the changes in abundance, both Bacteria and Archaea showed the same pattern; as salinity increased, the number of different clusters decreased, and only one cluster became dominant. Both clusters, however, showed a considerable degree of microdiversity. The meaning of such microdiversity remains to be determined.