Nonpareil: a redundancy-based approach to assess the level of coverage in metagenomic datasets

Do coastal bacterioplankton communities hold the molecular key to the rapid biodegradation of Polycyclic Aromatic Hydrocarbons (PAHs) from shipping scrubber effluent?

Shipping scrubber effluents, containing a cocktail of Polycyclic Aromatic Hydrocarbons (PAHs), show undisputed effects at single-species experiments while PAHs fate in the marine environment after effluent discharge is still investigated. Bacterioplankton, composed of abundant diverse taxa with xenobiotic-degrading capabilities, are the first responders to scrubber emissions and can affect PAHs impacts on marine life. This work aims to examine the fate of scrubber effluent PAHs and alkyl-PAHs in mesocosms of coastal bacterioplankton communities from a pristine (phytoplankton carbon biomass was 8.16 μg C L-1) and a eutrophic (105.35 μg C L-1) coastal site. High-throughput 16S rRNA metabarcoding revealed differential responses of the bacterioplankton linked to their initial community structure and population abundances. Taxa known for their PAHs-degrading capacity were retrieved, including the genera Roseobacter, Porticoccus, Marinomonas, Arcobacter, Lentibacter, Lacinutrix, Pseudospirillum, Glaciecola, Vibrio, Marivita, and Mycobacterium, and were found to have increased roles in shifted communities by increasing their relative abundances at least 5-fold in treatments with high scrubber effluent additions. Additionally, metagenomic analysis of shotgun sequencing, indicated an increase on the number of Clusters of Orthologous Genes (COGs) associated with pathways involved in PAHs degradation. Up to 198 more COGs involved in signal transduction were retrieved in scrubber effluent enriched mesocosms compared to controls, while 15, 86, and 136 more COGs associated with naphthalene, aromatic compound, and benzoate degradation, respectively, were detected in the pristine mesocosms after effluent additions. In both experiments, bacterioplankton responses towards xenobiotic degradation under increased PAHs and alkyl-PAHs were coupled with a drop in their concentrations, below the limit of detection by Day 3 of the experiment in the eutrophic community, and by half in Day 6 in the pristine environment's community. Our findings indicate that PAHs and alkyl-PAHs impacts can be rapidly reduced in natural systems of high bacterial activity.

Long-term multi-meta-omics resolves the ecophysiological controls of seasonal N2O emissions during wastewater treatment

Nitrous oxide (N2O) is the third most important greenhouse gas and originates primarily from natural and engineered microbiomes. Effective emission mitigations are currently hindered by the largely unresolved ecophysiological controls of coexisting N2O-converting metabolisms in complex communities. To address this, we used biological wastewater treatment as a model ecosystem and combined long-term metagenome-resolved metaproteomics with ex situ kinetic and full-scale operational characterization over nearly 2 years. By leveraging the evidence independently obtained at multiple ecophysiological levels, from individual genetic potential to actual metabolism and emergent community phenotype, the cascade of environmental and operational triggers driving seasonal N2O emissions has ultimately been resolved. We identified nitrifier denitrification as the dominant N2O-producing pathway and dissolved O2 as the prime operational parameter, paving the way to the design and fostering of robust emission control strategies. This work exemplifies the untapped potential of multi-meta-omics in the mechanistic understanding and ecological engineering of microbiomes towards reducing anthropogenic impacts and advancing sustainable biotechnological developments.

Long-term climate establishes functional legacies by altering microbial traits

Long-term climate history can influence rates of soil carbon cycling but the microbial traits underlying these legacy effects are not well understood. Legacies may result if historical climate differences alter the traits of soil microbial communities, particularly those associated with carbon cycling and stress tolerance. However, it is also possible that contemporary conditions can overcome the influence of historical climate, particularly under extreme conditions. Using shotgun metagenomics, we assessed the composition of soil microbial functional genes across a mean annual precipitation gradient that previously showed evidence of strong climate legacies in soil carbon flux and extracellular enzyme activity. Sampling coincided with recovery from a regional, multi-year severe drought, allowing us to document how the strength of climate legacies varied with contemporary conditions. We found increased investment in genes associated with resource cycling with historically higher precipitation across the gradient, particularly in traits related to resource transport and complex carbon degradation. This legacy effect was strongest in seasons with the lowest soil moisture, suggesting that contemporary conditions-particularly, resource stress under water limitation-influences the strength of legacy effects. In contrast, investment in stress tolerance did not vary with historical precipitation, likely due to frequent periodic drought throughout the gradient. Differences in the relative abundance of functional genes explained over half of variation in microbial functional capacity-potential enzyme activity-more so than historical precipitation or current moisture conditions. Together, these results suggest that long-term climate can alter the functional potential of soil microbial communities, leading to legacies in carbon cycling.

Antarctic Geothermal Soils Exhibit an Absence of Regional Habitat Generalist Microorganisms

Active geothermal systems are relatively rare in Antarctica and represent metaphorical islands ideal to study microbial dispersal. In this study, we tested the macro-ecological concept that high dispersal rates result in communities being dominated by either habitat generalists or specialists by investigating the microbial communities on four geographically separated geothermal sites on three Antarctic volcanoes (Mts. Erebus, Melbourne, and Rittman). We found that the microbial communities at higher temperature (max 65°C) sites (Tramway Ridge on Erebus and Rittmann) were unique from each other and were dominated by a variety of novel Archaea from class Nitrososphaeria, while lower temperature (max 50°C) sites (Western Crater on Erebus and Melbourne) had characteristically mesophilic communities (Planctomycetes, Acidobacteriota, etc.) that were highly similar. We found that 97% of the detected microbial taxa were regional habitat specialists, with no generalists, with community assembly driven by high dispersal rates and drift (25% and 30% of community assembly, respectively), not environmental selection. Our results indicate that for microbial communities experiencing high dispersal rates between isolated communities, habitat specialists may tend to out-compete habitat generalists.

Increased antibiotic resistance gene abundance linked to intensive bacterial competition in the phyllosphere across an elevational gradient

Antibiotic resistance genes (ARGs) are ancient and widespread in natural habitats, providing survival advantages for microbiomes under challenging conditions. In mountain ecosystems, phyllosphere bacterial communities face multiple stress conditions, and the elevational gradients of mountains represent crucial environmental gradients for studying biodiversity distribution patterns. However, the distribution patterns of ARGs in the phyllosphere along elevational gradients, and their correlation with bacterial community structures, remain poorly understood. Here, we applied metagenomic analyses to investigate the abundance and diversity of ARGs in 88 phyllosphere samples collected from Mount Tianmu, a national natural reserve. Our results showed that the abundance of ARGs in the phyllosphere increased along elevational gradients and was dominated by multidrug resistance and efflux pumps. The composition of bacterial communities, rather than plant traits or abiotic factors, significantly affected ARG abundance. Moreover, increased ARG abundance was correlated with greater phylogenetic overdispersion and a greater proportion of negative associations in the bacterial co-occurrence networks, suggesting that bacterial competition primarily shapes phyllosphere resistomes. These findings constitute a major advance in the biodiversity of phyllosphere resistomes along elevations, emphasizing the significant impact of bacterial community structure and assembly on ARG distribution, and are essential for understanding the emergence of ARGs.

The coral Oculina patagonica holobiont and its response to confinement, temperature, and Vibrio infections

BACKGROUND

Extensive research on the diversity and functional roles of the microorganisms associated with reef-building corals has been promoted as a consequence of the rapid global decline of coral reefs attributed to climate change. Several studies have highlighted the importance of coral-associated algae (Symbiodinium) and bacteria and their potential roles in promoting coral host fitness and survival. However, the complex coral holobiont extends beyond these components to encompass other entities such as protists, fungi, and viruses. While each constituent has been individually investigated in corals, a comprehensive understanding of their collective roles is imperative for a holistic comprehension of coral health and resilience.

RESULTS

The metagenomic analysis of the microbiome of the coral Oculina patagonica has revealed that fungi of the genera Aspergillus, Fusarium, and Rhizofagus together with the prokaryotic genera Streptomyces, Pseudomonas, and Bacillus were abundant members of the coral holobiont. This study also assessed changes in microeukaryotic, prokaryotic, and viral communities under three stress conditions: aquaria confinement, heat stress, and Vibrio infections. In general, stress conditions led to an increase in Rhodobacteraceae, Flavobacteraceae, and Vibrionaceae families, accompanied by a decrease in Streptomycetaceae. Concurrently, there was a significant decline in both the abundance and richness of microeukaryotic species and a reduction in genes associated with antimicrobial compound production by the coral itself, as well as by Symbiodinium and fungi.

CONCLUSION

