Untangling genomes from metagenomes: revealing an uncultured class of marine Euryarchaeota

From nets to networks: tools for deciphering phytoplankton metabolic interactions within communities and their global significance

Our oceans are populated with a wide diversity of planktonic organisms that form complex dynamic communities at the base of marine trophic networks. Within such communities are phytoplankton, unicellular photosynthetic taxa that provide an estimated half of global primary production and support biogeochemical cycles, along with other essential ecosystem services. One of the major challenges for microbial ecologists has been to try to make sense of this complexity. While phytoplankton distributions can be well explained by abiotic factors such as temperature and nutrient availability, there is increasing evidence that their ecological roles are tightly linked to their metabolic interactions with other plankton members through complex mechanisms (e.g. competition and symbiosis). Therefore, unravelling phytoplankton metabolic interactions is the key for inferring their dependency on, or antagonism with, other taxa and better integrating them into the context of carbon and nutrient fluxes in marine trophic networks. In this review, we attempt to summarize the current knowledge brought by ecophysiology, organismal imaging, in silico predictions and co-occurrence networks using 'omics data, highlighting successful combinations of approaches that may be helpful for future investigations of phytoplankton metabolic interactions within their complex communities.This article is part of the theme issue 'Connected interactions: enriching food web research by spatial and social interactions'.

Non-redundant metagenome-assembled genomes of activated sludge reactors at different disturbances and scales

Metagenome-assembled genomes (MAGs) are microbial genomes reconstructed from metagenomic data and can be assigned to known taxa or lead to uncovering novel ones. MAGs can provide insights into how microbes interact with the environment. Here, we performed genome-resolved metagenomics on sequencing data from four studies using sequencing batch reactors at microcosm (~25 mL) and mesocosm (~4 L) scales inoculated with sludge from full-scale wastewater treatment plants. These studies investigated how microbial communities in such plants respond to two environmental disturbances: the presence of toxic 3-chloroaniline and changes in organic loading rate. We report 839 non-redundant MAGs with at least 50% completeness and 10% contamination (MIMAG medium-quality criteria). From these, 399 are of putative high-quality, while sixty-seven meet the MIMAG high-quality criteria. MAGs in this catalogue represent the microbial communities in sixty-eight laboratory-scale reactors used for the disturbance experiments, and in the full-scale wastewater treatment plant which provided the source sludge. This dataset can aid meta-studies aimed at understanding the responses of microbial communities to disturbances, particularly as ecosystems confront rapid environmental changes.

AFITbin: a metagenomic contig binning method using aggregate l-mer frequency based on initial and terminal nucleotides

BACKGROUND

Using next-generation sequencing technologies, scientists can sequence complex microbial communities directly from the environment. Significant insights into the structure, diversity, and ecology of microbial communities have resulted from the study of metagenomics. The assembly of reads into longer contigs, which are then binned into groups of contigs that correspond to different species in the metagenomic sample, is a crucial step in the analysis of metagenomics. It is necessary to organize these contigs into operational taxonomic units (OTUs) for further taxonomic profiling and functional analysis. For binning, which is synonymous with the clustering of OTUs, the tetra-nucleotide frequency (TNF) is typically utilized as a compositional feature for each OTU.

RESULTS

In this paper, we present AFIT, a new l-mer statistic vector for each contig, and AFITBin, a novel method for metagenomic binning based on AFIT and a matrix factorization method. To evaluate the performance of the AFIT vector, the t-SNE algorithm is used to compare species clustering based on AFIT and TNF information. In addition, the efficacy of AFITBin is demonstrated on both simulated and real datasets in comparison to state-of-the-art binning methods such as MetaBAT 2, MaxBin 2.0, CONCOT, MetaCon, SolidBin, BusyBee Web, and MetaBinner. To further analyze the performance of the purposed AFIT vector, we compare the barcodes of the AFIT vector and the TNF vector.

CONCLUSION

The results demonstrate that AFITBin shows superior performance in taxonomic identification compared to existing methods, leveraging the AFIT vector for improved results in metagenomic binning. This approach holds promise for advancing the analysis of metagenomic data, providing more reliable insights into microbial community composition and function.

AVAILABILITY

A python package is available at: https://github.com/SayehSobhani/AFITBin .

Archaeal blooms and busts in an estuarine time series

Coastal bays, such as Delaware Bay, are highly productive, ecologically important transitions between rivers and the coastal ocean. They offer opportunities to investigate archaeal assemblages across seasons, with the exchange of water masses that occurs with tidal cycles, and in the context of variable organic matter quality. For a year-long estuarine, size-fractionated time series, we used amplicon sequencing, chemical measurements, and qPCR to follow archaeal groups through the seasons. We detected seasonally high abundances of Marine Group II archaea in summer months which correlate with indicators of phytoplankton production, although not phytoplankton biomass. Although previous studies have reported associations between Marine Group II archaea and particles, here they are almost entirely found in very small particles (0.22-0.7 μm), suggesting they are free-living cells. Populations of Nitrososphaeria did not vary with particle size or environmental conditions. Methanogens were significant fractions of archaeal sequences in large particles at low tide during winter months. Contrary to expectations, Nanoarchaeia were found predominantly in the free-living fraction despite the previous observation that they require an association with hosts. These results underscore the utility of time series studies in shallow, tidally mixed estuarine environments that capture variable conditions for understanding the ecology and biogeochemistry of planktic archaea.

Chronic industrial perturbation and seasonal change induces shift in the bacterial community from gammaproteobacteria to betaproteobacteria having catabolic potential for aromatic compounds at Amlakhadi canal

Escalating proportions of industrially contaminated sites are one of the major catastrophes faced at the present time due to the industrial revolution. The difficulties associated with culturing the microbes, has been circumvent by the direct use of metagenomic analysis of various complex niches. In this study, a metagenomic approach using next generation sequencing technologies was applied to exemplify the taxonomic abundance and metabolic potential of the microbial community residing in Amlakhadi canal, Ankleshwar at two different seasons. All the metagenomes revealed a predominance of Proteobacteria phylum. However, difference was observed within class level where Gammaproteobacteria was relatively high in polluted metagenome in Summer while in Monsoon the abundance shifted to Betaproteobacteria. Similarly, significant statistical differences were obtained while comparing the genera amongst contaminated sites where Serratia, Achromobacter, Stenotrophomonas and Pseudomonas were abundant in summer season and the dominance changed to Thiobacillus, Thauera, Acidovorax, Nitrosomonas, Sulfuricurvum, Novosphingobium, Hyphomonas and Geobacter in monsoon. Further upon functional characterization, the microbiomes revealed the diverse survival mechanisms, in response to the prevailing ecological conditions (such as degradation of aromatic compounds, heavy metal resistance, oxidative stress responses and multidrug resistance efflux pumps, etc.). The results have important implications in understanding and predicting the impacts of human-induced activities on microbial communities inhabiting natural niche and their responses in coping with the fluctuating pollution load.

Dancing the Nanopore limbo - Nanopore metagenomics from small DNA quantities for bacterial genome reconstruction

BACKGROUND

While genome-resolved metagenomics has revolutionized our understanding of microbial and genetic diversity in environmental samples, assemblies of short-reads often result in incomplete and/or highly fragmented metagenome-assembled genomes (MAGs), hampering in-depth genomics. Although Nanopore sequencing has increasingly been used in microbial metagenomics as long reads greatly improve the assembly quality of MAGs, the recommended DNA quantity usually exceeds the recoverable amount of DNA of environmental samples. Here, we evaluated lower-than-recommended DNA quantities for Nanopore library preparation by determining sequencing quality, community composition, assembly quality and recovery of MAGs.

RESULTS

We generated 27 Nanopore metagenomes using the commercially available ZYMO mock community and varied the amount of input DNA from 1000 ng (the recommended minimum) down to 1 ng in eight steps. The quality of the generated reads remained stable across all input levels. The read mapping accuracy, which reflects how well the reads match a known reference genome, was consistently high across all libraries. The relative abundance of the species in the metagenomes was stable down to input levels of 50 ng. High-quality MAGs (> 95% completeness, ≤ 5% contamination) could be recovered from metagenomes down to 35 ng of input material. When combined with publicly available Illumina reads for the mock community, Nanopore reads from input quantities as low as 1 ng improved the quality of hybrid assemblies.

CONCLUSION

Our results show that the recommended DNA amount for Nanopore library preparation can be substantially reduced without any adverse effects to genome recovery and still bolster hybrid assemblies when combined with short-read data. We posit that the results presented herein will enable studies to improve genome recovery from low-biomass environments, enhancing microbiome understanding.

Features of the Mechanism of Proton Transport in ESR, Retinal Protein from Exiguobacterium sibiricum

