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.

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.

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.

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.

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.

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 cbb₃- 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.

BinSanity: unsupervised clustering of environmental microbial assemblies using coverage and affinity propagation

Metagenomics has become an integral part of defining microbial diversity in various environments. Many ecosystems have characteristically low biomass and few cultured representatives. Linking potential metabolisms to phylogeny in environmental microorganisms is important for interpreting microbial community functions and the impacts these communities have on geochemical cycles. However, with metagenomic studies there is the computational hurdle of ‘binning’ contigs into phylogenetically related units or putative genomes. Binning methods have been implemented with varying approaches such as k-means clustering, Gaussian mixture models, hierarchical clustering, neural networks, and two-way clustering; however, many of these suffer from biases against low coverage/abundance organisms and closely related taxa/strains. We are introducing a new binning method, BinSanity, that utilizes the clustering algorithm affinity propagation (AP), to cluster assemblies using coverage with compositional based refinement (tetranucleotide frequency and percent GC content) to optimize bins containing multiple source organisms. This separation of composition and coverage based clustering reduces bias for closely related taxa. BinSanity was developed and tested on artificial metagenomes varying in size and complexity. Results indicate that BinSanity has a higher precision, recall, and Adjusted Rand Index compared to five commonly implemented methods. When tested on a previously published environmental metagenome, BinSanity generated high completion and low redundancy bins corresponding with the published metagenome-assembled genomes.

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

The rock-hosted, oceanic crustal aquifer is one of the largest ecosystems on Earth, yet little is known about its indigenous microorganisms. Here we provide the first phylogenetic and functional description of an active microbial community residing in the cold oxic crustal aquifer. Using subseafloor observatories, we recovered crustal fluids and found that the geochemical composition is similar to bottom seawater, as are cell abundances. However, based on relative abundances and functional potential of key bacterial groups, the crustal fluid microbial community is heterogeneous and markedly distinct from seawater. Potential rates of autotrophy and heterotrophy in the crust exceeded those of seawater, especially at elevated temperatures (25 °C) and deeper in the crust. Together, these results reveal an active, distinct, and diverse bacterial community engaged in both heterotrophy and autotrophy in the oxygenated crustal aquifer, providing key insight into the role of microbial communities in the ubiquitous cold dark subseafloor biosphere.

Comparative genomics of planktonic Flavobacteriaceae from the Gulf of Maine using metagenomic data

BACKGROUND

The Gulf of Maine is an important biological province of the Northwest Atlantic with high productivity year round. From an environmental Sanger-based metagenome, sampled in summer and winter, we were able to assemble and explore the partial environmental genomes of uncultured members of the class Flavobacteria. Each of the environmental genomes represents organisms that compose less than 1% of the total microbial metagenome.

RESULTS

Four partial environmental genomes were assembled with varying degrees of estimated completeness (37%-84% complete) and were analyzed from a perspective of gathering information regarding niche partitioning between co-occurring organisms. Comparative genomics revealed potentially important niche partitioning genomic variations, including iron transporters and genes associated with cell attachment and polymer degradation. Analysis of large syntenic regions helped reveal potentially ecologically relevant variations for Flavobacteriaceae in the Gulf of Maine, such as arginine biosynthesis, and identify a putative genomic island incorporating novel exogenous genes from the environment.

CONCLUSIONS

Biogeographic analysis revealed flavobacteria species with distinct abundance patterns suggesting the presence of local blooms relative to the other species, as well as seasonally selected organisms. The analysis of genomic content for the Gulf of Maine Flavobacteria supports the hypothesis of a particle-associated lifestyle and specifically highlights a number of putative coding sequences that may play a role in the remineralization of particulate organic matter. And lastly, analysis of the underlying sequences for each assembled genome revealed seasonal and nonseasonal variants of specific genes implicating a dynamic interaction between individuals within the species.

Microbial communities associated with ferromanganese nodules and the surrounding sediments