Our findings suggest that the interplay between microeukaryotic and prokaryotic components of the coral holobiont may be disrupted by stress conditions, such as confinement, increase of seawater temperature, or Vibrio infection, leading to a dysbiosis in the global microbial community that may increase coral susceptibility to diseases. Further, microeukaryotic community seems to exert influence on the prokaryotic community dynamics, possibly through predation or the production of secondary metabolites with anti-bacterial activity. Video Abstract.

The breath volatilome is shaped by the gut microbiota

The gut microbiota is widely implicated in host health and disease, inspiring translational efforts to implement our growing body of knowledge in clinical settings. However, the need to characterize gut microbiota by its genomic content limits the feasibility of rapid, point-of-care diagnostics. The microbiota produces a diverse array of xenobiotic metabolites that disseminate into tissues, including volatile organic compounds (VOCs) that may be excreted in breath. We hypothesize that breath contains gut microbe-derived VOCs that inform the composition and metabolic state of the microbiota. To explore this idea, we compared the breath volatilome and fecal gut microbiomes of 27 healthy children and found that breath VOC composition is correlated with gut microbiomes. To experimentally interrogate this finding, we devised a method for capturing exhaled breath from gnotobiotic mice. Breath volatiles are then profiled by gas-chromatography mass-spectrometry (GC-MS). Using this novel methodology, we found that the murine breath profile is markedly shaped by the composition of the gut microbiota. We also find that VOCs produced by gut microbes in pure culture can be identified in vivo in the breath of mice monocolonized with the same bacteria. Altogether, our studies identify microbe-derived VOCs excreted in breath and support a mechanism by which gut bacterial metabolism directly contributes to the mammalian breath VOC profiles.

Extremely halophilic brine community manipulation shows higher robustness of microbiomes inhabiting human-driven solar saltern than naturally driven lake

Hypersaline ecosystems display taxonomically similar assemblages with low diversities and highly dense accompanying viromes. The ecological implications of viral infection on natural microbial populations remain poorly understood, especially at finer scales of diversity. Here, we sought to investigate the influence of changes in environmental physicochemical conditions and viral predation pressure by autochthonous and allochthonous viruses on host dynamics. For this purpose, we transplanted two microbiomes coming from distant hypersaline systems (solar salterns of Es Trenc in Spain and the thalassohaline lake of Aran-Bidgol lake in Iran), by exchanging the cellular fractions with the sterile-filtered accompanying brines with and without the free extracellular virus fraction. The midterm exposure (1 month) of the microbiomes to the new conditions showed that at the supraspecific taxonomic range, the assemblies from the solar saltern brine more strongly resisted the environmental changes and viral predation than that of the lake. The metagenome-assembled genomes (MAGs) analysis revealed an intraspecific transition at the ecotype level, mainly driven by changes in viral predation pressure, by both autochthonous and allochthonous viruses.

IMPORTANCE

Viruses greatly influence succession and diversification of their hosts, yet the effects of viral infection on the ecological dynamics of natural microbial populations remain poorly understood, especially at finer scales of diversity. By manipulating the viral predation pressure by autochthonous and allochthonous viruses, we uncovered potential phage-host interaction, and their important role in structuring the prokaryote community at an ecotype level.

Selection and horizontal gene transfer underlie microdiversity-level heterogeneity in resistance gene fate during wastewater treatment

Activated sludge is the centerpiece of biological wastewater treatment, as it facilitates removal of sewage-associated pollutants, fecal bacteria, and pathogens from wastewater through semi-controlled microbial ecology. It has been hypothesized that horizontal gene transfer facilitates the spread of antibiotic resistance genes within the wastewater treatment plant, in part because of the presence of residual antibiotics in sewage. However, there has been surprisingly little evidence to suggest that sewage-associated antibiotics select for resistance at wastewater treatment plants via horizontal gene transfer or otherwise. We addressed the role of sewage-associated antibiotics in promoting antibiotic resistance using lab-scale sequencing batch reactors fed field-collected wastewater, metagenomic sequencing, and our recently developed bioinformatic tool Kairos. Here, we found confirmatory evidence that fluctuating levels of antibiotics in sewage are associated with horizontal gene transfer of antibiotic resistance genes, microbial ecology, and microdiversity-level differences in resistance gene fate in activated sludge.

Unlocking the biosynthetic potential and taxonomy of the Antarctic microbiome along temporal and spatial gradients

Extreme environments, such as Antarctica, select microbial communities that display a range of evolutionary strategies to survive and thrive under harsh environmental conditions. These include a diversity of specialized metabolites, which have the potential to be a source for new natural product discovery. Efforts using (meta)genome mining approaches to identify and understand biosynthetic gene clusters in Antarctica are still scarce, and the extent of their diversity and distribution patterns in the environment have yet to be discovered. Herein, we investigated the biosynthetic gene diversity of the biofilm microbial community of Whalers Bay, Deception Island, in the Antarctic Peninsula and revealed its distribution patterns along spatial and temporal gradients by applying metagenome mining approaches and multivariable analysis. The results showed that the Whalers Bay microbial community harbors a great diversity of biosynthetic gene clusters distributed into seven classes, with terpene being the most abundant. The phyla Proteobacteria and Bacteroidota were the most abundant in the microbial community and contributed significantly to the biosynthetic gene abundances in Whalers Bay. Furthermore, the results highlighted a significant correlation between the distribution of biosynthetic genes and taxonomic diversity, emphasizing the intricate interplay between microbial taxonomy and their potential for specialized metabolite production.IMPORTANCEThis research on antarctic microbial biosynthetic diversity in Whalers Bay, Deception Island, unveils the hidden potential of extreme environments for natural product discovery. By employing metagenomic techniques, the research highlights the extensive diversity of biosynthetic gene clusters and identifies key microbial phyla, Proteobacteria and Bacteroidota, as significant contributors. The correlation between taxonomic diversity and biosynthetic gene distribution underscores the intricate interplay governing specialized metabolite production. These findings are crucial for understanding microbial adaptation in extreme environments and hold significant implications for bioprospecting initiatives. The study opens avenues for discovering novel bioactive compounds with potential applications in medicine and industry, emphasizing the importance of preserving and exploring these polyextreme ecosystems to advance biotechnological and pharmaceutical research.

Uncovering novel bacterial and archaeal diversity: genomic insights from metagenome-assembled genomes in Cuatro Cienegas, Coahuila

A comprehensive study was conducted in the Cuatro Ciénegas Basin (CCB) in Coahuila, Mexico, which is known for its diversity of microorganisms and unique physicochemical properties. The study focused on the "Archaean Domes" (AD) site in the CCB, which is characterized by an abundance of hypersaline, non-lithifying microbial mats. In AD, we analyzed the small domes and circular structures using metagenome assembly genomes (MAGs) with the aim of expanding our understanding of the prokaryotic tree of life by uncovering previously unreported lineages, as well as analyzing the diversity of bacteria and archaea in the CCB. A total of 325 MAGs were identified, including 48 Archaea and 277 Bacteria. Remarkably, 22 archaea and 104 bacteria could not be classified even at the genus level, highlighting the remarkable novel diversity of the CCB. Besides, AD site exhibited significant diversity at the phylum level, with Proteobacteria being the most abundant, followed by Desulfobacteria, Spirochaetes, Bacteroidetes, Nanoarchaeota, Halobacteriota, Cyanobacteria, Planctomycetota, Verrucomicrobiota, Actinomycetes and Chloroflexi. In Archaea, the monophyletic groups of MAGs belonged to the Archaeoglobi, Aenigmarchaeota, Candidate Nanoarchaeota, and Halobacteriota. Among Bacteria, monophyletic groups were also identified, including Spirochaetes, Proteobacteria, Planctomycetes, Actinobacteria, Verrucomicrobia, Bacteroidetes, Candidate Bipolaricaulota, Desulfobacteria, and Cyanobacteria. These monophyletic groups were possibly influenced by geographic isolation, as well as the extreme and fluctuating environmental conditions in the pond AD, such as stoichiometric imbalance of C:N:P of 122:42:1, fluctuating pH (5-9.8) and high salinity (5.28% to saturation).

New avenues for potentially seeking microbial responses to climate change beneath Antarctic ice shelves

The signs of climate change are undeniable, and the impact of these changes on ecosystem function heavily depends on the response of microbes that underpin the food web. Antarctic ice shelf is a massive mass of floating ice that extends from the continent into the ocean, exerting a profound influence on global carbon cycles. Beneath Antarctic ice shelves, marine ice stores valuable genetic information, where marine microbial communities before the industrial revolution are archived. Here, in this proof-of-concept, by employing a combination of single-cell technologiesand metagenomics, we have been able to sequence frozen microbial DNA (≈300 years old) stored in the marine ice core B15 collected from the Filchnner-Ronne Ice Shelf. Metagenomic data indicated that Proteobacteria and Thaumarchaeota (e.g., Nitrosopumilus spp.), followed by Actinobacteria (e.g., Actinomarinales), were abundant. Remarkably, our data allow us to "travel to the past" and calibrate genomic and genetic evolutionary changes for ecologically relevant microbes and functions, such as Nitrosopumilus spp., preserved in the marine ice (≈300 years old) with those collected recently in seawater under an ice shelf (year 2017). The evolutionary divergence for the ammonia monooxygenase gene amoA involved in chemolithoautotrophy was about 0.88 amino acid and 2.8 nucleotide substitution rate per 100 sites in a century, while the accumulated rate of genomic SNPs was 2,467 per 1 Mb of genome and 100 years. Whether these evolutionary changes remained constant over the last 300 years or accelerated during post-industrial periods remains an open question that will be further elucidated.

