Charting the Complexity of the Marine Microbiome through Single-Cell Genomics.

Marine bacteria and archaea play key roles in global biogeochemistry. To improve our understanding of this complex microbiome, we employed single-cell genomics and a randomized, hypothesis-agnostic cell selection strategy to recover 12,715 partial genomes from the tropical and subtropical euphotic ocean. A substantial fraction of known prokaryoplankton coding potential was recovered from a single, 0.4 mL ocean sample, which indicates that genomic information disperses effectively across the globe. Yet, we found each genome to be unique, implying limited clonality within prokaryoplankton populations. Light harvesting and secondary metabolite biosynthetic pathways were numerous across lineages, highlighting the value of single-cell genomics to advance the identification of ecological roles and biotechnology potential of uncultured microbial groups. This genome collection enabled functional annotation and genus-level taxonomic assignments for >80% of individual metagenome reads from the tropical and subtropical surface ocean, thus offering a model to improve reference genome databases for complex microbiomes.

Microbes and microplastics: Community shifts along an urban coastal contaminant gradient.

Plastic substrates introduced to the environment during the Anthropocene have introduced new pathways for microbial selection and dispersal. Some plastic-colonising microorganisms have adapted phenotypes for plastic degradation (selection), while the spatial transport (dispersal) potential of plastic colonisers remains controlled by polymer-specific density, hydrography and currents. Plastic-degrading enzyme abundances have recently been correlated with concentrations of plastic debris in open ocean environments, making it critical to better understand colonisation of hydrocarbon degraders with plastic degradation potential in urbanised watersheds where plastic pollution often originates. We found that microbial colonisation by reputed hydrocarbon degraders on microplastics (MPs) correlated with a spatial contaminant gradient (New York City/Long Island waterways), polymer types, temporal scales, microbial domains and putative cell activity (DNA vs. RNA). Hydrocarbon-degrading taxa enriched on polyethylene and polyvinyl chloride substrates relative to other polymers and were more commonly recovered in samples proximal to New York City. These differences in MP colonisation could indicate phenotypic adaptation processes resulting from increased exposure to urban plastic runoff as well as differences in carbon bioavailability across polymer types. Shifts in MP community potential across urban coastal contaminant gradients and polymer types improve our understanding of environmental plastic discharge impacts toward biogeochemical cycling across the global ocean.

Protistan grazing impacts microbial communities and carbon cycling at deep-sea hydrothermal vents.

Microbial eukaryotes (or protists) in marine ecosystems are a link between primary producers and all higher trophic levels, and the rate at which heterotrophic protistan grazers consume microbial prey is a key mechanism for carbon transport and recycling in microbial food webs. At deep-sea hydrothermal vents, chemosynthetic bacteria and archaea form the base of a food web that functions in the absence of sunlight, but the role of protistan grazers in these highly productive ecosystems is largely unexplored. Here, we pair grazing experiments with a molecular survey to quantify protistan grazing and to characterize the composition of vent-associated protists in low-temperature diffuse venting fluids from Gorda Ridge in the northeast Pacific Ocean. Results reveal protists exert higher predation pressure at vents compared to the surrounding deep seawater environment and may account for consuming 28 to 62% of the daily stock of prokaryotic biomass within discharging hydrothermal vent fluids. The vent-associated protistan community was more species rich relative to the background deep sea, and patterns in the distribution and co-occurrence of vent microbes provide additional insights into potential predator-prey interactions. Ciliates, followed by dinoflagellates, Syndiniales, rhizaria, and stramenopiles, dominated the vent protistan community and included bacterivorous species, species known to host symbionts, and parasites. Our findings provide an estimate of protistan grazing pressure within hydrothermal vent food webs, highlighting the important role that diverse protistan communities play in deep-sea carbon cycling.

Parasitic infections by Group II Syndiniales target selected dinoflagellate host populations within diverse protist assemblages in a model coastal pond.

