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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08906-w.

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 (Hg^II) 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 Hg^II-sulfide species were strong predictors of both the Hg^II 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 Hg^II methylators in the water column of the Baltic Sea. In establishing relationships between hgc transcripts and the Hg^II 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.

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).

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

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

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.

Development of gut microbiota during the first 2 years of life

Although development of microbiota in childhood has been linked to chronic immune-related conditions, early childhood determinants of microbiota development have not been fully elucidated. We used 16S rRNA sequencing to analyse faecal and saliva samples from 83 children at four time-points during their first 2 years of life and from their mothers. Our findings confirm that gut microbiota in infants have low diversity and highlight that some properties are shared with the oral microbiota, although inter-individual differences are present. A considerable convergence in gut microbiota composition was noted across the first 2 years of life, towards a more diverse adult-like microbiota. Mode of delivery accounted for some of the inter-individual variation in early childhood, but with a pronounced attenuation over time. Our study extends previous research with further characterization of the major shift in gut microbiota composition during the first 2 years of life.

Norway spruce postglacial recolonization of Fennoscandia

Contrasting theories exist regarding how Norway spruce (Picea abies) recolonized Fennoscandia after the last glaciation and both early Holocene establishments from western microrefugia and late Holocene colonization from the east have been postulated. Here, we show that Norway spruce was present in southern Fennoscandia as early as 14.7 ± 0.1 cal. kyr BP and that the millennia-old clonal spruce trees present today in central Sweden likely arrived with an early Holocene migration from the east. Our findings are based on ancient sedimentary DNA from multiple European sites (N = 15) combined with nuclear and mitochondrial DNA analysis of ancient clonal (N = 135) and contemporary spruce forest trees (N = 129) from central Sweden. Our other findings imply that Norway spruce was present shortly after deglaciation at the margins of the Scandinavian Ice Sheet, and support previously disputed finds of pollen in southern Sweden claiming spruce establishment during the Lateglacial.

Contrasting theories exist regarding how Norway spruce recolonized Fennoscandia after the last glaciation. Here, the authors provide evidences from sedimentary ancient DNA and modern population genomics to support that Norway spruce was present in southern Fennoscandia shortly after deglaciation and the early Holocene migration from the east.

Environmental paleomicrobiology: using DNA preserved in aquatic sediments to its full potential

In‐depth knowledge about spatial and temporal variation in microbial diversity and function is needed for a better understanding of ecological and evolutionary responses to global change. In particular, the study of microbial ancient DNA preserved in sediment archives from lakes and oceans can help us to evaluate the responses of aquatic microbes in the past and make predictions about future biodiversity change in those ecosystems. Recent advances in molecular genetic methods applied to the analysis of historically deposited DNA in sediments have not only allowed the taxonomic identification of past aquatic microbial communities but also enabled tracing their evolution and adaptation to episodic disturbances and gradual environmental change. Nevertheless, some challenges remain for scientists to take full advantage of the rapidly developing field of paleo‐genetics, including the limited ability to detect rare taxa and reconstruct complete genomes for evolutionary studies. Here, we provide a brief review of some of the recent advances in the field of environmental paleomicrobiology and discuss remaining challenges related to the application of molecular genetic methods to study microbial diversity, ecology, and evolution in sediment archives. We anticipate that, in the near future, environmental paleomicrobiology will shed new light on the processes of microbial genome evolution and microbial ecosystem responses to quaternary environmental changes at an unprecedented level of detail. This information can, for example, aid geological reconstructions of biogeochemical cycles and predict ecosystem responses to environmental perturbations, including in the context of human‐induced global changes.

Connectivity of Fennoscandian Shield terrestrial deep biosphere microbiomes with surface communities

The deep biosphere is an energy constrained ecosystem yet fosters diverse microbial communities that are key in biogeochemical cycling. Whether microbial communities in deep biosphere groundwaters are shaped by infiltration of allochthonous surface microorganisms or the evolution of autochthonous species remains unresolved. In this study, 16S rRNA gene amplicon analyses showed that few groups of surface microbes infiltrated deep biosphere groundwaters at the Äspö Hard Rock Laboratory, Sweden, but that such populations constituted up to 49% of the microbial abundance. The dominant persisting phyla included Patescibacteria, Proteobacteria, and Epsilonbacteraeota. Despite the hydrological connection of the Baltic Sea with the studied groundwaters, infiltrating microbes predominantly originated from deep soil groundwater. Most deep biosphere groundwater populations lacked surface representatives, suggesting that they have evolved from ancient autochthonous populations. We propose that deep biosphere groundwater communities in the Fennoscandian Shield consist of selected infiltrated and indigenous populations adapted to the prevailing conditions.