Retinal-containing light-sensitive proteins - rhodopsins - are found in many microorganisms. Interest in them is largely explained by their role in light energy storage and photoregulation in microorganisms, as well as the prospects for their use in optogenetics to control neuronal activity, including treatment of various diseases. One of the representatives of microbial rhodopsins is ESR, the retinal protein of Exiguobacterium sibiricum. What distinguishes ESR from homologous proteins is the presence of a lysine residue (Lys96) as a proton donor for the Schiff base. This feature, along with the hydrogen bond of the proton acceptor Asp85 with the His57 residue, determines functional characteristics of ESR as a proton pump. This review examines the results of ESR studies conducted using various methods, including direct electrometry. Comparison of the obtained data with the results of structural studies and with other retinal proteins allows us to draw conclusions about the mechanisms of transport of hydrogen ions in ESR and similar retinal proteins.

A moderately thermophilic origin of a novel family of marine group II euryarchaeota from deep ocean

Marine group II (MGII) is the most abundant planktonic heterotrophic archaea in the ocean. The evolutionary history of MGII archaea is elusive. In this study, 13 new MGII metagenome-assembled genomes were recovered from surface to the hadal zone in Challenger Deep of the Mariana Trench; four of them from the deep ocean represent a novel group. The optimal growth temperature (OGT) of the common ancestor of MGII has been estimated to be at about 60°C and OGTs of MGIIc, MGIIb, and MGIIa at 47°C-50ºC, 37°C-44ºC, and 30°C-37ºC, respectively, suggesting the adaptation of these species to different temperatures during evolution. The estimated OGT range of MGIIc was supported by experimental measurements of cloned β-galactosidase that showed optimal enzyme activity around 50°C. These results indicate that MGIIc may have originated from a common ancestor that lived in warm or even hot marine environment, such as hydrothermal vents.

Vertical segregation and phylogenetic characterization of archaea and archaeal ammonia monooxygenase gene in the water column of the western Arctic Ocean

Archaea constitute a substantial fraction of marine microbial biomass and play critical roles in the biogeochemistry of oceans. However, studies on their distribution and ecology in the Arctic Ocean are relatively scarce. Here, we studied the distributions of archaea and archaeal ammonia monooxygenase (amoA) gene in the western Arctic Ocean, using the amplicon sequencing approach from the sea surface to deep waters up to 3040 m depth. A total of five archaeal phyla, Nitrososphaerota, "Euryarchaeota", "Halobacteriota," "Nanoarchaeota", and Candidatus Thermoplasmatota, were detected. We observed a clear, depth-dependent vertical segregation among archaeal communities. Ca. Thermoplasmatota (66.8%) was the most dominant phylum in the surface waters. At the same time, Nitrososphaerota (55.9%) was dominant in the deep waters. Most of the amoA gene OTUs (99%) belonged to the Nitrosopumilales and were further clustered into five subclades ("NP-Alpha", "NP-Delta", "NP-Epsilon", "NP-Gamma", and "NP-Theta"). "NP-Epsilon" was the most dominant clade throughout the water column and "NP_Alpha" showed higher abundance only in the deeper water. Salinity and inorganic nutrient concentrations were the major factors that determined the vertical segregation of archaea. We anticipate that the observed differences in the vertical distribution of archaea might contribute to the compartmentalization of dark carbon fixation and nitrification in deeper water and organic matter degradation in surface waters of the Arctic Ocean.

Compositional Changes in Sediment Microbiota Are Associated with Seasonal Variation of the Water Column in High-Altitude Hyperarid Andean Lake Systems

The lacustrine systems of La Brava and La Punta, located in the Tilopozo sector in the extreme south of Salar de Atacama, are pristine high-altitude Andean lakes found along the central Andes of South America. This shallow ecosystem suffers from permanent evaporation, leading to falling water levels, causing it to recede or disappear during the dry season. This dynamic causes physicochemical changes in lakes, such as low nutrient availability, pH change, and dissolved metals, which can influence the composition of the microbial community. In this study, we used a metataxonomic approach (16S rRNA hypervariable regions V3 to V4) to characterize the sedimentary microbiota of these lakes. To understand how the water column affects and is structured in the microbiota of these lakes, we combined the analysis of the persistence of the water column through satellite images and physicochemical characterization. Our results show a significant difference in abiotic factors and microbiota composition between La Punta and La Brava lakes. In addition, microbiota analysis revealed compositional changes in the ecological disaggregation (main and isolated bodies) and antagonistic changes in the abundance of certain taxa between lakes. These findings are an invaluable resource for understanding the microbiological diversity of high Andean lakes using a multidisciplinary approach that evaluates the microbiota behavior in response to abiotic factors. IMPORTANCE In this study, we analyzed the persistence of the water column through satellite images and physicochemical characterization to investigate the composition and diversity in High Andean Lake Systems in a hyperarid environment. In addition to the persistence of the water column, this approach can be used to analyze changes in the morphology of saline accumulations and persistence of snow or ice; for example, for establishing variable plant cover over time and evaluating the microbiota associated with soils with seasonal changes in plants. This makes it an ideal approach to search for novel extremophilic microorganisms with unique properties. In our case, it was used to study microorganisms capable of resisting desiccation and water restriction for a considerable period and adapting to survive in ecological niches, such as those with high UV irradiation, extreme drought, and high salt concentration.

A compendium of bacterial and archaeal single-cell amplified genomes from oxygen deficient marine waters

Oxygen-deficient marine waters referred to as oxygen minimum zones (OMZs) or anoxic marine zones (AMZs) are common oceanographic features. They host both cosmopolitan and endemic microorganisms adapted to low oxygen conditions. Microbial metabolic interactions within OMZs and AMZs drive coupled biogeochemical cycles resulting in nitrogen loss and climate active trace gas production and consumption. Global warming is causing oxygen-deficient waters to expand and intensify. Therefore, studies focused on microbial communities inhabiting oxygen-deficient regions are necessary to both monitor and model the impacts of climate change on marine ecosystem functions and services. Here we present a compendium of 5,129 single-cell amplified genomes (SAGs) from marine environments encompassing representative OMZ and AMZ geochemical profiles. Of these, 3,570 SAGs have been sequenced to different levels of completion, providing a strain-resolved perspective on the genomic content and potential metabolic interactions within OMZ and AMZ microbiomes. Hierarchical clustering confirmed that samples from similar oxygen concentrations and geographic regions also had analogous taxonomic compositions, providing a coherent framework for comparative community analysis.

Microbe, climate change and marine environment: Linking trends and research hotspots

Microbes, or microorganisms, have been the foundation of the biosphere for over 3 billion years and have played an essential role in shaping our planet. The available knowledge on the topic of microbes associated with climate change has the potential to reshape upcoming research trends globally. As climate change impacts the ocean or marine ecosystem, the responses of these "unseen life" will heavily influence the achievement of a sustainable evolutionary environment. The present study aims to identify microbial-related research under changing climate within the marine environment through the mapping of visualized graphs of the available literature. We used scientometric methods to retrieve documents from the Web of Science platform in the Core Collection (WOSCC) database, analyzing a total of 2767 documents based on scientometric indicators. Our findings show that this research area is growing exponentially, with the most influential keywords being "microbial diversity," "bacteria," and "ocean acidification," and the most cited being "microorganism" and "diversity." The identification of influential clusters in the field of marine science provides insight into the hot spots and frontiers of research in this area. Prominent clusters include "coral microbiome," "hypoxic zone," "novel Thermoplasmatota clade," "marine dinoflagellate bloom," and "human health." Analyzing emerging trends and transformative changes in this field can inform the creation of special issues or research topics in selected journals, thus increasing visibility and engagement among the scientific community.

Bacteria and Archaea Regulate Particulate Organic Matter Export in Suspended and Sinking Marine Particle Fractions

