Community-wide analysis of microbial genome sequence signatures

Genome-resolved metagenomics reveals the effect of nutrient availability on bacterial genomic properties across 44 European freshwater lakes

Understanding intricate microbial interactions in the environment is crucial. This is especially true for the relationships between nutrients and bacteria, as phosphorus, nitrogen and organic carbon availability are known to influence bacterial population dynamics. It has been suggested that low nutrient conditions prompt the evolutionary process of genome streamlining. This process helps conserve scarce nutrients and allows for proliferation. Genome streamlining is associated with genomic properties such as %GC content, genes encoding sigma factors, percent coding regions, gene redundancy, and functional shifts in processes like cell motility and ATP binding cassette transporters, among others. The current study aims to unveil the impact of nutrition on the genome size, %GC content, and functional properties of pelagic freshwater bacteria. We do this at finer taxonomic resolutions for many metagenomically characterized communities. Our study confirms the interplay of trophic level and genomic properties. It also highlights that different nutrient types, particularly phosphorus and nitrogen, impact these properties differently. We observed a covariation of functional traits with genome size. Larger genomes exhibit enriched pathways for motility, environmental interaction, and regulatory genes. ABC transporter genes reflect the availability of nutrients in the environment, with small genomes presumably relying more on metabolites from other organisms. We also discuss the distinct strategies different phyla adopt to adapt to oligotrophic environments. The findings contribute to our understanding of genomic adaptations within complex microbial communities.

Genomic comparison of deep-sea hydrothermal genera related to Aeropyrum, Thermodiscus and Caldisphaera, and proposed emended description of the family Acidilobaceae

Deep-sea hydrothermal vents host archaeal and bacterial thermophilic communities, including taxonomically and functionally diverse Thermoproteota. Despite their prevalence in high-temperature submarine communities, Thermoproteota are chronically under-represented in genomic databases and issues have emerged regarding their nomenclature, particularly within the Aeropyrum-Thermodiscus-Caldisphaera. To resolve some of these problems, we identified 47 metagenome-assembled genomes (MAGs) within this clade, from 20 previously published deep-sea hydrothermal vent and submarine volcano metagenomes, and 24 MAGs from public databases. Using phylogenomic analysis, Genome Taxonomy Database Toolkit (GTDB-Tk) taxonomic assessment, 16S rRNA gene phylogeny, average amino acid identity (AAI) and functional gene patterns, we re-evaluated of the taxonomy of the Aeropyrum-Thermodiscus-Caldisphaera. At least nine genus-level clades were identified with two or more MAGs. In accordance with SeqCode requirements and recommendations, we propose names for three novel genera, viz. Tiamatella incendiivivens, Hestiella acidicharens and Calypsonella navitae. A fourth genus was also identified related to Thermodiscus maritimus, for which no available sequenced genome exists. We propose the novel species Thermodiscus eudorianus to describe our high-quality Thermodiscus MAG, which represents the type genome for the genus. All three novel genera and T. eudorianus are likely anaerobic heterotrophs, capable of fermenting protein-rich carbon sources, while some Tiamatella, Calypsonella and T. eudorianus may also reduce polysulfides, thiosulfate, sulfur and/or selenite, and the likely acidophile, Hestiella, may reduce nitrate and/or perchlorate. Based on phylogenomic evidence, we also propose the family Acidilobaceae be amended to include Caldisphaera, Aeropyrum, Thermodiscus and Stetteria and the novel genera described here.

The ecological role of microbiome at community-, taxonomic - and genome-levels in black-odorous waters

Black-odorous waters (BOWs) are heavily polluted waters where microbial information remains elusive mechanistically. Based on gene amplicon and metagenomics sequencing, a comprehensive study was conducted to investigate the microbial communities in urban and rural BOWs. The results revealed that microbial communities' assembly in urban and rural BOWs was predominantly governed by stochastic factors at the community level. At the taxonomic level, there were 62 core species (58.48%) in water and 207 core species (44.56%) in sediment across urban and rural areas. Notably, significant differences were observed in the functional genetic composition of BOWs between urban and rural areas. Specifically, rural areas exhibited an enhanced abundance of genes involved in nitrogen fixation, Fe2+ transport, and sulfate reduction. Conversely, urban areas showed higher abundances of some genes associated with carbon fixation, nitrification and denitrification. A sulfur-centered ecological model of microbial communities was constructed by integrating data from the three levels of analysis, and 14 near-complete draft genomes were generated, representing a substantial portion of the microbial community (35.04% in rural BOWs and 29.97% in urban BOWs). This research provides significant insights into the sustainable management and preservation of aquatic ecosystems affected by BOWs.

Krumholzibacteriota and Deltaproteobacteria contain rare genetic potential to liberate carbon from monoaromatic compounds in subsurface coal seams

Biogenic methane in subsurface coal seam environments is produced by diverse consortia of microbes. Although this methane is useful for global energy security, it remains unclear which microbes can liberate carbon from the coal. Most of this carbon is relatively resistant to biodegradation, as it is contained within aromatic rings. Thus, to explore for coal-degrading taxa in the subsurface, this study reconstructed relevant metagenome-assembled genomes (MAGs) from coal seams by using a key genomic marker for the anaerobic degradation of monoaromatic compounds as a guide: the benzoyl-CoA reductase gene (bcrABCD). Three MAGs were identified with this genetic potential. The first represented a novel taxon from the Krumholzibacteriota phylum, which this study is the first to describe. This Krumholzibacteriota MAG contained a full set of genes for benzoyl-CoA dearomatization, in addition to other genes for anaerobic catabolism of monoaromatics. Analysis of Krumholzibacteriota MAGs from other environments revealed that this genetic potential may be common, and thus, Krumholzibacteriota may be important organisms for the liberation of recalcitrant carbon in a broad range of environments. Moreover, the assembly and characterization of two Syntrophorhabdus aromaticivorans MAGs from different continents and a Syntrophaceae sp. MAG implicate the Deltaproteobacteria class in coal seam monoaromatic degradation. Each of these taxa are potential rate-limiting organisms for subsurface coal-to-methane biodegradation. Their description here provides some understanding of their function within the coal seam microbiome and will help inform future efforts in coal bed methane stimulation, anoxic bioremediation of organic pollutants, and assessments of anoxic, subsurface carbon cycling and emissions.IMPORTANCESubsurface coal seams are highly anoxic, oligotrophic environments, where the main source of carbon is "locked away" within aromatic rings. Despite these challenges, many coal seams accumulate biogenic methane, implying that the coal seam microbiome is "unlocking" this carbon source in situ. For over two decades, researchers have endeavored to understand which organisms perform these processes. This study provides the first descriptions of organisms with this genetic potential from the coal seam environment. Here, we report metagenomic insights into carbon liberation from aromatic molecules and the degradation pathways involved and describe a Krumholzibacteriota, two Syntrophorhabdus aromaticivorans, and a Syntrophaceae MAG that contain this genetic potential. This is also the first time that the Krumholzibacteriota phylum has been implicated in anaerobic dearomatization of aromatic hydrocarbons. This potential is identified here in numerous MAGs from other terrestrial and marine subsurface habitats, implicating the Krumholzibacteriota in carbon-cycling processes across a broad range of environments.

Shifts in the Microbial Populations of Bioleach Reactors Are Determined by Carbon Sources and Temperature

In the present study, the effect of additional carbon sources (carbon dioxide and molasses) on the bio-oxidation of a pyrite-arsenopyrite concentrate at temperatures of 40-50 °C was studied, and novel data regarding the patterns of the bio-oxidation of gold-bearing sulfide concentrates and the composition of the microbial populations performing these processes were obtained. At 40 °C, additional carbon sources did not affect the bio-oxidation efficiency. At the same time, the application of additional carbon dioxide improved the bio-oxidation performance at temperatures of 45 and 50 °C and made it possible to avoid the inhibition of bio-oxidation due to an increase in the temperature. Therefore, the use of additional carbon dioxide may be proposed to prevent the negative effect of an increase in temperature on the bio-oxidation of sulfide concentrates. 16S rRNA gene profiling revealed archaea of the family Thermoplasmataceae (Acidiplasma, Ferroplasma, Cuniculiplasma, and A-plasma group) and bacteria of the genera Leptospirillum, with Sulfobacillus and Acidithiobacillus among the dominant groups in the community. Temperature influenced the composition of the communities to a greater extent than the additional sources of carbon and the mode of operation of the bioreactor. Elevating the temperature from 40 °C to 50 °C resulted in increases in the shares of Acidiplasma and Sulfobacillus and decreases in the relative abundances of Ferroplasma, Leptospirillum, and Acidithiobacillus, while Cuniculiplasma and A-plasma were more abundant at 45 °C. A metagenomic analysis of the studied population made it possible to characterize novel archaea belonging to an uncultivated, poorly-studied group of Thermoplasmatales which potentially plays an important role in the bio-oxidation process. Based on an analysis of the complete genome, we propose describing the novel species and novel genus as "Candidatus Carboxiplasma ferriphilum" gen. nov., spec. nov.

Metagenomic survey reveals global distribution and evolution of microbial sialic acid catabolism