IMPORTANCE

Several efforts have been undertaken to predict the response of microbes under climate change, mainly based on short-term microcosm experiments under forced conditions. A common concern is that manipulative experiments cannot properly simulate the response of microbes to climate change, which is a long-term evolutionary process. In this proof-of-concept study with a limited sample size, we demonstrate a novel approach yet to be fully explored in science for accessing genetic information from putative past marine microbes preserved under Antarctic ice shelves before the industrial revolution. This potentially allows us estimating evolutionary changes as exemplified in our study. We advocate for gathering a more comprehensive Antarctic marine ice core data sets across various periods and sites. Such a data set would enable the establishment of a robust baseline, facilitating a better assessment of the potential effects of climate change on key genetic signatures of microbes.

Sustained bacterial N2O reduction at acidic pH

Nitrous oxide (N2O) is a climate-active gas with emissions predicted to increase due to agricultural intensification. Microbial reduction of N2O to dinitrogen (N2) is the major consumption process but microbial N2O reduction under acidic conditions is considered negligible, albeit strongly acidic soils harbor nosZ genes encoding N2O reductase. Here, we study a co-culture derived from acidic tropical forest soil that reduces N2O at pH 4.5. The co-culture exhibits bimodal growth with a Serratia sp. fermenting pyruvate followed by hydrogenotrophic N2O reduction by a Desulfosporosinus sp. Integrated omics and physiological characterization revealed interspecies nutritional interactions, with the pyruvate fermenting Serratia sp. supplying amino acids as essential growth factors to the N2O-reducing Desulfosporosinus sp. Thus, we demonstrate growth-linked N2O reduction between pH 4.5 and 6, highlighting microbial N2O reduction potential in acidic soils.

Metagenomics of two aquifers with thermal anomalies in Mallorca Island, and proposal of new uncultivated taxa named following the rules of SeqCode

Groundwater offers an intriguing blend of distinctive physical and chemical conditions, constituting a challenge for microbial life. In Mallorca, the largest island of Balearic archipelago, harbours a variety of thermal anomalies (i.e., geothermal manifestation where surface aquifers exhibiting temperatures exceeding the regional average). The metagenomes of two aquifers in the centre and southern of the island showed Pseudomonadota to be the most represented phylum when using extracted 16S rRNA gene sequences. However, the microbial structures within and between aquifers were remarkably diverse but similar in their metabolic profiles as revealed by the metagenome-assembled genomes (MAGs) pointing to a prevalence of aerobic chemolithoautotrophic and heterotrophic metabolisms, especially in the Llucmajor aquifer. Also, some evidences of anaerobic lifestyles were detected, which would indicate that these environments either could suffer episodes of oxygen depletion or the anaerobes had been transported from deeper waters. We believe that the local environmental factors (temperature, external inputs or chemistry) seem to be more relevant than the connection and, eventually, transport of microbial cells within the aquifer in determining the highly divergent structures. Notably, most of the reconstructed genomes belonged to undescribed bacterial lineages and from them two high-quality MAGs could be classified as novel taxa named following the rules of the Code for Nomenclature of Prokaryotes Described from Sequence Data (SeqCode). Accordingly, we propose the new species and genus Costitxia debesea gen. nov., sp. nov., affiliated with the novel family Costitxiaceae fam. nov., order Costitxiales ord. nov. and class Costitxiia class. nov.; and the new new species and genus Lloretia debesea gen. nov. sp. nov. affiliated with the novel family Lloretiaceae fam. nov.

A metagenomic approach to unveil the association between fecal gut microbiota and short-chain fatty acids in diarrhea caused by diarrheagenic Escherichia coli in children

Diarrheagenic Escherichia coli (DEC) is the main cause of diarrhea in children under five years old. The virulence of DEC is tightly regulated by environmental signals influenced by the gut microbiota and its metabolites. Short-chain fatty acids (SCFAs) are the main metabolic product of anaerobic fermentation in the gut, but their role in DEC diarrhea has not yet been established. In this study, we determine the levels of acetate, propionate, and butyrate in stool samples from children with diarrhea caused by DEC, and we identify bacteria from the fecal gut microbiota associated with the production of SCFAs. The microbiota and SCFAs levels in stool samples obtained from 40 children with diarrhea and 43 healthy children were determined by 16S rRNA gene sequencing and HPLC, respectively. Additionally, shotgun metagenomics was used to identify metagenome-assembled genomes (MAGs) in a subgroup of samples. The results showed significantly higher levels of all SCFAs tested in diarrheal samples than in healthy controls. The abundance of Streptococcus sp., Limosilactobacillus, Blautia, Escherichia, Bacteroides, Megamonas, and Roseburia was higher in the DEC group than in healthy individuals. Functional analysis of bacteria and their main metabolic pathways made it possible to identify species MAGs that could be responsible for the detected SCFAs levels in DEC-positive diarrhea. In conclusion, based on our results and published data, we suggest that SCFAs may be important in the crosstalk between the microbiota and DEC pathogens in the gut.

The gut metagenome harbors metabolic and antibiotic resistance signatures of moderate-to-severe asthma

Asthma is a common allergic airway disease that has been associated with the development of the human microbiome early in life. Both the composition and function of the infant gut microbiota have been linked to asthma risk, but functional alterations in the gut microbiota of older patients with established asthma remain an important knowledge gap. Here, we performed whole metagenomic shotgun sequencing of 95 stool samples from a cross-sectional cohort of 59 healthy and 36 subjects with moderate-to-severe asthma to characterize the metagenomes of gut microbiota in adults and children 6 years and older. Mapping of functional orthologs revealed that asthma contributes to 2.9% of the variation in metagenomic content even when accounting for other important clinical demographics. Differential abundance analysis showed an enrichment of long-chain fatty acid (LCFA) metabolism pathways, which have been previously implicated in airway smooth muscle and immune responses in asthma. We also observed increased richness of antibiotic resistance genes (ARGs) in people with asthma. Several differentially abundant ARGs in the asthma cohort encode resistance to macrolide antibiotics, which are often prescribed to patients with asthma. Lastly, we found that ARG and virulence factor (VF) richness in the microbiome were correlated in both cohorts. ARG and VF pairs co-occurred in both cohorts suggesting that virulence and antibiotic resistance traits are coselected and maintained in the fecal microbiota of people with asthma. Overall, our results show functional alterations via LCFA biosynthetic genes and increases in antibiotic resistance genes in the gut microbiota of subjects with moderate-to-severe asthma and could have implications for asthma management and treatment.

Genome-Resolved Metagenomics and Denitrifying Strain Isolation Reveal New Insights into Microbial Denitrification in the Deep Vadose Zone

The heavy use of nitrogen fertilizer in intensive agricultural areas often leads to nitrate accumulation in subsurface soil and nitrate contamination in groundwater, which poses a serious risk to public health. Denitrifying microorganisms in the subsoil convert nitrate to gaseous forms of nitrogen, thereby mitigating the leaching of nitrate into groundwater. Here, we investigated denitrifying microorganisms in the deep vadose zone of a typical intensive agricultural area in China through microcosm enrichment, genome-resolved metagenomic analysis, and denitrifying bacteria isolation. A total of 1000 metagenome-assembled genomes (MAGs) were reconstructed, resulting in 98 high-quality, dereplicated MAGs that contained denitrification genes. Among them, 32 MAGs could not be taxonomically classified at the genus or species level, indicating that a broader spectrum of taxonomic groups is involved in subsoil denitrification than previously recognized. A denitrifier isolate library was constructed by using a strategy combining high-throughput and conventional cultivation techniques. Assessment of the denitrification characteristics of both the MAGs and isolates demonstrated the dominance of truncated denitrification. Functional screening revealed the highest denitrification activity in two complete denitrifiers belonging to the genus Pseudomonas. These findings greatly expand the current knowledge of the composition and function of denitrifying microorganisms in subsoils. The constructed isolate library provided the first pool of subsoil-denitrifying microorganisms that could facilitate the development of microbe-based technologies for nitrate attenuation in groundwater.

Long-term conservation tillage with reduced nitrogen fertilization intensity can improve winter wheat health via positive plant-microorganism feedback in the rhizosphere

