Microbial community structure in the North Pacific ocean

Gradients of bacteria in the oceanic water column reveal finely-resolved vertical distributions

Bacterial communities directly influence ecological processes in the ocean, and depth has a major influence due to the changeover in primary energy sources between the sunlit photic zone and dark ocean. Here, we examine the abundance and diversity of bacteria in Monterey Bay depth profiles collected from the surface to just above the sediments (e.g., 2000 m). Bacterial abundance in these Pacific Ocean samples decreased by >1 order of magnitude, from 1.22 ±0.69 ×106 cells ml-1 in the variable photic zone to 1.44 ± 0.25 ×105 and 6.71 ± 1.23 ×104 cells ml-1 in the mesopelagic and bathypelagic, respectively. V1-V2 16S rRNA gene profiling showed diversity increased sharply between the photic and mesopelagic zones. Weighted Gene Correlation Network Analysis clustered co-occurring bacterial amplicon sequence variants (ASVs) into seven subnetwork modules, of which five strongly correlated with depth-related factors. Within surface-associated modules there was a clear distinction between a 'copiotrophic' module, correlating with chlorophyll and dominated by e.g., Flavobacteriales and Rhodobacteraceae, and an 'oligotrophic' module dominated by diverse Oceanospirillales (such as uncultured JL-ETNP-Y6, SAR86) and Pelagibacterales. Phylogenetic reconstructions of Pelagibacterales and SAR324 using full-length 16S rRNA gene data revealed several additional subclades, expanding known microdiversity within these abundant lineages, including new Pelagibacterales subclades Ia.B, Id, and IIc, which comprised 4-10% of amplicons depending on the subclade and depth zone. SAR324 and Oceanospirillales dominated in the mesopelagic, with SAR324 clade II exhibiting its highest relative abundances (17±4%) in the lower mesopelagic (300-750 m). The two newly-identified SAR324 clades showed highest relative abundances in the photic zone (clade III), while clade IV was extremely low in relative abundance, but present across dark ocean depths. Hierarchical clustering placed microbial communities from 900 m samples with those from the bathypelagic, where Marinimicrobia was distinctively relatively abundant. The patterns resolved herein, through high resolution and statistical replication, establish baselines for marine bacterial abundance and taxonomic distributions across the Monterey Bay water column, against which future change can be assessed.

Global biogeography of the smallest plankton across ocean depths

Tiny ocean plankton (picoplankton) are fundamental for the functioning of the biosphere, but the ecological mechanisms shaping their biogeography were partially understood. Comprehending whether these microorganisms are structured by niche versus neutral processes is relevant in the context of global change. We investigate the ecological processes (selection, dispersal, and drift) structuring global-ocean picoplanktonic communities inhabiting the epipelagic (0 to 200 meters), mesopelagic (200 to 1000 meters), and bathypelagic (1000 to 4000 meters) zones. We found that selection decreased, while dispersal limitation increased with depth, possibly due to differences in habitat heterogeneity and dispersal barriers such as water masses and bottom topography. Picoplankton β-diversity positively correlated with environmental heterogeneity and water mass variability, but this relationship tended to be weaker for eukaryotes than for prokaryotes. Community patterns were more pronounced in the Mediterranean Sea, probably because of its cross-basin environmental heterogeneity and deep-water isolation. We conclude that different combinations of ecological mechanisms shape the biogeography of the ocean microbiome across depths.

The Diversity and Metabolism of Culturable Nitrate-Reducing Bacteria from the Photic Zone of the Western North Pacific Ocean

To better understand bacterial communities and metabolism under nitrogen deficiency, 154 seawater samples were obtained from 5 to 200 m at 22 stations in the photic zone of the Western North Pacific Ocean. Total 634 nitrate-utilizing bacteria were isolated using selective media and culture-dependent methods, and 295 of them were positive for nitrate reduction. These nitrate-reducing bacteria belonged to 19 genera and 29 species and among them, Qipengyuania flava, Roseibium aggregatum, Erythrobacter aureus, Vibrio campbellii, and Stappia indica were identified from all tested seawater layers of the photic zone and at almost all stations. Twenty-nine nitrate-reducing strains representing different species were selected for further the study of nitrogen, sulfur, and carbon metabolism. All 29 nitrate-reducing isolates contained genes encoding dissimilatory nitrate reduction or assimilatory nitrate reduction. Six nitrate-reducing isolates can oxidize thiosulfate based on genomic analysis and activity testing, indicating that nitrate-reducing thiosulfate-oxidizing bacteria exist in the photic zone. Five nitrate-reducing isolates obtained near the chlorophyll a-maximum layer contained a dimethylsulfoniopropionate synthesis gene and three of them contained both dimethylsulfoniopropionate synthesis and cleavage genes. This suggests that nitrate-reducing isolates may participate in dimethylsulfoniopropionate synthesis and catabolism in photic seawater. The presence of multiple genes for chitin degradation and extracellular peptidases may indicate that almost all nitrate-reducing isolates (28/29) can use chitin and proteinaceous compounds as important sources of carbon and nitrogen. Collectively, these results reveal culturable nitrate-reducing bacterial diversity and have implications for understanding the role of such strains in the ecology and biogeochemical cycles of nitrogen, sulfur, and carbon in the oligotrophic marine photic zone.

Water Mass Controlled Vertical Stratification of Bacterial and Archaeal Communities in the Western Arctic Ocean During Summer Sea-Ice Melting

The environmental variations and their interactions with the biosphere are vital in the Arctic Ocean during the summer sea-ice melting period in the current scenario of climate change. Hence, we analysed the vertical distribution of bacterial and archaeal communities in the western Arctic Ocean from sea surface melt-ponds to deep water up to a 3040 m depth. The distribution of microbial communities showed a clear stratification with significant differences among different water depths, and the water masses in the Arctic Ocean - surface mixed layer, Atlantic water mass and deep Arctic water - appeared as a major factor explaining their distribution in the water column. A total of 34 bacterial phyla were detected in the seawater and 10 bacterial phyla in melt-ponds. Proteobacteria was the dominant phyla in the seawater irrespective of depth, whereas Bacteroidota was the dominant phyla in the melt-ponds. A fast expectation-maximization microbial source tracking analysis revealed that only limited dispersion of the bacterial community was possible across the stratified water column. The surface water mass contributed 21% of the microbial community to the deep chlorophyll maximum (DCM), while the DCM waters contributed only 3% of the microbial communities to the deeper water masses. Atlantic water mass contributed 37% to the microbial community of the deep Arctic water. Oligotrophic heterotrophic bacteria were dominant in the melt-ponds and surface waters, whereas chemoautotrophic and mixotrophic bacterial and archaeal communities were abundant in deeper waters. Chlorophyll and ammonium were the major environmental factors that determined the surface microbial communities, whereas inorganic nutrient concentrations controlled the deep-water communities.

Linkages and key factors between soil bacterial and fungal communities along an altitudinal gradient of different slopes on mount Segrila, Tibet, China

Soil microbes are of great significance to many energy flow and material circulation processes in alpine forest ecosystems. The distribution pattern of soil microbial community along altitudinal gradients is an essential research topic for the Tibetan Plateau. Yet our understanding of linkages between soil microbial communities and key factors along an altitudinal gradient of different slopes remains limited. Here, the diversity, composition and interaction of bacterial and fungal communities and in response to environmental factors were compared across five elevation sites (3,500 m, 3,700 m, 3,900 m, 4,100 m, 4,300 m) on the eastern and western slopes of Mount Segrila, by using Illumina MiSeq sequencing. Our results showed that microbial community composition and diversity were distinct at different elevations, being mainly influenced by soil total nitrogen and carbonate. Structural equation models indicated that elevation had a greater influence than slope upon the soil microbial community. Co-occurrence network analysis suggested that fungi were stable but bacteria contributed more to among interactions of bacterial and fungal communities. Ascomycota was identified as a key hub for the internal interactions of microbial community, which might affect the soil microbial co-occurrence network resilience of alpine forest ecosystems on the Tibetan Plateau.

Microbial communities on plastic particles in surface waters differ from subsurface waters of the North Pacific Subtropical Gyre

The long-term fate of plastics in the ocean and their interactions with marine microorganisms remain poorly understood. In particular, the role of sinking plastic particles as a transport vector for surface microbes towards the deep sea has not been investigated. Here, we present the first data on the composition of microbial communities on floating and suspended plastic particles recovered from the surface to the bathypelagic water column (0-2000 m water depth) of the North Pacific Subtropical Gyre. Microbial community composition of suspended plastic particles differed from that of plastic particles afloat at the sea surface. However, in both compartments, a diversity of hydrocarbon-degrading bacteria was identified. These findings indicate that microbial community members initially present on floating plastics are quickly replaced by microorganisms acquired from deeper water layers, thus suggesting a limited efficiency of sinking plastic particles to vertically transport microorganisms in the North Pacific Subtropical Gyre.

