Unmasking the physiology of mercury detoxifying bacteria from polluted sediments.

Marine sediments polluted from anthropogenic activities can be major reservoirs of toxic mercury species. Some microorganisms in these environments have the capacity to detoxify these pollutants, by using the mer operon. In this study, we characterized microbial cultures isolated from polluted marine sediments growing under diverse environmental conditions of salinity, oxygen availability and mercury tolerance. Specific growth rates and percentage of mercury removal were measured in batch cultures for a selection of isolates. A culture affiliated with Pseudomonas putida (MERCC_1942), which contained a mer operon as well as other genes related to metal resistances, was selected as the best candidate for mercury elimination. In order to optimize mercury detoxification conditions for strain MERCC_1942 in continuous culture, three different dilution rates were tested in bioreactors until the cultures achieved steady state, and they were subsequently exposed to a mercury spike; after 24 h, strain MERCC_1942 removed up to 76% of the total mercury. Moreover, when adapted to high growth rates in bioreactors, this strain exhibited the highest specific mercury detoxification rates. Finally, an immobilization protocol using the sol-gel technology was optimized. These results highlight that some sediment bacteria show capacity to detoxify mercury and could be used for bioremediation applications.

A pan-genomic approach reveals novel Sulfurimonas clade in the ferruginous meromictic Lake Pavin.

The permanently anoxic waters in meromictic lakes create suitable niches for the growth of bacteria using sulphur metabolisms like sulphur oxidation. In Lake Pavin, the anoxic water mass hosts an active cryptic sulphur cycle that interacts narrowly with iron cycling, however the metabolisms of the microorganisms involved are poorly known. Here we combined metagenomics, single-cell genomics, and pan-genomics to further expand our understanding of the bacteria and the corresponding metabolisms involved in sulphur oxidation in this ferruginous sulphide- and sulphate-poor meromictic lake. We highlighted two new species within the genus Sulfurimonas that belong to a novel clade of chemotrophic sulphur oxidisers exclusive to freshwaters. We moreover conclude that this genus holds a key-role not only in limiting sulphide accumulation in the upper part of the anoxic layer but also constraining carbon, phosphate and iron cycling.

Influence of rapid vertical mixing on bacterial community assembly in stratified water columns.

Water column mixing homogenizes thermal and chemical gradients which are known to define distribution of microbial communities and influence the prevailing biogeochemical processes. Little is however known about the effects of rapid water column mixing on the vertical distribution of microbial communities in stratified reservoirs. To address this knowledge gap, physicochemical properties and microbial community composition from 16 S rRNA amplicon sequencing were analyzed before and after mixing of vertically stratified water-column bioreactors. Our results showed that α-diversity of bacterial communities decreased from bottom to surface during periods of thermal stratification. After an experimental mixing event, bacterial community diversity experienced a significant decrease throughout the water column and network connectivity was disrupted, followed by slow recovery. Significant differences in composition were seen for both total (DNA) and active (RNA) bacterial communities when comparing surface and bottom layer during periods of stratification, and when comparing samples collected before mixing and after re-stratification. The dominant predicted community assembly processes for stratified conditions were deterministic while such processes were less important during recovery from episodic mixing. Water quality characteristics of stratified water were significantly correlated with bacterial community diversity and structure. Furthermore, structural equation modeling analyses showed that changes in sulfur may have the greatest direct effect on bacterial community composition. Our results imply that rapid vertical mixing caused by episodic weather extremes and hydrological operations may have a long-term effect on microbial communities and biogeochemical processes.

Revisiting the mercury cycle in marine sediments: A potential multifaceted role for Desulfobacterota.