Microbiome-based solutions are regarded key for sustainable agroecosystems. However, it is unclear how agricultural practices affect the rhizosphere microbiome, plant-microorganism interactions and crop performance under field conditions. Therefore, we installed root observation windows in a winter wheat field cultivated either under long-term mouldboard plough (MP) or cultivator tillage (CT). Each tillage practice was also compared at two nitrogen (N) fertilization intensities, intensive (recommended N-supply with pesticides/growth regulators) or extensive (reduced N-supply, no fungicides/growth regulators). Shoot biomass, root exudates and rhizosphere metabolites, physiological stress indicators, and gene expression were analyzed together with the rhizosphere microbiome (bacterial/archaeal 16S rRNA gene, fungal ITS amplicon, and shotgun metagenome sequencing) shortly before flowering. Compared to MP, the rhizosphere of CT winter wheat contained more primary and secondary metabolites, especially benzoxazinoid derivatives. Potential copiotrophic and plant-beneficial taxa (e.g. Bacillus, Devosia, and Trichoderma) as well as functional genes (e.g. siderophore production, trehalose synthase, and ACC deaminase) were enriched in the CT rhizosphere, suggesting that tillage affected belowground plant-microorganism interactions. In addition, physiological stress markers were suppressed in CT winter wheat compared to MP. In summary, tillage practice was a major driver of crop performance, root deposits, and rhizosphere microbiome interactions, while the N-fertilization intensity was also relevant, but less important.

Temporal metagenomic characterization of microbial community structure and nitrogen modification genes within an activated sludge bioreactor system

IMPORTANCE

Wastewater treatment plays an essential role in minimizing negative impacts on downstream aquatic environments. Microbial communities are known to play a vital role in the wastewater treatment process, particularly in the removal of nitrogen and phosphorus, which can be especially damaging to aquatic ecosystems. There is limited understanding of how these microbial communities may change in response to fluctuating temperatures or how seasonality may impact their ability to participate in the treatment process. The findings of this study indicate that the microbial communities of wastewater are relatively stable both compositionally and functionally across fluctuating temperatures.

Diversity of Microbial Mats in the Makgadikgadi Salt Pans, Botswana

The Makgadikgadi Salt Pans are the remnants of a mega paleo-lake system in the central Kalahari, Botswana. Today, the Makgadikgadi Basin is an arid to semi-arid area receiving water of meteoric origin during the short, wet season. Large microbial mats, which support primary production, are formed due to desiccation during the dry season. This study aimed to characterise the microbial diversity of the microbial mats and the underlying sediment. The focus was the Ntwetwe Pan, located west of the Makgadikgadi Basin. Metagenomic analyses demonstrated that the mats consisted of a high relative abundance of Cyanobacteriota (synonym Cyanobacteria) (20.50-41.47%), Pseudomonadota (synonym Proteobacteria) (15.71 to 32.18%), and Actinomycetota (synonym Actinobacteria) (8.53-32.56%). In the underlying sediments, Pseudomonadota, Actinomycetota, and Euryarchaeota represented over 70% of the community. Localised fluctuations in water content and pH did not significantly affect the microbial diversity of the sediment or the mats.

Aerobic denitrification as an N2O source from microbial communities

Nitrous oxide (N2O) is a potent greenhouse gas of primarily microbial origin. Oxic and anoxic emissions are commonly ascribed to autotrophic nitrification and heterotrophic denitrification, respectively. Beyond this established dichotomy, we quantitatively show that heterotrophic denitrification can significantly contribute to aerobic nitrogen turnover and N2O emissions in complex microbiomes exposed to frequent oxic/anoxic transitions. Two planktonic, nitrification-inhibited enrichment cultures were established under continuous organic carbon and nitrate feeding, and cyclic oxygen availability. Over a third of the influent organic substrate was respired with nitrate as electron acceptor at high oxygen concentrations (>6.5 mg/L). N2O accounted for up to one-quarter of the nitrate reduced under oxic conditions. The enriched microorganisms maintained a constitutive abundance of denitrifying enzymes due to the oxic/anoxic frequencies exceeding their protein turnover-a common scenario in natural and engineered ecosystems. The aerobic denitrification rates are ascribed primarily to the residual activity of anaerobically synthesised enzymes. From an ecological perspective, the selection of organisms capable of sustaining significant denitrifying activity during aeration shows their competitive advantage over other heterotrophs under varying oxygen availabilities. Ultimately, we propose that the contribution of heterotrophic denitrification to aerobic nitrogen turnover and N2O emissions is currently underestimated in dynamic environments.

Taxonomic and functional stability overrules seasonality in polar benthic microbiomes

Coastal shelf sediments are hot spots of organic matter mineralization. They receive up to 50% of primary production, which, in higher latitudes, is strongly seasonal. Polar and temperate benthic bacterial communities, however, show a stable composition based on comparative 16S rRNA gene sequencing despite different microbial activity levels. Here, we aimed to resolve this contradiction by identifying seasonal changes at the functional level, in particular with respect to algal polysaccharide degradation genes, by combining metagenomics, metatranscriptomics, and glycan analysis in sandy surface sediments from Isfjorden, Svalbard. Gene expressions of diverse carbohydrate-active enzymes changed between winter and spring. For example, β-1,3-glucosidases (e.g. GH30, GH17, GH16) degrading laminarin, an energy storage molecule of algae, were elevated in spring, while enzymes related to α-glucan degradation were expressed in both seasons with maxima in winter (e.g. GH63, GH13_18, and GH15). Also, the expression of GH23 involved in peptidoglycan degradation was prevalent, which is in line with recycling of bacterial biomass. Sugar extractions from bulk sediments were low in concentrations during winter but higher in spring samples, with glucose constituting the largest fraction of measured monosaccharides (84% ± 14%). In porewater, glycan concentrations were ~18-fold higher than in overlying seawater (1107 ± 484 vs. 62 ± 101 μg C l-1) and were depleted in glucose. Our data indicate that microbial communities in sandy sediments digest and transform labile parts of photosynthesis-derived particulate organic matter and likely release more stable, glucose-depleted residual glycans of unknown structures, quantities, and residence times into the ocean, thus modulating the glycan composition of marine coastal waters.

Experimental evolution at ecological scales allows linking of viral genotypes to specific host strains

Viruses shape microbial community structure and activity through the control of population diversity and cell abundances. Identifying and monitoring the dynamics of specific virus-host pairs in nature is hampered by the limitations of culture-independent approaches such as metagenomics, which do not always provide strain-level resolution, and culture-based analyses, which eliminate the ecological background and in-situ interactions. Here, we have explored the interaction of a specific "autochthonous" host strain and its viruses within a natural community. Bacterium Salinibacter ruber strain M8 was spiked into its environment of isolation, a crystallizer pond from a coastal saltern, and the viral and cellular communities were monitored for one month using culture, metagenomics, and microscopy. Metagenome sequencing indicated that the M8 abundance decreased sharply after being added to the pond, likely due to forces other than viral predation. However, the presence of M8 selected for two species of a new viral genus, Phoenicisalinivirus, for which 120 strains were isolated. During this experiment, an assemblage of closely related viral genomic variants was replaced by a single population with the ability to infect M8, a scenario which was compatible with the selection of a genomic variant from the rare biosphere. Further analysis implicated a viral genomic region putatively coding for a tail fiber protein to be responsible for M8 specificity. Our results indicate that low abundance viral genotypes provide a viral seed bank that allows for a highly specialized virus-host response within a complex ecological background.

Recruitment of complete crAss-like phage genomes reveals their presence in chicken viromes, few human-specific phages, and lack of universal detection

The order Crassvirales, which includes the prototypical crAssphage (p-crAssphage), is predominantly associated with humans, rendering it the most abundant and widely distributed group of DNA phages in the human gut. The reported human specificity and wide global distribution of p-crAssphage makes it a promising human fecal marker. However, the specificity for the human gut as well as the geographical distribution around the globe of other members of the order Crassvirales remains unknown. To determine this, a recruitment analysis using 91 complete, non-redundant genomes of crAss-like phages in human and animal viromes revealed that only 13 crAss-like phages among the 91 phages analyzed were highly specific to humans, and p-crAssphage was not in this group. Investigations to elucidate whether any characteristic of the phages was responsible for their prevalence in humans showed that the 13 human crAss-like phages do not share a core genome. Phylogenomic analysis placed them in three independent families, indicating that within the Crassvirales group, human specificity is likely not a feature of a common ancestor but rather was introduced on separate/independent occasions in their evolutionary history. The 13 human crAss-like phages showed variable geographical distribution across human metagenomes worldwide, with some being more prevalent in certain countries than in others, but none being universally identified. The varied geographical distribution and the absence of a phylogenetic relationship among the human crAss-like phages are attributed to the emergence and dissemination of their bacterial host, the symbiotic human strains of Bacteroides, across various human populations occupying diverse ecological niches worldwide.

