Re-evaluation of dioxygenase gene phylogeny for the development and validation of a quantitative assay for environmental aromatic hydrocarbon degraders

Differential effect of monoterpenes and flavonoids on the transcription of aromatic ring-hydroxylating dioxygenase genes in Rhodococcus opacus C1 and Rhodococcus sp. WAY2

Aromatic ring-hydroxylating dioxygenases (ARHDs) play a crucial role in the aerobic biodegradation of both natural and anthropogenic aromatic compounds. Although their ability to process contaminants is not entirely understood, it is thought to have evolved from the transformation of structurally similar secondary plant metabolites (SPMs). Hence, to investigate this connection, we tested a variety of SPMs from the monoterpene and flavonoid classes as carbon sources and transcriptional effectors of several phylogenetically distant ARHD genes involved in the degradation of aromatic pollutants. Specifically, we focused on bphA1, nahA1 and phtA1 in Rhodococcus opacus C1, whose genomic analysis is also presented hereinafter, and bphA1a, nahA1-bphA1b and etbA1ab in Rhodococcus sp. WAY2. Whilst induction was only observed with (R)-carvone for bphA1a and nahA1-bphA1b of strain WAY2, and with p-cymene for nahA1 and nahA1-bphA1b of strains C1 and WAY2, respectively, an extensive inhibition by flavonoids was observed for most of the genes in both strains. To the best of our knowledge, our study is the first to report the effect of flavonoids and monoterpenes on the transcription of nahA1, etbA1 and phtA1 genes. In addition, we show that, in contrast to pseudomonads, many flavonoids inhibit the transcription of the ARHD genes in rhodococci. Thus, our work provides a new perspective on flavonoids as the transcriptional effectors of ARHDs, highlighting the significant variability of these enzymes and the divergent responses that they elicit. Moreover, our results contribute to understanding the complex interactions between microorganisms and SPMs and provide insights into the molecular basis of a number of them.

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.

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.

DARHD: A sequence database for aromatic ring-hydroxylating dioxygenase analysis and primer evaluation

Biodegradation of aromatic compounds is ubiquitous in the environment and important for controlling organic pollutants. Aromatic ring-hydroxylating dioxygenases (ARHDs) are responsible for the first and rate-limiting step of aerobic biodegradation of aromatic compounds. The ARHD α subunit is a good biomarker for studying functional microorganisms in the environment, however their diversity and corresponding primer coverage are unclear, both of which require a comprehensive sequence database for the ARHD α subunit. Here amino acid sequences of the ARHD α subunit were collected, and a total of 103 sequences were selected as seed sequences that were distributed in 72 bacterial genera with 34 gene names. Based on both homolog search and keyword confirmation against the GenBank, a sequence database of ARHD (DARHD) has been established and 6367 highly credible sequences were retrieved. DARHD contained 407 bacterial genera capable of degrading 38 aromatic substrates, and intricate relationships among the gene name, aromatic substrate and microbial taxa were observed. Thereafter, a total of 136 pairs of primers were collected and assessed. Results showed coverages of most published primers were low. Our research provides new insights for understanding the diversity of ARHD α subunit, and gives guidance on the design and application of primers in the future.

Abundance and diversity of n-alkane and PAH-degrading bacteria and their functional genes - Potential for use in detection of marine oil pollution

Monitoring environmental status through molecular investigation of microorganisms in the marine environment is suggested as a potentially very effective method for biomonitoring, with great potential for automation. There are several hurdles to that approach with regards to primer design, variability across geographical locations, seasons, and type of environmental pollution. Here, qPCR analysis of genes involved in the initial activation of aliphatic and aromatic hydrocarbons were used in a laboratory setup mimicking realistic oil leakage at sea. Seawater incubation experiments were carried out under two different seasons with two different oil types. Degenerate primers targeting initial oxygenases (alkane 1-monooxygenase; alkB and aromatic-ring hydroxylating dioxygenase; ARHD) were employed in qPCR assays to quantify the abundance of genes essential for oil degradation. Shotgun metagenomics was used to map the overall community dynamics and the diversity of alkB and ARHD genes represented in the microbial community. The amplicons generated through the qPCR assays were sequenced to reveal the diversity of oil-degradation related genes captured by the degenerate primers. We identified a major mismatch between the taxonomic diversity of alkB and ARHD genes amplified by the degenerate primers and those identified through shotgun metagenomics. More specifically, the designed primers did not amplify the alkB genes of the two most abundant alkane degraders that bloomed in the experiments, Oceanobacter and Oleispira. The relative abundance of alkB sequences from shotgun metagenomics and 16S rRNA-based Oleispira-specific qPCR assay were better signals for oil in water than the tested qPCR alkB assay. The ARHD assay showed a good agreement with PAHs degradation despite covering only 25% of the top 100 ARHD genes and missing several abundant Cycloclasticus sequences that were present in the metagenome. We conclude that further improvement of the degenerate primer approach is needed to rely on the use of oxygenase-related qPCR assays for oil leakage detection.