Sialic acids comprise a varied group of nine-carbon amino sugars found mostly in humans and other higher metazoans, playing major roles in cell interactions with external environments as well as other cells. Microbial sialic acid catabolism (SAC) has long been considered a virulence determinant, and appears to be mainly the purview of pathogenic and commensal bacterial species associated with eukaryotic hosts. Here, we used 2,521 (pre-)assembled metagenomes to evaluate the distribution of SAC in microbial communities from diverse ecosystems and human body parts. Our results demonstrated that microorganisms possessing SAC globally existed in non-host associated environments, although much less frequently than in mammal hosts. We also showed that the ecological significance and taxonomic diversity of microbial SAC have so far been largely underestimated. Phylogenetic analysis revealed a strong signal of horizontal gene transfer among distinct taxa and habitats, and also suggested a specific ecological pressure and a relatively independent evolution history in environmental communities. Our study expanded the known diversity of microbial SAC, and has provided the backbone for further studies on its ecological roles and potential pathogenesis.

A Metagenome from a Steam Vent in Los Azufres Geothermal Field Shows an Abundance of Thermoplasmatales archaea and Bacteria from the Phyla Actinomycetota and Pseudomonadota

Los Azufres National Park is a geothermal field that has a wide number of thermal manifestations; nevertheless, the microbial communities in many of these environments remain unknown. In this study, a metagenome from a sediment sample from Los Azufres National Park was sequenced. In this metagenome, we found that the microbial diversity corresponds to bacteria (Actinomycetota, Pseudomonadota), archaea (Thermoplasmatales and Candidatus Micrarchaeota and Candidatus Parvarchaeota), eukarya (Cyanidiaceae), and viruses (Fussellovirus and Caudoviricetes). The functional annotation showed genes related to the carbon fixation pathway, sulfur metabolism, genes involved in heat and cold shock, and heavy-metal resistance. From the sediment, it was possible to recover two metagenome-assembled genomes from Ferrimicrobium and Cuniculiplasma. Our results showed that there are a large number of microorganisms in Los Azufres that deserve to be studied.

Evidence for nontraditional mcr-containing archaea contributing to biological methanogenesis in geothermal springs

Recent discoveries of methyl-coenzyme M reductase-encoding genes (mcr) in uncultured archaea beyond traditional euryarchaeotal methanogens have reshaped our view of methanogenesis. However, whether any of these nontraditional archaea perform methanogenesis remains elusive. Here, we report field and microcosm experiments based on ¹³C-tracer labeling and genome-resolved metagenomics and metatranscriptomics, revealing that nontraditional archaea are predominant active methane producers in two geothermal springs. Archaeoglobales performed methanogenesis from methanol and may exhibit adaptability in using methylotrophic and hydrogenotrophic pathways based on temperature/substrate availability. A five-year field survey found Candidatus Nezhaarchaeota to be the predominant mcr-containing archaea inhabiting the springs; genomic inference and mcr expression under methanogenic conditions strongly suggested that this lineage mediated hydrogenotrophic methanogenesis in situ. Methanogenesis was temperature-sensitive , with a preference for methylotrophic over hydrogenotrophic pathways when incubation temperatures increased from 65° to 75°C. This study demonstrates an anoxic ecosystem wherein methanogenesis is primarily driven by archaea beyond known methanogens, highlighting diverse nontraditional mcr-containing archaea as previously unrecognized methane sources.

uBin: A manual refining tool for genomes from metagenomes

Resolving bacterial and archaeal genomes from metagenomes has revolutionized our understanding of Earth's biomes yet producing high-quality genomes from assembled fragments has been an ever-standing problem. While automated binning software and their combination produce prokaryotic bins in high throughput, their manual refinement has been slow, sometimes difficult or missing entirely facilitating error propagation in public databases. Here, we present uBin, a GUI-based, standalone bin refiner that runs on all major operating platforms and was additionally designed for educational purposes. When applied to the public CAMI dataset, refinement of bins using GC content, coverage and taxonomy was able to improve 78.9% of bins by decreasing their contamination. We also applied the bin refiner as a standalone binner to public metagenomes from the International Space Station and demonstrate the recovery of near-complete genomes, whose replication indices indicate the active proliferation of microbes in Earth's lower orbit. uBin is an easy to instal software for bin refinement, binning of simple metagenomes and communication of metagenomic results to other scientists and in classrooms. The software and its helper scripts are open source and available under https://github.com/ProbstLab/uBin.

Inference and reconstruction of the heimdallarchaeial ancestry of eukaryotes

In the ongoing debates about eukaryogenesis-the series of evolutionary events leading to the emergence of the eukaryotic cell from prokaryotic ancestors-members of the Asgard archaea play a key part as the closest archaeal relatives of eukaryotes1. However, the nature and phylogenetic identity of the last common ancestor of Asgard archaea and eukaryotes remain unresolved2-4. Here we analyse distinct phylogenetic marker datasets of an expanded genomic sampling of Asgard archaea and evaluate competing evolutionary scenarios using state-of-the-art phylogenomic approaches. We find that eukaryotes are placed, with high confidence, as a well-nested clade within Asgard archaea and as a sister lineage to Hodarchaeales, a newly proposed order within Heimdallarchaeia. Using sophisticated gene tree and species tree reconciliation approaches, we show that analogous to the evolution of eukaryotic genomes, genome evolution in Asgard archaea involved significantly more gene duplication and fewer gene loss events compared with other archaea. Finally, we infer that the last common ancestor of Asgard archaea was probably a thermophilic chemolithotroph and that the lineage from which eukaryotes evolved adapted to mesophilic conditions and acquired the genetic potential to support a heterotrophic lifestyle. Our work provides key insights into the prokaryote-to-eukaryote transition and a platform for better understanding the emergence of cellular complexity in eukaryotic cells.

Baseline metagenome-assembled genome (MAG) data of Sikkim hot springs from Indian Himalayan geothermal belt (IHGB) showcasing its potential CAZymes, and sulfur-nitrogen metabolic activity

Here we present the construction and characterization of metagenome assembled genomes (MAGs) from two hot springs residing in the vicinity of Indian Himalayan Geothermal Belt (IHGB). A total of 78 and 7 taxonomic bins were obtained for Old Yume Samdong (OYS) and New Yume Samdong (NYS) hot springs respectively. After passing all the criteria only 21 and 4 MAGs were further studied based on the successful prediction of their 16 S rRNA. Various databases were used such as GTDB, Kaiju, EzTaxon, BLAST XY Plot and NCBI BLAST to get the taxonomic classification of various 16 S rRNA predicted MAGs. The bacterial genomes found were from both thermophilic and mesophilic bacteria among which Proteobacteria, Chloroflexi, Bacteroidetes and Firmicutes were the abundant phyla. However, in case of OYS, two genomes belonged to archaeal Methanobacterium and Methanocaldococcus. Functional characterization revealed the richness of CAZymes such as Glycosyl Transferase (GT) (56.7%), Glycoside Hydrolase (GH) (37.4%), Carbohydrate Esterase family (CE) (8.2%), and Polysaccharide Lyase (PL) (1.9%). There were negligible antibiotic resistance genes in the MAGs however, a significant heavy metal tolerance gene was found in the MAGs. Thus, it may be assumed that there is no coexistence of antibiotic and heavy metal resistance genes in these hot spring microbiomes. Since the selected hot springs possess good sulfur content thus, we also checked the presence of genes for sulfur and nitrogen metabolism. It was found that MAGs from both the hot springs possess significant number of genes related to sulfur and nitrogen metabolism.

Anaerobic degradation of organic carbon supports uncultured microbial populations in estuarine sediments

BACKGROUND

A large proportion of prokaryotic microbes in marine sediments remains uncultured, hindering our understanding of their ecological functions and metabolic features. Recent environmental metagenomic studies suggested that many of these uncultured microbes contribute to the degradation of organic matter, accompanied by acetogenesis, but the supporting experimental evidence is limited.

RESULTS

Estuarine sediments were incubated with different types of organic matters under anaerobic conditions, and the increase of uncultured bacterial populations was monitored. We found that (1) lignin stimulated the increase of uncultured bacteria within the class Dehalococcoidia. Their ability to metabolize lignin was further supported by the presence of genes associated with a nearly complete degradation pathway of phenolic monomers in the Dehalococcoidia metagenome-assembled genomes (MAGs). (2) The addition of cellulose stimulated the increase of bacteria in the phylum Ca. Fermentibacterota and family Fibrobacterales, a high copy number of genes encoding extracellular endoglucanase or/and 1,4-beta-cellobiosidase for cellulose decomposition and multiple sugar transporters were present in their MAGs. (3) Uncultured lineages in the order Bacteroidales and the family Leptospiraceae were enriched by the addition of casein and oleic acid, respectively, a high copy number of genes encoding extracellular peptidases, and the complete β-oxidation pathway were found in those MAGs of Bacteroidales and Leptospiraceae, respectively. (4) The growth of unclassified bacteria of the order Clostridiales was found after the addition of both casein and cellulose. Their MAGs contained multiple copies of genes for extracellular peptidases and endoglucanase. Additionally, ¹³C-labeled acetate was produced in the incubations when ¹³C-labeled dissolved inorganic carbon was provided.