Why are so many enteric pathogen infections asymptomatic? Pathogen and gut microbiome characteristics associated with diarrhea symptoms and carriage of diarrheagenic E. coli in northern Ecuador

A high proportion of enteric infections, including those caused by diarrheagenic Escherichia coli (DEC), are asymptomatic for diarrhea. The factors responsible for the development of diarrhea symptoms, or lack thereof, remain unclear. Here, we used DEC isolate genome and whole stool microbiome data from a case-control study of diarrhea in Ecuador to examine factors associated with diarrhea symptoms accompanying DEC carriage. We investigated i) pathogen abundance, ii) gut microbiome characteristics, and iii) strain-level pathogen characteristics from DEC infections with diarrhea symptoms (symptomatic infections) and without diarrhea symptoms (asymptomatic infections). We also included data from individuals with and without diarrhea who were not infected with DEC (uninfected cases and controls). i) E. coli relative abundance in the gut microbiome was highly variable, but higher on-average in individuals with symptomatic compared to asymptomatic DEC infections. Similarly, the number and relative abundances of virulence genes in the gut were higher in symptomatic than asymptomatic DEC infections. ii) Measures of microbiome diversity were similar regardless of diarrhea symptoms or DEC carriage. Proteobacterial families that have been described as pathobionts were enriched in symptomatic infections and uninfected cases, whereas potentially beneficial taxa, including the Bacteroidaceae and Bifidobacteriaceae, were more abundant in individuals without diarrhea. An analysis of high-level gene functions recovered in metagenomes revealed that genes that were differentially abundant by diarrhea and DEC infection status were more abundant in symptomatic than asymptomatic DEC infections. iii) DEC isolates from symptomatic versus asymptomatic individuals showed no significant differences in virulence or accessory gene content, and there was no phylogenetic signal associated with diarrhea symptoms. Together, these data suggest signals that distinguish symptomatic from asymptomatic DEC infections. In particular, the abundance of E. coli, the virulence gene content of the gut microbiome, and the taxa present in the gut microbiome have an apparent role.

Seasonal variations of biofilm C, N and S cycling genes in a pilot-scale chlorinated drinking water distribution system

Biofilms in drinking water distribution systems (DWDS) host diverse microorganisms. However, the functional attributes of DWDS biofilms and their associations with seasonality remain unclear. This study aims to characterize variations in the microbial metabolic traits of DWDS biofilms collected during different seasons, using a pilot-scale DWDS in dark under plug-flow conditions during one-year operation period. Network analysis was used to predict the functional gene hosts. The overall functional attributes determined by shotgun metagenomics exhibited significant differences among seasons. Genes associated with aromatic metabolism, fatty acid biosynthesis and degradation, and capsular extracellular polymeric substance (EPS) were significantly upregulated in summer owing to the higher temperatures and chlorine in the influent of the DWDS. Moreover, the pathways associated with nitrogen, sulfur, glycolysis, and tricarboxylic acid (TCA) cycling, as well as carbon fixation were reconstructed and displayed according to the sampling season. Nitrogen reduction pathways [dissimilatory nitrate reduction to ammonium (DNRA) 73 %, assimilatory nitrate reduction to ammonium (ANRA) 21 %] were identified in DWDS biofilms, but nitrogen oxidation pathways were not. Sulfur cycling were involved in diverse pathways and genes. Glycolysis and TCA cycling offered electron donors and energy sources for nitrogen and sulfur reduction in biofilms. Carbon fixation was observed in DWDS biofilms, with the predominant pathway for fixing carbon dioxide being the reductive citrate cycle (38 %). Constructed functional gene networks composed of carbon, nitrogen, and sulfur cycling-related genes demonstrated synergistic effects (Positive proportion: 63.52-71.09 %). In addition, from spring to autumn, the network complexity decreased and network modularity increased. The assembly mechanism of carbon, nitrogen and sulfur cycling-related genes was driven by stochastic processes for all samples. These results highlight the diverse functional genes in DWDS biofilms, their synergetic interrelationships, and the seasonality effect on functional attributes.

Impact of dilution on stochastically driven methanogenic microbial communities of hypersaline anoxic sediments

Sediments underlying the solar salterns of S'Avall are anoxic hypersaline ecosystems dominated by anaerobic prokaryotes, and with the especial relevance of putative methanogenic archaea. Slurries from salt-saturated sediments, diluted in a gradient of salinity and incubated for > 4 years revealed that salt concentration was the major selection force that deterministically structured microbial communities. The dominant archaea in the original communities showed a decrease in alpha diversity with dilution accompanied by the increase of bacterial alpha diversity, being highest at 5% salts. Correspondingly, methanogens decreased and in turn sulfate reducers increased with decreasing salt concentrations. Methanogens especially dominated at 25%. Different concentrations of litter of Posidonia oceanica seagrass added as a carbon substrate, did not promote any clear relevant effect. However, the addition of ampicillin as selection pressure exerted important effects on the assemblage probably due to the removal of competitors or enhancers. The amended antibiotic enhanced methanogenesis in the concentrations ≤ 15% of salts, whereas it was depleted at salinities ≥ 20% revealing key roles of ampicillin-sensitive bacteria.

A glimpse of the paleome in endolithic microbial communities

BACKGROUND

The terrestrial subsurface is home to a significant proportion of the Earth's microbial biomass. Our understanding about terrestrial subsurface microbiomes is almost exclusively derived from groundwater and porous sediments mainly by using 16S rRNA gene surveys. To obtain more insights about biomass of consolidated rocks and the metabolic status of endolithic microbiomes, we investigated interbedded limestone and mudstone from the vadose zone, fractured aquifers, and deep aquitards.

RESULTS

By adapting methods from microbial archaeology and paleogenomics, we could recover sufficient DNA for downstream metagenomic analysis from seven rock specimens independent of porosity, lithology, and depth. Based on the extracted DNA, we estimated between 2.81 and 4.25 × 105 cells × g-1 rock. Analyzing DNA damage patterns revealed paleome signatures (genetic records of past microbial communities) for three rock specimens, all obtained from the vadose zone. DNA obtained from deep aquitards isolated from surface input was not affected by DNA decay indicating that water saturation and not flow is controlling subsurface microbial survival. Decoding the taxonomy and functional potential of paleome communities revealed increased abundances for sequences affiliated with chemolithoautotrophs and taxa such as Cand. Rokubacteria. We also found a broader metabolic potential in terms of aromatic hydrocarbon breakdown, suggesting a preferred utilization of sedimentary organic matter in the past.

CONCLUSIONS

Our study suggests that limestones function as archives for genetic records of past microbial communities including those sensitive to environmental stress at modern times, due to their specific conditions facilitating long-term DNA preservation. Video Abstract.

Recovery of the gut microbiome following enteric infection and persistence of antimicrobial resistance genes in specific microbial hosts

Enteric pathogens cause widespread foodborne illness and are increasingly resistant to important antibiotics yet their ecological impact on the gut microbiome and resistome is not fully understood. Herein, shotgun metagenome sequencing was applied to stool DNA from 60 patients (cases) during an enteric bacterial infection and after recovery (follow-ups). Overall, the case samples harbored more antimicrobial resistance genes (ARGs) with greater resistome diversity than the follow-up samples (p < 0.001), while follow-ups had more diverse gut microbiota (p < 0.001). Although cases were primarily defined by genera Escherichia, Salmonella, and Shigella along with ARGs for multi-compound and multidrug resistance, follow-ups had a greater abundance of Bacteroidetes and Firmicutes phyla and resistance genes for tetracyclines, macrolides, lincosamides, and streptogramins, and aminoglycosides. A host-tracking analysis revealed that Escherichia was the primary bacterial host of ARGs in both cases and follow-ups, with a greater abundance occurring during infection. Eleven distinct extended spectrum beta-lactamase (ESBL) genes were identified during infection, with some detectable upon recovery, highlighting the potential for gene transfer within the community. Because of the increasing incidence of disease caused by foodborne pathogens and their role in harboring and transferring resistance determinants, this study enhances our understanding of how enteric infections impact human gut ecology.

Seasonality Determines the Variations of Biofilm Microbiome and Antibiotic Resistome in a Pilot-Scale Chlorinated Drinking Water Distribution System Deciphered by Metagenome Assembly