Moroccan Lagoon Microbiomes

Lagoons are fragile marine ecosystems that are considerably affected by anthropogenic pollutants. We performed a spatiotemporal characterization of the microbiome of two Moroccan lagoons, Marchica and Oualidia, both classified as Ramsar sites, the former on the Mediterranean coast and the latter on the Atlantic coast. We investigated their microbial diversity and abundance using 16S rRNA amplicon- and shotgun-based metagenomics approaches during the summers of 2014 and 2015. The bacterial microbiome was composed primarily of Proteobacteria (25–53%, 29–29%), Cyanobacteria (34–12%, 11–0.53%), Bacteroidetes (24–16%, 23–43%), Actinobacteria (7–11%, 13–7%), and Verrucomicrobia (4–1%, 15–14%) in Marchica and Oualidia in 2014 and 2015, respectively. Interestingly, 48 strains were newly reported in lagoon ecosystems, while eight unknown viruses were detected in Mediterranean Marchica only. Statistical analysis showed higher microbial diversity in the Atlantic lagoon than in the Mediterranean lagoon and a robust relationship between alpha diversity and geographic sampling locations. This first-ever metagenomics study on Moroccan aquatic ecosystems enriched the national catalog of marine microorganisms. They will be investigated as candidates for bioindication properties, biomonitoring potential, biotechnology valorization, biodiversity protection, and lagoon health assessment.

Contrasting diversity patterns of prokaryotes and protists over time and depth at the San-Pedro Ocean Time series

Community dynamics are central in microbial ecology, yet we lack studies comparing diversity patterns among marine protists and prokaryotes over depth and multiple years. Here, we characterized microbes at the San-Pedro Ocean Time series (2005-2018), using SSU rRNA gene sequencing from two size fractions (0.2-1 and 1-80 μm), with a universal primer set that amplifies from both prokaryotes and eukaryotes, allowing direct comparisons of diversity patterns in a single set of analyses. The 16S + 18S rRNA gene composition in the small size fraction was mostly prokaryotic (>92%) as expected, but the large size fraction unexpectedly contained 46-93% prokaryotic 16S rRNA genes. Prokaryotes and protists showed opposite vertical diversity patterns; prokaryotic diversity peaked at mid-depth, protistan diversity at the surface. Temporal beta-diversity patterns indicated prokaryote communities were much more stable than protists. Although the prokaryotic communities changed monthly, the average community stayed remarkably steady over 14 years, showing high resilience. Additionally, particle-associated prokaryotes were more diverse than smaller free-living ones, especially at deeper depths, contributed unexpectedly by abundant and diverse SAR11 clade II. Eukaryotic diversity was strongly correlated with the diversity of particle-associated prokaryotes but not free-living ones, reflecting that physical associations result in the strongest interactions, including symbioses, parasitism, and decomposer relationships.

Patterns of Relative Bacterial Richness and Community Composition in Seawater and Marine Sediment Are Robust for Both Operational Taxonomic Units and Amplicon Sequence Variants

To understand the relative influences of operational taxonomic units (OTUs) and amplicon sequence variants (ASVs) on patterns of marine microbial diversity and community composition, we examined bacterial diversity and community composition of seawater from 12 sites in the North Atlantic Ocean and Canadian Arctic and sediment from two sites in the North Atlantic. For the seawater analyses, we included samples from three to six zones in the water column of each site. For the sediment analyses, we included over 20 sediment horizons at each of two sites. For all samples, we amplified the V4-V5 hypervariable region of the 16S ribosomal RNA (rRNA) gene. We analyzed each sample in two different ways: (i) by clustering its reads into 97%-similar OTUs and (ii) by assigning sequences to unique ASVs. OTU richness is much higher than ASV richness for every sample, but both OTUs and ASVs exhibit similar vertical patterns of relative diversity in both the water column and the sediment. Bacterial richness is highest just below the photic zone in the water column and at the seafloor in the sediment. For both OTUs and ASVs, richness estimates depend on the number of sequences analyzed. Both methods yield broadly similar community compositions for each sample at the taxonomic levels of phyla to families. While the two methods yield different richness values, broad-scale patterns of relative richness and community composition are similar with both methods.

Biogeography of Southern Ocean prokaryotes: a comparison of the Indian and Pacific sectors

We investigated the Southern Ocean (SO) prokaryote community structure via zero‐radius operational taxonomic unit (zOTU) libraries generated from 16S rRNA gene sequencing of 223 full water column profiles. Samples reveal the prokaryote diversity trend between discrete water masses across multiple depths and latitudes in Indian (71–99°E, summer) and Pacific (170–174°W, autumn‐winter) sectors of the SO. At higher taxonomic levels (phylum‐family) we observed water masses to harbour distinct communities across both sectors, but observed sectorial variations at lower taxonomic levels (genus‐zOTU) and relative abundance shifts for key taxa such as Flavobacteria, SAR324/Marinimicrobia, Nitrosopumilus and Nitrosopelagicus at both epi‐ and bathy‐abyssopelagic water masses. Common surface bacteria were abundant in several deep‐water masses and vice‐versa suggesting connectivity between surface and deep‐water microbial assemblages. Bacteria from same‐sector Antarctic Bottom Water samples showed patchy, high beta‐diversity which did not correlate well with measured environmental parameters or geographical distance. Unconventional depth distribution patterns were observed for key archaeal groups: Crenarchaeota was found across all depths in the water column and persistent high relative abundances of common epipelagic archaeon Nitrosopelagicus was observed in deep‐water masses. Our findings reveal substantial regional variability of SO prokaryote assemblages that we argue should be considered in wide‐scale SO ecosystem microbial modelling.

Planktonic Archaeal Ether Lipid Origins in Surface Waters of the North Pacific Subtropical Gyre

Thaumarchaeota and Thermoplasmatota are the most abundant planktonic archaea in the sea. Thaumarchaeota contain tetraether lipids as their major membrane lipids, but the lipid composition of uncultured planktonic Thermoplasmatota representatives remains unknown. To address this knowledge gap, we quantified archaeal cells and ether lipids in open ocean depth profiles (0–200 m) of the North Pacific Subtropical Gyre. Planktonic archaeal community structure and ether lipid composition in the water column partitioned into two separate clusters: one above the deep chlorophyll maximum, the other within and below it. In surface waters, Thermoplasmatota densities ranged from 2.11 × 10⁶ to 6.02 × 10⁶ cells/L, while Thaumarchaeota were undetectable. As previously reported for Thaumarchaeota, potential homologs of archaeal tetraether ring synthases were present in planktonic Thermoplasmatota metagenomes. Despite the absence of Thaumarchaeota in surface waters, measurable amounts of intact polar ether lipids were found there. Based on cell abundance estimates, these surface water archaeal ether lipids contributed only 1.21 × 10^–9 ng lipid/Thermoplasmatota cell, about three orders of magnitude less than that reported for Thaumarchaeota cells. While these data indicate that even if some tetraether and diether lipids may be derived from Thermoplasmatota, they would only comprise a small fraction of Thermoplasmatota total biomass. Therefore, while both MGI Thaumarchaeota and MGII/III Thermoplasmatota are potential biological sources of archaeal GDGTs, the Thaumarchaeota appear to be the major contributors of archaeal tetraether lipids in planktonic marine habitats. These results extend and confirm previous reports of planktonic archaeal lipid sources, and further emphasize the need for Thermoplasmatota cultivation, to better characterize the membrane lipid constituents of marine planktonic Thermoplasmatota, and more precisely define the sources and patterns of archaeal tetraether lipid distributions in marine plankton.

Dynamic prokaryotic communities in the dark western Mediterranean Sea

Dark ocean microbial dynamics are fundamental to understand ecosystem metabolism and ocean biogeochemical processes. Yet, the ecological response of deep ocean communities to environmental perturbations remains largely unknown. Temporal and spatial dynamics of the meso- and bathypelagic prokaryotic communities were assessed throughout a 2-year seasonal sampling across the western Mediterranean Sea. A common pattern of prokaryotic communities' depth stratification was observed across the different regions and throughout the seasons. However, sporadic and drastic alterations of the community composition and diversity occurred either at specific water masses or throughout the aphotic zone and at a basin scale. Environmental changes resulted in a major increase in the abundance of rare or low abundant phylotypes and a profound change of the community composition. Our study evidences the temporal dynamism of dark ocean prokaryotic communities, exhibiting long periods of stability but also drastic changes, with implications in community metabolism and carbon fluxes. Taken together, the results highlight the importance of monitoring the temporal patterns of dark ocean prokaryotic communities.

