The Amazon continuum dataset: quantitative metagenomic and metatranscriptomic inventories of the Amazon River plume, June 2010

Microbial ecology of northern Gulf of Mexico estuarine waters

Estuarine and coastal ecosystems are of high economic and ecological importance, owing to their diverse communities and the disproportionate role they play in carbon cycling, particularly in carbon sequestration. Organisms inhabiting these environments must overcome strong natural fluctuations in salinity, nutrients, and turbidity, as well as numerous climate change-induced disturbances such as land loss, sea level rise, and, in some locations, increasingly severe tropical cyclones that threaten to disrupt future ecosystem health. The northern Gulf of Mexico (nGoM) along the Louisiana coast contains dozens of estuaries, including the Mississippi-Atchafalaya River outflow, which dramatically influence the region due to their vast upstream watershed. Nevertheless, the microbiology of these estuaries and surrounding coastal environments has received little attention. To improve our understanding of microbial ecology in the understudied coastal nGoM, we conducted a 16S rRNA gene amplicon survey at eight sites and multiple time points along the Louisiana coast and one inland swamp spanning freshwater to high brackish salinities, totaling 47 duplicated Sterivex (0.2-2.7 µm) and prefilter (>2.7 µm) samples. We cataloged over 13,000 Amplicon Sequence ariants (ASVs) from common freshwater and marine clades such as SAR11 (Alphaproteobacteria), Synechococcus (Cyanobacteria), and acI and Candidatus Actinomarina (Actinobacteria). We observed correlations with freshwater or marine habitats in many organisms and characterized a group of taxa with specialized distributions across brackish water sites, supporting the hypothesis of an endogenous brackish-water community. Additionally, we observed brackish-water associations for several aquatic clades typically considered marine or freshwater taxa, such as SAR11 subclade II, SAR324, and the acI Actinobacteria. The data presented here expand the geographic coverage of microbial ecology in estuarine communities, help delineate the native and transitory members of these environments, and provide critical aquatic microbiological baseline data for coastal and estuarine sites in the nGoM.IMPORTANCEEstuarine and coastal waters are diverse ecosystems influenced by tidal fluxes, interconnected wetlands, and river outflows, which are of high economic and ecological importance. Microorganisms play a pivotal role in estuaries as "first responders" and ecosystem architects, yet despite their ecological importance, they remain underrepresented in microbial studies compared to open ocean environments. This leads to substantial knowledge gaps that are important for understanding global biogeochemical cycling and making decisions about conservation and management strategies in these environments. Our study makes key contributions to the microbial ecology of estuarine and coastal habitats in the northern Gulf of Mexico. Our microbial community data support the concept of a globally distributed, core brackish microbiome and emphasize previously underrecognized brackish-water taxa. Given the projected worsening of land loss, oil spills, and natural disasters in this region, our results will serve as important baseline data for researchers investigating the microbial communities found across estuaries.

Gene inversion led to the emergence of brackish archaeal heterotrophs in the aftermath of the Cryogenian Snowball Earth

Land-ocean interactions greatly impact the evolution of coastal life on earth. However, the ancient geological forces and genetic mechanisms that shaped evolutionary adaptations and allowed microorganisms to inhabit coastal brackish waters remain largely unexplored. In this study, we infer the evolutionary trajectory of the ubiquitous heterotrophic archaea Poseidoniales (Marine Group II archaea) presently occurring across global aquatic habitats. Our results show that their brackish subgroups had a single origination, dated to over 600 million years ago, through the inversion of the magnesium transport gene corA that conferred osmotic-stress tolerance. The subsequent loss and gain of corA were followed by genome-wide adjustment, characterized by a general two-step mode of selection in microbial speciation. The coastal family of Poseidoniales showed a rapid increase in the evolutionary rate during and in the aftermath of the Cryogenian Snowball Earth (∼700 million years ago), possibly in response to the enhanced phosphorus supply and the rise of algae. Our study highlights the close interplay between genetic changes and ecosystem evolution that boosted microbial diversification in the Neoproterozoic continental margins, where the Cambrian explosion of animals soon followed.

Bacterioplankton Communities in Dissolved Organic Carbon-Rich Amazonian Black Water

The Amazon River basin sustains dramatic hydrochemical gradients defined by three water types: white, clear, and black waters. In black water, important loads of allochthonous humic dissolved organic matter (DOM) result from the bacterioplankton degradation of plant lignin. However, the bacterial taxa involved in this process remain unknown, since Amazonian bacterioplankton has been poorly studied. Its characterization could lead to a better understanding of the carbon cycle in one of the Earth's most productive hydrological systems. Our study characterized the taxonomic structure and functions of Amazonian bacterioplankton to better understand the interplay between this community and humic DOM. We conducted a field sampling campaign comprising 15 sites distributed across the three main Amazonian water types (representing a gradient of humic DOM), and a 16S rRNA metabarcoding analysis based on bacterioplankton DNA and RNA extracts. Bacterioplankton functions were inferred using 16S rRNA data in combination with a tailored functional database from 90 Amazonian basin shotgun metagenomes from the literature. We discovered that the relative abundances of fluorescent DOM fractions (humic-, fulvic-, and protein-like) were major drivers of bacterioplankton structure. We identified 36 genera for which the relative abundance was significantly correlated with humic DOM. The strongest correlations were found in the Polynucleobacter, Methylobacterium, and Acinetobacter genera, three low abundant but omnipresent taxa that possessed several genes involved in the main steps of the β-aryl ether enzymatic degradation pathway of diaryl humic DOM residues. Overall, this study identified key taxa with DOM degradation genomic potential, the involvement of which in allochthonous Amazonian carbon transformation and sequestration merits further investigation. IMPORTANCE The Amazon basin discharge carries an important load of terrestrially derived dissolved organic matter (DOM) to the ocean. The bacterioplankton from this basin potentially plays important roles in transforming this allochthonous carbon, which has consequences on marine primary productivity and global carbon sequestration. However, the structure and function of Amazonian bacterioplanktonic communities remain poorly studied, and their interactions with DOM are unresolved. In this study, we (i) sampled bacterioplankton in all the main Amazon tributaries, (ii) combined information from the taxonomic structure and functional repertory of Amazonian bacterioplankton communities to understand their dynamics, (iii) identified the main physicochemical parameters shaping bacterioplanktonic communities among a set of >30 measured environmental parameters, and (iv) characterized how bacterioplankton structure varies according to the relative abundance of humic compounds, a by-product from the bacterial degradation process of allochthonous DOM.

