Original Research ARTICLE
Molecular and Microbiological Insights on the Enrichment Procedures for the Isolation of Petroleum Degrading Bacteria and Fungi
- 1Department for Sustainable Food Processes, Università Cattolica del Sacro Cuore, Piacenza, Italy
- 2Department of Life Sciences and Systems Biology, Mycotheca Universitatis Taurinensis, University of Turin, Turin, Italy
- 3Satt Grand-EST, Maison Règionale de l’Innovation, Dijon, France
- 4ARPAE-Emilia Romagna, Bologna, Italy
Autochthonous bioaugmentation, by exploiting the indigenous microorganisms of the contaminated environment to be treated, can represent a successful bioremediation strategy. In this perspective, we have assessed by molecular methods the evolution of bacterial and fungal communities during the selective enrichment on different pollutants of a soil strongly polluted by mixtures of aliphatic and polycyclic hydrocarbons. Three consecutive enrichments were carried out on soil samples from different soil depths (0–1, 1–2, 2–3 m), and analyzed at each step by means of high-throughput sequencing of bacterial and fungal amplicons biomarkers. At the end of the enrichments, bacterial and fungal contaminants degrading strains were isolated and identified in order to (i) compare the composition of enriched communities by culture-dependent and culture-independent molecular methods and to (ii) obtain a collection of hydrocarbon degrading microorganisms potentially exploitable for soil bioremediation. Molecular results highlighted that for both bacteria and fungi the pollutant had a partial shaping effect on the enriched communities, with paraffin creating distinct enriched bacterial community from oil, and polycyclic aromatic hydrocarbons generally overlapping; interestingly neither the soil depth or the enrichment step had significant effects on the composition of the final enriched communities. Molecular analyses well-agreed with culture-dependent analyses in terms of most abundant microbial genera. A total of 95 bacterial and 94 fungal strains were isolated after selective enrichment procedure on different pollutants. On the whole, isolated bacteria where manly ascribed to Pseudomonas genus followed by Sphingobacterium, Bacillus, Stenothrophomonas, Achromobacter, and Serratia. As for fungi, Fusarium was the most abundant genus followed by Trichoderma and Aspergillus. The species comprising more isolates, such as Pseudomonas putida, Achromobacter xylosoxidans and Ochromobactrum anthropi for bacteria, Fusarium oxysporum and Fusarium solani for fungi, were also the dominant OTUs assessed in Illumina.
Oil hydrocarbons are the most widespread environmental pollutants, including n-alkanes, cycloalkanes, and polycyclic aromatic hydrocarbons (PAHs) that have been regarded as serious ecological and public health concerns (Bao et al., 2012). Oil contamination of ecosystems is a serious issue associated with crude oil drilling, transportation, refining, and related activities which demands immediate attention for restoration. Crude petroleum oil is a complex mixture of hydrocarbons mainly composed of saturated and aromatic hydrocarbons, asphaltenes and resins (Moubasher et al., 2015). Aliphatic hydrocarbons consist of readily biodegradable n-alkanes followed by less biodegradable branched and cyclic alkanes (Das and Chandran, 2011). Similarly, PAHs are compounds with two or more aromatic rings that are the most recalcitrant components present at high percentages in crude oil (Haritash and Kaushik, 2009).
Due to the public awareness and to the strict legal constraints on the release of pollutants into the environment, it is necessary to find effective and affordable technologies for the treatment of oil industrial wastes. Bioremediation is a biological approach that relies on the metabolic potential of microorganisms to remove contaminants (Maiti et al., 2008; Megharaj et al., 2011; Hara et al., 2013). The use of bioremediation techniques proved economical, environmentally friendly and flexible (Obi et al., 2016). In this prospect, bioremediation is gaining more and more importance as constructive approach for the remediation of polluted sites. Several studies reported the catabolic abilities of indigenous microorganisms such as fungi, bacteria and algae to degrade hydrocarbons (Dean-Ross et al., 2002; Bundy et al., 2004; Maiti et al., 2008; Wang et al., 2011; Badr El-Din et al., 2014). These microorganisms, adapted to the contaminated environments, are equipped with specific enzyme systems that enable them to use hydrocarbons as sole carbon source. Different hydrocarbon degrading microorganisms such as bacteria and archaea have been found in hydrocarbons contaminated environments (Alonso-Gutierrez et al., 2009; An et al., 2013; Hazen et al., 2013; Head et al., 2014; Fowler et al., 2016). The structure of the microbial community in a soil, influences deeply the degree of oil hydrocarbons degradation. Liu et al. (2011) observed that, at the early stage of remediation, the bacterial community was responsible for the degradation of the saturated and partially aromatic hydrocarbons; the fungal community instead became dominant in decomposing the polar hydrocarbons fraction in post-remediation. Generally, thanks to the variety of extracellular enzymes and fungal hyphae, fungi are the first key players in degrading available contaminants and recalcitrant polymers (Wick et al., 2007; Fernandez-Luqueno et al., 2010; Deshmukh et al., 2016). Fungal mobilization and degradation of contaminants contribute to release bioavailable intermediates on which, in a later stage, the bacterial community can act more easily (Smits et al., 2005; Scullion, 2006; Leonardi et al., 2008; Yuan et al., 2018). Therefore, in order to characterize and to monitor the native microbial community, the dynamics and the functional potential of both bacteria and fungi in polluted ecosystems is essential for the development of a successful bioremediation strategy (Atlas, 1981; Head et al., 2014; Covino et al., 2016).
The success of on-site bioremediation, employing native microbial populations, could be voided by imbalanced nutrients and/or adverse factors (temperature, moisture content, pH, availability of electron donor and/or acceptor, high pollutant concentration, etc.), which are common in contaminated sites (Smith et al., 2015). Bioremediation can be accomplished by either boosting the growth of the indigenous microbial community through biostimulation or by introducing naturally occurring microorganisms with excellent catabolic abilities (bioagumentation) that are adapted to the ecological conditions of the site (Agnello et al., 2016).
Many investigators have reported successful bioremediation, where biostimulation with the addition of appropriate nutrients (N and/or P), to avoid metabolic limitations, resulted in an improved metabolic activity of indigenous microorganisms (Yu et al., 2005; Ghaly and Yusran, 2013; Suja et al., 2014; Smith et al., 2015). On the other hand, autochthonous bioagumentation, based on the re-inoculation in polluted sites of indigenous microorganisms previously enriched under laboratory conditions, enhanced the microbial activities, thus improving the degradation of hydrocarbons (Dueholm et al., 2015). In order to provide an inoculum for bioaugmentation, the isolation of microorganisms in pure culture from these contaminated environments is fundamental. However, due to the complexity of crude oil, a microbial consortium composed of microorganisms endowed with diverse metabolic capacities and syntrophic interactions would work better than a pure culture. Several reports demonstrated the better metabolic versatility of mixed cultures in using hydrocarbon pollutants as sole carbon source in comparison to pure cultures (Cerqueira et al., 2011; Das and Chandran, 2011). In laboratory conditions, bacterial, and fungal co-culture(s) showed improved degradation rates of diesel oil and of polyacyclic romatic hydrocarbons (PAHs) (Wang et al., 2011). Hence, catabolic interactions among different microbial groups during biodegradation is extremely important (Atlas, 1981; Varjani et al., 2015). Although identifications and characterizations of the microorganisms involved in the degrading processes are available (Desai and Banat, 1997; Rojo, 2009), less is known on the biodiversity and dynamics of the native hydrocarbons-degrading microbial community of a contaminated soil, especially during the enrichment process applied to isolate the most effective strains (Omrani et al., 2018). The development of effective bioremediation strategies requires an extensive understanding of the resident microorganisms of these habitats. Recent techniques such as high-throughput sequencing (HTS) have greatly facilitated the advancement of microbial ecological studies in oil-polluted sites.
In the present work we have assessed by molecular methods the evolution of bacterial and fungal communities during the enrichment on different pollutants of a soil strongly polluted by mixtures of aliphatic and polycyclic hydrocarbons. Three consecutive enrichments were carried out on soil samples from different soil depths (0–1, 1–2, 2–3 m), and analyzed at each step by means of HTS of bacterial and fungal amplicons biomarkers. At the end of the enrichment, bacterial, and fungal contaminants degrading strains were cultivated and identified.
Main aims of the work were: (i) to assess the effect of different pollutants used as sole carbon during a selective enrichment procedure source on the diversity of the microbial community of crude oil contaminated soil (ii) to compare the composition of enriched communities by culture-dependent and culture-independent molecular methods and (iii) to obtain a collection of hydrocarbon degrading fungi and bacteria potentially exploitable for soil bioremediation.
Materials and Methods
Sampling and Soil Characteristic
Crude oil contaminated soil samples from an area of the industrial SIN (Site of National Interest) located in Fidenza (Emilia-Romagna, Italy) were collected at three different depths (0–1, 1–2, and 2–3 m). The polluted site (about 80,000 m2 wide) has a long history of industrial exploitation and it is contaminated by a variety of pollutants such as Benzene-Toluene-Ethylbenzene-Xylene (BTEX), n-alkanes and PAHs.
The chemical composition of the soil samples is described in Supplementary Table S1.
Target Organic Contaminants
Target organic contaminants were chosen based on literature data and on the chemical characterization of the contaminated area of Fidenza SIN. All chemicals were purchased by Sigma-Aldrich (Germany). Benzene was selected as representative of BTEX. Pyrene, phenanthrene, and naphthalene were selected as representatives of 4, 3, and 2 rings PAHs. Stock solutions were prepared in methanol for naphthalene (20 mg/mL) and in ethanol 95% for phenanthrene (15 mg/mL) and for pyrene (5 mg/mL). Paraffin oil, representative of alkanes, and crude oil mixture from the contaminated Fidenza site were used at a final concentration of 1% v/v.
Microcosm Enrichments and Microbial Isolation
Different enrichment cultures were set up using the polluted soil as inoculum and the above-mentioned crude oil hydrocarbons components as the sole carbon source. One hundred g of each soil sample at the three depths (namely, S1, S2, and S3) were added to 900 mL sterile Mineral Medium (MM) supplemented with the target analyte as sole carbon source. To provide the proper microelement amount for bacterial and fungal growth, two mineral media were used: M9 mineral medium (Difco, Sparks, MD, United States) for bacteria and Czapek for fungi.
The flasks were incubated on a rotary shaker at 30°C and 180 rpm for bacteria and at 24°C and 120 rpm for fungi. After 7 days, 5 mL of the culture was transferred to another flask with MM and the corresponding target analyte: benzene (50 ppm); naphthalene (200 ppm), phenanthrene (200 ppm), pyrene (200 ppm); paraffin oil (1% v/v), and crude oil (1% v/v). Three consecutive subcultures were performed in the same conditions.
Each fungal and bacterial enrichment culture was plate on solid MM containing each pollutant as sole carbon source. To isolate pure microbial cultures, morphologically different colonies were selected and transferred to Malt Extract Agar (MEA) plates for fungi and Tryptone Soy Agar (TSA) plates for bacteria.
Molecular Characterization of Fungal and Bacterial Isolates
DNA from isolate purified colonies was extracted using Microlysis kit (Labogen, London, United Kingdom) according to the manufacturer’s protocol. The isolates were screened by randomly amplified polymorphic DNA-polymerase chain reaction (RAPD-PCR) using the single stranded oligonucleotides primer RAPD2 (5′-AGCAGCGTGG-3′) (Rebecchi et al., 2015) and GTG-5 (5′-GTG GTG GTG GTG GTG-3′) (Versalovic et al., 1994). The PCR fragment profiles were digitally captured using the BioImaging System Gene Genius and pattern analysis was performed with the Fingerprinting II software (Bio-Rad Laboratories, Hercules, CA, United States). The similarity in the profiles of bands was based on the Pearson correlation coefficient and the cluster analyses were performed by unweighted pair group method with arithmetic mean (UPGMA). A correlation coefficient of 70% was arbitrarily selected to distinguish the clusters, and one representative for each cluster was amplified using the primers P0 (5′-GAG AGT TTG ATC CTG GCT-3′) and P6 (5′-CTA CGG CTA CCT TGT TAC-3′) (Di Cello and Fani, 1996). PCR products were visualized by electrophoresis on 2.5% agarose gel in Tris-Acetate-EDTA (TAE). The PCR amplicons of approximately 1.5 kb, corresponding to the size of the full 16S rRNA gene were purified using the NucleoSpin gel and PCR clean-up according to the package insert (Macherey-Nagel, DE) and sequenced at the GATC Biotech (Germany). The taxonomical identification of sequences was performed using BLAST (Basic Local Alignment Search Tool) (Altschul et al., 1997) and by alignment against the Ribosomal Database Project (RDP) database using the Naïve Bayesian Classifier (Wang et al., 2007).
Genomic DNA of each strains was extracted from about 100 mg of mycelium scraped from the MEA petri dishes using the NucleoSpin® Plant II kit (Macherey-Nagel), according to the manufacturer’s instruction. The quality and quantity of extracted DNA was measured spectrophotometrically by Infinite M200 (TECAN Trading, Austria). De-replication of fungal isolates was performed by using the minisatellite core sequence derived from the wildtype phage M13 (5′-GAG GGT GGC GGT TCT-3′) as specific primer to amplify variable number tandem repeat (VNTR) (Poli et al., 2016). Molecular identification of each fungal isolate was carried out by amplification of specific markers (White et al., 1990; Glass and Donaldson, 1995; Bensch et al., 2012). PCR products were visualized by electrophoresis on a 1.5% agarose gel stained with ethidium bromide in Tris-Borate-EDTA (TBE). PCR products were purified and sequenced at Macrogen Europe (Amsterdam, Netherlands). Consensus sequences were obtained by using Sequencer 5.0 (Gene Code Corporation). Taxonomic assignments were inferred by querying with the Blastn algorithm (default setting), hosted at NCBI (National Center for Biotechnology Information), the newly generated sequences against the nucleotide database of NCBI (GenBank). Pairwise alignments were also performed against the CBS-Knaw Fungal Biodiversity Centre (Centraalbureau voor Schimmel cultures) database. Similarity values equal or higher than 98% (e-value > e-100) were considered reliable; results were confirmed morphologically.
DNA Extraction From Original Soils and Enrichments
Total microbial DNA was extracted from original soil samples and from each enrichment step of fungal and bacterial culture/microcosms with the PowerLyzer PowerSoil kit (MoBIO Laboratories, Inc., Carlsbad, CA, United States) according to the manufacturer’s instructions. DNA purity was checked with electrophoresis on a 0.8% agarose gel, while quantification was performed with the Quant-iTTM HS ds-DNA assay kit (Invitrogen, Paisley, United Kingdom) method in combination with the QuBitTM fluorometer.
Molecular Analyses of Microbial Diversity
Illumina Sequencing of 16S PCR Amplicons
For bacteria, PCR amplicons covering the V3–V4 regions of the 16S rRNA were analyzed in Illumina MiSeq with V3 chemistry in 300 bp paired-reads mode. PCR reactions were performed using indexed primer pairs 343F (5′-TACGGRAGGCAGCAG-3′) and 802R (5′-TACNVGGGTWTCTAATCC-3′), as described in Vasileiadis et al. (2015). A multiplexing strategy was employed to analyze several amplicon samples simultaneously in the same sequencing run. A nine nucleic acids extension was added to the 5′ end of the forward primer, where the first seven bases served as a tag to identify to each sample, and the following two bases were a linker designed not to match bacterial sequences in the same position according to RDP entries. In order to reduce possible biases related to the primer extension, the two step-PCR approach described in Berry et al. (2011) was adopted. The PCR conditions were set as described by Vasileiadis et al. (2015). The final PCR products were checked on 0.8% agarose gel and pooled in equimolar amounts according to QuBit measurements. The final pool was cleaned with the SPRI (Solid Phase Reverse Immobilization Method) using the Agencourt® AMPure® XP kit (Beckman Coulter, Milan, Italy). Finally, the pool was sequenced by BioFab Company (Rome, Italy) with a MiSeq Illumina instrument (Illumina, Inc., San Diego, CA, United States) operating with V3 chemistry and producing 300 bp paired-reads.
Illumina Sequencing of ITS PCR Amplicons
The nr ITS2 region was amplified from all DNA by means of a semi-nested PCR approach. In the first PCR, the nr ITS (ITS1-5.8S-ITS2) was amplified with universal primers ITS1F-ITS4 (White et al., 1990). For the second PCR, ITS3 and ITS4 (White et al., 1990) tagged primers were used to amplify the ITS2 region of each DNA sample (Voyron et al., 2017).
PCR products were pooled and purified using Wizard SV Gel and PCR Clean-Up System (Promega) following the manufacturer’s instructions. After quantification with Qubit 2.0 (Thermo Fisher Scientific, Waltham, MA, United States), the purified PCR products were mixed in equimolar amounts to prepare sequencing libraries. The libraries were paired-end sequenced using the Illumina MiSeq technology (2 bp × 250 bp) by IGA Technology Services S.r.l. Unipersonale (Udine, Italy).
Sequences Data Preparation, Bioinformatics, and Statistical Analyses
The first steps for sequences processing and filtering were the same for both 16S and ITS amplicons. Raw paired Illumina sequences were merged with the “pandaseq” script (Masella et al., 2012) with a minimum overlap of 30 bp between read pairs and 2 maximum allowed mismatches. Sequences were multiplexed according to sample indexes and primers with the fastx-toolkit1. Both bacterial and fungal amplicons were analyzed with taxonomy-based and OTU-based analyses: in the first case, all sequences were individually classified at taxonomical level against relevant database (GreenGenes for bacteria, UNITE for fungi), while in the OTU-based analyses, sequences were grouped at 97% similarities. For 16S amplicons, both operational taxonomic units (OTUs) and taxonomy-based matrixes were produced with a pipeline in Mothur (Schloss et al., 2009). For ITS amplicons, taxonomy-based analyses were also performed in Mothur, whereas OTUs were determined in UPARSE (Edgar, 2013). The reason for this discrepancy was that no aligned databases are available for ITS, and the OTU-clustering method implemented in Mothur does not allow analysis of sequences that have dissimilar lengths, which was the case with the ITS amplicons.
For bacterial sequences, Mothur v.1.39.5 (Schloss et al., 2009) was applied in order to remove sequences with large homopolymers (≥10), sequences that did not align within the targeted V3–V4 region, chimeric sequences (Edgar et al., 2011) and sequences that were not classified as bacterial after alignment against the Mothur version of the RDP training data set. The resulting high- quality sequences were analyzed with Mothur and R2 following the OTU and the taxonomy-based approach. Sequences were first aligned against the SILVA reference aligned database for bacteria (Pruesse et al., 2007) using the NAST algorithm and a kmer approach (Schloss, 2010) and then clustered at the 3% distance using the average linkage algorithm. OTUs were classified into taxa by alignment against the Greengenes database (McDonald et al., 2012).
ITS taxonomy-based analyses were conducted in Mothur: sequences shorter than 120 bp were discarded. We discarded homopolymers > 10 bp and chimeras, which were identified with the UCHIME algorithm implemented in Mothur, with the UNITE database version 6 as reference. The same database was used to classify the retained sequences and to eliminate non-fungal sequences. OTUs were produced in USEARCH with the –fastx_uniques and - cluster_otus commands. Sequences that did not belong to fungi were identified with the sintax command against the Utax reference database and discarded.
The OTU- and taxonomy-based matrixes obtained were analyzed in R to estimate the associated α and β diversity of the samples. The Good’s coverage estimate was calculated to assess the “percentage diversity” captured by sequencing. The most abundant OTUs identified were confirmed with BLAST (Basic Local Alignment Search) searches against the GenBank and the RDP database.
Hydrocarbon Degrading Bacterial and Fungal Strains
A collection of hydrocarbon-degrading bacterial and fungal strains was isolated from microcosms enriched on each pollutant (benzene, paraffin, pyrene, naphthalene, phenanthrene, and crude oil) and from three depths of the contaminated soil of the SIN in Fidenza. Molecular de-replication of isolated colonies resulted in 95 and 94 unique colonies for bacteria and fungi, respectively. Molecular identification of the strains was inferred by 16S rDNA and ITS sequencing and is reported in Table 1.
TABLE 1. Identification of bacteria and fungi isolated from benzene, paraffin, crude oil, naphthalene, pyrene, and phenanthrene enrichment microcosms.
For bacteria, sequence similarity searches in NCBI GenBank and RDP database revealed the presence of cultivable members of different genera under the phyla Proteobacteria, Firmicutes, and Sphingobacteria. On the whole, the bacteria isolated where ascribed mainly to 12 genera belonging to both Gram-negative and Gram-positive. Pseudomonas was the most abundant genera (56.8%) followed by Sphingobacterium (12%), Bacillus (6%), Stenothrophomonas (6%), Achromobacter (6%), and Serratia (3%). Most of the bacterial strains (46) were isolated from enrichments inoculated with soil S2 (Supplementary Table S2): the most represented bacterial genera in soil S2 were Pseudomonas (23 isolates), Serratia (3 isolates), and Sphingobacterium (9 isolates). A lower number of bacterial strains was isolated from enrichments of soil S1 and S3 (Supplementary Table S2). A total of 26 isolates were obtained from soil S1 and were mainly affiliated to the genera Bacillus, Pseudomonas and Pseudoxanthomonas, while the 23 bacterial strains isolated from the deepest soil S3 were mainly affiliated to the genus Pseudomonas (Supplementary Table S2). The influence of pollutants used in the enrichments on the distribution of bacterial species is reported in Table 1. The highest number of bacterial isolates was obtained from paraffin (29 isolates), followed by naphthalene (19 isolates), benzene (15 isolates), and phenanthrene (15 isolates) microcosms; 10 strains were isolated in the presence of pyrene and only six strains were isolated in the presence of crude oil (Table 1). The genus Pseudomonas dominated all the microcosms enriched on pollutants. Benzene and paraffin enrichments showed the highest diversity of bacterial genera. In particular, benzene enrichment selected also bacterial species of the genera Sphingobacterium and Stenotrophomonas, while paraffin enrichment resulted also in the isolation of species of Bacillus, Achromobacter, and Pseudoxanthomonas. Beside Pseudomonas, Serratia, and Ochromobactrum were the most abundant genera in crude oil and pyrene microcosms (Table 1). At species level, Pseudomonas putida was dominant and was isolated in the presence of all the tested pollutants. Other frequently isolated species were Pseudomonas fluorescence (four pollutants), Stenotrophomonas maltophilia (two pollutants), Stenotrophomonas acidaminiphila (two pollutants), Bacillus subtilis (two pollutants), and Sphingobacterium multivorum (two pollutants). Three strains of Serratia marcescens were isolated only from the enriched microcosms with the crude oil of the contaminated site (Table 1).
Fungal species with their relative abundance are reported in Supplementary Table S3, according to the soil depth and in Table 1 according to the pollutant used in the enrichments. Isolated fungi belonged mostly to the phylum Ascomycota, and only 5% to Basidiomycota; Mucoromycota were not isolated. Most of the strains were isolated from enrichments inoculated with contaminated soils S1 (35 strains) and S3 (36 strains); a lower number of isolates (23) were retrieved from the contaminated soil S2 (Supplementary Table S3). In general, Fusarium was the most abundant genus (43.8%) followed by Trichoderma (and its teleomorph Hypocrea lixii) (13.8%) and Aspergillus (11.0%). The majority of the isolates detected in soils S1 and S3 were indeed ascribable to these three genera, while in S2 only 44.4% of the strains were affiliated to the genus Fusarium and 16.7% belonged to the genus Cladosporium (Supplementary Table S3).
The retrieval of fungal species was influenced by the pollutant used as sole carbon source in the enrichments, as reported in Table 1. The highest number of isolates was obtained from pyrene and phenanthrene (23-22 isolates) followed by benzene and crude oil (15 isolates); 11 strains were isolated in the presence of naphthalene and only eight strains were isolated in the presence of paraffin oil (Table 1). With the only exception of crude oil, the genus Fusarium was the most abundant in each pollutant. Benzene and paraffin enrichments revealed the highest variety of genus and species.
Molecular Analyses of Bacterial and Fungal Communities
Illumina HTS of PCR amplicons of the 16S V3–V4 region resulted in 1,991,866 sequences, which were reduced to 1,742,447 after the exclusion of homopolymers, of sequences <380 bp, sequences not aligning to the targeted V3–V4 regions, chimera and sequences not classified as bacterial after alignment against the RDP 16S training set. In order to reduce biases in diversity estimates related to the analyses of samples with different number of sequences, rarefaction to a common minimum number of 10,065 sequences per sample was carried out (Gihring et al., 2012). The rarefaction step resulted in the loss of only 5 out of 82 samples, while retaining a sequencing depth able to depict most of the bacterial diversity in the samples, as indicated by an average Good’s coverage of 90.4%.
As for fungi, Illumina sequencing of ITS2 amplicons resulted in 795,742 high quality sequences. Taxonomy-based analyses were performed in Mothur and consisted in a screening to remove homopolymers, sequences shorter than 120 bp, chimera, and sequences that were not correctly classified as fungal after alignment against the UNITE v6 database. Rarefaction was performed with a common number of 5,001 sequences per sample. The coverage of fungal diversity was not affected by this rarefaction step: average Good’s coverage was 99.7%, thus showing that the sequencing effort covered the totality of fungal biodiversity in the samples. OTU-based analyses for fungi were performed in USEARCH and resulted in a very similar number of sequences per sample (4,942) and a total of 188 fungal OTUs.
Multivariate canonical correspondence analysis (CCA) was performed on bacterial and fungal OTUs abundance tables in order to assess how the microbial communities responded to contaminants, depth of soil collection and enrichment time (Figure 1). Among the tested factors, in the case of bacteria, pollutant was the most significant: the percentage of variance explained (36.1%, Figure 1A) was much higher than the variance explained by the soil (6.1%, Figure 1B) and by the time (8.3%, Figure 1C) factors. The bacterial communities differed among the six tested contaminants: pollutants that shared the chemical structure (e.g., the three PAHs) grouped together, while oil, benzene, and paraffin formed separate groups (Figure 1A). A similar picture was obtained by the CCA analyses performed on fungal communities: only the pollutant factor was significant (25.6% of variance explained, Figure 1D) contrary to soil and time factors. The pollutant used in the enrichment significantly affected the evolution of fungal communities: phenanthrene and pyrene enriched fungal communities formed groups that separated from those enriched from oil (Figure 1D). For both fungi and bacteria, the groups formed by soil-depth (Figures 1B–E) and time (Figures 1C–F) were mainly overlapping: this indicates that the depth, as well as the time of enrichment, did not affect the evolution of enriched microbial communities. Grouping of samples was also evaluated by UPGMA clustering of sequencing classified at genus level. Results are reported in Figure 2 for bacteria and in Figure 3 for fungi: samples are labeled according to the soil depth, followed by the pollutant and the enrichment step. Results confirmed that the relative composition in bacterial genera was partly dependent on the pollutant: several pyrene and phenanthrene samples formed two sub-clusters in a common cluster, with Azospirillum, Achromobacter, and Pseudomonas being the dominant genera; naphthalene enrichments were dominated by Achromobacter and Pseudomonas in similar proportions, while oil enrichments showed a dominance of Pseudomonas followed by Achromobacter. Finally, benzene and paraffin enrichments revealed a more diverse community, with high relative presences of Acinetobacter, Pseudoxanthomonas, and Pseudomonas. Neither the soil depth nor the time of enrichment influenced the samples clustering, confirming the CCA results presented above.
FIGURE 1. Canonical correspondence analyses (CCAs) to test the significance of the effects of pollutant (A,D), soil depth (B,E), and time (C,F) on the total structure of bacterial and fungal communities as determined by the relative abundances of all the OTUs resulting by Illumina sequencing of bacterial 16S and fungal ITS amplicons. (A–C) Graphs in the first line are for bacteria, (D–F) are for fungi.
FIGURE 2. Hierarchical clustering of bacterial 16S sequences classified at the genus level. Only taxa participating with >20% in at least one sample are shown, while taxa with lower participation are grouped in the “other” sequence group.
FIGURE 3. Hierarchical clustering of bacterial ITS sequences classified at the genus level. Only taxa participating with >5% in at least one sample are shown, while taxa with lower participation are grouped in the “other” sequence group.
Concerning fungi, no clear clustering of genera was shown according to pollutant, soil depth, or enrichment time (Figure 3). All samples were dominated by Fusarium (more than 80% in average), followed by Aspergillus, Penicillium, Trichoderma, and Arthrinium.
Comparison Between Culture-Dependent and Culture-Independent Analyses
A number of analyses were performed in order to compare the results on the diversity of bacterial and fungal strains isolated with culture-dependent and culture-independent methods. In Table 2 the relative percentages of bacterial genera retrieved from each of the three soil depths are compared. Results generally agree: Pseudomonas was the dominant genus in both cases, with 24.2, 25.2, and 63.5% of relative abundances in Illumina and 57.7, 48.9, and 69.6% as isolates for S1, S2, and S3, respectively. The other genera, in terms of relative percentages, were Bacillus, Achromobacter, Sphingobacterium, and Pseudotrophomonas for the culture-dependent data; Achromobacer, Agrobacterium, Azospirilllum, Shinella, Sphingobacterium, and Pseudotrophomonas for the Illumina data. It is worth noting that for some genera, the soil depth with the highest relative percentage was the same with both approaches. For instance, Pseudomonas reached 60% of abundance in soil three from both the approaches. There are a number of discrepancies: Agrobacterium sequences were above 9% in all soils while no strains belonging to this genus was isolated; and 2.2% of isolates (i.e., two strains) from S2 belong to Helicobacter while no sequences affiliated to this genus were detected with Illumina.
TABLE 2. Relative abundances of bacterial genera determined by Illumina sequencing of 16S amplicns or by isolation on selective media; the comparison was performed at the end of the enrichment (step III).
The same analyses were conducted on fungi (Table 3). Sordariomycetes was the most abundant class, with the order Hypocreales (i.e., Acremonium, Clonostachys, Fusarium, Hypocrea/Trichoderma) and Microascales (i.e., Ceriosporopsis, Pseudoallescheria/Scedosporium). Xylariales were instead a minor component of the detected fungi (i.e., Arthrinium, Eutypella). Results from culture-dependent and culture-indepentent methods confirmed each other in the case of the genus Fusarium, whose relative percentages were highly similar (ranging 40 and 69%). Clonotachys, Pseudallescheria, and Trichoderma were mostly isolated from S1 soils where also Illumina data showed the highest relative percentage. In the case of fungi too, there are some contrasting results, as fungal genera isolated in pure culture that were not detected with Illumina sequencing of ITS2 amplicons (e.g., Bjerkandera, Epicoccum, Eutypella, Irpex, Polyporus, and Sulcatispora). As for the genera Ceriosporopsis and Phanerochaete, they were detected with Illumina in significant amounts (above 5% of total sequences) but no isolates were obtained.
TABLE 3. Relative abundances of fungal genera as determined by Illumina sequencing of 16S PCR of amplicons or by isolation on selective media; the comparison was performed at the end of the enrichment (step III).
The most abundant OTUs (i.e., those who had a relative abundance higher than 1% in at least one sample) are reported in Tables 4, 5 for bacteria and fungi, respectively. This screening criterion retained 19 OTUs for bacteria, and 21 OTUs for fungi: in both tables, the highest taxonomic affiliation is reported, together with the relative abundances in each enrichment microcosm. In the case of bacteria, most OTUs (16 out of 19) were classified at species level, two at genus and only one at the family level. Nine out of 19 OTUs were retrieved in the cultivated isolates. By comparing the data on the isolates in Table 1 with the molecular data in Table 4, it is possible to notice that the species comprising more isolates (Pseudomonas putida, Achromobacter xylosoxidans, and Pseudomonas veronii) were also the most abundant in terms of Illumina OTUs. As above reported for genera, also for OTUs classified at genus or species level there was a good agreement between molecular and microbiological data.
TABLE 4. Relative percentages of most abundant bacterial OTUs after the last enrichment step for each pollutant and at each of the three soil depths tested.
TABLE 5. Relative percentages of most abundant fungal OTUs after the last enrichment step for each pollutant and at each of the three soil depths tested. OTUs that represent more than 1% of total fungal community in at least one sample are presented.
Regarding fungi, the rate of OTUs that could be classified at species level is similar to the one for bacteria: 14 out of 21. The number of OTUs also detected among isolates was instead slightly lower, 9 isolates (Table 5). It is worth noting that, as for bacteria, the taxa comprising more isolates were also the most abundant OTUs: that was the case for Fusarium oxysporum, Cladosporium spp., and Fusarium solani. These data highlighted the different adaptation abilities developed by fungi: while species of the genus Fusarium were abundant and widely distributed among pollutants and depth, other fungi seemed to be particularly adapted to individual contaminant. For instance, the relative percentage of Cladosporium spp. and Pseudoallescheria spp. was significantly high only in the presence of oil and benzene, respectively.
Since the discovery of the great plate count anomaly in the middle ‘80s of the last century (Staley and Konopka, 1985), microbiologists are facing the issue of non-cultivable microorganisms and the underestimation of the actual microbial community by means of direct isolation techniques. This has become more and more relevant with the advent of HTS techniques, and it is of utmost importance in the environmental microbiology sector of bioremediation, where studies must rely on the isolation and identification of bacterial and fungal strains to be inoculated in soil for efficient degradation of target pollutants (Megharaj et al., 2011). These strains are usually obtained through enrichment procedures, aimed at isolating autochthonous microorganisms capable of using the pollutants as sole carbon sources (Thompson et al., 2005). Although many data are available in literature about the assessment of microbial community dynamics following the addition of enriched microbial strains (MacNaughton et al., 1999; Aburto-Medina et al., 2012; Fuentes et al., 2014), little is known on the developments that occur within the bacterial and fungal communities during the enrichment steps. Furthermore, to the best of our knowledge, no attempt was made to assess how soil samples with a high petroleum hydrocarbons content and collected at different depths affect the organization of the bacterial and fungal communities.
In this work, we assessed whether the pollutant and the collection depth had significant effects on the microbial community arrangements. Enrichments were performed using five pollutants (three PAHs, benzene, and paraffin) and a complex mixture (i.e., crude oil collected from the polluted site). The results obtained (Figure 1) highlighted that for both bacteria and fungi, the pollutant partially shaped the enriched communities. Although it was difficult to highlight separate clusters according to the pollutants, some differences were observed as bacteria grown in the presence of paraffin and benzene or with PAHs, and fungi grown with paraffin or PAHs (Figure 1). This is not surprising, if we consider that different metabolic pathways are involved in the degradation of aliphatic or aromatic hydrocarbons, but that single strains can be also equipped with both pathways (Whyte et al., 1997). The microbial community responded uniquely to the presence of crude oil: only those strains, highly adapted to this extreme and toxic environment, were capable of colonizing this ecological niche.
According to literature, a relation between pollutants load and microbial communities can be drawn. Several studies assessed a different distribution of PAHs and other oil pollutants along soil depths (Cousins et al., 1999; Ping et al., 2007). It is well-known that microbial community organization is affected by soil depths, commonly causing a decline in biomass and diversity from the surface downward (Blume et al., 2002; Fierer et al., 2003; Oehl et al., 2005). As far as we know, no experiments were previously conducted to assess the impact of the soil samples depths on the community structure and enrichment outcomes. Contrary to this general background, in the present study, neither the soil depth nor the enrichment step had significant effects on the community structure (Figures 1–3). Regarding fungi, this could be explained by their physiological growing capabilities. Thanks to their hyphae growth, they can widely colonize the soil, creating a homogenous mycoflora in the space (Hesham et al., 2017).
The total load of bacteria and fungi was lower than unpolluted soils (Supplementary Table S1), although no differences were observed among the collection depths: a strong and long-lasting pollution of this site may shape a microbial community adapted to the pollutants pressure mostly composed by oil-utilizing strains. As illustrated in Figures 2, 3, differences among samples are more evident for bacteria. However, it is not possible to define clear clusters based on soil depths; the structure of the microbial community was not dependent by the collection point. Since the microbial community is selected by the total content of the pollutants, we can speculate that, in this soil, a pollutants threshold has been reached in the surface topsoil: beyond this level, only the selected and adapted microflora can survive. Indeed, in a less polluted site, the presence of both aerobic and anaerobic microorganisms decreased with the vertical soil profile (Biró et al., 2014). Regarding the enrichment steps, it was surprising to find that already after 1 week of enrichment, the communities were similar to those found after 4 weeks. This indicates that when dealing with strong polluted soils, as in this study, a single enrichment step may be enough.
A total of 95 bacterial and 94 fungal strains were identified. These strains were well-distributed along soil depths and pollutant (Supplementary Tables S2, S3). The taxonomical identification of bacteria confirmed the selection of specific bacterial communities following the addition of different hydrocarbons, indicating a high specialization of bacterial taxa involved in their degradation. Proteobacteria, Firmicutes, and Sphingobacteria were the prevailing bacterial phyla: all these phyla comprised several phylogenetic groups involved in the aerobic degradation of hydrocarbons (Vinas et al., 2005; Head et al., 2006; Yang et al., 2012; Yergeau et al., 2012). Proteobacteria was also found as the most abundant phylum distributed in an environment contaminated with petroleum muck and in activated biomass from a petrochemical industry wastewater sample (Joshi et al., 2014; Yadav et al., 2015). This phylum is divided into five major classes, all with oil-degrading genera (Gao et al., 2014). On the whole, isolated bacteria where ascribed mainly to genera belonging to Gamma-proteobacteria class, whose dominance within the enriched microcosms can be expected since it comprises many of the hydrocarbonoclastic bacteria responsible for the first steps of hydrocarbons degradation (Head et al., 2006; Yakimov et al., 2007; Barbato et al., 2016).
Bacterial species of the genera Sphingobacterium, Achromobacter, Pseudomonas, Stenotrophomonas, and Bacillus, isolated from this crude-oil contaminated soil, showed close lineage with previously reported hydrocarbons degrading members of same genera and most of them have been reported earlier as hydrocarbon degraders (Widdel and Rabus, 2001; Roy et al., 2002; Foght, 2008; Chandra et al., 2013; Varjani et al., 2015; Varjani and Upasani, 2016). Members of the genus Pseudomonas are widespread environmental microorganisms isolated from a variety of natural sources (Velmurugan et al., 2011) and well-known biodegraders, capable of metabolizing a range of compounds (Kostka et al., 2011; Mulet et al., 2011; Mahjoubi et al., 2013).
According to our results, the PAHs-degrading isolates in soil belong to the Sphingomonas and Pseudomonas bacteria. Pseudomonas usually contain nah-like genes that encode for the dioxygenase subfamily; the products of these genes are capable of degrading low-molecular- weight PAHs (Johnsen et al., 2007) and producing a rhamnolipid biosurfactant in soils contaminated with petroleum hydrocarbons (Sadouk et al., 2008; Xia et al., 2014). Unlike other gram-negative bacteria strains, members of the genus Sphingomonas are able to degrade a wide range of natural and xenobiotic compounds (Eguchi et al., 1996; Pinhassi and Hagstrom, 2000). Three strains of Serratia marcescens were also isolated from the enriched microcosms with the crude oil of the contaminated site, supporting previous evidences of degradation of crude oil and petroleum products by species of Serratia (De La Fuente et al., 1991; Rojas-Avelizapa et al., 2002; Wongsa et al., 2004; Rajasekar et al., 2007).
Generally, fungi involved in the degradation of PAHs include ligninolytic and non-ligninolytic fungi. Most of the strains isolated after the enrichment steps belonged to the phylum Ascomycota, although few Basidiomycota were isolated. The dominance of Ascomycota in polluted soil has been extensively acknowledged: Zhang et al. (2018) reported that Ascomycota represented up to 73–96% of the total 18S sequences of a coking area soil; eight of the 10 strains isolated from a PAHs contaminated pond were Ascomycetes and were also able to remove anthracene (Aranda et al., 2017). Finally, an oilfield with a history of 50 years pollution influenced mainly the bacterial community, whereas fungi, mostly related to this phylum, were abundant, demonstrating to be less sensitive to soil PAHs (Zhou et al., 2017).
Basidiomycota represented only the 5.3% of the total isolates. Similarly, in a toxic coking area, only a small fraction of the isolates belonged to Basidiomycota (0.34–7.0%) (Zhang et al., 2018). Besides Bjerkandera adusta, Irpex lacteus, Wallemia mellicola, Polyporus gayanus, Trametes gibbosa, the presence of species of Lopharia, Phanerochaete, and Trichosporon was assessed by HTS. Among them, T. gibbosa, P. gayanus, and W. mellicola were isolated for the first time in hydrocarbons polluted soil.
On the other hand, Mucoromycota were absent in this contaminated soil. Unlike Morais et al. (2016) who found more than 50% of Mucoromycota in a petroleum contaminated costal site, we did not identified any organism belonging to this phylum. Mucoromycota represented a small fraction of the isolates also in an estuarine sediment contaminated mainly by PAHs (Da Silva et al., 2003) and in heavy crude oil-contaminated soil (Zafra et al., 2014). In the present study, Mucoromycota could be included in the “unclassified” group (Figure 3), even though they would be a negligible fraction of the mycoflora.
All samples were dominated by Fusarium, followed by Aspergillus, Penicillium, and Trichoderma. The identification of these genera as polluted soil inhabitants has been already discussed elsewhere. Sordariomycetes is the most represented class with Acremonium, Arthrinium, Clonostachys, Eutypella, Fusarium, Hypocrea/Trichoderma, and Pseudoallescheria/Scedosporium. Fungal taxonomy analysis of an oil-contaminated soil detected almost 25% of Sordariomyecetes (Morais et al., 2016). Trichoderma harzianum (Da Silva et al., 2003), Scedosporium apiospermum and Acremonium sp. (Chaıneau et al., 1999; Zafra et al., 2014), Fusarium oxysporum (Marchand et al., 2017), and Fusarium solani (Chaıneau et al., 1999) were also isolated from contaminated soil or sediments. Fusarium (three strains), Trichoderma (three strains), and Pseudallescheria (two strains) represented the 40% of the isolates from heavy and extra-heavy crude oil polluted soils (Naranjo et al., 2007). Zhou et al. (2017) detected more than 20% of Gibberella species in an oilfield site. Since Gibberella is the teleomorph of the majority of Fusarium species (Leslie and Summerell, 2008) we can state that this finding is in agreement with our observations. The isolation of 42 strains of Fusarium, among which 64% were ascribable to F. solani, opens to intriguing solutions for their application in bioremediation processes. Indeed F. solani isolated from oil polluted soil by enrichment method using phenanthrene as the sole source of carbon was capable of degrading a mixture of low and high molecular weight PAHs (up to 84.8%) (Hesham et al., 2017). To the best of our knowledge, this is the first time that a strain of Eutypella scoparia has been isolated from a polluted soil. Previous researches reported the isolation of this species only in a non-polluted Arctic soil (Bergero et al., 1999; Liu et al., 2014).
Surprisingly, Morais et al. (2016) did not found Eurotiomycetes in a contaminated soil. On the contrary, in the present study, the fungal community was rich of Aspergillus and Penicillium. Aspergillus and Penicillium strains have been already isolated in several polluted samples (Chaıneau et al., 1999; Da Silva et al., 2003; Zafra et al., 2014; Al-Hawash et al., 2018). Data on soil isolation of Aspergillus niger (Chaıneau et al., 1999), Aspergillus fumigatus (Chaıneau et al., 1999; Zafra et al., 2014), Aspergillus terreus (Da Silva et al., 2003), Aspergillus flavus, Aspergillus niger, Aspergillus nomius (Zafra et al., 2014), Aspergillus sydowii (Tian et al., 2016), Aspergillus versicolor (Šimonovièová et al., 2017) have been reported although few information deals with the here-detected species of Aspergillus (i.e., A. creber, A. jensenii, A. protuberus). For instance, A. creber and A. jensenii, described as new species in Jurjevic et al. (2012), are mostly known as air-born contaminants or of clinical relevance (Micheluz et al., 2016; Siqueira et al., 2016); A. protuberus has been isolated from air and the mussel Mytilus edulis galloprovincialis (Siqueira et al., 2016; Klas et al., 2018). There is no evidence in literature of these species as a terrestrial inhabitant. Although already detected in soil, Aspergillus sclerotiorum and Aspergillus waskmanii have been only associated to not-contaminated lands (Nesci et al., 2006; Phainuphong et al., 2017) or as part of the rhizosphere of tussock (Di Menna et al., 2007). Although this is the first time a Penicillium catenatum is found in a polluted soil, it must be considered that the identification of soil fungi is not always accurate at species level (Al-Hawash et al., 2018).
A thorough comparison was performed between the relative abundances of bacterial and fungal strains. The comparisons carried out at genus level (Table 2 for bacteria, Table 3 for fungi) and at OTU level (Table 4 for bacteria, Table 5 for fungi) showed an agreement between the culture-dependent and the culture-independent approaches. In most cases, species level was reached, thus confirming that the Illumina sequencing of the two hypervariable V3–V4 regions and of the partial ITS2 are sufficient to reach the species level for bacteria and fungi respectively. The species with more isolates (such as P. putida, A. xylosoxidans and Ochromobactrum anthropi for bacteria, F. oxysporum and F. solani for fungi) were also the dominant OTUs, thus showing that the culturing methods reflect the relative composition of total communities. Furthermore, the two approaches also agree for the less abundant genera: this is expected since enrichment procedures indeed select for the cultivable species among the total bacterial and fungal communities that is originally present in the polluted soil. A novel aspect shown in our work is that the two approaches had detection levels that were generally in good agreement. Composition of the actual microbial community and are not much subjected to biases.
In the present study, we evaluated the evolution of bacterial and fungal communities enriched from polluted soil by culture-independent and dependent methods. The results showed that already after 1 week of enrichment, both bacterial and fungal communities resembled those found after 4 weeks, and that the soil depth did not influence the evolution of microbial communities, contrary to the pollutant used. Isolation results were in agreement with HTS molecular data, indicating that the strains isolated reflected the microbial composition of the enriched consortia. Future studies will focus on the biodegradation abilities of the isolates microorganisms.
EP and GV ideated and supervised the study. GS, FS, and AP performed the experiments and wrote the paper. A-LB and TR performed the experiments. CG carried out chemical analyses.
This research was conducted within the LIFE15 ENV/IT/000396 project “Bioremediation and revegetation to restore the public use of contaminated land.”
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors are grateful to Sotirios Vasileiadis for a number of scripts used for HTS data elaboration, and to Prof Mariangela Girlanda and Dr. Samuele Voyron for helping in Fungal metabarcoding analysis.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.02543/full#supplementary-material
Aburto-Medina, A., Adetutu, E. M., Aleer, S., Weber, J., Patil, S. S., Sheppard, P. J., et al. (2012). Comparison of indigenous and exogenous microbial populations during slurry phase biodegradation of long-term hydrocarbon-contaminated soil. Biodegradation 23, 813–822. doi: 10.1007/s10532-012-9563-8
Agnello, A. C., Bagard, M., Van Hullebusch, E. D., Esposito, G., and Huguenot, D. (2016). Comparative bioremediation of heavy metals and petroleum hydrocarbons co-contaminated soil by natural attenuation, phytoremediation, bioaugmentation and bioaugmentation-assisted phytoremediation. Sci. Total. Environ. 563-564, 693–703. doi: 10.1016/j.scitotenv.2015.10.061
Al-Hawash, A. B., Alkooranee, J. T., Abbood, H. A., Zhang, J., Sun, J., Zhang, X., et al. (2018). Isolation and characterization of two crude oil-degrading fungi strains from Rumaila oil field. Iraq. Biotechnol. Rep. 17, 104–109. doi: 10.1016/j.btre.2017.12.006
Alonso-Gutierrez, J., Figueras, A., Albaiges, J., Jimenez, N., Vinas, M., Solanas, A. M., et al. (2009). Bacterial communities from shoreline environments (Costa da morte, northwestern Spain) affected by the prestige oil spill. Appl. Environ. Microbiol. 75, 3407–3418. doi: 10.1128/AEM.01776-08
Altschul, S. F., Madden, T. L., Schaffer, A., Zhang, J., Zhang, Z., Miller, W., et al. (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic. Acid. Res. 25, 3389–3404. doi: 10.1093/nar/25.17.3389
An, D., Caffrey, S. M., Soh, J., Agrawal, A., Brown, D., Budwill, K., et al. (2013). Metagenomics of hydrocarbon resource environments indicates aerobic taxa and genes to be unexpectedly common. Environ. Sci. Technol. 47, 10708–10717. doi: 10.1021/es4020184
Aranda, E., Godoy, P., Reina, R., Badia-Fabregat, M., Rosell, M., Marco-Urrea, E., et al. (2017). Isolation of Ascomycota fungi with capability to transform PAHs: insights into the biodegradation mechanisms of Penicillium oxalicum. Int. Biodeter. Biodegr. 122, 141–150. doi: 10.1016/j.ibiod.2017.05.015
Badr El-Din, S. M., Moussa, T. A., Moawad, H., and Sharaf, O. A. (2014). Isolation and characterization of polyaromatic hydrocarbons degrading bacteria from compost leachate. J. Adv. Biol. 5, 651–660.
Bao, M. T., Wang, L. N., Sun, P. Y., Cao, L. X., Zou, J., and Li, Y. M. (2012). Biodegradation of crude oil using an efficient microbial consortium in a simulated marine environment. Mar. Pollut. Bull. 64, 1177–1185. doi: 10.1016/j.marpolbul.2012.03.020
Barbato, M., Mapelli, F., Magagnini, M., Chouaia, B., Armeni, M., Marasco, R., et al. (2016). Hydrocarbon pollutants shape bacterial community assembly of harbor sediments. Mar. Pollut. Bull. 104, 211–220. doi: 10.1016/j.marpolbul.2016.01.029
Berry, D., Mahfoudh, K. B., Wagner, M., and Loy, A. (2011). Barcoded primers used in multiplex amplicon pyrosequencing bias amplification. Appl. Environ. Microb. 77, 7846–7849. doi: 10.1128/AEM.05220-11
Biró, B., Toscano, G., Horváth, N., Matics, H., Domonkos, M., Scotti, R., et al. (2014). Vertical and horizontal distributions of microbial abundances and enzymatic activities in propylene-glycol-affected soils. Environ. Sci. Pollut. R. 21, 9095–9108. doi: 10.1007/s11356-014-2686-1
Blume, E., Bischoff, M., Reichert, J., Moorman, T., Konopka, A., and Turco, R. (2002). Surface and subsurface microbial biomass, community structure and metabolic activity as a function of soil depth and season. Appl. Soil Ecol. 20, 171–181. doi: 10.1016/S0929-1393(02)00025-2
Bundy, J. G., Paton, G. I., and Campbell, C. D. (2004). Combined microbial community level and single species biosensor responses to monitor recovery of oil polluted soil. Soil Biol. Biochem. 36, 1149–1159. doi: 10.1016/j.soilbio.2004.02.025
Cerqueira, V. S., Hollenbach, E. B., Maboni, F., Vainstein, M. H., Camargo, F. A., Do Carmo, R. P. M., et al. (2011). Biodegradation potential of oily sludge by pure and mixed bacterial cultures. Bioresour. Technol. 102, 11003–11010. doi: 10.1016/j.biortech.2011.09.074
Chaıneau, C., Morel, J., Dupont, J., Bury, E., and Oudot, J. (1999). Comparison of the fuel oil biodegradation potential of hydrocarbon-assimilating microorganisms isolated from a temperate agricultural soil. Sci. Tot. Environ. 227, 237–247. doi: 10.1016/S0048-9697(99)00033-9
Chandra, S., Sharma, R., and Singh, K. (2013). Application of bioremediation technology in the environment contaminated with petroleum hydrocarbon. Ann. Microbiol. 63, 417–431. doi: 10.1007/s13213-012-0543-3
Cousins, I. T., Gevao, B., and Jones, K. C. (1999). Measuring and modelling the vertical distribution of semi-volatile organic compounds in soils. I: PCB and PAH soil core data. Chemosphere 39, 2507–2518. doi: 10.1016/S0045-6535(99)00164-2
Covino, S., Stella, T., and Cajthaml, T. (2016). “Mycoremediation of organic pollutants: principles, opportunities, and pitfalls,” in Fungal Applications in Sustainable Environmental Biotechnology, ed. D. Purchase (Berlin: Springer), 185–231. doi: 10.1007/978-3-319-42852-9_8
Da Silva, M., Umbuzeiro, G. A., Pfenning, L. H., Canhos, V. P., and Esposito, E. (2003). Filamentous fungi isolated from estuarine sediments contaminated with industrial discharges. Soil Sediment Contam. 12, 345–356. doi: 10.1080/713610976
Dean-Ross, D., Moody, J., and Cerniglia, C. E. (2002). Utilization of mixtures of polycyclic aromatic hydrocarbons by bacteria isolated from contaminated sediment. FEMS Microbiol. Ecol. 41, 1–7. doi: 10.1111/j.1574-6941.2002.tb00960.x
Di Menna, M., Sayer, S., Barratt, B., Bell, N., Ferguson, C., and Townsend, R. (2007). Biodiversity of indigenous tussock grassland sites in Otago, Canterbury and the central North Island. V. Penicillia and Aspergilli. J. Roy. Soc. New Zeal. 37, 131–137. doi: 10.1080/03014220709510541
Dueholm, M. S., Marques, I. G., Karst, S. M., D’imperio, S., Tale, V. P., Lewis, D., et al. (2015). Survival and activity of individual bioaugmentation strains. Bioresour. Technol. 186, 192–199. doi: 10.1016/j.biortech.2015.02.111
Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C., and Knight, R. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200. doi: 10.1093/bioinformatics/btr381
Eguchi, M., Nishikawa, T., Macdonald, K., Cavicchioli, R., Gottschal, J., and Kjelleberg, S. (1996). Responses to stress and nutrient availability by the marine ultramicrobacterium Sphingomonas sp. strain RB2256. App. Environ. Microb. 62, 1287–1294.
Fernandez-Luqueno, F., Valenzuela-Encinas, C., Marsch, R., Martinez-Suarez, C., Vazquez-Nunez, E., and Dendooven, L. (2010). Microbial communities to mitigate contamination of PAHs in soil-possibilities and challenge: a review. Environ. Sci. Pollut. R. 18, 12–30. doi: 10.1007/s11356-010-0371-6
Fierer, N., Schimel, J. P., and Holden, P. A. (2003). Variations in microbial community composition through two soil depth profiles. Soil Biol. Biochem. 35, 167–176. doi: 10.1016/S0038-0717(02)00251-1
Fowler, S. J., Toth, C. R., and Gieg, L. M. (2016). Community structure in methanogenic enrichments provides insight into syntrophic interactions in hydrocarbon-impacted environments. Front. Microbiol. 7:562. doi: 10.3389/fmicb.2016.00562
Fuentes, S., Méndez, V., Aguila, P., and Seeger, M. (2014). Bioremediation of petroleum hydrocarbons: catabolic genes, microbial communities, and applications. Appl. Microbiol. Biotechnol. 98, 4781–4794. doi: 10.1007/s00253-014-5684-9
Gao, Y. C., Guo, S. H., Wang, J. N., Li, D., Wang, H., and Zeng, D. H. (2014). Effects of different remediation treatments on crude oil contaminated saline soil. Chemosphere 117, 486–493. doi: 10.1016/j.chemosphere.2014.08.070
Gihring, T. M., Green, S. J., and Schadt, C. W. (2012). Massively parallel rRNA gene sequencing exacerbates the potential for biased community diversity comparisons due to variable library sizes. Environ. Microbiol. 14, 285–290. doi: 10.1111/j.1462-2920.2011.02550.x
Hara, E., Kurihara, M., Nomura, N., Nakajima, T., and Uchiyama, H. (2013). Bioremediation field trial of oil-contaminated soil with food-waste compost. J. JSCE 1, 125–132. doi: 10.2208/journalofjsce.1.1_125
Head, I. M., Gray, N. D., and Larter, S. R. (2014). Life in the slow lane; biogeochemistry of biodegraded petroleum containing reservoirs and implications for energy recovery and carbon management. Front. Microbiol. 5:566. doi: 10.3389/fmicb.2014.00566
Hesham, A. E.-L., Mohamed, E. A., Mawad, A. M., Elfarash, A., Abd El-Fattah, B. S., and El-Rawy, M. A. (2017). Molecular characterization of Fusarium solani degrades a mixture of low and high molecular weight polycyclic aromatic hydrocarbons. Biotechnol. J. 11, 27–35.
Johnsen, A. R., Schmidt, S., Hybholt, T. K., Henriksen, S., Jacobsen, C. S., and Andersen, O. (2007). Strong impact on the polycyclic aromatic hydrocarbon (PAH)-degrading community of a PAH-polluted soil but marginal effect on PAH degradation when priming with bioremediated soil dominated by mycobacteria. Appl. Environ. Microbiol. 73, 1474–1480. doi: 10.1128/AEM.02236-06
Joshi, M. N., Dhebar, S. V., Dhebar, S. V., Bhargava, P., Pandit, A., Patel, R. P., et al. (2014). Metagenomics of petroleum muck: revealing microbial diversity and depicting microbial syntrophy. Arch. Microbiol. 196, 531–544. doi: 10.1007/s00203-014-0992-0
Jurjevic, Z., Peterson, S. W., and Horn, B. W. (2012). Aspergillus section versicolores: nine new species and multilocus DNA sequence based phylogeny. IMA Fungus 3, 59–59. doi: 10.5598/imafungus.2012.03.01.07
Klas, K. R., Kato, H., Frisvad, J. C., Yu, F., Newmister, S. A., Fraley, A. E., et al. (2018). Structural and stereochemical diversity in prenylated indole alkaloids containing the bicyclo-[2.2.2]-diazaoctane ring system from marine and terrestrial fungi. Nat. Prod. Rep. 35, 532–558. doi: 10.1039/c7np00042a
Kostka, J. E., Prakash, O., Overholt, W. A., Green, S. J., Freyer, G., Canion, A., et al. (2011). Hydrocarbon-degrading bacteria and the bacterial community response in gulf of Mexico beach sands impacted by the deepwater horizon oil spill. Appl. Environ. Microbiol. 77, 7962–7974. doi: 10.1128/AEM.05402-11
Leonardi, V., Giubilei, M. A., Federici, E., Spaccapelo, R., Sasek, V., Novotny, C., et al. (2008). Mobilizing agents enhance fungal degradation of polycyclic aromatic hydrocarbons and affect diversity of indigenous bacteria in soil. Biotechnol. Bioeng. 101, 273–285. doi: 10.1002/bit.21909
Liu, G. P.-W., Chang, T. C., Whang, L.-M., Kao, C.-H., Pan, P.-T., and Cheng, S.-S. (2011). Bioremediation of petroleum hydrocarbon contaminated soil: effects of strategies and microbial community shift. Int. Biodeter. Biodegr. 65, 1119–1127. doi: 10.1016/j.ibiod.2011.09.002
Liu, J. T., Hu, B., Gao, Y., Zhang, J. P., Jiao, B. H., Lu, X. L., et al. (2014). Bioactive tyrosine-derived cytochalasins from fungus Eutypella sp. D-1. Chem. Biodivers. 11, 800–806. doi: 10.1002/cbdv.201300218
MacNaughton, S. J., Stephen, J. R., Venosa, A. D., Davis, G. A., Chang, Y.-J., and White, D. C. (1999). Microbial population changes during bioremediation of an experimental oil spill. Appl. Environ. Microbiol. 65, 3566–3574.
Mahjoubi, M., Jaouani, A., Guesmi, A., Ben Amor, S., Jouini, A., Cherif, H., et al. (2013). Hydrocarbonoclastic bacteria isolated from petroleum contaminated sites in Tunisia: isolation, identification and characterization of the biotechnological potential. N. Biotechnol. 30, 723–733. doi: 10.1016/j.nbt.2013.03.004
Maiti, D., Chandra, K., Mondal, S., Ojha, A. K., Das, D., Roy, S. K., et al. (2008). Isolation and characterization of a heteroglycan from the fruits of Astraeus hygrometricus. Carbohyd. Res. 343, 817–824. doi: 10.1016/j.carres.2007.12.003
Marchand, C., St-Arnaud, M., Hogland, W., Bell, T. H., and Hijri, M. (2017). Petroleum biodegradation capacity of bacteria and fungi isolated from petroleum-contaminated soil. Int. Biodeter. Biodegr. 116, 48–57. doi: 10.1016/j.ibiod.2016.09.030
Masella, A. P., Bartram, A. K., Truszkowski, J. M., Brown, D. G., and Neufeld, J. D. (2012). PANDAseq: paired-end assembler for illumina sequences. BMC Bioinformatics 13:31. doi: 10.1186/1471-2105-13-31
McDonald, D., Price, M. N., Goodrich, J., Nawrocki, E. P., Desantis, T. Z., Probst, A., et al. (2012). An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 6, 610–618. doi: 10.1038/ismej.2011.139
Megharaj, M., Ramakrishnan, B., Venkateswarlu, K., Sethunathan, N., and Naidu, R. (2011). Bioremediation approaches for organic pollutants: a critical perspective. Environ. Int. 37, 1362–1375. doi: 10.1016/j.envint.2011.06.003
Micheluz, A., Sulyok, M., Manente, S., Krska, R., Varese, G., and Ravagnan, G. (2016). Fungal secondary metabolite analysis applied to cultural heritage: the case of a contaminated library in Venice. World Mycotoxin J. 9, 397–407. doi: 10.3920/WMJ2015.1958
Morais, D., Pylro, V., Clark, I. M., Hirsch, P. R., and Tótola, M. R. (2016). Responses of microbial community from tropical pristine coastal soil to crude oil contamination. PeerJ 4:e1733. doi: 10.7717/peerj.1733
Moubasher, H. A., Hegazy, A. K., Mohamed, N. H., Moustafa, Y. M., Kabiel, H. F., and Hamad, A. A. (2015). Phytoremediation of soils polluted with crude petroleum oil using Bassia scoparia and its associated rhizosphere microorganisms. Int. Biodeter. Biodegr. 98, 113–120. doi: 10.1016/j.ibiod.2014.11.019
Mulet, M., David, Z., Nogales, B., Bosch, R., Lalucat, J., and Garcia-Valdes, E. (2011). Pseudomonas diversity in crude-oil-contaminated intertidal sand samples obtained after the Prestige oil spill. Appl. Environ. Microbiol. 77, 1076–1085. doi: 10.1128/AEM.01741-10
Naranjo, L., Urbina, H., Sisto, A. D., and Leon, V. (2007). Isolation of autochthonous non-white rot fungi with potential for enzymatic upgrading of Venezuelan extra-heavy crude oil. Biocatal. Biotransfor. 25, 341–349. doi: 10.1080/10242420701379908
Nesci, A., Barros, G., Castillo, C., and Etcheverry, M. (2006). Soil fungal population in preharvest maize ecosystem in different tillage practices in Argentina. Soil Till. Res. 91, 143–149. doi: 10.1016/j.still.2005.11.014
Oehl, F., Sieverding, E., Ineichen, K., Ris, E. A., Boller, T., and Wiemken, A. (2005). Community structure of arbuscular mycorrhizal fungi at different soil depths in extensively and intensively managed agroecosystems. New Phytol. 165, 273–283. doi: 10.1111/j.1469-8137.2004.01235.x
Omrani, R., Spini, G., Puglisi, E., and Saidane, D. (2018). Modulation of microbial consortia enriched from different polluted environments during petroleum biodegradation. Biodegradation 29, 187–209. doi: 10.1007/s10532-018-9823-3
Phainuphong, P., Rukachaisirikul, V., Tadpetch, K., Sukpondma, Y., Saithong, S., Phongpaichit, S., et al. (2017). γ-butenolide and furanone derivatives from the soil-derived fungus Aspergillus sclerotiorum PSU-RSPG178. Phytochemistry 137, 165–173. doi: 10.1016/j.phytochem.2017.02.008
Ping, L., Luo, Y., Zhang, H., Li, Q., and Wu, L. (2007). Distribution of polycyclic aromatic hydrocarbons in thirty typical soil profiles in the Yangtze river delta region, east China. Environ. Poll. 147, 358–365. doi: 10.1016/j.envpol.2006.05.027
Poli, A., Lazzari, A., Prigione, V., Voyron, S., Spadaro, D., and Varese, G. C. (2016). Influence of plant genotype on the cultivable fungi associated to tomato rhizosphere and roots in differnet soils. Fungal Biol. U.K. 120, 862–872. doi: 10.1016/j.funbio.2016.03.008
Pruesse, E., Quast, C., Knittel, K., Fuchs, B. M., Ludwig, W., Peplies, J., et al. (2007). SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res. 35, 7188–7196. doi: 10.1093/nar/gkm864
Rajasekar, A., Babu, T. G., Pandian, S. T., Maruthamuthu, S., Palaniswamy, N., and Rajendran, A. (2007). Role of Serratia marcescens ACE2 on diesel degradation and its influence on corrosion. J. Ind. Microbiol. Biotechnol. 34, 589–598. doi: 10.1007/s10295-007-0225-5
Rebecchi, A., Pisacane, V., Miragoli, F., Polka, J., Falasconi, I., Morelli, L., et al. (2015). High-throughput assessment of bacterial ecology in hog, cow and ovine casings used in sausages production. Int. J. Food. Microbiol. 212, 49–59. doi: 10.1016/j.ijfoodmicro.2015.04.047
Rojas-Avelizapa, N. G., Cervantes-Gonzalez, E., Cruz-Camarillo, R., and Rojas-Avelizapa, L. I. (2002). Degradation of aromatic and asphaltenic fractions by Serratia liquefasciens and Bacillus sp. Bull. Environ. Contam. Toxicol. 69, 835–842. doi: 10.1007/s00128-002-0135-1
Roy, S., Hens, D., Biswas, D., Biswas, D., and Kumar, R. (2002). Survey of petroleum-degrading bacteria in coastal waters of Sunderban biosphere reserve. World J. Microb. Biotech. 18, 575–581. doi: 10.1023/A:1016362819746
Sadouk, Z., Tazerouti, A., and Hacene, H. (2008). Biodegradation of diesel oil and production of fatty acid esters by a newly isolated Pseudomonas citronellolis KHA. World J. Microb. Biotech. 25, 65–70. doi: 10.1007/s11274-008-9863-7
Schloss, P. D. (2010). The effects of alignment quality, distance calculation method, sequence filtering, and region on the analysis of 16S rRNA gene-based studies. PLoS Comput. Biolog. 6:e1000844. doi: 10.1371/journal.pcbi.1000844
Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541. doi: 10.1128/AEM.01541-09
Šimonovièová, A., Ferianc, P., Vojtková, H., Pangallo, D., Hanajík, P., Kraková, L., et al. (2017). Alkaline technosol contaminated by former mining activity and its culturable autochthonous microbiota. Chemosphere 171, 89–96. doi: 10.1016/j.chemosphere.2016.11.131
Siqueira, J. P. Z., Sutton, D. A., García, D., Gené, J., Thomson, P., Wiederhold, N., et al. (2016). Species diversity of Aspergillus section versicolores in clinical samples and antifungal susceptibility. Fungal Biol. U.K. 120, 1458–1467. doi: 10.1016/j.funbio.2016.02.006
Smith, E., Thavamani, P., Ramadass, K., Naidu, R., Srivastava, P., and Megharaj, M. (2015). Remediation trials for hydrocarbon-contaminated soils in arid environments: evaluation of bioslurry and biopiling techniques. Int. Biodeter.Biodegr. 101, 56–65. doi: 10.1016/j.ibiod.2015.03.029
Smits, T. H. M., Ford, R. M., Keel, C., Harms, H., Wick, L. Y., and Kohlmeier, S. (2005). Taking the fungal highway: mobilization of pollutant-degrading bacteria by fungi. Environ. Sci. Tech. 39, 4640–4646. doi: 10.1021/es047979z
Staley, J. T., and Konopka, A. (1985). Measurement of in situ activities of nonphotosynthetic microorganisms in aquatic and terrestrial habitats. Ann. Rev. Microbiol. 39, 321–346. doi: 10.1146/annurev.mi.39.100185.001541
Suja, F., Rahim, F., Taha, M. R., Hambali, N., Rizal Razali, M., Khalid, A., et al. (2014). Effects of local microbial bioaugmentation and biostimulation on the bioremediation of total petroleum hydrocarbons (TPH) in crude oil contaminated soil based on laboratory and field observations. Int. Biodeter.Biodegr. 90, 115–122. doi: 10.1016/j.ibiod.2014.03.006
Thompson, I. P., Van Der Gast, C. J., Ciric, L., and Singer, A. C. (2005). Bioaugmentation for bioremediation: the challenge of strain selection. Environ. Microbiol. 7, 909–915. doi: 10.1111/j.1462-2920.2005.00804.x
Tian, J., Dong, Q., Yu, C., Zhao, R., Wang, J., and Chen, L. (2016). Biodegradation of the organophosphate trichlorfon and its major degradation products by a novel Aspergillus sydowii PA F-2. J. Agr. Food Chem. 64, 4280–4287. doi: 10.1021/acs.jafc.6b00909
Varjani, S. J., Rana, D. P., Jain, A. K., Bateja, S., and Upasani, V. N. (2015). Synergistic ex-situ biodegradation of crude oil by halotolerant bacterial consortium of indigenous strains isolated from on shore sites of Gujarat, India. Int. Biodeter. Biodegr. 103, 116–124. doi: 10.1016/j.ibiod.2015.03.030
Varjani, S. J., and Upasani, V. N. (2016). Biodegradation of petroleum hydrocarbons by oleophilic strain of Pseudomonas aeruginosa NCIM 5514. Bioresour. Technol. 222, 195–201. doi: 10.1016/j.biortech.2016.10.006
Vasileiadis, S., Puglisi, E., Trevisa, M., Scheckel, K. G., Langdon, K., McLaughlin, M., et al. (2015). Changes in soil bacterial communities and diversity in response to long-term silver exposure. FEMS Microbiol. Ecol. 91:fiv114. doi: 10.1093/femsec/fiv114
Velmurugan, N., Kalpana, D., Cho, J.-Y., Lee, G.-H., Park, S.-H., and Lee, Y.-S. (2011). Phylogenetic analysis of culturable marine bacteria in sediments from South Korean Yellow Sea. Microbiology 80, 261–272. doi: 10.1134/S0026261711010188
Versalovic, J., Schneider, M., De Bruijn, F. J., and Lupski, J. R. (1994). Genomic fingerprinting of bacteria using repetitive sequence-based polymerase chain reaction. Method Mol. Cell. Bio. 5, 25–40.
Vinas, M., Sabate, J., Espuny, M. J., and Solanas, A. M. (2005). Bacterial community dynamics and polycyclic aromatic hydrocarbon degradation during bioremediation of heavily creosote-contaminated soil. Appl. Environ. Microbiol. 71, 7008–7018. doi: 10.1128/AEM.71.11.7008-7018.2005
Voyron, S., Ercole, E., Ghignone, S., Perotto, S., and Girlanda, M. (2017). Fine-scale spatial distribution of orchid mycorrhizal fungi in the soil of host-rich grasslands. New Phytol. 213, 1428–1439. doi: 10.1111/nph.14286
Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbio. 73, 5261–5267. doi: 10.1128/AEM.00062-07
Wang, X. B., Chi, C. Q., Nie, Y., Tang, Y. Q., Tan, Y., Wu, G., et al. (2011). Degradation of petroleum hydrocarbons (C6-C40) and crude oil by a novel Dietzia strain. Bioresour. Technol. 102, 7755–7761. doi: 10.1016/j.biortech.2011.06.009
Whyte, L. G., Bourbonniere, L., and Greer, C. W. (1997). Biodegradation of petroleum hydrocarbons by psychrotrophic Pseudomonas strains possessing both alkane (alk) and naphthalene (nah) catabolic pathways. Appl. Environ. Microbiol. 63, 3719–3723.
Wick, L. Y., Remer, R., Würz, B., Reichenbach, J., Braun, S., Schäfer, F., and Harms, H. (2007). Effect of fungal hyphae on the access of bacteria to phenanthrene in soil. Environ. Sci. Technol. 41, 500–505. doi: 10.1021/es061407s
Wongsa, P., Tanaka, M., Ueno, A., Hasanuzzaman, M., Yumoto, I., and Okuyama, H. (2004). Isolation and characterization of novel strains of Pseudomonas aeruginosa and Serratia marcescens possessing high efficiency to degrade gasoline, kerosene, diesel oil, and lubricating oil. Curr. Microbiol. 49, 415–422. doi: 10.1007/s00284-004-4347-y
Xia, W., Du, Z., Cui, Q., Dong, H., Wang, F., He, P., et al. (2014). Biosurfactant produced by novel Pseudomonas sp. WJ6 with biodegradation of n-alkanes and polycyclic aromatic hydrocarbons. J. Hazard. Mater. 276, 489–498. doi: 10.1016/j.jhazmat.2014.05.062
Yadav, T. C., Pal, R. R., Shastri, S., Jadeja, N. B., and Kapley, A. (2015). Comparative metagenomics demonstrating different degradative capacity of activated biomass treating hydrocarbon contaminated wastewater. Bioresour. Technol. 188, 24–32. doi: 10.1016/j.biortech.2015.01.141
Yang, S., Wen, X., Jin, H., and Wu, Q. (2012). Pyrosequencing investigation into the bacterial community in permafrost soils along the China-Russia crude oil pipeline (CRCOP). PLoS One 7:e52730. doi: 10.1371/journal.pone.0052730
Yergeau, E., Sanschagrin, S., Beaumier, D., and Greer, C. W. (2012). Metagenomic analysis of the bioremediation of diesel-contaminated Canadian high arctic soils. PLoS One 7:e30058. doi: 10.1371/journal.pone.0030058
Yu, K. S., Wong, A. H., Yau, K. W., Wong, Y. S., and Tam, N. F. (2005). Natural attenuation, biostimulation and bioaugmentation on biodegradation of polycyclic aromatic hydrocarbons (PAHs) in mangrove sediments. Mar. Pollut. Bull. 51, 1071–1077. doi: 10.1016/j.marpolbul.2005.06.006
Yuan, X., Zhang, X., Chen, X., Kong, D., Liu, X., and Shen, S. (2018). Synergistic degradation of crude oil by indigenous bacterial consortium and exogenous fungus Scedosporium boydii. Bioresour. Technol. 264, 190–197. doi: 10.1016/j.biortech.2018.05.072
Zafra, G., Absalón,ÁE., Cuevas, M. D. C., and Cortés-Espinosa, D. V. (2014). Isolation and selection of a highly tolerant microbial consortium with potential for PAH biodegradation from heavy crude oil-contaminated soils. Water Air Soil Poll. 225, 1826. doi: 10.1007/s11270-013-1826-4
Zhang, G., Guo, X., Zhu, Y., Liu, X., Han, Z., Sun, K., et al. (2018). The effects of different biochars on microbial quantity, microbial community shift, enzyme activity, and biodegradation of polycyclic aromatic hydrocarbons in soil. Geoderma 328, 100–108. doi: 10.1016/j.geoderma.2018.05.009
Zhou, Z.-F., Wang, M.-X., Zuo, X.-H., and Yao, Y.-H. (2017). Comparative investigation of bacterial, fungal, and archaeal community structures in soils in a typical oilfield in Jianghan, China. Arch. Environ. Con. Tox. 72, 65–77. doi: 10.1007/s00244-016-0333-1
Keywords: bioremediation, crude oil, soil contamination, enrichment culture, metagenomics, bacteria, fungi
Citation: Spini G, Spina F, Poli A, Blieux A-L, Regnier T, Gramellini C, Varese GC and Puglisi E (2018) Molecular and Microbiological Insights on the Enrichment Procedures for the Isolation of Petroleum Degrading Bacteria and Fungi. Front. Microbiol. 9:2543. doi: 10.3389/fmicb.2018.02543
Received: 31 July 2018; Accepted: 05 October 2018;
Published: 30 October 2018.
Edited by:Stéphane Pesce, Institut National de Recherche en Sciences et Technologies pour l’Environnement et l’Agriculture (IRSTEA), France
Reviewed by:Mariusz Cycoń, Medical University of Silesia, Poland
Maurizio Petruccioli, Università degli Studi della Tuscia, Italy
Copyright © 2018 Spini, Spina, Poli, Blieux, Regnier, Gramellini, Varese and Puglisi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Edoardo Puglisi, firstname.lastname@example.org