Marine sediments impacted by urban and industrial pollutants are typically exposed to reducing conditions and represent major reservoirs of toxic mercury species. Mercury methylation mediated by anaerobic microorganisms is favored under such conditions, yet little is known about potential microbial mechanisms for mercury detoxification. We used culture-independent (metagenomics, metabarcoding) and culture-dependent approaches in anoxic marine sediments to identify microbial indicators of mercury pollution and analyze the distribution of genes involved in mercury reduction (merA) and demethylation (merB). While none of the isolates featured merB genes, 52 isolates, predominantly affiliated with Gammaproteobacteria, were merA positive. In contrast, merA genes detected in metagenomes were assigned to different phyla, including Desulfobacterota, Actinomycetota, Gemmatimonadota, Nitrospirota, and Pseudomonadota. This indicates a widespread capacity for mercury reduction in anoxic sediment microbiomes. Notably, merA genes were predominately identified in Desulfobacterota, a phylum previously associated only with mercury methylation. Marker genes involved in the latter process (hgcAB) were also mainly assigned to Desulfobacterota, implying a potential central and multifaceted role of this phylum in the mercury cycle. Network analysis revealed that Desulfobacterota were associated with anaerobic fermenters, methanogens and sulfur-oxidizers, indicating potential interactions between key players of the carbon, sulfur and mercury cycling in anoxic marine sediments.

Shotgun metagenomes from productive lakes in an urban region of Sweden.

Urban lakes provide multiple benefits to society while influencing life quality. Moreover, lakes and their microbiomes are sentinels of anthropogenic impact and can be used for natural resource management and planning. Here, we release original metagenomic data from several well-characterized and anthropogenically impacted eutrophic lakes in the vicinity of Stockholm (Sweden). Our goal was to collect representative microbial community samples and use shotgun sequencing to provide a broad view on microbial diversity of productive urban lakes. Our dataset has an emphasis on Lake Mälaren as a major drinking water reservoir under anthropogenic impact. This dataset includes short-read sequence data and metagenome assemblies from each of 17 samples collected from eutrophic lakes near the greater Stockholm area. We used genome-resolved metagenomics and obtained 2378 metagenome assembled genomes that de-replicated into 514 species representative genomes. This dataset adds new datapoints to previously sequenced lakes and it includes the first sequenced set of metagenomes from Lake Mälaren. Our dataset serves as a baseline for future monitoring of drinking water reservoirs and urban lakes.

Sulfate-reduction and methanogenesis are coupled to Hg(II) and MeHg reduction in rice paddies.

Methylmercury (MeHg) produced in rice paddies is the main source of MeHg accumulation in rice, resulting in high risk of MeHg exposure to humans and wildlife. Net MeHg production is affected by Hg(II) reduction and MeHg demethylation, but it remains unclear to what extent these processes influence net MeHg production, as well as the role of the microbial guilds involved. We used isotopically labeled Hg species and specific microbial inhibitors in microcosm experiments to simultaneously investigate the rates of Hg(II) and MeHg transformations, as well as the key microbial guilds controlling these processes. Results showed that Hg(II) and MeHg reduction rate constants significantly decreased with addition of molybdate or BES, which inhibit sulfate-reduction and methanogenesis, respectively. This suggests that both sulfate-reduction and methanogenesis are important processes controlling Hg(II) and MeHg reduction in rice paddies. Meanwhile, up to 99% of MeHg demethylation was oxidative demethylation (OD) under the incubation conditions, suggesting that OD was the main MeHg degradative pathway in rice paddies. In addition, [202Hg(0)/Me202Hg] from the added 202Hg(NO3)2 was up to 13.9%, suggesting that Hg(II) reduction may constrain Hg(II) methylation in rice paddies at the abandoned Hg mining site. This study improves our understanding of Hg cycling pathways in rice paddies, and more specifically how reduction processes affect net MeHg production and related microbial metabolisms.

Redox gradient shapes the abundance and diversity of mercury-methylating microorganisms along the water column of the Black Sea.

