PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of ambexSpringerOpen.comThis journalSubmit a manuscriptRegisterSpringerOpen.comAMB Express
 
AMB Express. 2012; 2: 18.
Published online Mar 27, 2012. doi:  10.1186/2191-0855-2-18
PMCID: PMC3366876
Phylogenetic and functional diversity of metagenomic libraries of phenol degrading sludge from petroleum refinery wastewater treatment system
Cynthia C Silva,corresponding author1 Helen Hayden,2 Tim Sawbridge,2 Pauline Mele,2 Ricardo H Kruger,3 Marili VN Rodrigues,4 Gustavo GL Costa,5 Ramon O Vidal,5 Maíra P Sousa,6 Ana Paula R Torres,6 Vânia MJ Santiago,6 and Valéria M Oliveira1
1Microbial Resources Division, Research Center for Chemistry, Biology and Agriculture (CPQBA), Campinas University - UNICAMP, CP 6171, CEP 13081-970, Campinas, SP, Brazil
2Department of Primary Industries, BioSciences Research Division, 1 Park Dr, Bundoora, Melbourne, Vic 3083, Australia
3Department of Enzymology, University of Brasilia (UnB), Brasilia, Brazil
4Division of Organic and Pharmaceutical Chemistry, Research Center for Chemistry, Biology and Agriculture (CPQBA), Campinas University, Campinas, Brazil
5Laboratory of Genomics and Expression, Department of Genetics and Evolution, Institute of Biology, Campinas University - UNICAMP, CEP 13083-970, Campinas, SP, Brazil
6PETROBRAS R&D Center, University City, Square 7, Ilha do Fundão, Rio de Janeiro, RJ, Brazil
corresponding authorCorresponding author.
Cynthia C Silva: ccanedosilva/at/yahoo.com.br; Helen Hayden: helen.hayden/at/dpi.vic.gov.au; Tim Sawbridge: Tim.Sawbridge/at/dpi.vic.gov.au; Pauline Mele: pauline.mele/at/dpi.vic.gov.au; Ricardo H Kruger: kruger/at/unb.br; Marili VN Rodrigues: marili/at/cpqba.unicamp.br; Gustavo GL Costa: glacerda/at/lge.ibi.unicamp.br; Ramon O Vidal: ramon.vidal/at/gmail.com; Maíra P Sousa: mpsousa/at/petrobras.com.br; Ana Paula R Torres: aptorres/at/petrobras.com.br; Vânia MJ Santiago: vaniamjs/at/petrobras.com.br; Valéria M Oliveira: vmaia/at/cpqba.unicamp.br
Received February 15, 2012; Accepted March 27, 2012.
In petrochemical refinery wastewater treatment plants (WWTP), different concentrations of pollutant compounds are received daily in the influent stream, including significant amounts of phenolic compounds, creating propitious conditions for the development of particular microorganisms that can rapidly adapt to such environment. In the present work, the microbial sludge from a refinery WWTP was enriched for phenol, cloned into fosmid vectors and pyrosequenced. The fosmid libraries yielded 13,200 clones and a comprehensive bioinformatic analysis of the sequence data set revealed a complex and diverse bacterial community in the phenol degrading sludge. The phylogenetic analyses using MEGAN in combination with RDP classifier showed a massive predominance of Proteobacteria, represented mostly by the genera Diaphorobacter, Pseudomonas, Thauera and Comamonas. The functional classification of phenol degrading sludge sequence data set generated by MG-RAST showed the wide metabolic diversity of the microbial sludge, with a high percentage of genes involved in the aerobic and anaerobic degradation of phenol and derivatives. In addition, genes related to the metabolism of many other organic and xenobiotic compounds, such as toluene, biphenyl, naphthalene and benzoate, were found. Results gathered herein demonstrated that the phenol degrading sludge has complex phylogenetic and functional diversities, showing the potential of such community to degrade several pollutant compounds. This microbiota is likely to represent a rich resource of versatile and unknown enzymes which may be exploited for biotechnological processes such as bioremediation.
Keywords: Membrane bioreactor, Fosmid library, Pyrosequencing, Microbial diversity, Metabolic profile
Wastewater treatment plants (WWTP) represent habitats of continuous change in chemical composition (Szczepanowski et al. 2008). In oil refineries, the wastewater influent differs daily in terms of concentrations and composition of pollutant compounds, including light fraction aliphatic and aromatic petroleum hydrocarbons, organochlorines originated from cooling liquids used in the industrial process (Stepnowski et al. 2002) and other compounds such as phenol, chlorides, sulphides, sodium hydroxide, ammonia and heavy metals (Braile 1979; Mariano 2001). Phenols and derivatives are prominent pollutants in these wastes. These compounds are widely used as raw materials in the petrochemical industry and in oil refineries, for example in the washing and conditioning of alkaline or acid products. The increasing presence of phenols in the environment represents a serious ecological problem due to toxicity hazard for living creatures, including micro-organisms (Ojumu et al. 2005; Barrios-Martinez et al. 2006). Besides, the presence of phenols reduces significantly the biological degradation of the other compounds.
Several processes are used to eliminate phenolic compounds from industrial wastewater, but the biological treatments have been preferred for large-scale removal. However, this is not an easy task because of the proper toxicity of phenol towards microorganisms (Barrios-Martinez et al. 2006). In this sense, the monitoring of the microbiota is very important for efficient performance of biological treatment systems.
Traditionally, the diversity of microbial communities has been accessed by means of cultivation-based techniques or optical microscopy (Henze et al. 1997; Ojumu et al. 2005; Chang et al. 2005). Although very useful for taxonomic, physiological and genetic studies, culture-based techniques are insufficient for a more precise characterization of the functional and phylogenetic diversity of microbial communities, since it is now well known that only a small fraction (0.1 to 10%) of the microbial diversity in nature can be recovered in the laboratory (Torsvik et al. 1990; Amann et al. 1995).
In the last decade molecular cultivation-independent techniques have allowed the access to yet uncultivated microorganisms in several environmental niches (Handelsman 2004; Sleator et al. 2008), providing significant insights into bacterial communities in wastewater treatment processes (Sercu et al. 2006; Miura et al. 2007; Zang et al. 2008; Silva et al. 2010a; Silva et al. 2010b). Additionally, with the development of metagenomic approaches, the discovery and exploration of new microbial groups and functions have been accelerated (Handelsman 2004; Steele et al. 2009).
Nonetheless, molecular studies of microbial communities using Sanger sequencing have been limited by the number of sequences that can be obtained (Zhang et al. 2009). Microbial ecology studies based on 16S rRNA gene libraries have generally shown an underestimated bacterial diversity (Zhou et al. 2001; Zhang et al. 2009; Silva et al. 2010b). The recent application of next-generation sequencing technologies, such as pyrosequencing, has allowed one to obtain a huge number of sequences, usually sufficient to reveal the complexity of a microbial community in a given sample.
There are a few studies reporting the use of pyrosequencing for the phylogenetic and functional analysis of microbial communities in sludge from wastewater treatment processes. Some authors used this tool to look for plasmids with antibiotic resistance in sewage wastewater (Szczepanowski et al. 2008; Schlüter et al. 2008). Other studies employed 454-FLX pyrosequencing to investigate the microbial community from activated sludge of a domestic sewage wastewater treatment plant (Sanapareddy et al. 2009) or for a comprehensive understanding of the pathogen richn ess and abundance in metagenomic DNA samples derived from sewage sludge, composting and agricultural soil (Billy et al. 2010).
In this work we used the pyrosequencing approach followed by bioinformatics analysis to carry out an investigation of metagenomic libraries of phenol degrading sludge derived from petroleum refinery WWTP, allowing us to have a snapshot of the phylogenetic and functional microbial diversity selected after phenol enrichment of sludge.
Sampling
Sludge samples were collected from two different refinery wastewater treatment plants of the petroleum producer Petrobras (Brazil). The first sludge sample, named MBR1 and originated from a laboratory-scale membrane bioreactor (2 L) (MBR), was performed with continued aeration and collected after a 30-day period of high phenolic load feeding (68.5 mg L-1), as previously described by Viero and collaborators (Viero et al. 2008).
The second sample, named MBR2, was collected from a pilot submerged aerobic membrane bioreactor, previously described by Silva et al. (2010b), that has been in continuous operation for 18 months and is part of the industrial WWTP of Petrobras Refinery (Brazil). However, a phenol acclimation step was done to increase the possibility of finding genes related with phenol degradation. The sludge sample MBR2 was acclimatized in batch-culture for a 30 day-period to 1.0 g L-1 of phenol (Merck, USA) and then used for the metagenomic library construction. The acclimation step was performed, in triplicate, using 2.0 g L-1 of sludge, which was centrifuged and placed into an Erlenmeyer flask containing 300 mL of an initial rich nutrient medium (2.75 g L-1 K2HPO4, 2.25 g L-1 KH2PO4, 0.1 g L-1 NaCl, 1.0 g L-1 (NH4)2SO4, 0.2 g L-1 MgCl2.6H2O, 0.01 g L-1 CaCl2 and 1 g L-1 yeast extract as carbon source). These flasks were incubated at ambient temperature at 150 rpm. The initial carbon source was gradually diminished and replaced with phenol, in the proportion of 0.5 g L-1 less of yeast extract for each 0.2 g L-1 increment of phenol, until the yeast extract was totally eliminated. The modification of the initial medium composition was performed every three days, where an aliquot of the previous culture medium (10%) was added as inoculum to a new culture medium containing more phenol. The acclimatizing process was monitored through phenol quantification using gas chromatography. The sludge was collected after the microorganisms were considered totally adapted to 1.0 g L-1 phenol, when 100% phenol was removed in less than 24 hours.
Nucleic acid extraction and metagenomic fosmid library construction
High molecular weight DNA extraction from sludge samples was carried out using the protocol previously described by Silva et al. (2010a). Fosmid libraries were constructed using the CopyControl™ HTP Fosmid Library Production Kit (Epicentre, USA), according to the manufacturer's instructions. The DNA fragment size of interest was selected prior to the construction of the metagenomic libraries. The extracted DNA was run on a 1% low-melting agarose gel (Sigma, USA) and submitted to pulsed field gel electrophoresis (Pulsed-field CHEF DRIII System - BioRad- USA) at 0.5 s switch time, 9 Vcm-1, 120° included angle, 5 h at 14°C. DNA fragments of about 25-50 kb were isolated from the agarose gel, submitted to the end-repair reaction and ligated into the pCC2Fos fosmid vectors. These vectors were packed and transfected into Escherichia coli EPI300TM-T1R cells, which were plated onto Luria-Bertani (LB) agar medium containing chloramphenicol (12.5 μg/mL) and incubated overnight at 37°C. Two metagenomic libraries were constructed, one for each sludge sample (MBR1 and MBR2).
For the validation of the metagenomic libraries, six fosmid clones were randomly selected from each library. The fosmid DNA of each clone was extracted using the FosmidMax DNA Purification kit (Epicentre, USA), according to the manufacturer's protocol, and then digested using 10 U NotI (Promega, USA) at 37°C overnight. The band profiles of fosmid clones were checked in preparative pulsed field gel electrophoresis (Pulsed-field CHEF DRIII System - BioRad- USA) at angle 120°, 6 Vcm-1, 1 s - 12 s switch time, 10.5 h at 14°C.
Extraction of fosmid DNA pools and pyrosequencing
Clones derived from both metagenomic libraries were grown in 96-well plates containing 1 mL of LB medium and chloramphenicol (12.5 μg/mL) for 17 h at 37°C at 180 rpm. For fosmid extractions, an aliquot of 500 μL was aliquoted from each well of ten 96 well-plates and pooled, totaling a final volume of approximately 500 mL for each extraction. Each pool, containing 960 clones, was extracted using the Large-Construction Kit (Qiagen, USA), according to the manufacturer's protocol.
Fosmid DNA extracted from all pools for both libraries was then combined to produce a single DNA sample, a master pool of 5 μg DNA, for subsequent pyrosequencing. This master pool was firstly nebulized, for 2 minutes, to produce DNA fragments with average size of 500 bp. The DNA fragment ends were polished and purified to remove fragments < 500 bp and the Roche 454 GS FLX Titanium sequencing (454 Life Sciences Branford, CT, USA) was performed according to 454/Roche GS-FLX protocols.
Assembly and sequence analysis
The sorting and trimming of the metagenomic data, based on quality and size of the reads, as well as the contig assembling were done using the 454 Newbler assembler (version 2.0.01.14) for Genome Sequencer FLX (Roche, USA).
For the phylogenetic classification, the metagenomic data were compared using BLASTn (E-value ≤ 1e-2) against all 16S rRNA sequences deposited in the European Ribosomal RNA Database (http://bioinformatics.psb.ugent.be/webtools/rRNA/). Then, all reads present in the metagenomic data that matched with any rRNA sequence were selected and submitted to a phylogenetic classification using the RDP naïve Bayesian rRNA Classifier tool, at 80% confidence threshold, and E-value of ≤ 1e-2, available on the Ribosomal Database Project's Pyrosequencing Pipeline website.
The contig and singlet sequences were analyzed with MEGAN software and the MG-RAST platform. These sequence sets were submitted to protein BLAST for sequence matching using the NR database (http://blast.ncbi.nlm.nih.gov/Blast.cgi). The resulting BLASTp file was then analyzed using MEGAN metagenomic software (Huson et al. 2007). MEGAN software uses a homology-matching algorithm to generate a phylogenetic tree based on the NCBI taxonomic database.
The metabolic pathway classification of contig and singlet sequences was done using the MG-RAST platform (http://metagenomics.nmpdr.org/metagenomics.cgi?page=Logout). The data were submitted as text file (.txt) for annotation using the subsystems technology. In this approach, reads are classified in a hierarchical structure in which all genes required for a specific function are grouped into subsystems (Aziz et al. 2008). An overview of the experimental approach used in this study is illustrated in Figure Figure11.
Figure 1
Figure 1
Overview of the experimental approach employed in this work.
Nucleotide sequence accession number
The raw data set from the fosmid libraries is available from the NCBI Short Read Archive (SRX038779.2).
GC-MS analysis of phenol of acclimated sample
To verify if the MBR2 sample microorganisms were degrading phenol after acclimation process, the phenol concentration was determined by 1:1 extraction with CH2Cl2 according to the procedures described in Standard Methods (1998). Phenol quantitation was done using a calibration curve constructed from a fortified matrix in three levels. Growth culture samples were filtered through a 0.45 μm membrane filter to remove suspended biomass before analysis. One microliter of extracted sample was subsequently injected into a gas chromatography/mass spectrometry (model 6890 N, Agilent, USA) equipped with a 5975 selective mass detector and 7683B automatic injector. The sample components were separated in the HP-5 MS (30 m × 0.25 mm × 0.25 μm) capillary column using the following temperature programming: 60°C (3°/min), 240°C/10 min. The injector temperature was 220°C and the interface temperature was 250°C. Split ratio was 1/100 and the monitoring fragments were 94 (quantification ion), 66 and 65 (identification ions).
MBR2 sludge acclimation
The acclimation step was done only for the MBR2 sample, since this sludge sample originated from a membrane bioreactor operating without high phenol concentration. The biomass was progressively acclimated to a maximum phenol concentration of 1.0 g L-1 for 30 days, and after this period the biomass was able to degrade 1.0 g L-1 of phenol in less than 24 hours.
The biomass acclimation was checked by monitoring phenol degradation using visual and chromatographic analysis during 30 days. The visual analysis consisted of observation of microbial aggregates (flocculation) and the production of yellow pigments in the culture medium, which, according to Harayama et al. (1992), is indicative of phenol biodegradation through the meta-cleavage pathway. The yellow color comes from the product of cathecol cleavage reaction, 2-hydroxymuconic semialdehyde. Additionally, the phenol degradation was monitored by chromatographic experiments, which demonstrated that 90% of phenol was degraded in the first 12 hours, and in less than 24 hours 1.0 g L-1 of phenol was totally removed by the acclimated biomass.
Metagenomic fosmid libraries and pyrosequencing
Two metagenomic fosmid libraries were made using high molecular weight DNA extracted from the total microbial community of MBR1 and MBR2 sludge samples, yielding 10,000 and 3,200 clones, respectively. The libraries were validated by enzyme digestion in order to check the diversity of clones obtained. Pulsed field gel electrophoresis revealed distinct band profiles for the six randomly selected clones from each library, with the vector band present in the same position for all clones, confirming the occurrence of different insertion events.
The 13,200 clones obtained from both libraries were then submitted to fosmid DNA extraction in pools and subsequent pyrosequencing. A single pyrosequencing run was done for two libraries, yielding 322,742 reads with an average sequence length of 263 bp, the largest sequence being 735 bp. The metagenomic data was assembled using Newbler (version 2.0.01.14), which employs algorithm attempts to combine individual sequence reads into longer contigs. The assembly recruited 108,000 sequences, which resulted in 22,267 contigs and 99,786 singletons.
Phylogenetic composition of phenol degrading sludge by MEGAN
Community taxonomic analysis of metagenomic data was performed using MEGAN (Metagenome analysis software). This software generated a phylogenetic tree based on the NCBI taxonomic database, in which the size of each circular node is proportional to the number of assignments at the particular taxonomic level. The phylogenetic tree showed that the majority of the metagenomic sequences were affiliated to the Bacteria domain (64%) and the remaining sequences were distributed between the Archaea and Eukarya domains (Figure (Figure2).2). The eukaryotes were represented mainly by fungi and protozoa, common organisms present in the sludge of wastewater treatment system (Jiang et al. 2005; Oliveira et al. 2009). A considerable percentage of metagenomic sequences (34%) were unaffiliated. Possible these sequences are distantly related with any known sequence deposited in the public databases.
Figure 2
Figure 2
Phylogenetic tree of all reads derived from the enriched refinery sludge community by MEGAN analysis.
The Bacteria domain was represented by sixteen different phyla encompassing Proteobacteria, Candidate Division WWE1, Bacteroidetes/Chlorobi, Planctomycetes, Firmicutes, Deinococcus-Thermus, Fusobacteria, Chloroflexi, Cyanobacteria, Actinobacteria, Thermotogae, Synergistetes, Acidobacteria, Verrucomicrobia, Lentisphaerae and Gemmatimonadetes (Figure (Figure2).2). Among these, the most abundant were Proteobacteria (66%), followed by Bacteroidetes (7%), Actinobacteria (5.5%), Planctomycetes (4.9%) and Verrucomicrobia (1.3%). This finding supports previous reports of our research group that proteobacteria is a predominant phylum in sludge samples (Silva et al. 2010a; Silva et al. 2010b), suggesting that these organisms play key roles in biodegradation processes that take place in petroleum refinery wastewater treatment systems.
The metagenomic data analysis using MEGAN also yielded a profile of attributes for the bacterial community, in which the majority of the bacterial species were shown to be gram-negative, aerobic, motile, non-pathogenic, aquatic and mesophilic (Figure (Figure3).3). These characteristics are in accordance with those displayed by the vast majority of proteobacteria detected in the sludge of WWTPs (Miura et al. 2007; Kraigher et al. 2008; Sanapareddy et al. 2009, Silva et al. 2010a, Silva et al. 2010b).
Figure 3
Figure 3
Attributes of the metagenomic library of phenol degrading sludge from membrane bioreactor obtained by MEGAN analysis of all reads.
Phylogenetic classification of metagenomic sequences
The comparison of all metagenomic sequences with the European Database allowed the selection of 6,099 sequences presenting some 16S rRNA gene similarity, this dataset was submitted to automatic annotation using the RDP Classifier. This algorithm platform uses Bayesian statistics to assign 16S rRNA genes sequences to known taxa. A set of 247 reads was annotated at least to the phylum level based on known 16S rRNA genes sequences from RDP database at an E-value cutoff ≤ 0.01. The selected size of sort reads in the pyrosequencing data was 100 nt at least. This sequence length was chosen based on the work of Liu et al. (2007) showing that reads as short as 100 bases are long enough to accurately characterize taxa. The annotation results revealed the presence of six different phyla in the sludge samples under study, named Proteobacteria, Actinobacteria, Chloroflexi, Planctomycetes, Verrucomicrobia and Bacteroidetes (Figure (Figure4).4). Despite the fact that this analysis uses only 16S rRNA gene sequences for the classification, most of the results were consistent with the MEGAN analysis (Figure (Figure22).
Figure 4
Figure 4
Classification of the 16S rRNA gene sequences of metagenomic library of phenol degrading sludge from membrane bioreactor using the RDP database. (a) Phylum, (b) Class, (c) Order and (d) Genus. Sequences that were not classified at the Phylum level were (more ...)
At the phylum level, the observed taxa were dominated by representatives of Proteobacteria, which accounted for 91% of the classified sequences (Figure (Figure4a).4a). From these, about 45% were affiliated to class Betaproteobacteria (Figure (Figure4b).4b). The majority of the 16S rRNA sequences were classified at the order level and, distributed across 11 different orders with Burkholderiales (Betaproteobacteria) the most abundant, followed by Rhizobiales (Alphaproteobacteria) and Pseudomonadales (Gammaproteobacteria) (Figure (Figure4c).4c). At the genus level, the sequences could be classified in 17 known genera of which Diaphorobacter, Pseudomonas, Hyphomicrobium, Comamonas, Planctomyces and Thauera were the most prevalent. As classification of the sequences moved from phylum to genus, fewer sequences could be classified with a RDP confidence score of at least 80%, and 67% of the sequences were considered unassigned (Figure (Figure4d),4d), confirming previous reports on the huge bacterial diversity of wastewater treatment sludges (Miura et al. 2007; Sanapareddy et al. 2009; Silva et al. 2010a; Silva et al. 2010b).
Functional analysis of the metagenome phenol degrading sludge
The sequences were functionally annotated based upon MG-RAST plataform using the SEED database. The functional analysis allowed the classification of 42.8% of the metagenomic sequences into several subsystems with the majority of them related to the metabolism of carbohydrates, aromatic compounds, proteins, DNA, amino acids and derivatives, virulence, respiration and cofactors. By contrast, some subsystems were under-represented in the metagenomic data, including sequences involved in photosynthesis, secondary metabolism, macromolecular synthesis and nitrogen metabolism (Figure (Figure5).5). Despite of the low number of reads assigned to the nitrogen metabolism, the majority of these sequences represented genes coding enzymes linked to processes such as nitrate and nitrite ammonification and denitrification, which are considered important processes for the removal of complex nitrogen compounds in wastewater treatment. An interesting result was observed for the virulence subsystem (Figure (Figure5),5), which accounted for a high number of reads. Annotation of such reads revealed that almost 50% of them were related to antibiotic and toxic compounds resistance genes, which are extremely important features for microbial survival and adaptation in polluted environments.
Figure 5
Figure 5
Metabolic profile of phenol degrading sludge metagenomic libraries datasets using MG-RAST platform.
Analysis of metagenome for phenol degradation
The phenol can be aerobically degraded for two different pathways, ortho- or meta-pathway. The aromatic ring is initially monohydroxilated in the adjacent carbon of a hydroxyl group by the enzyme phenol hydroxylase (EC 1.14.13.7) resulting in catechol, which is in turn cleaved by either ortho- or meta-cleavage pathway. In case of the ortho-pathway, the ring is cleaved by the catechol 1,2-dioxygenase enzyme (EC 1.13.11.1), leading to the initial formation of succinyl-CoA and acetyl-CoA. In the meta- pathway, the catechol is cleaved by the catechol 2,3-dioxygenase enzyme (C23O), leading to the formation of pyruvate and acetaldehyde (Marimaa et al. 2006).
In MG-RAST analysis, the aromatic compound utilization profile of phenol degrading sludge was dominated by proteins annotated as subsystems of peripheral pathway for catabolism of aromatic compounds and metabolism of central aromatic intermediates, besides benzoate degradation, cresol degradation and toluene 4-mooxygenase, that are also important pathways involved in the catabolism of pollutant compounds found in refinery wastewater (Figure (Figure66).
Figure 6
Figure 6
Percentage of sequences associated to aromatic compound metabolism of phenol degrading sludge metagenomic libraries from membrane bioreactor by MG-RAST platform. Total of aromatic compounds metabolism reads = 3800.
The subsystem of peripheral pathway for catabolism of aromatic compounds encompassed the proteins involved with the degradation of phenol and derivatives, such as phenol hydroxylase (EC 1.14.13.7) (78 reads) and enzymes involved in the biphenyl degradation (75 reads). The subsystem of metabolism of central aromatic intermediates, although mainly represented by the homogentisate pathway, involved in the catabolism of aromatic rings, contained the proteins involved in the ortho- and meta-pathway of phenol degradation, such as catechol 1,2-dioxygenase (EC 1.13.11.1) (9 reads) and catechol 2,3- dioxygenase (EC 1.13.11.2) (23 reads), respectively.
These analyses showed that the genes related to phenol degradation have been well sampled in the phenol degrading sludge, revealing that the phenol enrichment was successful. Further studies will be carried out to verify if these genes represent new sequences coming from uncultivated microorganisms.
Analysis of the metabolic potential for organic matter removal
This analysis was done to investigate the presence of some important functions related to the biological treatment performance in the microbial sludge under study after the phenol enrichment.
The potential metabolic functions for organic matter removal were established by comparison of all pyrosequencing-derived sequences from phenol degrading sludge to the KEGG maps using functional assignments from SEED analysis (Meyer et al. 2008). Within the KEGG categories, matches were separated into different subcategories and several sequences were assigned to more than one subcategory. The most represented categories were metabolism of energy, amino acids, carbohydrates, cofactors and nucleotides. These categories showed a high number of distinct Enzyme Commissions (EC) belonging to a specific metabolic pathway. The EC number is the common name for IUBMB's Enzyme Nomenclature System, and plays a key role in classifying enzymatic reactions and linking the enzyme or protein genes to reactions in metabolic pathways (Yamanishi et al. 2009).
Among the KEGG categories found, the energy metabolism and the biodegradation of xenobiotics are emphasized, since they constitute important processes for wastewater biological treatment, such as removal of organic matter, including pollutant compounds, from the petroleum refining process. About 4,000 metagenomic sequences were related to the energy metabolism category, which involved homology to known genes of carbon fixation, oxidative phosphorylation, nitrogen metabolism and CO2 fixation (Table (Table1).1). On the other hand, the number of genes that matched with methane metabolism subcategories was low, which is consistent with the aerobic biological treatment. In the sulfur metabolism subcategory, most sequences showed similarities with enzymes related to sulfur assimilation pathways, such as cystathionine synthase, homoserine acetyltransferase, homoserine succinyltransferase and sulfate kinases. In the nitrogen metabolism subcategory, the majority of the sequences were related with nitrate and nitrite reductase, enzymes involved in nitrogen assimilation pathways. A total of 2,495 metagenomic sequences were assigned to the xenobiotic biodegradation category, which showed, in general, a more even distribution among the subcategories when compared to the previous category. Exceptions were the benzoate degradation via CoA ligation subcategory, which presented the highest number of matches, and the carbazole degradation subcategory, which in turn showed the lowest number of matches (Table (Table11).
Table 1
Table 1
Number of sequences showing homology to genes associated with KEGG pathways in the categories "carbohydrate metabolism", "biodegradation of xenobiotics" and "energy metabolism"
In addition, the analysis of the carbohydrate metabolism category showed a high number of matches (9,387 sequences) related to genes involved in several subcategories for utilization of carbohydrates (Table (Table1).1). These results demonstrated the broad capability of microorganisms from refinery sludge to metabolize distinct sugar sources.
This is the first report of a pyrosequencing approach for a broad phylogenetic and metabolic diversity analysis of metagenomic fosmid libraries derived from phenol degrading sludge samples of petroleum refinery WWTP. Sequencing of the metagenomic libraries has allowed us to have deeper insight to the complex metagenome of the phenol degrading sludge samples, with thousands of reads assigned to different taxa and metabolic categories driving the functioning of the microbial community in the membrane bioreactors.
Based on previous studies of the use of short sequence reads (~90 bp) to accurately classify microbial communities (Liu et al. 2007; Sanapareddy et al. 2009), this work used MEGAN and RDP Classifier tools for obtaining a microbial phylogenetic profile of the metagenomic data. The results obtained using both tools were consistent and showed that the most abundant phyla in the metagenomic data were similar, despite the fact that the MEGAN analysis is based on all metagenomic reads, while the RDP Classifier is based on only the 16S rRNA gene reads.
Results of RDP classifier were compared to a previous report by Silva and collaborators (2010a), where the authors used the phenol degrading sludge sample (MBR1) to construct a 16S rRNA gene library. The comparison of the two libraries showed that Proteobacteria was the predominant phylum, but the richness of phyla was higher in the metagenomic dataset, which contained 6 different phyla, whereas the 16S rRNA library revealed only 3 phyla for the MBR 1 sample. However, the comparison between deeper levels, such as order and genus, showed that the richness between both libraries is similar, although the groups found were different. Probably, these differences can be explained by the bias inherent to each method, e.g. pyrosequencing and 16S rRNA gene library, used to survey the bacterial diversity. Nonetheless, in this case, the data obtained using both approaches can be considered complementary in order to depict a bacterial diversity scenario of the phenol degrading sludge.
The high abundance of the Proteobacteria group in microbial communities from wastewater treatment samples was also observed in several other studies using PCR-based experiments targeting the 16S rRNA gene (Miura et al. 2007; Ahmed et al. 2007; Li et al. 2009; Sanapareddy et al. 2009; Silva et al. 2010a; Silva et al. 2010b). The predominance of Proteobacteria in such environments could be explained by the fact that this phylum comprises the most phylogenetically diverse group in the Bacteria Domain, known to be metabolically versatile, including aerobic and facultative aerobic bacteria (Madigan et al. 2008). These are quite interesting characteristics for microorganisms inhabiting wastewater treatment plant, an environment showing great daily variations in the composition and concentration of pollutant compounds. These results demonstrate that phylogenetic studies based on pyrosequencing of metagenomic fosmid libraries can give broad and reliable information about the predominant microbial groups present in the microbial community sampled.
Within the Proteobacteria, a few genera containing important species, e.g. Pseudomonas putida (Marques and Ramos 1993, Gonzales et al. 2001), Thauera aromatica (Breinig et al. 2000), Thauera aminoaromatica, Thauera phenylacetica (Mechichi et al. 2002), Thauera sp. DNT-1 (Shinoda et al. 2004) and Comamonas testosterone, are able to utilize different kinds of aromatic compounds, including phenol, polyphenol, toluene and halobenzoate as carbon sources. Functional studies based on SIP-RNA have revealed that members of Thauera genus dominated the phenol degradation process in bioreactor sludges (Manefield et al. 2002; Valle et al. 2004). Additionally, studies based on culturing analyses have showed that Comamonas testosteroni can be involved with the metabolism of aromatic compounds, such as phenol and 4-clorophenol (Bae et al. 1996; Arai et al. 1998). Basu et al. (2006) verified that Pseudomonas putida CSV86 is able to degrade preferentially naphthalene over glucose. Agarry and Solomon (2008) studied batch culture using synthetic phenol and observed that Pseudomonas fluorescens was able to degrade phenol in the concentration range of 100-500 mg/L.
A considerable fraction of the metagenomic sequence data was not assigned at the genus level, demonstrating the astonishing microbial diversity present in the sludge from wastewater treatment plants. Similar findings were reported by Sanapareddy et al. (2009) when analyzing sewage sludge from biologic treatment system, and these data corroborate other previous studies of complex environments, such as soil and oceans, in which the sequence classification becomes worse as one moves into deeper taxonomic levels (Roesch et al. 2007; Brown et al. 2009).
The phylogenetic richness observed reflected the wide metabolic diversity present in the metagenomic data from the WWTP sludge. As expected, genes assigned to the metabolism of carbohydrates, amino acids and proteins were more numerous since they are related to the housekeeping functions of all living organisms. The presence of sequences assigned to functions like nitrogen, phosphor, sulfur and aromatic compounds metabolisms are essential for the high performance of wastewater treatment plant, since they are indicative that the microorganisms from the sludge are degrading and/or assimilating such compounds. High concentrations of these organic compounds, such as ammonia, sulfate, phosphate, phenol may be toxic to human health and cause negative environmental impacts (Nair et al. 2008).
The metabolic profile revealed a broad set of important genes related to the utilization and mineralization of aromatic and/or xenobiotics compounds, including some key enzymes related to aerobic phenol degradation, such as phenol hydroxylase, catechol 1,2 dioxygenase and catechol 2,3 dioxygenase. A high number of sequences related to benzoate degradation pathway was observed as well, what is coherent with the fact that benzoate is a central intermediary compound in the anaerobic and aerobic metabolism of various aromatic compounds, such as toluene, xylene, fluorene, carbazole and biphenyl (Kim and Harwood 1991). In addition, the anaerobic phenol degradation via carboxylation of phenol to 4-hydroxybenzoate ends in the anaerobic benzoate pathway (Lack and Fuchs 1994). These results suggest that the microorganisms of the sludge under study have potential to degrade phenol and derivatives via aerobic and anaerobic pathway. This versatility is actually a great benefit to the performance of the biological treatment.
Finally, the phylogenetic and metabolic diversities observed are an indicative that the phenol enrichment of the sludge did not affect other important functions besides phenol degradation, which are necessary for the efficient performance of biological treatment systems. Additionally, the metagenome dataset generated by pirosequencing may provide useful sequence information for the characterization of whole catabolic pathways, particularly phenol degrading pathway, that support fundamental key processes occurring in the wastewater treatment plants. Future studies will be conducted aiming at the design and use of probes or primers to detect fosmid clones bearing specific target new genes and pathways related to pollutant compound degradation, thus offering efficient tools for the improvement of bioremediation technologies.
Competing interests
The authors declare that they have no competing interests.
Acknowledgements
The authors are grateful to Petrobras for technical and financial support. C.C.S. was supported by grants from FAPESP (Fundação de Amparo à Pesquisa do Estado de São Paulo). RK acknowledges FAPDF and CNPq.
  • Agarry SE, Solomon BO. Kinetics of batch microbial degradation of phenols by indigenous Pseudomonas fluorescenc. Int J Environ Sci Tech. 2008;5:223–232.
  • Ahmed Z, Cho J, Lim B-R, Song K-G, Ahn K-H. Effects of sludge retention time on membrane fouling and microbial community structure in a membrane bioreactor. J Membrane Sci. 2007;287:211–218. doi: 10.1016/j.memsci.2006.10.036. [Cross Ref]
  • Amann RI, Ludwig W, Schleifer K. Phylogenetic identification and in sit detection of individual microbial cells without cultivation. Microbiol Rev. 1995;59:143–169. [PMC free article] [PubMed]
  • Arai H, Akahira S, Ohishi T, Maeda M, Kudo T. Adaptation of Comamonas testosteron TA441 to utilize phenol: organization and regulation of the genes involved in phenol degradation. Microbiology. 1998;144:2895–2903. doi: 10.1099/00221287-144-10-2895. [PubMed] [Cross Ref]
  • Aziz AK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, Formsma K, Gerdes S, Glass EM, Kubal M, Meyer F, Olsen GJ, Olson R, Osterman AL, Overbeek RA, McNeil LK, Paarmann D, Paczian T, Parrello B, Pusch GD, Reich C, Stevens R, Vassieva O, Vonstein V, Wilke A, Zagnitko O. The RAST server: Rapid annotation using subsystems technology. BMC Genomics. 2008;9:75. doi: 10.1186/1471-2164-9-75. [PMC free article] [PubMed] [Cross Ref]
  • Bae HS, Lee JM, Kim YB, Lee ST. Biodegradation of the mixtures of 4-chlorophenol and phenol by Comamonas testosteron CPW301. Biodegradation. 1996;7:463–469. [PubMed]
  • Basu A, Apte AK, Phale PS. Preferential utilization of aromatic compounds over glucose by Pseudomonas putid CV86. Appl Environ Microbiol. 2006;72:2226–2230. doi: 10.1128/AEM.72.3.2226-2230.2006. [PMC free article] [PubMed] [Cross Ref]
  • Barrios-Martinez A, Barbot E, Marrot B, Moulin P, Roche N. Degradation of synthetic phenol-containing wastewaters by MBR. Journal Membrane Science. 2006;281:288–296. doi: 10.1016/j.memsci.2006.03.048. [Cross Ref]
  • Billy K, Viau E, Peccia J. Pyrosequencing of the 16S rRNA gene to reveal bacterial pathogen diversity in biosolids. Water Res. 2010;44:4252–4260. doi: 10.1016/j.watres.2010.05.039. [PubMed] [Cross Ref]
  • Braile PM. Manual of Industrial Wastewater Treatment. 1. CETESB, São Paulo; 1979.
  • Breinig S, Schiltz E, Fuchs G. Genes Involved in Anaerobic Metabolism of Phenol in the Bacterium Thauera aromatic. J Bacteriol. 2000;182:5849–5863. doi: 10.1128/JB.182.20.5849-5863.2000. [PMC free article] [PubMed] [Cross Ref]
  • Brown MV, Philip GK, Bunge JA, Smith MC, Bissett A, Lauro FM, Fuhrman JA, Donachie SP. Microbial community structure in the North Pacific Ocean. ISME J. 2009;3:1374–1386. doi: 10.1038/ismej.2009.86. [PubMed] [Cross Ref]
  • Chang BV, Chiang F, Yuan SY. Biodegradation of nonylphenol in sewage sludge. Chemosphere. 2005;60:1652–1659. doi: 10.1016/j.chemosphere.2005.02.042. [PubMed] [Cross Ref]
  • Gonzalez G, Herrera G, García MT, Peña M. Biodegradation of phenolic industrial wastewater in a fluidized bed bioreactor with immobilized cells of Pseudomonas putid. Biores Technol. 2001;80:17–142. [PubMed]
  • Handelsman J. Metagenomics: Application of Genomics to Uncultured microorganisms. Microbiol Mol Biol Rev. 2004;68:669–685. doi: 10.1128/MMBR.68.4.669-685.2004. [PMC free article] [PubMed] [Cross Ref]
  • Harayama S, Kok M, Neidle EL. Functional and evolutionary relationships among diverse oxygenases. Annu Rev Microbiol. 1992;46:565–601. doi: 10.1146/annurev.mi.46.100192.003025. [PubMed] [Cross Ref]
  • Henze M, Harremoes P, la Cour JJ, Arvin E. In: Wastewater treatment. 2. Forstner U, Murphy RJ, Rulkens WH, editor. Berlin: Springer; 1997. Environmental Engineering.
  • Huson DH, Auch AF, Qi J, Schuster SC. MEGAN analysis of metagenomic data. Genome Res. 2007;17:377–386. doi: 10.1101/gr.5969107. [PubMed] [Cross Ref]
  • Jiang Y, Wen JP, Li HM, Yang SL, Hu ZD. The biodegradation of phenol at high initial concentration by the yeast Candida tropicali. Biochem Eng J. 2005;24:243–247. doi: 10.1016/j.bej.2005.02.016. [Cross Ref]
  • Kim M-K, Harwood CS. Regulation of benzoate-CoA ligase in Rhodopseudomonas palustri. FEMS Microbiol Let. 1991;83:199–203.
  • Kraigher B, Kosjek T, Heath E, Kompare B, Mandic-Muleca I. Influence of pharmaceutical residues on the structure of activated sludge bacterial communities in wastewater treatment bioreactors. Water Res. 2008;42:4578–4588. doi: 10.1016/j.watres.2008.08.006. [PubMed] [Cross Ref]
  • Lack A, Fuchs G. Evidence that phenol phosphorylation to phenylphosphate is the first step in anaerobic phenol metabolism in a denitrifying Pseudomona sp. Arch Microbiol. 1994;161:132–139. [PubMed]
  • Li X, Ma H, Wang Q, Matsumoto S, Maeda T, Ogawa HI. Isolation, identification of sludge-lysing strain and its utilization in thermophilic aerobic digestion for waste activated sludge. Bioresource Technology. 2009;100:2475–2481. doi: 10.1016/j.biortech.2008.12.019. [PubMed] [Cross Ref]
  • Liu Z, Lozupone C, Hamady M, Bushman FD, Knight R. Short pyrosequencing reads suffice for accurate microbial community analysis. Nucleic Acids Res. 2007;35:e120. doi: 10.1093/nar/gkm541. [PMC free article] [PubMed] [Cross Ref]
  • Madigan M, Martinko J, Dunlap P, Clark D. Brock Biology of Microorganisms. 12. New York, USA: Pearson Higher Education; 2008.
  • Manefield M, Whiteley AS, Griffiths GI, Bailey MJ. RNA stable isotope probing, a novel means of linking microbial community function to phylogeny. Applied Environment and Microbiology. 2002;68:5367–5373. doi: 10.1128/AEM.68.11.5367-5373.2002. [PMC free article] [PubMed] [Cross Ref]
  • Mariano JB. Environmental impacts of petroleum refining. COPPE/UFRJ, Rio de Janeiro, RJ, Brazil: MSc These; 2001.
  • Merimaa M, Heinaru E, Liivak M, Vedler E, Heinaru A. Grouping of phenol hydroxylase and catechol 2,3-dioxygenase genes among phenol- and p-cresol-degrading Pseudomonas species and biotypes. Arch Microbiol. 2006;186:287–296. doi: 10.1007/s00203-006-0143-3. [PubMed] [Cross Ref]
  • Marques S, Ramos JL. Transcriptional control of the Pseudomonas putid TOL plasmid catabolic pathways. Mol Microbiol. 1993;9:923–929. doi: 10.1111/j.1365-2958.1993.tb01222.x. [PubMed] [Cross Ref]
  • Mechichi T, Stackebrandt E, Gadon N, Fuchs G. Phylogenetic and metabolic diversity of bacteria degrading aromatic compounds under denitrifying conditions, and description of Thauera phenylacetic sp. nov., Thauera aminoaromatic sp. nov., and Azoarcus buckeli sp. nov. Arch Microbiol. 2002;178:26–35. doi: 10.1007/s00203-002-0422-6. [PubMed] [Cross Ref]
  • Meyer F, Paarmann D, D'Souza M, Olson R, Glass EM, Kubal M, Paczian T, Rodriguez A, Stevens R, Wilke A, Wilkening J, Edwards RA. The Metagenomics RAST server - A public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinformatics. 2008;9:386. doi: 10.1186/1471-2105-9-386. [PMC free article] [PubMed] [Cross Ref]
  • Miura Y, Hiraiwa MN, Ito T, Itonaga T, Watanabe Y, Okabe S. Bacterial community structures in MBRs treating municipal wastewater: Relationship between community stability and reactor performance. Water Res. 2007;41:627–637. doi: 10.1016/j.watres.2006.11.005. [PubMed] [Cross Ref]
  • Nair CI, Jayachandran K, Shashidhar S. Biodegradation of phenol. African J Biotechnol. 2008;7:4951–4958.
  • Ojumu TV, Bello OO, Sonibare JA, Solomon BO. Evaluation of microbial systems for bioremediation of petroleum refinery effluents in Nigeria. African J Biotechnol. 2005;4:31–35.
  • Oliveira GSS, Araújo CVM, Fernandes JGS. Activated sludge system microbiology and its relationship with the industrial effluents treatment: the experience of the Cetrel. Eng Sanit Ambient. 2009;14:183–192.
  • Roesch LF, Fulthorpe RR, Riva A, Casella G, Hadwin AKM, Kent AD, Daroub SH, Camargo FAO, Farmerie WG, Triplett EW. Pyrosequencing enumerates and contrasts soil microbial diversity. ISME J. 2007;1:283–290. [PMC free article] [PubMed]
  • Sanapareddy N, Hamp TJ, Gonzalez LC, Hilger HA, Fodor AA, Clinton SM. Molecular diversity of a north carolina wastewater treatment plant as revealed by pyrosequencing. Appl Environ Microbiol. 2009;175:1688–1696. [PMC free article] [PubMed]
  • Schlüter A, Krause L, Szczepanowski R, Goesmann A, Pühler A. Genetic diversity and composition of a plasmid metagenome from a wastewater treatment plant. J Biotechnol. 2008;136:65–72. doi: 10.1016/j.jbiotec.2008.03.017. [PubMed] [Cross Ref]
  • Sercu B, Boon N, Verstraete W, van Langenhove H. H2S degradation is reflected by both the activity and composition of the microbial community in a compost biofilter. Appl Microbiol Biotechnol. 2006;72:1090–1098. doi: 10.1007/s00253-006-0382-x. [PubMed] [Cross Ref]
  • Shinoda Y, Sakai Y, Uenishi H, Uchihashi Y, Hiraishi A, Yukawa H, YurimotoH KN. Aerobic and anaerobic toluene degradation by a newly isolated denitrifying bacterium, Thauer sp. strain DNT-1. Appl Environ Microbiol. 2004;70:1385–1392. doi: 10.1128/AEM.70.3.1385-1392.2004. [PMC free article] [PubMed] [Cross Ref]
  • Silva CC, Jesus EC, Torres APR, Sousa MP, Santiago VMJ, Oliveira VM. Investigation of bacterial diversity in membrane bioreactor and conventional activated sludge processes from petroleum refineries using phylogenetic and statistical approaches. J Microbiol Biotechnol. 2010;20:447–459. [PubMed]
  • Silva CC, Viero AF, Dias ACF, Andreote FD, Jesus EC, De Paula SO, Torres APR, Santiago VMJ, Oliveira VM. Monitoring the bacterial community dynamics in a petroleum refinery wastewater membrane bioreactor fed with a high phenolic load. J Microbiol Biotechnol. 2010;20:17–25. [PubMed]
  • Sleator RD, Shortall C, Hill C. Under the microscope: Metagenomics. Let Appl Microbiol. 2008;47:361–366. doi: 10.1111/j.1472-765X.2008.02444.x. [PubMed] [Cross Ref]
  • Standard Methods for the Examination of Water and Wastewater. American Public Health Association and American Water Works Association and Water Pollution Control Federation. 20. USA: Washington; 1998.
  • Steele HL, Jaeger K-E, Daniel R, Streit WR. Advances in recovery of novel biocatalysts from metagenomes. J Mol Microbiol Biotechnol. 2009;16:25–37. doi: 10.1159/000142892. [PubMed] [Cross Ref]
  • Stepnowski P, Siedlecka EM, Behrend P, Jastorff B. Enhanced photo-degradation of contaminants in petroleum refinery wastewater. Water Res. 2002;36:2167–2172. doi: 10.1016/S0043-1354(01)00450-X. [PubMed] [Cross Ref]
  • Szczepanowski R, Bekel T, Goesmann A, Krause L, Krömeke H, Kaiser O, Eichler W, Pühler A, Schlüter A. Insight into the plasmid metagenome of wastewater treatment plant bacteria showing reduced susceptibility to antimicrobial drugs analyzed by the 454-pyrosequencing technology. J Biotechnol. 2008;136:54–64. doi: 10.1016/j.jbiotec.2008.03.020. [PubMed] [Cross Ref]
  • Torsvik V, Goksoyr J, Daae FL. High diversity in DNA of soil bacteria. Appl Environ Microbiol. 1990;56:782–787. [PMC free article] [PubMed]
  • Valle A, Bailey MJ, Whiteley AS, Manefield M. N-acyl- L -homoserine lactones (AHLs) affect microbial community composition and function in activated sludge. Environmental Microbiology. 2004;6:424–433. doi: 10.1111/j.1462-2920.2004.00581.x. [PubMed] [Cross Ref]
  • Viero AF, Melo TM, Torres APR, Ferreira NR, Sant'Anna GL Jr, Borges CP, Santiago VMJ. The effects of long-term feeding of high organic loading in a submerged membrane bioreactor treating oil refinery wastewater. J Membrane Sci. 2008;319:223–230. doi: 10.1016/j.memsci.2008.03.038. [Cross Ref]
  • Yamanishi Y, Hattori M, Kotera M, Goto S, Kanehisa M. E-zyme: predicting potential EC numbers from the chemical transformation pattern of substrate-product pairs. Bioinformatics. 2009;25:179–182. doi: 10.1093/bioinformatics/btp223. [PMC free article] [PubMed] [Cross Ref]
  • Zang K, Kurisu F, Kasuga I, Furumai H, Yagi O. Analysis of the phylogenetic diversity of estrone-degrading bacteria in activated sewage sludge using microautoradiography-fluorescence in situ hybridization. Systematic Appl Microbiol. 2008;31:206–214. doi: 10.1016/j.syapm.2008.03.005. [PubMed] [Cross Ref]
  • Zhang H, Banaszak JE, Parameswaran P, Alder J, Krajmalnik-Brown R, Rittmann BE. Focused-pulsed sludge pre-treatment increases the bacterial diversity and relative abundance of acetoclastic methanogens in a full-scale anaerobic digester. Water Res. 2009;43:4517–4526. doi: 10.1016/j.watres.2009.07.034. [PubMed] [Cross Ref]
  • Zhou J, Xia B, Treves DS, Wu L-Y, Marsh TL, O'Neill RV, Palumbo AV, Tiedje JM. Spatial and resource factors influencing high microbial diversity in soil. Appl Environ Microbiol. 2001;68:326–334. [PMC free article] [PubMed]
Articles from AMB Express are provided here courtesy of
Springer-Verlag