The biological carbon pump (BCP) in the Southern Ocean is driven by phytoplankton productivity and is a significant organic matter sink. However, the role of particle-attached (PA) and free-living (FL) prokaryotes (bacteria and archaea) and their diversity in influencing the efficiency of the BCP is still unclear. To investigate this, we analyzed the metagenomes linked to suspended and sinking marine particles from the Sub-Antarctic Southern Ocean Time Series (SOTS) by deploying a Marine Snow Catcher (MSC), obtaining suspended and sinking particulate material, determining organic carbon and nitrogen flux, and constructing metagenome-assembled genomes (MAGs). The suspended and sinking particle-pools were dominated by bacteria with the potential to degrade organic carbon. Bacterial communities associated with the sinking fraction had more genes related to the degradation of complex organic carbon than those in the suspended fraction. Archaea had the potential to drive nitrogen metabolism via nitrite and ammonia oxidation, altering organic nitrogen concentration. The data revealed several pathways for chemoautotrophy and the secretion of recalcitrant dissolved organic carbon (RDOC) from CO₂, with bacteria and archaea potentially sequestering particulate organic matter (POM) via the production of RDOC. These findings provide insights into the diversity and function of prokaryotes in suspended and sinking particles and their role in organic carbon/nitrogen export in the Southern Ocean. IMPORTANCE The biological carbon pump is crucial for the export of particulate organic matter in the ocean. Recent studies on marine microbes have shown the profound influence of bacteria and archaea as regulators of particulate organic matter export. Yet, despite the importance of the Southern Ocean as a carbon sink, we lack comparable insights regarding microbial contributions. This study provides the first insights regarding prokaryotic contributions to particulate organic matter export in the Southern Ocean. We reveal evidence that prokaryotic communities in suspended and sinking particle fractions harbor widespread genomic potential for mediating particulate organic matter export. The results substantially enhance our understanding of the role played by microorganisms in regulating particulate organic matter export in suspended and sinking marine fractions in the Southern Ocean.

Ecological divergence of syntopic marine bacterial species is shaped by gene content and expression

Identifying mechanisms by which bacterial species evolve and maintain genomic diversity is particularly challenging for the uncultured lineages that dominate the surface ocean. A longitudinal analysis of bacterial genes, genomes, and transcripts during a coastal phytoplankton bloom revealed two co-occurring, highly related Rhodobacteraceae species from the deeply branching and uncultured NAC11-7 lineage. These have identical 16S rRNA gene amplicon sequences, yet their genome contents assembled from metagenomes and single cells indicate species-level divergence. Moreover, shifts in relative dominance of the species during dynamic bloom conditions over 7 weeks confirmed the syntopic species' divergent responses to the same microenvironment at the same time. Genes unique to each species and genes shared but divergent in per-cell inventories of mRNAs accounted for 5% of the species' pangenome content. These analyses uncover physiological and ecological features that differentiate the species, including capacities for organic carbon utilization, attributes of the cell surface, metal requirements, and vitamin biosynthesis. Such insights into the coexistence of highly related and ecologically similar bacterial species in their shared natural habitat are rare.

To assemble or not to assemble: metagenomic profiling of microbially mediated biogeochemical pathways in complex communities

High-throughput profiling of microbial functional traits involved in various biogeochemical cycling pathways using shotgun metagenomic sequencing has been routinely applied in microbial ecology and environmental science. Multiple bioinformatics data processing approaches are available, including assembly-based (single-sample assembly and multi-sample assembly) and read-based (merged reads and raw data). However, it remains not clear how these different approaches may differ in data analyses and affect result interpretation. In this study, using two typical shotgun metagenome datasets recovered from geographically distant coastal sediments, the performance of different data processing approaches was comparatively investigated from both technical and biological/ecological perspectives. Microbially mediated biogeochemical cycling pathways, including nitrogen cycling, sulfur cycling and B12 biosynthesis, were analyzed. As a result, multi-sample assembly provided the most amount of usable information for targeted functional traits, at a high cost of computational resources and running time. Single-sample assembly and read-based analysis were comparable in obtaining usable information, but the former was much more time- and resource-consuming. Critically, different approaches introduced much stronger variations in microbial profiles than biological differences. However, community-level differences between the two sampling sites could be consistently observed despite the approaches being used. In choosing an appropriate approach, researchers shall balance the trade-offs between multiple factors, including the scientific question, the amount of usable information, computational resources and time cost. This study is expected to provide valuable technical insights and guidelines for the various approaches used for metagenomic data analysis.

Diversity, distribution, and functional potentials of magroviruses from marine and brackish waters

Marine group II (MGII) archaea (Ca. Poseidoniales) are among the most abundant microbes in global oceanic surface waters and play an important role in driving marine biogeochemical cycles. Magroviruses - the viruses of MGII archaea have been recently found to occur ubiquitously in surface ocean. However, their diversity, distribution, and potential ecological functions in coastal zones especially brackish waters are unknown. Here we obtained 234 non-redundant magroviral genomes from brackish surface waters by using homology searches for viral signature proteins highlighting the uncovered vast diversity of this novel viral group. Phylogenetic analysis based on these brackish magroviruses along with previously reported marine ones identified six taxonomic groups with close evolutionary connection to both haloviruses and the viruses of Marine Group I archaea. Magroviruses were present abundantly both in brackish and open ocean samples with some showing habitat specification and others having broad spectrums of distribution between different habitats. Genome annotation suggests they may be involved in regulating multiple metabolic pathways of MGII archaea. Our results uncover the previously overlooked diversity and ecological potentials of a major archaeal virial group in global ocean and brackish waters and shed light on the cryptic evolutionary history of archaeal viruses.

Long-Read Metagenomics and CAZyme Discovery

Microorganisms play a primary role in regulating biogeochemical cycles and are a valuable source of enzymes that have biotechnological applications, such as carbohydrate-active enzymes (CAZymes). However, the inability to culture the majority of microorganisms that exist in natural ecosystems restricts access to potentially novel bacteria and beneficial CAZymes. While commonplace molecular-based culture-independent methods such as metagenomics enable researchers to study microbial communities directly from environmental samples, recent progress in long-read sequencing technologies are advancing the field. We outline key methodological stages that are required as well as describe specific protocols that are currently used for long-read metagenomic projects dedicated to CAZyme discovery.

Persistent memory as an effective alternative to random access memory in metagenome assembly

BACKGROUND

The assembly of metagenomes decomposes members of complex microbe communities and allows the characterization of these genomes without laborious cultivation or single-cell metagenomics. Metagenome assembly is a process that is memory intensive and time consuming. Multi-terabyte sequences can become too large to be assembled on a single computer node, and there is no reliable method to predict the memory requirement due to data-specific memory consumption pattern. Currently, out-of-memory (OOM) is one of the most prevalent factors that causes metagenome assembly failures.

RESULTS

In this study, we explored the possibility of using Persistent Memory (PMem) as a less expensive substitute for dynamic random access memory (DRAM) to reduce OOM and increase the scalability of metagenome assemblers. We evaluated the execution time and memory usage of three popular metagenome assemblers (MetaSPAdes, MEGAHIT, and MetaHipMer2) in datasets up to one terabase. We found that PMem can enable metagenome assemblers on terabyte-sized datasets by partially or fully substituting DRAM. Depending on the configured DRAM/PMEM ratio, running metagenome assemblies with PMem can achieve a similar speed as DRAM, while in the worst case it showed a roughly two-fold slowdown. In addition, different assemblers displayed distinct memory/speed trade-offs in the same hardware/software environment.

CONCLUSIONS

We demonstrated that PMem is capable of expanding the capacity of DRAM to allow larger metagenome assembly with a potential tradeoff in speed. Because PMem can be used directly without any application-specific code modification, these findings are likely to be generalized to other memory-intensive bioinformatics applications.

Whole-genome scanning reveals environmental selection mechanisms that shape diversity in populations of the epipelagic diatom Chaetoceros

Diatoms form a diverse and abundant group of photosynthetic protists that are essential players in marine ecosystems. However, the microevolutionary structure of their populations remains poorly understood, particularly in polar regions. Exploring how closely related diatoms adapt to different environments is essential given their short generation times, which may allow rapid adaptations, and their prevalence in marine regions dramatically impacted by climate change, such as the Arctic and Southern Oceans. Here, we address genetic diversity patterns in Chaetoceros, the most abundant diatom genus and one of the most diverse, using 11 metagenome-assembled genomes (MAGs) reconstructed from Tara Oceans metagenomes. Genome-resolved metagenomics on these MAGs confirmed a prevalent distribution of Chaetoceros in the Arctic Ocean with lower dispersal in the Pacific and Southern Oceans as well as in the Mediterranean Sea. Single-nucleotide variants identified within the different MAG populations allowed us to draw a landscape of Chaetoceros genetic diversity and revealed an elevated genetic structure in some Arctic Ocean populations. Gene flow patterns of closely related Chaetoceros populations seemed to correlate with distinct abiotic factors rather than with geographic distance. We found clear positive selection of genes involved in nutrient availability responses, in particular for iron (e.g., ISIP2a, flavodoxin), silicate, and phosphate (e.g., polyamine synthase), that were further supported by analysis of Chaetoceros transcriptomes. Altogether, these results highlight the importance of environmental selection in shaping diatom diversity patterns and provide new insights into their metapopulation genomics through the integration of metagenomic and environmental data.

Seasonal microbial dynamics in the ocean inferred from assembled and unassembled data: a view on the unknown biosphere

