ORIGINAL RESEARCH article

Front. Microbiol., 20 June 2018

Sec. Microbiotechnology

Volume 9 - 2018 | https://doi.org/10.3389/fmicb.2018.01332

Functional Genetic Diversity and Culturability of Petroleum-Degrading Bacteria Isolated From Oil-Contaminated Soils

  • 1. Department of Energy & Resources Engineering, College of Engineering, Peking University, Beijing, China

  • 2. School of Environment, Tsinghua University, Beijing, China

Article metrics

View details

29

Citations

5k

Views

2k

Downloads

Abstract

In this study, we compared the culturability of aerobic bacteria isolated from long-term oil-contaminated soils via enrichment and direct-plating methods; bacteria were cultured at 30°C or ambient temperatures. Two soil samples were collected from two sites in the Shengli oilfield located in Dongying, China. One sample (S0) was close to the outlet of an oil-production water treatment plant, and the other sample (S1) was located 500 m downstream of the outlet. In total, 595 bacterial isolates belonging to 56 genera were isolated, distributed in Actinobacteria, Firmicutes, Bacterioidetes, and Proteobacteria. It was interesting that Actinobacteria and Firmicutes were not detected from the 16S rRNA gene clone library. The results suggested the activation of rare species during culture. Using the enrichment method, 239 isolates (31 genera) and 96 (22 genera) isolates were obtained at ambient temperatures and 30°C, respectively, from S0 soil. Using the direct-plating method, 97 isolates (15 genera) and 163 isolates (20 genera) were obtained at ambient temperatures and 30°C, respectively, from two soils. Of the 595 isolates, 244 isolates (41.7% of total isolates) could degrade n-hexadecane. A greater number of alkane-degraders was isolated at ambient temperatures using the enrichment method, suggesting that this method could significantly improve bacterial culturability. Interestingly, the proportion of alkane degrading isolates was lower in the isolates obtained using enrichment method than that obtained using direct-plating methods. Considering the greater species diversity of isolates obtained via the enrichment method, this technique could be used to increase the diversity of the microbial consortia. Furthermore, phenol hydroxylase genes (pheN), medium-chain alkane monooxygenases genes (alkB and CYP153A), and long-chain alkane monooxygenase gene (almA) were detected in 60 isolates (11 genotypes), 91 isolates (27 genotypes) and 93 isolates (24 genotypes), and 34 isolates (14 genotypes), respectively. This study could provide new insights into microbial resources from oil fields or other environments, and this information will be beneficial for bioremediation of petroleum contamination and for other industrial applications.

Introduction

Microbes not only have crucial roles in global carbon and nitrogen cycles, but also play an important role in industrial processes, such as bioremediation. Culture-independent methods can give us detailed information regarding the diversity of microbial communities across the world, and they can help to define functional genes via metagenomic, metatranscriptomic, and metaproteomic approaches. However, it is still important to study the functions and ecological roles of these microbes in their ecosystems; these studies are usually performed using pure cultures. Currently, the majority of microbes remain uncultivated (Hobbie et al., 1977; Roszak and Colwell, 1987). Therefore, improvements in the culturability of microbes are becoming increasingly important for understanding their functions in the environments, and their potential applications.

Bacterial alkane and aromatic compounds degradation are essential for bioremediation of petroleum-contaminated environments. More than 60 genera of aerobic bacteria and 5 genera of anaerobic bacteria that can degrade n-alkanes and more than 100 microbial genera than can degrade aromatic compounds have been reported. Alkane hydroxylases (AHs), which catalyze the hydroxylation of alkanes, are key enzymes in aerobic degradation of alkanes by bacteria. The integral-membrane alkane monooxygenase (AlkB)-related AHs and the cytochrome P450 CYP153 family are two common AHs in bacteria (Nie et al., 2014a). Furthermore, flavin-binding LC-alkane monooxygenase (almA) and thermophilic soluble LC-alkane monooxygenase (ladA) genes have also been found to be involved in hydroxylation of alkanes (Wentzel et al., 2007; Wang and Shao, 2012). Phenol hydroxylation is a vital step in biodegradation of phenolics, and is catalyzed by phenol hydroxylase. Among the phenol hydroxylase systems, two main classes have been reported, the “two-gene-component” and the “six-gene-component” systems (Powlowski and Shingler, 1994).