Westmeijer et al. employ high-throughput sequencing to investigate the connection between deep biosphere groundwaters and surface microbial communities. They suggest that the microbial communities of deep biosphere groundwaters in the Fennoscandian Shield are mostly comprised of autochthonous species, rather than migratory surface representatives.

Distinct B cell subsets in Peyer's patches convey probiotic effects by Limosilactobacillus reuteri

BACKGROUND

Intestinal Peyer's patches (PPs) form unique niches for bacteria-immune cell interactions that direct host immunity and shape the microbiome. Here we investigate how peroral administration of probiotic bacterium Limosilactobacillus reuteri R2LC affects B lymphocytes and IgA induction in the PPs, as well as the downstream consequences on intestinal microbiota and susceptibility to inflammation.

RESULTS

The B cells of PPs were separated by size to circumvent activation-dependent cell identification biases due to dynamic expression of markers, which resulted in two phenotypically, transcriptionally, and spatially distinct subsets: small IgD⁺/GL7^-/S1PR1⁺/Bcl6, CCR6-expressing pre-germinal center (GC)-like B cells with innate-like functions located subepithelially, and large GL7⁺/S1PR1^-/Ki67⁺/Bcl6, CD69-expressing B cells with strong metabolic activity found in the GC. Peroral L. reuteri administration expanded both B cell subsets and enhanced the innate-like properties of pre-GC-like B cells while retaining them in the sub-epithelial compartment by increased sphingosine-1-phosphate/S1PR1 signaling. Furthermore, L. reuteri promoted GC-like B cell differentiation, which involved expansion of the GC area and autocrine TGFβ-1 activation. Consequently, PD-1-T follicular helper cell-dependent IgA induction and production was increased by L. reuteri, which shifted the intestinal microbiome and protected against dextran-sulfate-sodium induced colitis and dysbiosis.

CONCLUSIONS

The Peyer's patches sense, enhance and transmit probiotic signals by increasing the numbers and effector functions of distinct B cell subsets, resulting in increased IgA production, altered intestinal microbiota, and protection against inflammation. Video abstract.

Integrating experiments with system-level biogeochemical modeling to understand nitrogen cycling of reservoir sediments at elevated hydrostatic pressure

Impoundment of rivers to construct reservoirs for hydropower and irrigation greatly increase the hydrostatic pressure acting on river sediments with potential repercussions for ecosystem-level microbial activity and metabolism. Understanding the functioning and responses of key biogeochemical cycles such as that of nitrogen cycling to shifting hydrostatic pressure is needed to estimate and predict the systemic nutrient dynamics in deep-water reservoirs. We studied the functioning of bacterial communities involved in nitrogen transformation in bioreactors maintained under contrasting hydrostatic pressures (0.5 MPa-3.0 MPa) and complemented the experimental approach with a functional gene-informed biogeochemical model. The model predictions were broadly consistent with observations from the experiment, suggesting that the rates of N₂O production decreased while the sediment concentration of nitrite increased significantly with increasing pressure, at least when exceeding 1.0 MPa. Changes in nitrite reduction (nirS) and aerobic ammonia oxidation (amoA) genes abundances were in accordance with the observed changes in N₂O production and nitrite levels. Moreover, the model predicted that the higher pressures (P > 1.5 MPa) would intensify the inhibition of N₂ production via denitrification and result in an accumulation of ammonia in the sediment along with a decrease in dissolved oxygen. The results imply that increased hydrostatic pressure caused by dam constructions may have a strong effect on microbial nitrogen conversion, and that this may result in lower nitrogen removal.

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

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.

Freshwater Chlorobia Exhibit Metabolic Specialization among Cosmopolitan and Endemic Populations

Photosynthetic bacteria from the class Chlorobia (formerly phylum Chlorobi) sustain carbon fixation in anoxic water columns. They harvest light at extremely low intensities and use various inorganic electron donors to fix carbon dioxide into biomass. Until now, most information on the functional ecology and local adaptations of Chlorobia members came from isolates and merely 26 sequenced genomes that may not adequately represent natural populations. To address these limitations, we analyzed global metagenomes to profile planktonic Chlorobia cells from the oxyclines of 42 freshwater bodies, spanning subarctic to tropical regions and encompassing all four seasons. We assembled and compiled over 500 genomes, including metagenome-assembled genomes (MAGs), single-amplified genomes (SAGs), and reference genomes from cultures, clustering them into 71 metagenomic operational taxonomic units (mOTUs or "species"). Of the 71 mOTUs, 57 were classified within the genus Chlorobium, and these mOTUs represented up to ∼60% of the microbial communities in the sampled anoxic waters. Several Chlorobium-associated mOTUs were globally distributed, whereas others were endemic to individual lakes. Although most clades encoded the ability to oxidize hydrogen, many lacked genes for the oxidation of specific sulfur and iron substrates. Surprisingly, one globally distributed Scandinavian clade encoded the ability to oxidize hydrogen, sulfur, and iron, suggesting that metabolic versatility facilitated such widespread colonization. Overall, these findings provide new insight into the biogeography of the Chlorobia and the metabolic traits that facilitate niche specialization within lake ecosystems.IMPORTANCE The reconstruction of genomes from metagenomes has helped explore the ecology and evolution of environmental microbiota. We applied this approach to 274 metagenomes collected from diverse freshwater habitats that spanned oxic and anoxic zones, sampling seasons, and latitudes. We demonstrate widespread and abundant distributions of planktonic Chlorobia-associated bacteria in hypolimnetic waters of stratified freshwater ecosystems and show they vary in their capacities to use different electron donors. Having photoautotrophic potential, these Chlorobia members could serve as carbon sources that support metalimnetic and hypolimnetic food webs.