Experimental and Genomic Evaluation of the Oestrogen Degrading Bacterium Rhodococcus equi ATCC13557

Rhodococcus equi ATCC13557 was selected as a model organism to study oestrogen degradation based on its previous ability to degrade 17α-ethinylestradiol (EE2). Biodegradation experiments revealed that R. equi ATCC13557 was unable to metabolise EE2. However, it was able to metabolise E2 with the major metabolite being E1 with no further degradation of E1. However, the conversion of E2 into E1 was incomplete, with 11.2 and 50.6% of E2 degraded in mixed (E1-E2-EE2) and E2-only conditions, respectively. Therefore, the metabolic pathway of E2 degradation by R. equi ATCC13557 may have two possible pathways. The genome of R. equi ATCC13557 was sequenced, assembled, and mapped for the first time. The genome analysis allowed the identification of genes possibly responsible for the observed biodegradation characteristics of R. equi ATCC13557. Several genes within R. equi ATCC13557 are similar, but not identical in sequence, to those identified within the genomes of other oestrogen degrading bacteria, including Pseudomonas putida strain SJTE-1 and Sphingomonas strain KC8. Homologous gene sequences coding for enzymes potentially involved in oestrogen degradation, most commonly a cytochrome P450 monooxygenase (oecB), extradiol dioxygenase (oecC), and 17β-hydroxysteroid dehydrogenase (oecA), were identified within the genome of R. equi ATCC13557. These searches also revealed a gene cluster potentially coding for enzymes involved in steroid/oestrogen degradation; 3-carboxyethylcatechol 2,3-dioxygenase, 2-hydroxymuconic semialdehyde hydrolase, 3-alpha-(or 20-beta)-hydroxysteroid dehydrogenase, 3-(3-hydroxy-phenyl)propionate hydroxylase, cytochrome P450 monooxygenase, and 3-oxosteroid 1-dehydrogenase. Further, the searches revealed steroid hormone metabolism gene clusters from the 9, 10-seco pathway, therefore R. equi ATCC13557 also has the potential to metabolise other steroid hormones such as cholesterol.

Comparative Genomics of the Rhodococcus Genus Shows Wide Distribution of Biodegradation Traits

The genus Rhodococcus exhibits great potential for bioremediation applications due to its huge metabolic diversity, including biotransformation of aromatic and aliphatic compounds. Comparative genomic studies of this genus are limited to a small number of genomes, while the high number of sequenced strains to date could provide more information about the Rhodococcus diversity. Phylogenomic analysis of 327 Rhodococcus genomes and clustering of intergenomic distances identified 42 phylogenomic groups and 83 species-level clusters. Rarefaction models show that these numbers are likely to increase as new Rhodococcus strains are sequenced. The Rhodococcus genus possesses a small “hard” core genome consisting of 381 orthologous groups (OGs), while a “soft” core genome of 1253 OGs is reached with 99.16% of the genomes. Models of sequentially randomly added genomes show that a small number of genomes are enough to explain most of the shared diversity of the Rhodococcus strains, while the “open” pangenome and strain-specific genome evidence that the diversity of the genus will increase, as new genomes still add more OGs to the whole genomic set. Most rhodococci possess genes involved in the degradation of aliphatic and aromatic compounds, while short-chain alkane degradation is restricted to a certain number of groups, among which a specific particulate methane monooxygenase (pMMO) is only found in Rhodococcus sp. WAY2. The analysis of Rieske 2Fe-2S dioxygenases among rhodococci genomes revealed that most of these enzymes remain uncharacterized.

Plant-microorganisms interaction promotes removal of air pollutants in Milan (Italy) urban area

Plants and phyllosphere microorganisms may effectively contribute to reducing air pollution in cities through the adsorption and biodegradation of pollutants onto leaves. In this work, during all seasons, we sampled atmospheric particulate matter (PM₁₀) and leaves of southern magnolia Magnolia grandiflora and deodar cedar Cedrus deodara, two evergreen plant species widespread in the urban area of Milan where the study was carried out. We then quantified Polycyclic Aromatic Hydrocarbons (PAHs) both in PM₁₀ and on leaves and used sequencing of 16S rRNA gene, shotgun metagenomics and qPCR analyses to investigate the microbial communities hosted by the sampled leaves. Taxonomic and functional profiles of epiphytic bacterial communities differed between host plant species and seasons and the microbial communities on leaves harboured genes involved in the degradation of hydrocarbons. Evidence collected in this work also suggested that the abundance of hydrocarbon-degrading microorganisms on evergreen leaves increased with the concentration of hydrocarbons when atmospheric pollutants were deposited at high concentration on leaves, and that the biodegradation on the phyllosphere can contribute to the removal of PAHs from the urban air.