Thermal discharge-induced seawater warming alters richness, community composition and interactions of bacterioplankton assemblages in a coastal ecosystem

Despite accumulating evidence on the impact of global climate warming on marine microbes, how increasing seawater temperature influences the marine bacterioplankton communities is elusive. As temperature gradient created by thermal discharges provides a suitable in situ model to study the influence of warming on marine microorganisms, surface seawater were sampled consecutively for one year (September-2016 to August-2017) from the control (unimpacted) and thermal discharge-impacted areas of a coastal power plant, located in India. The bacterioplankton community differences between control (n = 16) and thermal discharge-impacted (n = 26) areas, as investigated using 16S rRNA gene tag sequencing revealed reduced richness and varied community composition at thermal discharge-impacted areas. The relative proportion of Proteobacteria was found to be higher (average ~ 15%) while, Bacteroidetes was lower (average ~ 10%) at thermal discharge-impacted areas. Intriguingly, thermal discharge-impacted areas were overrepresented by several potential pathogenic bacterial genera (e.g. Pseudomonas, Acinetobacter, Sulfitobacter, Vibrio) and other native marine genera (e.g. Marinobacter, Pseudoalteromonas, Alteromonas, Pseudidiomarina, Halomonas). Further, co-occurrence networks demonstrated that complexity and connectivity of networks were altered in warming condition. Altogether, results indicated that increasing temperature has a profound impact on marine bacterioplankton richness, community composition, and inter-species interactions. Our findings are immensely important in forecasting the consequences of future climate changes especially, ocean warming on marine microbiota.

Temporal dynamics of carbon sequestration in coastal North Atlantic fjord system as seen through dissolved organic matter characterisation

Fjord systems in higher latitudes are unique coastal water ecosystems that facilitate the study of dissolved organic matter (DOM) dynamics from surface to deeper waters. The current work was undertaken in the Trondheim fjord characterized by North Atlantic waters, and compared DOM fractions from three depths - surface (3 m), intermediate (225 m) and deep (440 m) in four seasons, from late spring to winter in 2017. The high-resolution mass spectrometry data showed that DOM composition varies significantly in different seasons rather than in different depths in the fjord systems. The bacterial community composition was comparable except at spring surface and summer intermediate depths. Bacterial production was minimal below the euphotic layer, even with sufficient availability of inorganic nutrients. The bacterial production rate in the surface waters was about 7 times and over 50 times higher than that of the aphotic zone in the winter and the summer seasons, respectively. The surface heterotrophic microbial communities might have rapidly consumed the available labile DOM, with the production of more refractory DOM limiting bacterial production in aphotic layers. The greater number of CRAM-like formulas determined in the surface waters compared to other depths supports our hypothesis. The refractory DOM sequestered in the water column may either be exported into sediments attached to particulate matter and marine gels, or may escape into the atmosphere as carbon dioxide/monoxide during the photochemical oxidation pathways, suggesting that it is involved in climate change scenarios.

Microbial community transcriptional patterns vary in response to mesoscale forcing in the North Pacific Subtropical Gyre

The physical and biological dynamics that influence phytoplankton communities in the oligotrophic ocean are complex, changing across broad temporal and spatial scales. Eukaryotic phytoplankton (e.g., diatoms), despite their relatively low abundance in oligotrophic waters, are responsible for a large component of the organic matter flux to the ocean interior. Mesoscale eddies can impact both microbial community structure and function, enhancing primary production and carbon export, but the mechanisms that underpin these dynamics are still poorly understood. Here, mesoscale eddy influences on the taxonomic diversity and expressed functional profiles of surface communities of microeukaryotes and particle-associated heterotrophic bacteria from the North Pacific Subtropical Gyre were assessed over 2 years (spring 2016 and summer 2017). The taxonomic diversity of the microeukaryotes significantly differed by eddy polarity (cyclonic versus anticyclonic) and between sampling seasons/years and was significantly correlated with the taxonomic diversity of particle-associated heterotrophic bacteria. The expressed functional profile of these taxonomically distinct microeukaryotes varied consistently as a function of eddy polarity, with cyclones having a different expression pattern than anticyclones, and between sampling seasons/years. These data suggest that mesoscale forcing, and associated changes in biogeochemistry, could drive specific physiological responses in the resident microeukaryote community, independent of species composition.

The Problem with ‘Microbiome’

The term “microbiome” is currently applied predominantly to assemblages of organisms with 16S rRNA genes. In this context, “microbiome” is a misnomer that has been conferred a wide-ranging primacy over terms for community members lacking such genes, e.g., mycobiome, eukaryome, and virome, yet these are also important subsets of microbial communities. Widespread convenient and affordable 16S rRNA sequencing pipelines have accelerated continued use of such a “microbiome”, but at what intellectual and practical costs? Here we show that the use of “microbiome” in ribosomal gene-based studies has been egregiously misapplied, and discuss potential impacts. We argue that the current focus of “microbiome” research, predominantly on only ‘bacteria’, presents a dangerous narrowing of scope which encourages dismissal and even ignorance of other organisms’ contributions to microbial diversity, sensu stricto, and as etiologic agents; we put this in context by discussing cases in both marine microbial diversity and the role of pathogens in global amphibian decline. Fortunately, the solution is simple. We must use descriptive nouns that strictly reflect the outcomes attainable by the methods used. “Microbiome”, as a descriptive noun, should only be used when diversity in the three recognized domains is explored.

Shift and Metabolic Potentials of Microbial Eukaryotic Communities Across the Full Depths of the Mariana Trench

Microbial eukaryotes are widespread and play important roles in marine ecosystems. However, their ecological characteristics in the deep sea (>1,000 m), especially hadal trenches, were largely unknown. Here, we investigated the diversity and metabolic potentials of microbial eukaryotes along the whole water column of the Mariana Trench by metagenomics. Our results showed clear depth-related distribution of microbial eukaryotic community and associated metabolic potentials. Surface seawater was dominated by phototrophic/mixotrophic groups (e.g., Dinoflagellata) and genes involved in biosynthesis (photosynthesis and fatty acid biosynthesis), while deep (bathypelagic and/or hadal) seawaters were enriched with heterotrophic groups (e.g., Bicoecea) and genes related to digestion (lysosomal enzymes and V-type ATPase) and carbohydrate metabolism. Co-occurrence analysis revealed high intra-domain connectivity, indicating that microbial eukaryotic composition was more influenced by microbial eukaryotes themselves than bacteria. Increased abundance of genes associated with unsaturated fatty acid biosynthesis likely plays a role in resisting high hydrostatic pressure. Top1 and hupB genes, responsible for the formation and stabilization of DNA structure, were unique and abundant in the hadal zone and thus may be helpful to stabilize DNA structure in the deep sea. Overall, our results provide insights into the distribution and potential adaptability of microbial eukaryotes in the hadal zone.

Differential recruitment of opportunistic taxa leads to contrasting abilities in carbon processing by bathypelagic and surface microbial communities

Different factors affect the way dissolved organic matter (DOM) is processed in the ocean water column, including environmental conditions and the functional capabilities of the communities. Recent studies have shown that bathypelagic prokaryotes are metabolically flexible, but whether this versatility translates into a higher ability to process DOM has been barely explored. Here we performed a multifactorial transplant experiment to compare the growth, activity and changes in DOM quality in surface and bathypelagic waters inoculated with either surface or bathypelagic prokaryotic communities. The effect of nutrient additions to surface waters was also explored. Despite no differences in the cell abundance of surface and deep ocean prokaryotes were observed in any of the treatments, in surface waters with nutrients the heterotrophic production of surface prokaryotes rapidly decreased. Conversely, bathypelagic communities displayed a sustained production throughout the experiment. Incubations with surface prokaryotes always led to a significant accumulation of recalcitrant compounds, which did not occur with bathypelagic prokaryotes, suggesting they have a higher ability to process DOM. These contrasting abilities could be explained by the recruitment of a comparatively larger number of opportunistic taxa within the bathypelagic assemblages, which likely resulted in a broader community capability of substrate utilization.

Diversity of Pelagic and Benthic Bacterial Assemblages in the Western Pacific Ocean

