Increasing transposase abundance with ocean depth correlates with a particle-associated lifestyle.

Transposases are mobile genetic elements that move within and between genomes, promoting genomic plasticity in microorganisms. In marine microbial communities, the abundance of transposases increases with depth, but the reasons behind this trend remain unclear. Our analysis of metagenomes from the Tara Oceans and Malaspina Expeditions suggests that a particle-associated lifestyle is the main covariate for the high occurrence of transposases in the deep ocean, and this trend holds true for individual genomes as well as in a community-wide sense. We observed a strong and depth-independent correlation between transposase abundance and the presence of biofilm-associated genes, as well as the prevalence of secretory enzymes. This suggests that mobile genetic elements readily propagate among microbial communities within crowded biofilms. Furthermore, we show that particle association positively correlates with larger genome size, which is in turn associated with higher transposase abundance. Cassette sequences associated with transposons are enriched with genes related to defense mechanisms, which are more highly expressed in the deep sea. Thus, while transposons spread at the expense of their microbial hosts, they also introduce novel genes and potentially benefit the hosts in helping to compete for limited resources. Overall, our results suggest a new understanding of deep ocean particles as highways for gene sharing among defensively oriented microbial genomes.IMPORTANCEGenes can move within and between microbial genomes via mobile genetic elements, which include transposases and transposons. In the oceans, there is a puzzling increase in transposase abundance in microbial genomes as depth increases. To gain insight into this trend, we conducted an extensive analysis of marine microbial metagenomes and metatranscriptomes. We found a significant correlation between transposase abundance and a particle-associated lifestyle among marine microbes at both the metagenome and genome-resolved levels. We also observed a link between transposase abundance and genes related to defense mechanisms. These results suggest that as microbes become densely packed into crowded particles, mobile genes are more likely to spread and carry genetic material that provides a competitive advantage in crowded habitats. This may enable deep sea microbes to effectively compete in such environments.

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.

Identification of microbial metabolic functional guilds from large genomic datasets.

Heterotrophic microbes play an important role in the Earth System as key drivers of major biogeochemical cycles. Specifically, the consumption rate of organic matter is set by the interaction between diverse microbial communities and the chemical and physical environment in which they reside. Modeling these dynamics requires reducing the complexity of microbial communities and linking directly with biogeochemical functions. Microbial metabolic functional guilds provide one approach for reducing microbial complexity and incorporating microbial biogeochemical functions into models. However, we lack a way to identify these guilds. In this study, we present a method for defining metabolic functional guilds from annotated genomes, which are derived from both uncultured and cultured organisms. This method utilizes an Aspect Bernoulli (AB) model and was tested on three large genomic datasets with 1,733-3,840 genomes each. Ecologically relevant microbial metabolic functional guilds were identified including guilds related to DMSP degradation, dissimilatory nitrate reduction to ammonia, and motile copiotrophy. This method presents a way to generate hypotheses about functions co-occurring within individual microbes without relying on cultured representatives. Applying the concept of metabolic functional guilds to environmental samples will provide new insight into the role that heterotrophic microbial communities play in setting rates of carbon cycling.

Microbial Populations Are Shaped by Dispersal and Recombination in a Low Biomass Subseafloor Habitat.

The subseafloor is a vast habitat that supports microorganisms that have a global scale impact on geochemical cycles. Many of the endemic microbial communities inhabiting the subseafloor consist of small populations under growth-limited conditions. For small populations, stochastic evolutionary events can have large impacts on intraspecific population dynamics and allele frequencies. These conditions are fundamentally different from those experienced by most microorganisms in surface environments, and it is unknown how small population sizes and growth-limiting conditions influence evolution and population structure in the subsurface. Using a 2-year, high-resolution environmental time series, we examine the dynamics of microbial populations from cold, oxic crustal fluids collected from the subseafloor site North Pond, located near the mid-Atlantic ridge. Our results reveal rapid shifts in overall abundance, allele frequency, and strain abundance across the time points observed, with evidence for homologous recombination between coexisting lineages. We show that the subseafloor aquifer is a dynamic habitat that hosts microbial metapopulations that disperse frequently through the crustal fluids, enabling gene flow and recombination between microbial populations. The dynamism and stochasticity of microbial population dynamics in North Pond suggest that these forces are important drivers in the evolution of microbial populations in the vast subseafloor habitat. IMPORTANCE The cold, oxic subseafloor is an understudied habitat that is difficult to access, yet important to global biogeochemical cycles and starkly different compared to microbial habitats on the surface of the Earth. Our understanding of microbial evolution and population dynamics is largely molded by studies of microbes living in surface habitats that can host 10 to 1,000 times more microbial biomass than is frequently observed in the subsurface. This study provides an opportunity to observe population dynamics within a low biomass, growth-limited environment and reveals that microbial populations in the subseafloor are influenced by changes in selection pressure and gene sweeps. In addition, recombination between strains that have dispersed from elsewhere within the aquifer has an important impact on the evolution of microbial populations. Much of the microbial life on the planet exists under growth-limited conditions, and the subseafloor provides a natural laboratory to explore how life evolves in such environments.

