Metagenomes and metagenome assembled genomes from anaerobic digesters at three Canadian pulp and paper mills.

We present a dataset of six metagenomes and 323 metagenome assembled genomes (MAGs) describing the microbial community of anaerobic digesters at three Canadian pulp and paper mills. Our objective was to assess the coding potential of the microbial community and obtain draft genomes of key organisms in the digesters.

Metagenomic and genomic sequences from a nitrate-reducing benzene-degrading enrichment culture.

Metagenome-assembled genomes (MAGs) were recovered from metagenomic assemblies from a nitrate-reducing benzene-degrading enrichment culture. Ten MAGs of high quality or functional interest to benzene degradation are reported, seven of which are single contig genomes.

Draft genome of an anaerobic nitrate-reducing, benzene-degrading member of the order Thermincolales.

We present a metagenome assembled genome (MAG) of an anaerobic bacterium from a nitrate-reducing, benzene-degrading enrichment culture (NRBC). The draft Thermincolales genome consists of 20 contigs with a total length of 4.09 Mbp and includes putative carboxylase genes likely involved in benzene activation.

Metagenomes and metagenome-assembled genomes from a sequentially fed anaerobic digester treating solid organic municipal waste.

We present a data set of four metagenomes and 281 metagenome-assembled genomes describing the microbial community of a laboratory-scale high solids anaerobic digester. Our objective was to obtain information on the coding potential of the microbial community and draft genomes of the most abundant organisms in the digester.

Early Divergence and Gene Exchange Highways in the Evolutionary History of Mesoaciditogales.

The placement of a nonhyperthermophilic order Mesoaciditogales as the earliest branching clade within the Thermotogota phylum challenges the prevailing hypothesis that the last common ancestor of Thermotogota was a hyperthermophile. Yet, given the long branch leading to the only two Mesoaciditogales described to date, the phylogenetic position of the order may be due to the long branch attraction artifact. By testing various models and applying data recoding in phylogenetic reconstructions, we observed that early branching of Mesoaciditogales within Thermotogota is strongly supported by the conserved marker genes assumed to be vertically inherited. However, based on the taxonomic content of 1,181 gene families and a phylogenetic analysis of 721 gene family trees, we also found that a substantial number of Mesoaciditogales genes are more closely related to species from the order Petrotogales. These genes contribute to coenzyme transport and metabolism, fatty acid biosynthesis, genes known to respond to heat and cold stressors, and include many genes of unknown functions. The Petrotogales comprise moderately thermophilic and mesophilic species with similar temperature tolerances to that of Mesoaciditogales. Our findings hint at extensive horizontal gene transfer (HGT) between, or parallel independent gene gains by, the two ecologically similar lineages and suggest that the exchanged genes may be important for adaptation to comparable temperature niches.

Metagenomic and Genomic Sequences from a Methanogenic Benzene-Degrading Consortium.

Draft and complete metagenome assembled genomes (MAGs) were created from multiple metagenomic assemblies of DGG-B, a strictly anaerobic, stable mixed microbial consortium that degrades benzene completely to methane and CO2. Our objective was to obtain closed genome sequences of benzene-fermenting bacteria to enable the elucidation of their elusive anaerobic benzene degradation pathway.

"Candidatus Nealsonbacteria" Are Likely Biomass Recycling Ectosymbionts of Methanogenic Archaea in a Stable Benzene-Degrading Enrichment Culture.