Despite numerous studies on marine prokaryotes, the vertical distribution patterns of bacterial community, either on the taxonomic composition or the functional structure, remains relatively unexplored. Using HiSeq-derived 16S rRNA data, the depth-related distribution patterns of taxonomic diversity and functional structure predicted from diversity data in the water column and sediments of the Western Pacific Ocean were explored. The OTU richness declined along the water column after peaking between 100 to 200 m deep. Relative abundance of Cyanobacteria and SAR11 decreased significantly with depth, while Actinobacteria and Gammaproteobacteria increased. This clearly mirrors the vertical distribution pattern of the predicted functional composition with the shift between phototrophic to chemoheterotrophic groups from the surface to the deeper layers. In terms of community composition and functional structure, the epipelagic zone differed from other deeper ones (i.e., meso-, bathy-, and abyssopelagic zones) where no obvious differences were detected. For the epipelagic zone, temperature, dissolved oxygen, and salinity were recognized as the crucial factors shaping both community composition and the functional structure of bacteria. Compared with water samples, benthic sediment samples harbored unexpectedly higher read abundance of Proteobacteria, presenting distinguishable taxonomic and functional compositions. This study provides novel knowledge on the vertical distribution of bacterial taxonomic and functional compositions in the western Pacific.

Vertical shifts of particle-attached and free-living prokaryotes in the water column above the cold seeps of the South China Sea

Marine particle-attached (PA) and free-living (FL) microbes play important roles in the biogeochemical cycling of organic matter along the water column. Deep-sea cold seeps are highly productive and chemosynthetic ecosystems, their continuous emission of CH₄, CO₂, and H₂S can reach up to 100 m in the above water, therefore would influence the distribution and potential metabolic functions of deep-sea prokaryotes. In this study, the vertical distribution profiles of both PA and FL microbes in the water column above two cold seeps of the South China Sea were investigated using Illumina sequencing and quantitative PCR (qPCR) based on 16S rRNA gene. Photosynthetic and heterotrophic prokaryotes were predominant in respective surface and deep layers below the photic zone. The typical cold seep chemosynthetic microbes, such as methanotrophs and sulfate-reducing bacteria were observed with low proportions in the two cold seeps as well. Distinct PA and FL microbial fractions were found in terms of abundance and diversity. FL fraction exposed to the bulk water was significantly affected by temperature and inorganic nutrients, whereas PA fraction relied more on the organic matter of the particles and less susceptible to the environmental variability. Our study highlights the importance of vertical geochemical gradients on the distribution and potential metabolic choice of marine microbes and extends our current knowledge of depth-associated microbial distribution patterns.

A multidimensional perspective on microbial interactions

Beyond being simply positive or negative, beneficial or inhibitory, microbial interactions can involve a diverse set of mechanisms, dependencies and dynamical properties. These more nuanced features have been described in great detail for some specific types of interactions, (e.g. pairwise metabolic cross-feeding, quorum sensing or antibiotic killing), often with the use of quantitative measurements and insight derived from modeling. With a growing understanding of the composition and dynamics of complex microbial communities for human health and other applications, we face the challenge of integrating information about these different interactions into comprehensive quantitative frameworks. Here, we review the literature on a wide set of microbial interactions, and explore the potential value of a formal categorization based on multidimensional vectors of attributes. We propose that such an encoding can facilitate systematic, direct comparisons of interaction mechanisms and dependencies, and we discuss the relevance of an atlas of interactions for future modeling and rational design efforts.

Depth-Differentiation and Seasonality of Planktonic Microbial Assemblages in the Monterey Bay Upwelling System

Coastal upwelling regions are hotspots of biological productivity, supporting diverse communities of microbial life and metabolisms. Monterey Bay (MB), a coastal ocean embayment in central California, experiences seasonal upwelling of cold, nutrient-rich waters that sustain episodes of high phytoplankton production in surface waters. While productivity in surface waters is intimately linked to metabolisms of diverse communities of Archaea and Bacteria, a comprehensive understanding of the microbial community in MB is missing thus far, particularly in relation to the distinct hydrographic seasons characteristic of the MB system. Here we present the results of a 2-year microbial time-series survey in MB, investigating community composition and structure across spatiotemporal gradients. In deciphering these patterns, we used unique sequence variants (SVs) of the 16S rRNA gene (V4–V5 region), complemented with metagenomes and metatranscriptomes representing multiple depth profiles. We found clear depth-differentiation and recurring seasonal abundance patterns within planktonic communities, particularly when analyzed at finer taxonomic levels. Compositional changes were more pronounced in the upper 0–40 m of the water column, whereas deeper depths were characterized by temporally stable populations. In accordance with the dynamic nutrient profiles, the system appears to change from a Bacteroidetes- and Rhodobacterales-dominated upwelling period to an oceanic season dominated by oligotrophic groups such as SAR11 and picocyanobacteria. The cascade of environmental changes brought about by upwelling and relaxation events thus impacts microbial community structure in the bay, with important implications for the temporal variability of nutrient and energy fluxes within the MB ecosystem. Our observations emphasize the need for continued monitoring of planktonic microbial communities in order to predict and manage the behavior of this sensitive marine sanctuary ecosystem, over projected intensification of upwelling in the region.

Diatoms Are Selective Segregators in Global Ocean Planktonic Communities

Diatoms are key phytoplankton in the modern ocean that are involved in numerous biotic interactions, ranging from symbiosis to predation and viral infection, which have considerable effects on global biogeochemical cycles. However, despite recent large-scale studies of plankton, we are still lacking a comprehensive picture of the diversity of diatom biotic interactions in the marine microbial community. Through the ecological interpretation of both inferred microbial association networks and available knowledge on diatom interactions compiled in an open-access database, we propose an ecosystems approach for exploring diatom interactions in the ocean.

Diatoms are a major component of phytoplankton, believed to be responsible for around 20% of the annual primary production on Earth. As abundant and ubiquitous organisms, they are known to establish biotic interactions with many other members of plankton. Through analyses of cooccurrence networks derived from the Tara Oceans expedition that take into account both biotic and abiotic factors in shaping the spatial distributions of species, we show that only 13% of diatom pairwise associations are driven by environmental conditions; the vast majority are independent of abiotic factors. In contrast to most other plankton groups, on a global scale, diatoms display a much higher proportion of negative correlations with other organisms, particularly toward potential predators and parasites, suggesting that their biogeography is constrained by top-down pressure. Genus-level analyses indicate that abundant diatoms are not necessarily the most connected and that species-specific abundance distribution patterns lead to negative associations with other organisms. In order to move forward in the biological interpretation of cooccurrence networks, an open-access extensive literature survey of diatom biotic interactions was compiled, of which 18.5% were recovered in the computed network. This result reveals the extent of what likely remains to be discovered in the field of planktonic biotic interactions, even for one of the best-known organismal groups.

IMPORTANCE Diatoms are key phytoplankton in the modern ocean that are involved in numerous biotic interactions, ranging from symbiosis to predation and viral infection, which have considerable effects on global biogeochemical cycles. However, despite recent large-scale studies of plankton, we are still lacking a comprehensive picture of the diversity of diatom biotic interactions in the marine microbial community. Through the ecological interpretation of both inferred microbial association networks and available knowledge on diatom interactions compiled in an open-access database, we propose an ecosystems approach for exploring diatom interactions in the ocean.

Marked changes in diversity and relative activity of picoeukaryotes with depth in the world ocean

Microbial eukaryotes are key components of the ocean plankton. Yet, our understanding of their community composition and activity in different water layers of the ocean is limited, particularly for picoeukaryotes (0.2–3 µm cell size). Here, we examined the picoeukaryotic communities inhabiting different vertical zones of the tropical and subtropical global ocean: surface, deep chlorophyll maximum, mesopelagic (including the deep scattering layer and oxygen minimum zones), and bathypelagic. Communities were analysed by high-tthroughput sequencing of the 18S rRNA gene (V4 region) as represented by DNA (community structure) and RNA (metabolism), followed by delineation of Operational Taxonomic Units (OTUs) at 99% similarity. We found a stratification of the picoeukaryotic communities along the water column, with assemblages corresponding to the sunlit and dark ocean. Specific taxonomic groups either increased (e.g., Chrysophyceae or Bicosoecida) or decreased (e.g., Dinoflagellata or MAST-3) in abundance with depth. We used the rRNA:rDNA ratio of each OTU as a proxy of metabolic activity. The highest relative activity was found in the mesopelagic layer for most taxonomic groups, and the lowest in the bathypelagic. Altogether, we characterize the change in community structure and metabolic activity of picoeukaryotes with depth in the global ocean, suggesting a hotspot of activity in the mesopelagic.

Biodiversity patterns across taxonomic groups along a lake water-depth gradient: Effects of abiotic and biotic drivers