The formation and maintenance of deep-sea ferromanganese/polymetallic nodules still remains a mystery 140 years after their discovery. The wealth of rare metals concentrated in these nodules has spurred global interest in exploring the mining potential of these resources. The prevailing theory of abiotic formation has been called into question and the role of microbial metabolisms in nodule development is now an area of active research. To understand the community structure of microbes associated with nodules and their surrounding sediment, we performed targeted sequencing of the V4 hypervariable region of the 16S rRNA gene from three nodules collected from the central South Pacific. Results have shown that the microbial communities of the nodules are significantly distinct from the communities in the surrounding sediments, and that the interiors of the nodules harbor communities different from the exterior. This suggests not only differences in potential metabolisms between the nodule and sediment communities, but also differences in the dominant metabolisms of interior and exterior communities. We identified several operational taxonomic units (OTUs) unique to both the nodule and sediment environments. The identified OTUs were assigned putative taxonomic identifications, including two OTUs only found associated with the nodules, which were assigned to the α-Proteobacteria. Finally, we explored the diversity of the most assigned taxonomic group, the Thaumarchaea MG-1, which revealed novel OTUs compared to previous research from the region and suggests a potential role as a source of fixed carbon for ammonia oxidizing archaea in the environment.

Prospects for the study of evolution in the deep biosphere

Since the days of Darwin, scientists have used the framework of the theory of evolution to explore the interconnectedness of life on Earth and adaptation of organisms to the ever-changing environment. The advent of molecular biology has advanced and accelerated the study of evolution by allowing direct examination of the genetic material that ultimately determines the phenotypes upon which selection acts. The study of evolution has been furthered through examination of microbial evolution, with large population numbers, short generation times, and easily extractable DNA. Such work has spawned the study of microbial biogeography, with the realization that concepts developed in population genetics may be applicable to microbial genomes (Martiny et al., 2006; Manhes and Velicer, 2011). Microbial biogeography and adaptation has been examined in many different environments. Here we argue that the deep biosphere is a unique environment for the study of evolution and list specific factors that can be considered and where the studies may be performed. This publication is the result of the NSF-funded Center for Dark Energy Biosphere Investigations (C-DEBI) theme team on Evolution (www.darkenergybiosphere.org).

Metagenomic analysis of a complex marine planktonic thaumarchaeal community from the Gulf of Maine

Thaumarchaea, which represent as much as 20% of prokaryotic biomass in the open ocean, have been linked to environmentally relevant biogeochemical processes, such as ammonia oxidation (nitrification) and inorganic carbon fixation. We have used culture-independent methods to study this group because current cultivation limitations have proved a hindrance in studying these organisms. From a metagenomic data set obtained from surface waters from the Gulf of Maine, we have identified 36,111 sequence reads (containing 30 Mbp) likely derived from environmental planktonic Thaumarchaea. Metabolic analysis of the raw sequences and assemblies identified copies of the catalytic subunit required in aerobic ammonia oxidation. In addition, genes that comprise a nearly complete carbon assimilation pathway in the form of the 3-hyroxypropionate/4-hydroxybutyrate cycle were identified. Comparative genomics contrasting the putative environmental thaumarchaeal sequences and 'Candidatus Nitrosopumilus maritimus SCM1' revealed a number of genomic islands absent in the Gulf of Maine population. Analysis of these genomic islands revealed an integrase-associated island also found in distantly related microbial species, variations in the abundance of genes predicted to be important in thaumarchaeal respiratory chain, and the absence of a high-affinity phosphate uptake operon. Analysis of the underlying sequence diversity suggests the presence of at least two dominant environmental populations. Attempts to assemble complete environmental genomes were unsuccessful, but analysis of scaffolds revealed two diverging populations, including a thaumarchaeal-related scaffold with the full urease operon. Ultimately, the analysis revealed a number of insights into the metabolic potential of a predominantly uncultivated lineage of organisms. The predicted functions in the thaumarchaeal metagenomic sequences are directly supported by historic measurements of nutrient concentrations and provide new avenues of research in regards to understanding the role Thaumarchaea play in the environment.