The Candidate Phyla Radiation (CPR), also referred to as superphylum Patescibacteria, is a very large group of bacteria with no pure culture representatives discovered by 16S rRNA sequencing or genome-resolved metagenomic analyses of environmental samples. Within the CPR, candidate phylum Parcubacteria, previously referred to as OD1, is prevalent in anoxic sediments and groundwater. Previously, we had identified a specific member of the Parcubacteria (referred to as DGGOD1a) as an important member of a methanogenic benzene-degrading consortium. Phylogenetic analyses herein place DGGOD1a within the clade "Candidatus Nealsonbacteria." Because of its persistence over many years, we hypothesized that "Ca. Nealsonbacteria" DGGOD1a must play an important role in sustaining anaerobic benzene metabolism in the consortium. To try to identify its growth substrate, we amended the culture with a variety of defined compounds (pyruvate, acetate, hydrogen, DNA, and phospholipid), as well as crude culture lysate and three subfractions thereof. We observed the greatest (10-fold) increase in the absolute abundance of "Ca. Nealsonbacteria" DGGOD1a only when the consortium was amended with crude cell lysate. These results implicate "Ca. Nealsonbacteria" in biomass recycling. Fluorescence in situ hybridization and cryogenic transmission electron microscope images revealed that "Ca. Nealsonbacteria" DGGOD1a cells were attached to larger archaeal Methanothrix cells. This apparent epibiont lifestyle was supported by metabolic predictions from a manually curated complete genome. This is one of the first examples of bacterial-archaeal episymbiosis and may be a feature of other "Ca. Nealsonbacteria" found in anoxic environments. IMPORTANCE An anaerobic microbial enrichment culture was used to study members of candidate phyla that are difficult to grow in the lab. We were able to visualize tiny "Candidatus Nealsonbacteria" cells attached to a large Methanothrix cell, revealing a novel episymbiosis.

Complete Genome Sequence of Oceanotoga sp. Strain T3B (DSM 15011).

Oceanotoga sp. strain T3B was isolated from an estuarine sinkhole in the Bahamas. Here, we report its complete genome, which is currently the only sequenced genome from the genus Oceanotoga. The genome sequence provides new data for the genus Oceanotoga.

Microbial community dynamics of a sequentially fed anaerobic digester treating solid organic waste.

A 50-kg scale, high solids anaerobic digester (AD) comprising six sequentially fed leach beds with a leachate recirculation system was operated at 37°C for 88 weeks. The solid feedstock contained a constant fibre fraction (a mix of cardboard, boxboard, newsprint, and fine paper) and varying proportions of food waste. Previously, we reported on the stable operation of this digestion system, where significantly enhanced methane production from the fibre fraction was observed as the proportion of food waste increased. The objective of this study was to identify relationships between process parameters and the microbial community. Increasing food waste led to a large increase in the absolute microbial abundance in the circulating leachate. While 16S rRNA amplicons for Clostridium butyricum were most abundant and correlated with the amount of FW in the system and with the overall methane yield, it was more cryptic Candidatus Roizmanbacteria and Spirochaetaceae that correlated specifically with enhanced methane from the fiber fraction. A faulty batch of bulking agent led to hydraulic channeling, which was reflected in the leachate microbial profiles matching that of the incoming food waste. The system performance and microbial community re-established rapidly after reverting to better bulking agent, illustrating the robustness of the system.

Plasmid pMO1 from Marinitoga okinawensis, first non-cryptic plasmid reported within Thermotogota.

Mobile genetic elements (MGEs), such as viruses and plasmids, drive the evolution and adaptation of their cellular hosts from all three domains of life. This includes microorganisms thriving in the most extreme environments, like deep-sea hydrothermal vents. However, our knowledge about MGEs still remains relatively sparse in these abyssal ecosystems. Here we report the isolation, sequencing, assembly, and functional annotation of pMO1, a 28.2 kbp plasmid associated with the reference strain Marinitoga okinawensis. Carrying restriction/modification and chemotaxis protein-encoding genes, pMO1 likely affects its host's phenotype and represents the first non-cryptic plasmid described among the phylum Thermotogota.

Complex impacts of hydraulic fracturing return fluids on soil microbial community respiration, structure and functional potentials.

The consequences of soils exposed to hydraulic fracturing (HF) return fluid, often collectively termed flowback and produced water (FPW), are poorly understood, even though soils are a common receptor of FPW spills. Here, we investigate the impacts on soil microbiota exposed to FPW collected from the Montney Formation of western Canada. We measured soil respiration, microbial community structure and functional potentials under FPW exposure across a range of concentrations, exposure time and soil types (luvisol and chernozem). We find that soil type governs microbial community response upon FPW exposure. Within each soil, FPW exposure led to reduced biotic soil respiration, and shifted microbial community structure and functional potentials. We detect substantially higher species richness and more unique functional genes in FPW-exposed soils than in FPW-unexposed soils, with metagenome-assembled genomes (e.g. Marinobacter persicus) from luvisol soil exposed to concentrated FPW being most similar to genomes from HF/FPW sites. Our data demonstrate the complex impacts of microbial communities following FPW exposure and highlight the site-specific effects in evaluation of spills and agricultural reuse of FPW on the normal soil functions.

