Metagenomics and Quantitative Stable Isotope Probing Offer Insights into Metabolism of Polycyclic Aromatic Hydrocarbon Degraders in Chronically Polluted Seawater

Genetic redundancy in the naphthalene-degradation pathway of Cycloclasticus pugetii strain PS-1 enables response to varying substrate concentrations

Polycyclic aromatic hydrocarbon (PAH) contamination in marine environments range from low-diffusive inputs to high loads. The influence of PAH concentration on the expression of functional genes [e.g. those encoding ring-hydroxylating dioxygenases (RHDs)] has been overlooked in PAH biodegradation studies. However, understanding marker-gene expression under different PAH loads can help to monitor and predict bioremediation efficiency. Here, we followed the expression (via RNA sequencing) of Cycloclasticus pugetii strain PS-1 in cell suspension experiments under different naphthalene (100 and 30 mg L-1) concentrations. We identified genes encoding previously uncharacterized RHD subunits, termed rhdPS1α and rhdPS1β, that were highly transcribed in response to naphthalene-degradation activity. Additionally, we identified six RHD subunit-encoding genes that responded to naphthalene exposure. By contrast, four RHD subunit genes were PAH-independently expressed and three other RHD subunit genes responded to naphthalene starvation. Cycloclasticus spp. could, therefore, use genetic redundancy in key PAH-degradation genes to react to varying PAH loads. This genetic redundancy may restrict the monitoring of environmental hydrocarbon-degradation activity using single-gene expression. For Cycloclasticus pugetii strain PS-1, however, the newly identified rhdPS1α and rhdPS1β genes might be potential target genes to monitor its environmental naphthalene-degradation activity.

Utilization of lasso peptides for biodegradation of polycyclic aromatic hydrocarbons

Many microbial genes involved in degrading recalcitrant environmental contaminants such as polycyclic aromatic hydrocarbons (PAHs) have been identified and characterized. However, all molecular mechanisms required for PAH utilization have not yet been elucidated. In this work, we demonstrate the proposed involvement of lasso peptides in the utilization of the PAH phenanthrene in Sphingomonas BPH. Transpositional mutagenesis of Sphingomonas BPH with the miniTn5 transposon yielded 3 phenanthrene utilization deficient mutants, #257, #1778, and #1782. In mutant #1782, Tn5 had inserted into the large subunit of the naph/bph dioxygenase gene. In mutant #1778, Tn5 had inserted into the B2 protease gene of a lasso peptide cluster. This finding is the first report on the role of lasso peptides in PAH utilization. Our studies also demonstrate that interruption of the lasso peptide cluster resulted in a significant increase in the amount of biosurfactant produced in the presence of glucose when compared to the wild-type strain. Collectively, these results suggest that the mechanisms Sphingomonas BPH utilizes to degrade phenanthrene are far more complex than previously understood and that the #1778 mutant may be a good candidate for bioremediation when glucose is applied as an amendment due to its higher biosurfactant production.

Evidence for novel polycyclic aromatic hydrocarbon degradation pathways in culturable marine isolates

IMPORTANCE

Polycyclic aromatic hydrocarbon (PAH) pollution is widespread throughout marine environments and significantly affects native flora and fauna. Investigating microbes responsible for degrading PAHs in these environments provides a greater understanding of natural attenuation in these systems. In addition, the use of culture-based approaches to inform bioinformatic and omics-based approaches is useful in identifying novel mechanisms of PAH degradation that elude genetic biomarker-based investigations. Furthermore, culture-based approaches allow for the study of PAH co-metabolism, which increasingly appears to be a prominent mechanism for PAH degradation in marine microbes.

Exploring novel alkane-degradation pathways in uncultured bacteria from the North Atlantic Ocean

IMPORTANCE

Petroleum pollution in the ocean has increased because of rapid population growth and modernization, requiring urgent remediation. Our understanding of the metabolic response of native microbial communities to oil spills is not well understood. Here, we explored the baseline hydrocarbon-degrading communities of a subarctic Atlantic region to uncover the metabolic potential of the bacteria that inhabit the surface and subsurface water. We conducted enrichments with a 13C-labeled hydrocarbon to capture the fraction of the community actively using the hydrocarbon. We then combined this approach with metagenomics to identify the metabolic potential of this hydrocarbon-degrading community. This revealed previously undescribed uncultured bacteria with unique metabolic mechanisms involved in aerobic hydrocarbon degradation, indicating that temperature may be pivotal in structuring hydrocarbon-degrading baseline communities. Our findings highlight gaps in our understanding of the metabolic complexity of hydrocarbon degradation by native marine microbial communities.

Bacterial cell sensing and signaling pathway for external polycyclic aromatic hydrocarbons (PAHs)