Evaluation of computational phage detection tools for metagenomic datasets

Introduction

As new computational tools for detecting phage in metagenomes are being rapidly developed, a critical need has emerged to develop systematic benchmarks.

Methods

In this study, we surveyed 19 metagenomic phage detection tools, 9 of which could be installed and run at scale. Those 9 tools were assessed on several benchmark challenges. Fragmented reference genomes are used to assess the effects of fragment length, low viral content, phage taxonomy, robustness to eukaryotic contamination, and computational resource usage. Simulated metagenomes are used to assess the effects of sequencing and assembly quality on the tool performances. Finally, real human gut metagenomes and viromes are used to assess the differences and similarities in the phage communities predicted by the tools.

Results

We find that the various tools yield strikingly different results. Generally, tools that use a homology approach (VirSorter, MARVEL, viralVerify, VIBRANT, and VirSorter2) demonstrate low false positive rates and robustness to eukaryotic contamination. Conversely, tools that use a sequence composition approach (VirFinder, DeepVirFinder, Seeker), and MetaPhinder, have higher sensitivity, including to phages with less representation in reference databases. These differences led to widely differing predicted phage communities in human gut metagenomes, with nearly 80% of contigs being marked as phage by at least one tool and a maximum overlap of 38.8% between any two tools. While the results were more consistent among the tools on viromes, the differences in results were still significant, with a maximum overlap of 60.65%. Discussion: Importantly, the benchmark datasets developed in this study are publicly available and reusable to enable the future comparability of new tools developed.

Genomic and Environmental Factors Shape the Active Gill Bacterial Community of an Amazonian Teleost Holobiont

Fish bacterial communities provide functions critical for their host's survival in contrasting environments. These communities are sensitive to environmental-specific factors (i.e., physicochemical parameters, bacterioplankton), and host-specific factors (i.e., host genetic background). The relative contribution of these factors shaping Amazonian fish bacterial communities is largely unknown. Here, we investigated this topic by analyzing the gill bacterial communities of 240 wild flag cichlids (Mesonauta festivus) from 4 different populations (genetic clusters) distributed across 12 sites in 2 contrasting water types (ion-poor/acidic black water and ion-rich/circumneutral white water). Transcriptionally active gill bacterial communities were characterized by a 16S rRNA metabarcoding approach carried on RNA extractions. They were analyzed using comprehensive data sets from the hosts genetic background (Genotyping-By-Sequencing), the bacterioplankton (16S rRNA) and a set of 34 environmental parameters. Results show that the taxonomic structure of 16S rRNA gene transcripts libraries were significantly different between the 4 genetic clusters and also between the 2 water types. However, results suggest that the contribution of the host's genetic background was relatively weak in comparison to the environment-related factors in structuring the relative abundance of different active gill bacteria species. This finding was also confirmed by a mixed-effects modeling analysis, which indicated that the dissimilarity between the taxonomic structure of bacterioplanktonic communities possessed the best explicative power regarding the dissimilarity between gill bacterial communities' structure, while pairwise fixation indexes (FST) from the hosts' genetic data only had a weak explicative power. We discuss these results in terms of bacterial community assembly processes and flag cichlid fish ecology. IMPORTANCE Host-associated microbial communities respond to factors specific to the host physiology, genetic backgrounds, and life history. However, these communities also show different degrees of sensitivity to environment-dependent factors, such as abiotic physico-chemical parameters and ecological interactions. The relative importance of host- versus environment-associated factors in shaping teleost bacterial communities is still understudied and is paramount for their conservation and aquaculture. Here, we studied the relative importance of host- and environment-associated factors structuring teleost bacterial communities using gill samples from a wild Amazonian teleost model (Mesonauta festivus) sampled in contrasting habitats along a 1500 km section of the Amazonian basin, thus ensuring high genetic diversity. Results showed that the contribution of the host's genetic background was weak compared to environment-related bacterioplanktonic communities in shaping gill bacterial assemblages, thereby suggesting that our understanding of teleost microbiome assembly could benefit from further studies focused on the ecological interplay between host-associated and free-living communities.

An improved workflow for accurate and robust healthcare environmental surveillance using metagenomics

BACKGROUND

Effective surveillance of microbial communities in the healthcare environment is increasingly important in infection prevention. Metagenomics-based techniques are promising due to their untargeted nature but are currently challenged by several limitations: (1) they are not powerful enough to extract valid signals out of the background noise for low-biomass samples, (2) they do not distinguish between viable and nonviable organisms, and (3) they do not reveal the microbial load quantitatively. An additional practical challenge towards a robust pipeline is the inability to efficiently allocate sequencing resources a priori. Assessment of sequencing depth is generally practiced post hoc, if at all, for most microbiome studies, regardless of the sample type. This practice is inefficient at best, and at worst, poor sequencing depth jeopardizes the interpretation of study results. To address these challenges, we present a workflow for metagenomics-based environmental surveillance that is appropriate for low-biomass samples, distinguishes viability, is quantitative, and estimates sequencing resources.

RESULTS