In environmental metagenomic experiments, a very high proportion of the microbial sequencing data (> 70%) remains largely unexploited because rare and closely related genomes are missed in short-read assemblies. The identity and the potential metabolisms of a large fraction of natural microbial communities thus remain inaccessible to researchers. The purpose of this study was to explore the genomic content of unassembled metagenomic data and test their level of novelty. We used data from a three-year microbial metagenomic time series of the NW Mediterranean Sea, and conducted reference-free and database-guided analysis. The results revealed a significant genomic difference between the assembled and unassembled reads. The unassembled reads had a lower mean identity against public databases, and fewer metabolic pathways could be reconstructed. In addition, the unassembled fraction presented a clear temporal pattern, unlike the assembled ones, and a specific community composition that was similar to the rare communities defined by metabarcoding using the 16S rRNA gene. The rare gene pool was characterised by keystone bacterial taxa, and the presence of viruses, suggesting that viral lysis could maintain some taxa in a state of rarity. Our study demonstrates that unassembled metagenomic data can provide important information on the structure and functioning of microbial communities.

Proton transfer reactions in donor site mutants of ESR, a retinal protein from Exiguobacterium sibiricum

Light-driven proton transport by microbial retinal proteins such as archaeal bacteriorhodopsin involves carboxylic residues as internal proton donors to the catalytic center which is a retinal Schiff base (SB). The proton donor, Asp96 in bacteriorhodopsin, supplies a proton to the transiently deprotonated Schiff base during the photochemical cycle. Subsequent proton uptake resets the protonated state of the donor. This two step process became a distinctive signature of retinal based proton pumps. Similar steps are observed also in many natural variants of bacterial proteorhodopsins and xanthorhodopsins where glutamic acid residues serve as a proton donor. Recently, however, an exception to this rule was found. A retinal protein from Exiguobacterium sibiricum, ESR, contains a Lys residue in place of Asp or Glu, which facilitates proton transfer from the bulk to the SB. Lys96 can be functionally replaced with the more common donor residues, Asp or Glu. Proton transfer to the SB in the mutants containing these replacements (K96E and K96D/A47T) is much faster than in the proteins lacking the proton donor (K96A and similar mutants), and in the case of K96D/A47T, comparable with that in the wild type, indicating that carboxylic residues can replace Lys96 as proton donors in ESR. We show here that there are important differences in the functioning of these residues in ESR from the way Asp96 functions in bacteriorhodopsin. Reprotonation of the SB and proton uptake from the bulk occur almost simultaneously during the M to N transition (as in the wild type ESR at neutral pH), whereas in bacteriorhodopsin these two steps are well separated in time and occur during the M to N and N to O transitions, respectively.

Composition and Biogeography of Planktonic Pro- and Eukaryotic Communities in the Atlantic Ocean: Primer Choice Matters

Basin-scale biogeographic observations of marine pelagic pro- and eukaryotic communities are necessary to understand forces driving community composition and for providing a baseline to monitor global change. Deep sequencing of rRNA genes provides community composition at high resolution; yet, it is unclear how the choice of primers affects biogeographic patterns. Here, we re-amplified 16S rRNA genes from DNA sampled during R/V Polarstern Cruise ANT28-5 over a latitudinal transect across the Atlantic Ocean from 52°S to 47°N using universal V4-V5 primers and compared the results with those obtained previously with V5-V6 bacteria-specific primers. For validation of our results, we inferred community composition based on 16S rRNA genes of metagenomes from the same stations and single amplified genomes (SAGs) from the Global Ocean Reference Genome (GORG) database. We found that the universal V4-V5 primers retrieved SAR11 clades with similar relative proportions as those found in the GORG database while the V5-V6 primers recovered strongly diverging clade abundances. We confirmed an inverse bell-shaped distance-decay relationship and a latitudinal diversity gradient that did not decline linearly with absolute latitude in the Atlantic Ocean. Patterns were modified by sampling depth, sequencing depth, choice of primers, and abundance filtering. Especially richness patterns were not robust to methodological change. This study offers a detailed picture of the Atlantic Ocean microbiome using a universal set of PCR primers that allow for the conjunction of biogeographical patterns among organisms from different domains of life.

Availability of vitamin B12 and its lower ligand intermediate α-ribazole impact prokaryotic and protist communities in oceanic systems

Genome analyses predict that the cofactor cobalamin (vitamin B₁₂, called B₁₂ herein) is produced by only one-third of all prokaryotes but almost all encode at least one B₁₂-dependent enzyme, in most cases methionine synthase. This implies that the majority of prokaryotes relies on exogenous B₁₂ supply and interacts with producers. B₁₂ consists of a corrin ring centred around a cobalt ion and the lower ligand 5’6-dimethylbenzimidazole (DMB). It has never been tested whether availability of this pivotal cofactor, DMB or its intermediate α-ribazole affect growth and composition of prokaryotic microbial communities. Here we show that in the subtropical, equatorial and polar frontal Pacific Ocean supply of B₁₂ and α-ribazole enhances heterotrophic prokaryotic production and alters the composition of prokaryotic and heterotrophic protist communities. In the polar frontal Pacific, the SAR11 clade and Oceanospirillales increased their relative abundances upon B₁₂ supply. In the subtropical Pacific, Oceanospirillales increased their relative abundance upon B₁₂ supply as well but also downregulated the transcription of the btuB gene, encoding the outer membrane permease for B₁₂. Surprisingly, Prochlorococcus, known to produce pseudo-B₁₂ and not B₁₂, exhibited significant upregulation of genes encoding key proteins of photosystem I + II, carbon fixation and nitrate reduction upon B₁₂ supply in the subtropical Pacific. These findings show that availability of B₁₂ and α-ribazole affect growth and composition of prokaryotic and protist communities in oceanic systems thus revealing far-reaching consequences of methionine biosynthesis and other B₁₂-dependent enzymatic reactions on a community level.

A deep siamese neural network improves metagenome-assembled genomes in microbiome datasets across different environments

Metagenomic binning is the step in building metagenome-assembled genomes (MAGs) when sequences predicted to originate from the same genome are automatically grouped together. The most widely-used methods for binning are reference-independent, operating de novo and enable the recovery of genomes from previously unsampled clades. However, they do not leverage the knowledge in existing databases. Here, we introduce SemiBin, an open source tool that uses deep siamese neural networks to implement a semi-supervised approach, i.e. SemiBin exploits the information in reference genomes, while retaining the capability of reconstructing high-quality bins that are outside the reference dataset. Using simulated and real microbiome datasets from several different habitats from GMGCv1 (Global Microbial Gene Catalog), including the human gut, non-human guts, and environmental habitats (ocean and soil), we show that SemiBin outperforms existing state-of-the-art binning methods. In particular, compared to other methods, SemiBin returns more high-quality bins with larger taxonomic diversity, including more distinct genera and species.

Biogeographic, Driving Factors, Assembly, and Co-occurrence Patterns of Archaeal Community in Biocrusts

Archaea exhibit strong community heterogeneity with microhabitat gradients and are a non-negligible part of biocrust's microorganisms. The study on archaeal biogeography in biocrusts could provide new insights for its application in environmental restoration. However, only a few studies on assembly processes and co-occurrence patterns of the archaeal community in patchy biocrusts have been reported, especially considering the number of species pools (SPs). Here, we comprehensively collected biocrusts across 3,500 km of northern China. Different successional biocrusts from various regions contain information of local climate and microenvironments, which can shape multiple unique archaeal SPs. The archaeal community differences in the same successional stage exceeded the variations between successional stages, which was due to the fact that the heterogeneous taxa tended to exchange between unknown patches driven by drift. We also comparatively studied the driving forces of community heterogeneity across three to ten SPs, and assembly and co-occurrence patterns were systematically analyzed. The results revealed that the impact of spatial factors on biogeographic patterns was greater than that of environmental and successional factors and that impact decreased with the number of SPs considered. Meanwhile, community heterogeneity at the phylogenetic facet was more sensitive to these driving factors than the taxonomic facet. Subgroups 1 (SG1) and 2 (SG2) of the archaeal communities in biocrusts were dominated by Nitrososphaeraceae and Haloarchaea, respectively. The former distribution pattern was associated with non-salinity-related variables and primarily assembled by drift, whereas the latter was associated with salinity-related variables and primarily assembled by homogeneous selection. Finally, network analysis indicated that the SG1 network had a higher proportion of competition and key taxa than the SG2 network, but the network of SG2 was more complex. Our study suggested that the development of the archaeal community was not consistent with biocrusts succession. The dominant taxa may determine the patterns of community biogeography, assembly, and co-occurrence.

Biodiversity and Biogeography of Abundant and Rare Microbial Assemblages in the Western Subtropical Pacific Ocean

