mOTUpan: a robust Bayesian approach to leverage metagenome-assembled genomes for core-genome estimation

Recent advances in sequencing and bioinformatics have expanded the tree of life by providing genomes for uncultured environmentally relevant clades, either through metagenome-assembled genomes or through single-cell genomes. While this expanded diversity can provide novel insights into microbial population structure, most tools available for core-genome estimation are sensitive to genome completeness. Consequently, a major portion of the huge phylogenetic diversity uncovered by environmental genomic approaches remains excluded from such analyses. We present mOTUpan, a novel iterative Bayesian method for computing the core genome for sets of genomes of highly diverse completeness range. The likelihood for each gene cluster to belong to core or accessory genome is estimated by computing the probability of its presence/absence pattern in the target genome set. The core-genome prediction is computationally efficient and can be scaled up to thousands of genomes. It has shown comparable estimates to state-of-the-art tools Roary and PPanGGOLiN for high-quality genomes and is capable of using genomes at lower completeness thresholds. mOTUpan wraps a bootstrapping procedure to estimate the quality of a specific core-genome prediction, as the accuracy of each run will depend on the specific completeness distribution and the number of genomes in the dataset under scrutiny. mOTUpan is implemented in the mOTUlizer software package, and available at github.com/moritzbuck/mOTUlizer, under GPL 3.0 license

Exploring environmental intra-species diversity through non-redundant pangenome assemblies

At the genome level, microorganisms are highly adaptable both in terms of allele and gene composition. Such heritable traits emerge in response to different environmental niches and can have a profound influence on microbial community dynamics. As a consequence, any individual genome or population will contain merely a fraction of the total genetic diversity of any operationally defined "species", whose ecological potential can thus be only fully understood by studying all of their genomes and the genes therein. This concept, known as the pangenome, is valuable for studying microbial ecology and evolution, as it partitions genomes into core (present in all the genomes from a species, and responsible for housekeeping and species-level niche adaptation among others) and accessory regions (present only in some, and responsible for intra-species differentiation). Here we present SuperPang, an algorithm producing pangenome assemblies from a set of input genomes of varying quality, including metagenome-assembled genomes (MAGs). SuperPang runs in linear time and its results are complete, non-redundant, preserve gene ordering and contain both coding and non-coding regions. Our approach provides a modular view of the pangenome, identifying operons and genomic islands, and allowing to track their prevalence in different populations. We illustrate this by analysing intra-species diversity in Polynucleobacter, a bacterial genus ubiquitous in freshwater ecosystems, characterized by their streamlined genomes and their ecological versatility. We show how SuperPang facilitates the simultaneous analysis of allelic and gene content variation under different environmental pressures, allowing us to study the drivers of microbial diversification at unprecedented resolution

16S rRNA gene sequences of Candidatus Methylumidiphilus (Methylococcales), a putative methanotrophic genus in lakes and ponds

acceptedVersionPeer reviewe

Streamlined and Abundant Bacterioplankton Thrive in Functional Cohorts

While fastidious microbes can be abundant and ubiquitous in their natural communities, many fail to grow axenically in laboratories due to auxotrophies or other dependencies. To overcome auxotrophies, these microbes rely on their surrounding cohort. A cohort may consist of kin (ecotypes) or more distantly related organisms (community) with the cooperation being reciprocal or nonreciprocal and expensive (Black Queen hypothesis) or costless (by-product). These metabolic partnerships (whether at single species population or community level) enable dominance by and coexistence of these lineages in nature. Here we examine the relevance of these cooperation models to explain the abundance and ubiquity of the dominant fastidious bacterioplankton of a dimictic mesotrophic freshwater lake. Using both culture-dependent (dilution mixed cultures) and culture-independent (small subunit [SSU] rRNA gene time series and environmental metagenomics) methods, we independently identified the primary cohorts of actinobacterial genera "Candidatus Planktophila" (acI-A) and "Candidatus Nanopelagicus" (acI-B) and the proteobacterial genus "Candidatus Fonsibacter" (LD12). While "Ca. Planktophila" and "Ca. Fonsibacter" had no correlation in their natural habitat, they have the potential to be complementary in laboratory settings. We also investigated the bifunctional catalase-peroxidase enzyme KatG (a common good which "Ca. Planktophila" is dependent upon) and its most likely providers in the lake. Further, we found that while ecotype and community cooperation combined may explain "Ca. Planktophila" population abundance, the success of "Ca. Nanopelagicus" and "Ca. Fonsibacter" is better explained as a community by-product. Ecotype differentiation of "Ca. Fonsibacter" as a means of escaping predation was supported but not for overcoming auxotrophies.IMPORTANCE This study examines evolutionary and ecological relationships of three of the most ubiquitous and abundant freshwater bacterial genera: "Ca. Planktophila" (acI-A), "Ca. Nanopelagicus" (acI-B), and "Ca. Fonsibacter" (LD12). Due to high abundance, these genera might have a significant influence on nutrient cycling in freshwaters worldwide, and this study adds a layer of understanding to how seemingly competing clades of bacteria can coexist by having different cooperation strategies. Our synthesis ties together network and ecological theory with empirical evidence and lays out a framework for how the functioning of populations within complex microbial communities can be studied