The workflow was developed using a representative microbiome sample, which was created by aggregating 120 surface swabs collected from a medical intensive care unit. Upon evaluating and optimizing techniques as well as developing new modules, we recommend best practices and introduce a well-structured workflow. We recommend adopting liquid-liquid extraction to improve DNA yield and only incorporating whole-cell filtration when the nonbacterial proportion is large. We suggest including propidium monoazide treatment coupled with internal standards and absolute abundance profiling for viability assessment and involving cultivation when demanding comprehensive profiling. We further recommend integrating internal standards for quantification and additionally qPCR when we expect poor taxonomic classification. We also introduce a machine learning-based model to predict required sequencing effort from accessible sample features. The model helps make full use of sequencing resources and achieve desired outcomes. Video Abstract CONCLUSIONS: This workflow will contribute to more accurate and robust environmental surveillance and infection prevention. Lessons gained from this study will also benefit the continuing development of methods in relevant fields.

Two Metatranscriptomic Profiles through Low-Dissolved-Oxygen Waters (DO, 0 to 33 µM) in the Eastern Tropical North Pacific Ocean

We present 16 seawater metatranscriptomes collected from a marine oxygen-deficient zone (ODZ) in the eastern tropical North Pacific (ETNP). This data set will be useful for identifying shifts in microbial community structure and function through oxic/anoxic transition zones, where overlapping aerobic and anaerobic microbial processes impact marine biogeochemical cycling.

Microbial population genomes from the Amazon River reveal possible modulation of the organic matter degradation process in tropical freshwaters

Rivers connect the carbon cycle in land with that in aquatic ecosystems by transporting and transforming terrestrial organic matter (TeOM). The Amazon River receives huge loads of TeOM from the surrounding rainforest, promoting a substantial microbial heterotrophic activity and consequently, CO₂ outgassing. In the Amazon River, microbes degrade up to 55% of the lignin present in the TeOM. Yet, the main microbial genomes involved in TeOM degradation were unknown. Here, we characterize 51 population genomes (PGs) representing some of the most abundant microbes in the Amazon River deriving from 106 metagenomes. The 51 reconstructed PGs are among the most abundant microbes in the Amazon River, and 53% of them are not able to degrade TeOM. Among the PGs capable of degrading TeOM, 20% were exclusively cellulolytic, while the others could also oxidize lignin. The transport and consumption of lignin oxidation byproducts seemed to be decoupled from the oxidation process, being apparently performed by different groups of microorganisms. By connecting the genomic features of abundant microbes in the Amazon River with the degradation machinery of TeOM, we suggest that a complex microbial consortium could explain the quick turnover of TeOM previously observed in this ecosystem.

Parasites, pathogens, and other symbionts of copepods

There is a large diversity of eukaryotic symbionts of copepods, dominated by epizootic protists such as ciliates, and metazoan parasites. Eukaryotic endoparasites, copepod-associated bacteria, and viruses are less well known, partly due to technical limitations. However, new molecular techniques, combined with a range of other approaches, provide a complementary toolkit for understanding the complete symbiome of copepods and how the symbiome relates to their ecological roles, relationships with other biota, and responses to environmental change. In this review we provide the most complete overview of the copepod symbiome to date, including microeukaryotes, metazoan parasites, bacteria, and viruses, and provide extensive literature databases to inform future studies.

Sweet spheres: succession and CAZyme expression of marine bacterial communities colonizing a mix of alginate and pectin particles

Polysaccharide particles are important substrates and microhabitats for marine bacteria. However, substrate-specific bacterial dynamics in mixtures of particle types with different polysaccharide composition, as likely occurring in natural habitats, are undescribed. Here, we studied the composition, functional diversity and gene expression of marine bacterial communities colonizing a mix of alginate and pectin particles. Amplicon, metagenome and metatranscriptome sequencing revealed that communities on alginate and pectin particles significantly differed from their free-living counterparts. Unexpectedly, microbial dynamics on alginate and pectin particles were similar, with predominance of amplicon sequence variants (ASVs) from Tenacibaculum, Colwellia, Psychrobium and Psychromonas. Corresponding metagenome-assembled genomes (MAGs) expressed diverse alginate lyases, several colocalized in polysaccharide utilization loci. Only a single, low-abundant MAG showed elevated transcript abundances of pectin-degrading enzymes. One specific Glaciecola ASV dominated the free-living fraction, possibly persisting on particle-derived oligomers through different glycoside hydrolases. Elevated ammonium uptake and metabolism signified nitrogen as an important factor for degrading carbon-rich particles, whereas elevated methylcitrate and glyoxylate cycles suggested nutrient limitation in surrounding waters. The bacterial preference for alginate, whereas pectin primarily served as colonization scaffold, illuminates substrate-driven dynamics within mixed polysaccharide pools. These insights expand our understanding of bacterial niche specialization and the biological carbon pump in macroalgae-rich habitats.

Abundance and Diversity of Phages, Microbial Taxa, and Antibiotic Resistance Genes in the Sediments of the River Ganges Through Metagenomic Approach

In this study, we have analyzed the metagenomic DNA from the pooled sediment sample of the river Ganges to explore the abundance and diversity of phages, microbial community, and antibiotic resistance genes (ARGs). Utilizing data from Illumina platform, 4,174 (∼0.0013%) reads were classified for the 285 different DNA viruses largely dominated by the group of 260 distinctive phages (3,602 reads, ∼86.3%). Among all, Microcystis (782 hits), Haemophilus (403), Synechococcus (386), Pseudomonas (279), Enterococcus (232), Bacillus (196), Rhodococcus (166), Caulobacter (163), Salmonella (146), Enterobacteria (143), Mycobacterium and (128) phages show the highest abundance and account for ∼90% of the total identified phages. In addition, we have also identified corresponding host pertaining to these phages. Mainly, Proteobacteria (∼69.3%) dominates the microbial population structure. Primarily, orders such as Caulobacterales (∼28%), Burkholderiales (∼13.9%), Actinomycetales (∼13.7%), and Pseudomonadales (∼7.5%) signify the core section. Furthermore, 21,869 (∼0.00695%) reads were classified in 20 ARG types (classes) and 240 ARGs (subtypes), among which 4 ARG types, namely multidrug resistance (12,041 reads, ∼55%), bacitracin (3,202 reads, ∼15%), macrolide-lincosamide-streptogramin (1,744 reads, ∼7.98%), and fosmidomycin (990 reads, ∼4.53%), have the highest abundance. Simultaneously, six resistance mechanisms were also recognized with the dominance of antibiotic efflux (72.8%, 15,919 reads). The results unveil the distribution of (pro)-phages; microbial community; and various ARGs in the Ganges river sediments.