Newly identified proviruses in Thermotogota suggest that viruses are the vehicles on the highways of interphylum gene sharing.

Phylogenomic analyses of bacteria from the phylum Thermotogota have shown extensive lateral gene transfer with distantly related organisms, particularly with Firmicutes. One likely mechanism of such DNA transfer is viruses. However, to date, only three temperate viruses have been characterized in this phylum, all infecting bacteria from the Marinitoga genus. Here we report 17 proviruses integrated into genomes of bacteria belonging to eight Thermotogota genera and induce viral particle production from one of the proviruses. All except an incomplete provirus from Mesotoga fall into two groups based on sequence similarity, gene synteny and taxonomic classification. Proviruses of Group 1 are found in the genera Geotoga, Kosmotoga, Marinitoga, Thermosipho and Mesoaciditoga and are similar to the previously characterized Marinitoga viruses, while proviruses from Group 2 are distantly related to the Group 1 proviruses, have different genome organization and are found in Petrotoga and Defluviitoga. Genes carried by both groups are closely related to Firmicutes and Firmicutes (pro)viruses in phylogenetic analyses. Moreover, one of the groups show evidence of recent gene exchange and may be capable of infecting cells from both phyla. We hypothesize that viruses are responsible for a large portion of the observed gene flow between Firmicutes and Thermotogota.

The first head-tailed virus, MFTV1, infecting hyperthermophilic methanogenic deep-sea archaea.

Deep-sea hydrothermal vents are inhabited by complex communities of microbes and their viruses. Despite the importance of viruses in controlling the diversity, adaptation and evolution of their microbial hosts, to date, only eight bacterial and two archaeal viruses isolated from abyssal ecosystems have been described. Thus, our efforts focused on gaining new insights into viruses associated with deep-sea autotrophic archaea. Here, we provide the first evidence of an infection of hyperthermophilic methanogenic archaea by a head-tailed virus, Methanocaldococcus fervens tailed virus 1 (MFTV1). MFTV1 has an isometric head of 50 nm in diameter and a 150 nm-long non-contractile tail. Virions are released continuously without causing a sudden drop in host growth. MFTV1 infects Methanocaldococcus species and is the first hyperthermophilic head-tailed virus described thus far. The viral genome is a double-stranded linear DNA of 31 kb. Interestingly, our results suggest potential strategies adopted by the plasmid pMEFER01, carried by M. fervens, to spread horizontally in hyperthermophilic methanogens. The data presented here open a new window of understanding on how the abyssal mobilome interacts with hyperthermophilic marine archaea.

Response of aquatic microbial communities and bioindicator modelling of hydraulic fracturing flowback and produced water.

The response of microbial communities to releases of hydraulic fracturing flowback and produced water (PW) may influence ecosystem functions. However, knowledge of the effects of PW spills on freshwater microbiota is limited. Here, we conducted two separate experiments: 16S rRNA gene sequencing combined with random forests modelling was used to assess freshwater community changes in simulated PW spills by volume from 0.05% to 50%. In a separate experiment, live/dead cell viability in a freshwater community was tested during exposure to 10% PW by volume. Three distinct patterns of microbial community shifts were identified: (i) indigenous freshwater genera remained dominant in <2.5% PW, (ii) from 2.5% to 5% PW, potential PW organic degraders such as Pseudomonas, Rheinheimera and Brevundimonas became dominant, and (iii) no significant change in the relative abundance of taxa was observed in >5% PW. Microbial taxa including less abundant genera such as Cellvibrio were potential bioindicators for the degree of contamination with PW. Additionally, live cells were quickly damaged by adding 10% PW, but cell counts recovered in the following days. Our study shows that the responses of freshwater microbiota vary by spill size, and these responses show promise as effective fingerprints for PW spills in aquatic environments.