In the global context of seawater deoxygenation triggered by climate change and anthropogenic activities, changes in redox gradients impacting biogeochemical transformations of pollutants, such as mercury, become more likely. Being the largest anoxic basin worldwide, with high concentrations of the potent neurotoxic methylmercury (MeHg), the Black Sea is an ideal natural laboratory to provide new insights about the link between dissolved oxygen concentration and hgcAB gene-carrying (hgc+) microorganisms involved in the formation of MeHg. We combined geochemical and microbial approaches to assess the effect of vertical redox gradients on abundance, diversity, and metabolic potential of hgc+ microorganisms in the Black Sea water column. The abundance of hgcA genes [congruently estimated by quantitative PCR (qPCR) and metagenomics] correlated with MeHg concentration, both maximal in the upper part of the anoxic water. Besides the predominant Desulfobacterales, hgc+ microorganisms belonged to a unique assemblage of diverse-previously underappreciated-anaerobic fermenters from Anaerolineales, Phycisphaerae (characteristic of the anoxic and sulfidic zone), Kiritimatiellales, and Bacteroidales (characteristic of the suboxic zone). The metabolic versatility of Desulfobacterota differed from strict sulfate reduction in the anoxic water to reduction of various electron acceptors in the suboxic water. Linking microbial activity and contaminant concentration in environmental studies is rare due to the complexity of biological pathways. In this study, we disentangle the role of oxygen in shaping the distribution of Hg-methylating microorganisms consistently with MeHg concentration, and we highlight their taxonomic and metabolic niche partitioning across redox gradients, improving the prediction of the response of marine communities to the expansion of oxygen-deficient zones. IMPORTANCE Methylmercury (MeHg) is a neurotoxin detected at high concentrations in certain marine ecosystems, posing a threat to human health. MeHg production is mainly mediated by hgcAB gene-carrying (hgc+) microorganisms. Oxygen is one of the main factors controlling Hg methylation; however, its effect on the diversity and ecology of hgc+ microorganisms remains unknown. Under the current context of seawater deoxygenation, mercury cycling is expected to be disturbed. Here, we show the strong effect of oxygen gradients on the distribution of potential Hg methylators. In addition, we show for the first time the significant contribution of a unique assemblage of potential fermenters from Anaerolineales, Phycisphaerae, and Kiritimatiellales to Hg methylation, stratified in different redox niches along the Black Sea gradient. Our results considerably expand the known taxonomic diversity and ecological niches prone to the formation of MeHg and contribute to better apprehend the consequences of oxygen depletion in seawater.

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.

Riverine microbial communities impacted by per- and polyfluoroalkyl substances (PFAS) emissions from a fluoropolymer-manufacturing plant.

Per- and polyfluoroalkyl substances (PFAS) are widespread pollutants that can influence microorganisms. To unveil the effects of PFAS in natural microecosystems, a study that focused on the bacterial, fungal, and microeukaryotic communities around the PFAS point source was conducted in China. A total of 255 specific taxa were significantly different between the upstream and downstream samples, 54 of which were directly correlated with PFAS concentration. Stenotrophomonas (99.2 %), Ralstonia (90.7 %), Phoma (21.9 %), and Alternaria (97.6 %) were the dominant genera in sediment samples from the downstream communities. In addition, most of the dominant taxa were significantly correlated with PFAS concentration. Furthermore, the type of microorganism (bacteria, fungi, and microeukaryotes) and habitat (sediment or pelagic) also influence the microbial community responses to PFAS exposure. Pelagic microorganisms featured more PFAS-correlated biomarker taxa (36 pelagic microeukaryotic biomarkers and 8 pelagic bacteria biomarkers) than the sediments (9 sediment fungi biomarkers and 5 sediment bacteria biomarker). In general, around the factory, the microbial community was more variable in pelagic, summer, and microeukaryotic conditions than in other types. Attention needs to be paid to these variables in the future effect of PFAS on microorganisms.

Thiobacillus as a key player for biofilm formation in oligotrophic groundwaters of the Fennoscandian Shield.

Biofilm formation is a common adaptation for microbes in energy-limited conditions such as those prevalent in the vast deep terrestrial biosphere. However, due to the low biomass and the inaccessible nature of subsurface groundwaters, the microbial populations and genes involved in its formation are understudied. Here, a flow-cell system was designed to investigate biofilm formation under in situ conditions in two groundwaters of contrasting age and geochemistry at the Äspö Hard Rock Laboratory, Sweden. Metatranscriptomes showed Thiobacillus, Sideroxydans, and Desulforegula to be abundant and together accounted for 31% of the transcripts in the biofilm communities. Differential expression analysis highlighted Thiobacillus to have a principal role in biofilm formation in these oligotrophic groundwaters by being involved in relevant processes such as the formation of extracellular matrix, quorum sensing, and cell motility. The findings revealed an active biofilm community with sulfur cycling as a prominent mode of energy conservation in the deep biosphere.

Large-scale phylogenomics of aquatic bacteria reveal molecular mechanisms for adaptation to salinity.