CONCLUSIONS

Our results provide new insights into the roles of microorganisms during organic carbon degradation in anaerobic estuarine sediments and suggest that these macro and single molecular organic carbons support the persistence and increase of uncultivated bacteria. Acetogenesis is an additional important microbial process alongside organic carbon degradation. Video Abstract.

Hot spring distribution and survival mechanisms of thermophilic comammox Nitrospira

The recent discovery of Nitrospira species capable of complete ammonia oxidation (comammox) in non-marine natural and engineered ecosystems under mesothermal conditions has changed our understanding of microbial nitrification. However, little is known about the occurrence of comammox bacteria or their ability to survive in moderately thermal and/or hyperthermal habitats. Here, we report the wide distribution of comammox Nitrospira in five terrestrial hot springs at temperatures ranging from 36 to 80°C and provide metagenome-assembled genomes of 11 new comammox strains. Interestingly, the identification of dissimilatory nitrate reduction to ammonium (DNRA) in thermophilic comammox Nitrospira lineages suggests that they have versatile ecological functions as both sinks and sources of ammonia, in contrast to the described mesophilic comammox lineages, which lack the DNRA pathway. Furthermore, the in situ expression of key genes associated with nitrogen metabolism, thermal adaptation, and oxidative stress confirmed their ability to survive in the studied hot springs and their contribution to nitrification in these environments. Additionally, the smaller genome size and higher GC content, less polar and more charged amino acids in usage profiles, and the expression of a large number of heat shock proteins compared to mesophilic comammox strains presumably confer tolerance to thermal stress. These novel insights into the occurrence, metabolic activity, and adaptation of comammox Nitrospira in thermal habitats further expand our understanding of the global distribution of comammox Nitrospira and have significant implications for how these unique microorganisms have evolved thermal tolerance strategies.

Previously uncharacterized rectangular bacterial structures in the dolphin mouth

Much remains to be explored regarding the diversity of uncultured, host-associated microbes. Here, we describe rectangular bacterial structures (RBSs) in the mouths of bottlenose dolphins. DNA staining revealed multiple paired bands within RBSs, suggesting the presence of cells dividing along the longitudinal axis. Cryogenic transmission electron microscopy and tomography showed parallel membrane-bound segments that are likely cells, encapsulated by an S-layer-like periodic surface covering. RBSs displayed unusual pilus-like appendages with bundles of threads splayed at the tips. We present multiple lines of evidence, including genomic DNA sequencing of micromanipulated RBSs, 16S rRNA gene sequencing, and fluorescence in situ hybridization, suggesting that RBSs are bacterial and distinct from the genera Simonsiella and Conchiformibius (family Neisseriaceae), with which they share similar morphology and division patterning. Our findings highlight the diversity of novel microbial forms and lifestyles that await characterization using tools complementary to genomics such as microscopy.

Back to Basics: A Simplified Improvement to Multiple Displacement Amplification for Microbial Single-Cell Genomics

Microbial single-cell genomics (SCG) provides access to the genomes of rare and uncultured microorganisms and is a complementary method to metagenomics. Due to the femtogram-levels of DNA in a single microbial cell, sequencing the genome requires whole genome amplification (WGA) as a preliminary step. However, the most common WGA method, multiple displacement amplification (MDA), is known to be costly and biased against specific genomic regions, preventing high-throughput applications and resulting in uneven genome coverage. Thus, obtaining high-quality genomes from many taxa, especially minority members of microbial communities, becomes difficult. Here, we present a volume reduction approach that significantly reduces costs while improving genome coverage and uniformity of DNA amplification products in standard 384-well plates. Our results demonstrate that further volume reduction in specialized and complex setups (e.g., microfluidic chips) is likely unnecessary to obtain higher-quality microbial genomes. This volume reduction method makes SCG more feasible for future studies, thus helping to broaden our knowledge on the diversity and function of understudied and uncharacterized microorganisms in the environment.

Genomic Signature in Evolutionary Biology: A Review

Organisms are unique physical entities in which information is stored and continuously processed. The digital nature of DNA sequences enables the construction of a dynamic information reservoir. However, the distinction between the hardware and software components in the information flow is crucial to identify the mechanisms generating specific genomic signatures. In this work, we perform a bibliometric analysis to identify the different purposes of looking for particular patterns in DNA sequences associated with a given phenotype. This study has enabled us to make a conceptual breakdown of the genomic signature and differentiate the leading applications. On the one hand, it refers to gene expression profiling associated with a biological function, which may be shared across taxa. This signature is the focus of study in precision medicine. On the other hand, it also refers to characteristic patterns in species-specific DNA sequences. This interpretation plays a key role in comparative genomics, identifying evolutionary relationships. Looking at the relevant studies in our bibliographic database, we highlight the main factors causing heterogeneities in genome composition and how they can be quantified. All these findings lead us to reformulate some questions relevant to evolutionary biology.

MetaBinner: a high-performance and stand-alone ensemble binning method to recover individual genomes from complex microbial communities

Binning aims to recover microbial genomes from metagenomic data. For complex metagenomic communities, the available binning methods are far from satisfactory, which usually do not fully use different types of features and important biological knowledge. We developed a novel ensemble binner, MetaBinner, which generates component results with multiple types of features by k-means and uses single-copy gene information for initialization. It then employs a two-stage ensemble strategy based on single-copy genes to integrate the component results efficiently and effectively. Extensive experimental results on three large-scale simulated datasets and one real-world dataset demonstrate that MetaBinner outperforms the state-of-the-art binners significantly.

Uncovering Lasonolide A Biosynthesis Using Genome-Resolved Metagenomics

Invertebrates, particularly sponges, have been a dominant source of new marine natural products. For example, lasonolide A (LSA) is a potential anticancer molecule isolated from the marine sponge Forcepia sp., with nanomolar growth inhibitory activity and a unique cytotoxicity profile against the National Cancer Institute 60-cell-line screen. Here, we identified the putative biosynthetic pathway for LSA. Genomic binning of the Forcepia sponge metagenome revealed a Gram-negative bacterium belonging to the phylum Verrucomicrobia as the candidate producer of LSA. Phylogenetic analysis showed that this bacterium, here named "Candidatus Thermopylae lasonolidus," only has 88.78% 16S rRNA identity with the closest relative, Pedosphaera parvula Ellin514, indicating that it represents a new genus. The lasonolide A (las) biosynthetic gene cluster (BGC) was identified as a trans-acyltransferase (AT) polyketide synthase (PKS) pathway. Compared with its host genome, the las BGC exhibits a significantly different GC content and pentanucleotide frequency, suggesting a potential horizontal acquisition of the gene cluster. Furthermore, three copies of the putative las pathway were identified in the candidate producer genome. Differences between the three las repeats were observed, including the presence of three insertions, two single-nucleotide polymorphisms, and the absence of a stand-alone acyl carrier protein in one of the repeats. Even though the verrucomicrobial producer shows signs of genome reduction, its genome size is still fairly large (about 5 Mbp), and, compared to its closest free-living relative, it contains most of the primary metabolic pathways, suggesting that it is in the early stages of reduction. IMPORTANCE While sponges are valuable sources of bioactive natural products, a majority of these compounds are produced in small quantities by uncultured symbionts, hampering the study and clinical development of these unique compounds. Lasonolide A (LSA), isolated from marine sponge Forcepia sp., is a cytotoxic molecule active at nanomolar concentrations, which causes premature chromosome condensation, blebbing, cell contraction, and loss of cell adhesion, indicating a novel mechanism of action and making it a potential anticancer drug lead. However, its limited supply hampers progression to clinical trials. We investigated the microbiome of Forcepia sp. using culture-independent DNA sequencing, identified genes likely responsible for LSA synthesis in an uncultured bacterium, and assembled the symbiont's genome. These insights provide future opportunities for heterologous expression and cultivation efforts that may minimize LSA's supply problem.

Insights into the microbial response mechanisms to ciprofloxacin during sulfur-mediated biological wastewater treatment using a metagenomics approach

The fate and removal of ciprofloxacin, a class of fluoroquinolone antibiotic, during sulfur-mediated biological wastewater treatment has been recently well documented. However, little is known regarding the genetic response of microorganisms to ciprofloxacin. Here, a lab-scale anaerobic sulfate-reducing bioreactor was continuously operated over a long term for ciprofloxacin-contaminated wastewater treatment to investigate the response of the microorganisms to ciprofloxacin by adopting a metagenomics approach. It was found that total organic carbon (TOC) removal and sulfate reduction were promoted by approximately 10% under ciprofloxacin stress, along with the enrichment of functional genera (e.g., Desulfobacter, Geobacter) involved in carbon and sulfur metabolism. The metagenomic analytical results demonstrated that ciprofloxacin triggered the microbial SOS response, as demonstrated by the up-regulation of the multidrug efflux pump genes (8-125-fold higher than that of the control) and ciprofloxacin-degrading genes (4-33-fold higher than that of the control). Moreover, the contents of ATP, NADH, and cytochrome C, as well as related functional genes (including genes involved in energy generation, electron transport, carbon metabolism, and sulfur metabolism) were markedly increased under ciprofloxacin stress. This demonstrated that the carbon and sulfur metabolisms were enhanced for energy (ATP) generation and electron transport in response to ciprofloxacin-induced stress. Interestingly, the microbes tended to cooperate while being subjected exposure to exogenous ciprofloxacin according to the reconstructed metabolic network using the NetSeed model. Particularly, the species with higher complementarity indices played more pivotal roles in strengthening microbial metabolism and the SOS response under long-term ciprofloxacin stress. This study characterized the response mechanisms of microorganisms to ciprofloxacin at the genetic level in sulfur-mediated biological wastewater treatment. These new understandings will contribute the scientific basis for improving and optimizing the sulfur-mediated bioprocess for antibiotics-laden wastewater treatment.