Complete Genome Sequence of Acidithiobacillus ferridurans JAGS, Isolated from Acidic Mine Drainage.

We report a complete genome sequence of Acidithiobacillus ferridurans JAGS, determined using PacBio single-molecule real-time (SMRT) sequencing. The circular genome of JAGS (2,933,811 bp; GC content, 58.57%) contains 3,001 protein-coding sequences, 46 tRNAs, and 6 rRNAs. Predicted genes indicate the potential to fix CO2 and N2 and to utilize Fe2+, S0, and H2 as energy sources.

Insights into origins and function of the unexplored majority of the reductive dehalogenase gene family as a result of genome assembly and ortholog group classification.

Organohalide respiring bacteria (OHRB) express reductive dehalogenases for energy conservation and growth. Some of these enzymes catalyze the reductive dehalogenation of chlorinated and brominated pollutants in anaerobic subsurface environments, providing a valuable ecosystem service. Dehalococcoides mccartyi strains have been most extensively studied owing to their ability to dechlorinate all chlorinated ethenes - most notably carcinogenic vinyl chloride - to ethene. The genomes of OHRB, particularly obligate OHRB, often harbour multiple putative reductive dehalogenase genes (rdhA), most of which have yet to be characterized. We recently sequenced and closed the genomes of eight new strains, increasing the number of available D. mccartyi genomes in NCBI from 16 to 24. From all available OHRB genomes, we classified predicted translations of reductive dehalogenase genes using a previously established 90% amino acid pairwise identity cut-off to identify Ortholog Groups (OGs). Interestingly, the majority of D. mccartyi dehalogenase gene sequences, once classified into OGs, exhibited a remarkable degree of synteny (gene order) in all genomes sequenced to date. This organization was not apparent without the classification. A high degree of synteny indicates that differences arose from rdhA gene loss rather than recombination. Phylogenetic analysis suggests that most rdhA genes have a long evolutionary history in the Dehalococcoidia with origin prior to speciation of Dehalococcoides and Dehalogenimonas. We also looked for evidence of synteny in the genomes of other species of OHRB. Unfortunately, there are too few closed Dehalogenimonas genomes to compare at this time. There is some partial evidence for synteny in the Dehalobacter restrictus genomes, but here too more closed genomes are needed for confirmation. Interestingly, we found that the rdhA genes that encode enzymes that catalyze dehalogenation of industrial pollutants are the only rdhA genes with strong evidence of recent lateral transfer - at least in the genomes examined herein. Given the utility of the RdhA sequence classification to comparative analyses, we are building a public web server () for the community to use, which allows users to add and classify new sequences, and download the entire curated database of reductive dehalogenases.

Antibiotic resistome and microbial community structure during anaerobic co-digestion of food waste, paper and cardboard.

Solid organic waste is a significant source of antibiotic resistance genes (ARGs) and effective treatment strategies are urgently required to limit the spread of antimicrobial resistance. Here, we studied ARG diversity and abundance as well as the relationship between antibiotic resistome and microbial community structure within a lab-scale solid-state anaerobic digester treating a mixture of food waste, paper and cardboard. A total of 10 samples from digester feed and digestion products were collected for microbial community analysis including small subunit rRNA gene sequencing, total community metagenome sequencing and high-throughput quantitative PCR. We observed a significant shift in microbial community composition and a reduction in ARG diversity and abundance after 6 weeks of digestion. ARGs were identified in all samples with multidrug resistance being the most abundant ARG type. Thirty-two per cent of ARGs detected in digester feed were located on plasmids indicating potential for horizontal gene transfer. Using metagenomic assembly and binning, we detected potential bacterial hosts of ARGs in digester feed, which included Erwinia, Bifidobacteriaceae, Lactococcus lactis and Lactobacillus. Our results indicate that the process of sequential solid-state anaerobic digestion of food waste, paper and cardboard tested herein provides a significant reduction in the relative abundance of ARGs per 16S rRNA gene.

Genomic analysis of the mesophilic Thermotogae genus Mesotoga reveals phylogeographic structure and genomic determinants of its distinct metabolism.