The levels of chlorophyll a and nutrient concentrations in the surface waters of the western subtropical Pacific Ocean are among the lowest globally. In addition, our knowledge of basin-scale diversity and biogeography of microbial communities in this vast extremely oligotrophic environment is still rather limited. Here, high-throughput sequencing was used to examine the biodiversity and biogeography of abundant and rare microbial assemblages throughout the water column from the surface to a depth of 3,000 m across a horizontal distance of 1,100 km in the western Pacific Ocean. Microbial alpha diversity in the 200-m layer was higher than at other depths, with Gammaproteobacteria, Alphaproteobacteria, and Clostridia as the dominant classes in all samples. Distinctly vertical distributions within the microbial communities were revealed, with no difference horizontally. Some microbes exhibited depth stratification. For example, the relative abundances of Cyanobacteria and Alphaproteobacteria decreased with depth, while Nitrososphaeria, Actinobacteria, and Gammaproteobacteria increased with depth in the aphotic layers. Furthermore, we found that environmental (selective process) and spatial (neutral process) factors had different effects on abundant and rare taxa. Geographical distance showed little effect on the dispersal of all and abundant taxa, while statistically significant distance-decay relationships were observed among the rare taxa. Temperature and chlorophyll a were strongly associated with all, abundant, and rare taxa in the photic layers, while total inorganic nitrogen was recognized as the crucial factor in the aphotic layers. Variance partitioning analysis indicated that environmental selection played a relatively important role in shaping all and abundant taxa, while the variation in rare taxa explained by environmental and spatial processes was relatively low, as more than 70% of the variation remained unexplained. This study provides novel knowledge related to microbial community diversity in the western subtropical Pacific Ocean, and the analyzes biogeographical patterns among abundant and rare taxa.

METABOLIC: high-throughput profiling of microbial genomes for functional traits, metabolism, biogeochemistry, and community-scale functional networks

BACKGROUND

Advances in microbiome science are being driven in large part due to our ability to study and infer microbial ecology from genomes reconstructed from mixed microbial communities using metagenomics and single-cell genomics. Such omics-based techniques allow us to read genomic blueprints of microorganisms, decipher their functional capacities and activities, and reconstruct their roles in biogeochemical processes. Currently available tools for analyses of genomic data can annotate and depict metabolic functions to some extent; however, no standardized approaches are currently available for the comprehensive characterization of metabolic predictions, metabolite exchanges, microbial interactions, and microbial contributions to biogeochemical cycling.

RESULTS

We present METABOLIC (METabolic And BiogeOchemistry anaLyses In miCrobes), a scalable software to advance microbial ecology and biogeochemistry studies using genomes at the resolution of individual organisms and/or microbial communities. The genome-scale workflow includes annotation of microbial genomes, motif validation of biochemically validated conserved protein residues, metabolic pathway analyses, and calculation of contributions to individual biogeochemical transformations and cycles. The community-scale workflow supplements genome-scale analyses with determination of genome abundance in the microbiome, potential microbial metabolic handoffs and metabolite exchange, reconstruction of functional networks, and determination of microbial contributions to biogeochemical cycles. METABOLIC can take input genomes from isolates, metagenome-assembled genomes, or single-cell genomes. Results are presented in the form of tables for metabolism and a variety of visualizations including biogeochemical cycling potential, representation of sequential metabolic transformations, community-scale microbial functional networks using a newly defined metric "MW-score" (metabolic weight score), and metabolic Sankey diagrams. METABOLIC takes ~ 3 h with 40 CPU threads to process ~ 100 genomes and corresponding metagenomic reads within which the most compute-demanding part of hmmsearch takes ~ 45 min, while it takes ~ 5 h to complete hmmsearch for ~ 3600 genomes. Tests of accuracy, robustness, and consistency suggest METABOLIC provides better performance compared to other software and online servers. To highlight the utility and versatility of METABOLIC, we demonstrate its capabilities on diverse metagenomic datasets from the marine subsurface, terrestrial subsurface, meadow soil, deep sea, freshwater lakes, wastewater, and the human gut.

CONCLUSION

METABOLIC enables the consistent and reproducible study of microbial community ecology and biogeochemistry using a foundation of genome-informed microbial metabolism, and will advance the integration of uncultivated organisms into metabolic and biogeochemical models. METABOLIC is written in Perl and R and is freely available under GPLv3 at https://github.com/AnantharamanLab/METABOLIC . Video abstract.

Interactions Between Marine Group II Archaea and Phytoplankton Revealed by Population Correlations in the Northern Coast of South China Sea

Marine Group II (MGII) archaea (Poseidoniales) are the most abundant surface marine planktonic archaea and are widely distributed in both coastal and pelagic waters. The factors affecting their distribution and activity are poorly understood. MGII archaea have the metabolic potential to utilize algae-derived organic matter and are frequently observed in high abundance during or following phytoplankton blooms, suggesting that they are key players of the marine food web. In this study, we studied interactions between MGII archaea and the diverse taxa of phytoplankton in the northern coast of South China Sea. Non-metric multidimensional scaling and cluster analyses demonstrated distinct MGII community patterns in the Pearl River plume (PRP) and the open regions of the northern South China Sea (ONSCS), with MGIIb dominating the former and MGIIa and MGIIb showing remarkable variations in the latter for the same sampling season. Nevertheless, positive correlations (Pearson correlation: R > 0.8 and P < 0.01) in absolute abundances of ribosomal RNA (rRNA)-derived complementary DNA and rRNA genes from network analyses were found between MGII archaea and phytoplankton (cyanobacteria, haptophytes, and stramenopiles in both PRP and ONSCS) among different particle size fractions, indicating their intrinsic relationships under changing environmental conditions. The results of this study may shed light on the multiple interactions between co-existing species in the micro-niches of different oceanic regions.

Classification and Identification of Archaea Using Single-Cell Raman Ejection and Artificial Intelligence: Implications for Investigating Uncultivated Microorganisms

Archaea can produce special cellular components such as polyhydroxyalkanoates, carotenoids, rhodopsin, and ether lipids, which have valuable applications in medicine and green energy production. Most of the archaeal species are uncultivated, posing challenges to investigating their biomarker components and biochemical properties. In this study, we applied Raman spectroscopy to examine the biological characteristics of nine archaeal isolates, including halophilic archaea (Haloferax larsenii, Haloarcula argentinensis, Haloferax mediterranei, Halomicrobium mukohataei, Halomicrobium salinus, Halorussus sp., Natrinema gari), thermophilic archaea (Sulfolobus acidocaldarius), and marine group I (MGI) archaea (Nitrosopumilus maritimus). Linear discriminant analysis of the Raman spectra allowed visualization of significant separations among the nine archaeal isolates. Machine-learning classification models based on support vector machine achieved accuracies of 88-100% when classifying the nine archaeal species. The predicted results were validated by DNA sequencing analysis of cells isolated from the mixture by Raman-activated cell sorting. Raman spectra of uncultured archaea (MGII) were also obtained based on Raman spectroscopy and fluorescence in situ hybridization. The results combining multiple Raman-based techniques indicated that MGII may have the ability to produce lipids distinct from other archaeal species. Our study provides a valuable approach for investigating and classifying archaea, especially uncultured species, at the single-cell level.

Planktonic Archaeal Ether Lipid Origins in Surface Waters of the North Pacific Subtropical Gyre

Thaumarchaeota and Thermoplasmatota are the most abundant planktonic archaea in the sea. Thaumarchaeota contain tetraether lipids as their major membrane lipids, but the lipid composition of uncultured planktonic Thermoplasmatota representatives remains unknown. To address this knowledge gap, we quantified archaeal cells and ether lipids in open ocean depth profiles (0–200 m) of the North Pacific Subtropical Gyre. Planktonic archaeal community structure and ether lipid composition in the water column partitioned into two separate clusters: one above the deep chlorophyll maximum, the other within and below it. In surface waters, Thermoplasmatota densities ranged from 2.11 × 10⁶ to 6.02 × 10⁶ cells/L, while Thaumarchaeota were undetectable. As previously reported for Thaumarchaeota, potential homologs of archaeal tetraether ring synthases were present in planktonic Thermoplasmatota metagenomes. Despite the absence of Thaumarchaeota in surface waters, measurable amounts of intact polar ether lipids were found there. Based on cell abundance estimates, these surface water archaeal ether lipids contributed only 1.21 × 10^–9 ng lipid/Thermoplasmatota cell, about three orders of magnitude less than that reported for Thaumarchaeota cells. While these data indicate that even if some tetraether and diether lipids may be derived from Thermoplasmatota, they would only comprise a small fraction of Thermoplasmatota total biomass. Therefore, while both MGI Thaumarchaeota and MGII/III Thermoplasmatota are potential biological sources of archaeal GDGTs, the Thaumarchaeota appear to be the major contributors of archaeal tetraether lipids in planktonic marine habitats. These results extend and confirm previous reports of planktonic archaeal lipid sources, and further emphasize the need for Thermoplasmatota cultivation, to better characterize the membrane lipid constituents of marine planktonic Thermoplasmatota, and more precisely define the sources and patterns of archaeal tetraether lipid distributions in marine plankton.