Understanding biodiversity patterns and the role of biotic attributes in governing these patterns remains one of the most important challenges in ecology. Here, taking water depth in Lake Lugu as a typical geographical gradient, we studied how these different taxa, that is bacteria, diatoms and chironomids, respond to the water depth and environmental gradients using molecular and morphological methods. We further evaluated the relative importance of water depth, environmental variables and biotic attributes in explaining biological characteristics, such as biomass, species richness, and community composition. The biomass of chironomids and the richness of bacteria and chironomids showed a nonlinearly decreasing pattern associated with increased water depth, while biomass and species richness of diatoms showed U-shaped and hump-shaped patterns, respectively. The three taxonomic groups all showed increasing dissimilarity with water depth changes, and there was clear cross-taxon congruence among the variations in community composition. Abiotic variables were pivotal in structuring biological characteristics; however, the biotic attributes also explained a unique portion of their variations. This suggests that biotic interactions significantly influenced the patterns of biomass, species richness, and community compositions along the water depth gradient for the three taxonomic groups studied. Our results provide new evidence that biotic attributes could help in predicting the biodiversity of aquatic communities along geographical gradients, such as water depth.

Distinct temporal dynamics of planktonic archaeal and bacterial assemblages in the bays of the Yellow Sea

The Yellow Sea features unique characteristics due to strong tides and nutrient-enriched freshwater outflows from China and Korea. The coupling of archaeal and bacterial assemblages associated with environmental factors at two bay areas in the Yellow Sea was investigated. Temporal variations of the archaeal and bacterial assemblages were shown to be greater than the spatial variations based on an analysis of the 16S rRNA gene sequences. Distinct temporal dynamics of both planktonic archaeal and bacterial assemblages was associated with temperature, NO₂^-, and chlorophyll a ([chl-a]) concentrations in the bays of the Yellow Sea. The [chl-a] was the prime predictor of bacterial abundance, and some taxa were clearly correlated with [chl-a]. Bacteroidetes and Alpha-proteobacteria dominated at high [chl-a] stations while Gamma-proteobacteria (esp. SAR86 clade) and Actinobacteria (Candidatus Actinomarina clade) were abundant at low [chl-a] stations. The archaeal abundance was comparable with the bacterial abundance in most of the October samples. Co-dominance of Marine Group II (MGII) and Candidatus Nitrosopumilus suggests that the assimilation of organic nitrogen by MGII could be coupled with nitrification by ammonia-oxidizing archaea. The distinct temporal dynamics of the archaeal and bacterial assemblages might be attributable to the strong tides and the inflow of nutrient-rich freshwater.

Influence of 16S rRNA Hypervariable Region on Estimates of Bacterial Diversity and Community Composition in Seawater and Marine Sediment

To assess the influence of 16S ribosomal RNA (rRNA) tag choice on estimates of microbial diversity and/or community composition in seawater and marine sediment, we examined bacterial diversity and community composition from a site in the Central North Atlantic and a site in the Equatorial Pacific. For each site, we analyzed samples from four zones in the water column, a seafloor sediment sample, and two subseafloor sediment horizons (with stratigraphic ages of 1.5 and 5.5 million years old). We amplified both the V4 and V6 hypervariable regions of the 16S rRNA gene and clustered the sequences into operational taxonomic units (OTUs) of 97% similarity to analyze for diversity and community composition. OTU richness is much higher with the V6 tag than with the V4 tag, and subsequently OTU-level community composition is quite different between the two tags. Vertical patterns of relative diversity are broadly the same for both tags, with maximum taxonomic richness in seafloor sediment and lowest richness in subseafloor sediment at both geographic locations. Genetic dissimilarity between sample locations is also broadly the same for both tags. Community composition is very similar for both tags at the class level, but very different at the level of 97% similar OTUs. Class-level diversity and community composition of water-column samples are very similar at each water depth between the Atlantic and Pacific. However, sediment communities differ greatly from the Atlantic site to the Pacific site. Finally, for relative patterns of diversity and class-level community composition, deep sequencing and shallow sequencing provide similar results.

The Relative Abundance of Benthic Bacterial Phyla Along a Water-Depth Gradient in a Plateau Lake: Physical, Chemical, and Biotic Drivers

Water-depth biodiversity gradient, one of the typical biogeographical patterns on Earth, is understudied for bacteria in freshwater ecosystems, and thus left the underlying mechanisms poorly understood especially for benthic bacteria. Here, we investigated the water-depth distribution of surface sediment bacterial phyla and their driving factors in Lake Lugu, a plateau lake in Southwest China. Our results revealed that the relative abundance of 11 dominant bacterial phyla showed various water-depth patterns, such as increasing, decreasing, hump-shaped, and U-shaped patterns. These patterns across phyla were consistent with their different niche positions of water depth, while the occupancy-abundance relationships were not dependent on phylum attributes. Consistently, phylum abundance was best explained by water depth; other physical and chemical factors, such as metal ion concentrations, SiO₂, and pH, can also explain the variations in some bacterial phyla. Chemical variables were the main drivers of the dominant bacterial phyla. However, biotic variables also showed substantial importance for some phyla, such as Planctomycetes, Actinobacteria, and WS3. This work could provide new insights into the general water-depth patterns and underlying mechanisms of the relative abundance of bacterial phyla in freshwater ecosystems.

Differential co-occurrence relationships shaping ecotype diversification within Thaumarchaeota populations in the coastal ocean water column

Ecological factors contributing to depth-related diversification of marine Thaumarchaeota populations remain largely unresolved. To investigate the role of potential microbial associations in shaping thaumarchaeal ecotype diversification, we examined co-occurrence relationships in a community composition dataset (16S rRNA V4-V5 region) collected as part of a 2-year time series in coastal Monterey Bay. Ecotype groups previously defined based on functional gene diversity-water column A (WCA), water column B (WCB) and Nitrosopumilus-like clusters-were recovered in the thaumarchaeal 16S rRNA gene phylogeny. Networks systematically reflected depth-related patterns in the abundances of ecotype populations, suggesting thaumarchaeal ecotypes as keystone members of the microbial community below the euphotic zone. Differential environmental controls on the ecotype populations were further evident in subnetwork modules showing preferential co-occurrence of OTUs belonging to the same ecotype cluster. Correlated abundances of Thaumarchaeota and heterotrophic bacteria (e.g., Bacteroidetes, Marinimicrobia and Gammaproteobacteria) indicated potential reciprocal interactions via dissolved organic matter transformations. Notably, the networks recovered ecotype-specific associations between thaumarchaeal and Nitrospina OTUs. Even at depths where WCB-like Thaumarchaeota dominated, Nitrospina OTUs were found to preferentially co-occur with WCA-like and Nitrosopumilus-like thaumarchaeal OTUs, highlighting the need to investigate the ecological implications of the composition of nitrifier assemblages in marine waters.

Heterogeneity of Microbial Communities on Deep-Sea Ferromanganese Crusts in the Takuyo-Daigo Seamount

Rock outcrops of aged deep-sea seamounts are generally covered with Fe and Mn oxides, known as ferromanganese (Fe-Mn) crusts. Although the presence of microorganisms in Fe-Mn crusts has been reported, limited information is currently available on intra- and inter-variations in crust microbial communities. Therefore, we collected several Fe-Mn crusts in bathyal and abyssal zones (water depths of 1,150-5,520 m) in the Takuyo-Daigo Seamount in the northwestern Pacific, and examined microbial communities on the crusts using culture-independent molecular and microscopic analyses. Quantitative PCR showed that microbial cells were abundant (10⁶-10⁸ cells g^-1) on Fe-Mn crust surfaces through the water depths. A comparative 16S rRNA gene analysis revealed community differences among Fe-Mn crusts through the water depths, which may have been caused by changes in dissolved oxygen concentrations. Moreover, community differences were observed among positions within each Fe-Mn crust, and potentially depended on the availability of sinking particulate organic matter. Microscopic and elemental analyses of thin Fe-Mn crust sections revealed the accumulation of microbial cells accompanied by the depletion of Mn in valleys of bumpy crust surfaces. Our results suggest that heterogeneous and abundant microbial communities play a role in the biogeochemical cycling of Mn, in addition to C and N, on crusts and contribute to the extremely slow growth of Fe-Mn crusts.

Metabolic Modeling of Microbial Community Interactions for Health, Environmental and Biotechnological Applications