The crossing of environmental barriers poses major adaptive challenges. Rareness of freshwater-marine transitions separates the bacterial communities, but how these are related to brackish counterparts remains elusive, as do the molecular adaptations facilitating cross-biome transitions. We conducted large-scale phylogenomic analysis of freshwater, brackish, and marine quality-filtered metagenome-assembled genomes (11,248). Average nucleotide identity analyses showed that bacterial species rarely existed in multiple biomes. In contrast, distinct brackish basins cohosted numerous species, but their intraspecific population structures displayed clear signs of geographic separation. We further identified the most recent cross-biome transitions, which were rare, ancient, and most commonly directed toward the brackish biome. Transitions were accompanied by systematic changes in amino acid composition and isoelectric point distributions of inferred proteomes, which evolved over millions of years, as well as convergent gains or losses of specific gene functions. Therefore, adaptive challenges entailing proteome reorganization and specific changes in gene content constrains the cross-biome transitions, resulting in species-level separation between aquatic biomes.

Out of the blue: the independent activity of sulfur-oxidizers and diatoms mediate the sudden color shift of a tropical river.

BACKGROUND: Río Celeste ("Sky-Blue River") is a river located in the Tenorio National Park (Costa Rica) that has become an important hotspot for eco-tourism due to its striking sky-blue color. A previous study indicated that this color is not caused by dissolved chemical species, but by formation of light-scattering aluminosilicate particles at the mixing point of two colorless streams, the acidic Quebrada Agria and the neutral Río Buenavista. RESULTS: We now present microbiological information on Río Celeste and its two tributaries, as well as a more detailed characterization of the particles that occur at the mixing point. Our results overturn the previous belief that the light scattering particles are formed by the aggregation of smaller particles coming from Río Buenavista, and rather point to chemical formation of hydroxyaluminosilicate colloids when Quebrada Agria is partially neutralized by Río Buenavista, which also contributes silica to the reaction. The process is mediated by the activities of different microorganisms in both streams. In Quebrada Agria, sulfur-oxidizing bacteria generate an acidic environment, which in turn cause dissolution and mobilization of aluminum and other metals. In Río Buenavista, the growth of diatoms transforms dissolved silicon into colloidal biogenic forms which may facilitate particle precipitation. CONCLUSIONS: We show how the sky-blue color of Río Celeste arises from the tight interaction between chemical and biological processes, in what constitutes a textbook example of emergent behavior in environmental microbiology.

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

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

Microcosm experiment to test bacterial responses to perfluorooctanoate exposure.

The impact of perfluoroalkyl and polyfluoroalkyl substances on microbial communities is challenging to investigate in situ because of the complexity and dynamics of natural ecosystems. In the present study, four microcosms were established to explore the impact of perfluorooctanoate (PFOA) on bacterial communities in riverine and marine settings. PFOA distribution between the aqueous and sedimentary phases fluctuated in both PFOA-amended and unamended control systems. PFOA was more rapidly partitioned into the sediment in marine than in riverine microcosms. Differences in iron concentration and salinity may influence PFOA exchange between water and sediment. In marine microcosms, the alpha diversity of bacterial communities was significantly correlated to PFOA concentration. PFOA tended to correlate more strongly with bacterial community composition in water than in sediment. At the whole system level, Lefse's analysis indicated Algoriphagus halophilus as biomarkers for PFOA exposure in both riverine and marine systems, and the family Flavobacteriaceae were also more abundant in the exposed systems. In terms of temporal variation (comparison between three time points in the systems), metastat analysis showed great variability of potential PFOA-sensitive bacteria at the genus level. As such, most PFOA-sensitive genera were transitory and variable and existed for a short term in different systems (river, sea, blank, and experiment) and phases. Compared with other PFOA-sensitive genera, we suggest that further research is carried out to explore the use of Limnobacter as a bioindicator for temporal monitoring of PFOA pollution.

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.

Multiomics in the central Arctic Ocean for benchmarking biodiversity change.

Multiomics approaches need to be applied in the central Arctic Ocean to benchmark biodiversity change and to identify novel species and their genes. As part of MOSAiC, EcoOmics will therefore be essential for conservation and sustainable bioprospecting in one of the least explored ecosystems on Earth.

Genome-resolved analyses show an extensive diversification in key aerobic hydrocarbon-degrading enzymes across bacteria and archaea.

