Microbial Functional Diversity Correlates with Species Diversity along a Temperature Gradient.

Microbial community diversity is often correlated with physical environmental stresses like acidity, salinity, and temperature. For example, species diversity usually declines with increasing temperature above 20°C. However, few studies have examined whether the genetic functional diversity of community metagenomes varies in a similar way as species diversity along stress gradients. Here, we investigated bacterial communities in thermal spring sediments ranging from 21 to 88°C, representing communities of 330 to 3,800 bacterial and archaeal species based on 16S rRNA gene amplicon analysis. Metagenomes were sequenced, and Pfam abundances were used as a proxy for metagenomic functional diversity. Significant decreases in both species diversity and Pfam diversity were observed with increasing temperatures. The relationship between Pfam diversity and species diversity followed a power function with the steepest slopes in the high-temperature, low-diversity region of the gradient. Species additions to simple thermophilic communities added many new Pfams, while species additions to complex mesophilic communities added relatively fewer new Pfams, indicating that species diversity does not approach saturation as rapidly as Pfam diversity does. Many Pfams appeared to have distinct temperature ceilings of 60 to 80°C. This study suggests that temperature stress limits both taxonomic and functional diversity of microbial communities, but in a quantitatively different manner. Lower functional diversity at higher temperatures is probably due to two factors, including (i) the absence of many enzymes not adapted to thermophilic conditions, and (ii) the fact that high-temperature communities are comprised of fewer species with smaller average genomes and, therefore, contain fewer rare functions. IMPORTANCE Only recently have microbial ecologists begun to assess quantitatively how microbial species diversity correlates with environmental factors like pH, temperature, and salinity. However, still, very few studies have examined how the number of distinct biochemical functions of microbial communities, termed functional diversity, varies with the same environmental factors. Our study examined 18 microbial communities sampled across a wide temperature gradient and found that increasing temperature reduced both species and functional diversity, but in different ways. Initially, functional diversity increased sharply with increasing species diversity but eventually plateaued, following a power function. This pattern has been previously predicted in theoretical models, but our study validates this predicted power function with field metagenomic data. This study also presents a unique overview of the distribution of metabolic functions along a temperature gradient, demonstrating that many functions have temperature "ceilings" above which they are no longer found.

Seasonal Dynamics of Methanotrophic Bacteria in a Boreal Oil Sands End Pit Lake.

Base Mine Lake (BML) is the first full-scale demonstration end pit lake for the oil sands mining industry in Canada. We examined aerobic methanotrophic bacteria over all seasons for 5 years in this dimictic lake. Methanotrophs comprised up to 58% of all bacterial reads in 16S rRNA gene amplicon sequencing analyses (median 2.8%), and up to 2.7 × 104 cells mL-1 of water (median 0.5 × 103) based on qPCR of pmoA genes. Methanotrophic activity and populations in the lake water were highest during fall turnover and remained high through the winter ice-covered period into spring turnover. They declined during summer stratification, especially in the epilimnion. Three methanotroph genera (Methylobacter, Methylovulum, and Methyloparacoccus) cycled seasonally, based on both relative and absolute abundance measurements. Methylobacter and Methylovulum populations peaked in winter/spring, when methane oxidation activity was psychrophilic. Methyloparacoccus populations increased in the water column through summer and fall, when methane oxidation was mesophilic, and also predominated in the underlying tailings sediment. Other, less abundant genera grew primarily during summer, possibly due to distinct CH4/O2 microniches created during thermal stratification. These data are consistent with temporal and spatial niche differentiation based on temperature, CH4 and O2. This pit lake displays methane cycling and methanotroph population dynamics similar to natural boreal lakes. IMPORTANCE The study examined methanotrophic bacteria in an industrial end pit lake, combining molecular DNA methods (both quantitative and descriptive) with biogeochemical measurements. The lake was sampled over 5 years, in all four seasons, as often as weekly, and included sub-ice samples. The resulting multiseason and multiyear data set is unique in its size and intensity, and allowed us to document clear and consistent seasonal patterns of growth and decline of three methanotroph genera (Methylobacter, Methylovulum, and Methyloparacoccus). Laboratory experiments suggested that one major control of this succession was niche partitioning based on temperature. The study helps to understand microbial dynamics in engineered end pit lakes, but we propose that the dynamics are typical of boreal stratified lakes and widely applicable in microbial ecology and limnology. Methane-oxidizing bacteria are important model organisms in microbial ecology and have implications for global climate change.