Metagenomic methylation patterns resolve bacterial genomes of unusual size and structural complexity

The plasticity of bacterial and archaeal genomes makes examining their ecological and evolutionary dynamics both exciting and challenging. The same mechanisms that enable rapid genomic change and adaptation confound current approaches for recovering complete genomes from metagenomes. Here, we use strain-specific patterns of DNA methylation to resolve complex bacterial genomes from long-read metagenomic data of a marine microbial consortium, the “pink berries” of the Sippewissett Marsh (USA). Unique combinations of restriction-modification (RM) systems encoded by the bacteria produced distinctive methylation profiles that were used to accurately bin and classify metagenomic sequences. Using this approach, we finished the largest and most complex circularized bacterial genome ever recovered from a metagenome (7.9 Mb with >600 transposons), the finished genome of Thiohalocapsa sp. PB-PSB1 the dominant bacteria in the consortia. From genomes binned by methylation patterns, we identified instances of horizontal gene transfer between sulfur-cycling symbionts (Thiohalocapsa sp. PB-PSB1 and Desulfofustis sp. PB-SRB1), phage infection, and strain-level structural variation. We also linked the methylation patterns of each metagenome-assembled genome with encoded DNA methyltransferases and discovered new RM defense systems, including novel associations of RM systems with RNase toxins.

Heterotrophic Bacteria Dominate Catalase Expression during Microcystis Blooms

In the oligotrophic oceans, key autotrophs depend on "helper" bacteria to reduce oxidative stress from hydrogen peroxide (H₂O₂) in the extracellular environment. H₂O₂ is also a ubiquitous stressor in freshwaters, but the effects of H₂O₂ on autotrophs and their interactions with bacteria are less well understood in freshwaters. Naturally occurring H₂O₂ in freshwater systems is proposed to impact the proportion of microcystin-producing (toxic) and non-microcystin-producing (nontoxic) Microcystis in blooms, which influences toxin concentrations and human health impacts. However, how different strains of Microcystis respond to naturally occurring H₂O₂ concentrations and the microbes responsible for H₂O₂ decomposition in freshwater cyanobacterial blooms are unknown. To address these knowledge gaps, we used metagenomics and metatranscriptomics to track the presence and expression of genes for H₂O₂ decomposition by microbes during a cyanobacterial bloom in western Lake Erie in the summer of 2014. katG encodes the key enzyme for decomposing extracellular H₂O₂ but was absent in most Microcystis cells. katG transcript relative abundance was dominated by heterotrophic bacteria. In axenic Microcystis cultures, an H₂O₂ scavenger (pyruvate) significantly improved growth rates of one toxic strain while other toxic and nontoxic strains were unaffected. These results indicate that heterotrophic bacteria play a key role in H₂O₂ decomposition in Microcystis blooms and suggest that their activity may affect the fitness of some Microcystis strains and thus the strain composition of Microcystis blooms but not along a toxic versus nontoxic dichotomy. IMPORTANCE Cyanobacterial harmful algal blooms (CHABs) threaten freshwater ecosystems globally through the production of toxins. Toxin production by cyanobacterial species and strains during CHABs varies widely over time and space, but the ecological drivers of the succession of toxin-producing species remain unclear. Hydrogen peroxide (H₂O₂) is ubiquitous in natural waters, inhibits microbial growth, and may determine the relative proportions of Microcystis strains during blooms. However, the mechanisms and organismal interactions involved in H₂O₂ decomposition are unexplored in CHABs. This study shows that some strains of bloom-forming freshwater cyanobacteria benefit from detoxification of H₂O₂ by associated heterotrophic bacteria, which may impact bloom development.

Uptake of Phytoplankton-Derived Carbon and Cobalamins by Novel Acidobacteria Genera in Microcystis Blooms Inferred from Metagenomic and Metatranscriptomic Evidence

Interactions between bacteria and phytoplankton can influence primary production, community composition, and algal bloom development. However, these interactions are poorly described for many consortia, particularly for freshwater bloom-forming cyanobacteria. Here, we assessed the gene content and expression of two uncultivated Acidobacteria from Lake Erie Microcystis blooms. These organisms were targeted because they were previously identified as important catalase producers in Microcystis blooms, suggesting that they protect Microcystis from H₂O₂. Metatranscriptomics revealed that both Acidobacteria transcribed genes for uptake of organic compounds that are known cyanobacterial products and exudates, including lactate, glycolate, amino acids, peptides, and cobalamins. Expressed genes for amino acid metabolism and peptide transport and degradation suggest that use of amino acids and peptides by Acidobacteria may regenerate nitrogen for cyanobacteria and other organisms. The Acidobacteria genomes lacked genes for biosynthesis of cobalamins but expressed genes for its transport and remodeling. This indicates that the Acidobacteria obtained cobalamins externally, potentially from Microcystis, which has a complete gene repertoire for pseudocobalamin biosynthesis; expressed them in field samples; and produced pseudocobalamin in axenic culture. Both Acidobacteria were detected in Microcystis blooms worldwide. Together, the data support the hypotheses that uncultured and previously unidentified Acidobacteria taxa exchange metabolites with phytoplankton during harmful cyanobacterial blooms and influence nitrogen available to phytoplankton. Thus, novel Acidobacteria may play a role in cyanobacterial physiology and bloom development. IMPORTANCE Interactions between heterotrophic bacteria and phytoplankton influence competition and successions between phytoplankton taxa, thereby influencing ecosystem-wide processes such as carbon cycling and algal bloom development. The cyanobacterium Microcystis forms harmful blooms in freshwaters worldwide and grows in buoyant colonies that harbor other bacteria in their phycospheres. Bacteria in the phycosphere and in the surrounding community likely influence Microcystis physiology and ecology and thus the development of freshwater harmful cyanobacterial blooms. However, the impacts and mechanisms of interaction between bacteria and Microcystis are not fully understood. This study explores the mechanisms of interaction between Microcystis and uncultured members of its phycosphere in situ with population genome resolution to investigate the cooccurrence of Microcystis and freshwater Acidobacteria in blooms worldwide.

Conserved and lineage-specific hypothetical proteins may have played a central role in the rise and diversification of major archaeal groups

BACKGROUND

Archaea play fundamental roles in the environment, for example by methane production and consumption, ammonia oxidation, protein degradation, carbon compound turnover, and sulfur compound transformations. Recent genomic analyses have profoundly reshaped our understanding of the distribution and functionalities of Archaea and their roles in eukaryotic evolution.

RESULTS

Here, 1179 representative genomes were selected from 3197 archaeal genomes. The representative genomes clustered based on the content of 10,866 newly defined archaeal protein families (that will serve as a community resource) recapitulates archaeal phylogeny. We identified the co-occurring proteins that distinguish the major lineages. Those with metabolic roles were consistent with experimental data. However, two families specific to Asgard were determined to be new eukaryotic signature proteins. Overall, the blocks of lineage-specific families are dominated by proteins that lack functional predictions.

CONCLUSIONS

Given that these hypothetical proteins are near ubiquitous within major archaeal groups, we propose that they were important in the origin of most of the major archaeal lineages. Interestingly, although there were clearly phylum-specific co-occurring proteins, no such blocks of protein families were shared across superphyla, suggesting a burst-like origin of new lineages early in archaeal evolution.

Potential metabolic and genetic interaction among viruses, methanogen and methanotrophic archaea, and their syntrophic partners

The metabolism of methane in anoxic ecosystems is mainly mediated by methanogens and methane-oxidizing archaea (MMA), key players in global carbon cycling. Viruses are vital in regulating their host fate and ecological function. However, our knowledge about the distribution and diversity of MMA viruses and their interactions with hosts is rather limited. Here, by searching metagenomes containing mcrA (the gene coding for the α-subunit of methyl-coenzyme M reductase) from a wide variety of environments, 140 viral operational taxonomic units (vOTUs) that potentially infect methanogens or methane-oxidizing archaea were retrieved. Four MMA vOTUs (three infecting the order Methanobacteriales and one infecting the order Methanococcales) were predicted to cross-domain infect sulfate-reducing bacteria. By facilitating assimilatory sulfur reduction, MMA viruses may increase the fitness of their hosts in sulfate-depleted anoxic ecosystems and benefit from synthesis of the sulfur-containing amino acid cysteine. Moreover, cell-cell aggregation promoted by MMA viruses may be beneficial for both the viruses and their hosts by improving infectivity and environmental stress resistance, respectively. Our results suggest a potential role of viruses in the ecological and environmental adaptation of methanogens and methane-oxidizing archaea.