Landscape Setting Drives the Microbial Eukaryotic Community Structure in Four Swedish Mountain Lakes over the Holocene

On the annual and interannual scales, lake microbial communities are known to be heavily influenced by environmental conditions both in the lake and in its terrestrial surroundings. However, the influence of landscape setting and environmental change on shaping these communities over a longer (millennial) timescale is rarely studied. Here, we applied an 18S metabarcoding approach to DNA preserved in Holocene sediment records from two pairs of co-located Swedish mountain lakes. Our data revealed that the microbial eukaryotic communities were strongly influenced by catchment characteristics rather than location. More precisely, the microbial communities from the two bedrock lakes were largely dominated by unclassified Alveolata, while the peatland lakes showed a more diverse microbial community, with Ciliophora, Chlorophyta and Chytrids among the more predominant groups. Furthermore, for the two bedrock-dominated lakes-where the oldest DNA samples are dated to only a few hundred years after the lake formation-certain Alveolata, Chlorophytes, Stramenopiles and Rhizaria taxa were found prevalent throughout all the sediment profiles. Our work highlights the importance of species sorting due to landscape setting and the persistence of microbial eukaryotic diversity over millennial timescales in shaping modern lake microbial communities.

Functionally reversible impacts of disturbances on lake food webs linked to spatial and seasonal dependencies

Increasing human impact on the environment is causing drastic changes in disturbance regimes and how they prevail over time. Of increasing relevance is to further our understanding on biological responses to pulse disturbances (short duration) and how they interact with other ongoing press disturbances (constantly present). Because the temporal and spatial contexts of single experiments often limit our ability to generalize results across space and time, we conducted a modularized mesocosm experiment replicated in space (five lakes along a latitudinal gradient in Scandinavia) and time (two seasons, spring and summer) to generate general predictions on how the functioning and composition of multitrophic plankton communities (zoo-, phyto- and bacterioplankton) respond to pulse disturbances acting either in isolation or combined with press disturbances. As pulse disturbance, we used short-term changes in fish presence, and as press disturbance, we addressed the ongoing reduction in light availability caused by increased cloudiness and lake browning in many boreal and subarctic lakes. First, our results show that the top-down pulse disturbance had the strongest effects on both functioning and composition of the three trophic levels across sites and seasons, with signs for interactive impacts with the bottom-up press disturbance on phytoplankton communities. Second, community composition responses to disturbances were highly divergent between lakes and seasons: temporal accumulated community turnover of the same trophic level either increased (destabilization) or decreased (stabilization) in response to the disturbances compared to control conditions. Third, we found functional recovery from the pulse disturbances to be frequent at the end of most experiments. In a broader context, these results demonstrate that top-down, pulse disturbances, either alone or with additional constant stress upon primary producers caused by bottom-up disturbances, can induce profound but often functionally reversible changes across multiple trophic levels, which are strongly linked to spatial and temporal context dependencies. Furthermore, the identified dichotomy of disturbance effects on the turnover in community composition demonstrates the potential of disturbances to either stabilize or destabilize biodiversity patterns over time across a wide range of environmental conditions.

Methanogenesis Is an Important Process in Controlling MeHg Concentration in Rice Paddy Soils Affected by Mining Activities

Rice paddies are agricultural sites of special concern because the potent toxin methylmercury (MeHg), produced in rice paddy soils, accumulates in rice grains. MeHg cycling is mostly controlled by microbes but their importance in MeHg production and degradation in paddy soils and across a Hg concentration gradient remains unclear. Here we used surface and rhizosphere soil samples in a series of incubation experiments in combination with stable isotope tracers to investigate the relative importance of different microbial groups on MeHg production and degradation across a Hg contamination gradient. We showed that sulfate reduction was the main driver of MeHg formation and concentration at control sites, and that methanogenesis had an important and complex role in MeHg cycling as Hg concentrations increased. The inhibition of methanogenesis at the mining sites led to an increase in MeHg production up to 16.6-fold and a decrease in MeHg degradation by up to 77%, suggesting that methanogenesis is associated with MeHg degradation as Hg concentrations increased. This study broadens our understanding of the roles of microbes in MeHg cycling and highlights methanogenesis as a key control of MeHg concentrations in rice paddies, offering the potential for mitigation of Hg contamination and for the safe production of rice in Hg-contaminated areas.