In nature, microbes do not exist in isolation but co-exist in a variety of ecological and biological environments and on various host organisms. Due to their close proximity, these microbes interact among themselves, and also with the hosts in both positive and negative manners. Moreover, these interactions may modulate dynamically upon external stimulus as well as internal community changes. This demands systematic techniques such as mathematical modeling to understand the intrinsic community behavior. Here, we reviewed various approaches for metabolic modeling of microbial communities. If detailed species-specific information is available, segregated models of individual organisms can be constructed and connected via metabolite exchanges; otherwise, the community may be represented as a lumped network of metabolic reactions. The constructed models can then be simulated to help fill knowledge gaps, and generate testable hypotheses for designing new experiments. More importantly, such community models have been developed to study microbial interactions in various niches such as host microbiome, biogeochemical and bioremediation, waste water treatment and synthetic consortia. As such, the metabolic modeling efforts have allowed us to gain new insights into the natural and synthetic microbial communities, and design interventions to achieve specific goals. Finally, potential directions for future development in metabolic modeling of microbial communities were also discussed.

Winter picoplankton diversity in an oligotrophic marginal sea

Marine picoplankton, unicellular organisms with cell sizes up to 3 μm in diameter, numerically dominate marine ecosystems, encompassing Archaea, Bacteria, Eukarya (protists and fungi) as well as viruses. Autotrophic and heterotrophic picoplankton abundance and community composition with a focus on picoeukaryotes (PEs) were investigated in the winter of 2016 at three stations along a coast-to-offshore transect in the southern Adriatic Sea. Abundances were estimated by flow cytometry, while community composition by Illumina High-Throughput Sequencing (HTS) of 16S and 18S rRNA genes. The photosynthetic picoplankton diversity was also investigated by High-Performance Liquid Chromatography (HPLC) of liposoluble pigments. Heterotrophic bacteria and cyanobacteria (Prochlorococcus and Synechococcus) accounted for up to 7 × 10⁵; 2.3 × 10⁴ and 2.5 × 10⁴ cells mL^-1, respectively, while photosynthetic picoeukaryotes peaked with 3 × 10³ cells mL^-1. Prokaryotes, as revealed by HTS were dominated by Alphaproteobacteria (mainly SAR11, 44.91% of total 16S sequence reads), followed by Gammaproteobacteria (Oceanospirillales and Pseudomonadales, 14.96%), Bacteroidetes (mainly Flavobacteriales, 13%), Cyanobacteria (Prochlorococcus and Synechococcus, 9.52%), Marinimicrobia (SAR406, 7.97%), Deltaproteobacteria (SAR324, 3.83%), Actinobacteria (2.24%) and Chloroflexi (SAR202, 1.90%). Photosynthetic pigment concentrations were very low (12.12 μgL^-1 at the most) and taxonomic pigments could be attributed to Prochlorococcus, Synechococcus, Prymnesiophyceae, Bacillariophyceae, Chrysophyceae, and Prasinophyceae. HTS data revealed that PEs were dominated by heterotrophs, such as Syndiniophyceae, parasitic dinoflagellates (79.67% of total 18S sequence reads), Dinophyceae (8.7%) and the radiolarians Collodaria belonging to Sphaerozoidae (22.1%) and Spumellaria (5.0%). On the other hand, photoautotrophs, including Chlorophyta (Mamiellophyceae, Prasinophyceae, Trebouxiophyceae, and Ulvophyceae), Stramenopiles (Bacillariophyta, Chrysophyceae, Dictyochophyceae, Pelagophyceae), photoautotrophic Cryptophyta and some Haptophyta (Prymnesiophyceae), did not exceed 5% of total sequence reads. This study provides the first snapshot of the PEs diversity in oligotrophic euphotic waters of the southern Adriatic Sea, hence setting the stage for large-scale surveying and characterization of the eukaryotic diversity in the entire basin.

Vertical Distribution of Microbial Eukaryotes From Surface to the Hadal Zone of the Mariana Trench

Marine microbial eukaryotes are ubiquitous, comprised of phylogenetically diverse groups and play key roles in microbial food webs and global biogeochemical cycling. However, their vertical distribution in the deep sea has received little attention. In this study, we investigated the composition and diversity of the eukaryotes of both 0.2–3 μm and >3 μm size fractions from the surface to the hadal zone (8727 m) of the Mariana Trench using Illumina MiSeq sequencing for the 18S rDNA. The microbial eukaryotic community structure differed substantially across size fractions and depths. Operational taxonomic unit (OTU) richness in the >3 μm fraction was higher than that in the 0.2–3 μm fraction at the same depth. For the 0.2–3 μm fraction, sequences of Retaria (Rhizaria) were most abundant in the surface water (53.5%). Chrysophyceae (Stramenopiles) sequences dominated mostly in the samples from water depths below 1795 m. For the >3 μm fraction, sequences of Dinophyceae (Alveolata) were most abundant in surface waters (49.3%) and remained a significant proportion of total sequences at greater depths (9.8%, on average). Retaria sequences were abundant in samples of depths ≥1000 m. Amoebozoa and Apusozoa sequences were enriched in the hadal sample, comprising 38 and 20.4% of total sequences, respectively. Fungi (Opisthokonta) sequences were most abundant at 1759 m in both size fractions. Strong positive associations were found between Syndiniales (mainly MALV-I and MALV-II) and Retaria while negative associations were shown between MALV-II and Fungi in a co-occurrence analysis. This study compared the community structure of microbial eukaryotes in different zones in the deep sea and identified a distinct hadal community in the larger size fraction, suggesting the uniqueness of the eukaryotes in the biosphere in the Mariana Trench.

Oceanographic boundaries constrain microbial diversity gradients in the South Pacific Ocean

Marine microbes along with microeukaryotes are key regulators of oceanic biogeochemical pathways. Here we present a high-resolution (every 0.5° of latitude) dataset describing microbial pro- and eukaryotic richness in the surface and just below the thermocline along a 7,000-km transect from 66°S at the Antarctic ice edge to the equator in the South Pacific Ocean. The transect, conducted in austral winter, covered key oceanographic features including crossing of the polar front (PF), the subtropical front (STF), and the equatorial upwelling region. Our data indicate that temperature does not determine patterns of marine microbial richness, complementing the global model data from Ladau et al. [Ladau J, et al. (2013) ISME J 7:1669-1677]. Rather, NH₄⁺, nanophytoplankton, and primary productivity were the main drivers for archaeal and bacterial richness. Eukaryote richness was highest in the least-productive ocean region, the tropical oligotrophic province. We also observed a unique diversity pattern in the South Pacific Ocean: a regional increase in archaeal and bacterial diversity between 10°S and the equator. Rapoport's rule describes the tendency for the latitudinal ranges of species to increase with latitude. Our data showed that the mean latitudinal ranges of archaea and bacteria decreased with latitude. We show that permanent oceanographic features, such as the STF and the equatorial upwelling, can have a significant influence on both alpha-diversity and beta-diversity of pro- and eukaryotes.

Agulhas Current properties shape microbial community diversity and potential functionality

Understanding the impact of oceanographic features on marine microbial ecosystems remains a major ecological endeavour. Here we assess microbial diversity, community structure and functional capacity along the Agulhas Current system and the Subtropical Front in the South Indian Ocean (SIO). Samples collected from the epipelagic, oxygen minimum and bathypelagic zones were analysed by 16S rRNA gene amplicon and metagenomic sequencing. In contrast to previous studies, we found high taxonomic richness in surface and deep water samples, but generally low richness for OMZ communities. Beta-diversity analysis revealed significant dissimilarity between the three water depths. Most microbial communities were dominated by marine Gammaproteobacteria, with strikingly low levels of picocyanobacteria. Community composition was strongly influenced by specific environmental factors including depth, salinity, and the availability of both oxygen and light. Carbon, nitrogen and sulfur cycling capacity in the SIO was linked to several autotrophic and copiotrophic Alphaproteobacteria and Gammaproteobacteria. Taken together, our data suggest that the environmental conditions in the Agulhas Current system, particularly depth-related parameters, substantially influence microbial community structure. In addition, the capacity for biogeochemical cycling of nitrogen and sulfur is linked primarily to the dominant Gammaproteobacteria taxa, whereas ecologically rare taxa drive carbon cycling.

Sinking particles promote vertical connectivity in the ocean microbiome

Prokaryotes dominate the living biomass and the biological diversity of the ocean, one of the largest ecosystems on earth. The sinking of particles is a widespread mechanism that transports materials to the deep ocean, with a significant role in the global carbon cycle. Whether this process constitutes a global dispersal pathway for prokaryotic diversity connecting surface communities to those in the dark ocean has never been tested. Here we show that surface and deep-sea prokaryotic communities are strongly connected, constituting a vast oceanic metacommunity where local assemblages are linked through the transport of sinking particles. This vertical dispersal, mediated mainly by the largest sinking particles, emerges as a fundamental process shaping the assembly and biogeography of deep ocean prokaryotic communities.