Understanding the biofilm microbiome and antibiotic resistome evolution in drinking water distribution systems (DWDSs) is crucial to ensure the safety of drinking water. We explored the 10 month evolution of the microbial community, antibiotic resistance genes (ARGs), mobile gene elements (MGEs) co-existing with ARGs and pathogenic ARG hosts, and the ARG driving factors in DWDS biofilms using metagenomics assembly. Sampling season was critical in determining the microbial community and antibiotic resistome shift. Pseudomonas was the primary biofilm colonizer, and biofilms diversified more as the formation time increased. Most genera tended to cooperate to adapt to an oligotrophic environment with disinfectant stress. Biofilm microbial community and antibiotic resistome assembly were mainly determined by stochastic processes and changed with season. Metagenome assembly provided the occurrence and fates of MGEs co-existing with ARGs and ARG hosts in DWDS biofilms. The abundance of ARG- and MGE-carrying pathogen Stenotrophomonas maltophilia was high in summer. It primarily harbored the aph(3)-IIb, multidrug transporter, smeD, and metallo-beta-lactamase ARGs, which were transferred via recombination. The microbial community was the most crucial factor driving the antibiotic resistance shift. We provide novel insights about the evolution of pathogens and ARGs and their correlations in DWDS biofilms to ensure the safety of drinking water.

Altered infective competence of the human gut microbiome in COVID-19

BACKGROUND

Infections with SARS-CoV-2 have a pronounced impact on the gastrointestinal tract and its resident microbiome. Clear differences between severe cases of infection and healthy individuals have been reported, including the loss of commensal taxa. We aimed to understand if microbiome alterations including functional shifts are unique to severe cases or a common effect of COVID-19. We used high-resolution systematic multi-omic analyses to profile the gut microbiome in asymptomatic-to-moderate COVID-19 individuals compared to a control group.

RESULTS

We found a striking increase in the overall abundance and expression of both virulence factors and antimicrobial resistance genes in COVID-19. Importantly, these genes are encoded and expressed by commensal taxa from families such as Acidaminococcaceae and Erysipelatoclostridiaceae, which we found to be enriched in COVID-19-positive individuals. We also found an enrichment in the expression of a betaherpesvirus and rotavirus C genes in COVID-19-positive individuals compared to healthy controls.

CONCLUSIONS

Our analyses identified an altered and increased infective competence of the gut microbiome in COVID-19 patients. Video Abstract.

Deduplication Improves Cost-Efficiency and Yields of De Novo Assembly and Binning of Shotgun Metagenomes in Microbiome Research

In the last decade, metagenomics has greatly revolutionized the study of microbial communities. However, the presence of artificial duplicate reads raised mainly from the preparation of metagenomic DNA sequencing libraries and their impacts on metagenomic assembly and binning have never been brought to attention. Here, we explicitly investigated the effects of duplicate reads on metagenomic assemblies and binning based on analyses of five groups of representative metagenomes with distinct microbiome complexities. Our results showed that deduplication considerably increased the binning yields (by 3.5% to 80%) for most of the metagenomic data sets examined thanks to the improved contig length and coverage profiling of metagenome-assembled contigs, whereas it slightly decreased the binning yields of metagenomes with low complexity (e.g., human gut metagenomes). Specifically, 411 versus 397, 331 versus 317, 104 versus 88, and 9 versus 5 metagenome-assembled genomes (MAGs) were recovered from MEGAHIT assemblies of bioreactor sludge, surface water, lake sediment, and forest soil metagenomes, respectively. Noticeably, deduplication significantly reduced the computational costs of the metagenomic assembly, including the elapsed time (9.0% to 29.9%) and the maximum memory requirement (4.3% to 37.1%). Collectively, we recommend the removal of duplicate reads in metagenomes with high complexity before assembly and binning analyses, for example, the forest soil metagenomes examined in this study. IMPORTANCE Duplicated reads in shotgun metagenomes are usually considered technical artifacts. Their presence in metagenomes would theoretically not only introduce bias into the quantitative analysis but also result in mistakes in the coverage profile, leading to adverse effects on or even failures in metagenomic assembly and binning, as the widely used metagenome assemblers and binners all need coverage information for graph partitioning and assembly binning, respectively. However, this issue was seldom noticed, and its impacts on downstream essential bioinformatic procedures (e.g., assembly and binning) remained unclear. In this study, we comprehensively evaluated for the first time the implications of duplicate reads for the de novo assembly and binning of real metagenomic data sets by comparing the assembly qualities, binning yields, and requirements for computational resources with and without the removal of duplicate reads. It was revealed that deduplication considerably increased the binning yields of metagenomes with high complexity and significantly reduced the computational costs, including the elapsed time and the maximum memory requirement, for most of the metagenomes studied. These results provide empirical references for more cost-efficient metagenomic analyses in microbiome research.

The gut metagenome harbors metabolic and antibiotic resistance signatures of moderate-to-severe asthma

Asthma is a common allergic airway disease that develops in association with the human microbiome early in life. Both the composition and function of the infant gut microbiota have been linked to asthma risk, but functional alterations in the gut microbiota of older patients with established asthma remain an important knowledge gap. Here, we performed whole metagenomic shotgun sequencing of 95 stool samples from 59 healthy and 36 subjects with moderate-to-severe asthma to characterize the metagenomes of gut microbiota in children and adults 6 years and older. Mapping of functional orthologs revealed that asthma contributes to 2.9% of the variation in metagenomic content even when accounting for other important clinical demographics. Differential abundance analysis showed an enrichment of long-chain fatty acid (LCFA) metabolism pathways which have been previously implicated in airway smooth muscle and immune responses in asthma. We also observed increased richness of antibiotic resistance genes (ARGs) in people with asthma. One differentially abundant ARG was a macrolide resistance marker, ermF, which significantly co-occurred with the Bacteroides fragilis toxin, suggesting a possible relationship between enterotoxigenic B. fragilis, antibiotic resistance, and asthma. Lastly, we found multiple virulence factor (VF) and ARG pairs that co-occurred in both cohorts suggesting that virulence and antibiotic resistance traits are co-selected and maintained in the fecal microbiota of people with asthma. Overall, our results show functional alterations via LCFA biosynthetic genes and increases in antibiotic resistance genes in the gut microbiota of subjects with moderate-to-severe asthma and could have implications for asthma management and treatment.

Salt tolerance evolution facilitates antibiotic resistome in soil microbiota: Evidences from dissemination evaluation, hosts identification and co-occurrence exploration

Salinity is considered as one of the vital factors affecting the profiles of antibiotic resistance genes (ARGs) in soils, whereby its roles in shaping the antibiotic resistome were still poorly understood. Here, metagenomic analysis was conducted to track the ARGs distributions and dissemination in soils during salt accumulation and desalinization processes. Neutral-salt accumulation for 45 and 90 days significantly increased the relative abundances of ARGs and mobile genetic elements (MGEs) carrying antibiotic resistance contigs (ARCs). The ARGs within antibiotic efflux and target protection families primarily carried by Streptomyces, Nocardioides, Rhodanobacter and Monashia were largely enriched by salinity. The ARGs subtypes of the resistance-nodulation-division (RND) family, ATP-binding cassette (ABC) family, rRNA methyltransferase and other efflux were closely associated with MGEs, contributing to the enrichment of ARGs. Moreover, the ARGs subtypes and transposons were genetically linked with the salt-tolerance mechanisms of organic osmolyte transporters and K⁺ uptake proteins on the same ARC, demonstrating the coselection of ARGs and halotolerant genes. Furthermore, the antibiotic resistome could recover to a normal state after the prolonged incubation by alleviating salt stress. Nevertheless, the acquisition of ARGs by opportunistic pathogens after salt treatment was increased, serving to prioritize further efforts on the health risks correlated with resistance propagation and human exposure in saline soils.

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.

Microbial carbon, sulfur, iron, and nitrogen cycling linked to the potential remediation of a meromictic acidic pit lake

Cueva de la Mora is a permanently stratified acidic pit lake and a model system for extreme acid mine drainage (AMD) studies. Using a combination of amplicon sequencing, metagenomics and metatranscriptomics we performed a taxonomically resolved analysis of microbial contributions to carbon, sulfur, iron, and nitrogen cycling. We found that active green alga Coccomyxa onubensis dominated the upper layer and chemocline. The chemocline had activity for iron(II) oxidation carried out by populations of Ca. Acidulodesulfobacterium, Ferrovum, Leptospirillium, and Armatimonadetes. Predicted activity for iron(III) reduction was only detected in the deep layer affiliated with Proteobacteria. Activity for dissimilatory nitrogen cycling including nitrogen fixation and nitrate reduction was primarily predicted in the chemocline. Heterotrophic archaeal populations with predicted activity for sulfide oxidation related to uncultured Thermoplasmatales dominated in the deep layer. Abundant sulfate-reducing Desulfomonile and Ca. Acidulodesulfobacterium populations were active in the chemocline. In the deep layer, uncultured populations from the bacterial phyla Actinobacteria, Chloroflexi, and Nitrospirae contributed to both sulfate reduction and sulfide oxidation. Based on this information we evaluated the potential for sulfide mineral precipitation in the deep layer as a tool for remediation. We argue that sulfide precipitation is not limited by microbial genetic potential but rather by the quantity and quality of organic carbon reaching the deep layer as well as by oxygen additions to the groundwater enabling sulfur oxidation. Addition of organic carbon and elemental sulfur should stimulate sulfate reduction and limit reoxidation of sulfide minerals.

