Direct observations of microbial community succession on sinking marine particles.

Microbial community dynamics on sinking particles control the amount of carbon that reaches the deep ocean and the length of time that carbon is stored, with potentially profound impacts on Earth's climate. A mechanistic understanding of the controls on sinking particle distributions has been hindered by limited depth- and time-resolved sampling and methods that cannot distinguish individual particles. Here, we analyze microbial communities on nearly 400 individual sinking particles in conjunction with more conventional composite particle samples to determine how particle colonization and community assembly might control carbon sequestration in the deep ocean. We observed community succession with corresponding changes in microbial metabolic potential on the larger sinking particles transporting a significant fraction of carbon to the deep sea. Microbial community richness decreased as particles aged and sank; however, richness increased with particle size and the attenuation of carbon export. This suggests that the theory of island biogeography applies to sinking marine particles. Changes in POC flux attenuation with time and microbial community composition with depth were reproduced in a mechanistic ecosystem model that reflected a range of POC labilities and microbial growth rates. Our results highlight microbial community dynamics and processes on individual sinking particles, the isolation of which is necessary to improve mechanistic models of ocean carbon uptake.

Synechococcus nitrogen gene loss in iron-limited ocean regions.

Synechococcus are the most abundant cyanobacteria in high latitude regions and are responsible for an estimated 17% of annual marine net primary productivity. Despite their biogeochemical importance, Synechococcus populations have been unevenly sampled across the ocean, with most studies focused on low-latitude strains. In particular, the near absence of Synechococcus genomes from high-latitude, High Nutrient Low Chlorophyll (HNLC) regions leaves a gap in our knowledge of picocyanobacterial adaptations to iron limitation and their influence on carbon, nitrogen, and iron cycles. We examined Synechococcus populations from the subarctic North Pacific, a well-characterized HNLC region, with quantitative metagenomics. Assembly with short and long reads produced two near complete Synechococcus metagenome-assembled genomes (MAGs). Quantitative metagenome-derived abundances of these populations matched well with flow cytometry counts, and the Synechococcus MAGs were estimated to comprise >99% of the Synechococcus at Station P. Whereas the Station P Synechococcus MAGs contained multiple genes for adaptation to iron limitation, both genomes lacked genes for uptake and assimilation of nitrate and nitrite, suggesting a dependence on ammonium, urea, and other forms of recycled nitrogen leading to reduced iron requirements. A global analysis of Synechococcus nitrate reductase abundance in the TARA Oceans dataset found nitrate assimilation genes are also lower in other HNLC regions. We propose that nitrate and nitrite assimilation gene loss in Synechococcus may represent an adaptation to severe iron limitation in high-latitude regions where ammonium availability is higher. Our findings have implications for models that quantify the contribution of cyanobacteria to primary production and subsequent carbon export.

Carbon content, carbon fixation yield and dissolved organic carbon release from diverse marine nitrifiers.

Nitrifying microorganisms, including ammonia-oxidizing archaea, ammonia-oxidizing bacteria, and nitrite-oxidizing bacteria, are the most abundant chemoautotrophs in the ocean and play an important role in the global carbon cycle by fixing dissolved inorganic carbon (DIC) into biomass. The release of organic compounds by these microbes is not well quantified, but may represent an as-yet unaccounted source of dissolved organic carbon (DOC) available to marine food webs. Here, we provide measurements of cellular carbon and nitrogen quotas, DIC fixation yields and DOC release of 10 phylogenetically diverse marine nitrifiers. All investigated strains released DOC during growth, representing on average 5-15% of the fixed DIC. Changes in substrate concentration and temperature did not affect the proportion of fixed DIC released as DOC, but release rates varied between closely related species. Our results also indicate previous studies may have underestimated DIC fixation yields of marine nitrite oxidizers due to partial decoupling of nitrite oxidation from CO2 fixation, and due to lower observed yields in artificial compared to natural seawater medium. The results of this study provide critical values for biogeochemical models of the global carbon cycle, and help to further constrain the implications of nitrification-fueled chemoautotrophy for marine food-web functioning and the biological sequestration of carbon in the ocean.