The sinking of organic particles formed in the photic layer is a main vector of carbon export into the deep ocean. Although sinking particles are heavily colonized by microbes, so far it has not been explored whether this process plays a role in transferring prokaryotic diversity from surface to deep oceanic layers. Using Illumina sequencing of the 16S rRNA gene, we explore here the vertical connectivity of the ocean microbiome by characterizing marine prokaryotic communities associated with five different size fractions and examining their compositional variability from surface down to 4,000 m across eight stations sampled in the Atlantic, Pacific, and Indian Oceans during the Malaspina 2010 Expedition. Our results show that the most abundant prokaryotes in the deep ocean are also present in surface waters. This vertical community connectivity seems to occur predominantly through the largest particles because communities in the largest size fractions showed the highest taxonomic similarity throughout the water column, whereas free-living communities were more isolated vertically. Our results further suggest that particle colonization processes occurring in surface waters determine to some extent the composition and biogeography of bathypelagic communities. Overall, we postulate that sinking particles function as vectors that inoculate viable particle-attached surface microbes into the deep-sea realm, determining to a considerable extent the structure, functioning, and biogeography of deep ocean communities.

Bacterial Diversity and Nitrogen Utilization Strategies in the Upper Layer of the Northwestern Pacific Ocean

Nitrogen (N) is a primary limiting nutrient for bacterial growth and productivity in the ocean. To better understand bacterial community and their N utilization strategy in different N regimes of the ocean, we examined bacterial diversity, diazotrophic diversity, and N utilization gene expressions in the northwestern Pacific Ocean (NWPO) using a combination of high-throughput sequencing and real-time qPCR methods. 521 and 204 different operational taxonomic units (OTUs) were identified in the 16s rRNA and nifH libraries from nine surface samples. Of the 16s rRNA gene OTUs, 11.9% were observed in all samples while 3.5 and 15.9% were detected only in N-sufficient and N-deficient samples. Proteobacteria, Cyanobacteria and Bacteroidetes dominated the bacterial community. Prochlorococcus and Pseudoalteromonas were the most abundant at the genus level in N-deficient regimes, while SAR86, Synechococcus and SAR92 were predominant in the Kuroshio-Oyashio confluence region. The distribution of the nifH gene presented great divergence among sampling stations: Cyanobacterium_UCYN-A dominated the N-deficient stations, while clusters related to the Alpha-, Beta-, and Gamma-Proteobacteria were abundant in other stations. Temperature was the main factor that determined bacterial community structure and diversity while concentration of NO_X-N was significantly correlated with structure and distribution of N₂-fixing microorganisms. Expression of the ammonium transporter was much higher than that of urea transporter subunit A (urtA) and ferredoxin-nitrate reductase, while urtA had an increased expression in N-deficient surface water. The predicted ammonium transporter and ammonium assimilation enzymes were most abundant in surface samples while urease and nitrogenase were more abundant in the N-deficient regions. These findings underscore the fact that marine bacteria have evolved diverse N utilization strategies to adapt to different N habitats, and that urea metabolism is of vital ecological importance in N-deficient regimes.

Vertically distinct microbial communities in the Mariana and Kermadec trenches

Hadal trenches, oceanic locations deeper than 6,000 m, are thought to have distinct microbial communities compared to those at shallower depths due to high hydrostatic pressures, topographical funneling of organic matter, and biogeographical isolation. Here we evaluate the hypothesis that hadal trenches contain unique microbial biodiversity through analyses of the communities present in the bottom waters of the Kermadec and Mariana trenches. Estimates of microbial protein production indicate active populations under in situ hydrostatic pressures and increasing adaptation to pressure with depth. Depth, trench of collection, and size fraction are important drivers of microbial community structure. Many putative hadal bathytypes, such as members related to the Marinimicrobia, Rhodobacteraceae, Rhodospirilliceae, and Aquibacter, are similar to members identified in other trenches. Most of the differences between the two trench microbiomes consists of taxa belonging to the Gammaproteobacteria whose distributions extend throughout the water column. Growth and survival estimates of representative isolates of these taxa under deep-sea conditions suggest that some members may descend from shallower depths and exist as a potentially inactive fraction of the hadal zone. We conclude that the distinct pelagic communities residing in these two trenches, and perhaps by extension other trenches, reflect both cosmopolitan hadal bathytypes and ubiquitous genera found throughout the water column.

Depth-Resolved Distribution of Particle-Attached and Free-Living Bacterial Communities in the Water Column of the New Britain Trench

Particle-attached (PA) and free-living (FL) microorganisms play significant but different roles in mineralization of organic matter (OM) in the ocean. Currently, little is known about PA and FL microbial communities in bathyal and abyssal pelagic waters, and understanding of their diversity and distribution in the water column and their interactions with environmental factors in the trench area is limited. We investigated for the first time the variations of abundance and diversities of the PA and FL bacterial communities in the epi-, bathy-, and abyssopelagic zones of the New Britain Trench (NBT). The PA communities showed decreasing species richness but increasing relative abundance with depth, suggesting the increasing ecological significance of the PA bacteria in the deep ocean. The abundance and diversity of PA and FL bacterial communities in the NBT water column appeared to be shaped by different sets of environment factors, which might be related to different micro-niches of the two communities. Analysis on species distribution suggested that the differences between PA and FL bacteria communities mainly resulted from the different relative abundance of the “shared taxa” in the two types of communities. These findings provide valuable information for understanding the relative ecological roles of the PA and FL bacterial communities and their interactions with environmental factors in different pelagic zones along the vertical profile of the NBT water column.

Extracellular DNA as a genetic recorder of microbial diversity in benthic deep-sea ecosystems

Extracellular DNA in deep-sea sediments represents a major repository of genes, which previously belonged to living organisms. However, the extent to which these extracellular genes influence current estimates of prokaryotic biodiversity is unknown. We investigated the abundance and diversity of 16S rDNA sequences contained within extracellular DNA from continental margins of different biogeographic regions. We also compared the taxonomic composition of microbial assemblages through the analysis of extracellular DNA and DNA associated with living cells. 16S rDNA contained in the extracellular DNA pool contributed up to 50% of the total 16S rDNA copy number determined in the sediments. Ca. 4% of extracellular Operational Taxonomic Units (OTUs) were shared among the different biogeographic regions revealing the presence of a core of preserved OTUs. A higher fraction of OTUs was exclusive of each region potentially due to its geographic and thermohaline characteristics. Ca. one third of the OTUs identified in the extracellular DNA were absent from living prokaryotic assemblages, possibly representing the signatures of past assemblages. Our findings expand the knowledge of the contribution of extracellular microbial sequences to current estimates of prokaryotic diversity obtained through the analyses of “environmental DNA”, and open new perspectives for understanding microbial successions in benthic ecosystems.

Metagenomics of Bacterial Diversity in Villa Luz Caves with Sulfur Water Springs

New biotechnology applications require in-depth preliminary studies of biodiversity. The methods of massive sequencing using metagenomics and bioinformatics tools offer us sufficient and reliable knowledge to understand environmental diversity, to know new microorganisms, and to take advantage of their functional genes. Villa Luz caves, in the southern Mexican state of Tabasco, are fed by at least 26 groundwater inlets, containing 300-500 mg L-1 H2S and <0.1 mg L-1 O2. We extracted environmental DNA for metagenomic analysis of collected samples in five selected Villa Luz caves sites, with pH values from 2.5 to 7. Foreign organisms found in this underground ecosystem can oxidize H2S to H2SO4. These include: biovermiculites, a bacterial association that can grow on the rock walls; snottites, that are whitish, viscous biofilms hanging from the rock walls, and sacks or bags of phlegm, which live within the aquatic environment of the springs. Through the emergency food assistance program (TEFAP) pyrosequencing, a total of 20,901 readings of amplification products from hypervariable regions V1 and V3 of 16S rRNA bacterial gene in whole and pure metagenomic DNA samples were generated. Seven bacterial phyla were identified. As a result, Proteobacteria was more frequent than Acidobacteria. Finally, acidophilic Proteobacteria was detected in UJAT5 sample.

Paddy soil microbial communities driven by environment- and microbe-microbe interactions: A case study of elevation-resolved microbial communities in a rice terrace