The genus Mesotoga, the only described mesophilic Thermotogae lineage, is common in mesothermic anaerobic hydrocarbon-rich environments. Besides mesophily, Mesotoga displays lineage-specific phenotypes, such as no or little H2 production and dependence on sulfur-compound reduction, which may influence its ecological role. We used comparative genomics of 18 Mesotoga strains (pairwise 16S rRNA identity >99%) and a transcriptome of M. prima to investigate how life at moderate temperatures affects phylogeography and to interrogate the genomic features of its lineage-specific metabolism. We propose that Mesotoga accomplish H2 oxidation and thiosulfate reduction using a sulfide dehydrogenase and a hydrogenase-complex and that a pyruvate:ferredoxin oxidoreductase acquired from Clostridia is responsible for oxidizing acetate. Phylogenetic analysis revealed three distinct Mesotoga lineages (89.6%-99.9% average nucleotide identity [ANI] within lineages, 79.3%-87.6% ANI between lineages) having different geographic distribution patterns and high levels of intra-lineage recombination but little geneflow between lineages. Including data from metagenomes, phylogeographic patterns suggest that geographical separation historically has been more important for Mesotoga than hyperthermophilic Thermotoga and we hypothesize that distribution of Mesotoga is constrained by their anaerobic lifestyle. Our data also suggest that recent anthropogenic activities and environments (e.g., wastewater treatment, oil exploration) have expanded Mesotoga habitats and dispersal capabilities.

Thermosipho spp. Immune System Differences Affect Variation in Genome Size and Geographical Distributions.

Thermosipho species inhabit thermal environments such as marine hydrothermal vents, petroleum reservoirs, and terrestrial hot springs. A 16S rRNA phylogeny of available Thermosipho spp. sequences suggested habitat specialists adapted to living in hydrothermal vents only, and habitat generalists inhabiting oil reservoirs, hydrothermal vents, and hotsprings. Comparative genomics of 15 Thermosipho genomes separated them into three distinct species with different habitat distributions: The widely distributed T. africanus and the more specialized, T. melanesiensis and T. affectus. Moreover, the species can be differentiated on the basis of genome size (GS), genome content, and immune system composition. For instance, the T. africanus genomes are largest and contained the most carbohydrate metabolism genes, which could explain why these isolates were obtained from ecologically more divergent habitats. Nonetheless, all the Thermosipho genomes, like other Thermotogae genomes, show evidence of genome streamlining. GS differences between the species could further be correlated to differences in defense capacities against foreign DNA, which influence recombination via HGT. The smallest genomes are found in T. affectus that contain both CRISPR-cas Type I and III systems, but no RM system genes. We suggest that this has caused these genomes to be almost devoid of mobile elements, contrasting the two other species genomes that contain a higher abundance of mobile elements combined with different immune system configurations. Taken together, the comparative genomic analyses of Thermosipho spp. revealed genetic variation allowing habitat differentiation within the genus as well as differentiation with respect to invading mobile DNA.

Electron donor-driven bacterial and archaeal community patterns along forest ring edges in Ontario, Canada.

Forest rings are 50-1600 m diameter circular structures found in boreal forests around the globe. They are believed to be chemically reducing chimney features, having an accumulation of reduced species in the middle of the ring and oxidation processes occurring at the ring's edges. It has been suggested that microorganisms could be responsible for charge transfer from the inside to the outside of the ring. To explore this, we focused on the changes in bacterial and archaeal communities in the ring edges of two forest rings, the 'Bean' and the 'Thorn North' ring, in proximity to each other in Ontario, Canada. The drier samples from the methane-sourced Bean ring were characterized by the abundance of bacteria from the classes Deltaproteobacteria and Gemmatimonadetes. Geobacter spp. and methanotrophs, such as Candidatus Methylomirabilis and Methylobacter, were highly abundant in these samples. The Thorn North ring, centred on an H2 S accumulation in groundwater, had wetter samples and its communities were dominated by the classes Alphaproteobacteria and Anaerolineae. This ring's microbial communities showed an overall higher microbial diversity supported by higher available free energy. For both rings, the species diversity was highest near the borders of the 20-30 m broad ring edges.