Systematic comparison of unilamellar vesicles reveals that archaeal core lipid membranes are more permeable than bacterial membranes.

One of the deepest branches in the tree of life separates the Archaea from the Bacteria. These prokaryotic groups have distinct cellular systems including fundamentally different phospholipid membrane bilayers. This dichotomy has been termed the lipid divide and possibly bestows different biophysical and biochemical characteristics on each cell type. Classic experiments suggest that bacterial membranes (formed from lipids extracted from Escherichia coli, for example) show permeability to key metabolites comparable to archaeal membranes (formed from lipids extracted from Halobacterium salinarum), yet systematic analyses based on direct measurements of membrane permeability are absent. Here, we develop a new approach for assessing the membrane permeability of approximately 10 µm unilamellar vesicles, consisting of an aqueous medium enclosed by a single lipid bilayer. Comparing the permeability of 18 metabolites demonstrates that diether glycerol-1-phosphate lipids with methyl branches, often the most abundant membrane lipids of sampled archaea, are permeable to a wide range of compounds useful for core metabolic networks, including amino acids, sugars, and nucleobases. Permeability is significantly lower in diester glycerol-3-phosphate lipids without methyl branches, the common building block of bacterial membranes. To identify the membrane characteristics that determine permeability, we use this experimental platform to test a variety of lipid forms bearing a diversity of intermediate characteristics. We found that increased membrane permeability is dependent on both the methyl branches on the lipid tails and the ether bond between the tails and the head group, both of which are present on the archaeal phospholipids. These permeability differences must have had profound effects on the cell physiology and proteome evolution of early prokaryotic forms. To explore this further, we compare the abundance and distribution of transmembrane transporter-encoding protein families present on genomes sampled from across the prokaryotic tree of life. These data demonstrate that archaea tend to have a reduced repertoire of transporter gene families, consistent with increased membrane permeation. These results demonstrate that the lipid divide demarcates a clear difference in permeability function with implications for understanding some of the earliest transitions in cell origins and evolution.

Constraining the composition and quantity of organic matter used by abundant marine Thaumarchaeota.