Candidatus Methylumidiphilus Drives Peaks in Methanotrophic Relative Abundance in Stratified Lakes and Ponds Across Northern Landscapes

Boreal lakes and ponds produce two-thirds of the total natural methane emissions above the latitude of 50 degrees North. These lake emissions are regulated by methanotrophs which can oxidize up to 99% of the methane produced in the sediments and the water column. Despite their importance, the diversity and distribution of the methanotrophs in lakes are still poorly understood. Here, we used shotgun metagenomic data to explore the diversity and distribution of methanotrophs in 40 oxygen-stratified water bodies in boreal and subarctic areas in Europe and North America. In our data, gammaproteobacterial methanotrophs (order Methylococcales) generally dominated the methanotrophic communities throughout the water columns. A recently discovered lineage of Methylococcales, Candidatus Methylumidiphilus, was present in all the studied water bodies and dominated the methanotrophic community in lakes with a high relative abundance of methanotrophs. Alphaproteobacterial methanotrophs were the second most abundant group of methanotrophs. In the top layer of the lakes, characterized by low CH4 concentration, their abundance could surpass that of the gammaproteobacterial methanotrophs. These results support the theory that the alphaproteobacterial methanotrophs have a high affinity for CH4 and can be considered stress-tolerant strategists. In contrast, the gammaproteobacterial methanotrophs are competitive strategists. In addition, relative abundances of anaerobic methanotrophs, Candidatus Methanoperedenaceae and Candidatus Methylomirabilis, were strongly correlated, suggesting possible co-metabolism. Our data also suggest that these anaerobic methanotrophs could be active even in the oxic layers. In non-metric multidimensional scaling, alpha- and gammaproteobacterial methanotrophs formed separate clusters based on their abundances in the samples, except for the gammaproteobacterial Candidatus Methylumidiphilus, which was separated from these two clusters. This may reflect similarities in the niche and environmental requirements of the different genera within alpha- and gammaproteobacterial methanotrophs. Our study confirms the importance of O-2 and CH4 in shaping the methanotrophic communities and suggests that one variable cannot explain the diversity and distribution of the methanotrophs across lakes. Instead, we suggest that the diversity and distribution of freshwater methanotrophs are regulated by lake-specific factors

Yersinia Virulence Factor YopM Induces Sustained RSK Activation by Interfering with Dephosphorylation

Background: Pathogenic yersiniae inject several effector proteins (Yops) into host cells, which subverts immune functions and enables the bacteria to survive within the host organism. YopM, whose deletion in enteropathogenic yersiniae results in a dramatic loss of virulence, has previously been shown to form a complex with and activate the multifunctional kinases PKN2 and RSK1 in transfected cells. Methodology/Principal Findings: In a near physiological approach with double-affinity-tagged YopM being translocated into the macrophage cell line J774A.1 via the natural type three secretion system of Yersinia we verified the interaction of YopM with PKN2 and RSK1 and detected association with additional PKN and RSK isoforms. In transfected and infected cells YopM induced sustained phosphorylation of RSK at its activation sites serine-380 and serine-221 even in the absence of signalling from its upstream kinase ERK1/2, suggesting inhibition of dephosphorylation. ATP-depletion and in vitro assays using purified components directly confirmed that YopM shields RSK isoforms from phosphatase activity towards serines 380 and 221. Conclusions/Significance: Our study suggests that during Yersinia infection YopM induces sustained activation of RSK by blocking dephosphorylation of its activatory phosphorylation sites. This may represent a novel mode of action of a bacterial virulence factor