The environment drives microbial trait variability in aquatic habitats

A prerequisite to improve the predictability of microbial community dynamics is to understand the mechanisms of microbial assembly. To study factors that contribute to microbial community assembly, we examined the temporal dynamics of genes in five aquatic metagenome time-series, originating from marine offshore or coastal sites and one lake. With this trait-based approach we expected to find gene-specific patterns of temporal allele variability that depended on the seasonal metacommunity size of carrier-taxa and the variability of the milieu and the substrates to which the resulting proteins were exposed. In more detail, we hypothesized that a larger seasonal metacommunity size would result in increased temporal variability of functional units (i.e., gene alleles), as shown previously for taxonomic units. We further hypothesized that multicopy genes would feature higher temporal variability than single-copy genes, as gene multiplication can result from high variability in substrate quality and quantity. Finally, we hypothesized that direct exposure of proteins to the extracellular environment would result in increased temporal variability of the respective gene compared to intracellular proteins that are less exposed to environmental fluctuations. The first two hypotheses were confirmed in all data sets, while significant effects of the subcellular location of gene products was only seen in three of the five time-series. The gene with the highest allele variability throughout all data sets was an iron transporter, also representing a target for phage infection. Previous work has emphasized the role of phage-prokaryote interactions as a major driver of microbial diversity. Our finding therefore points to a potentially important role of iron transporter-mediated phage infections for the assembly and maintenance of diversity in aquatic prokaryotes.

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 (the hgcAB gene cluster). We determined the relative abundance of the hgcAB genes and their taxonomic identity in 81 brackish metagenomes that cover spatial, seasonal and redox variability in the Baltic Sea water column. The hgcAB genes were predominantly detected in anoxic water, but some hgcAB genes 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 of hgcAB genes 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.

Diazotroph Genomes and Their Seasonal Dynamics in a Stratified Humic Bog Lake

Aquatic N-fixation is generally associated with the growth and mass development of Cyanobacteria in nitrogen-deprived photic zones. However, sequenced genomes and environmental surveys suggest active aquatic N-fixation also by many non-cyanobacterial groups. Here, we revealed the seasonal variation and genomic diversity of potential N-fixers in a humic bog lake using metagenomic data and nif gene clusters analysis. Groups with diazotrophic operons were functionally divergent and included Cholorobi, Geobacter, Desulfobacterales, Methylococcales, and Acidobacteria. In addition to nifH (a gene that encodes the dinitrogenase reductase component of the molybdenum nitrogenase), we also identified sequences corresponding to vanadium and iron-only nitrogenase genes. Within the Chlorobi population, the nitrogenase (nifH) cluster was included in a well-structured retrotransposon. Furthermore, the presence of light-harvesting photosynthesis genes implies that anoxygenic photosynthesis may fuel nitrogen fixation under the prevailing low-irradiance conditions. The presence of rnf genes (related to the expression of H⁺/Na⁺-translocating ferredoxin: NAD+ oxidoreductase) in Methylococcales and Desulfobacterales suggests that other energy-generating processes may drive the costly N-fixation in the absence of photosynthesis. The highly reducing environment of the anoxic bottom layer of Trout Bog Lake may thus also provide a suitable niche for active N-fixers and primary producers. While future studies on the activity of these potential N-fixers are needed to clarify their role in freshwater nitrogen cycling, the metagenomic data presented here enabled an initial characterization of previously overlooked diazotrophs in freshwater biomes.

Opposing spatial trends in methylmercury and total mercury along a peatland chronosequence trophic gradient

Peatlands are abundant elements of boreal landscapes where inorganic mercury (IHg) can be transformed into bioaccumulating and highly toxic methylmercury (MeHg). We studied fifteen peatlands divided into three age classes (young, intermediate and old) along a geographically constrained chronosequence to determine the role of biogeochemical factors and nutrient availability in controlling the formation of MeHg. In the 10 cm soil layer just below the average annual growing season water table, concentrations of MeHg and %MeHg (of total Hg) were higher in younger, more mesotrophic peatlands than in older, more oligotrophic peatlands. In contrast, total mercury (THg) concentrations were higher in the older peatlands. Partial least squares (PLS) analysis indicates that the net MeHg production was positively correlated to trophic demands of vegetation and an increased availability of potential electron acceptors and donors for Hg methylating microorganisms. An important question for further studies will be to elucidate why there is less THg in the younger peatlands compared to the older peatlands, even though the age of the superficial peat itself is similar for all sites. We hypothesize that ecosystem features which enhance microbial processes involved in Hg methylation also promote Hg reduction that makes previously deposited Hg more available for evasion back to the atmosphere.