Dissecting the dominant hot spring microbial populations based on community-wide sampling at single-cell genomic resolution.

With advances in DNA sequencing and miniaturized molecular biology workflows, rapid and affordable sequencing of single-cell genomes has become a reality. Compared to 16S rRNA gene surveys and shotgun metagenomics, large-scale application of single-cell genomics to whole microbial communities provides an integrated snapshot of community composition and function, directly links mobile elements to their hosts, and enables analysis of population heterogeneity of the dominant community members. To that end, we sequenced nearly 500 single-cell genomes from a low diversity hot spring sediment sample from Dewar Creek, British Columbia, and compared this approach to 16S rRNA gene amplicon and shotgun metagenomics applied to the same sample. We found that the broad taxonomic profiles were similar across the three sequencing approaches, though several lineages were missing from the 16S rRNA gene amplicon dataset, likely the result of primer mismatches. At the functional level, we detected a large array of mobile genetic elements present in the single-cell genomes but absent from the corresponding same species metagenome-assembled genomes. Moreover, we performed a single-cell population genomic analysis of the three most abundant community members, revealing differences in population structure based on mutation and recombination profiles. While the average pairwise nucleotide identities were similar across the dominant species-level lineages, we observed differences in the extent of recombination between these dominant populations. Most intriguingly, the creek's Hydrogenobacter sp. population appeared to be so recombinogenic that it more closely resembled a sexual species than a clonally evolving microbe. Together, this work demonstrates that a randomized single-cell approach can be useful for the exploration of previously uncultivated microbes from community composition to population structure.

Genomic Analysis of the Yet-Uncultured Binatota Reveals Broad Methylotrophic, Alkane-Degradation, and Pigment Production Capacities.

The recent leveraging of genome-resolved metagenomics has generated an enormous number of genomes from novel uncultured microbial lineages yet left many clades undescribed. Here, we present a global analysis of genomes belonging to Binatota (UBP10), a globally distributed, yet-uncharacterized bacterial phylum. All orders in Binatota encoded the capacity for aerobic methylotrophy using methanol, methylamine, sulfomethanes, and chloromethanes as the substrates. Methylotrophy in Binatota was characterized by order-specific substrate degradation preferences, as well as extensive metabolic versatility, i.e., the utilization of diverse sets of genes, pathways, and combinations to achieve a specific metabolic goal. The genomes also encoded multiple alkane hydroxylases and monooxygenases, potentially enabling growth on a wide range of alkanes and fatty acids. Pigmentation is inferred from a complete pathway for carotenoids (lycopene, ß- and γ-carotenes, xanthins, chlorobactenes, and spheroidenes) production. Further, the majority of genes involved in bacteriochlorophyll a, c, and d biosynthesis were identified, although absence of key genes and failure to identify a photosynthetic reaction center preclude proposing phototrophic capacities. Analysis of 16S rRNA databases showed the preferences of Binatota to terrestrial and freshwater ecosystems, hydrocarbon-rich habitats, and sponges, supporting their potential role in mitigating methanol and methane emissions, breakdown of alkanes, and their association with sponges. Our results expand the lists of methylotrophic, aerobic alkane-degrading, and pigment-producing lineages. We also highlight the consistent encountering of incomplete biosynthetic pathways in microbial genomes, a phenomenon necessitating careful assessment when assigning putative functions based on a set-threshold of pathway completion.IMPORTANCE A wide range of microbial lineages remain uncultured, yet little is known regarding their metabolic capacities, physiological preferences, and ecological roles in various ecosystems. We conducted a thorough comparative genomic analysis of 108 genomes belonging to the Binatota (UBP10), a globally distributed, yet-uncharacterized bacterial phylum. We present evidence that members of the order Binatota specialize in methylotrophy and identify an extensive repertoire of genes and pathways mediating the oxidation of multiple one-carbon (C1) compounds in Binatota genomes. The occurrence of multiple alkane hydroxylases and monooxygenases in these genomes was also identified, potentially enabling growth on a wide range of alkanes and fatty acids. Pigmentation is inferred from a complete pathway for carotenoids production. We also report on the presence of incomplete chlorophyll biosynthetic pathways in all genomes and propose several evolutionary-grounded scenarios that could explain such a pattern. Assessment of the ecological distribution patterns of the Binatota indicates preference of its members to terrestrial and freshwater ecosystems characterized by high methane and methanol emissions, as well as multiple hydrocarbon-rich habitats and marine sponges.