Toward Accurate and Robust Environmental Surveillance Using Metagenomics

Environmental surveillance is a critical tool for combatting public health threats represented by the global COVID-19 pandemic and the continuous increase of antibiotic resistance in pathogens. With its power to detect entire microbial communities, metagenomics-based methods stand out in addressing the need. However, several hurdles remain to be overcome in order to generate actionable interpretations from metagenomic sequencing data for infection prevention. Conceptually and technically, we focus on viability assessment, taxonomic resolution, and quantitative metagenomics, and discuss their current advancements, necessary precautions and directions to further development. We highlight the importance of building solid conceptual frameworks and identifying rational limits to facilitate the application of techniques. We also propose the usage of internal standards as a promising approach to overcome analytical bottlenecks introduced by low biomass samples and the inherent lack of quantitation in metagenomics. Taken together, we hope this perspective will contribute to bringing accurate and consistent metagenomics-based environmental surveillance to the ground.

Planet Microbe: a platform for marine microbiology to discover and analyze interconnected 'omics and environmental data

In recent years, large-scale oceanic sequencing efforts have provided a deeper understanding of marine microbial communities and their dynamics. These research endeavors require the acquisition of complex and varied datasets through large, interdisciplinary and collaborative efforts. However, no unifying framework currently exists for the marine science community to integrate sequencing data with physical, geological, and geochemical datasets. Planet Microbe is a web-based platform that enables data discovery from curated historical and on-going oceanographic sequencing efforts. In Planet Microbe, each ‘omics sample is linked with other biological and physiochemical measurements collected for the same water samples or during the same sample collection event, to provide a broader environmental context. This work highlights the need for curated aggregation efforts that can enable new insights into high-quality metagenomic datasets. Planet Microbe is freely accessible from https://www.planetmicrobe.org/.

Uncovering the genomic potential of the Amazon River microbiome to degrade rainforest organic matter

BACKGROUND

The Amazon River is one of the largest in the world and receives huge amounts of terrestrial organic matter (TeOM) from the surrounding rainforest. Despite this TeOM is typically recalcitrant (i.e. resistant to degradation), only a small fraction of it reaches the ocean, pointing to a substantial TeOM degradation by the river microbiome. Yet, microbial genes involved in TeOM degradation in the Amazon River were barely known. Here, we examined the Amazon River microbiome by analysing 106 metagenomes from 30 sampling points distributed along the river.

RESULTS

We constructed the Amazon River basin Microbial non-redundant Gene Catalogue (AMnrGC) that includes ~ 3.7 million non-redundant genes, affiliating mostly to bacteria. We found that the Amazon River microbiome contains a substantial gene-novelty compared to other relevant known environments (rivers and rainforest soil). Genes encoding for proteins potentially involved in lignin degradation pathways were correlated to tripartite tricarboxylates transporters and hemicellulose degradation machinery, pointing to a possible priming effect. Based on this, we propose a model on how the degradation of recalcitrant TeOM could be modulated by labile compounds in the Amazon River waters. Our results also suggest changes of the microbial community and its genomic potential along the river course.

CONCLUSIONS

Our work contributes to expand significantly our comprehension of the world's largest river microbiome and its potential metabolism related to TeOM degradation. Furthermore, the produced gene catalogue (AMnrGC) represents an important resource for future research in tropical rivers. Video abstract.

Metabolic flexibility of SUP05 under low DO growth conditions

Chemoautotrophic bacteria from the SUP05 clade often dominate anoxic waters within marine oxygen minimum zones (OMZs) where they use energy gained from the oxidation of reduced sulfur to fuel carbon fixation. Some of these SUP05 bacteria are facultative aerobes that can use either nitrate or oxygen as a terminal electron acceptor making them ideally suited to thrive at the boundaries of OMZs where they experience fluctuations in dissolved oxygen (DO). SUP05 metabolism in these regions, and therefore the biogeochemical function of SUP05, depends largely on their sensitivity to oxygen. We evaluated growth and quantified differences in gene expression in Ca. T. autotrophicus strain EF1 from the SUP05 clade under high DO (22 μM), anoxic, and low DO (3.8 μM) concentrations. We show that strain EF1 cells respire oxygen and nitrate and that cells have higher growth rates, express more genes, and fix more carbon when oxygen becomes available for aerobic respiration. Evidence that facultatively aerobic SUP05 are more active and respire nitrate when oxygen becomes available at low concentrations suggests that they are an important source of nitrite across marine OMZ boundary layers.

BMC3C: binning metagenomic contigs using codon usage, sequence composition and read coverage

Motivation

Metagenomics investigates the DNA sequences directly recovered from environmental samples. It often starts with reads assembly, which leads to contigs rather than more complete genomes. Therefore, contig binning methods are subsequently used to bin contigs into genome bins. While some clustering-based binning methods have been developed, they generally suffer from problems related to stability and robustness.

Results

We introduce BMC3C, an ensemble clustering-based method, to accurately and robustly bin contigs by making use of DNA sequence Composition, Coverage across multiple samples and Codon usage. BMC3C begins by searching the proper number of clusters and repeatedly applying the k-means clustering with different initializations to cluster contigs. Next, a weight graph with each node representing a contig is derived from these clusters. If two contigs are frequently grouped into the same cluster, the weight between them is high, and otherwise low. BMC3C finally employs a graph partitioning technique to partition the weight graph into subgraphs, each corresponding to a genome bin. We conduct experiments on both simulated and real-world datasets to evaluate BMC3C, and compare it with the state-of-the-art binning tools. We show that BMC3C has an improved performance compared to these tools. To our knowledge, this is the first time that the codon usage features and ensemble clustering are used in metagenomic contig binning.