A Genomic Perspective Across Earth's Microbiomes Reveals That Genome Size in Archaea and Bacteria Is Linked to Ecosystem Type and Trophic Strategy

Our view of genome size in Archaea and Bacteria has remained skewed as the data has been dominated by genomes of microorganisms that have been cultivated under laboratory settings. However, the continuous effort to catalog Earth's microbiomes, specifically propelled by recent extensive work on uncultivated microorganisms, provides an opportunity to revise our perspective on genome size distribution. We present a meta-analysis that includes 26,101 representative genomes from 3 published genomic databases; metagenomic assembled genomes (MAGs) from GEMs and stratfreshDB, and isolates from GTDB. Aquatic and host-associated microbial genomes present on average the smallest estimated genome sizes (3.1 and 3.0 Mbp, respectively). These are followed by terrestrial microbial genomes (average 3.7 Mbp), and genomes from isolated microorganisms (average 4.3 Mbp). On the one hand, aquatic and host-associated ecosystems present smaller genomes sizes in genera of phyla with genome sizes above 3 Mbp. On the other hand, estimated genome size in phyla with genomes under 3 Mbp showed no difference between ecosystems. Moreover, we observed that when using 95% average nucleotide identity (ANI) as an estimator for genetic units, only 3% of MAGs cluster together with genomes from isolated microorganisms. Although there are potential methodological limitations when assembling and binning MAGs, we found that in genome clusters containing both environmental MAGs and isolate genomes, MAGs were estimated only an average 3.7% smaller than isolate genomes. Even when assembly and binning methods introduce biases, estimated genome size of MAGs and isolates are very similar. Finally, to better understand the ecological drivers of genome size, we discuss on the known and the overlooked factors that influence genome size in different ecosystems, phylogenetic groups, and trophic strategies

Microbial communities mediating net methylmercury formation along a trophic gradient in a peatland chronosequence

Peatlands are generally important sources of methylmercury (MeHg) to adjacent aquatic ecosystems, increasing the risk of human and wildlife exposure to this highly toxic compound. While microorganisms play important roles in mercury (Hg) geochemical cycles where they directly and indirectly affect MeHg formation in peatlands, potential linkages between net MeHg formation and microbial communities involving these microorganisms remain unclear. To address this gap, microbial community composition and specific marker gene transcripts were investigated along a trophic gradient in a geographically constrained peatland chronosequence. Our results showed a clear spatial pattern in microbial community composition along the gradient that was highly driven by peat soil properties and significantly associated with net MeHg formation as approximated by MeHg concentration and %MeHg of total Hg concentration. Known fermentative, syntrophic, methanogenic and iron-reducing metabolic guilds had the strong positive correlations to net MeHg formation, while methanotrophic and methylotrophic microorganisms were negatively correlated. Our results indicated that sulfate reducers did not have a key role in net MeHg formation. Microbial activity as interpreted from 16S rRNA sequences was significantly correlated with MeHg and %MeHg. Our findings shed new light on the role of microbial community in net MeHg formation of peatlands that undergo ontogenetic change

Energy efficiency and biological interactions define the core microbiome of deep oligotrophic groundwater