Shifts in mercury methylation across a peatland chronosequence: From sulfate reduction to methanogenesis and syntrophy

Peatlands are globally important ecosystems where inorganic mercury is converted to bioaccumulating and highly toxic methylmercury, resulting in high risks of methylmercury exposure in adjacent aquatic ecosystems. Although biological mercury methylation has been known for decades, there is still a lack of knowledge about the organisms involved in mercury methylation and the drivers controlling their methylating capacity. In order to investigate the metabolisms responsible for mercury methylation and methylmercury degradation as well as the controls of both processes, we studied a chronosequence of boreal peatlands covering fundamentally different biogeochemical conditions. Potential mercury methylation rates decreased with peatland age, being up to 53 times higher in the youngest peatland compared to the oldest. Methylation in young mires was driven by sulfate reduction, while methanogenic and syntrophic metabolisms became more important in older systems. Demethylation rates were also highest in young wetlands, with a gradual shift from biotic to abiotic methylmercury degradation along the chronosequence. Our findings reveal how metabolic shifts drive mercury methylation and its ratio to demethylation as peatlands age.

Oral Microbiota Development in Early Childhood

Early life determinants of the oral microbiota have not been thoroughly elucidated. We studied the association of birth and early childhood characteristics with oral microbiota composition using 16 S ribosomal RNA (rRNA) gene sequencing in a population-based Swedish cohort of 59 children sampled at 6, 12 and 24 months of age. Repeated-measurement regression models adjusted for potential confounders confirmed and expanded previous knowledge about the profound shift of oral microbiota composition in early life. These alterations included increased alpha diversity, decreased beta diversity and alteration of bacterial composition with changes in relative abundance of 14 of the 20 most common operational taxonomic units (OTUs). We also found that birth characteristics, breastfeeding and antibiotic use were associated with overall phyla distribution and/or with the relative abundance of specific OTUs. Further, we detected a novel link between morning salivary cortisol level, a physiological marker of neuroendocrine activity and stress, and overall phyla distribution as well as with decreased abundance of the most common OTU mapped to the Streptococcaceae family. In conclusion, a major part of the maturation of the oral microbiome occurs during the first two years of life, and this development may be influenced by early life circumstances.

Uncovering microbial inter-domain interactions in complex communities

Interactions between unicellular eukaryotes and bacteria are difficult to characterize in the environment owing to their large number and inherently microscopic scale. Although particular co-occurrences can be recovered through targeted approaches, e.g. single-cell sequencing or fluorescence in situ hybridization, the vast majority of the interactions remain unseen. Here, we discuss Emulsion, Paired Isolation and Concatenation polymerase chain reaction (epicPCR) as a tool to uncover these interactions in very high throughput. Originally developed for taxonomy-to-function linkage in bacterial communities, epicPCR has the potential to recover the complete interaction network in a given environment at single-cell resolution. This approach relies on the encapsulation of protistan single cells in emulsion droplets that can subsequently be gelified into beads. In this way, encapsulated cells can be exposed to lysis reagents and further phylogenetic paired marker amplification. A bacterium that physically co-occurs with the eukaryote will be jointly trapped, and the amplification will generate a concatenated PCR product containing physically coupled taxonomic markers from both partners, creating a link. Further amplification and sequencing enable the construction of an association pattern with statistically verified physical co-occurrences. Here, we discuss the potential, challenges and limitations of epicPCR. We argue that the microscopic scale at which epicPCR operates, the high throughput it delivers and its exploratory nature make it an unparalleled approach to unravel associations between microbes directly from environmental samples. This article is part of a discussion meeting issue 'Single cell ecology'.

Dietary Fiber in Bilberry Ameliorates Pre-Obesity Events in Rats by Regulating Lipid Depot, Cecal Short-Chain Fatty Acid Formation and Microbiota Composition