FRAGTE2: An Enhanced Algorithm to Pre-Select Closely Related Genomes for Bacterial Species Demarcation

We previously reported on FRAGTE (hereafter termed FRAGTE1), a promising algorithm for sieving (pre-selecting genome pairs for whole-genome species demarcation). However, the overall amount of pairs sieved by FRAGTE1 is still large, requiring seriously unaffordable computing cost, especially for large datasets. Here, we present FRAGTE2. Tests on simulated genomes, real genomes, and metagenome-assembled genomes revealed that (i) FRAGTE2 outstandingly reduces ~50–60.10% of the overall amount of pairs sieved by FRAGTE1, dramatically decreasing the computing cost required for whole-genome species demarcation afterward; (ii) FRAGTE2 shows superior sensitivity than FRAGTE1; (iii) FRAGTE2 shows higher specificity than FRAGTE1; and (iv) FRAGTE2 is faster than or comparable with FRAGTE1. Besides, FRAGTE2 is independent of genome completeness, the same as FRAGTE1. We therefore recommend FRAGTE2 tailored for sieving to facilitate species demarcation in prokaryotes.

Host-Associated Phages Disperse across the Extraterrestrial Analogue Antarctica

Extreme Antarctic conditions provide one of the closest analogues of extraterrestrial environments. Since air and snow samples, especially from polar regions, yield DNA amounts in the lower picogram range, binning of prokaryotic genomes is challenging and renders studying the dispersal of biological entities across these environments difficult. Here, we hypothesized that dispersal of host-associated bacteriophages (adsorbed, replicating, or prophages) across the Antarctic continent can be tracked via their genetic signatures, aiding our understanding of virus and host dispersal across long distances. Phage genome fragments (PGFs) reconstructed from surface snow metagenomes of three Antarctic stations were assigned to four host genomes, mainly Betaproteobacteria, including Ralstonia spp. We reconstructed the complete genome of a temperate phage with nearly complete alignment to a prophage in the reference genome of Ralstonia pickettii 12D. PGFs from different stations were related to each other at the genus level and matched similar hosts. Metagenomic read mapping and nucleotide polymorphism analysis revealed a wide dispersal of highly identical PGFs, 13 of which were detected in seawater from the Western Antarctic Peninsula at a distance of 5,338 km from the snow sampling stations. Our results suggest that host-associated phages, especially of Ralstonia sp., disperse over long distances despite the harsh conditions of the Antarctic continent. Given that 14 phages associated with two R. pickettii draft genomes isolated from space equipment were identified, we conclude that Ralstonia phages are ideal mobile genetic elements to track dispersal and contamination in ecosystems relevant for astrobiology. IMPORTANCE Host-associated phages of the bacterium Ralstonia identified in snow samples can be used to track microbial dispersal over thousands of kilometers across the Antarctic continent, which functions as an extraterrestrial analogue because of its harsh environmental conditions. Due to the presence of these bacteria carrying genome-integrated prophages on space-related equipment and the potential for dispersal of host-associated phages demonstrated here, our work has implications for planetary protection, a discipline in astrobiology interested in preventing contamination of celestial bodies with alien biomolecules or forms of life.

Dataset for Genome Sequencing and De Novo Assembly of the Candidate Phyla Radiation in Supragingival Plaque

The Candidate Phyla Radiation (CPR), as a newly discovered and difficult-to-culture bacterium, accounts for the majority of the bacterial domain, which may be related to various oral diseases, including dental caries. Restricted by laboratory culture conditions, there is limited knowledge about oral CPR. Advances in metagenomics provide a new way to study CPR through molecular biology. Here, we used metagenomic assembly and binning to reconstruct more and higher quality metagenome-assembled genomes (MAGs) of CPR from oral dental plaque. These MAGs represent novel CPR species, which differed from all known CPR organisms. Relative abundance of different CPR MAGs in the caries and caries-free group was estimated by mapping metagenomic reads to newly constructed MAGs. The relative abundance of two CPR MAGs was significantly increased in the caries group, indicating that there might be a relationship with caries activity. The detection of a large number of unclassified CPR MAGs in the dataset implies that the phylogenetic diversity of CPR is enormous. The results provide a reference value for exploring the ecological distribution and function of uncultured or difficult-to-culture microorganisms.

Anaerobic derivates of mitochondria and peroxisomes in the free-living amoeba Pelomyxa schiedti revealed by single-cell genomics

Background

Mitochondria and peroxisomes are the two organelles that are most affected during adaptation to microoxic or anoxic environments. Mitochondria are known to transform into anaerobic mitochondria, hydrogenosomes, mitosomes, and various transition stages in between, collectively called mitochondrion-related organelles (MROs), which vary in enzymatic capacity. Anaerobic peroxisomes were identified only recently, and their putatively most conserved function seems to be the metabolism of inositol. The group Archamoebae includes anaerobes bearing both anaerobic peroxisomes and MROs, specifically hydrogenosomes in free-living Mastigamoeba balamuthi and mitosomes in the human pathogen Entamoeba histolytica, while the organelles within the third lineage represented by Pelomyxa remain uncharacterized.

Results

We generated high-quality genome and transcriptome drafts from Pelomyxa schiedti using single-cell omics. These data provided clear evidence for anaerobic derivates of mitochondria and peroxisomes in this species, and corresponding vesicles were tentatively identified in electron micrographs. In silico reconstructed MRO metabolism harbors respiratory complex II, electron-transferring flavoprotein, a partial TCA cycle running presumably in the reductive direction, pyruvate:ferredoxin oxidoreductase, [FeFe]-hydrogenases, a glycine cleavage system, a sulfate activation pathway, and an expanded set of NIF enzymes for iron-sulfur cluster assembly. When expressed in the heterologous system of yeast, some of these candidates localized into mitochondria, supporting their involvement in the MRO metabolism. The putative functions of P. schiedti peroxisomes could be pyridoxal 5′-phosphate biosynthesis, amino acid and carbohydrate metabolism, and hydrolase activities. Unexpectedly, out of 67 predicted peroxisomal enzymes, only four were also reported in M. balamuthi, namely peroxisomal processing peptidase, nudix hydrolase, inositol 2-dehydrogenase, and d-lactate dehydrogenase. Localizations in yeast corroborated peroxisomal functions of the latter two.

Conclusions

This study revealed the presence and partially annotated the function of anaerobic derivates of mitochondria and peroxisomes in P. schiedti using single-cell genomics, localizations in yeast heterologous systems, and transmission electron microscopy. The MRO metabolism resembles that of M. balamuthi and most likely reflects the state in the common ancestor of Archamoebae. The peroxisomal metabolism is strikingly richer in P. schiedti. The presence of myo-inositol 2-dehydrogenase in the predicted peroxisomal proteome corroborates the situation in other Archamoebae, but future experimental evidence is needed to verify additional functions of this organelle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01247-w.

Genetic diversity in terrestrial subsurface ecosystems impacted by geological degassing

Earth’s mantle releases 38.7 ± 2.9 Tg/yr CO₂ along with other reduced and oxidized gases to the atmosphere shaping microbial metabolism at volcanic sites across the globe, yet little is known about its impact on microbial life under non-thermal conditions. Here, we perform comparative metagenomics coupled to geochemical measurements of deep subsurface fluids from a cold-water geyser driven by mantle degassing. Key organisms belonging to uncultivated Candidatus Altiarchaeum show a global biogeographic pattern and site-specific adaptations shaped by gene loss and inter-kingdom horizontal gene transfer. Comparison of the geyser community to 16 other publicly available deep subsurface sites demonstrate a conservation of chemolithoautotrophic metabolism across sites. In silico replication measures suggest a linear relationship of bacterial replication with ecosystems depth with the exception of impacted sites, which show near surface characteristics. Our results suggest that subsurface ecosystems affected by geological degassing are hotspots for microbial life in the deep biosphere.

Geological degassing can impact subsurface metabolism. Here, the authors describe microbial communities from a cold-water geyser are described and compared with other deep subsurface sites, finding a key role for an uncultivated archaeon.

Reconstruction of Archaeal Genomes from Short-Read Metagenomes

As the majority of biological diversity remains unexplored and uncultured, investigating it requires culture-independent approaches. Archaea in particular suffer from a multitude of issues that make their culturing problematic, from them being frequently members of the rare biosphere, to low growth rates, to them thriving under very specific and often extreme environmental and community conditions that are difficult to replicate. OMICs techniques are state of the art approaches that allow direct high-throughput investigations of environmental samples at all levels from nucleic acids to proteins, lipids, and secondary metabolites. Metagenomics, as the foundation for other OMICs techniques, facilitates the identification and functional characterization of the microbial community members and can be combined with other methods to provide insights into the microbial activities, both on the RNA and protein levels. In this chapter, we provide a step-by-step workflow for the recovery of archaeal genomes from metagenomes, starting from raw short-read sequences. This workflow can be applied to recover bacterial genomes as well.