GAL08, an Uncultivated Group of Acidobacteria, Is a Dominant Bacterial Clade in a Neutral Hot Spring.

GAL08 are bacteria belonging to an uncultivated phylogenetic cluster within the phylum Acidobacteria. We detected a natural population of the GAL08 clade in sediment from a pH-neutral hot spring located in British Columbia, Canada. To shed light on the abundance and genomic potential of this clade, we collected and analyzed hot spring sediment samples over a temperature range of 24.2-79.8°C. Illumina sequencing of 16S rRNA gene amplicons and qPCR using a primer set developed specifically to detect the GAL08 16S rRNA gene revealed that absolute and relative abundances of GAL08 peaked at 65°C along three temperature gradients. Analysis of sediment collected over multiple years and locations revealed that the GAL08 group was consistently a dominant clade, comprising up to 29.2% of the microbial community based on relative read abundance and up to 4.7 × 105 16S rRNA gene copy numbers per gram of sediment based on qPCR. Using a medium quality threshold, 25 single amplified genomes (SAGs) representing these bacteria were generated from samples taken at 65 and 77°C, and seven metagenome-assembled genomes (MAGs) were reconstructed from samples collected at 45-77°C. Based on average nucleotide identity (ANI), these SAGs and MAGs represented three separate species, with an estimated average genome size of 3.17 Mb and GC content of 62.8%. Phylogenetic trees constructed from 16S rRNA gene sequences and a set of 56 concatenated phylogenetic marker genes both placed the three GAL08 bacteria as a distinct subgroup of the phylum Acidobacteria, representing a candidate order (Ca. Frugalibacteriales) within the class Blastocatellia. Metabolic reconstructions from genome data predicted a heterotrophic metabolism, with potential capability for aerobic respiration, as well as incomplete denitrification and fermentation. In laboratory cultivation efforts, GAL08 counts based on qPCR declined rapidly under atmospheric levels of oxygen but increased slightly at 1% (v/v) O2, suggesting a microaerophilic lifestyle.

Ecological and genomic analyses of candidate phylum WPS-2 bacteria in an unvegetated soil.

Members of the bacterial candidate phylum WPS-2 (or Eremiobacterota) are abundant in several dry, bare soil environments. In a bare soil deposited by an extinct iron-sulfur spring, we found that WPS-2 comprised up to 24% of the bacterial community and up to 108 cells per g of soil based on 16S rRNA gene sequencing and quantification. A single genus-level cluster (Ca. Rubrimentiphilum) predominated in bare soils but was less abundant in adjacent forest. Nearly complete genomes of Ca. Rubrimentiphilum were recovered as single amplified genomes (SAGs) and metagenome-assembled genomes (MAGs). Surprisingly, given the abundance of WPS-2 in bare soils, the genomes did not indicate any capacity for autotrophy, phototrophy, or trace gas metabolism. Instead, they suggest a predominantly aerobic organoheterotrophic lifestyle, perhaps based on scavenging amino acids, nucleotides, and complex oligopeptides, along with lithotrophic capacity on thiosulfate. Network analyses of the entire community showed that some species of Chloroflexi, Actinobacteria, and candidate phylum AD3 (or Dormibacterota) co-occurred with Ca. Rubrimentiphilum and may represent ecological or metabolic partners. We propose that Ca. Rubrimentiphilum act as efficient heterotrophic scavengers. Combined with previous studies, these data suggest that the phylum WPS-2 includes bacteria with diverse metabolic capabilities.

Novel copper-containing membrane monooxygenases (CuMMOs) encoded by alkane-utilizing Betaproteobacteria.