The mechanism by which a bacterial cell senses external nutrients remains largely unknown. In this study, we identified a bacterial cell sensing system for polycyclic aromatic hydrocarbons (PAHs) in a common marine PAH-using bacterium, Cycloclasticus. It consists of an outer membrane receptor (PahS) and a periplasmic protein (PahP) in combination with a two-component sensing system (TCS) that ensures a rapid response to PAH occurrence by directly controlling serial reactions including chemotactic sensing and movement, PAH uptake and intracellular PAH metabolism. PahS protrudes from the cell and acts as a PAH sensor, transducing the PAH signal across the outer membrane to its periplasmic partner PahP, which in turn transduces the PAH signal across the periplasm to a specialized TCS. This sensing system plays a critical role in sensing and promoting the metabolism of PAHs, which can be scavenged by various hydrocarbon-degrading bacteria.

Advances and perspectives of using stable isotope probing (SIP)-based technologies in contaminant biodegradation

Stable isotope probing (SIP) is a powerful tool to study microbial community structure and function in both nature and engineered environments. Coupling with advanced genomics and other techniques, SIP studies have generated substantial information to allow researchers to draw a clearer picture of what is occurring in complex microbial ecosystems. This review provides an overview of the advances of SIP-based technologies over time, summarizes the status of SIP applications to contaminant biodegradation, provides critical perspectives on ecological interactions within the community, and important factors (controllable and non-controllable) to be considered in SIP experimental designs and data interpretation. Current trend and perspectives of adapting SIP techniques for environmental applications are also discussed.

A standardized quantitative analysis strategy for stable isotope probing metagenomics

Stable isotope probing (SIP) facilitates culture-independent identification of active microbial populations within complex ecosystems through isotopic enrichment of nucleic acids. Many DNA-SIP studies rely on 16S rRNA gene sequences to identify active taxa, but connecting these sequences to specific bacterial genomes is often challenging. Here, we describe a standardized laboratory and analysis framework to quantify isotopic enrichment on a per-genome basis using shotgun metagenomics instead of 16S rRNA gene sequencing. To develop this framework, we explored various sample processing and analysis approaches using a designed microbiome where the identity of labeled genomes and their level of isotopic enrichment were experimentally controlled. With this ground truth dataset, we empirically assessed the accuracy of different analytical models for identifying active taxa and examined how sequencing depth impacts the detection of isotopically labeled genomes. We also demonstrate that using synthetic DNA internal standards to measure absolute genome abundances in SIP density fractions improves estimates of isotopic enrichment. In addition, our study illustrates the utility of internal standards to reveal anomalies in sample handling that could negatively impact SIP metagenomic analyses if left undetected. Finally, we present SIPmg, an R package to facilitate the estimation of absolute abundances and perform statistical analyses for identifying labeled genomes within SIP metagenomic data. This experimentally validated analysis framework strengthens the foundation of DNA-SIP metagenomics as a tool for accurately measuring the in situ activity of environmental microbial populations and assessing their genomic potential. IMPORTANCE Answering the questions, "who is eating what?" and "who is active?" within complex microbial communities is paramount for our ability to model, predict, and modulate microbiomes for improved human and planetary health. These questions can be pursued using stable isotope probing to track the incorporation of labeled compounds into cellular DNA during microbial growth. However, with traditional stable isotope methods, it is challenging to establish links between an active microorganism's taxonomic identity and genome composition while providing quantitative estimates of the microorganism's isotope incorporation rate. Here, we report an experimental and analytical workflow that lays the foundation for improved detection of metabolically active microorganisms and better quantitative estimates of genome-resolved isotope incorporation, which can be used to further refine ecosystem-scale models for carbon and nutrient fluxes within microbiomes.

Substrate-independent expression of key functional genes in Cycloclasticus pugetii strain PS-1 limits their use as markers for PAH biodegradation

Microbial degradation of petroleum hydrocarbons is a crucial process for the clean-up of oil-contaminated environments. Cycloclasticus spp. are well-known polycyclic aromatic hydrocarbon (PAH) degraders that possess PAH-degradation marker genes including rhd3α, rhd2α, and pahE. However, it remains unknown if the expression of these genes can serve as an indicator for active PAH degradation. Here, we determined transcript-to-gene (TtG) ratios with (reverse transcription) qPCR in cultures of Cycloclasticus pugetii strain PS-1 grown with naphthalene, phenanthrene, a mixture of these PAHs, or alternate substrates (i.e., no PAHs). Mean TtG ratios of 1.99 × 10^-2, 1.80 × 10^-3, and 3.20 × 10^-3 for rhd3α, rhd2α, and pahE, respectively, were measured in the presence or absence of PAHs. The TtG values suggested that marker-gene expression is independent of PAH degradation. Measurement of TtG ratios in Arctic seawater microcosms amended with water-accommodated crude oil fractions, and incubated under in situ temperature conditions (i.e., 1.5°C), only detected Cycloclasticus spp. rhd2α genes and transcripts (mean TtG ratio of 4.15 × 10^-1). The other marker genes-rhd3α and pahE-were not detected, suggesting that not all Cycloclasticus spp. carry these genes and a broader yet-to-be-identified repertoire of PAH-degradation genes exists. The results indicate that the expression of PAH marker genes may not correlate with PAH-degradation activity, and transcription data should be interpreted cautiously.