Rice paddies are a significant source of the greenhouse gas methane, which mainly originates from microbial activity. Methane generation in anaerobic systems involves complex interactions of multiple functional microbial groups. Rice paddies installed in hilly terrain are often terraced, providing multiple quasi-independent plots differing primarily in their elevation up a hillside. This represents an excellent study site to explore the influence of environmental factors on microbial communities and interactions among microbial populations. In this study, we used a combination of geochemical analyses, high-throughput amplicon sequencing, and statistical methods to elucidate these interactions. Sulfate, total nitrogen, total iron, and total organic carbon were determined to be critical factors in steering the ecosystem composition and function. Sulfate-reducing bacteria predominated in the rice terrace microbial communities, and Fe(III)-reducing and methane-oxidizing bacteria were abundant as well. Biotic interactions indicated by co-occurrence network analysis suggest mutualistic interactions among these three functional groups. Paddy-scale methane production may be affected by competition among methanogens and sulfate- and Fe(III)-reducing bacteria, or by direct methane oxidation by methane-oxidizing bacteria.

CAPSULE

Microbial communities were characterized in rice terrace. The environment- and microbe-microbe interactions indicated the mitigation of sulfate and Fe on methane production.

Diversity of Microbial Communities and Quantitative Chemodiversity in Layers of Marine Sediment Cores from a Causeway (Kaichu-Doro) in Okinawa Island, Japan

Microbial community diversity and chemodiversity were investigated in marine sediments adjacent to the Okinawan “Kaichu-Doro” Causeway, which was constructed 46 years ago to connect a group of four islands (Henza-jima, Miyagi-jima, Ikei-jima, Hamahiga-jima) to the Okinawan main island. This causeway was not built on pilings, but by land reclamation; hence, it now acts as a long, thin peninsula. The construction of this causeway was previously shown to have influenced the surrounding marine ecosystem, causing ecosystem fragmentation and loss of water circulation. In this study, we collected sediment cores (n = 10) from five paired sites in 1 m water depths. Each pair of sites consisted of one site each on the immediate north and south sides of the causeway. Originally the members of each pair were much closer to each other (<150 m) than to other pairs, but now the members of each pair are isolated by the causeway. Each core was 60–80 cm long and was divided into 15-cm layers. We examined the vertical diversity of microbial communities and chemical compounds to determine the correlation between chemodiversity and microbial communities among marine sediment cores and layers. Principal coordinate analyses (PCoA) of detected compounds and of bacterial and archaeal operational taxonomic units (OTUs) revealed that the north and south sides of the causeway are relatively isolated, with each side having unique microbial OTUs. Additionally, some bacterial families (e.g., Acidaminobacteraceae, Rhizobiaceae, and Xanthomonadaceae) were found only on the south side of Kaichu-Doro. Interestingly, we found that the relative abundance of OTUs for some microbial families increased from top to bottom, but this was reversed in some other families. We conclude that the causeway has altered microbial community composition and metabolite profiles in marine sediments.

Censusing marine eukaryotic diversity in the twenty-first century

The ocean constitutes one of the vastest and richest biomes on our planet. Most recent estimations, all based on indirect approaches, suggest that there are millions of marine eukaryotic species. Moreover, a large majority of these are small (less than 1 mm), cryptic and still unknown to science. However, this knowledge gap, caused by the lack of diagnostic morphological features in small organisms and the limited sampling of the global ocean, is currently being filled, thanks to new DNA-based approaches. The molecular technique of PCR amplification of homologous gene regions combined with high-throughput sequencing, routinely used to census unculturable prokaryotes, is now also being used to characterize whole communities of marine eukaryotes. Here, we review how this methodological advancement has helped to better quantify the magnitude and patterns of marine eukaryotic diversity, with an emphasis on taxonomic groups previously largely overlooked. We then discuss obstacles remaining to achieve a global understanding of marine eukaryotic diversity. In particular, we argue that 18S variable regions do not provide sufficient taxonomic resolution to census marine life, and suggest combining broad eukaryotic surveys targeting the 18S rRNA region with more taxon-focused analyses of hypervariable regions to improve our understanding of the diversity of species, the functional units of marine ecosystems.

This article is part of the themed issue ‘From DNA barcodes to biomes’.

Spatial variability of marine bacterial and archaeal communities along the particulate matter continuum

Biotic and abiotic particles shape the microspatial architecture that defines the microbial aquatic habitat, being particles highly variable in size and quality along oceanic horizontal and vertical gradients. We analysed the prokaryotic (bacterial and archaeal) diversity and community composition present in six distinct particle size classes ranging from the pico- to the microscale (0.2 to 200 μm). Further, we studied their variations along oceanographic horizontal (from the coast to open oceanic waters) and vertical (from the ocean surface into the meso- and bathypelagic ocean) gradients. In general, prokaryotic community composition was more variable with depth than in the transition from the coast to the open ocean. Comparing the six size-fractions, distinct prokaryotic communities were detected in each size-fraction, and whereas bacteria were more diverse in the larger size-fractions, archaea were more diverse in the smaller size-fractions. Comparison of prokaryotic community composition among particle size-fractions showed that most, but not all, taxonomic groups have a preference for a certain size-fraction sustained with depth. Species sorting, or the presence of diverse ecotypes with distinct size-fraction preferences, may explain why this trend is not conserved in all taxa.

Basin Scale Variation on the Composition and Diversity of Archaea in the Pacific Ocean

The Archaea are a widely distributed group of prokaryotes that inhabit and thrive in many different environments. In the sea, they play key roles in various global biogeochemical processes. Here, in order to investigate the vertical profiles of archaeal community across a large geographic distance, the compositions of archaeal communities in seven seawater columns in the Pacific Ocean were investigated using high throughput 454 pyrosequencing of the 16S rRNA gene. The surface archaeal communities showed lower diversity and greater variability than those in the deeper layers. Two of the major archaeal phyla that displayed different depth preferences were Thaumarchaeota and Euryarchaeota. The majority of Thaumarchaeota belonged to Marine Group I (MGI), which had high relative abundance in deep water. In contrast, Euryarchaeota, which mainly consisted of Marine Group II (MGII) and III (MGIII), were dominant in the surface layer. Compared with MGI and MGII, MGIII were less abundant in seawater and generally absent from the surface water of the subarctic Pacific. In addition, niche separation in the MGI, MGII, and MGIII subgroups was also observed. For example, MGI.C and MGII.A (the major subgroups of MGI and MGII, respectively) displayed a strong negative correlation with each other. The highest level of archaeal diversity was found in the core of an oxygen minimum zone (OMZ) located off Costa Rica, which resulted from the co-occurrence of both anaerobic and aerobic archaea. For example, methanotrophic archaea ANME-2, methanogenic archaea and several sediment origin archaea, such as Marine Benthic Group A (MBGA) and Bathyarchaeota, were all detected at relatively high abundance in the OMZ. Together, our findings indicate that vertical heterogeneities along water columns and latitudinal differentiation in the surface waters are ubiquitous features of archaeal communities in the Pacific Ocean, and the OMZ off Costa Rica is an archaeal biodiversity hot-spot.

From the Surface to the Deep-Sea: Bacterial Distributions across Polymetallic Nodule Fields in the Clarion-Clipperton Zone of the Pacific Ocean

Marine bacteria regulate fluxes of matter and energy essential for pelagic and benthic organisms and may also be involved in the formation and maintenance of commercially valuable abyssal polymetallic nodules. Future mining of these nodule fields is predicted to have substantial effects on biodiversity and physicochemical conditions in mined areas. Yet, the identity and distributions of bacterial populations in deep-sea sediments and associated polymetallic nodules has received relatively little attention. We examined bacterial communities using high-throughput sequencing of bacterial 16S rRNA gene fragments from samples collected in the water column, sediment, and polymetallic nodules in the Pacific Ocean (bottom depth ≥4,000 m) in the eastern Clarion-Clipperton Zone. Operational taxonomic units (OTUs; defined at 99% 16S rRNA gene identity) affiliated with JTB255 (Gammaproteobacteria) and Rhodospirillaceae (Alphaproteobacteria) had higher relative abundances in the nodule and sediment habitats compared to the water column. Rhodobiaceae family and Vibrio OTUs had higher relative abundance in nodule samples, but were less abundant in sediment and water column samples. Bacterial communities in sediments and associated with nodules were generally similar; however, 5,861 and 6,827 OTUs found in the water column were retrieved from sediment and nodule habitats, respectively. Cyanobacterial OTUs clustering among Prochlorococcus and Synechococcus were detected in both sediments and nodules, with greater representation among nodule samples. Such results suggest that vertical export of typically abundant photic-zone microbes may be an important process in delivery of water column microorganisms to abyssal habitats, potentially influencing the structure and function of communities in polymetallic nodule fields.