Petrobactin, a siderophore produced by Alteromonas, mediates community iron acquisition in the global ocean.

It is now widely accepted that siderophores play a role in marine iron biogeochemical cycling. However, the mechanisms by which siderophores affect the availability of iron from specific sources and the resulting significance of these processes on iron biogeochemical cycling as a whole have remained largely untested. In this study, we develop a model system for testing the effects of siderophore production on iron bioavailability using the marine copiotroph Alteromonas macleodii ATCC 27126. Through the generation of the knockout cell line ΔasbB::kmr, which lacks siderophore biosynthetic capabilities, we demonstrate that the production of the siderophore petrobactin enables the acquisition of iron from mineral sources and weaker iron-ligand complexes. Notably, the utilization of lithogenic iron, such as that from atmospheric dust, indicates a significant role for siderophores in the incorporation of new iron into marine systems. We have also detected petrobactin, a photoreactive siderophore, directly from seawater in the mid-latitudes of the North Pacific and have identified the biosynthetic pathway for petrobactin in bacterial metagenome-assembled genomes widely distributed across the global ocean. Together, these results improve our mechanistic understanding of the role of siderophore production in iron biogeochemical cycling in the marine environment wherein iron speciation, bioavailability, and residence time can be directly influenced by microbial activities.

Time-series transcriptomics from cold, oxic subseafloor crustal fluids reveals a motile, mixotrophic microbial community.

The oceanic crustal aquifer is one of the largest habitable volumes on Earth, and it harbors a reservoir of microbial life that influences global-scale biogeochemical cycles. Here, we use time series metagenomic and metatranscriptomic data from a low-temperature, ridge flank environment representative of the majority of global hydrothermal fluid circulation in the ocean to reconstruct microbial metabolic potential, transcript abundance, and community dynamics. We also present metagenome-assembled genomes from recently collected fluids that are furthest removed from drilling disturbances. Our results suggest that the microbial community in the North Pond aquifer plays an important role in the oxidation of organic carbon within the crust. This community is motile and metabolically flexible, with the ability to use both autotrophic and organotrophic pathways, as well as function under low oxygen conditions by using alternative electron acceptors such as nitrate and thiosulfate. Anaerobic processes are most abundant in subseafloor horizons deepest in the aquifer, furthest from connectivity with the deep ocean, and there was little overlap in the active microbial populations between sampling horizons. This work highlights the heterogeneity of microbial life in the subseafloor aquifer and provides new insights into biogeochemical cycling in ocean crust.

Marine Dadabacteria exhibit genome streamlining and phototrophy-driven niche partitioning.

The remineralization of organic material via heterotrophy in the marine environment is performed by a diverse and varied group of microorganisms that can specialize in the type of organic material degraded and the niche they occupy. The marine Dadabacteria are cosmopolitan in the marine environment and belong to a candidate phylum for which there has not been a comprehensive assessment of the available genomic data to date. Here in, we assess the functional potential of the marine pelagic Dadabacteria in comparison to members of the phylum that originate from terrestrial, hydrothermal, and subsurface environments. Our analysis reveals that the marine pelagic Dadabacteria have streamlined genomes, corresponding to smaller genome sizes and lower nitrogen content of their DNA and predicted proteome, relative to their phylogenetic counterparts. Collectively, the Dadabacteria have the potential to degrade microbial dissolved organic matter, specifically peptidoglycan and phospholipids. The marine Dadabacteria belong to two clades with apparent distinct ecological niches in global metagenomic data: a clade with the potential for photoheterotrophy through the use of proteorhodopsin, present predominantly in surface waters up to 100 m depth; and a clade lacking the potential for photoheterotrophy that is more abundant in the deep photic zone.