Seasonal variation of temperature affects HMW-PAH accumulation in fishery species by bacterially mediated LMW-PAH degradation

Currently, the specific mechanism generating seasonal variation in polycyclic aromatic hydrocarbons (PAHs) via bacterial biodegradation remains unclear, and whether this alteration affects PAH bioaccumulation is unknown. Therefore, we performed a study between 2015 and 2020 to investigate the effects of seasonal variation on bacterial communities and PAH bioaccumulation in the Pearl River Estuary. Significantly high PAH concentrations in both aquatic and fishery species were determined in dry seasons (the mean ∑₁₆PAH concentration: water, 37.24 ng/L (2015), 30.83 ng/L (2020); fish, 51.01 ng/L (2015) and 72.60 ng/L (2020)) compared to wet seasons (the mean ∑₁₆PAH concentration: water, 22.38 ng/L (2015), 19.40 ng/L(2020); fish, 25.28 ng/L (2015) and 32.59 ng/L (2020)). Distinct differences in taxonomic and functional composition of bacterial communities related to biodegradation of low molecular weight PAHs (LMW-PAHs) were observed between seasons, and the concentrations of PAHs were negatively correlated with seasonal variation in temperature. Temperature-related specific bacterial taxa (e.g., Stenotrophomonas) directly or indirectly participated in LMW-PAH degradation via encoding PAH degradation enzymes (e.g., protocatechuate 4,5-dioxygenase) that subsequently led to bioaccumulation of high molecular weight PAHs (HMW-PAHs) in wild and fishery species due to LMW-PAHs in the water. Based on this alteration, the ecological risk posed by PAHs decreased in wet seasons, and an unbalanced spatio-temporal distribution of PAHs was observed in this estuary. These results suggest that seasonal variation of temperature affects HMW-PAH accumulation in fishery species via bacterially mediated LMW-PAH biodegradation.

Quantitative Stable-Isotope Probing (qSIP) with Metagenomics Links Microbial Physiology and Activity to Soil Moisture in Mediterranean-Climate Grassland Ecosystems

The growth and physiology of soil microorganisms, which play vital roles in biogeochemical cycling, are shaped by both current and historical soil environmental conditions. Here, we developed and applied a genome-resolved metagenomic implementation of quantitative stable isotope probing (qSIP) with an H₂¹⁸O labeling experiment to identify actively growing soil microorganisms and their genomic capacities. qSIP enabled measurement of taxon-specific growth because isotopic incorporation into microbial DNA requires production of new genome copies. We studied three Mediterranean grassland soils across a rainfall gradient to evaluate the hypothesis that historic precipitation levels are an important factor controlling trait selection. We used qSIP-informed genome-resolved metagenomics to resolve the active subset of soil community members and identify their characteristic ecophysiological traits. Higher year-round precipitation levels correlated with higher activity and growth rates of flagellar motile microorganisms. In addition to heavily isotopically labeled bacteria, we identified abundant isotope-labeled phages, suggesting phage-induced cell lysis likely contributed to necromass production at all three sites. Further, there was a positive correlation between phage activity and the activity of putative phage hosts. Contrary to our expectations, the capacity to decompose the diverse complex carbohydrates common in soil organic matter or oxidize methanol and carbon monoxide were broadly distributed across active and inactive bacteria in all three soils, implying that these traits are not highly selected for by historical precipitation. IMPORTANCE Soil moisture is a critical factor that strongly shapes the lifestyle of soil organisms by changing access to nutrients, controlling oxygen diffusion, and regulating the potential for mobility. We identified active microorganisms in three grassland soils with similar mineral contexts, yet different historic rainfall inputs, by adding water labeled with a stable isotope and tracking that isotope in DNA of growing microbes. By examining the genomes of active and inactive microorganisms, we identified functions that are enriched in growing organisms, and showed that different functions were selected for in different soils. Wetter soil had higher activity of motile organisms, but activity of pathways for degradation of soil organic carbon compounds, including simple carbon substrates, were comparable for all three soils. We identified many labeled, and thus active bacteriophages (viruses that infect bacteria), implying that the cells they killed contributed to soil organic matter. The activity of these bacteriophages was significantly correlated with activity of their hosts.

Site-tropism of streptococci in the oral microbiome