pahE, a Functional Marker Gene for Polycyclic Aromatic Hydrocarbon-Degrading Bacteria

The characterization of native polycyclic aromatic hydrocarbon (PAH)-degrading bacteria is significant for understanding the PAH degradation process in the natural environment and developing effective remediation technologies. Most previous investigations of PAH-degrading bacteria in environmental samples employ pahAc, which encodes the α-subunit of PAH ring-hydroxylating dioxygenase, as a functional marker gene. However, the poor phylogenetic resolution and nonspecificity of pahAc result in a misestimation of PAH-degrading bacteria. Here, we propose a PAH hydratase-aldolase-encoding gene, pahE, as a superior biomarker for PAH-degrading bacteria. Comparative phylogenetic analysis of the key enzymes involved in the upper pathway of PAH degradation indicated that pahE evolved dependently from a common ancestor. A phylogenetic tree constructed based on PahE is largely congruent with PahAc-based phylogenies, except for the dispersion of several clades of other non-PAH-degrading aromatic hydrocarbon dioxygenases present in the PahAc tree. Analysis of pure strains by PCR confirmed that pahE can specifically distinguish PAH-degrading bacteria, while pahAc cannot. Illumina sequencing of pahE and pahAc amplicons showed more genotypes and higher specificity and resolution for pahE Novel reads were also discovered among the pahE amplicons, suggesting the presence of novel PAH-degrading populations. These results suggest that pahE is a more powerful biomarker for exploring the ecological role and degradation potential of PAH-degrading bacteria in ecosystems, which is significant to the bioremediation of PAH pollution and environmental microbial ecology.IMPORTANCE PAH contamination has become a worldwide environmental issue because of the potential toxic effects on natural ecosystems and human health. Biotransformation and biodegradation are considered the main natural elimination forms of PAHs from contaminated sites. Therefore, the knowledge of the degradation potential of the microbial community in contaminated sites is crucial for PAH pollution bioremediation. However, the nonspecificity of pahAc as a functional marker of PAH-degrading bacteria has resulted neither in a reliable prediction of PAH degradation potential nor an accurate assessment of degradation. Here, we introduced pahE encoding the PAH hydratase-aldolase as a new and better functional marker gene of PAH-degrading bacteria. This study provides a powerful molecular tool to more effectively explore the ecological role and degradation potential of PAH-degrading bacteria in ecosystems, which is significant to the bioremediation of PAH pollution.

Diversity of aromatic hydroxylating dioxygenase genes in mangrove microbiome and their biogeographic patterns across global sites

Aromatic hydrocarbons (AH), such as polycyclic aromatic hydrocarbons, are compounds largely found in nature. Aromatic‐ring‐hydroxylating dioxygenases (ARHD) are proteins involved in AH degradation pathways. We used ARHD functional genes from an oil‐impacted mangrove area and compared their diversity with other sites around the world to understand the ARHD biogeographic distribution patterns. For this, a comprehensive database was established with 166 operational protein families (OPFs) from 1,758 gene sequences obtained from 15 different sites worldwide, of which twelve are already published studies and three are unpublished. Based on a deduced ARHD peptide sequences consensus phylogeny, we examined trends and divergences in the sequence phylogenetic clustering from the different sites. The taxonomic affiliation of the OPF revealed that Pseudomonas, Streptomyces, Variovorax, Bordetella and Rhodococcus were the five most abundant genera, considering all sites. The functional diversity analysis showed the enzymatic prevalence of benzene 1,2‐dioxygenase, 3‐phenylpropionate dioxygenase and naphthalene 1,2‐dioxygenase, in addition to 10.98% of undefined category ARHDs. The ARHD gene correlation analysis among different sites was essentially important to gain insights on spatial distribution patterns, genetic congruence and ecological coherence of the bacterial groups found. This work revealed the genetic potential from the mangrove sediment for AH biodegradation and a considerable evolutionary proximity among the dioxygenase OPFs found in Antarctica and South America sites, in addition to high level of endemism in each continental region.

Benzene Degradation by a Variovorax Species within a Coal Tar-Contaminated Groundwater Microbial Community