Obesity is linked to non-alcoholic fatty liver disease and risk factors associated to metabolic syndrome. Bilberry (Vaccinium myrtillus) that contains easily fermentable fiber may strengthen the intestinal barrier function, attenuate inflammation and modulate gut microbiota composition, thereby prevent obesity development. In the current study, liver lipid metabolism, fat depot, cecal and serum short-chain fatty acids (SCFAs) and gut microbiome were evaluated in rats fed bilberries in a high-fat (HFD + BB) or low-fat (LFD + BB) setting for 8 weeks and compared with diets containing equal amount of fiber resistant to fermentation (cellulose, HFD and LFD). HFD fed rats did not obtain an obese phenotype but underwent pre-obesity events including increased liver index, lipid accumulation and increased serum cholesterol levels. This was linked to shifts of cecal bacterial community and reduction of major SCFAs. Bilberry inclusion improved liver metabolism and serum lipid levels. Bilberry inclusion under either LFD or HFD, maintained microbiota homeostasis, stimulated interscapular-brown adipose tissue depot associated with increased mRNA expression of uncoupling protein-1; enhanced SCFAs in the cecum and circulation; and promoted butyric acid and butyrate-producing bacteria. These findings suggest that bilberry may serve as a preventative dietary measure to optimize microbiome and associated lipid metabolism during or prior to HFD.

Presence of Aedes and Anopheles mosquito larvae is correlated to bacteria found in domestic water-storage containers

Water-storage containers are common in households where access to water is scarce and often act as breeding sites for vector mosquitoes. Bacteria in these containers may be important for attracting or repelling ovipositing mosquitoes. We hypothesized that bacterial community composition in water-storage containers would represent either inhibitory or suitable environmental conditions for mosquito larvae. To investigate this, we characterized the bacterial community composition in water-storage containers and correlated these communities to Aedes and Anopheles larval densities. Water samples were collected over two years from 13 containers in an Indian village and analyzed by high throughput 16S rRNA gene amplicon sequencing. Comparisons of bacterial community composition between water with and without mosquito larvae showed that Xanthomonadaceae, Comamonadaceae and Burkholderiaceae were more common (P < 0.05) in absence of larvae, while Lachnospiraceae, Synechococcaceae, Alcaligenaceae and Cryomorphaceae were more common (P < 0.05) in presence of larvae. Indicator analysis identified operational taxonomic units designated as CL500-29 marine group (Acidimicrobiaceae) and FukuN101 (Microbacteriaceae) for absence and presence of larvae, respectively. These results contribute to the understanding of which bacteria, directly or indirectly, can be linked to absence or presence of mosquitoes around households and set the basis for potential measures to be taken against these vector mosquitoes.

Freshwater carbon and nutrient cycles revealed through reconstructed population genomes

Although microbes mediate much of the biogeochemical cycling in freshwater, the categories of carbon and nutrients currently used in models of freshwater biogeochemical cycling are too broad to be relevant on a microbial scale. One way to improve these models is to incorporate microbial data. Here, we analyze both genes and genomes from three metagenomic time series and propose specific roles for microbial taxa in freshwater biogeochemical cycles. Our metagenomic time series span multiple years and originate from a eutrophic lake (Lake Mendota) and a humic lake (Trout Bog Lake) with contrasting water chemistry. Our analysis highlights the role of polyamines in the nitrogen cycle, the diversity of diazotrophs between lake types, the balance of assimilatory vs. dissimilatory sulfate reduction in freshwater, the various associations between types of phototrophy and carbon fixation, and the density and diversity of glycoside hydrolases in freshwater microbes. We also investigated aspects of central metabolism such as hydrogen metabolism, oxidative phosphorylation, methylotrophy, and sugar degradation. Finally, by analyzing the dynamics over time in nitrogen fixation genes and Cyanobacteria genomes, we show that the potential for nitrogen fixation is linked to specific populations in Lake Mendota. This work represents an important step towards incorporating microbial data into ecosystem models and provides a better understanding of how microbes may participate in freshwater biogeochemical cycling.

Depth and Dissolved Organic Carbon Shape Microbial Communities in Surface Influenced but Not Ancient Saline Terrestrial Aquifers

The continental deep biosphere is suggested to contain a substantial fraction of the earth's total biomass and microorganisms inhabiting this environment likely have a substantial impact on biogeochemical cycles. However, the deep microbial community is still largely unknown and can be influenced by parameters such as temperature, pressure, water residence times, and chemistry of the waters. In this study, 21 boreholes representing a range of deep continental groundwaters from the Äspö Hard Rock Laboratory were subjected to high-throughput 16S rRNA gene sequencing to characterize how the different water types influence the microbial communities. Geochemical parameters showed the stability of the waters and allowed their classification into three groups. These were (i) waters influenced by infiltration from the Baltic Sea with a "modern marine (MM)" signature, (ii) a "thoroughly mixed (TM)" water containing groundwaters of several origins, and (iii) deep "old saline (OS)" waters. Decreasing microbial cell numbers positively correlated with depth. In addition, there was a stronger positive correlation between increased cell numbers and dissolved organic carbon for the MM compared to the OS waters. This supported that the MM waters depend on organic carbon infiltration from the Baltic Sea while the ancient saline waters were fed by "geogases" such as carbon dioxide and hydrogen. The 16S rRNA gene relative abundance of the studied groundwaters revealed different microbial community compositions. Interestingly, the TM water showed the highest dissimilarity compared to the other two water types, potentially due to the several contrasting water types contributing to this groundwater. The main identified microbial phyla in the groundwaters were Gammaproteobacteria, unclassified sequences, Campylobacterota (formerly Epsilonproteobacteria), Patescibacteria, Deltaproteobacteria, and Alphaproteobacteria. Many of these taxa are suggested to mediate ferric iron and nitrate reduction, especially in the MM waters. This indicated that nitrate reduction may be a neglected but important process in the deep continental biosphere. In addition to the high number of unclassified sequences, almost 50% of the identified phyla were archaeal or bacterial candidate phyla. The percentage of unknown and candidate phyla increased with depth, pointing to the importance and necessity of further studies to characterize deep biosphere microbial populations.