Availability and implementation

The codes of BMC3C are available at http://mlda.swu.edu.cn/codes.php?name=BMC3C.

Supplementary information

Supplementary data are available at Bioinformatics online.

Identification of discriminatory antibiotic resistance genes among environmental resistomes using extremely randomized tree algorithm

BACKGROUND

The interconnectivities of built and natural environments can serve as conduits for the proliferation and dissemination of antibiotic resistance genes (ARGs). Several studies have compared the broad spectrum of ARGs (i.e., "resistomes") in various environmental compartments, but there is a need to identify unique ARG occurrence patterns (i.e., "discriminatory ARGs"), characteristic of each environment. Such an approach will help to identify factors influencing ARG proliferation, facilitate development of relative comparisons of the ARGs distinguishing various environments, and help pave the way towards ranking environments based on their likelihood of contributing to the spread of clinically relevant antibiotic resistance. Here we formulate and demonstrate an approach using an extremely randomized tree (ERT) algorithm combined with a Bayesian optimization technique to capture ARG variability in environmental samples and identify the discriminatory ARGs. The potential of ERT for identifying discriminatory ARGs was first evaluated using in silico metagenomic datasets (simulated metagenomic Illumina sequencing data) with known variability. The application of ERT was then demonstrated through analyses using publicly available and in-house metagenomic datasets associated with (1) different aquatic habitats (e.g., river, wastewater influent, hospital effluent, and dairy farm effluent) to compare resistomes between distinct environments and (2) different river samples (i.e., Amazon, Kalamas, and Cam Rivers) to compare resistome characteristics of similar environments.

RESULTS

The approach was found to readily identify discriminatory ARGs in the in silico datasets. Also, it was not found to be biased towards ARGs with high relative abundance, which is a common limitation of feature projection methods, and instead only captured those ARGs that elicited significant profiles. Analyses of publicly available metagenomic datasets further demonstrated that the ERT approach can effectively differentiate real-world environmental samples and identify discriminatory ARGs based on pre-defined categorizing schemes.

CONCLUSIONS

Here a new methodology was formulated to characterize and compare variances in ARG profiles between metagenomic data sets derived from similar/dissimilar environments. Specifically, identification of discriminatory ARGs among samples representing various environments can be identified based on factors of interest. The methodology could prove to be a particularly useful tool for ARG surveillance and the assessment of the effectiveness of strategies for mitigating the spread of antibiotic resistance. The python package is hosted in the Git repository: https://github.com/gaarangoa/ExtrARG.

Marine-freshwater prokaryotic transitions require extensive changes in the predicted proteome

BACKGROUND

The adaptation of a marine prokaryote to live in freshwater environments or vice versa is generally believed to be an unusual and evolutionary demanding process. However, the reasons are not obvious given the similarity of both kinds of habitats.

RESULTS

We have found major differences at the level of the predicted metaproteomes of marine and freshwater habitats with more acidic values of the isoelectric points (pI) in marine microbes. Furthermore, by comparing genomes of marine-freshwater phylogenetic relatives, we have found higher pI values (basic shift) in the freshwater ones. This difference was sharper in secreted > cytoplasmic > membrane proteins. The changes are concentrated on the surface of soluble proteins. It is also detectable at the level of total amino acid composition and involves similarly core and flexible genome- encoded proteins.

CONCLUSIONS

The marked changes at the level of protein amino acid composition and pI provide a tool to predict the preferred habitat of a culture or a metagenome-assembled genome (MAG). The exact physiological explanation for such variations in the pIs and electrostatic surface potentials is not known yet. However, these changes might reflect differences in membrane bioenergetics derived from the absence of significant Na⁺ concentrations in most freshwater habitats. In any case, the changes in amino acid composition in most proteins imply that a long evolutionary time is required to adapt from one type of habitat to the other.

The MAR databases: development and implementation of databases specific for marine metagenomics