Marine Group I (MGI) Thaumarchaeota were originally described as chemoautotrophic nitrifiers, but molecular and isotopic evidence suggests heterotrophic and/or mixotrophic capabilities. Here, we investigated the quantity and composition of organic matter assimilated by individual, uncultured MGI cells from the Pacific Ocean to constrain their potential for mixotrophy and heterotrophy. We observed that most MGI cells did not assimilate carbon from any organic substrate provided (glucose, pyruvate, oxaloacetate, protein, urea, and amino acids). The minority of MGI cells that did assimilate it did so exclusively from nitrogenous substrates (urea, 15% of MGI and amino acids, 36% of MGI), and only as an auxiliary carbon source (<20% of that subset's total cellular carbon was derived from those substrates). At the population level, MGI assimilation of organic carbon comprised just 0.5%-11% of total biomass carbon. We observed extensive assimilation of inorganic carbon and urea- and amino acid-derived nitrogen (equal to that from ammonium), consistent with metagenomic and metatranscriptomic analyses performed here and previously showing a widespread potential for MGI to perform autotrophy and transport and degrade organic nitrogen. Our results constrain the quantity and composition of organic matter used by MGI and suggest they use it primarily to meet nitrogen demands for anabolism and nitrification.

Microbial functional diversity across biogeochemical provinces in the central Pacific Ocean.

Enzymes catalyze key reactions within Earth's life-sustaining biogeochemical cycles. Here, we use metaproteomics to examine the enzymatic capabilities of the microbial community (0.2 to 3 µm) along a 5,000-km-long, 1-km-deep transect in the central Pacific Ocean. Eighty-five percent of total protein abundance was of bacterial origin, with Archaea contributing 1.6%. Over 2,000 functional KEGG Ontology (KO) groups were identified, yet only 25 KO groups contributed over half of the protein abundance, simultaneously indicating abundant key functions and a long tail of diverse functions. Vertical attenuation of individual proteins displayed stratification of nutrient transport, carbon utilization, and environmental stress. The microbial community also varied along horizontal scales, shaped by environmental features specific to the oligotrophic North Pacific Subtropical Gyre, the oxygen-depleted Eastern Tropical North Pacific, and nutrient-rich equatorial upwelling. Some of the most abundant proteins were associated with nitrification and C1 metabolisms, with observed interactions between these pathways. The oxidoreductases nitrite oxidoreductase (NxrAB), nitrite reductase (NirK), ammonia monooxygenase (AmoABC), manganese oxidase (MnxG), formate dehydrogenase (FdoGH and FDH), and carbon monoxide dehydrogenase (CoxLM) displayed distributions indicative of biogeochemical status such as oxidative or nutritional stress, with the potential to be more sensitive than chemical sensors. Enzymes that mediate transformations of atmospheric gases like CO, CO2, NO, methanethiol, and methylamines were most abundant in the upwelling region. We identified hot spots of biochemical transformation in the central Pacific Ocean, highlighted previously understudied metabolic pathways in the environment, and provided rich empirical data for biogeochemical models critical for forecasting ecosystem response to climate change.

The microbiome of a bacterivorous marine choanoflagellate contains a resource-demanding obligate bacterial associate.

Microbial predators such as choanoflagellates are key players in ocean food webs. Choanoflagellates, which are the closest unicellular relatives of animals, consume bacteria and also exhibit marked biological transitions triggered by bacterial compounds, yet their native microbiomes remain uncharacterized. Here we report the discovery of a ubiquitous, uncultured bacterial lineage we name Candidatus Comchoanobacterales ord. nov., related to the human pathogen Coxiella and physically associated with the uncultured marine choanoflagellate Bicosta minor. We analyse complete 'Comchoano' genomes acquired after sorting single Bicosta cells, finding signatures of obligate host-dependence, including reduction of pathways encoding glycolysis, membrane components, amino acids and B-vitamins. Comchoano encode the necessary apparatus to import energy and other compounds from the host, proteins for host-cell associations and a type IV secretion system closest to Coxiella's that is expressed in Pacific Ocean metatranscriptomes. Interactions between choanoflagellates and their microbiota could reshape the direction of energy and resource flow attributed to microbial predators, adding complexity and nuance to marine food webs.

Exaggerated trans-membrane charge of ammonium transporters in nutrient-poor marine environments.

Transporter proteins are a vital interface between cells and their environment. In nutrient-limited environments, microbes with transporters that are effective at bringing substrates into their cells will gain a competitive advantage over variants with reduced transport function. Microbial ammonium transporters (Amt) bring ammonium into the cytoplasm from the surrounding periplasm space, but diagnosing Amt adaptations to low nutrient environments solely from sequence data has been elusive. Here, we report altered Amt sequence amino acid distribution from deep marine samples compared to variants sampled from shallow water in two important microbial lineages of the marine water column community-Marine Group I Archaea (Thermoproteota) and the uncultivated gammaproteobacterial lineage SAR86. This pattern indicates an evolutionary pressure towards an increasing dipole in Amt for these clades in deep ocean environments and is predicted to generate stronger electric fields facilitating ammonium acquisition. This pattern of increasing dipole charge with depth was not observed in lineages capable of accessing alternative nitrogen sources, including the abundant alphaproteobacterial clade SAR11. We speculate that competition for ammonium in the deep ocean drives transporter sequence evolution. The low concentration of ammonium in the deep ocean is therefore likely due to rapid uptake by Amts concurrent with decreasing nutrient flux.

Complete Genome Sequences of Two Phylogenetically Distinct Nitrospina Strains Isolated from the Atlantic and Pacific Oceans.

The complete genome sequences of two chemoautotrophic nitrite-oxidizing bacteria of the genus Nitrospina are reported. Nitrospina gracilis strain Nb-211 was isolated from the Atlantic Ocean, and Nitrospina sp. strain Nb-3 was isolated from the Pacific Ocean. We report two highly similar ~3.07-Mbp genome sequences that differ by the presence of ferric iron chelator (siderophore) biosynthesis genes.

Microbial metabolites in the marine carbon cycle.

One-quarter of photosynthesis-derived carbon on Earth rapidly cycles through a set of short-lived seawater metabolites that are generated from the activities of marine phytoplankton, bacteria, grazers and viruses. Here we discuss the sources of microbial metabolites in the surface ocean, their roles in ecology and biogeochemistry, and approaches that can be used to analyse them from chemistry, biology, modelling and data science. Although microbial-derived metabolites account for only a minor fraction of the total reservoir of marine dissolved organic carbon, their flux and fate underpins the central role of the ocean in sustaining life on Earth.

Metabolic versatility of the nitrite-oxidizing bacterium Nitrospira marina and its proteomic response to oxygen-limited conditions.

The genus Nitrospira is the most widespread group of nitrite-oxidizing bacteria and thrives in diverse natural and engineered ecosystems. Nitrospira marina Nb-295T was isolated from the ocean over 30 years ago; however, its genome has not yet been analyzed. Here, we investigated the metabolic potential of N. marina based on its complete genome sequence and performed physiological experiments to test genome-derived hypotheses. Our data confirm that N. marina benefits from additions of undefined organic carbon substrates, has adaptations to resist oxidative, osmotic, and UV light-induced stress and low dissolved pCO2, and requires exogenous vitamin B12. In addition, N. marina is able to grow chemoorganotrophically on formate, and is thus not an obligate chemolithoautotroph. We further investigated the proteomic response of N. marina to low (∼5.6 µM) O2 concentrations. The abundance of a potentially more efficient CO2-fixing pyruvate:ferredoxin oxidoreductase (POR) complex and a high-affinity cbb3-type terminal oxidase increased under O2 limitation, suggesting a role in sustaining nitrite oxidation-driven autotrophy. This putatively more O2-sensitive POR complex might be protected from oxidative damage by Cu/Zn-binding superoxide dismutase, which also increased in abundance under low O2 conditions. Furthermore, the upregulation of proteins involved in alternative energy metabolisms, including Group 3b [NiFe] hydrogenase and formate dehydrogenase, indicate a high metabolic versatility to survive conditions unfavorable for aerobic nitrite oxidation. In summary, the genome and proteome of the first marine Nitrospira isolate identifies adaptations to life in the oxic ocean and provides insights into the metabolic diversity and niche differentiation of NOB in marine environments.

Single-Cell Transcriptomics of Abedinium Reveals a New Early-Branching Dinoflagellate Lineage.

Dinoflagellates possess many cellular characteristics with unresolved evolutionary histories. These include nuclei with greatly expanded genomes and chromatin packaged using histone-like proteins and dinoflagellate-viral nucleoproteins instead of histones, highly reduced mitochondrial genomes with extensive RNA editing, a mix of photosynthetic and cryptic secondary plastids, and tertiary plastids. Resolving the evolutionary origin of these traits requires understanding their ancestral states and early intermediates. Several early-branching dinoflagellate lineages are good candidates for such reconstruction, however these cells tend to be delicate and environmentally sparse, complicating such analyses. Here, we employ transcriptome sequencing from manually isolated and microscopically documented cells to resolve the placement of two cells of one such genus, Abedinium, collected by remotely operated vehicle in deep waters off the coast of Monterey Bay, CA. One cell corresponds to the only described species, Abedinium dasypus, whereas the second cell is distinct and formally described as Abedinium folium, sp. nov. Abedinium has classically been assigned to the early-branching dinoflagellate subgroup Noctilucales, which is weakly supported by phylogenetic analyses of small subunit ribosomal RNA, the single characterized gene from any member of the order. However, an analysis based on 221 proteins from the transcriptome places Abedinium as a distinct lineage, separate from and basal to Noctilucales and the rest of the core dinoflagellates. The transcriptome also contains evidence of a cryptic plastid functioning in the biosynthesis of isoprenoids, iron-sulfur clusters, and heme, a mitochondrial genome with all three expected protein-coding genes (cob, cox1, and cox3), and the presence of some but not all dinoflagellate-specific chromatin packaging proteins.

Microbial communities can predict the ecological condition of headwater streams.

Humanity's reliance on clean water and the ecosystem services provided makes identifying efficient and effective ways to assess the ecological condition of streams ever more important. We used high throughput sequencing of the 16S rRNA region to explore relationships between stream microbial communities, environmental attributes, and assessments of stream ecological condition. Bacteria and archaea in microbial community samples collected from the water column and from stream sediments during spring and summer were used to replicate standard assessments of ecological condition performed with benthic macroinvertebrate collections via the Benthic Index of Biotic Integrity (BIBI). Microbe-based condition assessments were generated at different levels of taxonomic resolution from phylum to OTU (Operational Taxonomic Units) in order to understand appropriate levels of taxonomic aggregation. Stream sediment microbial communities from both spring and summer were much better than the water column at replicating BIBI condition assessment results. Accuracies were as high as 100% on training data used to build the models and up to 80% on validation data used to assess predictions. Assessments using all OTUs usually had the highest accuracy on training data, but were lower on validation data due to overfitting. In contrast, assessments at the order-level had similar performance accuracy for validation data, and a reduced subset of orders also performed well, suggesting the method could be generalized to other watersheds. Subsets of the important orders responded similarly to environmental gradients compared to the entire community, where strong shifts in community structure occurred for known aquatic stressors such as pH, dissolved organic carbon, and nitrate nitrogen. The results suggest the stream microbes may be useful for assessing the ecological condition of streams and especially useful for stream restorations where many eukaryotic taxa have been eliminated due to prior degradation and are unable to recolonize.

Heterotrophic Thaumarchaea with Small Genomes Are Widespread in the Dark Ocean.

The Thaumarchaeota is a diverse archaeal phylum comprising numerous lineages that play key roles in global biogeochemical cycling, particularly in the ocean. To date, all genomically characterized marine thaumarchaea are reported to be chemolithoautotrophic ammonia oxidizers. In this study, we report a group of putatively heterotrophic marine thaumarchaea (HMT) with small genome sizes that is globally abundant in the mesopelagic, apparently lacking the ability to oxidize ammonia. We assembled five HMT genomes from metagenomic data and show that they form a deeply branching sister lineage to the ammonia-oxidizing archaea (AOA). We identify this group in metagenomes from mesopelagic waters in all major ocean basins, with abundances reaching up to 6% of that of AOA. Surprisingly, we predict the HMT have small genomes of ∼1 Mbp, and our ancestral state reconstruction indicates this lineage has undergone substantial genome reduction compared to other related archaea. The genomic repertoire of HMT indicates a versatile metabolism for aerobic chemoorganoheterotrophy that includes a divergent form III-a RuBisCO, a 2M respiratory complex I that has been hypothesized to increase energetic efficiency, and a three-subunit heme-copper oxidase complex IV that is absent from AOA. We also identify 21 pyrroloquinoline quinone (PQQ)-dependent dehydrogenases that are predicted to supply reducing equivalents to the electron transport chain and are among the most highly expressed HMT genes, suggesting these enzymes play an important role in the physiology of this group. Our results suggest that heterotrophic members of the Thaumarchaeota are widespread in the ocean and potentially play key roles in global chemical transformations.IMPORTANCE It has been known for many years that marine Thaumarchaeota are abundant constituents of dark ocean microbial communities, where their ability to couple ammonia oxidation and carbon fixation plays a critical role in nutrient dynamics. In this study, we describe an abundant group of putatively heterotrophic marine Thaumarchaeota (HMT) in the ocean with physiology distinct from those of their ammonia-oxidizing relatives. HMT lack the ability to oxidize ammonia and fix carbon via the 3-hydroxypropionate/4-hydroxybutyrate pathway but instead encode a form III-a RuBisCO and diverse PQQ-dependent dehydrogenases that are likely used to conserve energy in the dark ocean. Our work expands the scope of known diversity of Thaumarchaeota in the ocean and provides important insight into a widespread marine lineage.

Global reconstruction reduces the uncertainty of oceanic nitrous oxide emissions and reveals a vigorous seasonal cycle.

Assessment of the global budget of the greenhouse gas nitrous oxide ([Formula: see text]O) is limited by poor knowledge of the oceanic [Formula: see text]O flux to the atmosphere, of which the magnitude, spatial distribution, and temporal variability remain highly uncertain. Here, we reconstruct climatological [Formula: see text]O emissions from the ocean by training a supervised learning algorithm with over 158,000 [Formula: see text]O measurements from the surface ocean-the largest synthesis to date. The reconstruction captures observed latitudinal gradients and coastal hot spots of [Formula: see text]O flux and reveals a vigorous global seasonal cycle. We estimate an annual mean [Formula: see text]O flux of 4.2 ± 1.0 Tg N[Formula: see text], 64% of which occurs in the tropics, and 20% in coastal upwelling systems that occupy less than 3% of the ocean area. This [Formula: see text]O flux ranges from a low of 3.3 ± 1.3 Tg N[Formula: see text] in the boreal spring to a high of 5.5 ± 2.0 Tg N[Formula: see text] in the boreal summer. Much of the seasonal variations in global [Formula: see text]O emissions can be traced to seasonal upwelling in the tropical ocean and winter mixing in the Southern Ocean. The dominant contribution to seasonality by productive, low-oxygen tropical upwelling systems (>75%) suggests a sensitivity of the global [Formula: see text]O flux to El Niño-Southern Oscillation and anthropogenic stratification of the low latitude ocean. This ocean flux estimate is consistent with the range adopted by the Intergovernmental Panel on Climate Change, but reduces its uncertainty by more than fivefold, enabling more precise determination of other terms in the atmospheric [Formula: see text]O budget.

Diversity, ecology and evolution of Archaea.

Compared to bacteria, our knowledge of archaeal biology is limited. Historically, microbiologists have mostly relied on culturing and single-gene diversity surveys to understand Archaea in nature. However, only six of the 27 currently proposed archaeal phyla have cultured representatives. Advances in genomic sequencing and computational approaches are revolutionizing our understanding of Archaea. The recovery of genomes belonging to uncultured groups from the environment has resulted in the description of several new phyla, many of which are globally distributed and are among the predominant organisms on the planet. In this Review, we discuss how these genomes, together with long-term enrichment studies and elegant in situ measurements, are providing insights into the metabolic capabilities of the Archaea. We also debate how such studies reveal how important Archaea are in mediating an array of ecological processes, including global carbon and nutrient cycles, and how this increase in archaeal diversity has expanded our view of the tree of life and early archaeal evolution, and has provided new insights into the origin of eukaryotes.

Author Correction: Diversity, ecology and evolution of Archaea.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Headwater Stream Microbial Diversity and Function across Agricultural and Urban Land Use Gradients.

Anthropogenic activity impacts stream ecosystems, resulting in a loss of diversity and ecosystem function; however, little is known about the response of aquatic microbial communities to changes in land use. Here, microbial communities were characterized in 82 headwater streams across a gradient of urban and agricultural land uses using 16S rRNA gene amplicon sequencing and compared to a rich data set of physicochemical variables and traditional benthic invertebrate indicators. Microbial diversity and community structures differed among watersheds with high agricultural, urban, and forested land uses, and community structure differed in streams classified as being in good, fair, poor, and very poor condition using benthic invertebrate indicators. Microbial community similarity decayed with geodesic distance across the study region but not with environmental distance. Stream community respiration rates ranged from 21.7 to 1,570 mg O2 m-2 day-1 and 31.9 to 3,670 mg O2 m-2 day-1 for water column and sediments, respectively, and correlated with nutrients associated with anthropogenic influence and microbial community structure. Nitrous oxide (N2O) concentrations ranged from 0.22 to 4.41 µg N2O liter-1; N2O concentration was negatively correlated with forested land use and was positively correlated with dissolved inorganic nitrogen concentrations. Our findings suggest that stream microbial communities are impacted by watershed land use and can potentially be used to assess ecosystem health.IMPORTANCE Stream ecosystems are frequently impacted by changes in watershed land use, resulting in altered hydrology, increased pollutant and nutrient loads, and habitat degradation. Macroinvertebrates and fish are strongly affected by changes in stream conditions and are commonly used in biotic indices to assess ecosystem health. Similarly, microbes respond to environmental stressors, and changes in community composition alter key ecosystem processes. The response of microbes to habitat degradation and their role in global biogeochemical cycles provide an opportunity to use microbes as a monitoring tool. Here, we identify stream microbes that respond to watershed urbanization and agricultural development and demonstrate that microbial diversity and community structure can be used to assess stream conditions and ecosystem functioning.

Controlled sampling of ribosomally active protistan diversity in sediment-surface layers identifies putative players in the marine carbon sink.

Marine sediments are one of the largest carbon reservoir on Earth, yet the microbial communities, especially the eukaryotes, that drive these ecosystems are poorly characterised. Here, we report implementation of a sampling system that enables injection of reagents into sediments at depth, allowing for preservation of RNA in situ. Using the RNA templates recovered, we investigate the 'ribosomally active' eukaryotic diversity present in sediments close to the water/sediment interface. We demonstrate that in situ preservation leads to recovery of a significantly altered community profile. Using SSU rRNA amplicon sequencing, we investigated the community structure in these environments, demonstrating a wide diversity and high relative abundance of stramenopiles and alveolates, specifically: Bacillariophyta (diatoms), labyrinthulomycetes and ciliates. The identification of abundant diatom rRNA molecules is consistent with microscopy-based studies, but demonstrates that these algae can also be exported to the sediment as active cells as opposed to dead forms. We also observe many groups that include, or branch close to, osmotrophic-saprotrophic protists (e.g. labyrinthulomycetes and Pseudofungi), microbes likely to be important for detrital decomposition. The sequence data also included a diversity of abundant amplicon-types that branch close to the Fonticula slime moulds. Taken together, our data identifies additional roles for eukaryotic microbes in the marine carbon cycle; where putative osmotrophic-saprotrophic protists represent a significant active microbial-constituent of the upper sediment layer.

Unexpected mitochondrial genome diversity revealed by targeted single-cell genomics of heterotrophic flagellated protists.

Most eukaryotic microbial diversity is uncultivated, under-studied and lacks nuclear genome data. Mitochondrial genome sampling is more comprehensive, but many phylogenetically important groups remain unsampled. Here, using a single-cell sorting approach combining tubulin-specific labelling with photopigment exclusion, we sorted flagellated heterotrophic unicellular eukaryotes from Pacific Ocean samples. We recovered 206 single amplified genomes, predominantly from underrepresented branches on the tree of life. Seventy single amplified genomes contained unique mitochondrial contigs, including 21 complete or near-complete mitochondrial genomes from formerly under-sampled phylogenetic branches, including telonemids, katablepharids, cercozoans and marine stramenopiles, effectively doubling the number of available samples of heterotrophic flagellate mitochondrial genomes. Collectively, these data identify a dynamic history of mitochondrial genome evolution including intron gain and loss, extensive patterns of genetic code variation and complex patterns of gene loss. Surprisingly, we found that stramenopile mitochondrial content is highly plastic, resembling patterns of variation previously observed only in plants.