NCBI's Virus Discovery Codeathon: Building "FIVE" -The Federated Index of Viral Experiments API Index.

Viruses represent important test cases for data federation due to their genome size and the rapid increase in sequence data in publicly available databases. However, some consequences of previously decentralized (unfederated) data are lack of consensus or comparisons between feature annotations. Unifying or displaying alternative annotations should be a priority both for communities with robust entry representation and for nascent communities with burgeoning data sources. To this end, during this three-day continuation of the Virus Hunting Toolkit codeathon series (VHT-2), a new integrated and federated viral index was elaborated. This Federated Index of Viral Experiments (FIVE) integrates pre-existing and novel functional and taxonomy annotations and virus-host pairings. Variability in the context of viral genomic diversity is often overlooked in virus databases. As a proof-of-concept, FIVE was the first attempt to include viral genome variation for HIV, the most well-studied human pathogen, through viral genome diversity graphs. As per the publication of this manuscript, FIVE is the first implementation of a virus-specific federated index of such scope. FIVE is coded in BigQuery for optimal access of large quantities of data and is publicly accessible. Many projects of database or index federation fail to provide easier alternatives to access or query information. To this end, a Python API query system was developed to enhance the accessibility of FIVE.

Biases in genome reconstruction from metagenomic data.

BACKGROUND: Advances in sequencing, assembly, and assortment of contigs into species-specific bins has enabled the reconstruction of genomes from metagenomic data (MAGs). Though a powerful technique, it is difficult to determine whether assembly and binning techniques are accurate when applied to environmental metagenomes due to a lack of complete reference genome sequences against which to check the resulting MAGs. METHODS: We compared MAGs derived from an enrichment culture containing ~20 organisms to complete genome sequences of 10 organisms isolated from the enrichment culture. Factors commonly considered in binning software-nucleotide composition and sequence repetitiveness-were calculated for both the correctly binned and not-binned regions. This direct comparison revealed biases in sequence characteristics and gene content in the not-binned regions. Additionally, the composition of three public data sets representing MAGs reconstructed from the Tara Oceans metagenomic data was compared to a set of representative genomes available through NCBI RefSeq to verify that the biases identified were observable in more complex data sets and using three contemporary binning software packages. RESULTS: Repeat sequences were frequently not binned in the genome reconstruction processes, as were sequence regions with variant nucleotide composition. Genes encoded on the not-binned regions were strongly biased towards ribosomal RNAs, transfer RNAs, mobile element functions and genes of unknown function. Our results support genome reconstruction as a robust process and suggest that reconstructions determined to be >90% complete are likely to effectively represent organismal function; however, population-level genotypic heterogeneity in natural populations, such as uneven distribution of plasmids, can lead to incorrect inferences.

MetaSanity: an integrated microbial genome evaluation and annotation pipeline.

SUMMARY: As the importance of microbiome research continues to become more prevalent and essential to understanding a wide variety of ecosystems (e.g. marine, built, host associated, etc.), there is a need for researchers to be able to perform highly reproducible and quality analysis of microbial genomes. MetaSanity incorporates analyses from 11 existing and widely used genome evaluation and annotation suites into a single, distributable workflow, thereby decreasing the workload of microbiologists by allowing for a flexible, expansive data analysis pipeline. MetaSanity has been designed to provide separate, reproducible workflows that (i) can determine the overall quality of a microbial genome, while providing a putative phylogenetic assignment, and (ii) can assign structural and functional gene annotations with varying degrees of specificity to suit the needs of the researcher. The software suite combines the results from several tools to provide broad insights into overall metabolic function. Importantly, this software provides built-in optimization for 'big data' analysis by storing all relevant outputs in an SQL database, allowing users to query all the results for the elements that will most impact their research. AVAILABILITY AND IMPLEMENTATION: MetaSanity is provided under the GNU General Public License v.3.0 and is available for download at https://github.com/cjneely10/MetaSanity. This application is distributed as a Docker image. MetaSanity is implemented in Python3/Cython and C++. Instructions for its installation and use are available within the GitHub wiki page at https://github.com/cjneely10/MetaSanity/wiki, and additional instructions are available at https://cjneely10.github.io/year-archive/. MetaSanity is optimized for users with limited programing experience. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Author Correction: A distinct and active bacterial community in cold oxygenated fluids circulating beneath the western flank of the Mid-Atlantic ridge.

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