We introduce the marine databases; MarRef, MarDB and MarCat (https://mmp.sfb.uit.no/databases/), which are publicly available resources that promote marine research and innovation. These data resources, which have been implemented in the Marine Metagenomics Portal (MMP) (https://mmp.sfb.uit.no/), are collections of richly annotated and manually curated contextual (metadata) and sequence databases representing three tiers of accuracy. While MarRef is a database for completely sequenced marine prokaryotic genomes, which represent a marine prokaryote reference genome database, MarDB includes all incomplete sequenced prokaryotic genomes regardless level of completeness. The last database, MarCat, represents a gene (protein) catalog of uncultivable (and cultivable) marine genes and proteins derived from marine metagenomics samples. The first versions of MarRef and MarDB contain 612 and 3726 records, respectively. Each record is built up of 106 metadata fields including attributes for sampling, sequencing, assembly and annotation in addition to the organism and taxonomic information. Currently, MarCat contains 1227 records with 55 metadata fields. Ontologies and controlled vocabularies are used in the contextual databases to enhance consistency. The user-friendly web interface lets the visitors browse, filter and search in the contextual databases and perform BLAST searches against the corresponding sequence databases. All contextual and sequence databases are freely accessible and downloadable from https://s1.sfb.uit.no/public/mar/.

Metatranscriptomic exploration of microbial functioning in clouds

Clouds constitute the uppermost layer of the biosphere. They host diverse communities whose functioning remains obscure, although biological activity potentially participates to atmospheric chemical and physical processes. In order to gain information on the metabolic functioning of microbial communities in clouds, we conducted coordinated metagenomics/metatranscriptomics profiling of cloud water microbial communities. Samples were collected from a high altitude atmospheric station in France and examined for biological content after untargeted amplification of nucleic acids. Living microorganisms, essentially bacteria, maintained transcriptional and translational activities and expressed many known complementary physiological responses intended to fight oxidants, osmotic variations and cold. These included activities of oxidant detoxification and regulation, synthesis of osmoprotectants/cryoprotectants, modifications of membranes, iron uptake. Consistently these energy-demanding processes were fueled by central metabolic routes involved in oxidative stress response and redox homeostasis management, such as pentose phosphate and glyoxylate pathways. Elevated binding and transmembrane ion transports demonstrated important interactions between cells and their cloud droplet chemical environments. In addition, polysaccharides, potentially beneficial for survival like exopolysaccharides, biosurfactants and adhesins, were synthesized. Our results support a biological influence on cloud physical and chemical processes, acting notably on the oxidant capacity, iron speciation and availability, amino-acids distribution and carbon and nitrogen fates.

Impacts of land-use and land-cover change on stream hydrochemistry in the Cerrado and Amazon biomes

Studies on the impacts of land-use and land-cover change on stream hydrochemistry in active deforestation zones of the Amazon agricultural frontier are limited and have often used low-temporal-resolution datasets. Moreover, these impacts are not concurrently assessed in well-established agricultural areas and new deforestations hotspots. We aimed to identify these impacts using an experimental setup to collect high-temporal-resolution hydrological and hydrochemical data in two pairs of low-order streams in catchments under contrasting land use and land cover (native vegetation vs. pasture) in the Amazon and Cerrado biomes. Our results indicate that the conversion of natural landscapes to pastures increases carbon and nutrient fluxes via streamflow in both biomes. These changes were the greatest in total inorganic carbon in the Amazon and in potassium in the Cerrado, representing a 5.0- and 5.5-fold increase in the fluxes of each biome, respectively. We found that stormflow, which is often neglected in studies on stream hydrochemistry in the tropics, plays a substantial role in the carbon and nutrient fluxes, especially in the Amazon biome, as its contributions to hydrochemical fluxes are mostly greater than the volumetric contribution to the total streamflow. These findings demonstrate that assessments of the impacts of deforestation in the Amazon and Cerrado biomes should also take into account rapid hydrological pathways; however, this can only be achieved through collection of high-temporal-resolution data.

Virioplankton Assemblage Structure in the Lower River and Ocean Continuum of the Amazon

The Amazon River forms a vast plume in the Atlantic Ocean that can extend for more than 1,000 km. Microbial communities promote a globally relevant carbon sink system in the plume. Despite the importance of viruses for the global carbon cycle, the diversity and the possible roles of viruses in the Amazonia are poorly understood. The present work assesses, for the first time, the abundance and diversity of viruses simultaneously in the river and ocean in order to elucidate their possible roles. DNA sequence assembly yielded 29,358 scaffolds, encoding 82,546 viral proteins, with 15 new complete viral genomes from the 12 river and ocean locations. Viral diversity was clearly distinguished by river and ocean. Bacteriophages were the most abundant and occurred throughout the continuum. Viruses that infect eukaryotes were more abundant in the river, whereas phages appeared to have strong control over the host prokaryotic populations in the plume.

The Amazon River watershed and its associated plume comprise a vast continental and oceanic area. The microbial activities along this continuum contribute substantially to global carbon and nutrient cycling, and yet there is a dearth of information on the diversity, abundance, and possible roles of viruses in this globally important river. The aim of this study was to elucidate the diversity and structure of virus assemblages of the Amazon River-ocean continuum. Environmental viral DNA sequences were obtained for 12 locations along the river’s lower reach (n = 5) and plume (n = 7). Sequence assembly yielded 29,358 scaffolds, encoding 82,546 viral proteins, with 15 new complete viral genomes. Despite the spatial connectivity mediated by the river, virome analyses and physical-chemical water parameters clearly distinguished river and plume ecosystems. Bacteriophages were ubiquitous in the continuum and were more abundant in the transition region. Eukaryotic viruses occurred mostly in the river, while the plume had more viruses of autotrophic organisms (Prochlorococcus, Synechococcus) and heterotrophic bacteria (Pelagibacter). The viral families Microviridae and Myoviridae were the most abundant and occurred throughout the continuum. The major functions of the genes in the continuum involved viral structures and life cycles, and viruses from plume locations and Tapajós River showed the highest levels of functional diversity. The distribution patterns of the viral assemblages were defined not only by the occurrence of possible hosts but also by water physical and chemical parameters, especially salinity. The findings presented here help to improve understanding of the possible roles of viruses in the organic matter cycle along the river-ocean continuum.

IMPORTANCE The Amazon River forms a vast plume in the Atlantic Ocean that can extend for more than 1,000 km. Microbial communities promote a globally relevant carbon sink system in the plume. Despite the importance of viruses for the global carbon cycle, the diversity and the possible roles of viruses in the Amazon are poorly understood. The present work assesses, for the first time, the abundance and diversity of viruses simultaneously in the river and ocean in order to elucidate their possible roles. DNA sequence assembly yielded 29,358 scaffolds, encoding 82,546 viral proteins, with 15 new complete viral genomes from the 12 river and ocean locations. Viral diversity was clearly distinguished by river and ocean. Bacteriophages were the most abundant and occurred throughout the continuum. Viruses that infect eukaryotes were more abundant in the river, whereas phages appeared to have strong control over the host prokaryotic populations in the plume.

1,003 reference genomes of bacterial and archaeal isolates expand coverage of the tree of life

We present 1,003 reference genomes that were sequenced as part of the Genomic Encyclopedia of Bacteria and Archaea (GEBA) initiative, selected to maximize sequence coverage of phylogenetic space. These genomes double the number of existing type strains and expand their overall phylogenetic diversity by 25%. Comparative analyses with previously available finished and draft genomes reveal a 10.5% increase in novel protein families as a function of phylogenetic diversity. The GEBA genomes recruit 25 million previously unassigned metagenomic proteins from 4,650 samples, improving their phylogenetic and functional interpretation. We identify numerous biosynthetic clusters and experimentally validate a divergent phenazine cluster with potential new chemical structure and antimicrobial activity. This Resource is the largest single release of reference genomes to date. Bacterial and archaeal isolate sequence space is still far from saturated, and future endeavors in this direction will continue to be a valuable resource for scientific discovery.

Bacterial Biogeography across the Amazon River-Ocean Continuum

Spatial and temporal patterns in microbial biodiversity across the Amazon river-ocean continuum were investigated along ∼675 km of the lower Amazon River mainstem, in the Tapajós River tributary, and in the plume and coastal ocean during low and high river discharge using amplicon sequencing of 16S rRNA genes in whole water and size-fractionated samples (0.2–2.0 μm and >2.0 μm). River communities varied among tributaries, but mainstem communities were spatially homogeneous and tracked seasonal changes in river discharge and co-varying factors. Co-occurrence network analysis identified strongly interconnected river assemblages during high (May) and low (December) discharge periods, and weakly interconnected transitional assemblages in September, suggesting that this system supports two seasonal microbial communities linked to river discharge. In contrast, plume communities showed little seasonal differences and instead varied spatially tracking salinity. However, salinity explained only a small fraction of community variability, and plume communities in blooms of diatom-diazotroph assemblages were strikingly different than those in other high salinity plume samples. This suggests that while salinity physically structures plumes through buoyancy and mixing, the composition of plume-specific communities is controlled by other factors including nutrients, phytoplankton community composition, and dissolved organic matter chemistry. Co-occurrence networks identified interconnected assemblages associated with the highly productive low salinity near-shore region, diatom-diazotroph blooms, and the plume edge region, and weakly interconnected assemblages in high salinity regions. This suggests that the plume supports a transitional community influenced by immigration of ocean bacteria from the plume edge, and by species sorting as these communities adapt to local environmental conditions. Few studies have explored patterns of microbial diversity in tropical rivers and coastal oceans. Comparison of Amazon continuum microbial communities to those from temperate and arctic systems suggest that river discharge and salinity are master variables structuring a range of environmental conditions that control bacterial communities across the river-ocean continuum.

Expression patterns of elemental cycling genes in the Amazon River Plume

Metatranscriptomics and metagenomics data sets benchmarked with internal standards were used to characterize the expression patterns for biogeochemically relevant bacterial and archaeal genes mediating carbon, nitrogen, phosphorus and sulfur uptake and metabolism through the salinity gradient of the Amazon River Plume. The genes were identified in 48 metatranscriptomic and metagenomic data sets summing to >500 million quality-controlled reads from six locations in the plume ecosystem. The ratio of transcripts per gene copy (a direct measure of expression made possible by internal standard additions) showed that the free-living bacteria and archaea exhibited only small changes in the expression levels of biogeochemically relevant genes through the salinity and nutrient zones of the plume. In contrast, the expression levels of genes in particle-associated cells varied over orders of magnitude among the stations, with the largest differences measured for genes mediating aspects of nitrogen cycling (nifH, amtB and amoA) and phosphorus acquisition (pstC, phoX and phoU). Taxa varied in their baseline gene expression levels and extent of regulation, and most of the spatial variation in the expression level could be attributed to changes in gene regulation after removing the effect of shifting taxonomic composition. We hypothesize that changes in microbial element cycling along the Amazon River Plume are largely driven by shifting activities of particle-associated cells, with most activities peaking in the mesohaline regions where N₂ fixation rates are elevated.

Patterns of Transcript Abundance of Eukaryotic Biogeochemically-Relevant Genes in the Amazon River Plume

The Amazon River has the largest discharge of all rivers on Earth, and its complex plume system fuels a wide array of biogeochemical processes, across a large area of the western tropical North Atlantic. The plume thus stimulates microbial processes affecting carbon sequestration and nutrient cycles at a global scale. Chromosomal gene expression patterns of the 2.0 to 156 μm size-fraction eukaryotic microbial community were investigated in the Amazon River Plume, generating a robust dataset (more than 100 million mRNA sequences) that depicts the metabolic capabilities and interactions among the eukaryotic microbes. Combining classical oceanographic field measurements with metatranscriptomics yielded characterization of the hydrographic conditions simultaneous with a quantification of transcriptional activity and identity of the community. We highlight the patterns of eukaryotic gene expression for 31 biogeochemically significant gene targets hypothesized to be valuable within forecasting models. An advantage to this targeted approach is that the database of reference sequences used to identify the target genes was selectively constructed and highly curated optimizing taxonomic coverage, throughput, and the accuracy of annotations. A coastal diatom bloom highly expressed nitrate transporters and carbonic anhydrase presumably to support high growth rates and enhance uptake of low levels of dissolved nitrate and CO₂. Diatom-diazotroph association (DDA: diatoms with nitrogen fixing symbionts) blooms were common when surface salinity was mesohaline and dissolved nitrate concentrations were below detection, and hence did not show evidence of nitrate utilization, suggesting they relied on ammonium transporters to aquire recently fixed nitrogen. These DDA blooms in the outer plume had rapid turnover of the photosystem D1 protein presumably caused by photodegradation under increased light penetration in clearer waters, and increased expression of silicon transporters as silicon became limiting. Expression of these genes, including carbonic anhydrase and transporters for nitrate and phosphate, were found to reflect the physiological status and biogeochemistry of river plume environments. These relatively stable patterns of eukaryotic transcript abundance occurred over modest spatiotemporal scales, with similarity observed in sample duplicates collected up to 2.45 km in space and 120 minutes in time. These results confirm the use of metatranscriptomics as a valuable tool to understand and predict microbial community function.

An extensive reef system at the Amazon River mouth

Large rivers create major gaps in reef distribution along tropical shelves. The Amazon River represents 20% of the global riverine discharge to the ocean, generating up to a 1.3 × 10(6)-km(2) plume, and extensive muddy bottoms in the equatorial margin of South America. As a result, a wide area of the tropical North Atlantic is heavily affected in terms of salinity, pH, light penetration, and sedimentation. Such unfavorable conditions were thought to imprint a major gap in Western Atlantic reefs. We present an extensive carbonate system off the Amazon mouth, underneath the river plume. Significant carbonate sedimentation occurred during lowstand sea level, and still occurs in the outer shelf, resulting in complex hard-bottom topography. A permanent near-bottom wedge of ocean water, together with the seasonal nature of the plume's eastward retroflection, conditions the existence of this extensive (~9500 km(2)) hard-bottom mosaic. The Amazon reefs transition from accretive to erosional structures and encompass extensive rhodolith beds. Carbonate structures function as a connectivity corridor for wide depth-ranging reef-associated species, being heavily colonized by large sponges and other structure-forming filter feeders that dwell under low light and high levels of particulates. The oxycline between the plume and subplume is associated with chemoautotrophic and anaerobic microbial metabolisms. The system described here provides several insights about the responses of tropical reefs to suboptimal and marginal reef-building conditions, which are accelerating worldwide due to global changes.

Next-generation sequencing (NGS) for assessment of microbial water quality: current progress, challenges, and future opportunities

Water quality is an emergent property of a complex system comprised of interacting microbial populations and introduced microbial and chemical contaminants. Studies leveraging next-generation sequencing (NGS) technologies are providing new insights into the ecology of microbially mediated processes that influence fresh water quality such as algal blooms, contaminant biodegradation, and pathogen dissemination. In addition, sequencing methods targeting small subunit (SSU) rRNA hypervariable regions have allowed identification of signature microbial species that serve as bioindicators for sewage contamination in these environments. Beyond amplicon sequencing, metagenomic and metatranscriptomic analyses of microbial communities in fresh water environments reveal the genetic capabilities and interplay of waterborne microorganisms, shedding light on the mechanisms for production and biodegradation of toxins and other contaminants. This review discusses the challenges and benefits of applying NGS-based methods to water quality research and assessment. We will consider the suitability and biases inherent in the application of NGS as a screening tool for assessment of biological risks and discuss the potential and limitations for direct quantitative interpretation of NGS data. Secondly, we will examine case studies from recent literature where NGS based methods have been applied to topics in water quality assessment, including development of bioindicators for sewage pollution and microbial source tracking, characterizing the distribution of toxin and antibiotic resistance genes in water samples, and investigating mechanisms of biodegradation of harmful pollutants that threaten water quality. Finally, we provide a short review of emerging NGS platforms and their potential applications to the next generation of water quality assessment tools.

Metagenomic and metatranscriptomic inventories of the lower Amazon River, May 2011

BACKGROUND

The Amazon River runs nearly 6500 km across the South American continent before emptying into the western tropical North Atlantic Ocean. In terms of both volume and watershed area, it is the world's largest riverine system, affecting elemental cycling on a global scale.

RESULTS

A quantitative inventory of genes and transcripts benchmarked with internal standards was obtained at five stations in the lower Amazon River during May 2011. At each station, metagenomes and metatranscriptomes were obtained in duplicate for two microbial size fractions (free-living, 0.2 to 2.0 μm; particle-associated, 2.0 to 297 μm) using 150 × 150 paired-end Illumina sequencing. Forty eight sample datasets were obtained, averaging 15 × 10(6) potential protein-encoding reads each (730 × 10(6) total). Prokaryotic metagenomes and metatranscriptomes were dominated by members of the phyla Actinobacteria, Planctomycetes, Betaproteobacteria, Verrucomicrobia, Nitrospirae, and Acidobacteria. The actinobacterium SCGC AAA027-L06 reference genome recruited the greatest number of reads overall, with this single bin contributing an average of 50 billion genes and 500 million transcripts per liter of river water. Several dominant taxa were unevenly distributed between the free-living and particle-associated size fractions, such as a particle-associated bias for reads binning to planctomycete Schlesneria paludicola and a free-living bias for actinobacterium SCGC AAA027-L06. Gene expression ratios (transcripts to gene copy ratio) increased downstream from Óbidos to Macapá and Belém, indicating higher per cell activity of Amazon River bacteria and archaea as river water approached the ocean.

CONCLUSION

This inventory of riverine microbial genes and transcripts, benchmarked with internal standards for full quantitation, provides an unparalleled window into microbial taxa and functions in the globally important Amazon River ecosystem.

A new class of marine Euryarchaeota group II from the Mediterranean deep chlorophyll maximum

We have analyzed metagenomic fosmid clones from the deep chlorophyll maximum (DCM), which, by genomic parameters, correspond to the 16S ribosomal RNA (rRNA)-defined marine Euryarchaeota group IIB (MGIIB). The fosmid collections associated with this group add up to 4 Mb and correspond to at least two species within this group. From the proposed essential genes contained in the collections, we infer that large sections of the conserved regions of the genomes of these microbes have been recovered. The genomes indicate a photoheterotrophic lifestyle, similar to that of the available genome of MGIIA (assembled from an estuarine metagenome in Puget Sound, Washington Pacific coast), with a proton-pumping rhodopsin of the same kind. Several genomic features support an aerobic metabolism with diversified substrate degradation capabilities that include xenobiotics and agar. On the other hand, these MGIIB representatives are non-motile and possess similar genome size to the MGIIA-assembled genome, but with a lower GC content. The large phylogenomic gap with other known archaea indicates that this is a new class of marine Euryarchaeota for which we suggest the name Thalassoarchaea. The analysis of recruitment from available metagenomes indicates that the representatives of group IIB described here are largely found at the DCM (ca. 50 m deep), in which they are abundant (up to 0.5% of the reads), and at the surface mostly during the winter mixing, which explains formerly described 16S rRNA distribution patterns. Their uneven representation in environmental samples that are close in space and time might indicate sporadic blooms.