Protists are integral to marine food webs and biogeochemical cycles; however, there is a paucity of data describing specific ecological niches for some of the most abundant taxa in marker gene libraries. Syndiniales are one such group, often representing the majority of sequence reads recovered from picoplankton samples across the global ocean. However, the prevalence and impacts of syndinian parasitism in marine environments remain unclear. We began to address these critical knowledge gaps by generating a high-resolution time series (March-October 2018) in a productive coastal pond. Seasonal shifts in protist populations, including parasitic Syndiniales, were documented during periods of higher primary productivity and increased summer temperature-driven stratification. Elevated concentrations of infected hosts and free-living parasite spores occurred at nearly monthly intervals in July, August, and September. We suggest intensifying stratification during this period correlated with the increased prevalence of dinoflagellates that were parasitized by Group II Syndiniales. Infections in some protist populations were comparable to previously reported large single-taxon dinoflagellate blooms. Infection dynamics in Salt Pond demonstrated the propagation of syndinian parasites through mixed protist assemblages and highlighted patterns of host/parasite interactions that better reflect many other marine environments where single taxon blooms are uncommon.

Diverse nitrogen cycling pathways across a marine oxygen gradient indicate nitrogen loss coupled to chemoautotrophic activity.

Genetic markers and geochemical assays of microbial nitrogen cycling processes, including autotrophic and heterotrophic denitrification, anammox, ammonia oxidation, and nitrite oxidation, were examined across the oxycline, suboxic, and anoxic zones of the Cariaco Basin, Venezuela. Ammonia and nitrite oxidation genes were expressed through the entire gradient. Transcripts associated with autotrophic and heterotrophic denitrifiers were mostly confined to the suboxic zone and below but were also present in particles in the oxycline. Anammox genes and transcripts were detected over a narrow depth range near the bottom of the suboxic zone and coincided with secondary NO2 - maxima and available NH4 + . Dissolved inorganic nitrogen (DIN) amendment incubations and comparisons between our sampling campaigns suggested that denitrifier activity may be closely coupled with NO3 - availability. Expression of denitrification genes at depths of high rates of chemoautotrophic carbon fixation and phylogenetic analyses of nitrogen cycling genes and transcripts indicated a diverse array of denitrifiers, including chemoautotrophs capable of using NO3 - to oxidize reduced sulfur species. Thus, results suggest that the Cariaco Basin nitrogen cycle is influenced by autotrophic carbon cycling in addition to organic matter oxidation and anammox.

Organic matter processing by microbial communities throughout the Atlantic water column as revealed by metaproteomics.

The phylogenetic composition of the heterotrophic microbial community is depth stratified in the oceanic water column down to abyssopelagic layers. In the layers below the euphotic zone, it has been suggested that heterotrophic microbes rely largely on solubilized particulate organic matter as a carbon and energy source rather than on dissolved organic matter. To decipher whether changes in the phylogenetic composition with depth are reflected in changes in the bacterial and archaeal transporter proteins, we generated an extensive metaproteomic and metagenomic dataset of microbial communities collected from 100- to 5,000-m depth in the Atlantic Ocean. By identifying which compounds of the organic matter pool are absorbed, transported, and incorporated into microbial cells, intriguing insights into organic matter transformation in the deep ocean emerged. On average, solute transporters accounted for 23% of identified protein sequences in the lower euphotic and ∼39% in the bathypelagic layer, indicating the central role of heterotrophy in the dark ocean. In the bathypelagic layer, substrate affinities of expressed transporters suggest that, in addition to amino acids, peptides and carbohydrates, carboxylic acids and compatible solutes may be essential substrates for the microbial community. Key players with highest expression of solute transporters were Alphaproteobacteria, Gammaproteobacteria, and Deltaproteobacteria, accounting for 40%, 11%, and 10%, respectively, of relative protein abundances. The in situ expression of solute transporters indicates that the heterotrophic prokaryotic community is geared toward the utilization of similar organic compounds throughout the water column, with yet higher abundances of transporters targeting aromatic compounds in the bathypelagic realm.