However, compared to known isolates, there are still many more unknown petroleum-degrading bacteria in the environment (Nie et al., 2014a). Culture-independent approaches using both, 16S rRNA genes (Jiao et al., 2016) and functional genes (Yang et al., 2014), have been used for isolation of diverse bacteria from petroleum contaminated soils. There are more studies focusing on the succession of bacterial community structures in oil-contaminated soils during bioremediation (Ai-Kindi and Abed, 2016; Leewis et al., 2016; Pacwa-Plociniczak et al., 2016; Wu et al., 2016, 2017). It has been shown that bacterial compositions differ significantly between culture-dependent and -independent methods (Stefani et al., 2015). After years of enrichment, the diverse microbes in petroleum contaminated soils provide an ideal species pool for isolating petroleum-degrading bacteria (Bhattacharya et al., 2003; M'rassi et al., 2015; El Mahdi et al., 2016). However, compared to the microbial diversity in petroleum contaminated soils estimated using culture-independent studies, the number of isolates obtained from these environments is still limited and the functions of majority of the microbes in such environments are still unknown.

The limitations of culturing such microbes may be due to failure to mimic essential aspects of bacterial environments, such as nutrition, environmental conditions, and microbial interactions (Zengler et al., 2002). Previous studies have tried to improve the culturability of microbes isolated from the environment by varying incubating conditions, such as using different medium (Gartner et al., 2011; Sun et al., 2014b), temperature (Weichart et al., 1992; Junge et al., 2002), pressure (Gartner et al., 2011), as well as adding antioxidants (Muresu et al., 2013). Few studies have focused on the effects of oscillations in ambient-temperature on bacterial culturability. Interestingly, bacterial strains also synchronize their circadian rhythms to day/night cycles (Huang et al., 1990), and oscillations in ambient-temperature in nature may influence the growth of microbes. In this study in order to investigate the functions of microbes in petroleum-contaminated soils, we improved the culturability of microbes by simulating ambient temperatures in nature, and characterized new petroleum-degrading bacteria from two different petroleum contaminated soils. The distribution of genes associated with petroleum degradation was also investigated.

Materials and methods

Sample collection

Two soil samples were collected from two sites in the Shengli oilfield located in Dongying, Shandong Province, China. One sample (S0) was close to the outlet of an oil-production water treatment plant, and the other sample (S1) was located 500 m downstream of the outlet. The soil at these sites had been contaminated by crude oil for more than 10 years. Soil samples were transported to the laboratory at 4°C immediately following collection. The characteristics of the soil, analyzed by standard methods, are listed in Table S1.

The oil-production water used for cultivating the bacterial strains was sampled from the Gudao wastewater treatment plant in the same oil field. The physical-chemical properties of the oil-production water were determined using standard methods. The characteristics of the oil-production water were as follows: pH 8.5, 22.82 mg total nitrogen l−1, 0.18 mg total phosphorous l−1, 1,640 mg CODCr l−1. The heavy crude oil used as carbon and energy source was collected directly from the outlet of the production well.

16S rRNA gene clone library analysis

Extraction of DNA from oil-contaminated soil and amplification of 16S rRNA genes were carried out as previously described (Yu et al., 2011). PCR products were purified using the PCR-product purification kit (BioTech, China), and were ligated into pGEM T-easy vectors (Promega, USA) according to manufacturer's instruction. The ligated products were transformed into Escherichia coli DH5α, and ampicillin resistant colonies were picked up. The inserted fragments from each colony were amplified using the primers SP6 (ATT TAG GTG ACA CTA TAG) and T7 (TAA TAC GAC TCA CTA TAG GG). Amplified products were separated into different types using restriction fragment length polymorphism (RFLP) analysis with RsaI and HhaI in order to pick up clones with different 16S rRNA gene insertion. In total, 55 positive clones were picked up from the clone library; these clones were subdivided into 21 different types according to their RFLP results. Thirty clones were selected for DNA sequencing.

Bacteria isolation

Aerobic bacteria were isolated from soil sample S0 using two different strategies. In the first method, soil sample S0 was immediately transferred to the appropriate media for enrichment after being transported to the laboratory. Three kinds of media were used for enrichment. The G3 mineral salt medium consisted of 5.0 g NaCl l−1, 1.0 g K2HPO4 l−1, 1.0 g NH4H2PO4 l−1, 1.0 g NH4SO4 l−1, 0.2 g MgSO4·7H2O l−1, and 3.0 g KNO3 l−1 in ddH2O, adjusted to pH 8.0–8.5. The G2 medium was made using oil-production water supplemented with 1.0 g yeast extract l−1. The G5 mineral salt medium had the same mineral salt concentration and pH as the G3 medium, but the salts were dissolved in 0.9 l oil-production water and 0.1 l ddH2O. One gram of sampled soil was added into each medium supplemented with 1% sterilized heavy crude oil. Cultures were then incubated at ambient, or at a constant temperature of 30°C; cultures were shaken at 150 rpm. Cultures were passaged using fresh media supplemented with 2–10% crude oil every 20 days. After three passages, the enriched cultures were serially diluted and plated onto agars. Agar plates were incubated at the same temperatures as the enrichment temperature. After a 14–day incubation period, all colonies on agar plates were collected for further investigation.

For the second strategy, bacteria were directly isolated from soil samples without enrichment. Soil samples were diluted with the 10-fold dilution method, and the dilutions were plated on oil-production water agar (PW agar: 1000 ml of oil-production water, 15 g of agar) and mineral salt agar (MSM agar: 0.5 g NaCl l−1, 1 g NH4H2PO4 l−1, 1 g (NH4)2SO4 l−1, 0.2 g MgSO4·7H2O l−1, 3 g KNO3 l−1, 1 g K2HPO4 l−1, 15 g agar l−1). The agar plates were incubated at either ambient or constant temperature of 30°C in the dark. After a 1-week incubation period, all colonies present on the plates were collected for further investigation. Bacteria were also directly isolated from soil samples using the same method.

The prefixes of the isolate names were proposed according to sampling sites, the isolation medium, and the isolation temperature (Table S2). Purified isolates were primarily identified via 16S rRNA gene analysis, which was performed according to previously described protocol (Wang et al., 2007b). The 16S rRNA gene sequence obtained was analyzed using the BLAST tool in NCBI (http://blast.ncbi.nlm.nih.gov/Blast.cgi) and Classifier in RDP (http://rdp.cme.msu.edu/classifier/classifier.jsp). Reference sequences were retrieved from GenBank. After multiple sequence alignments using the CLUSTAL_X tool and manual correction, a phylogenetic tree based on the 16S rRNA gene was constructed using the neighbor-joining method (Saitou and Nei, 1987) contained in the MEGA 6.0 software (Tamura et al., 2013). Tree topology was evaluated with bootstrap analysis based on 1,000 resampling replicates.

Amplification of alkB, CYP153A, alma, and pheN genes

Genomic DNA was extracted from the isolates as previously described (Wang et al., 2007b). The phenol hydroxylase genes pheN as well as alkane monooxygenase genes including alkB, CYP153A, and almA were amplified and sequenced as previously described (Sun et al., 2014b). Representative isolates of all phylotypes were selected for target gene detection. The PCR products were purified and directly sequenced using amplification primers (Sun et al., 2014b). The obtained sequences were searched against the NR database using BLAST (http://blast.ncbi.nlm.nih.gov/Blast.cgi), and the nearest reference sequences were retrieved. The derived amino acid sequences were aligned by CLUSTAL_X and manual correction. The phylogenetic tree was constructed after alignment using the neighbor-joining method (Saitou and Nei, 1987) in the MEGA 6.0 software (Tamura et al., 2013). Tree topology was evaluated with bootstrap analysis based on 1,000 resampling replicates.

Detection of phenol- and alkane-degrading abilities

The isolates were first cultured in artificial sea water medium (ASW: 1.0 g yeast extract L−1, 5.0 g tryptone l−1, 4.0 g Na2SO4 l−1, 0.68 g KCl l−1, 0.1 g KBr L−1, 0.025 g H3BO3 l−1, 5.4 g MgCl2·H2O l−1, 1.5 g CaCl2·H2O l−1, 0.024 g SrCl2·6H2O l−1, 0.2 g NaHCO3 l−1, 0.2 g Na2HPO4 l−1, 0.5 g NH4Cl l−1, 0.002 g NaF l−1, 24 g NaCl l−1, 15 g agar l−1; pH 7.0-8.0) at 25°C with shaking at 150 rpm for 5 days. Cells were then harvested by centrifugation (2,000 g at 4°C for 10 min). Cell pellets were washed twice using sterile saline, and were suspended in an equal volume of saline. Then, 30 μl cells were inoculated into 3 ml MSM supplemented with 100 mg phenol l−1 (final concentration) or 100 mg hexadecane l−1 as sole carbon sources. Cell cultures were incubated at 30°C with shaking at 150 rpm in the dark, and cell growth was measured. After 5 days, cell cultures were sampled, and the concentrations of residual phenol and hexadecane were detected, as previously described (Sun et al., 2011).

Sequence accession numbers

The 16S rRNA gene sequences of the clone library from contaminated soils were deposited into GenBank under the accession numbers KC763859-KC763888. The 16S rRNA gene sequences of the isolated strains were under the accession numbers KP711432-KP711625; the phenol hydroxylase genes and the flavin-binding monooxygenase gene almA were under the numbers KP711626-KP711653; the alkB and CYP153A genes were deposited into GenBank under the numbers KT160363-KT160415.

Results

Isolation of aerobic bacterial strains from soils under ambient and constant temperatures

In this work, a total of 595 bacterial isolates were isolated from the soils, including 335 isolates via S0 soil enrichments and 260 isolates via direct isolation from both soil S0 and S1. Among these, 336 and 259 bacterial species were obtained under ambient and constant temperatures, respectively.

The isolates were distributed into four phyla (Actinobacteria, Firmicutes, Bacterioidetes, and Proteobacteria), and belonged to 56 genera; they were classified into 149 phylotypes. Proteobacteria was the most dominant phylum in the isolates. At the genus level, Pseudomonas, Marinobacter, and Bacillus were the dominant genera (Figure 1, Figure S1, Table 2). According to colony morphology and phylogenetic analysis, 335 isolates obtained from “S0” soil using the enrichment method were affiliated into 39 genera and were classified into 109 phylotypes. It was notable that enrichment and isolation at ambient temperatures enhanced bacterial culturability. (i) The number of isolates obtained when incubated at ambient temperatures was much higher as compared with that when incubated at a constant temperature of 30°C. Ambient temperature enrichment yielded 239 isolates that were distributed into 31 genera; 96 isolates distributed into 22 genera were obtained from enrichment under constant temperature. (ii) Isolates cultured in ambient temperatures exhibited higher α-diversity as compared with those cultured at 30°C. The Shannon-Wiener index values for bacteria cultured at ambient temperature were all higher than 2.0, while those cultured at a constant temperature were < 2.0 (Table 2); (iii) Higher number of unique genera (only obtained by one cultivating method) were obtained by incubating at ambient temperatures. We collected 17 unique genera, including Achromobacter, Agromyces, Bordetella, Brachybacterium, Dietzia, Enterobacter, Erythrobacter, Hydrogenophaga, Hyphomonas, Idiomarina, Kocuria, Mycobacterium, Providencia, Pusillimonas, Rheinheimera, Staphylococcus, and Vibrio, by incubating at ambient temperatures. Only seven unique genera were obtained by incubating at a constant temperature of 30°C. In addition, 15 genera were isolated under both culture conditions (Figure 1 and Figure S2). (iv) More novel or potential novel strains with lower 16S rRNA gene identities (< 98.6%) were obtained at ambient temperatures as compared with those obtained at 30°C (Table 1); two strains obtained at ambient temperatures have been published as type strains for novel species: Negadavirga shengliensis SLG210A2-5 (Hu et al., 2015) and Glycocaulis albus SLG210-30A1T (Lv et al., 2014). These results suggested that ambient temperature cultures improve the culturability of bacteria from environmental samples.

Figure 1

Figure 1

Composition fractions of isolates cultured in different conditions and from different soil samples on a (A) genus and (B) phylum levels. Bacterial composition obtained from of (C) cultured strains and (D) via the clone library method.

Table 1

S0S1
G2G3G5PWMSMPWMSM
ABABABABABABAB
Isolates103178842483773831343792282
Phylotypes357201022156221112217838
Genus176138151341288105519
n-Hexadecane degrader431297162452216132231528
n-Hexadecane degrader (%)41.75.933.016.733.364.971.457.951.638.259.533.368.234.1
Novel or potential strains71226101001007
Simpson index0.790.790.870.670.850.840.710.760.670.700.760.810.800.91
Shannon-Wiener index2.081.562.181.432.211.951.151.821.441.501.750.731.452.58

Basic information of strains isolated from different enrichment conditions.

Furthermore, the effects of ambient/constant temperature differ between various taxa. For example, more number of Gammaproteobacteria and Betaproteobacteria were obtained by cultivating under the ambient temperatures: 8 Betaproteobacteria phylotypes (35 isolates, 6 genera) were isolated at ambient temperature, while only one phylotype (including 23 isolates) were obtained at 30°C. Similarly, a higher number of Gammaproteobacteria isolates were obtained at ambient temperature as compared with those obtained at a constant temperature. However, at the order level, Pseudoxanthomonas belonging to Xanthomonadales was only obtained from cultures incubated at 30°C (Figure S1).

A total of 110 and 150 isolates were obtained from S0 and S1 soils, respectively, using the direct-plating (without enrichment) method. These 260 isolates (27 genera) were affiliated with Actinobacteria, Firmicutes, and Proteobacteria, and were classified into 32 genera, and 45 phylotypes. In contrast to the enrichment method, more strains were isolated at constant temperature (163 isolates, 20 genera), rather than at ambient temperatures (97 isolates, 15 genera) using this method. (i) The α-diversities (Simpson index) of isolates cultured under ambient temperatures were lower as compared with those cultured under 30°C constant temperature (Table 1). (ii) More number of novel strains were obtained by culturing at a constant temperature (8 novel species) than were obtained by culturing at ambient temperatures (1 novel species) (Table 1). The strains directly isolated from the soil have been identified as novel taxonomic units, and include Salinarimonas ramus SL014B-41A4T (Cai et al., 2011b), Tessaracoccus oleiagri SL014B-20A1T (Cai et al., 2011b), Halomonas shengliensis SL014B-85T (Wang et al., 2007a), Halomonas gudaonensis SL014B-69T (Wang et al., 2007b), Marinobacter gudaonensis SL014B-61AT (Gu et al., 2007), Rubrimonas shengliensis SL014B-28A2T (Cai et al., 2011a), Polymorphum gilvum SL003B-26A1T (Cai et al., 2011a) and Nitratireductor shengliensis SL014B-25AT (Pan et al., 2014), and Marinobacter shengliensis SL013A-34A2T (Luo et al., 2015).

Isolation conditions affect the proportion of petroleum degraders

Hexadecane was able to be used as the sole carbon and energy source for growth of 244 isolates (Tables 1, 2). Using the enrichment method, 88 (36.8% of isolates) and 32 (33.0% of isolates) hexadecane-degrading isolates were obtained by incubating at ambient temperatures and a constant temperature, respectively. Specifically, using the G2 medium for enriching strains, 43 and one hexadecane-degrading isolates were obtained by the ambient and constant temperature method; the ratio of the number of hexadecane-degrading isolates obtained under ambient temperatures to that obtained under a constant temperature method were 29:7 and 16:24 for G3 and G5 media, respectively. Using the direct-plating method, 21 and 35 hexadecane-degrading isolates were obtained from S0 soil at ambient and constant temperatures, respectively; 37 and 31 hexadecane-degrading isolates were isolated from S1 soil at ambient and constant temperatures, respectively. It is worth noting that the number of hexadecane-degrading isolates obtained by culturing at ambient temperature were slightly higher than those obtained by culturing at a constant temperature. In addition, the direct plating method led to a higher ratio of hexadecane degraders as compared with that obtained via the enrichment method using G2 and G3 as enrichment media. Specifically, some isolates, such as Pseudomonas spp. SLG510A3-14 and SL014A-33A1-2, Enterobacter sp. SLG510A3-6, Acinetobacter spp. SLG310A2-9 and SL003A-18A1, Nitratireductor sp. SLG510B6-18, Bacillus sp. SLG210A2-18-1, Rhizobium sp. SL014B-37A2, Labrenzia sp. SL014B-4A2, Paracoccus sp. SL013B-8A1, Planococcus sp. SL013A-41A-2, Halomonas sp. SL014B-63A1-2, and Marinobacter sp. SL013A-25A demonstrated robust growth in MSM supplemented with hexadecane as the sole carbon source. Among them, nine strains were obtained from ambient-temperature cultures, and only six were obtained from constant-temperature cultures; five strains were obtained via the enrichment method, while eight strains were obtained via the direct-plating method. Furthermore, alkane-degrading ability was first detected in the Enterobacter and Nitratireductor strains.

Table 2

StrainClosest type strainsAmountGenBank No.Identities (%)BiodegradationGenes
Phen-HexapheNalkBP450almA
SLG510A3-4A2Pseudomonas songnenensis1NR_14829599.9+
SLG210A2-34-3Pseudomonas songnenensis1NR_14829599.9+
SLG310A2-2Pseudomonas songnenensis2NR_14829599.9+
SLG310A2-40Pseudomonas songnenensis13LT62979799.0+
SLG510A3-14Pseudomonas aeruginosa2KT825518100+++++
SLG310A2-6A1-2Pseudomonas balearica3NR_02597299.7++
SLG210B1-9A1Marinobacter alkaliphilus1NR_112223100++
SLG510A10-4Marinobacter shengliensis3KP711478.1100++++
SLG210A2-10B2-2Providencia burhodogranariea2NR_10491499.4++
SLG510A3-6Enterobacter asburiae1CP01186399.6++++++
SLG310A2-9Acinetobacter venetianus16NR_042049100++++++
SLG510A10-22B1Acinetobacter venetianus1NR_042049100++++++
SLG510A3-11A1Acinetobacter venetianus2NR_042049100++++++
SLG510A6-3Acinetobacter venetianus2NR_042049100++++++
SLG510B10-21Acinetobacter venetianus2NR_042049100++++++
SLG310B1-12Pseudoxanthomonas mexicana2NR_113973100+
SLG310B2-40A1Pseudoxanthomonas mexicana1NR_113973100+
SLG210A2-35-2Hydrogenophaga palleronii1NR_11413299.1++
SLG210A2-28A1Bordetella bronchiseptica2NR_11362898.1+
SLG310A2-11A1Achromobacter pulmonis7NR_11764499.8++
SLG510A3-26Achromobacter pulmonis2NR_11764499.8++
SLG210A2-36-1Tistrella mobilis1NR_114036100+
SLG210A3-2-1Tistrella mobilis1NR_114036100+
SLG510A2-12-2Tistrella mobilis1NR_114036100+
SLG510A10-23Tistrella mobilis1NR_114036100+
SLG510B2-10Tistrella mobilis8NR_114036100+
SLG510B10-1Tistrella mobilis5NR_114036100+
SLG310A2-23Sphingopyxis chilensis1NR_02463199.6+
SLG310B1-4-2Sphingopyxis chilensis4NR_02463199.6+
SLG310B2-17Sphingopyxis chilensis4NR_02463199.6+
SLG310A2-18B3Erythrobacter citreus2NR_02874199.4+
SLG210A2-1-1Hyphomonas polymorpha3KF86314498.6+
SLG310B2-7A2-2Brevundimonas diminuta1NR_11360299.9++
SLG510B6-18Nitratireductor aquibiodomus1NR_02526299.6++++++
SLG310A2-16Nitratireductor indicus4NR_117518100+
SLG510B6-1Nitratireductor indicus2NR_117518100+
SLG310A2-20Nitratireductor soli3KP63957099.9+
SLG310B1-11Nitratireductor indicus1NR_11751897.8+
SLG310B2-11Nitratireductor indicus1NR_11751897.8+
SLG310B1-40A2Rhizobium halotolerans1NR_12563299.9+
SLG310B2-4-1Rhizobium halotolerans1NR_12563299.9+
SLG510A10-14B1Rhizobium halotolerans1NR_12563299.9+
SLG510A3-33Rhizobium halotolerans1NR_12563299.9+
SLG510A10-12Rhizobium helanshanense1NR_13301999.6++
SLG210A2-38-1Agrobacterium fabrum1NR_074266100++
SLG210B1-12Agrobacterium fabrum1NR_074266100++
SLG510A10-7B2Agrobacterium fabrum1NR_074266100++
SLG510B10-28Agrobacterium fabrum2NR_074266100++
SLG210A6-6Exiguobacterium aestuarii1NR_04300599.5+
SLG210A2-29-1Bacillus tianshenii11NR_13370499.8+++
SLG210A3-15-4Bacillus tianshenii4NR_13370499.8+++
SLG210A2-14C1Bacillus muralis5NR_04208399.9+
SLG210A3-9-1Bacillus muralis1NR_04208399.9+
SLG210A2-25-1Bacillus vietnamensis1NR_11399599.1+
SLG210A2-2Bacillus idriensis1NR_02902299.7+
SLG210A2-11Bacillus halosaccharovorans1NR_10911699.1+
SLG210A2-18-1Bacillus litoralis3NR_04301598.5++++
SLG210A3-8Bacillus litoralis3NR_04301598.5++++
SLG510B10-40B1-2Microbacterium thalassium2NR_04248199.4++--
SLG310A2-4-1Agromyces italicus1NR_04302699.5+
SLG510A2-27-2Kocuria turfanensis2DQ531634100+
SLG510B2-13A2Arthrobacter tumbae1NR_04207899.7++
SLG210A6-7Arthrobacter protophormiae2NR_026195100+
SLG210A2-38-2Arthrobacter protophormiae1NR_026195100+
SLG510B10-2-2Micrococcus terreus1NR_116649100+
SL003A-1A2Rhizobium rosettiformans1NR_116445100+
SL003B-10A2Rhizobium selenitireducens6NR_04421699.0+
SL004A-8A1Rhizobium rosettiformans14NR_116445100+
SL014A-30ARhizobium rosettiformans1NR_116445100+
SL014B-16A1-1Rhizobium selenitireducens1NR_04421699.0+
SL014B-37A2Rhizobium tarimense3NR_11785098.8++++
SL014B-87A3Stappia indica1NR_11643199.8++
SL003B-19A3Polymorphum gilvum4NR_074240100+
SL014B-4A2Labrenzia aggregata1NR_11386199.7+++
SL014B-40A2Stappia meyerae1NR_11594998.8++
SL004A-4AStappia meyerae1NR_11594998.8++----
SL014B-56A3-1Paracoccus homiensis4NR_04373399.5+
SL013B-8A1Paracoccus marcusii1NR_044922100++++
SL014B-90A-2Salinarimonas rosea3NR_11648798.6++
SL014B-17A3Porphyrobacter donghaensis1NR_02581699.8+
SL013B-2A1-3Porphyrobacter donghaensis1NR_02581699.8+
SL013A-50APorphyrobacter colymbi2NR_11432899.4+
SL014A-10A1Porphyrobacter colymbi9NR_11432899.4+
SL013A-40A1Arthrobacter tumbae1NR_04207899.9+
SL013A-41A-2Planococcus antarcticus3CP01653499.0+++++
SL003B-7A-2Planococcus antarcticus1CP01653499.0+++++
SL014B-56A1-1Planococcus antarcticus1CP01653499.0+++++
SL014B-4A1Lysinibacillus alkalisoli1KX25875799.6+
SL003B-20A2Lysinibacillus fusiformis1NR_11256999.9+
SL004B-24ALysinibacillus fusiformis1NR_11256999.9+
SL004B-4ABacillus megaterium5KJ47672199.7+
SL014B-7A1-2Bacillus megaterium1KJ47672199.7+
SL003B-7B1-1Bacillus subterraneus2NR_10474999.4+
SL013B-14A3Bacillus stratosphericus1NR_118441100++
SL003B-30B1Bacillus stratosphericus1NR_118441100++
SL014B-20A2Bacillus stratosphericus1NR_118441100++
SL003B-19A2Bacillus mojavensis5NR_118290100+
SL003B-11B2-2Bacillus aquimaris1NR_02524199.5-
SL013A-42A1Thauera selenatis1NR_02647498.5+
SL013A-36A-1Thauera chlorobenzoica2CP01883999.0++
SL004B-31A1Vibrio fluvialis2CP01403599.9+
SL014B-21B1-2Pseudomonas stutzeri1CP002881100++++
SL004A-44APseudomonas stutzeri6CP002881100+
SL004B-21B1-1Pseudomonas stutzeri6CP002881100+
SL003B-18B2-1Pseudomonas stutzeri2CP002881100+
SL003A-2A1-2Pseudomonas xanthomarina2NR_04104499.1+
SL003B-7A-3Pseudomonas japonica1KT82551999.9++
SL004B-28A-3Pseudomonas japonica1KT82551999.9+++
SL014B-42A3Pseudomonas japonica2KT82551999.9++++
SL003A-18A1Acinetobacter indicus1NR_117784100+++++
SL014B-63A1-2Halomonas ventosae1NR_04281299.7+++++
SL014A-11A2Halomonas ventosae3NR_04281299.9++
SL013A-32A1Halomonas ventosae1NR_04281299.9+++
SL004A-30AHalomonas ventosae1NR_04281299.9+++
SL014A-23AHalomonas campaniensis1NR_04215799.6++
SL013A-18AHalomonas meridiana1NR_04206699.3+
SL013A-7A1Halomonas zhaodongensis2NR_12561299.7++
SL004B-23A3-1Halomonas zhaodongensis1NR_12561299.7++
SL013B-12A1-3Halomonas zhaodongensis1NR_12561299.7++
SL014B-11A2*Marinobacter gudaonensis9DQ414419100+++
SL014A-2A1Marinobacter flavimaris2NR_02579999.5+++
SL013A-35AMarinobacter flavimaris2NR_02579999.5+++
SL014B-80A2Marinobacter hydrocarbonoclasticus2NR_074619100++++
SL014A-10A1-2Marinobacter hydrocarbonoclasticus2NR_074619100+++
SL013A-7A2Marinobacter hydrocarbonoclasticus2NR_074619100+++
SL004A-39AMarinobacter hydrocarbonoclasticus1NR_07461999.9++
SL014B-45A-2Marinobacter shengliensis1KP711478100+++++
SL013B-15A5Marinobacter shengliensis3KP711478100++++
SL013A-34A-2*Marinobacter shengliensis1KP711478100+++++
SL003A-19A-1Marinobacter shengliensis1KP711478100+++
SL013A-38AMarinobacter alkaliphilus1NR_11222399.7++
SL013A-24AMarinobacter shengliensis1KP711478100++++
SL013A-25AMarinobacter shengliensis10KP711478100+++++++
SL004B-10A-2Donghicola xiamenensis2NR_04356598.1+
SL014B-22A1Labrenzia alba1NR_04237899.0+
SL014B-55A2Rhizobium tarimense2NR_117850100+
SL004B-4ABacillus megaterium1KJ47672199.7+
SL014B-14A2Microbacterium saccharophilum1NR_11434299.5+
SL003B-19B2-1Leucobacter komagatae1NR_114929100+
SL003B-1A8Pseudoxanthomonas mexicana1NR_113973100+
SL004B-16A2-3Pseudomonas stutzeri1CP00288199.7++
SL013B-15A2Marinobacter shengliensis1KP71147899.9++
SL014B-81A3-1Marinobacter segnicrescens1NR_044124100+
SL014B-23A1-1Marinobacter flavimaris3NR_02579999.4++
SL014B-22A2-2Ornithinimicrobium pekingense1NR_04388699.7+++
SL014A-33A1-2Pseudomonas graminis1NR_02639598.9+++++++

Degrading characteristics of isolated bacterial strains and the gene amplification of phenol-degrading gene and alkane-degrading genes.

a

growth ability of the strains with n-hexadecane in MSM medium, +++++, ++++, +++, ++, + and –, indicating the growth capability from strong to weak with n-hexadecane as sole carbon and energy source, measured by optical density at 600 nm. +++++, growth (OD600 > 0.5) after a 7-day incubation at 30°C; ++++, growth (0.5 > OD600 > 0.4) after a 7-day incubation at 30°C; +++, growth (0.4 > OD600 > 0.3); ++, growth (0.3 > OD600 > 0.2) after a 7-day incubation at 30°C; +, growth (0.1 < OD600 < 0.2) after a 7-day incubation at 25°C; –, no growth.

*

strain is the type strain.

Results indicated that 27 isolates were capable of degrading phenols. Among them, 18 isolates were obtained from ambient temperature cultures, and 9 isolates were obtained from constant temperature cultures. All phenol-degrading bacterial strains were assigned to Proteobacteria: one Brevundimonas and two Stappia isolates belonging to Alphaproteobacteria, three Thauera and nine Achromobacter isolates belonging to Betaproteobacteria, as well as five Marinobacter and seven Pseudomonas isolates belonging to Gammaproteobacteria (Table 2). These were some differences in the taxonomy of phenol-degrading strains isolated from oil-production water; all phenol-degrading strains belonged to Betaproteobacteria (Sun et al., 2014b).

Distribution of pheN, alkB, CYP153A, and almA genes in isolates

Phenol hydroxylase, which catalyzes the initial step in most phenol catabolic pathway, is usually used to evaluate the phenol-degrading ability of bacterial strains (Qian et al., 1997; Sun et al., 2012, 2014a). The pheN gene was detected in 60 isolates, including 35 isolates (6 phylotypes; 10.7% of total isolates) in ambient-temperature cultures and 25 isolates (8 phylotypes; 9.6% of total isolates from constant temperature) in constant-temperature cultures. Conversely, more number of strains harboring the pheN gene were isolated by the direct-plating method (34 isolates, 11 phylotypes; 13.6% of the total strains) as compared with that isolated by the enrichment method (26 isolates, 4 phylotypes; 7.7% of the total strains). These pheN genes were categorized into two clusters according to their phylogenetic distance (Figure 2). Cluster I contained 52 sequences (4 genotypes) from 52 isolates (10 phylotypes) belonging to Marinobacter, Pseudomonas, Achromobacter, and Halomonas. Among these, Achromobacter belonged to Betaproteobacteria, while the others belonged to Gammaproteobacteria. Gene sequences in Cluster I were closely related to pheN from the Pseudomonas strains. Cluster II contained 8 sequences (3 genotypes) from 4 phylotypes (8 isolates) belonging to Pseudomonas, Paracoccus, Porphyrobacter, and Brevundimonas; this cluster was closely related to those belonging to Alphaproteobacteria and Actinobacteria. We noted that different cultivating methods may result in different genes; strains harboring the Cluster II pheN genes in the present study were only obtained from constant temperature cultures, while all strains obtained from ambient temperature cultures and a few strains obtained from constant temperature cultures contained Cluster I pheN genes. Additionally, although the pheN gene was detected in Marinobacter sp. SLG510A10-4, Paracoccus sp. SL013B-8A1, Porphyrobacter sp. SL014B-17A3, and Halomonas sp. SL014B-63A1-2 and SL013A-18A, phenol-degrading abilities were not detected in these strains (Figure 2).

Figure 2

Figure 2

Phylogenetic trees were constructed based on partial amino sequences of (A) phenol hydroxylase and (B) 16S rRNA gene sequence of bacterial strains obtained from the oil-contaminated soils. Tree topology was evaluated by bootstrap analysis based on 1,000 resampling replicates. Bootstrap values (%) are indicated at the nodes; only values > 50 are shown. The strains isolated at ambient temperature were indicated by bolded red letters; isolates cultured at 30°C were indicated by bold letters in blue; numbers in square brackets are the number of isolates within this phylotype. Scale bars represent 0.05 substitutions per site.

Terminal hydroxylation of alkanes yields the corresponding primary alcohols, and it is a crucial step that significantly enhances the hydrophilicity of substrates (Nie et al., 2011, 2014a,b). In bacteria, this reaction is catalyzed by different enzymes, including those encoded by almA, CYP153A, and alkB/alkM. A total of 27 unique CYP153A sequences were derived from 91 bacterial isolates that span over 13 genera; 24 unique alkB/alkM sequences were derived from 93 isolates in nine genera, and 14 almA sequences were obtained from 34 isolates in two genera.

The sequences of CYP153A, encoding a medium-chain-length alkane hydroxylase obtained from 91 isolates, were categorized into five clusters (Figure 3). As the biggest cluster, Cluster I contained 43 sequences (6 genotypes) from 36 isolates (11 phylotypes) belonging to Marinobacter and Rhizobium. Among them, three and two CYP153A sequences are present in Marinobacter spp. SL013B-15A2 and SL013A-34A-2, respectively. Cluster II included a unique sequence from 24 Acinetobacter isolates, which is closely related to Acinetobacter sp. NBRC100985 (ZP09221884). Cluster III consisted of two unique sequences from 11 Marinobacter isolates, which were closely related to Alvanivorax sp. DG881 (EDX90518) and Marinobacter hydrocarbonoclasticus ATCC 49840 (YP005431650) (Grimaud et al., 2012). Cluster IV was composed of P450 CYP153A genes from Amycolicicoccus subflavus DQS3-9A1T (Nie et al., 2013). This cluster consisted of only two unique sequences from two distinct genera: Nitratireductor belonging to Alphaproteobacteria, and Pseudomonas belonging to Gammaproteobacteria. Cluster V consisted of seven unique sequences from 25 isolates, which belonged to Hydrogenophaga (Betaproteobacteria), Erythrobacter, Hyphomonas, Sphingopyxis, and Porphyrobacter; these genes clustered with the CYP153 genes from strains belonging to Alphaproteobacteria. The ratio of the number of CYP153A-harboring isolates in ambient-temperature cultures to that in constant-temperature cultures was 24:13 (8:3 phylotypes), 22:2 (2:1), 0:7 (0:2), 1:1 (1:1), and 7:18 (4:3) in Cluster I, II, III, IV, and V, respectively, suggesting that ambient temperature may benefit isolation of stains containing the CYP153A gene. Furthermore, the method of isolation also somewhat affected the type of strains isolated. The ratio of the number of CYP153A-harboring isolates obtained by the enrichment method to that obtained by direct-plating method was 33:4 (8:4), 22:2 (2:1), 1:1 (1:1), 0:7 (0:2), and 14:11 (5:2) in cluster I, II, III, IV, and group II, respectively. Additionally, horizontal gene transfer was found in CYP153A genes. For example, CYP153A sequences from the Rhizobium sp. SLG510A10-12 and Marinobacter sp. SLG210B1-9A1 strains had identical CYP153A genes, and clustered together in the phylogenetic tree; this was inconsistent with the 16S rRNA gene tree (Figure 1). Furthermore, strains harboring Cluster I, Cluster IV, as well as Cluster V CYP153A genes did not grow in MSM supplemented with hexadecane as the sole carbon source. Interestingly, most of strains harboring the Cluster II and IV CYP153A gene were obtained by the enrichment method (Figure 3).

Figure 3

Figure 3

Phylogenetic trees were constructed based on the partial amino sequences of the (A) P450 (CYP153A) and (B) 16S rRNA gene sequences of bacterial strains which were obtained from the oil-contaminated soils. The strains isolated at ambient temperature were indicated by bolded red letters; isolates cultured at 30°C were indicated by bold letters in blue; numbers in square brackets are the number of isolates within this phylotype. The scale bar represents 0.1 substitutions per site.

In total, 13 unique medium-length-chain hydroxylase alkB/alkM genes, an integral-membrane non-heme diiron monooxygenase, were obtained from 93 isolates and categorized into four clusters (Figure 4). As the biggest cluster, cluster I comprised of six unique sequences clustered with alkB/alkM genes from Marinobacter hydrocarbonoclasticus VT8 (Singer et al., 2011). Among these sequences, only one sequence was obtained from Rhizobium sp. SLG510A10-12, the others were obtained from Marinobacter strains. Genes in cluster II, which clustered with Acinetobacter sp. ADP1, consisted only of one unique sequence from three distinct genera: Acinetobacter, Bacillus, and Microbacterium, which belonged to Gammaproteobacteria, Firmicutes, and Actinobacteria, respectively. Cluster III contained two unique sequences from 24 isolates belonging to Nitratireductor (Alphaproteobacteria), Microbacterium (Actinobacteria), and Pseudomonas (Gammaproteobacteria) strains. This cluster was usually detected in strains belonging to Actinobacteria, in the Dietzia genus (Sun et al., 2014b). Cluster IV had only a single unique sequence from three Marinobacter isolates obtained at a constant temperature, and it showed the highest identity to the alkB gene in Pseudomonas sp. GPo1 (van Beilen et al., 2001). Cluster V consisted of three unique sequences from 17 isolates belonging to two Gammaproteobacteria genera: Enterobacter and Pseudomonas. Similar to the case of CYP153A gene, cultivation mode also somewhat affected the clusters of alkB/alkM genes. (i) The ratio of the number of alkB/alkM-harboring isolates by direct-plating method to that by enrichment methods was 42:5 (11:3), 0:24 (0:4), 1:3 (1:3), 3:0 (1:0), and 0:17 (0:3) in Cluster I, II, III, IV, and V, respectively. The direct-plating method resulted in more isolates (strains) harboring sequences in Cluster I and IV, while the enrichment method resulted in strains harboring alkB/alkM in Cluster II and IV. (ii) The ratio of the number of strains obtained under ambient temperatures to that obtained under a constant temperature was 27:20 (9:5), 22:5 (3:2), 1:3 (1:2), 0:3 (0:1), and 14:6 (2:2) in Cluster I, II, III, IV, and V, respectively, suggesting that ambient temperatures promote isolation of strains harboring cluster II and Cluster V alkB/alkM genes, while a constant temperature results in higher number of isolates containing alkB/alkM genes in Cluster III and IV. For the first time, an alkane monooxygenase gene and hexadecane-degrading ability were detected in Nitratireductor species. (Figure 4) More interestingly, among the strains harboring Cluster I and V alkB/alkM genes, only strains Marinobacter sp. SL013A-34A-2 harboring alkB gene of Cluster I and Enterobacter sp. SLG510A3-6 could utilize alkanes as the sole carbon source for growth; no obvious growth was detected in other strains. In contrast, four of the six strains harboring Cluster II genes, and all three strains harboring Cluster III alkB/alkM demonstrated significant growth under the same culture conditions. It was interesting to note that most of the strains harboring Cluster II and III alkB/alkM were obtained via the enrichment method (Figure 4).

Figure 4

Figure 4

Phylogenetic tree was constructed based on the partial amino sequences of (A) alkane monooxygenase (AlkB) and (B) 16S rRNA gene sequences of bacterial strains obtained from oil-contaminated soils. The scale bar represents 0.1 substitutions per site. Acinetobacter venetianus 6A2 has the same 16S rRNA as Acinetobacter venetianus RAG-1 (Throne-Holst et al., 2006).

Another putative monooxygenase of the flavin-binding family AlmA catalyzes the metabolism of long-chain n-alkanes. In the present study, a total of 12 unique almA sequences were detected in 36 isolates belonging to two genera: Marinobacter and Rhizobium. These genes were all closely related to those from Marinobacter strains regardless of their distinct taxonomic position. We found only 5.5% (19 isolates) of ambient-temperature isolates were contained the almA genes, while 6.1% of constant-temperature (17 isolates) isolates harbored the same genes (Figure 5).

Figure 5

Figure 5

Phylogenetic trees were constructed based on the partial amino sequence of (A) alkane monooxygenase (AlmA) and (B) the 16S rRNA gene sequences of bacterial strains obtained from the oil-contaminated soils. The scale bar represents 0.5 and 0.05 substitutions per site.

It was observed that the Marinobacter strains can express multiple alkane hydroxylase genes. For example, three kinds of alkane monooxygenase genes including alkB, CYP153A, alma, and phenol hydroxylase genes were simultaneously detected in Marinobacter spp. SLG510A10-4, SL014B-80A2, and SL014B-45A-2. Marinobacter gudaonensis SL014B-11A2T, Marinobacter spp. SL013A-35A, SL004A-39A, SL013A-34A-2, SL013A-38A, SL013A-24A, SL013A-25A, and SL013B-15A2 all expressed alkB, CYP153A, and almA genes. Among them, strain SL013A-34A-2 had two copies of the CYP153A gene. In addition, Marinobacter sp. SLG210B1-9A1 harbored both alkB and CYP153A genes, and CYP153A and almA genes were simultaneously detected in Marinobacter sp. SL014B-23A1-1. Furthermore, three copies of the CYP153A gene were detected in Marinobacter sp. SL013-15A2. However, most of these strains did not grow in the MSM supplemented with hexadecane.

Discussion

The clone library result showed that the bacterial community in oil-contaminated soils is comprised of species belonging to Proteobacteria (69.1%), Cyanobacteria (1.8%), Bacteroidetes (1.8%), and several unknown bacterial species (27.0%); no strains belonging to Actinobacteria and Firmicutes were detected. Proteobacteria in oil-contaminated soil consisted of Gammaproteobacteria, Epsilonproteobacteria, Alphaproteobacteria, Betaproteobacteria, and Deltaproteobacteria. Among them, the proportion of Gammaproteobacteria and Epsilonproteobacteria were much greater as compared with that of the three other subdivisions of Proteobacteria (Figure 1D). It is surprising that some popular oil-degrading species in the environment such as Pseudomonas and Dietzia strains (Tang et al., 2012; Gharibzahedi et al., 2014) were not detected with the clone library method. Isolation of bacterial strains from the bacterial community in the soil yielded very different result as compared with that from the cultivated method. For example, Actinobacteria and Firmicutes strains were isolated from the two soil samples, but were not detected in the colony library. It has been reported that the rare biosphere taxa in a given environmental have important ecological roles and are reservoirs of genetic and functional diversity (Lynch and Neufeld, 2015). They may become dominant under optimal conditions. Our study suggested the activation and rapid growth of rare species during culture. Previous studies also demonstrated that the culture-dependent and culture-independent methods provided two extremely divergent microbial communities: nearly 40% isolated strains from oil-contaminated soil were unable to be detected by high-throughput sequencing method (M'rassi et al., 2015). It was found that most of the shared operational taxonomic units (OTUs) between the two isolation methods belonged to the Sphingobium and Sphingomonas (Sphingomonadales), and Pseudomonas (Pseudomonadales) and Arthrobacter (Actinomycetales) genera.

Along with application of culture-independent methods in microbial ecology, it is recognized that most of the strains in the environment are uncultured, and are difficult to culture with the usual methods; these are known as uncultured strains (Zengler et al., 2002). It is believed that cultural media, temperature, pressure, antioxidant, and carbon and nitrogen sources could affect the isolation of microbes from various environments (Gartner et al., 2011; Sun et al., 2014b). Results from the present study also suggested that different isolating conditions, including carbon source, incubation temperature, and concentration of carbon affect the diversities of isolates. Among these factors, temperature is one the most important factors affecting the culturability of bacteria from the environments. According to their optimal temperature, bacteria could be subdivided into psychrophilic, mesophilic, and thermophilic bacteria. In earlier studies, constant temperatures were used to incubate and isolate bacterial strains, such as at low temperature (Song et al., 2009), extreme high temperature (Gonzalez et al., 1996), and mesophilic temperature (Sun et al., 2014b); very few studies used ambient-temperatures. The present study suggested that using ambient-temperature during the enrichment step could significantly enhance bacterial culturability. This may be attributed to the circadian rhythm. Previous studies have demonstrated that prokaryotes also follow circadian rhythms, similar to those experience by eukaryotes (Huang et al., 1990). Furthermore, this rhythm is temperature-dependent (Sweeney and Borgese, 1989). The natural environment shows oscillating temperatures that are correlated with day/night cycle, such as the temperature of the top layer of soils. Because of this endogenous “circadian clocks,” indigenous microbes may have adapted to oscillating temperatures along with day/night cycles, thus forming the special bacterial circadian rhythm found in the natural environment. As a result, incubation at constant temperature may disturb the inherent rhythm of bacterial strains in the environment, which leads to reduced culturability. Additionally, each strain has an inherent optimal temperature for its metabolism. A constant incubation temperature may only be able to provide the optimal temperature for a small number of the strains, while the other strains are kept in the dormant phase or grow at much slower rates due to the unchanging temperature. Therefore, ambient temperatures that change with time may provide a wider growth window for more bacterial strains; however, they may grow slowly under this condition. The above two reasons may explain why using ambient temperatures could enhance culturability of strains isolated from the environment.

In earlier studies, bacterial strains harboring the alkB/alkM gene were found to affected by the isolation medium (Sun et al., 2014b). In the present study, isolation of bacteria containing various phenol and alkane hydroxylase genes was also affected by the isolating condition, namely the cultivating mode and temperature. The different kinds of alkane or phenol hydroxylase genes may bestow the host with different alkane- or phenol-degrading abilities (Nie et al., 2014b; Sun et al., 2014a,b). In fact, enrichment using specific carbon sources is a selective process of survival of the fittest, during which a large number of strains without alkane-degrading abilities or with inefficient alkane-degrading abilities are eliminated. The different alkane-hydroxylase catalysis ability could ultimately determine the survival of bacterial strains. More investigations need to be made on the mechanisms that underlie the relationship between the temperature and bacterial genes clusters.

Results from the present study also suggested that the direct-plating method could result in higher percentage of hexadecane- and phenol-degrading isolates as compared with that obtained via enrichment using mineral salt medium and oil-production water. The lower number of alkane- or phenol-degrading isolates in the enrichment method may be due to the abundance of rare microorganism that could not degrade crude oil. In this system, metabolites from crude oil degraders provided other bacteria with essential nutrients. In addition, crude oil contains diverse carbon sources, which can be easily used by many bacterial strains for growth. Because of the existence of the easily degradable carbon source, strains that are unable to degrade AHs also could grow vigorously in the enrichment culture. Therefore, the process of enrichment with crude oil may decline the percentage of n-hexadecane degraders.

It has been shown that pheN, CYP153A, alkB, and almA genes correlated with the phylogeny of their host. For example, CYP153A, alkB, and almA genes from Marinobacter strains were often tightly clustered (Figures 35). In contrast, many diverse strains also contained near similar genes. For example, the alkB genes from Acinetobacter sp. SLG310A2-9, Rhizobium sp. SLG510A10-12, and Microbacterium sp. SLG510B10-40B1-2 were found to be highly similar to those from Acinetobacter venetianus (AEW67043) (Mara et al., 2012). This may be attributed to horizontal gene transfer (HGT), which has been shown to be ubiquitous in the environment (de la Cruz and Davies, 2000). More interestingly, the less-oil-contaminated soil (S1) also possessed numerous alkane-degrading bacterial strains; some bacterial strains exhibited three different types of alkane-degrading genes. This suggested that alkane-degrading strains are ubiquitously distributed in the natural environments regardless of oil contamination.

It is worth noting that the occurrence of the multiple alkB, CYP153A, and almA alkane hydroxylases in one bacterium is quite frequent in Marinobacter strains. The coexistence of multiple alkane hydroxylases was also observed in Dietzia sp. DQ12-45-1b (Nie et al., 2011), Amycolicicoccus subflavus DQS3-9A1T (AEF42720) (Nie et al., 2013), Acinetobacter sp. ADP1 (Barbe et al., 2004), and many Alcanivorax isolates (Wang et al., 2010; Wang and Shao, 2012). These three alkane hydroxylases demonstrated complementary substrate ranges: CYP153A usually catalyzed the hydroxylation of the alkane with C5-C10, AlkB catalyzed C5-C16, while AlmA catalyzed C32-C36. The coexistence of alkane hydroxylase genes undoubtedly broadened the substrate range, thereby enhancing the adaptive ability of the host. While many strains showed tremendous growth in MSM containing n-hexadecane as the sole carbon and energy source, none of three abovementioned genes were amplified from them, including Rhizobium sp. SL014B-37A2, Stappia sp. SL014B-87A3, Labrenzia sp. SL014B-4A2, Paracoccus sp. SL013B-8A1, Planococcus sp. SL013A-41A-2, Pseudomonas spp. SLG510A3-14, Enterobacter sp. SLG510A3-6, and Bacillus spp. SLG210A2-18-1 and SLG210A2-29A1. It is possible that these isolates harbored other alkane hydroxylases, which should be investigated in the future.

In summary, under different incubation conditions, 595 bacterial strains were isolated from oil contaminated soils. These strains were affiliated into 56 genera belonging to Actinobacteria, Proteobacteria, Firmicutes, and Bacteroidetes. Using the enrichment method, 239 isolates distributed into 31 genera, and 96 isolates distributed into 22 genera were obtained by culturing at ambient temperatures and a constant temperature, respectively, from a single soil sample. With the direct-plating method, 97 isolates affiliated with 15 genera and 163 isolates affiliated with 20 genera were obtained by culturing under ambient temperatures and a constant temperature, respectively, from two soil samples. Furthermore, phenol hydroxylase genes, short- and medium-chain alkane monooxygenase AlkB and CYP153A, and long-chain alkane monooxygenase AlmA were detected in numerous bacterial strains, including some genera for the first time. We found that several strains could utilize phenol and n-alkanes as the sole carbon source for growth, but no known related genes were found in these strains, suggesting the existence of novel genes in these strains that need to be further studied in the future.

Statements

Author contributions

J-QS and X-LW designed the research. J-QS and LX performed the research. X-YL analyzed the 16S rRNA sequences. G-FZ and HC participated in isolation of bacterial strains. J-QS and YN wrote the paper.

Funding

This research was supported by the National Natural Science Foundation of China (31200100 and 31225001), the National Basic Research Program of China (973 Program, 2014CB846002), and the Inner Mongolia Science & Technology Plan (20130101).

Conflict of interest

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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.01332/full#supplementary-material

References

  • 1

    Ai-KindiS.AbedR. M. (2016). Effect of biostimulation using sewage sludge, soybean meal, and wheat straw on oil degradation and bacterial community composition in a contaminated desert soil. Front. Microbiol.7:240. 10.3389/fmicb.2016.00240

  • 2

    BarbeV.VallenetD.FonknechtenN.KreimeyerA.OztasS.LabarreL.et al. (2004). Unique features revealed by the genome sequence of Acinetobacter sp. ADP1, a versatile and naturally transformation competent bacterium. Nucleic Acids Res.32, 57665779. 10.1093/nar/gkh910

  • 3

    BhattacharyaD.SarmaP. M.KrishnanS.MishraS.LalB. (2003). Evaluation of genetic diversity among Pseudomonas citronellolis strains isolated from oily sludge-contaminated sites. Appl. Environ. Microbiol.69, 14351441. 10.1128/aem.69.3.1435-1441.2003

  • 4

    CaiM.WangL.CaiH.LiY.TangY. Q.WuX. L. (2011a). Rubrimonas shengliensis sp. nov. and Polymorphum gilvum gen. nov., sp. nov., novel members of Alphaproteobacteria from crude oil contaminated saline soil. Syst. Appl. Microbiol.34, 321327. 10.1016/j.syapm.2011.03.004

  • 5

    CaiM.WangL.CaiH.LiY.WangY. N.TangY. Q.et al. (2011b). Salinarimonas ramus sp nov and Tessaracoccus oleiagri sp nov., isolated from a crude oil-contaminated saline soil. Int. J. Syst. Evol. Microbiol.61, 17671775. 10.1099/ijs.0.025932-0

  • 6

    de la CruzF.DaviesJ. (2000). Horizontal gene transfer and the origin of species: lessons from bacteria. Trends Microbiol.8, 128133. 10.1016/S0966-842x(00)01703-0

  • 7

    El MahdiA. M.AzizH. A.Abu AmrS. S.El-GendyN. S.NassarH. N. (2016). Isolation and characterization of Pseudomonas sp NAF1 and its application in biodegradation of crude oil. Environ. Earth Sci.75:380. 10.1007/s12665-016-5296-z

  • 8

    GärtnerA.BlümelM.WieseJ.ImhoffJ. F. (2011). Isolation and characterisation of bacteria from the Eastern Mediterranean deep sea. Antonie Van Leeuwenhoek100, 421435. 10.1007/s10482-011-9599-5

  • 9

    GharibzahediS. M. T.RazaviS. H.MousaviS. M. (2014). Characterization of bacteria of the genus Dietzia: an updated review. Ann. Microbiol.64, 111. 10.1007/s13213-013-0603-3

  • 10

    GonzalezJ. M.KatoC.HorikoshiK. (1996). Culturability and survival of an extreme thermophile from deep-sea hydrothermal vents. Arch. Microbiol.166, 6467. 10.1007/s002030050356

  • 11

    GrimaudR.GhiglioneJ. F.CagnonC.LaugaB.VaysseP. J.Rodriguez-BlancoA.et al. (2012). Genome sequence of the marine bacterium Marinobacter hydrocarbonoclasticus sp17, which forms biofilms on hydrophobic organic compounds. J. Bacteriol.194, 35393540. 10.1128/Jb.00500-12

  • 12

    GuJ.CaiH.YuS. L.QuR.YinB.GuoY. F.et al. (2007). Marinobacter gudaonensis sp nov., isolated from an oil-polluted saline soil in a Chinese oilfield. Int. J. Syst. Evol. Microbiol.57, 250254. 10.1099/ijs.0.64522-0

  • 13

    HobbieJ. E.DaleyR. J.JasperS. (1977). Use of nuclepore filters for counting bacteria by fluorescence microscopy. Appl. Environ. Microbiol.33, 12251228.

  • 14

    HuB.YangQ.CaiM.TangY. Q.ZhaoG. F.WuX. L. (2015). Negadavirga shengliensis gen. nov., sp nov., a novel member of the family Cyclobacteriaceae isolated from oil-contaminated saline soil. Antonie Van Leeuwenhoek107, 663673. 10.1007/s10482-014-0361-7

  • 15

    HuangT. C.TuJ.ChowT. J.ChenT. H. (1990). Circadian-rhythm of the prokaryote Synechococcus sp. RF-1. Plant Physiol.92, 531533. 10.1104/pp.92.2.531

  • 16

    JiaoS.LiuZ. S.LinY. B.YangJ.ChenW. M.WeiG. H. (2016). Bacterial communities in oil contaminated soils: biogeography and co-occurrence patterns. Soil Biol. Biochem.98, 6473. 10.1016/j.soilbio.2016.04.005

  • 17

    JungeK.ImhoffF.StaleyT.DemingJ. W. (2002). Phylogenetic diversity of numerically important arctic sea-ice bacteria cultured at subzero temperature. Microb. Ecol.43, 315328. 10.1007/s00248-001-1026-4

  • 18

    LeewisM. C.UhlikO.FraraccioS.McFarlinK.KottaraA.GloverC.et al. (2016). Differential impacts of willow and mineral fertilizer on bacterial communities and biodegradation in diesel fuel oil-contaminated soil. Front. Microbiol.7:837. 10.3389/fmicb.2016.00837

  • 19

    LuoY. J.XieB. S.LvX. L.CaiM.WangY. N.CuiH. L.et al. (2015). Marinobacter shengliensis sp nov., a moderately halophilic bacterium isolated from oil-contaminated saline soil. Antonie Van Leeuwenhoek107, 10851094. 10.1007/s10482-015-0401-y

  • 20

    LvX. L.XieB. S.CaiM.GengS.TangY. Q.WangY. N.et al. (2014). Glycocaulis albus sp nov., a moderately halophilic dimorphic prosthecate bacterium isolated from petroleum-contaminated saline soil. Int. J. Syst. Evol. Microbiol.64, 31813187. 10.1099/ijs.0.063537-0

  • 21

    LynchM. D.NeufeldJ. D. (2015). Ecology and exploration of the rare biosphere. Nat. Rev. Microbiol.13, 217229. 10.1038/nrmicro3400

  • 22

    MaraK.DecorosiF.VitiC.GiovannettiL.PapaleoM. C.MaidaI.et al. (2012). Molecular and phenotypic characterization of Acinetobacter strains able to degrade diesel fuel. Res. Microbiol.163, 161172. 10.1016/j.resmic.2011.12.002

  • 23

    Guermouche M'rassiA.BensalahF.GuryJ.DuranR. (2015). Isolation and characterization of different bacterial strains for bioremediation of n-alkanes and polycyclic aromatic hydrocarbons. Environ. Sci. Pollut. Res.22, 1533215346. 10.1007/s11356-015-4343-8

  • 24

    MuresuR.TondelloA.PoloneE.SulasL.BaldanB.SquartiniA. (2013). Antioxidant treatments counteract the non-culturability of bacterial endophytes isolated from legume nodules. Arch. Microbiol.195, 385391. 10.1007/s00203-013-0886-6

  • 25

    NieY.ChiC. Q.FangH.LiangJ. L.LuS. L.LaiG. L.et al. (2014a). Diverse alkane hydroxylase genes in microorganisms and environments. Sci. Rep.4:4968. 10.1038/srep04968

  • 26

    NieY.FangH.LiY.ChiC. Q.TangY. Q.WuX. L. (2013). The genome of the moderate halophile Amycolicicoccus subflavus DQS3-9A1t reveals four alkane hydroxylation systems and provides some clues on the genetic basis for its adaptation to a petroleum environment. PLoS ONE8:e70986. 10.1371/journal.pone.0070986

  • 27

    NieY.LiangJ.FangH.TangY. Q.WuX. L. (2011). Two novel alkane hydroxylase-rubredoxin fusion genes isolated from a Dietzia bacterium and the functions of fused rubredoxin domains in long-chain n-alkane degradation. Appl. Environ. Microbiol.77, 72797288. 10.1128/Aem.00203-11

  • 28

    NieY.LiangJ. L.FangH.TangY. Q.WuX. L. (2014b). Characterization of a CYP153 alkane hydroxylase gene in a Gram-positive Dietzia sp. DQ12-45-1b and its “team role” with alkW1 in alkane degradation. Appl. Microb. Biotechnol.98, 163173. 10.1007/s00253-013-4821-1

  • 29

    Pacwa-PlociniczakM.PlazaG. A.Piotrowska-SegetZ. (2016). Monitoring the changes in a bacterial community in petroleum-polluted soil bioaugmented with hydrocarbon-degrading strains. Appl. Soil Ecol.105, 7685. 10.1016/j.apsoil.2016.04.005

  • 30

    PanX. C.GengS.MeiR.WangY. N.CaiH.LiuX. Y.et al. (2014). Nitratireductor shengliensis sp nov., isolated from an oil-polluted saline Soil. Curr. Microbiol.69, 561566. 10.1007/s00284-014-0624-6

  • 31

    PowlowskiJ.ShinglerV. (1994). Genetics and biochemistry of phenol degradation by Pseudomonas sp. CF600. Biodegradation5, 219236.

  • 32

    QianH.EdlundU.PowlowskiJ.ShinglerV.SethsonI. (1997). Solution structure of phenol hydroxylase protein component P2 determined by NMR spectroscopy. Biochemistry36, 495504. 10.1021/bi9619233

  • 33

    RoszakD. B.ColwellR. R. (1987). Metabolic activity of bacterial cells enumerated by direct viable count. Appl. Environ. Microbiol.53, 28892893.

  • 34

    SaitouN.NeiM. (1987). The neighbor-joining method - a new method for reconstructing phylogenetic trees. Mol. Biol. Evol.4, 406425.

  • 35

    SingerE.WebbE. A.NelsonW. C.HeidelbergJ. F.IvanovaN.PatiA.et al. (2011). Genomic potential of marinobacter aquaeolei, a biogeochemical “opportunitroph”. Appl. Environ. Microbiol.77, 27632771. 10.1128/aem.01866-10

  • 36

    SongJ.OhH. M.ChoJ. C. (2009). Improved culturability of SAR11 strains in dilution-to-extinction culturing from the East Sea, West Pacific Ocean. FEMS Microbiol. Lett.295, 141147. 10.1111/j.1574-6968.2009.01623.x

  • 37

    StefaniF. O.BellT. H.MarchandC.de la ProvidenciaI. E.El YassimiA.St-ArnaudM.et al. (2015). Culture-dependent and -independent methods capture different microbial community fractions in hydrocarbon-contaminated soils. PLoS ONE10:e0128272. 10.1371/journal.pone.0128272

  • 38

    SunJ. Q.XuL.TangY. Q.ChenF. M.WuX. L. (2012). Simultaneous degradation of phenol and n-hexadecane by Acinetobacter strains. Bioresour. Technol.123, 664668. 10.1016/j.biortech.2012.06.072

  • 39

    SunJ. Q.XuL.TangY. Q.ChenF. M.LiuW. Q.WuX. L. (2011). Degradation of pyridine by one Rhodococcus strain in the presence of chromium (VI) or phenol. J. Hazard. Mater.191, 6268. 10.1016/j.jhazmat.2011.04.034

  • 40

    SunJ. Q.XuL.TangY. Q.ChenF. M.ZhaoJ. J.WuX. L. (2014a). Bacterial pyridine hydroxylation is ubiquitous in environment. Appl. Microbiol. Biotechnol.98, 455464. 10.1007/s00253-013-4818-9

  • 41

    SunJ. Q.XuL.ZhangZ.LiY.TangY. Q.WuX. L. (2014b). Diverse bacteria isolated from microtherm oil-production water. Antonie Van Leeuwenhoek105, 401411. 10.1007/s10482-013-0088-x

  • 42

    SweeneyB. M.BorgeseM. B. (1989). A circadian-rhythm in cell-division in a prokaryote, the Cyanobacterium synechococcus Wh7803. J. Phycol.25, 183186. 10.1111/j.0022-3646.1989.00183.x

  • 43

    TamuraK.StecherG.PetersonD.FilipskiA.KumarS. (2013). MEGA6: molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol.30, 27252729. 10.1093/molbev/mst197

  • 44

    TangY. Q.LiY.ZhaoJ. Y.ChiC. Q.HuangL. X.DongH. P.et al. (2012). Microbial communities in long-term, water-flooded petroleum reservoirs with different in situ temperatures in the huabei oilfield, China. PLoS ONE7:e33535. 10.1371/journal.pone.0033535

  • 45

    Throne-HolstM.MarkussenS.WinnbergA.EllingsenT. E.KotlarH. K.ZotchevS. B. (2006). Utilization of n-alkanes by a newly isolated strain of Acinetobacter venetianus: the role of two AlkB-type alkane hydroxylases. Appl. Microbiol. Biotechnol.72, 353360. 10.1007/s00253-005-0262-9

  • 46

    van BeilenJ. B.PankeS.LucchiniS.FranchiniA. G.RothlisbergerM.WitholtB. (2001). Analysis of Pseudomonas putida alkane-degradation gene clusters and flanking insertion sequences: evolution and regulation of the alk genes. Microbiology147, 16211630. 10.1099/00221287-147-6-1621

  • 47

    WangL.WangW.LaiQ.ShaoZ. (2010). Gene diversity of CYP153A and AlkB alkane hydroxylases in oil-degrading bacteria isolated from the Atlantic Ocean. Environ. Microbiol.12, 12301242. 10.1111/j.1462-2920.2010.02165.x

  • 48

    WangW.ShaoZ. (2012). Diversity of flavin-binding monooxygenase genes (almA) in marine bacteria capable of degradation long-chain alkanes. FEMS Microbiol. Ecol.80, 523533. 10.1111/j.1574-6941.2012.01322.x

  • 49

    WangY. N.CaiH.ChiC. Q.LuA. H.LinX. G.JiangZ. F.et al. (2007a). Halomonas shengliensis sp nov., a moderately halophilic, denitrifying, crude-oil-utilizing bacterium. Int. J. Syst. Evol. Microbiol.57, 12221226. 10.1099/ijs.0.64973-0

  • 50

    WangY. N.CaiH.YuS. L.WangZ. Y.LiuJ.WuX. L. (2007b). Halomonas gudaonensis sp nov., isolated from a saline soil contaminated by crude oil. Int. J. Syst. Evol. Microbiol.57, 911915. 10.1099/ijs.0.64826-0

  • 51

    WeichartD.OliverJ. D.KjellebergS. (1992). Low-temperature induced nonculturability and killing of Vibrio vulnificus. FEMS Microbiol. Lett.100, 205210. 10.1016/0378-1097(92)90210-F

  • 52

    WentzelA.EllingsenT. E.KotlarH. K.ZotchevS. B.Throne-HolstM. (2007). Bacterial metabolism of long-chain n-alkanes. Appl. Microbiol. Biotechnol.76, 12091221. 10.1007/s00253-007-1119-1

  • 53

    WuM. L.DickW. A.LiW.WangX. C.YangQ.WangT. T.et al. (2016). Bioaugmentation and biostimulation of hydrocarbon degradation and the microbial community in a petroleum-contaminated soil. Int. Biodeterior. Biodegradation107, 158164. 10.1016/j.ibiod.2015.11.019

  • 54

    WuM. L.YeX. Q.ChenK. L.LiW.YuanJ.JiangX. (2017). Bacterial community shift and hydrocarbon transformation during bioremediation of short-term petroleum-contaminated soil. Environ. Pollut.223, 657664. 10.1016/j.envpol.2017.01.079

  • 55

    YangY.WangJ.LiaoJ.XieS.HuangY. (2014). Distribution of naphthalene dioxygenase genes in crude oil-contaminated soils. Microb. Ecol.68, 785793. 10.1007/s00248-014-0457-7

  • 56

    YuS.LiS.TangY.WuX. (2011). Succession of bacterial community along with the removal of heavy crude oil pollutants by multiple biostimulation treatments in the Yellow River Delta, China. J. Environ. Sci.23, 15331543. 10.1016/S1001-0742(10)60585-2

  • 57

    ZenglerK.ToledoG.RappeM.ElkinsJ.MathurE. J.ShortJ. M.et al. (2002). Cultivating the uncultured. Proc. Natl. Acad. Sci. U.S.A.99, 1568115686. 10.1073/pnas.252630999

Summary

Keywords

ambient temperature, bacteria, bioremediation, culturability, petroleum degradation

Citation

Sun J-Q, Xu L, Liu X-Y, Zhao G-F, Cai H, Nie Y and Wu X-L (2018) Functional Genetic Diversity and Culturability of Petroleum-Degrading Bacteria Isolated From Oil-Contaminated Soils. Front. Microbiol. 9:1332. doi: 10.3389/fmicb.2018.01332

Received

06 April 2018

Accepted

31 May 2018

Published

20 June 2018

Volume

9 - 2018

Edited by

Rajesh K. Sani, South Dakota School of Mines and Technology, United States

Reviewed by

Jiandong Jiang, Nanjing Agricultural University, China; Binbin Liu, Chinese Academy of Sciences, China

Updates

Copyright

*Correspondence: Xiao-Lei Wu

This article was submitted to Microbiotechnology, Ecotoxicology and Bioremediation, a section of the journal Frontiers in Microbiology

†These authors have contributed equally to this work.

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics