Climate change threatens to release abundant carbon that is sequestered at high latitudes, but the constraints on microbial metabolisms that mediate the release of methane and carbon dioxide are poorly understood. The role of viruses, which are known to affect microbial dynamics, metabolism and biogeochemistry in the oceans, remains largely unexplored in soil. Here, we aimed to investigate how viruses influence microbial ecology and carbon metabolism in peatland soils along a permafrost thaw gradient in Sweden. We recovered 1,907 viral populations (genomes and large genome fragments) from 197 bulk soil and size-fractionated metagenomes, 58% of which were detected in metatranscriptomes and presumed to be active. In silico predictions linked 35% of the viruses to microbial host populations, highlighting likely viral predators of key carbon-cycling microorganisms, including methanogens and methanotrophs. Lineage-specific virus/host ratios varied, suggesting that viral infection dynamics may differentially impact microbial responses to a changing climate. Virus-encoded glycoside hydrolases, including an endomannanase with confirmed functional activity, indicated that viruses influence complex carbon degradation and that viral abundances were significant predictors of methane dynamics. These findings suggest that viruses may impact ecosystem function in climate-critical, terrestrial habitats and identify multiple potential viral contributions to soil carbon cycling.

As global temperatures rise, large amounts of carbon sequestered in permafrost are becoming available for microbial degradation. Accurate prediction of carbon gas emissions from thawing permafrost is limited by our understanding of these microbial communities. Here we use metagenomic sequencing of 214 samples from a permafrost thaw gradient to recover 1,529 metagenome-assembled genomes, including many from phyla with poor genomic representation. These genomes reflect the diversity of this complex ecosystem, with genus-level representatives for more than sixty per cent of the community. Meta-omic analysis revealed key populations involved in the degradation of organic matter, including bacteria whose genomes encode a previously undescribed fungal pathway for xylose degradation. Microbial and geochemical data highlight lineages that correlate with the production of greenhouse gases and indicate novel syntrophic relationships. Our findings link changing biogeochemistry to specific microbial lineages involved in carbon processing, and provide key information for predicting the effects of climate change on permafrost systems.

No abstract

The fate of carbon sequestered in permafrost is a key concern for future global warming as this large carbon stock is rapidly becoming a net methane source due to widespread thaw. Methane release from permafrost is moderated by methanotrophs, which oxidise 20–60% of this methane before emission to the atmosphere. Despite the importance of methanotrophs to carbon cycling, these microorganisms are under-characterised and have not been studied across a natural permafrost thaw gradient. Here, we examine methanotroph communities from the active layer of a permafrost thaw gradient in Stordalen Mire (Abisko, Sweden) spanning three years, analysing 188 metagenomes and 24 metatranscriptomes paired with in situ biogeochemical data. Methanotroph community composition and activity varied significantly as thaw progressed from intact permafrost palsa, to partially thawed bog and fully thawed fen. Thirteen methanotroph population genomes were recovered, including two novel genomes belonging to the uncultivated upland soil cluster alpha (USCα) group and a novel potentially methanotrophic Hyphomicrobiaceae. Combined analysis of porewater δ13C-CH4 isotopes and methanotroph abundances showed methane oxidation was greatest below the oxic–anoxic interface in the bog. These results detail the direct effect of thaw on autochthonous methanotroph communities, and their consequent changes in population structure, activity and methane moderation potential.

Large-scale metagenomic datasets enable the recovery of hundreds of population genomes from environmental samples. However, these genomes do not typically represent the full diversity of complex microbial communities. Gene-centric approaches can be used to gain a comprehensive view of diversity by examining each read independently, but traditional pairwise comparison approaches typically over-classify taxonomy and scale poorly with increasing metagenome and database sizes. Here we introduce GraftM, a tool that uses gene specific packages to rapidly identify gene families in metagenomic data using hidden Markov models (HMMs) or DIAMOND databases, and classifies these sequences using placement into pre-constructed gene trees. The speed and accuracy of GraftM was benchmarked with in silico and in vitro mock communities using taxonomic markers, and was found to have higher accuracy at the family level with a processing time 2.0–3.7× faster than currently available software. Exploration of a wetland metagenome using 16S rRNA- and methyl-coenzyme M reductase (McrA)-specific gpkgs revealed taxonomic and functional shifts across a depth gradient. Analysis of the NCBI nr database using the McrA gpkg allowed the detection of novel sequences belonging to phylum-level lineages. A growing collection of gpkgs is available online (https://github.com/geronimp/graftM_gpkgs), where curated packages can be uploaded and exchanged.

Sponges underpin the productivity of coral reefs, yet few of their microbial symbionts have been functionally characterised. Here we present an analysis of ~1200 metagenome-assembled genomes (MAGs) spanning seven sponge species and 25 microbial phyla. Compared to MAGs derived from reef seawater, sponge-associated MAGs were enriched in glycosyl hydrolases targeting components of sponge tissue, coral mucus and macroalgae, revealing a critical role for sponge symbionts in cycling reef organic matter. Further, visualisation of the distribution of these genes amongst symbiont taxa uncovered functional guilds for reef organic matter degradation. Genes for the utilisation of sialic acids and glycosaminoglycans present in sponge tissue were found in specific microbial lineages that also encoded genes for attachment to sponge-derived fibronectins and cadherins, suggesting these lineages can utilise specific structural elements of sponge tissue. Further, genes encoding CRISPR and restriction-modification systems used in defence against mobile genetic elements were enriched in sponge symbionts, along with eukaryote-like gene motifs thought to be involved in maintaining host association. Finally, we provide evidence that many of these sponge-enriched genes are laterally transferred between microbial taxa, suggesting they confer a selective advantage within the sponge niche and therefore play a critical role in host ecology and evolution.

The methyl-coenzyme M reductase (MCR) complex is a key enzyme in archaeal methane generation and has recently been proposed to also be involved in the oxidation of short-chain hydrocarbons including methane, butane, and potentially propane. The number of archaeal clades encoding the MCR continues to grow, suggesting that this complex was inherited from an ancient ancestor, or has undergone extensive horizontal gene transfer. Expanding the representation of MCR-encoding lineages through metagenomic approaches will help resolve the evolutionary history of this complex. Here, a near-complete Archaeoglobi metagenome-assembled genome (MAG; Ca . Polytropus marinifundus gen. nov. sp. nov.) was recovered from the deep subseafloor along the Juan de Fuca Ridge flank that encodes two divergent McrABG operons similar to those found in Ca . Bathyarchaeota and Ca . Syntrophoarchaeum MAGs. Ca . P. marinifundus is basal to members of the class Archaeoglobi, and encodes the genes for β-oxidation, potentially allowing an alkanotrophic metabolism similar to that proposed for Ca . Syntrophoarchaeum. Ca . P. marinifundus also encodes a respiratory electron transport chain that can potentially utilize nitrate, iron, and sulfur compounds as electron acceptors. Phylogenetic analysis suggests that the Ca . P. marinifundus MCR operons were horizontally transferred, changing our understanding of the evolution and distribution of this complex in the Archaea.

Summary: Finding and translating stretches of DNA lacking stop codons is a task common in the analysis of sequence data. However, the computational tools for finding open reading frames are sufficiently slow that they are becoming a bottleneck as the volume of sequence data grows. This computational bottleneck is especially problematic in metagenomics when searching unassembled reads, or screening assembled contigs for genes of interest. Here, we present OrfM, a tool to rapidly identify open reading frames (ORFs) in sequence data by applying the Aho–Corasick algorithm to find regions uninterrupted by stop codons. Benchmarking revealed that OrfM finds identical ORFs to similar tools (‘GetOrf’ and ‘Translate’) but is four-five times faster. While OrfM is sequencing platform-agnostic, it is best suited to large, high quality datasets such as those produced by Illumina sequencers.Availability and Implementation: Source code and binaries are freely available for download at http://github.com/wwood/OrfM or through GNU Guix under the LGPL 3+ license. OrfM is implemented in C and supported on GNU/Linux and OSX.Contacts: b.woodcroft@uq.edu.auSupplementary information: Supplementary data are available at Bioinformatics online.