A detailed understanding of where bacteria localize is necessary to advance microbial ecology and microbiome-based therapeutics. The site-specialist hypothesis predicts that most microbes in the human oral cavity have a primary habitat type within the mouth where they are most abundant. We asked whether this hypothesis accurately describes the distribution of the members of the genus Streptococcus, a clinically relevant taxon that dominates most oral sites. Prior analysis of 16S rRNA gene sequencing data indicated that some oral Streptococcus clades are site-specialists while others may be generalists. However, within complex microbial populations composed of numerous closely related species and strains, such as the oral streptococci, genome-scale analysis is necessary to provide the resolution to discriminate closely related taxa with distinct functional roles. Here, we assess whether individual species within this genus are specialists using publicly available genomic sequence data that provide species-level resolution. We chose a set of high-quality representative genomes for human oral Streptococcus species. Onto these genomes, we mapped shotgun metagenomic sequencing reads from supragingival plaque, tongue dorsum, and other sites in the oral cavity. We found that every abundant Streptococcus species in the healthy human oral cavity showed strong site-tropism and that even closely related species such as S. mitis, S. oralis, and S. infantis specialized in different sites. These findings indicate that closely related bacteria can have distinct habitat distributions in the absence of dispersal limitation and under similar environmental conditions and immune regimes. Substantial overlap between the core genes of these three species suggests that site-specialization is determined by subtle differences in genomic content.

HT-SIP: a semi-automated stable isotope probing pipeline identifies cross-kingdom interactions in the hyphosphere of arbuscular mycorrhizal fungi

BACKGROUND

Linking the identity of wild microbes with their ecophysiological traits and environmental functions is a key ambition for microbial ecologists. Of many techniques that strive for this goal, Stable-isotope probing-SIP-remains among the most comprehensive for studying whole microbial communities in situ. In DNA-SIP, actively growing microorganisms that take up an isotopically heavy substrate build heavier DNA, which can be partitioned by density into multiple fractions and sequenced. However, SIP is relatively low throughput and requires significant hands-on labor. We designed and tested a semi-automated, high-throughput SIP (HT-SIP) pipeline to support well-replicated, temporally resolved amplicon and metagenomics experiments. We applied this pipeline to a soil microhabitat with significant ecological importance-the hyphosphere zone surrounding arbuscular mycorrhizal fungal (AMF) hyphae. AMF form symbiotic relationships with most plant species and play key roles in terrestrial nutrient and carbon cycling.

RESULTS

Our HT-SIP pipeline for fractionation, cleanup, and nucleic acid quantification of density gradients requires one-sixth of the hands-on labor compared to manual SIP and allows 16 samples to be processed simultaneously. Automated density fractionation increased the reproducibility of SIP gradients compared to manual fractionation, and we show adding a non-ionic detergent to the gradient buffer improved SIP DNA recovery. We applied HT-SIP to ¹³C-AMF hyphosphere DNA from a ¹³CO₂ plant labeling study and created metagenome-assembled genomes (MAGs) using high-resolution SIP metagenomics (14 metagenomes per gradient). SIP confirmed the AMF Rhizophagus intraradices and associated MAGs were highly enriched (10-33 atom% ¹³C), even though the soils' overall enrichment was low (1.8 atom% ¹³C). We assembled 212 ¹³C-hyphosphere MAGs; the hyphosphere taxa that assimilated the most AMF-derived ¹³C were from the phyla Myxococcota, Fibrobacterota, Verrucomicrobiota, and the ammonia-oxidizing archaeon genus Nitrososphaera.

CONCLUSIONS

Our semi-automated HT-SIP approach decreases operator time and improves reproducibility by targeting the most labor-intensive steps of SIP-fraction collection and cleanup. We illustrate this approach in a unique and understudied soil microhabitat-generating MAGs of actively growing microbes living in the AMF hyphosphere (without plant roots). The MAGs' phylogenetic composition and gene content suggest predation, decomposition, and ammonia oxidation may be key processes in hyphosphere nutrient cycling. Video Abstract.

Integrating Viral Metagenomics into an Ecological Framework

Viral metagenomics has expanded our knowledge of the ecology of uncultured viruses, within both environmental (e.g., terrestrial and aquatic) and host-associated (e.g., plants and animals, including humans) contexts. Here, we emphasize the implementation of an ecological framework in viral metagenomic studies to address questions in virology rarely considered ecological, which can change our perception of viruses and how they interact with their surroundings. An ecological framework explicitly considers diverse variants of viruses in populations that make up communities of interacting viruses, with ecosystem-level effects. It provides a structure for the study of the diversity, distributions, dynamics, and interactions of viruses with one another, hosts, and the ecosystem, including interactions with abiotic factors. An ecological framework in viral metagenomics stands poised to broadly expand our knowledge in basic and applied virology. We highlight specific fundamental research needs to capitalize on its potential and advance the field. Expected final online publication date for the Annual Review of Virology, Volume 8 is September 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.