Investigations of environmental microbial communities are crucial for the discovery of populations capable of degrading hazardous compounds and may lead to improved bioremediation strategies. The goal of this study was to identify microorganisms responsible for aerobic benzene degradation in coal tar-contaminated groundwater. Benzene degradation was monitored in laboratory incubations of well waters using gas chromatography mass spectrometry (GC-MS). Stable isotope probing (SIP) experiments using [¹³C]benzene enabled us to obtain ¹³C-labled community DNA. From this, 16S rRNA clone libraries identified Gammaproteobacteria and Betaproteobacteria as the active benzene-metabolizing microbial populations. Subsequent cultivation experiments yielded nine bacterial isolates that grew in the presence of benzene; five were confirmed in laboratory cultures to grow on benzene. The isolated benzene-degrading organisms were genotypically similar (>97% 16S rRNA gene nucleotide identities) to the organisms identified in SIP experiments. One isolate, Variovorax MAK3, was further investigated for the expression of a putative aromatic ring-hydroxylating dioxygenase (RHD) hypothesized to be involved in benzene degradation. Microcosm experiments using Variovorax MAK3 revealed a 10-fold increase in RHD (Vapar_5383) expression, establishing a link between this gene and benzene degradation. Furthermore, the addition of Variovorax MAK3 to microcosms prepared from site waters accelerated community benzene degradation and correspondingly increased RHD gene expression. In microcosms using uninoculated groundwater, quantitative (q)PCR assays (with 16S rRNA and RDH genes) showed that Variovorax was present and responsive to added benzene. These data demonstrate how the convergence of cultivation-dependent and -independent techniques can boost understandings of active populations and functional genes in complex benzene-degrading microbial communities.

IMPORTANCE

Benzene is a human carcinogen whose presence in contaminated groundwater drives environmental cleanup efforts. Although the aerobic biodegradation of benzene has long been established, knowledge of the identity of the microorganisms in complex naturally occurring microbial communities responsible for benzene biodegradation has evaded scientific inquiry for many decades. Here, we applied a molecular biology technique known as stable isotope probing (SIP) to the microbial communities residing in contaminated groundwater samples to identify the community members active in benzene biodegradation. We complemented this approach by isolating and growing in the laboratory a bacterium representative of the bacteria found using SIP. Further characterization of the isolated bacterium enabled us to track the expression of a key gene that attacks benzene both in pure cultures of the bacterium and in the naturally occurring groundwater microbial community. This work advances information regarding the documentation of microbial processes, especially the populations and genes that contribute to bioremediation.

In Silico Identification of Bioremediation Potential: Carbamazepine and Other Recalcitrant Personal Care Products

Emerging contaminants are principally personal care products not readily removed by conventional wastewater treatment and, with an increasing reliance on water recycling, become disseminated in drinking water supplies. Carbamazepine, a widely used neuroactive pharmaceutical, increasingly escapes wastewater treatment and is found in potable water. In this study, a mechanism is proposed by which carbamazepine resists biodegradation, and a previously unknown microbial biodegradation was predicted computationally. The prediction identified biphenyl dioxygenase from Paraburkholderia xenovorans LB400 as the best candidate enzyme for metabolizing carbamazepine. The rate of degradation described here is 40 times greater than the best reported rates. The metabolites cis-10,11-dihydroxy-10,11-dihydrocarbamazepine and cis-2,3-dihydroxy-2,3-dihydrocarbamazepine were demonstrated with the native organism and a recombinant host. The metabolites are considered nonharmful and mitigate the generation of carcinogenic acridine products known to form when advanced oxidation methods are used in water treatment. Other recalcitrant personal care products were subjected to prediction by the Pathway Prediction System and tested experimentally with P. xenovorans LB400. It was shown to biodegrade structurally diverse compounds. Predictions indicated hydrolase or oxygenase enzymes catalyzed the initial reactions. This study highlights the potential for using the growing body of enzyme-structural and genomic information with computational methods to rapidly identify enzymes and microorganisms that biodegrade emerging contaminants.

Role of environmental factors and microorganisms in determining the fate of polycyclic aromatic hydrocarbons in the marine environment

Polycyclic aromatic hydrocarbons (PAHs) are widespread in marine ecosystems and originate from natural sources and anthropogenic activities. PAHs enter the marine environment in two main ways, corresponding to chronic pollution or acute pollution by oil spills. The global PAH fluxes in marine environments are controlled by the microbial degradation and the biological pump, which plays a role in particle settling and in sequestration through bioaccumulation. Due to their low water solubility and hydrophobic nature, PAHs tightly adhere to sediments leading to accumulation in coastal and deep sediments. Microbial assemblages play an important role in determining the fate of PAHs in water and sediments, supporting the functioning of biogeochemical cycles and the microbial loop. This review summarises the knowledge recently acquired in terms of both chronic and acute PAH pollution. The importance of the microbial ecology in PAH-polluted marine ecosystems is highlighted as well as the importance of gaining further in-depth knowledge of the environmental services provided by microorganisms.