NCBI's Virus Discovery Hackathon: Engaging Research Communities to Identify Cloud Infrastructure Requirements.

A wealth of viral data sits untapped in publicly available metagenomic data sets when it might be extracted to create a usable index for the virological research community. We hypothesized that work of this complexity and scale could be done in a hackathon setting. Ten teams comprised of over 40 participants from six countries, assembled to create a crowd-sourced set of analysis and processing pipelines for a complex biological data set in a three-day event on the San Diego State University campus starting 9 January 2019. Prior to the hackathon, 141,676 metagenomic data sets from the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) were pre-assembled into contiguous assemblies (contigs) by NCBI staff. During the hackathon, a subset consisting of 2953 SRA data sets (approximately 55 million contigs) was selected, which were further filtered for a minimal length of 1 kb. This resulted in 4.2 million (Mio) contigs, which were aligned using BLAST against all known virus genomes, phylogenetically clustered and assigned metadata. Out of the 4.2 Mio contigs, 360,000 contigs were labeled with domains and an additional subset containing 4400 contigs was screened for virus or virus-like genes. The work yielded valuable insights into both SRA data and the cloud infrastructure required to support such efforts, revealing analysis bottlenecks and possible workarounds thereof. Mainly: (i) Conservative assemblies of SRA data improves initial analysis steps; (ii) existing bioinformatic software with weak multithreading/multicore support can be elevated by wrapper scripts to use all cores within a computing node; (iii) redesigning existing bioinformatic algorithms for a cloud infrastructure to facilitate its use for a wider audience; and (iv) a cloud infrastructure allows a diverse group of researchers to collaborate effectively. The scientific findings will be extended during a follow-up event. Here, we present the applied workflows, initial results, and lessons learned from the hackathon.

Genome Sequence of Mariprofundus sp. Strain EBB-1, a Novel Marine Autotroph Isolated from an Iron-Sulfur Mineral.

Mariprofundus sp. strain EBB-1 was isolated from a pyrrhotite biofilm incubated in seawater from East Boothbay (ME, USA). Strain EBB-1 is an autotrophic member of the class Zetaproteobacteria with the ability to form iron oxide biominerals. Here, we present the 2.88-Mb genome sequence of EBB-1, which contains 2,656 putative protein-coding sequences.

Metabolic diversity within the globally abundant Marine Group II Euryarchaea offers insight into ecological patterns.

Despite their discovery over 25 years ago, the Marine Group II Euryarchaea (MGII) remain a difficult group of organisms to study, lacking cultured isolates and genome references. The MGII have been identified in marine samples from around the world, and evidence supports a photoheterotrophic lifestyle combining phototrophy via proteorhodopsins with the remineralization of high molecular weight organic matter. Divided between two clades, the MGII have distinct ecological patterns that are not understood based on the limited number of available genomes. Here, I present a comparative genomic analysis of 250 MGII genomes, providing a comprehensive investigation of these mesophilic archaea. This analysis identifies 17 distinct subclades including nine subclades that previously lacked reference genomes. The metabolic potential and distribution of the MGII genera reveals distinct roles in the environment, identifying algal-saccharide-degrading coastal subclades, protein-degrading oligotrophic surface ocean subclades, and mesopelagic subclades lacking proteorhodopsins, common in all other subclades.

Genome Sequence of Geothermobacter sp. Strain HR-1, an Iron Reducer from the Lo'ihi Seamount, Hawai'i.

Geothermobacter sp. strain HR-1 was isolated from the Lo'ihi Seamount vent system in the Pacific Ocean at a depth of 1,000 m. Reported here is its 3.84-Mb genome sequence.