Meta-omics highlights the diversity, activity and adaptations of fungi in deep oceanic crust.

The lithified oceanic crust, lower crust gabbros in particular, has remained largely unexplored by microbiologists. Recently, evidence for heterogeneously distributed viable and transcriptionally active autotrophic and heterotrophic microbial populations within low-biomass communities was found down to 750 m below the seafloor at the Atlantis Bank Gabbro Massif, Indian Ocean. Here, we report on the diversity, activity and adaptations of fungal communities in the deep oceanic crust from ~10 to 780 mbsf by combining metabarcoding analyses with mid/high-throughput culturing approaches. Metabarcoding along with culturing indicate a low diversity of viable fungi, mostly affiliated to ubiquitous (terrestrial and aquatic environments) taxa. Ecophysiological analyses coupled with metatranscriptomics point to viable and transcriptionally active fungal populations engaged in cell division, translation, protein modifications and other vital cellular processes. Transcript data suggest possible adaptations for surviving in the nutrient-poor, lithified deep biosphere that include the recycling of organic matter. These active communities appear strongly influenced by the presence of cracks and veins in the rocks where fluids and resulting rock alteration create micro-niches.

Free-living chemoautotrophic and particle-attached heterotrophic prokaryotes dominate microbial assemblages along a pelagic redox gradient.

Using the anoxic Cariaco Basin as a natural laboratory, particle association of bacterial and archaeal taxa was assessed by iTag sequencing and qPCR gene assays of samples spanning an oxic-anoxic-euxinic gradient. A total of 10%-12% of all bacterial and archaeal cells were found in the particle-associated (PA) fraction, operationally defined as prokaryotes captured on 2.7 µm membranes. Both redox condition and size fraction segregated bacterial taxa. Archaeal taxa varied according to redox conditions, but were similar between size fractions. Taxa putatively associated with chemoautotrophic sulfur oxidation and nitrification dominated the free-living (FL) fraction throughout the oxycline (< 1-120 µM O2 ) and upper anoxic layer. Bacteria in the oxycline's PA fraction included taxa known to be aerobic and anaerobic chemoorganotrophs. At shallow anoxic depths, PA taxa were primarily affiliated with anaerobic sulfate ( SO42-)-reducing lineages. PA fractions in the most sulfidic samples were dominated by taxa affiliated with CH4 oxidizing, fermenting and SO42- reducing lineages. Prevalence of particle-associated SO42--reducing taxa and abundant sulfur-oxidizing taxa in both size fractions across the oxic-anoxic interface is consistent with the cryptic sulfur cycling concept. Bacterial assemblage diversity in the PA fraction always exceeded the FL fraction except in the most oxic samples, whereas Archaeal diversity was not consistently different between size fractions. Our results suggest that these particle-associated and free-living bacterial assemblages are functionally different and that the interplay between particle microhabitats and surrounding geochemical regimes is a strong selective force shaping microbial communities throughout the water column.

Metazoans of redoxcline sediments in Mediterranean deep-sea hypersaline anoxic basins.