Genome-resolved metagenomics using environmental and clinical samples

Recent advances in high-throughput sequencing technologies and computational methods have added a new dimension to metagenomic data analysis i.e. genome-resolved metagenomics. In general terms, it refers to the recovery of draft or high-quality microbial genomes and their taxonomic classification and functional annotation. In recent years, several studies have utilized the genome-resolved metagenome analysis approach and identified previously unknown microbial species from human and environmental metagenomes. In this review, we describe genome-resolved metagenome analysis as a series of four necessary steps: (i) preprocessing of the sequencing reads, (ii) de novo metagenome assembly, (iii) genome binning and (iv) taxonomic and functional analysis of the recovered genomes. For each of these four steps, we discuss the most commonly used tools and the currently available pipelines to guide the scientific community in the recovery and subsequent analyses of genomes from any metagenome sample. Furthermore, we also discuss the tools required for validation of assembly quality as well as for improving quality of the recovered genomes. We also highlight the currently available pipelines that can be used to automate the whole analysis without having advanced bioinformatics knowledge. Finally, we will highlight the most widely adapted and actively maintained tools and pipelines that can be helpful to the scientific community in decision making before they commence the analysis.

Quantifying the Impact of Human Leukocyte Antigen on the Human Gut Microbiota

The composition of the gut microbiota is affected by a number of factors, including the innate and adaptive immune system. The major histocompatibility complex (MHC), or the human leukocyte antigen (HLA) in humans, performs an essential role in vertebrate immunity and is very polymorphic in different populations. HLA determines the specificity of T lymphocyte and natural killer (NK) cell responses, including those against the commensal bacteria present in the human gut. Thus, it is likely that our HLA molecules, and thereby the adaptive immune response, can shape the composition of our microbiota. Here, we investigated the effect of HLA haplotype on the microbiota composition. We performed HLA typing and microbiota composition analyses on 3,002 public human gut microbiome data sets. We found that individuals with functionally similar HLA molecules are also similar in their microbiota composition. Our results show a statistical association between host HLA haplotype and gut microbiota composition. Because the HLA haplotype is a readily measurable parameter of the human immune system, these results open the door to incorporating the genetics of the immune system into predictive microbiome models. IMPORTANCE The microorganisms that live in the digestive tracts of humans, known as the gut microbiota, are essential for hosts' survival, as they support crucial functions. For example, they support the host in facilitating the uptake of nutrients and give colonization resistance against pathogens. The composition of the gut microbiota varies among humans. Studies have proposed multiple factors driving the observed variation, including diet, lifestyle, and health condition. Another major influence on the microbiota is the host's genetic background. We hypothesized the immune system to be one of the most important genetic factors driving the differences observed between gut microbiotas. Therefore, we searched for a link between the polymorphic molecules that shape human immune responses and the composition of the microbiota. HLA molecules are the most polymorphic molecules in our genome and therefore makes an excellent candidate to test such an association. To our knowledge for the first time, our results indicate a significant impact of the HLA on the human gut microbiota.

Expanding Archaeal Diversity and Phylogeny: Past, Present, and Future

The discovery of the Archaea is a major scientific hallmark of the twentieth century. Since then, important features of their cell biology, physiology, ecology, and diversity have been revealed. Over the course of some 40 years, the diversity of known archaea has expanded from 2 to about 30 phyla comprising over 20,000 species. Most of this archaeal diversity has been revealed by environmental 16S rRNA amplicon sequencing surveys using a broad range of universal and targeted primers. Of the few primers that target a large fraction of known archaeal diversity, all display a bias against recently discovered lineages, which limits studies aiming to survey overall archaeal diversity. Induced by genomic exploration of archaeal diversity, and improved phylogenomics approaches, archaeal taxonomic classification has been frequently revised. Due to computational limitations and continued discovery of new lineages, a stable archaeal phylogeny is not yet within reach. Obtaining phylogenetic and taxonomic consensus of archaea should be a high priority for the archaeal research community. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Bacterium Lacking a Known Gene for Retinal Biosynthesis Constructs Functional Rhodopsins

Microbial rhodopsins, comprising a protein moiety (rhodopsin apoprotein) bound to the light-absorbing chromophore retinal, function as ion pumps, ion channels, or light sensors. However, recent genomic and metagenomic surveys showed that some rhodopsin-possessing prokaryotes lack the known genes for retinal biosynthesis. Since rhodopsin apoproteins cannot absorb light energy, rhodopsins produced by prokaryotic strains lacking genes for retinal biosynthesis are hypothesized to be non-functional in cells. In the present study, we investigated whether Aurantimicrobium minutum KNC^T, which is widely distributed in terrestrial environments and lacks any previously identified retinal biosynthesis genes, possesses functional rhodopsin. We initially measured ion transport activity in cultured cells. A light-induced pH change in a cell suspension of rhodopsin-possessing bacteria was detected in the absence of exogenous retinal. Furthermore, spectroscopic analyses of the cell lysate and HPLC-MS/MS analyses revealed that this strain contained an endogenous retinal. These results confirmed that A. minutum KNC^T possesses functional rhodopsin and, hence, produces retinal via an unknown biosynthetic pathway. These results suggest that rhodopsin-possessing prokaryotes lacking known retinal biosynthesis genes also have functional rhodopsins.

Seasonality of archaeal proteorhodopsin and associated Marine Group IIb ecotypes (Ca. Poseidoniales) in the North Western Mediterranean Sea

The Archaea Marine Group II (MGII) is widespread in the world's ocean where it plays an important role in the carbon cycle. Despite recent discoveries on the group's metabolisms, the ecology of this newly proposed order (Candidatus Poseidoniales) remains poorly understood. Here we used a combination of time-series metagenome-assembled genomes (MAGs) and high-frequency 16S rRNA data from the NW Mediterranean Sea to test if the taxonomic diversity within the MGIIb family (Candidatus Thalassarchaeaceae) reflects the presence of different ecotypes. The MAGs' seasonality revealed a MGIIb family composed of different subclades that have distinct lifestyles and physiologies. The vitamin metabolisms were notably different between ecotypes with, in some, a possible link to sunlight's energy. Diverse archaeal proteorhodopsin variants, with unusual signature in key amino acid residues, had distinct seasonal patterns corresponding to changing day length. In addition, we show that in summer, archaea, as opposed to bacteria, disappeared completely from surface waters. Our results shed light on the diversity and the distribution of the euryarchaeotal proteorhodopsin, and highlight that MGIIb is a diverse ecological group. The work shows that time-series based studies of the taxonomy, seasonality, and metabolisms of marine prokaryotes is critical to uncover their diverse role in the ocean.

Discovering multiple types of DNA methylation from bacteria and microbiome using nanopore sequencing