Genomes from uncultivated prokaryotes: a comparison of metagenome-assembled and single-amplified genomes

BACKGROUND

Prokaryotes dominate the biosphere and regulate biogeochemical processes essential to all life. Yet, our knowledge about their biology is for the most part limited to the minority that has been successfully cultured. Molecular techniques now allow for obtaining genome sequences of uncultivated prokaryotic taxa, facilitating in-depth analyses that may ultimately improve our understanding of these key organisms.

RESULTS

We compared results from two culture-independent strategies for recovering bacterial genomes: single-amplified genomes and metagenome-assembled genomes. Single-amplified genomes were obtained from samples collected at an offshore station in the Baltic Sea Proper and compared to previously obtained metagenome-assembled genomes from a time series at the same station. Among 16 single-amplified genomes analyzed, seven were found to match metagenome-assembled genomes, affiliated with a diverse set of taxa. Notably, genome pairs between the two approaches were nearly identical (average 99.51% sequence identity; range 98.77-99.84%) across overlapping regions (30-80% of each genome). Within matching pairs, the single-amplified genomes were consistently smaller and less complete, whereas the genetic functional profiles were maintained. For the metagenome-assembled genomes, only on average 3.6% of the bases were estimated to be missing from the genomes due to wrongly binned contigs.

CONCLUSIONS

The strong agreement between the single-amplified and metagenome-assembled genomes emphasizes that both methods generate accurate genome information from uncultivated bacteria. Importantly, this implies that the research questions and the available resources are allowed to determine the selection of genomics approach for microbiome studies.

Temporal variability of bacterioplankton is habitat driven

Temporal dynamics of bacterioplankton are rarely investigated for multiple habitats and years within individual lakes, limiting our understanding of the variability of bacterioplankton community (BC) composition with respect to environmental factors. We assessed the BC composition of a littoral and two pelagic habitats (euphotic zone and hypolimnion) of Lake Tovel monthly from April 2014 to May 2017 by high-throughput sequencing of the V3-V4 hypervariable region of the 16S rRNA gene. The three habitats differed in temperature, light, oxygen and hydrology. In particular, the littoral was the most hydrologically unstable because it receives most of the lake inflow, the hypolimnion was the most stable because of its hydrologically sheltered position, and the pelagic euphotic habitat was intermediate. Consequently, we hypothesized different temporal patterns of BC composition for all three habitats according to their environmental differences. We applied PERMANOVA, nonmetric multidimensional scaling and source-sink analysis to characterize BC composition. Overall, BCs were different among habitats with the littoral showing the highest variability and the hypolimnion the highest stability. The BC of rainy 2014 was distinct from the BCs of other years irrespective of the habitats considered. Seasonal differences in BCs were limited to spring, probably linked to meltwater inflow and mixing. Thus, temporal effects related to year and season were linked to the hydrological gradient of habitats. We suggest that despite potential within-lake dispersal of bacterioplankton by water flow and mixing, local environmental conditions played a major role in Lake Tovel, fostering distinct BCs in the three habitats.

Investigation of viable taxa in the deep terrestrial biosphere suggests high rates of nutrient recycling

The deep biosphere is the largest 'bioreactor' on earth, and microbes inhabiting this biome profoundly influence global nutrient and energy cycles. An important question for deep biosphere microbiology is whether or not specific populations are viable. To address this, we used quantitative PCR and high throughput 16S rRNA gene sequencing of total and viable cells (i.e. with an intact cellular membrane) from three groundwaters with different ages and chemical constituents. There were no statistically significant differences in 16S rRNA gene abundances and microbial diversity between total and viable communities. This suggests that populations were adapted to prevailing oligotrophic conditions and that non-viable cells are rapidly degraded and recycled into new biomass. With higher concentrations of organic carbon, the modern marine and undefined mixed waters hosted a community with a larger range of predicted growth strategies than the ultra-oligotrophic old saline water. These strategies included fermentative and potentially symbiotic lifestyles by candidate phyla that typically have streamlined genomes. In contrast, the old saline waters had more 16S rRNA gene sequences in previously cultured lineages able to oxidize hydrogen and fix carbon dioxide. This matches the paradigm of a hydrogen and carbon dioxide-fed chemolithoautotrophic deep biosphere.

