Recent advances in understanding the ecology of marine systems have been greatly facilitated by the growing availability of metagenomic data, which provide information on the identity, diversity and functional potential of the microbial community in a particular place and time. Here we present a dataset comprising over 5 terabases of metagenomic data from 610 samples spanning diverse regions of the Atlantic and Pacific Oceans. One set of metagenomes, collected on GEOTRACES cruises, captures large geographic transects at multiple depths per station. The second set represents two years of time-series data, collected at roughly monthly intervals from 3 depths at two long-term ocean sampling sites, Station ALOHA and BATS. These metagenomes contain genomic information from a diverse range of bacteria, archaea, eukaryotes and viruses. The data’s utility is strengthened by the availability of extensive physical, chemical, and biological measurements associated with each sample. We expect that these metagenomes will facilitate a wide range of comparative studies that seek to illuminate new aspects of marine microbial ecosystems.

Mesophilic ammonia-oxidizing Archaea (AOA) are abundant in a diverse range of marine environments, including the deep ocean, as revealed by the quantification of the archaeal amoA gene encoding the alpha-subunit of the ammonia monooxygenase. Using two different amoA primer sets, two distinct ecotypes of marine Crenarchaeota Group I (MCGI) were detected in the waters of the tropical Atlantic and the coastal Arctic. The HAC-AOA ecotype (high ammonia concentration AOA) was ≍ 8000 times and 15 times more abundant in the coastal Arctic and the top 300 m layer of the open equatorial Atlantic, respectively, than the LAC-AOA (low ammonia concentration AOA) ecotype. In contrast, the LAC-AOA ecotype dominated the lower meso- and bathypelagic waters of the tropical Atlantic (≍ 50 times more abundant than the HAC-AOA) where ammonia concentrations are well below the detection limit using conventional spectrophotometric or fluorometric methods. Cluster analysis of the sequences from the clone libraries obtained by the two amoA primer sets revealed two phylogenetically distinct clusters. Taken together, our results suggest the presence of two ecotypes of archaeal ammonia oxidizers corresponding to the medium (1.24 µM on average in the coastal Arctic) and low ammonia concentration (< 0.01 µM) in the shallow and the deep waters respectively.

Viruses are an abundant, diverse and dynamic component of marine ecosystems and have a key role in the biogeochemical processes of the ocean by controlling prokaryotic and phytoplankton abundance and diversity. However, most of the studies on virus–prokaryote interactions in marine environments have been performed in nearshore waters. To assess potential variations in the relation between viruses and prokaryotes in different oceanographic provinces, we determined viral and prokaryotic abundance and production throughout the water column along a latitudinal transect in the North Atlantic. Depth-related trends in prokaryotic and viral abundance (both decreasing by one order of magnitude from epi- to abyssopelagic waters), and prokaryotic production (decreasing by three orders of magnitude) were observed along the latitudinal transect. The virus-to-prokaryote ratio (VPR) increased from ∼19 in epipelagic to ∼53 in the bathy- and abyssopelagic waters. Although the lytic viral production decreased significantly with depth, the lysogenic viral production did not vary with depth. In bathypelagic waters, pronounced differences in prokaryotic and viral abundance were found among different oceanic provinces with lower leucine incorporation rates and higher VPRs in the North Atlantic Gyre province than in the provinces further north and south. The percentage of lysogeny increased from subpolar regions toward the more oligotrophic lower latitudes. Based on the observed trends over this latitudinal transect, we conclude that the viral–host interactions significantly change among different oceanic provinces in response to changes in the biotic and abiotic variables.

Prochlorococcus and Synechococcus are the dominant primary producers in marine ecosystems and perform a significant fraction of ocean carbon fixation. These cyanobacteria interact with a diverse microbial community that coexists with them. Comparative genomics of cultivated isolates has helped address questions regarding patterns of evolution and diversity among microbes, but the fraction that can be cultivated is miniscule compared to the diversity in the wild. To further probe the diversity of these groups and extend the utility of reference sequence databases, we report a data set of single cell genomes for 489 Prochlorococcus, 50 Synechococcus, 9 extracellular virus particles, and 190 additional microorganisms from a diverse range of bacterial, archaeal, and viral groups. Many of these uncultivated single cell genomes are derived from samples obtained on GEOTRACES cruises and at well-studied oceanographic stations, each with extensive suites of physical, chemical, and biological measurements. The genomic data reported here greatly increases the number of available Prochlorococcus genomes and will facilitate studies on evolutionary biology, microbial ecology, and biological oceanography.