BACKGROUND: The deep-sea hypersaline anoxic basins (DHABs) of the Mediterranean (water depth ~3500 m) are some of the most extreme oceanic habitats known. Brines of DHABs are nearly saturated with salt, leading many to suspect they are uninhabitable for eukaryotes. While diverse bacterial and protistan communities are reported from some DHAB haloclines and brines, loriciferans are the only metazoan reported to inhabit the anoxic DHAB brines. Our goal was to further investigate metazoan communities in DHAB haloclines and brines. RESULTS: We report observations from sediments of three DHAB (Urania, Discovery, L'Atalante) haloclines, comparing these to observations from sediments underlying normoxic waters of typical Mediterranean salinity. Due to technical difficulties, sampling of the brines was not possible. Morphotype analysis indicates nematodes are the most abundant taxon; crustaceans, loriciferans and bryozoans were also noted. Among nematodes, Daptonema was the most abundant genus; three morphotypes were noted with a degree of endemicity. The majority of rRNA sequences were from planktonic taxa, suggesting that at least some individual metazoans were preserved and inactive. Nematode abundance data, in some cases determined from direct counts of sediments incubated in situ with CellTracker(TM) Green, was patchy but generally indicates the highest abundances in either normoxic control samples or in upper halocline samples; nematodes were absent or very rare in lower halocline samples. Ultrastructural analysis indicates the nematodes in L'Atalante normoxic control sediments were fit, while specimens from L'Atalante upper halocline were healthy or had only recently died and those from the lower halocline had no identifiable organelles. Loriciferans, which were only rarely encountered, were found in both normoxic control samples as well as in Discovery and L'Atalante haloclines. It is not clear how a metazoan taxon could remain viable under this wide range of conditions. CONCLUSIONS: We document a community of living nematodes in normoxic, normal saline deep-sea Mediterranean sediments and in the upper halocline portions of the DHABs. Occurrences of nematodes in mid-halocline and lower halocline samples did not provide compelling evidence of a living community in those zones. The possibility of a viable metazoan community in brines of DHABs is not supported by our data at this time.

Depth shapes α- and ß-diversities of microbial eukaryotes in surficial sediments of coastal ecosystems.

Little is known about the relative influence of historic processes and environmental gradients on shaping the diversity of single-celled eukaryotes in marine benthos. By combining pyrosequencing of 18S ribosomal RNA genes with data on multiple environmental factors, we investigated the diversity of microeukaryotes in surficial sediments of three basins of the Yellow Sea Large Marine Ecosystem. A considerable proportion (about 20%) of reads was affiliated with known parasitoid protists. Dinophyta and Ciliophora appeared dominant in terms of relative proportion of reads and operational taxonomic unit (OTU) richness. Overall, OTU richness of benthic microeukaryotes decreased with increasing water depth and decreasing pH. While community composition was significantly different among basins, partial Mantel tests indicated a depth-decay pattern of community similarity, whereby water depth, rather than geographic distance or environment, shaped ß-diversity of benthic microeukaryotes (including both the abundant and the rare biosphere) on a regional scale. Similar hydrographic and mineralogical factors contributed to the biogeography of both the abundant and the rare OTUs. The trace metal vanadium had a significant effect on the biogeography of the rare biosphere. Our study sheds new light on the composition, diversity patterns and underlying mechanisms of single-celled eukaryote distribution in surficial sediments of coastal oceans.

Inter-comparison of the potentially active prokaryotic communities in the halocline sediments of Mediterranean deep-sea hypersaline basins.

The sediment microbiota of the Mediterranean deep-sea anoxic hypersaline basins (DHABs) are understudied relative to communities in the brines and halocline waters. In this study, the active fraction of the prokaryotic community in the halocline sediments of L' Atalante, Urania, and Discovery DHABs was investigated based on extracted total RNA and 454 pyrosequencing of the 16S rRNA gene. Bacterial and archaeal communities were different in the sediments underlying the halocline waters of the three habitats, reflecting the unique chemical settings of each basin. The relative abundance of unique operational taxonomic units (OTUs) was also different between deep-sea control sediments and sediments underlying DHAB haloclines, suggesting adaptation to the steep DHAB chemical gradients. Only a few OTUs were affiliated to known bacterial halophilic and/or anaerobic groups. Many OTUs, including some of the dominant ones, were related to aerobic taxa. Archaea were detected only in few halocline samples, with lower OTU richness relative to Bacteria, and were dominated by taxa associated with methane cycling. This study suggests that, while metabolically active prokaryotic communities appear to be present in sediments underlying the three DHABs investigated, their diversity and activity are likely to be more reduced in sediments underlying the brines.

Ciliates along oxyclines of permanently stratified marine water columns.

Studies of microbial communities in areas of the world where permanent marine water column oxyclines exist suggest they are "hotspots" of microbial activity, and that these water features and the anoxic waters below them are inhabited by diverse protist taxa, including ciliates. These communities have minimal taxonomic overlap with those in overlying oxic water columns. Some ciliate taxa have been detected in multiple locations where these stable water column oxyclines exist; however, differences in such factors as hydrochemistry in the habitats that have been studied suggest local selection for distinct communities. We compare published data on ciliate communities from studies of deep marine water column oxyclines in Caricao Basin, Venezuela, and the Black Sea, with data from coastal, shallower oxycline waters in Framvaren and Mariager fjords, and from several deep-sea hypersaline anoxic basins in the Eastern Mediterranean Sea. Putative symbioses between Bacteria, Archaea, and ciliates observed along these oxyclines suggests a strategy of cooperative metabolism for survival that includes chemosynthetic autotrophy and exchanges of metabolic intermediates or end products between hosts and their prokaryotic partners.

Plasmid-Borne Biosynthetic Gene Clusters within a Permanently Stratified Marine Water Column.

Plasmids are mobile genetic elements known to carry secondary metabolic genes that affect the fitness and survival of microbes in the environment. Well-studied cases of plasmid-encoded secondary metabolic genes in marine habitats include toxin/antitoxin and antibiotic biosynthesis/resistance genes. Here, we examine metagenome-assembled genomes (MAGs) from the permanently-stratified water column of the Cariaco Basin for integrated plasmids that encode biosynthetic gene clusters of secondary metabolites (smBGCs). We identify 16 plasmid-borne smBGCs in MAGs associated primarily with Planctomycetota and Pseudomonadota that encode terpene-synthesizing genes, and genes for production of ribosomal and non-ribosomal peptides. These identified genes encode for secondary metabolites that are mainly antimicrobial agents, and hence, their uptake via plasmids may increase the competitive advantage of those host taxa that acquire them. The ecological and evolutionary significance of smBGCs carried by prokaryotes in oxygen-depleted water columns is yet to be fully elucidated.

Predicting metabolic modules in incomplete bacterial genomes with MetaPathPredict.

The reconstruction of complete microbial metabolic pathways using 'omics data from environmental samples remains challenging. Computational pipelines for pathway reconstruction that utilize machine learning methods to predict the presence or absence of KEGG modules in incomplete genomes are lacking. Here, we present MetaPathPredict, a software tool that incorporates machine learning models to predict the presence of complete KEGG modules within bacterial genomic datasets. Using gene annotation data and information from the KEGG module database, MetaPathPredict employs deep learning models to predict the presence of KEGG modules in a genome. MetaPathPredict can be used as a command line tool or as a Python module, and both options are designed to be run locally or on a compute cluster. Benchmarks show that MetaPathPredict makes robust predictions of KEGG module presence within highly incomplete genomes.

Microbial eukaryotic predation pressure and biomass at deep-sea hydrothermal vents.

Deep-sea hydrothermal vent geochemistry shapes the foundation of the microbial food web by fueling chemolithoautotrophic microbial activity. Microbial eukaryotes (or protists) play a critical role in hydrothermal vent food webs as consumers and hosts of symbiotic bacteria, and as a nutritional source to higher trophic levels. We measured microbial eukaryotic cell abundance and predation pressure in low-temperature diffuse hydrothermal fluids at the Von Damm and Piccard vent fields along the Mid-Cayman Rise in the Western Caribbean Sea. We present findings from experiments performed under in situ pressure that show cell abundances and grazing rates higher than those done at 1 atmosphere (shipboard ambient pressure); this trend was attributed to the impact of depressurization on cell integrity. A relationship between the protistan grazing rate, prey cell abundance, and temperature of end-member hydrothermal vent fluid was observed at both vent fields, regardless of experimental approach. Our results show substantial protistan biomass at hydrothermally fueled microbial food webs, and when coupled with improved grazing estimates, suggest an important contribution of grazers to the local carbon export and supply of nutrient resources to the deep ocean.

Diverse secondary metabolites are expressed in particle-associated and free-living microorganisms of the permanently anoxic Cariaco Basin.

Secondary metabolites play essential roles in ecological interactions and nutrient acquisition, and are of interest for their potential uses in medicine and biotechnology. Genome mining for biosynthetic gene clusters (BGCs) can be used for the discovery of new compounds. Here, we use metagenomics and metatranscriptomics to analyze BGCs in free-living and particle-associated microbial communities through the stratified water column of the Cariaco Basin, Venezuela. We recovered 565 bacterial and archaeal metagenome-assembled genomes (MAGs) and identified 1154 diverse BGCs. We show that differences in water redox potential and microbial lifestyle (particle-associated vs. free-living) are associated with variations in the predicted composition and production of secondary metabolites. Our results indicate that microbes, including understudied clades such as Planctomycetota, potentially produce a wide range of secondary metabolites in these anoxic/euxinic waters.

Eukaryotic genomes from a global metagenomic data set illuminate trophic modes and biogeography of ocean plankton.

Metagenomics is a powerful method for interpreting the ecological roles and physiological capabilities of mixed microbial communities. Yet, many tools for processing metagenomic data are neither designed to consider eukaryotes nor are they built for an increasing amount of sequence data. EukHeist is an automated pipeline to retrieve eukaryotic and prokaryotic metagenome-assembled genomes (MAGs) from large-scale metagenomic sequence data sets. We developed the EukHeist workflow to specifically process large amounts of both metagenomic and/or metatranscriptomic sequence data in an automated and reproducible fashion. Here, we applied EukHeist to the large-size fraction data (0.8-2,000 µm) from Tara Oceans to recover both eukaryotic and prokaryotic MAGs, which we refer to as TOPAZ (Tara Oceans Particle-Associated MAGs). The TOPAZ MAGs consisted of >900 environmentally relevant eukaryotic MAGs and >4,000 bacterial and archaeal MAGs. The bacterial and archaeal TOPAZ MAGs expand upon the phylogenetic diversity of likely particle- and host-associated taxa. We use these MAGs to demonstrate an approach to infer the putative trophic mode of the recovered eukaryotic MAGs. We also identify ecological cohorts of co-occurring MAGs, which are driven by specific environmental factors and putative host-microbe associations. These data together add to a number of growing resources of environmentally relevant eukaryotic genomic information. Complementary and expanded databases of MAGs, such as those provided through scalable pipelines like EukHeist, stand to advance our understanding of eukaryotic diversity through increased coverage of genomic representatives across the tree of life.IMPORTANCESingle-celled eukaryotes play ecologically significant roles in the marine environment, yet fundamental questions about their biodiversity, ecological function, and interactions remain. Environmental sequencing enables researchers to document naturally occurring protistan communities, without culturing bias, yet metagenomic and metatranscriptomic sequencing approaches cannot separate individual species from communities. To more completely capture the genomic content of mixed protistan populations, we can create bins of sequences that represent the same organism (metagenome-assembled genomes [MAGs]). We developed the EukHeist pipeline, which automates the binning of population-level eukaryotic and prokaryotic genomes from metagenomic reads. We show exciting insight into what protistan communities are present and their trophic roles in the ocean. Scalable computational tools, like EukHeist, may accelerate the identification of meaningful genetic signatures from large data sets and complement researchers' efforts to leverage MAG databases for addressing ecological questions, resolving evolutionary relationships, and discovering potentially novel biodiversity.

Microbial eukaryote life in the new hypersaline deep-sea basin Thetis.

Only recently, a novel anoxic hypersaline (thalassic) basin in the eastern Mediterranean was discovered at a depth of 3,258 m. The halite-saturated brine of this polyextreme basin revealed one of the highest salt concentrations ever reported for such an environment (salinity of 348‰). Using a eukaryote-specific probe and fluorescence in situ hybridization, we counted 0.6 × 10(4) protists per liter of anoxic brine. SSU rRNA sequence analyses, based on amplification of environmental cDNA identified fungi as the most diverse taxonomic group of eukaryotes in the brine, making deep-sea brines sources of unknown fungal diversity and hotspots for the discovery of novel metabolic pathways and for secondary metabolites. The second most diverse phylotypes are ciliates and stramenopiles (each 20%). The occurrence of closely related ciliate sequences exclusively in other Mediterranean brine basins suggests specific adaptations of the respective organisms to such habitats. Betadiversity-analyses confirm that microeukaryote communities in the brine and the interface are notably different. Several distinct morphotypes in brine samples suggest that the rRNA sequences detected in Thetis brine can be linked to indigenous polyextremophile protists. This contradicts previous assumptions that such extremely high salt concentrations are anathema to eukaryotic life. The upper salinity limits for eukaryotic life remain unidentified.

Diversity and spatial distribution of prokaryotic communities along a sediment vertical profile of a deep-sea mud volcano.

We investigated the top 30-cm sediment prokaryotic community structure in 5-cm spatial resolution, at an active site of the Amsterdam mud volcano, East Mediterranean Sea, based on the 16S rRNA gene diversity. A total of 339 and 526 sequences were retrieved, corresponding to 25 and 213 unique (≥98% similarity) phylotypes of Archaea and Bacteria, respectively, in all depths. The Shannon-Wiener diversity index H was higher for Bacteria (1.92-4.03) than for Archaea (0.99-1.91) and varied differently between the two groups. Archaea were dominated by anaerobic methanotrophs ANME-1, -2 and -3 groups and were related to phylotypes involved in anaerobic oxidation of methane from similar habitats. The much more complex Bacteria community consisted of 20 phylogenetic groups at the phylum/candidate division level. Proteobacteria, in particular δ-Proteobacteria, was the dominant group. In most sediment layers, the dominant phylotypes of both the Archaea and Bacteria communities were found in neighbouring layers, suggesting some overlap in species richness. The similarity of certain prokaryotic communities was also depicted by using four different similarity indices. The direct comparison of the retrieved phylotypes with those from the Kazan mud volcano of the same field revealed that 40.0% of the Archaea and 16.9% of the Bacteria phylotypes are common between the two systems. The majority of these phylotypes are closely related to phylotypes originating from other mud volcanoes, implying a degree of endemicity in these systems.

Eukaryotic Parasites Are Integral to a Productive Microbial Food Web in Oxygen-Depleted Waters.

Oxygen-depleted water columns (ODWCs) host a diverse community of eukaryotic protists that change dramatically in composition over the oxic-anoxic gradient. In the permanently anoxic Cariaco Basin, peaks in eukaryotic diversity occurred in layers where dark microbial activity (chemoautotrophy and heterotrophy) were highest, suggesting a link between prokaryotic activity and trophic associations with protists. Using 18S rRNA gene sequencing, parasites and especially the obligate parasitic clade, Syndiniales, appear to be particularly abundant, suggesting parasitism is an important, but overlooked interaction in ODWC food webs. Syndiniales were also associated with certain prokaryotic groups that are often found in ODWCs, including Marinimicrobia and Marine Group II archaea, evocative of feedbacks between parasitic infection events, release of organic matter, and prokaryotic assimilative activity. In a network analysis that included all three domains of life, bacterial and archaeal taxa were putative bottleneck and hub species, while a large proportion of edges were connected to eukaryotic nodes. Inclusion of parasites resulted in a more complex network with longer path lengths between members. Together, these results suggest that protists, and especially protistan parasites, play an important role in maintaining microbial food web complexity, particularly in ODWCs, where protist diversity and microbial productivity are high, but energy resources are limited relative to euphotic waters.