Bacterial DNA methylation occurs at diverse sequence contexts and plays important functional roles in cellular defense and gene regulation. Existing methods for detecting DNA modification from nanopore sequencing data do not effectively support de novo study of unknown bacterial methylomes. In this work, we observed that a nanopore sequencing signal displays complex heterogeneity across methylation events of the same type. To enable nanopore sequencing for broadly applicable methylation discovery, we generated a training dataset from an assortment of bacterial species and developed a method, named nanodisco ( https://github.com/fanglab/nanodisco ), that couples the identification and fine mapping of the three forms of methylation into a multi-label classification framework. We applied it to individual bacteria and the mouse gut microbiome for reliable methylation discovery. In addition, we demonstrated the use of DNA methylation for binning metagenomic contigs, associating mobile genetic elements with their host genomes and identifying misassembled metagenomic contigs.

Comparative Study of the Rhizosphere and Root Endosphere Microbiomes of Cholistan Desert Plants

Microbial communities associated with the rhizosphere and roots of desert halophytes play an important role in plants’ growth and development. Very limited information has been available on the microbial diversity of arid environments of Pakistan. Hence in the current study, the microbial diversity of rhizosphere and root endosphere of desert halophytes, Zygophyllum simplex, Haloxylon salicoricum, Aerva javanica, and Capparis decidua was evaluated. The rhizosphere and root endosphere samples of desert halophytes collected from the three geographic sites of Cholistan desert, Punjab, Pakistan were analyzed by using 16S rRNA based Illumina sequencing. The results showed that Proteobacteria were more abundant in the rhizospheric soils while Actinobacteria were more dominant in the root endosphere of halophytes. Bacteroidetes, Firmicutes, and Deinococcus-Thermus were identified from all rhizospheric soils and roots across the three sites, with variable percentage. Bacillus, Kocuria, Pseudomonas, Halomonas, and Flavobacterium were commonly identified from the rhizosphere and root endosphere of halophytes across all the three sites. At the genus level, microbial diversity from Haloxylon showed the greatest variations between the rhizosphere and root endosphere from the site 2. This study revealed that microbial diversity analysis can be used to study how changes in abiotic factors such as soil moisture content and salinity affect the microbial communities associated with the rhizospheric soils and root endosphere of halophytes across the three sites. This study will also help in the discovery of potential inoculants for crops growing in arid and semi-arid regions of Pakistan.

Metagenomic Data Assembly - The Way of Decoding Unknown Microorganisms

Metagenomics is a segment of conventional microbial genomics dedicated to the sequencing and analysis of combined genomic DNA of entire environmental samples. The most critical step of the metagenomic data analysis is the reconstruction of individual genes and genomes of the microorganisms in the communities using metagenomic assemblers - computational programs that put together small fragments of sequenced DNA generated by sequencing instruments. Here, we describe the challenges of metagenomic assembly, a wide spectrum of applications in which metagenomic assemblies were used to better understand the ecology and evolution of microbial ecosystems, and present one of the most efficient microbial assemblers, SPAdes that was upgraded to become applicable for metagenomics.

Transcriptional activity differentiates families of Marine Group II Euryarchaeota in the coastal ocean

Marine Group II Euryarchaeota (Candidatus Poseidoniales), abundant but yet-uncultivated members of marine microbial communities, are thought to be (photo)heterotrophs that metabolize dissolved organic matter (DOM), such as lipids and peptides. However, little is known about their transcriptional activity. We mapped reads from a metatranscriptomic time series collected at Sapelo Island (GA, USA) to metagenome-assembled genomes to determine the diversity of transcriptionally active Ca. Poseidoniales. Summer metatranscriptomes had the highest abundance of Ca. Poseidoniales transcripts, mostly from the O1 and O3 genera within Ca. Thalassarchaeaceae (MGIIb). In contrast, transcripts from fall and winter samples were predominantly from Ca. Poseidoniaceae (MGIIa). Genes encoding proteorhodopsin, membrane-bound pyrophosphatase, peptidase/proteases, and part of the ß-oxidation pathway were highly transcribed across abundant genera. Highly transcribed genes specific to Ca. Thalassarchaeaceae included xanthine/uracil permease and receptors for amino acid transporters. Enrichment of Ca. Thalassarchaeaceae transcript reads related to protein/peptide, nucleic acid, and amino acid transport and metabolism, as well as transcript depletion during dark incubations, provided further evidence of heterotrophic metabolism. Quantitative PCR analysis of South Atlantic Bight samples indicated consistently abundant Ca. Poseidoniales in nearshore and inshore waters. Together, our data suggest that Ca. Thalassarchaeaceae are important photoheterotrophs potentially linking DOM and nitrogen cycling in coastal waters.

Marine Group-II archaea dominate particle-attached as well as free-living archaeal assemblages in the surface waters of Kongsfjorden, Svalbard, Arctic Ocean

Marine archaea are a significant component of the global oceanic ecosystems, including the polar oceans. However, only a few attempts have been made to study archaea in the high Arctic fjords. Given the importance of Archaea in carbon and nitrogen cycling, it is imperative to explore their diversity and community composition in the high Arctic fjords, such as Kongsfjorden (Svalbard). In the present study, we evaluated archaeal diversity and community composition in the size-fractionated microbial population, viz-a-viz free-living (FL; 0.2-3 μm) and particle-attached (PA; > 3 μm) using archaeal V3-V4 16S rRNA gene amplicon sequencing. Our results indicate that the overall archaeal community in the surface water of Kongsfjorden was dominated by the members of the marine group-II (MGII) archaea, followed by the MGI group members, including Nitrosopumilaceae and Nitrososphaeraceae. Although a clear niche partitioning between PA and FL archaeal communities was not observed, 2 OTUs among 682 OTUs, and 3 ASVs out of 1932 ASVs were differentially abundant among the fractions. OTU001/ASV0002, classified as MGIIa, was differentially abundant in the PA fraction. OTU006/ASV0006/ASV0010 affiliated with MGIIb were differentially abundant in the FL fraction. Particulate organic nitrogen and C:N ratio were the most significant variables (P < 0.05) explaining the observed variation in the FL and PA archaeal communities, respectively. These results indicate an exchange between archaeal communities or a generalist lifestyle switching between FL and PA fractions. Besides, the particles' elemental composition (carbon and nitrogen) seems to play an essential role in shaping the PA archaeal communities in the surface waters of Kongsfjorden.

Exploring Marine Planktonic Archaea: Then and Now

In 1977, Woese and Fox leveraged molecular phylogenetic analyses of ribosomal RNAs and identified a new microbial domain of life on Earth, the Archaebacteria (now known as Archaea). At the time of their discovery, only one archaebacterial group, the strictly anaerobic methanogens, was known. But soon, other phenotypically unrelated microbial isolates were shown to belong to the Archaea, many originating from extreme habitats, including extreme halophiles, extreme thermophiles, and thermoacidophiles. Since most Archaea seemed to inhabit extreme or strictly anoxic habitats, it came as a surprise in 1992 when two new lineages of archaea were reported to be abundant in oxygen rich, temperate marine coastal waters and the deep ocean. Since that time, studies of marine planktonic archaea have revealed many more surprises, including their unexpected ubiquity, unusual symbiotic associations, unpredicted physiologies and biogeochemistry, and global abundance. In this Perspective, early work conducted on marine planktonic Archaea by my lab group and others is discussed in terms of the relevant historical context, some of the original research motivations, and surprises and discoveries encountered along the way.

Biases in genome reconstruction from metagenomic data

BACKGROUND

Advances in sequencing, assembly, and assortment of contigs into species-specific bins has enabled the reconstruction of genomes from metagenomic data (MAGs). Though a powerful technique, it is difficult to determine whether assembly and binning techniques are accurate when applied to environmental metagenomes due to a lack of complete reference genome sequences against which to check the resulting MAGs.

METHODS

We compared MAGs derived from an enrichment culture containing ~20 organisms to complete genome sequences of 10 organisms isolated from the enrichment culture. Factors commonly considered in binning software-nucleotide composition and sequence repetitiveness-were calculated for both the correctly binned and not-binned regions. This direct comparison revealed biases in sequence characteristics and gene content in the not-binned regions. Additionally, the composition of three public data sets representing MAGs reconstructed from the Tara Oceans metagenomic data was compared to a set of representative genomes available through NCBI RefSeq to verify that the biases identified were observable in more complex data sets and using three contemporary binning software packages.

RESULTS

Repeat sequences were frequently not binned in the genome reconstruction processes, as were sequence regions with variant nucleotide composition. Genes encoded on the not-binned regions were strongly biased towards ribosomal RNAs, transfer RNAs, mobile element functions and genes of unknown function. Our results support genome reconstruction as a robust process and suggest that reconstructions determined to be >90% complete are likely to effectively represent organismal function; however, population-level genotypic heterogeneity in natural populations, such as uneven distribution of plasmids, can lead to incorrect inferences.

Distinct Distribution of Archaea From Soil to Freshwater to Estuary: Implications of Archaeal Composition and Function in Different Environments

In addition to inhabiting extreme territories, Archaea are widely distributed in common environments spanning from terrestrial to aquatic environments. This study investigated and compared archaeal community structures from three different habitats (representing distinct environments): agriculture soils (from farming system trials FST, PA, United States), freshwater biofilms (from White Clay Creek, PA, United States), and estuary water (Chesapeake Bay, United States). High-throughput sequencing of 16S rRNA genes indicated that Thaumarchaeota, Euryarchaeota, Nanoarchaeota, Crenarchaeota, and Diapherotrites were the commonly found dominant phyla across these three environments. Similar to Bacteria, distinct community structure and distribution patterns for Archaea were observed in soils vs. freshwater vs. estuary. However, the abundance, richness, evenness, and diversity of archaeal communities were significantly greater in soils than it was in freshwater and estuarine environments. Indicator species (or amplicon sequence variants, ASVs) were identified from different nitrogen and carbon cycling archaeal groups in soils (Nitrososphaerales, Nitrosotaleales, Nitrosopumilales, Methanomassiliicoccales, Lainarchaeales), freshwater biofilms (Methanobacteria, Nitrososphaerales) and Chesapeake Bay (Marine Group II, Nitrosopumilales), suggesting the habitat-specificity of their biogeochemical contributions to different environments. Distinct functional aspects of Archaea were also confirmed by functional predictions (PICRUSt2 analysis). Further, co-occurrence network analysis indicated that only soil Archaea formed stable modules. Keystone species (ASVs) were identified mainly from Methanomassiliicoccales, Nitrososphaerales, Nitrosopumilales. Overall, these results indicate a strong habitat-dependent distribution of Archaea and their functional partitions within the local environments.

Evolutionary Mechanisms of Long-Term Genome Diversification Associated With Niche Partitioning in Marine Picocyanobacteria

Marine picocyanobacteria of the genera Prochlorococcus and Synechococcus are the most abundant photosynthetic organisms on Earth, an ecological success thought to be linked to the differential partitioning of distinct ecotypes into specific ecological niches. However, the underlying processes that governed the diversification of these microorganisms and the appearance of niche-related phenotypic traits are just starting to be elucidated. Here, by comparing 81 genomes, including 34 new Synechococcus, we explored the evolutionary processes that shaped the genomic diversity of picocyanobacteria. Time-calibration of a core-protein tree showed that gene gain/loss occurred at an unexpectedly low rate between the different lineages, with for instance 5.6 genes gained per million years (My) for the major Synechococcus lineage (sub-cluster 5.1), among which only 0.71/My have been fixed in the long term. Gene content comparisons revealed a number of candidates involved in nutrient adaptation, a large proportion of which are located in genomic islands shared between either closely or more distantly related strains, as identified using an original network construction approach. Interestingly, strains representative of the different ecotypes co-occurring in phosphorus-depleted waters (Synechococcus clades III, WPC1, and sub-cluster 5.3) were shown to display different adaptation strategies to this limitation. In contrast, we found few genes potentially involved in adaptation to temperature when comparing cold and warm thermotypes. Indeed, comparison of core protein sequences highlighted variants specific to cold thermotypes, notably involved in carotenoid biosynthesis and the oxidative stress response, revealing that long-term adaptation to thermal niches relies on amino acid substitutions rather than on gene content variation. Altogether, this study not only deciphers the respective roles of gene gains/losses and sequence variation but also uncovers numerous gene candidates likely involved in niche partitioning of two key members of the marine phytoplankton.

Simulation of Enhanced Growth of Marine Group II Euryarchaeota From the Deep Chlorophyll Maximum of the Western Pacific Ocean: Implication for Upwelling Impact on Microbial Functions in the Photic Zone

Mesoscale eddies can have a strong impact on regional biogeochemistry and primary productivity. To investigate the effect of the upwelling of seawater by western Pacific eddies on the composition of the active planktonic marine archaeal community composition of the deep chlorophyll maximum (DCM) layer, mesoscale cold-core eddies were simulated in situ by mixing western Pacific DCM layer water with mesopelagic layer (400 m) water. Illumina sequencing of the 16S rRNA gene and 16S rRNA transcripts indicated that the specific heterotrophic Marine Group IIb (MGIIb) taxonomic group of the DCM layer was rapidly stimulated after receiving fresh substrate from 400 m water, which was dominated by uncultured autotrophic Marine Group I (MGI) archaea. Furthermore, niche differentiation of autotrophic ammonia-oxidizing archaea (MGI) was demonstrated by deep sequencing of 16S rRNA, amoA, and accA genes, respectively. Similar distribution patterns of active Marine Group III (MGIII) were observed in the DCM layer with or without vertical mixing, indicating that they are inclined to utilize the substrates already present in the DCM layer. These findings underscore the importance of mesoscale cyclonic eddies in stimulating microbial processes involved in the regional carbon cycle.

Microbial single-cell omics: the crux of the matter

Abstract

Single-cell genomics and transcriptomics can provide reliable context for assembled genome fragments and gene expression activity on the level of individual prokaryotic genomes. These methods are rapidly emerging as an essential complement to cultivation-based, metagenomics, metatranscriptomics, and microbial community-focused research approaches by allowing direct access to information from individual microorganisms, even from deep-branching phylogenetic groups that currently lack cultured representatives. Their integration and binning with environmental ‘omics data already provides unprecedented insights into microbial diversity and metabolic potential, enabling us to provide information on individual organisms and the structure and dynamics of natural microbial populations in complex environments. This review highlights the pitfalls and recent advances in the field of single-cell omics and its importance in microbiological and biotechnological studies.

Key points

• Single-cell omics expands the tree of life through the discovery of novel organisms, genes, and metabolic pathways.

• Disadvantages of metagenome-assembled genomes are overcome by single-cell omics.

• Functional analysis of single cells explores the heterogeneity of gene expression.

• Technical challenges still limit this field, thus prompting new method developments.

Ecogenomics of the Marine Benthic Filamentous Cyanobacterium Adonisia

Turfs are among the major benthic components of reef systems worldwide. The nearly complete genome sequences, basic physiological characteristics, and phylogenomic reconstruction of two phycobiliprotein-rich filamentous cyanobacteria strains isolated from turf assemblages from the Abrolhos Bank (Brazil) are investigated. Both Adonisia turfae CCMR0081^T (= CBAS 745^T) and CCMR0082 contain approximately 8 Mbp in genome size and experiments identified that both strains exhibit chromatic acclimation. Whereas CCMR0081^T exhibits chromatic acclimation type 3 (CA3) regulating both phycocyanin (PC) and phycoerythrin (PE), CCMR0082 strain exhibits chromatic acclimation type 2 (CA2), in correspondence with genes encoding specific photosensors and regulators for PC and PE. Furthermore, a high number and diversity of secondary metabolite synthesis gene clusters were identified in both genomes, and they were able to grow at high temperatures (28 °C, with scant growth at 30 °C). These characteristics provide insights into their widespread distribution in reef systems.

Combining whole-genome shotgun sequencing and rRNA gene amplicon analyses to improve detection of microbe-microbe interaction networks in plant leaves

Microorganisms from all domains of life establish associations with plants. Although some harm the plant, others antagonize pathogens or prime the plant immune system, support the acquisition of nutrients, tune plant hormone levels, or perform additional services. Most culture-independent plant microbiome research has focused on amplicon sequencing of the 16S rRNA gene and/or the internal transcribed spacer (ITS) of rRNA genomic loci, which show the relative abundance of the microbes to each other. Here, we describe shotgun sequencing of 275 wild Arabidopsis thaliana leaf microbiomes from southwest Germany, with additional bacterial 16S and eukaryotic ITS1 rRNA amplicon data from 176 of these samples. Shotgun data, which unlike the amplicon data capture the ratio of microbe to plant DNA, enable scaling of microbial read abundances to reflect the microbial load on the host. In a more cost-effective hybrid strategy, we show they also allow a similar scaling of amplicon data to overcome compositionality problems. Our wild plants were dominated by bacterial sequences, with eukaryotes contributing only a minority of reads. Microbial membership showed weak associations with both site of origin and plant genotype, both of which were highly confounded in this dataset. There was large variation among microbiomes, with one extreme comprising samples of low complexity and a high load of microorganisms typical of infected plants, and the other extreme being samples of high complexity and a low microbial load. Critically, considering absolute microbial load led to fundamentally different conclusions about microbiome assembly and the interaction networks among major taxa.

Evaluating metagenomics tools for genome binning with real metagenomic datasets and CAMI datasets

Background

Shotgun metagenomics based on untargeted sequencing can explore the taxonomic profile and the function of unknown microorganisms in samples, and complement the shortage of amplicon sequencing. Binning assembled sequences into individual groups, which represent microbial genomes, is the key step and a major challenge in metagenomic research. Both supervised and unsupervised machine learning methods have been employed in binning. Genome binning belonging to unsupervised method clusters contigs into individual genome bins by machine learning methods without the assistance of any reference databases. So far a lot of genome binning tools have emerged. Evaluating these genome tools is of great significance to microbiological research. In this study, we evaluate 15 genome binning tools containing 12 original binning tools and 3 refining binning tools by comparing the performance of these tools on chicken gut metagenomic datasets and the first CAMI challenge datasets.

Results

For chicken gut metagenomic datasets, original genome binner MetaBat, Groopm2 and Autometa performed better than other original binner, and MetaWrap combined the binning results of them generated the most high-quality genome bins. For CAMI datasets, Groopm2 achieved the highest purity (> 0.9) with good completeness (> 0.8), and reconstructed the most high-quality genome bins among original genome binners. Compared with Groopm2, MetaBat2 had similar performance with higher completeness and lower purity. Genome refining binners DASTool predicated the most high-quality genome bins among all genomes binners. Most genome binner performed well for unique strains. Nonetheless, reconstructing common strains still is a substantial challenge for all genome binner.

Conclusions

In conclusion, we tested a set of currently available, state-of-the-art metagenomics hybrid binning tools and provided a guide for selecting tools for metagenomic binning by comparing range of purity, completeness, adjusted rand index, and the number of high-quality reconstructed bins. Furthermore, available information for future binning strategy were concluded.