BACKGROUND: Hydrocarbons (HCs) are organic compounds composed solely of carbon and hydrogen that are mainly accumulated in oil reservoirs. As the introduction of all classes of hydrocarbons including crude oil and oil products into the environment has increased significantly, oil pollution has become a global ecological problem. However, our perception of pathways for biotic degradation of major HCs and key enzymes in these bioconversion processes has mainly been based on cultured microbes and is biased by uneven taxonomic representation. Here we used Annotree to provide a gene-centric view of the aerobic degradation ability of aliphatic and aromatic HCs in 23,446 genomes from 123 bacterial and 14 archaeal phyla.  RESULTS: Apart from the widespread genetic potential for HC degradation in Proteobacteria, Actinobacteriota, Bacteroidota, and Firmicutes, genomes from an additional 18 bacterial and 3 archaeal phyla also hosted key HC degrading enzymes. Among these, such degradation potential has not been previously reported for representatives in the phyla UBA8248, Tectomicrobia, SAR324, and Eremiobacterota. Genomes containing whole pathways for complete degradation of HCs were only detected in Proteobacteria and Actinobacteriota. Except for several members of Crenarchaeota, Halobacterota, and Nanoarchaeota that have tmoA, ladA, and alkB/M key genes, respectively, representatives of archaeal genomes made a small contribution to HC degradation. None of the screened archaeal genomes coded for complete HC degradation pathways studied here; however, they contribute significantly to peripheral routes of HC degradation with bacteria. CONCLUSION: Phylogeny reconstruction showed that the reservoir of key aerobic hydrocarbon-degrading enzymes in Bacteria and Archaea undergoes extensive diversification via gene duplication and horizontal gene transfer. This diversification could potentially enable microbes to rapidly adapt to novel and manufactured HCs that reach the environment.

Regional factors as major drivers for microbial community turnover in tropical cascading reservoirs.

The turnover of microbial communities across space is dictated by local and regional factors. Locally, selection shapes community assembly through biological interactions between organisms and the environment, while regional factors influence microbial dispersion patterns. Methods used to disentangle the effects of local and regional factors typically do not aim to identify ecological processes underlying the turnover. In this paper, we identified and quantified these processes for three operational microbial subcommunities (cyanobacteria, particle-attached, and free-living bacteria) from a tropical cascade of freshwater reservoirs with decreasing productivity, over two markedly different dry and rainy seasons. We hypothesized that during the dry season communities would mainly be controlled by selection shaped by the higher environmental heterogeneity that results from low hydrological flow and connectivity between reservoirs. We expected highly similar communities shaped by dispersal and a more homogenized environment during the rainy season, enhanced by increased flow rates. Even if metacommunities were largely controlled by regional events in both periods, the selection had more influence on free-living communities during the dry period, possibly related to elevated dissolved organic carbon concentration, while drift as a purely stochastic factor, had more influence on cyanobacterial communities. Each subcommunity had distinct patterns of turnover along the cascade related to diversity (Cyanobacteria), lifestyle and size (Free-living), and spatial dynamics (particle-attached).

Expression Levels of hgcAB Genes and Mercury Availability Jointly Explain Methylmercury Formation in Stratified Brackish Waters.

Neurotoxic methylmercury (MeHg) is formed by microbial methylation of inorganic divalent Hg (HgII) and constitutes severe environmental and human health risks. The methylation is enabled by hgcA and hgcB genes, but it is not known if the associated molecular-level processes are rate-limiting or enable accurate prediction of MeHg formation in nature. In this study, we investigated the relationships between hgc genes and MeHg across redox-stratified water columns in the brackish Baltic Sea. We showed, for the first time, that hgc transcript abundance and the concentration of dissolved HgII-sulfide species were strong predictors of both the HgII methylation rate and MeHg concentration, implying their roles as principal joint drivers of MeHg formation in these systems. Additionally, we characterized the metabolic capacities of hgc+ microorganisms by reconstructing their genomes from metagenomes (i.e., hgc+ MAGs), which highlighted the versatility of putative HgII methylators in the water column of the Baltic Sea. In establishing relationships between hgc transcripts and the HgII methylation rate, we advance the fundamental understanding of mechanistic principles governing MeHg formation in nature and enable refined predictions of MeHg levels in coastal seas in response to the accelerating spread of oxygen-deficient zones.

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.