The work conducted by the U.S. Department of Energy Joint Genome Institute, a DOE Office of Science User Facility, is supported under Contract No. DE-AC02-05CH11231. The Swedish Research Council (contracts 2018-04311, 2017-04422, and 2014-4398) and The Swedish Nuclear Fuel and Waste Management Company (SKB) supported the study. M.D. thanks the Crafoord Foundation (contracts 20180599 and 20130557), the Nova Center for University Studies, Research and Development, and Familjen Hellmans Stiftelse for financial support. M.D. and D.S. thank the Carl Tryggers Foundation (grant KF16: 18) for financial support. S.B. and M.M. acknowledge financial support from the Swedish Research Council and Science for Life Laboratory. High-throughput sequencing was also carried out at the National Genomics Infrastructure hosted by the Science for Life Laboratory. Bioinformatics analyses were carried out utilizing the Uppsala Multidisciplinary Center for Advanced Computational Science (UPPMAX) at Uppsala University (projects b2013127, SNIC 2019/3-22, and SNIC 2020/5-19) with support from a SciLifeLab-WABI bioinformatics grant. We would also like to thank Mats Astrom for his comments on the Aspo HRL lithology. JS is financially supported by the Knut and Alice Wallenberg Foundation as part of the National Bioinformatics Infrastructure Sweden at SciLifeLab.While oligotrophic deep groundwaters host active microbes attuned to the low-end of the bioenergetics spectrum, the ecological constraints on microbial niches in these ecosystems and their consequences for microbiome convergence are unknown. Here, we provide a genome-resolved, integrated omics analysis comparing archaeal and bacterial communities in disconnected fracture fluids of the Fennoscandian Shield in Europe. Leveraging a dataset that combines metagenomes, single cell genomes, and metatranscriptomes, we show that groundwaters flowing in similar lithologies offer fixed niches that are occupied by a common core microbiome. Functional expression analysis highlights that these deep groundwater ecosystems foster diverse, yet cooperative communities adapted to this setting. We suggest that these communities stimulate cooperation by expression of functions related to ecological traits, such as aggregate or biofilm formation, while alleviating the burden on microorganisms producing compounds or functions that provide a collective benefit by facilitating reciprocal promiscuous metabolic partnerships with other members of the community. We hypothesize that an episodic lifestyle enabled by reversible bacteriostatic functions ensures the subsistence of the oligotrophic deep groundwater microbiome.United States Department of Energy (DOE) DE-AC02-05CH11231Swedish Research CouncilEuropean Commission 2018-04311 2017-04422 2014-4398Swedish Nuclear Fuel and Waste Management Company (SKB)Crafoord Foundation 20180599 20130557Nova Center for University Studies, Research and DevelopmentFamiljen Hellmans StiftelseCarl Tryggers Foundation KF16: 18Swedish Research Council European CommissionScience for Life LaboratoryUppsala Multidisciplinary Center for Advanced Computational Science (UPPMAX) at Uppsala University b2013127 SNIC 2019/3-22 SNIC 2020/5-19SciLifeLab-WABI bioinformatics grantKnut and Alice Wallenberg Foundation as part of the National Bioinformatics Infrastructure Sweden at SciLifeLa

Deltaproteobacteria and Spirochaetes-Like Bacteria Are Abundant Putative Mercury Methylators in Oxygen-Deficient Water and Marine Particles in the Baltic Sea

Methylmercury (MeHg), a neurotoxic compound biomagnifying in aquatic food webs, can be a threat to human health via fish consumption. However, the composition and distribution of the microbial communities mediating the methylation of mercury (Hg) to MeHg in marine systems remain largely unknown. In order to fill this knowledge gap, we used the Baltic Sea Reference Metagenome (BARM) dataset to study the abundance and distribution of the genes involved in Hg methylation (thehgcABgene cluster). We determined the relative abundance of thehgcABgenes and their taxonomic identity in 81 brackish metagenomes that cover spatial, seasonal and redox variability in the Baltic Sea water column. ThehgcABgenes were predominantly detected in anoxic water, but somehgcABgenes were also detected in hypoxic and normoxic waters. Phylogenetic analysis identified putative Hg methylators within Deltaproteobacteria, in oxygen-deficient water layers, but also Spirochaetes-like and Kiritimatiellaeota-like bacteria. Higher relative quantities ofhgcABgenes were found in metagenomes from marine particles compared to free-living communities in anoxic water, suggesting that such particles are hotspot habitats for Hg methylators in oxygen-depleted seawater. Altogether, our work unveils the diversity of the microorganisms with the potential to mediate MeHg production in the Baltic Sea and pinpoint the important ecological niches for these microorganisms within the marine water column