Until recently, ammonia oxidation, a key process in the global nitrogen cycle, was thought to be mediated exclusively by a few bacterial groups. It has been shown now, that also Crenarchaeota are capable to perform this initial nitrification step. The abundance of ammonia oxidizing Bacteria and Archaea was determined using the bacterial and archaeal ammonia monooxygenase-a subunit (amoA) gene as functional markers in a quantitative PCR approach and related to the abundance of Bacteria and Archaea in the Eastern Mediterranean Sea. Archaeal amoA copy numbers decreased from 4000-5000 copies ml À1 seawater from the 200-500 m depth layer to 20 copies ml À1 at 1000 m depth. b-Proteobacterial amoA genes were below the detection limit in all the samples. The archaeal amoA copy numbers were correlated with NO 2 À concentrations, suggesting that ammonia-oxidizing Archaea may play a significant role in the nitrification in the mesopelagic waters of the Eastern Mediterranean Sea. In the bathypelagic waters, however, archaeal amoA gene abundance was rather low although Crenarchaeota were abundant, indicating that Crenarchaeota might largely lack the amoA gene in these deep waters. Terminal restriction fragment length polymorphism analysis of the archaeal community revealed a distinct clustering with the mesopelagic community distinctly different from the archaeal communities of both, the surface waters and the 3000-4000 m layers. Hence, the archaeal community in the Eastern Mediterranean Sea appears to be highly stratified despite the absence of major temperature and density gradients between the meso-and bathypelagic waters of the Mediterranean Sea.

In marine ecosystems, Thaumarchaeota are most likely the major ammonia oxidizers. While ammonia concentrations vary by about two orders of magnitude in the oceanic water column, archaeal ammonia oxidizers (AOA) vary by only one order of magnitude from surface to bathypelagic waters. Thus, the question arises whether the key enzyme responsible for ammonia oxidation, ammonia monooxygenase (amo), exhibits different affinities to ammonia along the oceanic water column and consequently, whether there are different ecotypes of AOA present in the oceanic water column. We determined the abundance and phylogeny of AOA based on their amoA gene. Two ecotypes of AOA exhibited a distribution pattern reflecting the reported availability of ammonia and the physico-chemical conditions throughout the Atlantic, and from epi- to bathypelagic waters. The distinction between these two ecotypes was not only detectable at the nucleotide level. Consistent changes were also detected at the amino acid level. These changes include substitutions of polar to hydrophobic amino acid, and glycine substitutions that could have an effect on the configuration of the amo protein and thus, on its activity. Although we cannot identify the specific effect, the ratio of non-synonymous to synonymous substitutions (dN/dS) between the two ecotypes indicates a strong positive selection between them. Consequently, our results point to a certain degree of environmental selection on these two ecotypes that have led to their niche specialization.

SummaryZooplankton and microbes play a key role in the ocean's biological cycles by releasing and consuming copious amounts of particulate and dissolved organic matter. Additionally, zooplankton provide a complex microhabitat rich in organic and inorganic nutrients in which bacteria thrive. In this study, we assessed the phylogenetic composition and metabolic potential of microbial communities associated with crustacean zooplankton species collected in the North Atlantic. Using Illumina sequencing of the 16S rRNA gene, we found significant differences between the microbial communities associated with zooplankton and those inhabiting the surrounding seawater. Metagenomic analysis of the zooplankton‐associated microbial community revealed a highly specialized bacterial community able to exploit zooplankton as microhabitat and thus, mediating biogeochemical processes generally underrepresented in the open ocean. The zooplankton‐associated bacterial community is able to colonize the zooplankton's internal and external surfaces using a large set of adhesion mechanisms and to metabolize complex organic compounds released or exuded by the zooplankton such as chitin, taurine and other complex molecules. Moreover, the high number of genes involved in iron and phosphorus metabolisms in the zooplankton‐associated microbiome suggests that this zooplankton‐associated bacterial community mediates specific biogeochemical processes (through the proliferation of specific taxa) that are generally underrepresented in the ambient waters.

Viral and prokaryotic abundance, production and diversity were determined throughout the water column of the subtropical Atlantic Ocean to assess potential variations in the relation between viruses and prokaryotes. Prokaryotic abundance and heterotrophic activity decreased by one and three orders of magnitude, respectively, from the epi-to the abyssopelagic layer. Although the lytic viral production (VP) decreased with depth, lysogenic VP was variable throughout the water column and did not show any trend with depth. The bacterial, archaeal and viral community composition were depth-stratified as determined by the automated ribosomal intergenic spacer analysis, terminalrestriction fragment length polymorphism and randomly amplified polymorphic DNA-PCR, respectively. Generally, the number of operational taxonomic units (OTUs) did not reveal consistent trends throughout the water column. Viral and prokaryotic abundance were strongly related to heterotrophic prokaryotic production, suggesting similar linkage strength between the viral and prokaryotic communities from the lower epi-to the abyssopelagic layer in the Atlantic Ocean. Strikingly, the prokaryotic and viral parameters exhibited a similar variability throughout the water column down to the abyssopelagic layers, suggesting that the dark ocean is as dynamic a system as is the lower epipelagic layer. It also indicates that viruses are apparently having a similar role for prokaryotic mortality in the dark oceanic realm as in surface waters. The more than twofold increase in bacterial OTUs from 2750 m depth to 45000 m depth and the concurrent decrease in viral OTUs, however, suggests that viruses might exhibit a wider host range in deep waters than in surface waters.