High quality genome assembly of the amitochondriate eukaryote Monocercomonoides exilis

Monocercomonoides exilis is considered the first known eukaryote to completely lack mitochondria. This conclusion is based primarily on a genomic and transcriptomic study which failed to identify any mitochondrial hallmark proteins. However, the available genome assembly has limited contiguity and around 1.5 % of the genome sequence is represented by unknown bases. To improve the contiguity, we re-sequenced the genome and transcriptome of M. exilis using Oxford Nanopore Technology (ONT). The resulting draft genome is assembled in 101 contigs with an N50 value of 1.38 Mbp, almost 20 times higher than the previously published assembly. Using a newly generated ONT transcriptome, we further improve the gene prediction and add high quality untranslated region (UTR) annotations, in which we identify two putative polyadenylation signals present in the 3′UTR regions and characterise the Kozak sequence in the 5′UTR regions. All these improvements are reflected by higher BUSCO genome completeness values. Regardless of an overall more complete genome assembly without missing bases and a better gene prediction, we still failed to identify any mitochondrial hallmark genes, thus further supporting the hypothesis on the absence of mitochondrion.

The economical lifestyle of CPR bacteria in groundwater allows little preference for environmental drivers

BACKGROUND

The highly diverse Cand. Patescibacteria are predicted to have minimal biosynthetic and metabolic pathways, which hinders understanding of how their populations differentiate in response to environmental drivers or host organisms. Their mechanisms employed to cope with oxidative stress are largely unknown. Here, we utilized genome-resolved metagenomics to investigate the adaptive genome repertoire of Patescibacteria in oxic and anoxic groundwaters, and to infer putative host ranges.

RESULTS

Within six groundwater wells, Cand. Patescibacteria was the most dominant (up to 79%) super-phylum across 32 metagenomes sequenced from DNA retained on 0.2 and 0.1 µm filters after sequential filtration. Of the reconstructed 1275 metagenome-assembled genomes (MAGs), 291 high-quality MAGs were classified as Cand. Patescibacteria. Cand. Paceibacteria and Cand. Microgenomates were enriched exclusively in the 0.1 µm fractions, whereas candidate division ABY1 and Cand. Gracilibacteria were enriched in the 0.2 µm fractions. On average, Patescibacteria enriched in the smaller 0.1 µm filter fractions had 22% smaller genomes, 13.4% lower replication measures, higher proportion of rod-shape determining proteins, and of genomic features suggesting type IV pili mediated cell-cell attachments. Near-surface wells harbored Patescibacteria with higher replication rates than anoxic downstream wells characterized by longer water residence time. Except prevalence of superoxide dismutase genes in Patescibacteria MAGs enriched in oxic groundwaters (83%), no major metabolic or phylogenetic differences were observed. The most abundant Patescibacteria MAG in oxic groundwater encoded a nitrate transporter, nitrite reductase, and F-type ATPase, suggesting an alternative energy conservation mechanism. Patescibacteria consistently co-occurred with one another or with members of phyla Nanoarchaeota, Bacteroidota, Nitrospirota, and Omnitrophota. Among the MAGs enriched in 0.2 µm fractions,, only 8% Patescibacteria showed highly significant one-to-one correlation, mostly with Omnitrophota. Motility and transport related genes in certain Patescibacteria were highly similar to genes from other phyla (Omnitrophota, Proteobacteria and Nanoarchaeota).

CONCLUSION

Other than genes to cope with oxidative stress, we found little genomic evidence for niche adaptation of Patescibacteria to oxic or anoxic groundwaters. Given that we could detect specific host preference only for a few MAGs, we speculate that the majority of Patescibacteria is able to attach multiple hosts just long enough to loot or exchange supplies.

Leave no stone unturned: individually adapted xerotolerant Thaumarchaeota sheltered below the boulders of the Atacama Desert hyperarid core

BACKGROUND

The hyperarid core of the Atacama Desert is an extremely harsh environment thought to be colonized by only a few heterotrophic bacterial species. Current concepts for understanding this extreme ecosystem are mainly based on the diversity of these few species, yet a substantial area of the Atacama Desert hyperarid topsoil is covered by expansive boulder accumulations, whose underlying microbiomes have not been investigated so far. With the hypothesis that these sheltered soils harbor uniquely adapted microbiomes, we compared metagenomes and geochemistry between soils below and beside boulders across three distantly located boulder accumulations in the Atacama Desert hyperarid core.

RESULTS

Genome-resolved metagenomics of eleven samples revealed substantially different microbial communities in soils below and beside boulders, despite the presence of shared species. Archaea were found in significantly higher relative abundance below the boulders across all samples within distances of up to 205 km. These key taxa belong to a novel genus of ammonia-oxidizing Thaumarchaeota, Candidatus Nitrosodeserticola. We resolved eight mid-to-high quality genomes of this genus and used comparative genomics to analyze its pangenome and site-specific adaptations. Ca. Nitrosodeserticola genomes contain genes for ammonia oxidation, the 3-hydroxypropionate/4-hydroxybutyrate carbon fixation pathway, and acetate utilization indicating a chemolithoautotrophic and mixotrophic lifestyle. They also possess the capacity for tolerating extreme environmental conditions as highlighted by the presence of genes against oxidative stress and DNA damage. Site-specific adaptations of the genomes included the presence of additional genes for heavy metal transporters, multiple types of ATP synthases, and divergent genes for aquaporins.

CONCLUSION

We provide the first genomic characterization of hyperarid soil microbiomes below the boulders in the Atacama Desert, and report abundant and highly adapted Thaumarchaeaota with ammonia oxidation and carbon fixation potential. Ca. Nitrosodeserticola genomes provide the first metabolic and physiological insight into a thaumarchaeal lineage found in globally distributed terrestrial habitats characterized by various environmental stresses. We consequently expand not only the known genetic repertoire of Thaumarchaeota but also the diversity and microbiome functioning in hyperarid ecosystems. Video Abstract.

Prospects for multi-omics in the microbial ecology of water engineering

Advances in high-throughput sequencing technologies and bioinformatics approaches over almost the last three decades have substantially increased our ability to explore microorganisms and their functions - including those that have yet to be cultivated in pure isolation. Genome-resolved metagenomic approaches have enabled linking powerful functional predictions to specific taxonomical groups with increasing fidelity. Additionally, related developments in both whole community gene expression surveys and metabolite profiling have permitted for direct surveys of community-scale functions in specific environmental settings. These advances have allowed for a shift in microbiome science away from descriptive studies and towards mechanistic and predictive frameworks for designing and harnessing microbial communities for desired beneficial outcomes. Water engineers, microbiologists, and microbial ecologists studying activated sludge, anaerobic digestion, and drinking water distribution systems have applied various (meta)omics techniques for connecting microbial community dynamics and physiologies to overall process parameters and system performance. However, the rapid pace at which new omics-based approaches are developed can appear daunting to those looking to apply these state-of-the-art practices for the first time. Here, we review how modern genome-resolved metagenomic approaches have been applied to a variety of water engineering applications from lab-scale bioreactors to full-scale systems. We describe integrated omics analysis across engineered water systems and the foundations for pairing these insights with modeling approaches. Lastly, we summarize emerging omics-based technologies that we believe will be powerful tools for water engineering applications. Overall, we provide a framework for microbial ecologists specializing in water engineering to apply cutting-edge omics approaches to their research questions to achieve novel functional insights. Successful adoption of predictive frameworks in engineered water systems could enable more economically and environmentally sustainable bioprocesses as demand for water and energy resources increases.

Large-scale k-mer-based analysis of the informational properties of genomes, comparative genomics and taxonomy

Information theoretic approaches are ubiquitous and effective in a wide variety of bioinformatics applications. In comparative genomics, alignment-free methods, based on short DNA words, or k-mers, are particularly powerful. We evaluated the utility of varying k-mer lengths for genome comparisons by analyzing their sequence space coverage of 5805 genomes in the KEGG GENOME database. In subsequent analyses on four k-mer lengths spanning the relevant range (11, 21, 31, 41), hierarchical clustering of 1634 genus-level representative genomes using pairwise 21- and 31-mer Jaccard similarities best recapitulated a phylogenetic/taxonomic tree of life with clear boundaries for superkingdom domains and high subtree similarity for named taxons at lower levels (family through phylum). By analyzing ~14.2M prokaryotic genome comparisons by their lowest-common-ancestor taxon levels, we detected many potential misclassification errors in a curated database, further demonstrating the need for wide-scale adoption of quantitative taxonomic classifications based on whole-genome similarity.

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.

AI for the collective analysis of a massive number of genome sequences: various examples from the small genome of pandemic SARS-CoV-2 to the human genome