An in vitro evaluation of browser and grazer fermentation efficiency and microbiota using European moose spring and summer foods

Evolutionary morphological and physiological differences between browsers and grazers contribute to species‐specific digestion efficiency of food resources. Rumen microbial community structure of browsers is supposedly adapted to characteristic nutrient composition of the diet source. If this assumption is correct, domesticated ruminants, or grazers, are poor model animals for assessing the nutritional value of food consumed by browsing game species. In this study, typical spring and summer foods of the European moose (Alces alces) were combined with rumen fluid collected from both dairy cows (Bos taurus) and from moose, with the aim of comparing fermentation efficiency and microbial community composition. The nutritional value of the food resources was characterized by chemical analysis and advanced in vitro measurements. The study also addressed whether or not feed evaluation based on in vitro techniques with cattle rumen fluid as inoculum could be a practical alternative when evaluating the nutritional value of plants consumed by wild browsers. Our results suggest that the fermentation characteristics of moose spring and summer food are partly host‐specific and related to the contribution of the bacterial phyla Firmicutes and Bacteriodetes to the rumen microbial community. Host‐specific adaptations of the ruminal microbial community structure could be explained from the evolutionary adaptations related to feeding habitats and morphophysiological differences between browsers and grazers. However, the observed overall differences in microbial community structure could not be related to ruminal digestion parameters measured in vitro. The in vitro evaluation of digestion efficiency reveals that equal amounts of methane were produced across all feed samples regardless of whether the ruminal fluid was from moose or dairy cow. The results of this study suggested that the nutritional value of browsers' spring and summer food can be predicted using rumen fluid from domesticated grazers as inoculum in in vitro assessments of extent of digestion when excluding samples of the white water lily root, but not of fermentation characteristics as indicated by the proportions of individual fermentation fatty acids to the total of volatile fatty acids.

Formation of mercury methylation hotspots as a consequence of forestry operations

Earlier studies have shown that boreal forest logging can increase the concentration and export of methylmercury (MeHg) in stream runoff. Here we test whether forestry operations create soil environments of high MeHg net formation associated with distinct microbial communities. Furthermore, we test the hypothesis that Hg methylation hotspots are more prone to form after stump harvest than stem-only harvest, because of more severe soil compaction and soil disturbance. Concentrations of MeHg, percent MeHg of total Hg (THg), and bacterial community composition were determined at 200 soil sampling positions distributed across eight catchments. Each catchment was either stem-only harvested (n=3), stem- and stump-harvested (n=2) or left undisturbed (n=3). In support of our hypothesis, higher MeHg to THg ratios was observed in one of the stump-harvested catchments. While the effects of natural variation could not be ruled out, we noted that most of the highest % MeHg was observed in water-filled cavities created by stump removal or driving damage. This catchment also featured the highest bacterial diversity and highest relative abundance of bacterial families known to include Hg methylators. We propose that water-logged and disturbed soil environments associated with stump harvest can favor methylating microorganisms, which also enhance MeHg formation.

Geobacteraceae are important members of mercury-methylating microbial communities of sediments impacted by waste water releases

Microbial mercury (Hg) methylation in sediments can result in bioaccumulation of the neurotoxin methylmercury (MMHg) in aquatic food webs. Recently, the discovery of the gene hgcA, required for Hg methylation, revealed that the diversity of Hg methylators is much broader than previously thought. However, little is known about the identity of Hg-methylating microbial organisms and the environmental factors controlling their activity and distribution in lakes. Here, we combined high-throughput sequencing of 16S rRNA and hgcA genes with the chemical characterization of sediments impacted by a waste water treatment plant that releases significant amounts of organic matter and iron. Our results highlight that the ferruginous geochemical conditions prevailing at 1–2 cm depth are conducive to MMHg formation and that the Hg-methylating guild is composed of iron and sulfur-transforming bacteria, syntrophs, and methanogens. Deltaproteobacteria, notably Geobacteraceae, dominated the hgcA carrying communities, while sulfate reducers constituted only a minor component, despite being considered the main Hg methylators in many anoxic aquatic environments. Because iron is widely applied in waste water treatment, the importance of Geobacteraceae for Hg methylation and the complexity of Hg-methylating communities reported here are likely to occur worldwide in sediments impacted by waste water treatment plant discharges and in iron-rich sediments in general.