Metagenomic insights into the functions of microbial communities in sulfur-rich sediment of a shallow-water hydrothermal vent off Kueishan Island

Hydrothermal vent (HTV) systems are important habitats for understanding the biological processes of extremophiles on Earth and their relative contributions to material and energy cycles in the ocean. Current understanding on hydrothermal systems have been primarily focused on deep-sea HTVs, and little is known about the functions and metabolisms of microorganisms in shallow-water HTVs (SW-HTVs), which are distinguished from deep-sea HTVs by a depth limit of 200 m. In this study, we analyzed metagenomes of sulfur-rich sediment samples collected from a SW-HTV of Kueishan Island, located in a marginal sea of the western Pacific Ocean. Comparing with a previously published report of pelagic samples from the nearby sampling site, microbial communities in the SW-HTV sediments enriching with genes of both aerobic and anaerobic respiration inferred variable environments in the tested sediments. Abundant genes of energy metabolism encoding sulfur oxidation, H2 oxidation, and carbon fixation were detected from the sediment samples. Sixty-eight metagenome-assembled-genomes (MAGs) were reconstructed to further understand the metabolism and potential interactions between different microbial taxa in the SW-HTVs sediment. MAGs with the highest abundant were chemolithotrophic sulfur-oxidization bacteria, including Sulfurovum represented Campylobacteria involved sox multienzyme, sulfide oxidation genes and rTCA cycle, and Gammaproteobacteria involved dsr gene and CBB cycle. In addition, Desulfobacterota with the potential to participate in sulfur-disproportionating processes also had higher abundance than the sample’s overall mean value. The interaction of these bacterial groups allows the microbial communities to efficiently metabolize a large variety of sulfur compounds. In addition, the potential to use simple organic carbon, such as acetate, was found in chemolithotrophic Campylobacterial MAGs. Collectively, our results revealed the complexity of environmental conditions of the vent sediment and highlight the interactive relationships of the dominant microbial populations in driving sulfur cycles in the SW-HTV sediments off Kueishan Island.

Key players in syntrophic propionate oxidation revealed by metagenome-assembled genomes from anaerobic digesters bioaugmented with propionic acid enriched microbial consortia

Propionic acid (HPr) is frequently accumulated in anaerobic digesters due to its thermodynamically unfavorable degradation reaction. Here, we identify key players in HPr oxidation and organic overloading recovery from metagenome-assembled genomes (MAGs) recovered from anaerobic digesters inoculated with HPr-enriched microbial consortia before initiating organic overloading. Two independent HPr-enrichment cultures commonly selected two uncultured microorganisms represented with high relative abundance: Methanoculleus sp002497965 and JABUEY01 sp013314815 (a member of the Syntrophobacteraceae family). The relative abundance of JABUEY01 sp013314815 was 60 times higher in bioaugmented bioreactors compared to their unaugmented counterparts after recovery from organic overloading. Genomic analysis of JABUEY01 sp013314815 revealed its metabolic potential for syntrophic propionate degradation when partnered with hydrogenotrophic methanogens (e.g., Methanoculleus sp002497965) via the methylmalonyl-CoA pathway. Our results identified at least two key species that are responsible for efficient propionate removal and demonstrate their potential applications as microbial cocktails for stable AD operation.

Globally distributed mining-impacted environments are underexplored hotspots of multidrug resistance genes

Mining is among the human activities with widest environmental impacts, and mining-impacted environments are characterized by high levels of metals that can co-select for antibiotic resistance genes (ARGs) in microorganisms. However, ARGs in mining-impacted environments are still poorly understood. Here, we conducted a comprehensive study of ARGs in such environments worldwide, taking advantage of 272 metagenomes generated from a global-scale data collection and two national sampling efforts in China. The average total abundance of the ARGs in globally distributed studied mine sites was 1572 times per gigabase, being rivaling that of urban sewage but much higher than that of freshwater sediments. Multidrug resistance genes accounted for 40% of the total ARG abundance, tended to co-occur with multimetal resistance genes, and were highly mobile (e.g. on average 16% occurring on plasmids). Among the 1848 high-quality metagenome-assembled genomes (MAGs), 85% carried at least one multidrug resistance gene plus one multimetal resistance gene. These high-quality ARG-carrying MAGs considerably expanded the phylogenetic diversity of ARG hosts, providing the first representatives of ARG-carrying MAGs for the Archaea domain and three bacterial phyla. Moreover, 54 high-quality ARG-carrying MAGs were identified as potential pathogens. Our findings suggest that mining-impacted environments worldwide are underexplored hotspots of multidrug resistance genes.

Atmospheric chemosynthesis is phylogenetically and geographically widespread and contributes significantly to carbon fixation throughout cold deserts

Cold desert soil microbiomes thrive despite severe moisture and nutrient limitations. In Eastern Antarctic soils, bacterial primary production is supported by trace gas oxidation and the light-independent RuBisCO form IE. This study aims to determine if atmospheric chemosynthesis is widespread within Antarctic, Arctic and Tibetan cold deserts, to identify the breadth of trace gas chemosynthetic taxa and to further characterize the genetic determinants of this process. H₂ oxidation was ubiquitous, far exceeding rates reported to fulfill the maintenance needs of similarly structured edaphic microbiomes. Atmospheric chemosynthesis occurred globally, contributing significantly (p < 0.05) to carbon fixation in Antarctica and the high Arctic. Taxonomic and functional analyses were performed upon 18 cold desert metagenomes, 230 dereplicated medium-to-high-quality derived metagenome-assembled genomes (MAGs) and an additional 24,080 publicly available genomes. Hydrogenotrophic and carboxydotrophic growth markers were widespread. RuBisCO IE was discovered to co-occur alongside trace gas oxidation enzymes in representative Chloroflexota, Firmicutes, Deinococcota and Verrucomicrobiota genomes. We identify a novel group of high-affinity [NiFe]-hydrogenases, group 1m, through phylogenetics, gene structure analysis and homology modeling, and reveal substantial genetic diversity within RuBisCO form IE (rbcL1E), and high-affinity 1h and 1l [NiFe]-hydrogenase groups. We conclude that atmospheric chemosynthesis is a globally-distributed phenomenon, extending throughout cold deserts, with significant implications for the global carbon cycle and bacterial survival within environmental reservoirs.

A consensus protocol for the recovery of mercury methylation genes from metagenomes

Mercury (Hg) methylation genes (hgcAB) mediate the formation of the toxic methylmercury and have been identified from diverse environments, including freshwater and marine ecosystems, Arctic permafrost, forest and paddy soils, coal-ash amended sediments, chlor-alkali plants discharges and geothermal springs. Here we present the first attempt at a standardized protocol for the detection, identification and quantification of hgc genes from metagenomes. Our Hg-MATE (Hg-cycling Microorganisms in Aquatic and Terrestrial Ecosystems) database, a catalogue of hgc genes, provides the most accurate information to date on the taxonomic identity and functional/metabolic attributes of microorganisms responsible for Hg methylation in the environment. Furthermore, we introduce "marky-coco", a ready-to-use bioinformatic pipeline based on de novo single-metagenome assembly, for easy and accurate characterization of hgc genes from environmental samples. We compared the recovery of hgc genes from environmental metagenomes using the marky-coco pipeline with an approach based on co-assembly of multiple metagenomes. Our data show similar efficiency in both approaches for most environments except those with high diversity (i.e., paddy soils) for which a co-assembly approach was preferred. Finally, we discuss the definition of true hgc genes and methods to normalize hgc gene counts from metagenomes.

Uncovering Hidden Members and Functions of the Soil Microbiome Using De Novo Metaproteomics

Metaproteomics has been increasingly utilized for high-throughput characterization of proteins in complex environments and has been demonstrated to provide insights into microbial composition and functional roles. However, significant challenges remain in metaproteomic data analysis, including creation of a sample-specific protein sequence database. A well-matched database is a requirement for successful metaproteomics analysis, and the accuracy and sensitivity of PSM identification algorithms suffer when the database is incomplete or contains extraneous sequences. When matched DNA sequencing data of the sample is unavailable or incomplete, creating the proteome database that accurately represents the organisms in the sample is a challenge. Here, we leverage a de novo peptide sequencing approach to identify the sample composition directly from metaproteomic data. First, we created a deep learning model, Kaiko, to predict the peptide sequences from mass spectrometry data and trained it on 5 million peptide–spectrum matches from 55 phylogenetically diverse bacteria. After training, Kaiko successfully identified organisms from soil isolates and synthetic communities directly from proteomics data. Finally, we created a pipeline for metaproteome database generation using Kaiko. We tested the pipeline on native soils collected in Kansas, showing that the de novo sequencing model can be employed as an alternative and complementary method to construct the sample-specific protein database instead of relying on (un)matched metagenomes. Our pipeline identified all highly abundant taxa from 16S rRNA sequencing of the soil samples and uncovered several additional species which were strongly represented only in proteomic data.