In genetics and related fields, huge amounts of data, such as genome sequences, are accumulating, and the use of artificial intelligence (AI) suitable for big data analysis has become increasingly important. Unsupervised AI that can reveal novel knowledge from big data without prior knowledge or particular models is highly desirable for analyses of genome sequences, particularly for obtaining unexpected insights. We have developed a batch-learning self-organizing map (BLSOM) for oligonucleotide compositions that can reveal various novel genome characteristics. Here, we explain the data mining by the BLSOM: an unsupervised AI. As a specific target, we first selected SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) because a large number of viral genome sequences have been accumulated via worldwide efforts. We analyzed more than 0.6 million sequences collected primarily in the first year of the pandemic. BLSOMs for short oligonucleotides (e.g., 4-6-mers) allowed separation into known clades, but longer oligonucleotides further increased the separation ability and revealed subgrouping within known clades. In the case of 15-mers, there is mostly one copy in the genome; thus, 15-mers that appeared after the epidemic started could be connected to mutations, and the BLSOM for 15-mers revealed the mutations that contributed to separation into known clades and their subgroups. After introducing the detailed methodological strategies, we explain BLSOMs for various topics, such as the tetranucleotide BLSOM for over 5 million 5-kb fragment sequences derived from almost all microorganisms currently available and its use in metagenome studies. We also explain BLSOMs for various eukaryotes, including fishes, frogs and Drosophila species, and found a high separation ability among closely related species. When analyzing the human genome, we found enrichments in transcription factor-binding sequences in centromeric and pericentromeric heterochromatin regions. The tDNAs (tRNA genes) could be separated according to their corresponding amino acid.

Improving Fungal Cultivability for Natural Products Discovery

The pool of fungal secondary metabolites can be extended by activating silent gene clusters of cultured strains or by using sensitive biological assays that detect metabolites missed by analytical methods. Alternatively, or in parallel with the first approach, one can increase the diversity of existing culture collections to improve the access to new natural products. This review focuses on the latter approach of screening previously uncultured fungi for chemodiversity. Both strategies have been practiced since the early days of fungal biodiscovery, yet relatively little has been done to overcome the challenge of cultivability of as-yet-uncultivated fungi. Whereas earlier cultivability studies using media formulations and biological assays to scrutinize fungal growth and associated factors were actively conducted, the application of modern omics methods remains limited to test how to culture the fungal dark matter and recalcitrant groups of described fungi. This review discusses the development of techniques to increase the cultivability of filamentous fungi that include culture media formulations and the utilization of known chemical growth factors, in situ culturing and current synthetic biology approaches that build upon knowledge from sequenced genomes. We list more than 100 growth factors, i.e., molecules, biological or physical factors that have been demonstrated to induce spore germination as well as tens of inducers of mycelial growth. We review culturing conditions that can be successfully manipulated for growth of fungi and visit recent information from omics methods to discuss the metabolic basis of cultivability. Earlier work has demonstrated the power of co-culturing fungi with their host, other microorganisms or their exudates to increase their cultivability. Co-culturing of two or more organisms is also a strategy used today for increasing cultivability. However, fungi possess an increased risk for cross-contaminations between isolates in existing in situ or microfluidics culturing devices. Technological improvements for culturing fungi are discussed in the review. We emphasize that improving the cultivability of fungi remains a relevant strategy in drug discovery and underline the importance of ecological and taxonomic knowledge in culture-dependent drug discovery. Combining traditional and omics techniques such as single cell or metagenome sequencing opens up a new era in the study of growth factors of hundreds of thousands of fungal species with high drug discovery potential.

Desulfovibrio feeding Methanobacterium with electrons in conductive methanogenic aggregates from coastal zones

Geobacter, as a typical electroactive microorganism, is the "engine" of interspecies electron transfer (IET) between microorganisms. However, it does not have a dominant position in all natural environments. It is not known what performs a similar function as Geobacter in coastal zones. Metagenomic and metatranscriptomic analysis revealed that Desulfovibrio and Methanobacterium species were the most abundant in electrochemically active aggregates. Metatranscriptomic analysis showed that Desulfovibrio species highly expressed genes for ethanol metabolism and extracellular electron transfer involving cytochromes, pili and flagella. Methanobacterium species in the aggregates also expressed genes for enzymes involved in reducing carbon dioxide to methane. Pure cultures demonstrated that the isolated Desulfovibrio sp. strain JY contributed to aggregate conductivity and directly transferred electrons to Methanothrix harundinacea, which is unable to use H₂ or formate. Most importantly, further coculture studies indicated that Methanobacterium strain YSL might directly accept electrons from the Desulfovibrio strain JY for the reduction of carbon dioxide to methane in the aggregate. This finding suggested that the possibility of DIET by Desulfovibrio similar to Geobacter species in conductive methanogenic aggregates can not be excluded.

Deciphering Symbiotic Interactions of "Candidatus Aenigmarchaeota" with Inferred Horizontal Gene Transfers and Co-occurrence Networks

"Candidatus Aenigmarchaeota" ("Ca. Aenigmarchaeota") represents one of the earliest proposed evolutionary branches within the Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota, and Nanohaloarchaeota (DPANN) superphylum. However, their ecological roles and potential host-symbiont interactions are still poorly understood. Here, eight metagenome-assembled genomes (MAGs) were reconstructed from hot spring ecosystems, and further in-depth comparative and evolutionary genomic analyses were conducted on these MAGs and other genomes downloaded from public databases. Although with limited metabolic capacities, we reported that "Ca. Aenigmarchaeota" in thermal environments harbor more genes related to carbohydrate metabolism than "Ca. Aenigmarchaeota" in nonthermal environments. Evolutionary analyses suggested that members from the Thaumarchaeota, Aigarchaeota, Crenarchaeota, and Korarchaeota (TACK) superphylum and Euryarchaeota contribute substantially to the niche expansion of "Ca. Aenigmarchaeota" via horizontal gene transfer (HGT), especially genes related to virus defense and stress responses. Based on co-occurrence network results and recent genetic exchanges among community members, we conjectured that "Ca. Aenigmarchaeota" may be symbionts associated with one MAG affiliated with the genus Pyrobaculum, though host specificity might be wide and variable across different "Ca. Aenigmarchaeota" organisms. This study provides significant insight into possible DPANN-host interactions and ecological roles of "Ca. Aenigmarchaeota." IMPORTANCE Recent advances in sequencing technology promoted the blowout discovery of super tiny microbes in the Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota, and Nanohaloarchaeota (DPANN) superphylum. However, the unculturable properties of the majority of microbes impeded our investigation of their behavior and symbiotic lifestyle in the corresponding community. By integrating horizontal gene transfer (HGT) detection and co-occurrence network analysis on "Candidatus Aenigmarchaeota" ("Ca. Aenigmarchaeota"), we made one of the first attempts to infer their putative interaction partners and further decipher the potential functional and genetic interactions between the symbionts. We revealed that HGTs contributed by members from the Thaumarchaeota, Aigarchaeota, Crenarchaeota, and Korarchaeota (TACK) superphylum and Euryarchaeota conferred "Ca. Aenigmarchaeota" with the ability to survive under different environmental stresses, such as virus infection, high temperature, and oxidative stress. This study demonstrates that the interaction partners might be inferable by applying informatics analyses on metagenomic sequencing data.

Comparative Genomics Reveals Thermal Adaptation and a High Metabolic Diversity in "Candidatus Bathyarchaeia"

"Candidatus Bathyarchaeia" is a phylogenetically diverse and widely distributed lineage often in high abundance in anoxic submarine sediments; however, their evolution and ecological roles in terrestrial geothermal habitats are poorly understood. In the present study, 35 Ca. Bathyarchaeia metagenome-assembled genomes (MAGs) were recovered from hot spring sediments in Tibet and Yunnan, China. Phylogenetic analysis revealed all MAGs of Ca. Bathyarchaeia can be classified into 7 orders and 15 families. Among them, 4 families have been first discovered in the present study, significantly expanding the known diversity of Ca. Bathyarchaeia. Comparative genomics demonstrated Ca. Bathyarchaeia MAGs from thermal habitats to encode a large variety of genes related to carbohydrate degradation, which are likely a metabolic adaptation of these organisms to a lifestyle at high temperatures. At least two families are potential methanogens/alkanotrophs, indicating a potential for the catalysis of short-chain hydrocarbons. Three MAGs from Family-7.3 are identified as alkanotrophs due to the detection of an Mcr complex. Family-2 contains the largest number of genes relevant to alkyl-CoM transformation, indicating the potential for methylotrophic methanogenesis, although their evolutionary history suggests the ancestor of Ca. Bathyarchaeia was unable to metabolize alkanes. Subsequent lineages have acquired the ability via horizontal gene transfer. Overall, our study significantly expands our knowledge and understanding of the metabolic capabilities, habitat adaptations, and evolution of Ca. Bathyarchaeia in thermal environments. IMPORTANCE Ca. Bathyarchaeia MAGs from terrestrial hot spring habitats are poorly revealed, though they have been studied extensively in marine ecosystems. In this study, we uncovered the metabolic capabilities and ecological role of Ca. Bathyarchaeia in hot springs and give a comprehensive comparative analysis between thermal and nonthermal habitats to reveal the thermal adaptability of Ca. Bathyarchaeia. Also, we attempt to determine the evolutionary history of methane/alkane metabolism in Ca. Bathyarchaeia, since it appears to be the first archaea beyond Euryarchaeota which contains the mcrABG genes. The reclassification of Ca. Bathyarchaeia and significant genomic differences among different lineages largely expand our knowledge on these cosmopolitan archaea, which will be beneficial in guiding the future studies.