Copper-containing membrane monooxygenases (CuMMOs) are encoded by xmoCAB(D) gene clusters and catalyze the oxidation of methane, ammonia, or some short-chain alkanes and alkenes. In a metagenome constructed from an oilsands tailings pond we detected an xmoCABD gene cluster with <59% derived protein sequence identity to genes from known bacteria. Stable isotope probing experiments combined with a specific xmoA qPCR assay demonstrated that the bacteria possessing these genes were incapable of methane assimilation, but did grow on ethane and propane. Single-cell amplified genomes (SAGs) from propane-enriched samples were screened with the specific PCR assay to identify bacteria possessing the target gene cluster. Multiple SAGs of Betaproteobacteria belonging to the genera Rhodoferax and Polaromonas possessed homologues of the metagenomic xmoCABD gene cluster. Unexpectedly, each of these two genera also possessed other xmoCABD paralogs, representing two additional lineages in phylogenetic analyses. Metabolic reconstructions from SAGs predicted that neither bacterium encoded enzymes with the potential to support catabolic methane or ammonia oxidation, but that both were capable of higher n-alkane degradation. The involvement of the encoded CuMMOs in alkane oxidation was further suggested by reverse transcription PCR analyses, which detected elevated transcription of the xmoA genes upon enrichment of water samples with propane as the sole energy source. Enrichments, isotope incorporation studies, genome reconstructions, and gene expression studies therefore all agreed that the unknown xmoCABD operons did not encode methane or ammonia monooxygenases, but rather n-alkane monooxygenases. This study broadens the known diversity of CuMMOs and identifies these enzymes in non-nitrifying Betaproteobacteria.

Stable Isotope and Metagenomic Profiling of a Methanogenic Naphthalene-Degrading Enrichment Culture.

Polycyclic aromatic hydrocarbons (PAH) such as naphthalene are widespread, recalcitrant pollutants in anoxic and methanogenic environments. A mechanism catalyzing PAH activation under methanogenic conditions has yet to be discovered, and the microbial communities coordinating their metabolism are largely unknown. This is primarily due to the difficulty of cultivating PAH degraders, requiring lengthy incubations to yield sufficient biomass for biochemical analysis. Here, we sought to characterize a new methanogenic naphthalene-degrading enrichment culture using DNA-based stable isotope probing (SIP) and metagenomic analyses. 16S rRNA gene sequencing of fractionated DNA pinpointed an unclassified Clostridiaceae species as a putative naphthalene degrader after two months of SIP incubation. This finding was supported by metabolite and metagenomic evidence of genes predicted to encode for enzymes facilitating naphthalene carboxylic acid CoA-thioesterification and degradation of an unknown arylcarboxyl-CoA structure. Our findings also suggest a possible but unknown role for Desulfuromonadales in naphthalene degradation. This is the first reported functional evidence of PAH biodegradation by a methanogenic consortium, and we envision that this approach could be used to assess carbon flow through other slow growing enrichment cultures and environmental samples.

Benzene and Naphthalene Degrading Bacterial Communities in an Oil Sands Tailings Pond.

Oil sands process-affected water (OSPW), produced by surface-mining of oil sands in Canada, is alkaline and contains high concentrations of salts, metals, naphthenic acids, and polycyclic aromatic compounds (PAHs). Residual hydrocarbon biodegradation occurs naturally, but little is known about the hydrocarbon-degrading microbial communities present in OSPW. In this study, aerobic oxidation of benzene and naphthalene in the surface layer of an oil sands tailings pond were measured. The potential oxidation rates were 4.3 µmol L-1 OSPW d-1 for benzene and 21.4 µmol L-1 OSPW d-1 for naphthalene. To identify benzene and naphthalene-degrading microbial communities, metagenomics was combined with stable isotope probing (SIP), high-throughput sequencing of 16S rRNA gene amplicons, and isolation of microbial strains. SIP using 13C-benzene and 13C-naphthalene detected strains of the genera Methyloversatilis and Zavarzinia as the main benzene degraders, while strains belonging to the family Chromatiaceae and the genus Thauera were the main naphthalene degraders. Metagenomic analysis revealed a diversity of genes encoding oxygenases active against aromatic compounds. Although these genes apparently belonged to many phylogenetically diverse taxa, only a few of these taxa were predominant in the SIP experiments. This suggested that many members of the community are adapted to consuming other aromatic compounds, or are active only under specific conditions. 16S rRNA gene sequence datasets have been submitted to the Sequence Read Archive (SRA) under accession number SRP109130. The Gold Study and Project submission ID number in Joint Genome Institute IMG/M for the metagenome is Gs0047444 and Gp0055765.