Deep-Sea Sediments from the Southern Gulf of Mexico Harbor a Wide Diversity of PKS I Genes

The excessive use of antibiotics has triggered the appearance of new resistant strains, which is why great interest has been taken in the search for new bioactive compounds capable of overcoming this emergency in recent years. Massive sequencing tools have enabled the detection of new microorganisms that cannot be cultured in a laboratory, thus opening the door to the search for new biosynthetic genes. The great variety in oceanic environments in terms of pressure, salinity, temperature, and nutrients enables marine microorganisms to develop unique biochemical and physiological properties for their survival, enhancing the production of secondary metabolites that can vary from those produced by terrestrial microorganisms. We performed a search for type I PKS genes in metagenomes obtained from the marine sediments of the deep waters of the Gulf of Mexico using Hidden Markov Models. More than 2000 candidate genes were detected in the metagenomes that code for type I PKS domains, while biosynthetic pathways that may code for other secondary metabolites were also detected. Our research demonstrates the great potential use of the marine sediments of the Gulf of Mexico for identifying genes that code for new secondary metabolites.

Metagenomics indicate that public health risk may be higher from flooding following dry versus rainy periods

Urban floodwater could lead to significant risk for public and environmental health from mobilization of microbial pathogens and overflow of wastewater treatment systems. Here, we attempted to assess this risk by obtaining metagenomic profiles of antibiotic resistance genes (ARGs), virulence factors (VFs) and pathogens present in floodwater samples collected in urban Atlanta, GA that were categorized in two distinct groups: floods that occurred after periods of drought and those after regular (seasonal) rain events. Even though no major (known) pathogens were present at the limit of detection of our sequencing effort (~3 Gbp/sample), we observed that floodwaters after drought showed a 2.5-fold higher abundance of both ARGs and VFs compared to floodwater after rainy days. These differences were mainly derived by several novel species of the Pseudomonas genus, which were more dominant in the former versus the latter samples and carried several genes to cope with osmotic stress in addition to ARGs and VFs. These results revealed that there are previously undescribed species that become mobilized after flooding events in the Southeast US urban settings and could represent an increased public health risk, especially after periods of drought, which warrants further attention.

Identification of key microorganisms in facultative stabilization ponds from dairy industries, using metagenomics

Wastewater stabilization ponds are a natural form of wastewater treatment. Their low operation and maintenance costs have made them popular, especially in developing countries. In these systems, effluents are retained for long periods of time, allowing the microbial communities present in the ponds to degrade the organic matter present, using both aerobic and anaerobic processes. Even though these systems are widespread in low income countries, there are no studies about the microorganisms present in them and how they operate. In this study, we analised the microbial communities of two serial full-scale stabilization ponds systems using whole genome shotgun sequencing. First, a taxonomic profiling of the reads was performed, to estimate the microbial diversity. Then, the reads of each system were assembled and binned, allowing the reconstruction of 110 microbial genomes. A functional analysis of the genomes allowed us to find how the main metabolic pathways are carried out, and we propose several organisms that would be key to this kind of environment, since they play an important role in these metabolic pathways. This study represents the first genome-centred approach to understand the metabolic processes in facultative ponds. A better understanding of these microbial communities and how they stabilize the effluents of dairy industries is necessary to improve them and to minimize the environmental impact of dairy industries wastewater.

Clinical Metagenomics Is Increasingly Accurate and Affordable to Detect Enteric Bacterial Pathogens in Stool

Stool culture is the gold standard method to diagnose enteric bacterial infections; however, many clinical laboratories are transitioning to syndromic multiplex PCR panels. PCR is rapid, accurate, and affordable, yet does not yield subtyping information critical for foodborne disease surveillance. A metagenomics-based stool testing approach could simultaneously provide diagnostic and public health information. Here, we evaluated shotgun metagenomics to assess the detection of common enteric bacterial pathogens in stool. We sequenced 304 stool specimens from 285 patients alongside routine diagnostic testing for Salmonella spp., Campylobacter spp., Shigella spp., and shiga-toxin producing Escherichia coli. Five analytical approaches were assessed for pathogen detection: microbiome profiling, Kraken2, MetaPhlAn, SRST2, and KAT-SECT. Among analysis tools and databases compared, KAT-SECT analysis provided the best sensitivity and specificity for all pathogens tested compared to culture (91.2% and 96.2%, respectively). Where metagenomics detected a pathogen in culture-negative specimens, standard PCR was positive 85% of the time. The cost of metagenomics is approaching the current combined cost of PCR, reflex culture, and whole genome sequencing for pathogen detection and subtyping. As cost, speed, and analytics for single-approach metagenomics improve, it may be more routinely applied in clinical and public health laboratories.

Assessment of Hydrocarbon Degradation Potential in Microbial Communities in Arctic Sea Ice

The anthropogenic release of oil hydrocarbons into the cold marine environment is an increasing concern due to the elevated usage of sea routes and the exploration of new oil drilling sites in Arctic areas. The aim of this study was to evaluate prokaryotic community structures and the genetic potential of hydrocarbon degradation in the metagenomes of seawater, sea ice, and crude oil encapsulating the sea ice of the Norwegian fjord, Ofotfjorden. Although the results indicated substantial differences between the structure of prokaryotic communities in seawater and sea ice, the crude oil encapsulating sea ice (SIO) showed increased abundances of many genera-containing hydrocarbon-degrading organisms, including Bermanella, Colwellia, and Glaciecola. Although the metagenome of seawater was rich in a variety of hydrocarbon degradation-related functional genes (HDGs) associated with the metabolism of n-alkanes, and mono- and polyaromatic hydrocarbons, most of the normalized gene counts were highest in the clean sea ice metagenome, whereas in SIO, these counts were the lowest. The long-chain alkane degradation gene almA was detected from all the studied metagenomes and its counts exceeded ladA and alkB counts in both sea ice metagenomes. In addition, almA was related to the most diverse group of prokaryotic genera. Almost all 18 good- and high-quality metagenome-assembled genomes (MAGs) had diverse HDGs profiles. The MAGs recovered from the SIO metagenome belonged to the abundant taxa, such as Glaciecola, Bermanella, and Rhodobacteracea, in this environment. The genera associated with HDGs were often previously known as hydrocarbon-degrading genera. However, a substantial number of new associations, either between already known hydrocarbon-degrading genera and new HDGs or between genera not known to contain hydrocarbon degraders and multiple HDGs, were found. The superimposition of the results of comparing HDG associations with taxonomy, the HDG profiles of MAGs, and the full genomes of organisms in the KEGG database suggest that the found relationships need further investigation and verification.

Functional filtering and random processes affect the assembly of microbial communities of snow algae blooms at Maritime Antarctic

The increasing temperatures at the West Antarctic Peninsula (Maritime Antarctic) could lead to a higher occurrence of snow algal blooms which are ubiquitous events that change the snow coloration, reducing albedo and in turn exacerbating melting. However, there is a limited understanding of snow algae blooms biodiversity, composition, and their functional profiles, especially in one of the world's areas most affected by climate change. In this study we used 16S rRNA and 18S rRNA metabarcoding, and shotgun metagenomics to assess the diversity, composition, and functional potential of the snow algae blooms bacterial and eukaryotic communities at three different sites of Maritime Antarctic, between different colors of the algae blooms and between seasonal and semi-permanent snowfields. We tested the hypothesis that the functional potential of snow algae blooms is conserved despite a changing taxonomic composition. Furthermore, we determined taxonomic co-occurrence patterns of bacteria and eukaryotes and assessed the potential for the exchange of metabolites among bacterial taxa. Here, we tested the prediction that there are co-occurring taxa within snow algae whose biotic interactions are marked by the exchange of metabolites. Our results show that the composition of snow algae blooms vary significantly among sites. For instance, a higher abundance of fungi and protists were detected in Fildes Peninsula compared with Doumer Island and O'Higgins. Likewise, the composition varied between snow colors and snow types. However, the functional potential varied only among sampling sites with a higher abundance of genes involved in tolerance to environmental stress at O'Higgins. Co-occurrence patterns of dominant bacterial genera such as Pedobacter, Polaromonas, Flavobacterium and Hymenobacter were recorded, contrasting the absence of co-occurring patterns displayed by Chlamydomonadales algae with other eukaryotes. Finally, genome-scale metabolic models revealed that bacteria within snow algae blooms likely compete for resources instead of forming cooperative communities.