A Novel Microbialite-Associated Phototrophic Chloroflexi Lineage Exhibiting a Quasi-Clonal Pattern along Depth

Chloroflexales (Chloroflexi) are typical members of the anoxygenic photosynthesizing component of microbial mats and have mostly been characterized from communities associated to hot springs. Here, we report the assembly of five metagenome-assembled genomes (MAGs) of a novel lineage of Chloroflexales found in mesophilic lithifying microbial mats (microbialites) in Lake Alchichica (Mexico). Genomic and phylogenetic analyses revealed that the bins shared 92% of their genes, and these genes were nearly identical despite being assembled from samples collected along a depth gradient (1-15 m depth). We tentatively name this lineage Candidatus Lithoflexus mexicanus. Metabolic predictions based on the MAGs suggest that these chlorosome-lacking mixotrophs share features in central carbon metabolism, electron transport, and adaptations to life under oxic and anoxic conditions, with members of two related lineages, Chloroflexineae and Roseiflexineae. Contrasting with the other diverse microbialite community members, which display much lower genomic conservation along the depth gradient, Ca. L. mexicanus MAGs exhibit remarkable similarity. This might reflect a particular flexibility to acclimate to varying light conditions with depth or the capacity to occupy a very specific spatial ecological niche in microbialites from different depths. Alternatively, Ca. L. mexicanus may also have the ability to modulate its gene expression as a function of the local environmental conditions during diel cycles in microbialites along the depth gradient.

Printing Microbial Dark Matter: Using Single Cell Dispensing and Genomics to Investigate the Patescibacteria/Candidate Phyla Radiation

As of today, the majority of environmental microorganisms remain uncultured. They are therefore referred to as "microbial dark matter." In the recent past, cultivation-independent methods like single-cell genomics (SCG) enabled the discovery of many previously unknown microorganisms, among them the Patescibacteria/Candidate Phyla Radiation (CPR). This approach was shown to be complementary to metagenomics, however, the development of additional and refined sorting techniques beyond the most commonly used fluorescence-activated cell sorting (FACS) is still desirable to enable additional downstream applications. Adding image information on the number and morphology of sorted cells would be beneficial, as would be minimizing cell stress caused by sorting conditions such as staining or pressure. Recently, a novel cell sorting technique has been developed, a microfluidic single-cell dispenser, which assesses the number and morphology of the cell in each droplet by automated light microscopic processing. Here, we report for the first time the successful application of the newly developed single-cell dispensing system for label-free isolation of individual bacteria from a complex sample retrieved from a wastewater treatment plant, demonstrating the potential of this technique for single cell genomics and other alternative downstream applications. Genome recovery success rated above 80% with this technique-out of 880 sorted cells 717 were successfully amplified. For 50.1% of these, analysis of the 16S rRNA gene was feasible and led to the sequencing of 50 sorted cells identified as Patescibacteria/CPR members. Subsequentially, 27 single amplified genomes (SAGs) of 15 novel and distinct Patescibacteria/CPR members, representing yet unseen species, genera and families could be captured and reconstructed. This phylogenetic distinctness of the recovered SAGs from available metagenome-assembled genomes (MAGs) is accompanied by the finding that these lineages-in whole or in part-have not been accessed by genome-resolved metagenomics of the same sample, thereby emphasizing the importance and opportunities of SCGs.

Mercury methylation by metabolically versatile and cosmopolitan marine bacteria

Microbes transform aqueous mercury (Hg) into methylmercury (MeHg), a potent neurotoxin that accumulates in terrestrial and marine food webs, with potential impacts on human health. This process requires the gene pair hgcAB, which encodes for proteins that actuate Hg methylation, and has been well described for anoxic environments. However, recent studies report potential MeHg formation in suboxic seawater, although the microorganisms involved remain poorly understood. In this study, we conducted large-scale multi-omic analyses to search for putative microbial Hg methylators along defined redox gradients in Saanich Inlet, British Columbia, a model natural ecosystem with previously measured Hg and MeHg concentration profiles. Analysis of gene expression profiles along the redoxcline identified several putative Hg methylating microbial groups, including Calditrichaeota, SAR324 and Marinimicrobia, with the last the most active based on hgc transcription levels. Marinimicrobia hgc genes were identified from multiple publicly available marine metagenomes, consistent with a potential key role in marine Hg methylation. Computational homology modelling predicts that Marinimicrobia HgcAB proteins contain the highly conserved amino acid sites and folding structures required for functional Hg methylation. Furthermore, a number of terminal oxidases from aerobic respiratory chains were associated with several putative novel Hg methylators. Our findings thus reveal potential novel marine Hg-methylating microorganisms with a greater oxygen tolerance and broader habitat range than previously recognized.

Protein Family Content Uncovers Lineage Relationships and Bacterial Pathway Maintenance Mechanisms in DPANN Archaea

DPANN are small-celled archaea that are generally predicted to be symbionts, and in some cases are known episymbionts of other archaea. As the monophyly of the DPANN remains uncertain, we hypothesized that proteome content could reveal relationships among DPANN lineages, constrain genetic overlap with bacteria, and illustrate how organisms with hybrid bacterial and archaeal protein sets might function. We tested this hypothesis using protein family content that was defined in part using 3,197 genomes including 569 newly reconstructed genomes. Protein family content clearly separates the final set of 390 DPANN genomes from other archaea, paralleling the separation of Candidate Phyla Radiation (CPR) bacteria from all other bacteria. This separation is partly driven by hypothetical proteins, some of which may be symbiosis-related. Pacearchaeota with the most limited predicted metabolic capacities have Form II/III and III-like Rubisco, suggesting metabolisms based on scavenged nucleotides. Intriguingly, the Pacearchaeota and Woesearchaeota with the smallest genomes also tend to encode large extracellular murein-like lytic transglycosylase domain proteins that may bind and degrade components of bacterial cell walls, indicating that some might be episymbionts of bacteria. The pathway for biosynthesis of bacterial isoprenoids is widespread in Woesearchaeota genomes and is encoded in proximity to genes involved in bacterial fatty acids synthesis. Surprisingly, in some DPANN genomes we identified a pathway for synthesis of queuosine, an unusual nucleotide in tRNAs of bacteria. Other bacterial systems are predicted to be involved in protein refolding. For example, many DPANN have the complete bacterial DnaK-DnaJ-GrpE system and many Woesearchaeota and Pacearchaeota possess bacterial group I chaperones. Thus, many DPANN appear to have mechanisms to ensure efficient protein folding of both archaeal and laterally acquired bacterial proteins.

Mercury methylation by metabolically versatile and cosmopolitan marine bacteria

Microbes transform aqueous mercury (Hg) into methylmercury (MeHg), a potent neurotoxin that accumulates in terrestrial and marine food webs, with potential impacts on human health. This process requires the gene pair hgcAB, which encodes for proteins that actuate Hg methylation, and has been well described for anoxic environments. However, recent studies report potential MeHg formation in suboxic seawater, although the microorganisms involved remain poorly understood. In this study, we conducted large-scale multi-omic analyses to search for putative microbial Hg methylators along defined redox gradients in Saanich Inlet, British Columbia, a model natural ecosystem with previously measured Hg and MeHg concentration profiles. Analysis of gene expression profiles along the redoxcline identified several putative Hg methylating microbial groups, including Calditrichaeota, SAR324 and Marinimicrobia, with the last the most active based on hgc transcription levels. Marinimicrobia hgc genes were identified from multiple publicly available marine metagenomes, consistent with a potential key role in marine Hg methylation. Computational homology modelling predicts that Marinimicrobia HgcAB proteins contain the highly conserved amino acid sites and folding structures required for functional Hg methylation. Furthermore, a number of terminal oxidases from aerobic respiratory chains were associated with several putative novel Hg methylators. Our findings thus reveal potential novel marine Hg-methylating microorganisms with a greater oxygen tolerance and broader habitat range than previously recognized.

The hunt for a yet unknown: Common molecular signature in some genetically monomorphic enterobacteria

Mark Achtman introduced the term "genetically monomorphic bacteria" (GM bacteria) for some human and plant pathogens. They displayed a great uniformity in terms of their "genetic" properties. This "uniformity" poses a challenge to microbiologists. To address these problems, we used CodonW and IslandViewer 3 as analytical tools and took Escherichia coli, Salmonella, and Shigella strains as a model organisms. We hypothesized that GM bacterium contains a common molecular signature among them. We have found a significant correlation regarding the number of protein-coding genes, predicted highly expressed genes, and the highest length of gene in this regard. On the other hand, the correspondence analysis of pathogenicity-related genes identified by IslandViewer 3 displayed a somewhat unique pattern in GM bacteria. The probable pathogenic genes are clustered into two separate groups, which is a hallmark of some pattern. Similar genes of non-monomorphic pathogenic strain clustered almost similarly, but the clusters are joined together, they are not completely separated. These features, in our considered view, may be considered as codon usages signatures of these bacteria, and E. coli in particular.