Genome Sequence of Hydrogenovibrio sp. Strain SC-1, a Chemolithoautotrophic Sulfur and Iron Oxidizer.

Hydrogenovibrio sp. strain SC-1 was isolated from pyrrhotite incubated in situ in the marine surface sediment of Catalina Island, CA. Strain SC-1 has demonstrated autotrophic growth through the oxidation of thiosulfate and iron. Here, we present the 2.45-Mb genome sequence of SC-1, which contains 2,262 protein-coding genes.

The reconstruction of 2,631 draft metagenome-assembled genomes from the global oceans.

Microorganisms play a crucial role in mediating global biogeochemical cycles in the marine environment. By reconstructing the genomes of environmental organisms through metagenomics, researchers are able to study the metabolic potential of Bacteria and Archaea that are resistant to isolation in the laboratory. Utilizing the large metagenomic dataset generated from 234 samples collected during the Tara Oceans circumnavigation expedition, we were able to assemble 102 billion paired-end reads into 562 million contigs, which in turn were co-assembled and consolidated in to 7.2 million contigs ≥2 kb in length. Approximately 1 million of these contigs were binned to reconstruct draft genomes. In total, 2,631 draft genomes with an estimated completion of ≥50% were generated (1,491 draft genomes >70% complete; 603 genomes >90% complete). A majority of the draft genomes were manually assigned phylogeny based on sets of concatenated phylogenetic marker genes and/or 16S rRNA gene sequences. The draft genomes are now publically available for the research community at-large.

A dynamic microbial community with high functional redundancy inhabits the cold, oxic subseafloor aquifer.

The rock-hosted subseafloor crustal aquifer harbors a reservoir of microbial life that may influence global marine biogeochemical cycles. Here we utilized metagenomic libraries of crustal fluid samples from North Pond, located on the flanks of the Mid-Atlantic Ridge, a site with cold, oxic subseafloor fluid circulation within the upper basement to query microbial diversity. Twenty-one samples were collected during a 2-year period to examine potential microbial metabolism and community dynamics. We observed minor changes in the geochemical signatures over the 2 years, yet the microbial community present in the crustal fluids underwent large shifts in the dominant taxonomic groups. An analysis of 195 metagenome-assembled genomes (MAGs) were generated from the data set and revealed a connection between litho- and autotrophic processes, linking carbon fixation to the oxidation of sulfide, sulfur, thiosulfate, hydrogen, and ferrous iron in members of the Proteobacteria, specifically the Alpha-, Gamma- and Zetaproteobacteria, the Epsilonbacteraeota and the Planctomycetes. Despite oxic conditions, analysis of the MAGs indicated that members of the microbial community were poised to exploit hypoxic or anoxic conditions through the use of microaerobic cytochromes, such as cbb3- and bd-type cytochromes, and alternative electron acceptors, like nitrate and sulfate. Temporal and spatial trends from the MAGs revealed a high degree of functional redundancy that did not correlate with the shifting microbial community membership, suggesting functional stability in mediating subseafloor biogeochemical cycles. Collectively, the repeated sampling at multiple sites, together with the successful binning of hundreds of genomes, provides an unprecedented data set for investigation of microbial communities in the cold, oxic crustal aquifer.

290 metagenome-assembled genomes from the Mediterranean Sea: a resource for marine microbiology.

The Tara Oceans Expedition has provided large, publicly-accessible microbial metagenomic datasets from a circumnavigation of the globe. Utilizing several size fractions from the samples originating in the Mediterranean Sea, we have used current assembly and binning techniques to reconstruct 290 putative draft metagenome-assembled bacterial and archaeal genomes, with an estimated completion of ≥50%, and an additional 2,786 bins, with estimated completion of 0-50%. We have submitted our results, including initial taxonomic and phylogenetic assignments, for the putative draft genomes to open-access repositories for the scientific community to use in ongoing research.

Genome Sequence of Geothermobacter sp. Strain EPR-M, a Deep-Sea Hydrothermal Vent Iron Reducer.

Geothermobacter sp. strain EPR-M was isolated from a hydrothermal vent on the East Pacific Rise and has been shown to participate in the reduction of Fe(III) oxides. Here, we report its 3.73-Mb draft genome sequence.