Abstract
Microbial consortia consisting of cooperational strains exhibit biodegradation performance superior to that of single microbial strains and improved remediation efficiency by relieving the environmental stress. Tetrahydrofuran (THF), a universal solvent widely used in chemical and pharmaceutical synthesis, significantly affects the environment. As a refractory pollutant, THF can be degraded by some microbial strains under suitable conditions. There are often a variety of stresses, especially pH stress, that inhibit the THF-degradation efficiency of microbial consortia. Therefore, it is necessary to study the molecular mechanisms of microbial cooperational degradation of THF. In this study, under conditions of low pH (initial pH = 7.0) stress, a synergistic promotion of the THF degradation capability of the strain Rhodococcus ruber YYL was found in the presence of a non-THF degrading strain Bacillus cereus MLY1. Metatranscriptome analysis revealed that the low pH stress induced the strain YYL to up-regulate the genes involved in anti-oxidation, mutation, steroid and bile acid metabolism, and translation, while simultaneously down-regulating the genes involved in ATP production. In the co-culture system, strain MLY1 provides fatty acids, ATP, and amino acids for strain YYL in response to low pH stress during THF degradation. In return, YYL shares the metabolic intermediates of THF with MLY1 as carbon sources. This study provides the preliminary mechanism to understand how microbial consortia improve the degradation efficiency of refractory furan pollutants under environmental stress conditions.
Introduction
Bioremediation is an eco-friendly waste management technique that uses naturally-occurring organisms to break down hazardous substances into less toxic or non-toxic substances. Microbial bioremediation has been the most significant remediation method for the mineralization of pollutants in contaminated environments (Vilchez-Vargas et al., 2010; ). For this reason, many degrading microorganisms have been isolated from the polluted environments and tested for their degradation potential and mechanisms. However, degrading microorganisms do not often exist in mono-culture in the environment. Instead, many bacteria are found in close association with one another and exhibit a superior biodegradation performance than that by a single species (). For example, compared to single bacteria, a consortium of Pseudomonas sp. SUK1 and Aspergillus ochraceus NCIM-1146 shows an increased ability to break down the textile dye reactive navy blue HE2R in 24 h (). The biochemical cooperation of different microorganisms accelerates the complete degradation of some pollutants (Wu et al., 2010). In addition, some bacteria may display no degradation ability but protect the bacteria with degradation ability or change the structure of the contaminants to enhance the degradation process (; Zhao and Wong, 2009). Therefore, interactions between microbes play significant roles in hazardous substance degradation.
Tetrahydrofuran (THF) is a heterocyclic ether with the formula (CH2)4O. As a polar, versatile solvent, THF is widely used as a reaction media for the synthesis of polymers. Additionally, THF is easily detected in the groundwater () and can penetrate the skin, cause rapid dehydration, inhibit cytochrome P450, and induce central nervous system irritation, narcosis, edema, and colonic muscle spasms in animals (; ). As an eco-friendly and cost-effective strategy for pollutant removal, the application of THF-degrading microorganisms has been rarely investigated. To date, only a few strains of the genera Rhodococcus (; Yao et al., 2009; Tajima et al., 2012), Pseudonocardia (; ), and Pseudomonas () have been reported to bear the ability to use THF as the sole carbon source. In our previous studies, we isolated a THF-degrading strain, named Rhodococcus ruber YYL, which had a maximum THF degradation rate of 137.60 mg⋅h-1⋅g-1 YYL dry weight (Yao et al., 2009, 2013). In addition, the symbiotic non-THF-degrading Bacillus cereus MLY1 was isolated with strain YYL from the activated sludge and could successfully augment the strain YYL colonization in activated sludge and remarkably improve the THF removal in the reactor (Yao et al., 2013). However, the mechanism of the non-THF-degrading microorganism’s survival in the mineral medium and the interactions with strain YYL remain unclear.
The environmental pH has been one of the key factors affecting the microbial degradation of pollutants (; ; ). Strain YYL shows its optimal degradation efficiency at initial pH 8.3 (Yao et al., 2009). However, actual environment conditions cannot meet this optimal pH. Strain YYL produces acidic substances continuously during the degradation of THF and decreases the pH of its environment, negatively affecting its degradation efficiency. We found that the THF degradation efficiency of YYL could be improved when under low pH stress (initial pH = 7.0) by the addition of the strain MLY1.
To study the interactions between strains YYL and MLY1, these two strains were co-cultured under the initial pH of 7.0 or 8.3 and subjected to metatranscriptome analysis. Here we focus on the difference in THF removal and gene expression in strain YYL and MLY1 between the co-culture and mono-culture systems with different pH conditions.
Materials and Methods
Strains, Culture Conditions, and Co-culture Experiments
The THF-degrading strain R. ruber YYL was cultured and maintained in 100 mL liquid optimal base mineral medium (BMM) at initial pH of 8.3 with 20 mM THF (Yao et al., 2009); meanwhile, B. cereus MLY1 was also cultured and maintained in 100 mL BMM with 1.0 g/L yeast extract. For the initial phase of each co-culture experiment, the stationary phase cells of strain YYL or MLY1 were centrifuged (5 min, 8000 rpm), the growth media was decanted, and the cells were suspended in BMM at the initial pH of 8.3 or 7.0. Subsequently, 100 mL BMM (initial pH = 8.3 or 7.0) with 20 mM THF was inoculated with 2 mL of strain YYL (OD600 = 1.5) and 1 mL of MLY1 (OD600 = 1.5) to make a cell ratio of 2:1. The corresponding mono-culture systems were incubated with the same inoculum of YYL and 1 mL of BMM. One Erlenmeyer flask of culture was used for each sample, and four replicates were prepared for every sampling time point in the four treatment groups. All the flasks were cultured at 140 rpm and a temperature of 30°C.
Sample Collection and Detection
Samples were collected for the determination of the THF concentration and the pH during the incubation. THF concentrations were measured by GC-2014C gas chromatography equipped with a flame ionization detector (FID) and an AOC-20i auto injector (SHIMADZU, Japan). The injector, oven, and detector temperatures were set to 200, 160, and 200°C, respectively. The THF peak was observed at a retention time of 1.7 min. For the transcriptional analysis, the sampling time for the experiment was chosen based on the THF removal efficiency and strain growth stage in the four treatment groups. In this study, samples were collected at the initial stage, the first 4 days, during which the cells of all the treatments were at the exponential growth phase. Also, the transcriptional changes in strain YYL that could describe its response within the symbiotic system were identified. All the collected samples were immediately centrifuged at 4°C, flash-frozen in liquid nitrogen, and stored at -80°C. To detect strain MLY1 growth, the ratio of the DNA levels of the specific genes thm/GerM were quantified using quantitative PCR (qPCR) to detect the ratio of strain YYL to strain MLY1 in co-culture (; ); thm is the THF monooxygenase gene responsible for the first step of THF degradation as described in previous studies (), and the GerM gene encodes a lipoprotein that stabilizes the GerA-GerQ complex through an interaction with the remodeled cell wall during spore formation in bacilli (). The primers for the qPCR are shown in Supplementary Table S1.
RNA Isolations, Library Construction, and Sequencing
Total RNA was extracted from four biological replicates per growth condition using the RNeasy Mini Kit (Cat#74106, QIAGEN, Germany) according to the manufacturer’s instructions. The RNA quality and yield were measured with gel electrophoresis using 1% agarose gel and NanoDrop 2000 Spectrophotometer (Thermo Scientific, Wilmington, DE, United States). Subsequently, the extracted quadripartite RNA was pooled equally and purified with the RNeasy Micro Kit (Cat#74004, QIAGEN, Germany) and the RNase-Free DNase Set (Cat#79254, QIAGEN, Germany) (Zhang and Gant, 2005; ; ; Zhao et al., 2017). The success of DNA removal was confirmed by PCR of the 16S rRNA gene with primers 27F and 1492R () (Supplementary Table S1). The reaction contained 5 μL of the 10x LA Taq buffer (Mg+2), 4 μL of dNTPs (2.5 mM), 1 μL of primer 27F (10 μM), 1 μL of primer 1492R (10 μM), 1 μL of the template DNA/RNA, 0.5 μL of the LA Taq polymerase, and 37.5 μL of ddH2O. The PCR conditions used were 94°C for 30 s and 28 cycles of 55°C for 30 s and 72°C for 90 s. Again, the integrity and quantity were assessed using NanoDrop 2000 Spectrophotometer and Agilent Bioanalyzer 2100 (Agilent Technologies, United States). The depletion of ribosomal RNA before cDNA synthesis was performed using the Ribo-Zero kit (Epicentre Biotechnologies, United States) for meta-bacteria. The mRNA was fragmented using the Ambion RNA fragmentation kit (Ambion, United States), and double-stranded cDNA was generated using the Qubit dsDNA HS Assay Kit (Invitrogen, United States). Subsequently, four mRNA libraries were constructed, evaluated using the Agilent Bioanalyzer 2100, and sequenced using 2 × 125 paired-end reads on an Illumina HiSeq2500 sequencer using the HiSeq SBS Kit v4 (Illumina, United States). Approximately 5 Gbp of clean data was targeted for each library. There were 57,246,312; 38,944,960; 43,907,548; and 47,725,396 reads of total clean data obtained, and 57,189,768; 38,862,480; 43,867,792; and 47,697,494 reads of clean data were left after rRNA removal for the mono- and co-culture performed with initial pH of 7.0 and 8.3.
RNA-Seq Data Analysis
Raw reads were filtered to remove low quality reads using Seqtk1 with the following standards: (1) reads with adaptors were removed; (2) reads containing more than 50 bases with low quality (Q20) were removed; (3) reads with more than 3 N bases were removed; (4) low quality bases or N bases assigned at the 3′ tail were removed; and (5) reads shorter than 20 bp were removed. To eliminate all ribosomal RNA sequences, reads mapping to the rRNA (5S rRNA, 16S rRNA, and 23S rRNA) of YYL and MLY1 were removed, and the remaining clean reads were used for the subsequent analysis.
All the clean reads from the mono- and co-culture systems were pooled and assembled using the Trinity assembly algorithm for Primary UniGene and CAP3 EST for Final UniGene (First_contig and Second_contig) (; ). The acquired Final UniGene sequences were searched against the NCBI non-redundant (nr) database (March 2014) using BLASTx (version BLAST-2.2.28+) () for the detection of protein-coding genes with the parameter e-value < 10-30, identity >60%, and genes were identified using self-developed Perl scripts; the maximum overlap against adjacent genes was 100 bp. Subsequently, further redundancy was removed using CD-HIT-EST v4.6 () with a sequence identity threshold of 99% in every 1000 bp (Wong et al., 2015). A collection of 9,483 non-redundant genes were used as a reference genome for the differential gene analysis.
Clean reads from each condition were aligned to the whole-body transcriptome of YYL. The normalized output for each gene expression was calculated as Reads Per Kilobase per Million mapped reads (RPKM) (). For Differentially Expressed Gene (DEG) sets, hierarchical clustering analysis was performed using the complete linkage and Euclidean distance as a measure of similarity. All data analyses and visualization of DEGs were conducted using R3.0.22. Gene expression under different pH conditions in the mono- and co-culture systems were compared, and the differential expression was considered as a fold-change ≥2 and false discovery rate (FDR) < 0.05.
All the gene sequences were searched against the Gene Ontology (GO) (e-value < 10-5), Kyoto Encyclopedia of Genes and Genomes (KEGG), and evolutionary genealogy of genes Non-supervised Orthologous Groups (eggNOG) using BLASTx (version: BLAST-2.2.28+) (; ). To confirm which GO terms and metabolic pathways from the mono- and co-culture systems responded to the different pH conditions, GO and KEGG enrichment analyses were performed using DAVID (). The significantly enriched GO terms and KEGG pathways in the DEGs were identified with the hypergeometric test in the entire genome background. The calculated p-value was corrected by the Bonferroni method with a threshold of p-value ≤ 0.05 to define the GO terms and KEGG pathways as significantly enriched. Subsequently, the major biological functions, the most important biochemical metabolic pathways, and signal transduction pathways with differentially expressed genes (DEGs) were identified by GO and pathway significance enrichment.
Sequencing data are available in the NCBI database under the accession numbers SRR5723771, SRR5723772, SRR5723773, and SRR5723774.
Results
THF Degradation in the Incubation Systems
As previously mentioned, the pH is a significant variable affecting the THF degradation efficiency of strain YYL, and its optimal initial pH for THF degradation is 8.3 (Yao et al., 2009). The degradation efficiency of strain YYL at the sub-optimal initial pH of 7.0 can be improved by co-culture of YYL with MLY1. Therefore, strains YYL and MLY1 were mono- or co-cultured at initial pH of 7.0 and 8.3 and subjected to degrading physiological analysis (Figure 1). When the THF concentration and pH between the co-culture and mono-culture systems were compared, there was no significant difference found under the initial pH of 8.3, but a lower THF concentration and pH were detected under the initial pH of 7.0 in the co-culture than in the mono-culture (Figures 1A–D). Meanwhile, the ratio of strain YYL/strain MLY1 was lower when under the initial pH of 7.0 than that observed under the initial pH of 8.3 (Figure 1E), suggesting that strain MLY1 plays a more significant role under low pH conditions. Additionally, a more obvious synergistic promotion of the THF degradation was detected between strain YYL and MLY1 with low pH stress (initial pH = 7.0).
FIGURE 1
Transcriptional Response of Strain YYL to pH Stress in Mono-culture
As discussed above, the low initial pH of 7.0 was unfavorable for the THF degradation by strain YYL. To explore the effects of pH on the gene expression of strain YYL, transcriptional differences in the mono-culture systems of strain YYL under the initial pH of 7.0 and 8.3 were compared. As shown in Figure 2A, clear differences were observed in the general patterns of the transcriptional responses to the two-different pH. Growing strain YYL under the initial pH of 7.0 led to both increases and decreases in gene expression when compared to growing at the initial pH of 8.3. To explore the functions of the DEGs from strain YYL, KEGG, and GO annotation were analyzed.
FIGURE 2
The KEGG enrichment analysis revealed that steroid degradation and biosynthesis, degradation of the aromatic compounds, xylene, and dioxin, RNA transport, spliceosome, non-homologous end-joining, alpha-linoleic acid metabolism, and bile acid biosynthesis were up-regulated when strain YYL was grown under the initial pH of 7.0 (Figure 2B). Also, the GO functional analysis showed that pH stress (initial pH = 7.0) induced up-regulation of genes involved in the transposase and catalase activity, response to oxidative stress, transposition, and heterothallic cell-cell adhesion, with down-regulation of the ATPase activity, exopolyphosphatase activity, DNA N-glycosylase activity, branched-chain amino acid transport, and cellular carbohydrate metabolic processes occurring (Figure 3).
FIGURE 3
In summary, strain YYL responded to the low initial pH (7.0) stress by enhancing steroid concentration, anti-oxidation, additional compound degradation, translation, and mutation activity. Meanwhile, the low initial pH (7.0) stress decreased the ATPase activity, carbohydrate metabolic activity, and cell respiration of strain YYL.
Transcriptional Response of Strain YYL to Co-culture Under Different pH Values
According to the THF degradation curves (Figures 1A,B), strain MLY1 exhibited a synergistic promotion of the THF degradation by strain YYL when grown with initial pH of 7.0 but had no such effect when grown with initial pH of 8.3. Therefore, we hypothesized that strain MLY1 would induce different effects on the gene expression of strain YYL in co-cultured systems grxown at a different pH. As shown in Figures 4A,C, the correlation of the gene expression from strain YYL between the mono-culture and co-culture indicates that strain MLY1 influences the gene expression of strain YYL.
FIGURE 4
The DEGs from strain YYL derived from comparing mono- and co-culture under the initial pH of 7.0 were also analyzed by KEGG and GO annotation. As shown in Figures 4B,D, significantly different KEGG pathways were up-regulated in mono-culture when grown with initial pH of 7.0. The depleted KEGG pathways from the co-culture system included RNA transport, primary bile acid biosynthesis, glycolysis, 2-oxocarboxylic acid metabolism, and amino acid biosynthesis. The GO functional analysis revealed that the depleted GO functions from the co-culture system included 113 terms, such as the fatty acid biosynthesis process, response to acid, amino acid biosynthesis process, ATPase activity, and malate dehydrogenase activity (Supplementary Figure S1A); the enriched GO functions from the co-culture system included 25 terms, such as the preribosome, phosphoprotein phosphatase activity, and isopentenyl diphosphate biosynthetic process (Supplementary Figure S1B).
In the co-culture grown with initial pH of 8.3, the depleted KEGG pathways from strain YYL included metabolic pathways, oxidative phosphorylation, valine, leucine, and isoleucine degradation, geraniol degradation, and biosynthesis of secondary metabolites (Figure 5D). The GO functional analysis revealed that the depleted GO functions from the co-culture system included nine terms, such as RNA binding, proton-transporting ATP synthase activity, and growth (Supplementary Figure S2); the enriched GO function terms from the co-culture system included only the ribonucleo-protein complex and cytosolic large ribosomal subunit.
FIGURE 5
In summary, strain MLY1 imposed more effects on the gene expression of strain YYL when grown with initial pH of 7.0 than when grown with initial pH of 8.3. With initial pH of 7.0, strain YYL exhibited a lower level of amino acid biosynthesis, fatty acid biosynthesis, ATPase activity, response to acid, and translation activity because of the influence of strain MLY1; with initial pH of 8.3, strain YYL exhibited a lower level of the functions involved in cellular structure and biological function, such as cell wall synthesis, RNA binding, and cell growth.
Transcriptional Response to pH Stress in Co-culture
In the mono-culture system, pH stress (initial pH = 7.0) induced significant alterations in the gene expression of strain YYL. To explore the differences in gene expression between the two pH values in the co-culture systems, the transcriptomes of strain YYL and MLY1 were analyzed.
As shown in Figure 5A, most of the DEGs from strain YYL were up-regulated when YYL was grown with initial pH of 7.0 in the co-culture system. The KEGG enrichment analysis of DEGs indicated that oxidative phosphorylation was up-regulated in strain YYL because of the pH stress (initial pH = 7.0) (Figure 5B). The GO functional analysis showed the enrichment of GO function terms involved in the structural constituent of the ribosome, ATP synthesis, translation, and growth; the depleted GO function terms included rRNA binding, ribosome, ribonucleo-protein complex, and cytosolic large ribosomal subunit (Supplementary Figure S3).
As shown in Figure 5C, strain MLY1 in the co-culture system exhibited different gene expression when grown under the initial pH of 7.0 and 8.3. The KEGG analysis of differential genes showed the up-regulation of pathways from strain MLY1 grown with initial pH of 7.0 such as the biosynthesis of amino acids, 2-oxocarboxylic acid metabolism, sucrose metabolism, and ribosome (Figure 5D). The GO functional analysis revealed only two GO functional terms that were up-regulated in strain MLY1 grown with initial pH of 7.0 (Supplementary Figure S4).
In summary, the effects of low initial pH (7.0) stress on the gene expression of strain YYL were less notable in the co-culture system than in the mono-culture system. Meanwhile, strain MLY1 enhanced 2-oxocarboxylic acid metabolism and the synthesis of amino acids, ATP, and ribosomes to help strain YYL respond to the low initial pH (7.0) stress in the co-culture system.
Discussion
The response to pH stress plays a vital role in the survival of Gram-positive bacteria, and the mechanisms utilized by the gram-positive bacteria in response to pH stress operate in different ways (). However, the THF-degrading strain YYL exhibits weak resistance to low pH stress, with a low growth rate and THF degradation efficiency under low pH stress (Yao et al., 2009). In this study, two genes involved in anti-oxidation and mutation in strain YYL were unregulated when YYL was grown with initial pH of 7.0 suggesting that this pH renders an environmental stress on strain YYL. To deal with this stress, the synergistic relationships among microorganisms can be beneficial in the natural environment (; ). Based on our analysis of the THF degradation efficiency, a synergistic relationship between strain YYL and the non-THF-degrading strain MLY1 exists when strain YYL experiences pH stress (initial pH = 7.0) but not when strain YYL is not stressed (initial pH = 8.3).
Based on these results, we propose a model for the cooperation between strains YYL and MLY1, shown in Figure 6. In co-culture system with THF as the sole carbon source, THF is degraded by YYL, and the easily usable intermediates are utilized by MLY1 as carbon sources (Figure 6). The unfavorable pH environment causes a significant change in the expression of genes involved in the metabolism of fatty acids required for lipid synthesis, ATP for energy, and amino acids for translation in strain YYL (Figure 3 and Supplementary Figure S1A), which is lethal for YYL if the pH stress cannot be relieved. THF degradation is an acid-producing process, and the protons produced need be transported out of the cell, a process that is dependent upon ATP () (Figure 6). Fortunately, the symbiotic strain MLY1, surviving on the intermediates of the THF degradation, provides lipids, ATP, and amino acids for strain YYL to relieve the pH stress (Figure 6). Based on the above interactions between strain YYL and MLY1, the symbiotic system is stable to completely degrade the THF.
FIGURE 6
Strain MLY1 Could Contribute Fatty Acids and Lipids to Strain YYL
In response to pH stress (initial pH = 7.0), strain YYL up-regulates the genes involved in steroid biosynthesis and degradation and bile acid biosynthesis. Steroids and bile acids are both sterols, and interactive transformation exists between the steroids and bile acids in the genus Rhodococcus (; Yam et al., 2011). Up-regulation of steroid metabolism plays a significant role in the response to environmental stress by Rhodococcus sp. (; ). This may be because of the fact that steroid metabolism is related to energy supply and membrane structure in microbial cells (). In addition, previous literature has reported that most of the enzymes in charge of steroid transformation belong to the family of cytochrome P450 monooxygenases (P450s) in the Actinobacteria (); this enzyme activity can be inhibited by THF (Urlacher and Girhard, 2012). Therefore, the THF inhibition of the P450s in strain YYL grown with initial pH of 7.0 might induce the up-regulation of the genes encoding for the P450s to ensure adequate P450 function for steroid transformation, for instance.
At initial pH of 7.0, the co-culture of strain YYL and strain MLY1 causes down-regulation of the bile acid and fatty acid biosynthesis necessary for lipid production in strain YYL (Supplementary Figure S1A). Simultaneously, strain MLY1 up-regulates the genes involved in the biosynthesis of fatty acids and lipids, including cutin, suberine, and wax (Figure 5). We propose that strain MLY1 might contribute the fatty acids and lipids to help strain YYL respond to the low pH stress in the co-culture system.
Flow of ATP from Strain MLY1 to Strain YYL
Bacteria frequently encounter environmental stress that generates a severe demand for ATP (); for example, exposure to low pH in the case of Escherichia coli requires ATP (). However, strain YYL exposure to pH stress (initial pH = 7.0) causes down-regulation of the genes involved in ATP production, such as the ATPase, cell respiration, and carbohydrate metabolic functions. Generally, ATPase activity, a good marker for the condition of the cell membrane and its enzymatic activity (), is essential for growth; the ATPase counteracts variations in the cytoplasmic pH by pumping protons out of the cells during low pH stress (; ). The down-regulation of ATPases by strain YYLunder low pH stress (initial pH = 7.0) is not beneficial for responding to the pH stress. Therefore, the down-regulation of genes involved in ATP production likely explains the poor efficiency of THF degradation by strain YYL when grown at initial pH of 7.0.
The lack of ATPase activity and extrusion of protons during low pH stress in strain YYL leads YYL to require adequate support from strain MLY1. Comparing strain YYL co-cultured with strain MLY1 using the optimal pH (initial pH = 8.3), we noted the down-regulation of genes involved in ATP synthesis coupled with proton transport and the tricarboxylic acid cycle in the strain YYL co-cultured under pH stress (initial pH = 7.0) (Figure 3). Simultaneously, strain YYL up-regulated the genes involved in starch and sucrose metabolism, which could generate ATP in response to the low pH stress.
Strain MLY1 Could Provide Amino Acids and Synthetic Proteins to Strain YYL
During environmental stress, organisms must tightly regulate the activity of translation because of the high-energy consumption of protein synthesis (Sorensen and Loeschcke, 2007; Yamasaki and Anderson, 2008). To deal with pH stress (initial pH = 7.0), strain YYL up-regulates the genes responsible for RNA transport but down-regulates the genes responsible for amino acid transport. These alterations in gene expression in strain YYL might illustrate a significant reprogramming of protein translation for strain YYL to respond to the pH stress.
In a synergistic relationship of bacteria, the exchange of amino acids frequently occurs (). When under pH stress (initial pH = 7.0), co-culture with strain MLY1 causes the down-regulation of valine, leucine, and isoleucine biosynthesis in strain YYL (Supplementary Figure S1B). Simultaneously, strain MLY1 up-regulates the genes involved in the biosynthesis of amino acids and translation (Figure 5D and Supplementary Figure S4). It is likely that strain MLY1 provides amino acids and some synthetic proteins for strain YYL in the co-culture system.
In summary, this work demonstrates that pH stress imposes an inhibitory effect on the THF-degrading strain YYL. Transcriptome analysis reveals that strain YYL up-regulates the genes involved in anti-oxidation, mutation, steroid and bile acid metabolism, and translation while simultaneously down-regulating the genes involved in ATP production when experiencing the low pH stress. MLY1 has no THF degradation activity, but it could provide fatty acids, ATP, and amino acids for strain YYL in response to the pH stress in the co-culture systems which might relieve the inhibition of strain YYL.
Limitation Statement
The results reported in this article are only preliminary as four biological replicates were pooled for each condition investigated prior to sequencing.
Statements
Author contributions
Performed experiments: ZuL and ZH. Analyzed data: ZuL, ZH, HH, and XR. Conceived and designed experiments: ZuL, ZH, HH, XR, AO, and ZhL. All authors have agreed to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.
Funding
This work was financially supported by the National Natural Science Foundation of China (Nos. 41630637 and 31422003).
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.2017.02297/full#supplementary-material
References
1
AltschulS. F.GishW.MillerW.MyersE. W.LipmanD. J. (1990). Basic local alignment search tool.J. Mol. Biol.215403–410. 10.1016/S0022-2836(05)80360-2
2
AshburnerM.BallC. A.BlakeJ. A.BotsteinD.ButlerH.CherryJ. M.et al (2000). Gene ontology: tool for the unification of biology. The Gene Ontology Consortium.Nat. Genet.2525–29. 10.1038/75556
3
BearsonS.BearsonB.FosterJ. W. (1997). Acid stress responses in enterobacteria.FEMS. Microbiol. Lett.147173–180. 10.1111/j.1574-6968.1997.tb10238.x
4
BenomarS.RanavaD.CardenasM. L.TrablyE.RafrafiY.DucretA.et al (2015). Nutritional stress induces exchange of cell material and energetic coupling between bacterial species.Nat. Commun.6:6283. 10.1038/ncomms7283
5
Buch-PedersenM. J.PedersenB. P.VeierskovB.PalmgrenM. G. (2009). Protons and how they are transported by proton pumps.Pflugers Arch.457573–579. 10.1007/s00424-008-0503-8
6
ChenJ. M.ZhouY. Y.ChenD. Z.JinX. J. (2010). A newly isolated strain capable of effectively degrading tetrahydrofuran and its performance in a continuous flow system.Bioresour. Technol.1016461–6467. 10.1016/j.biortech.2010.03.064
7
ChhabraR. S.ElwellM. R.ChouB.MillerR. A.RenneR. A. (1990). Subchronic toxicity of tetrahydrofuran vapors in rats and mice.Fundam. Appl. Toxicol.14338–345. 10.1016/0272-0590(90)90213-4
8
CotterP. D.HillC. (2003). Surviving the acid test: responses of gram-positive bacteria to low pH.Microbiol. Mol. Biol. Rev.67429–453. 10.1128/MMBR.67.3.429-453.2003
9
DayeK. J.GroffJ. C.KirpekarA. C.MazumderR. (2003). High efficiency degradation of tetrahydrofuran (THF) using a membrane bioreactor: identification of THF-degrading cultures of Pseudonocardia sp. strain M1 and Rhodococcus ruber isolate M2.J. Ind. Microbiol. Biotechnol.30705–714. 10.1007/s10295-003-0103-8
10
DeLongE. F. (1992). Archaea in coastal marine environments.Proc. Natl. Acad. Sci. U.S.A.895685–5689. 10.1073/pnas.89.12.5685
11
DeviK.MishraS. K.SahuJ.PandaD.ModiM. K.SenP. (2016). Genome wide transcriptome profiling reveals differential gene expression in secondary metabolite pathway of Cymbopogon winterianus.Sci. Rep.6:21026. 10.1038/srep21026
12
Di GioiaD.FambriniL.CoppiniE.FavaF.BarberioC. (2004). Aggregation-based cooperation during bacterial aerobic degradation of polyethoxylated nonylphenols.Res. Microbiol.155761–769. 10.1016/j.resmic.2004.05.015
13
DonovaM. V. (2007). Transformation of steroids by actinobacteria: a review.Appl. Biochem. Microbiol.435–18. 10.1134/S0003683807010012
14
GarbayS.LonvaudfunelA. (1994). Characterization of membrane-bound ATPase activity of leuconostoc-oenos-growth-conditions.Appl. Microbiol. Biotechnol.41597–602. 10.1007/BF00178496
15
GrabherrM. G.HaasB. J.YassourM.LevinJ. Z.ThompsonD. A.AmitI.et al (2011). Full-length transcriptome assembly from RNA-Seq data without a reference genome.Nat. Biotechnol.29644–652. 10.1038/nbt.1883
16
Hall-StoodleyL.CostertonJ. W.StoodleyP. (2004). Bacterial biofilms: from the natural environment to infectious diseases.Nat. Rev. Microbiol.295–108. 10.1038/nrmicro821
17
HaussmannU.WoltersD. A.FranzelB.EltisL. D.PoetschA. (2013). Physiological adaptation of the Rhodococcus jostii RHA1 membrane proteome to steroids as growth substrates.J. Proteome Res.121188–1198. 10.1021/pr300816n
18
HeZ.YaoY.LuZ.YeY. (2014). Dynamic metabolic and transcriptional profiling of Rhodococcus sp. strain YYL during the degradation of tetrahydrofuran.Appl. Environ. Microbiol.802656–2664. 10.1128/AEM.04131-13
19
HuangD. W.ShermanB. T.LempickiR. A. (2009). Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources.Nat. Protoc.444–57. 10.1038/nprot.2008.211
20
HuangX.MadanA. (1999). CAP3: a DNA sequence assembly program.Genome Res.9868–877. 10.1101/gr.9.9.868
21
HuertasM.SáezL.RoldánM.Luque-AlmagroV.Martínez-LuqueM.BlascoR.et al (2010). Alkaline cyanide degradation by Pseudomonas pseudoalcaligenes CECT5344 in a batch reactor. Influence of pH.J. Hazard. Mater.17972–78. 10.1016/j.jhazmat.2010.02.059
22
IsaacsonC.MohrT. K.FieldJ. A. (2006). Quantitative determination of 14-dioxane and tetrahydrofuran in groundwater by solid phase extraction GC/MS/MS.Environ. Sci. Technol.407305–7311. 10.1021/es0615270
23
KadamA. A.TelkeA. A.JagtapA. S.GovindwarS. P. (2011). Decolorization of adsorbed textile dyes by developed consortium of Pseudomonas sp. SUK1 and Aspergillus ochraceus NCIM-1146 under solid state fermentation.J. Hazard. Mater.189486–494. 10.1016/j.jhazmat.2011.02.066
24
KainkaryamR. M.BruexA.GilbertA. C.SchiefelbeinJ.WoolfP. J. (2010). PoolMC: smart pooling of mRNA samples in microarray experiments.BMC Bioinformatics11:299. 10.1186/1471-2105-11-299
25
KobayashiH.SuzukiT.UnemotoT. (1986). Streptococcal cytoplasmic pH is regulated by changes in amount and activity of a proton-translocating ATPase.J. Biol. Chem.261627–630.
26
KohlweyerU.ThiemerB.SchraderT.AndreesenJ. R. (2000). Tetrahydrofuran degradation by a newly isolated culture of Pseudonocardia sp. strain K1.FEMS Microbiol. Lett.186301–306. 10.1111/j.1574-6968.2000.tb09121.x
27
KrastanovA.AlexievaZ.YemendzhievH. (2013). Microbial degradation of phenol and phenolic derivatives.Eng. Life Sci.1376–87. 10.1002/elsc.201100227
28
KuhnertW. L.ZhengG. L.FaustoferriR. C.QuiveyR. G. (2004). The F-ATPase operon promoter of Streptococcus mutans is transcriptionally regulated in response to external pH.J. Bacteriol.1868524–8528. 10.1128/Jb.186.24.8524-8528.2004
29
KullenM. J.KlaenhammerT. R. (1999). Identification of the pH-inducible, proton-translocating F1F0-ATPase (atpBEFHAGDC) operon of Lactobacillus acidophilus by differential display: gene structure, cloning and characterization.Mol. Microbiol.331152–1161. 10.1046/j.1365-2958.1999.01557.x
30
LarkinM. J.KulakovL. A.AllenC. C. R. (2006). Biodegradation by members of the genus Rhodococcus: biochemistry, physiology, and genetic adaptation.Adv. Appl. Microbiol.591–29. 10.1016/S0065-2164(06)59001-X
31
LiW.GodzikA. (2006). Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences.Bioinformatics221658–1659. 10.1093/bioinformatics/btl158
32
LoweS. E.JainM. K.ZeikusJ. G. (1993). Biology, ecology, and biotechnological applications of anaerobic bacteria adapted to environmental stresses in temperature, pH, salinity, or substrates.Microbiol. Rev.57451–509.
33
MaX. K.DingN.PetersonE. C.DaugulisA. J. (2016). Heavy metals species affect fungal-bacterial synergism during the bioremediation of fluoranthene.Appl. Microbiol. Biotechnol.1007741–7750. 10.1007/s00253-016-7595-4
34
MalleyL. A.ChristophG. R.StadlerJ. C.HansenJ. F.BiesemeierJ. A.JastiS. L. (2001). Acute and subchronic neurotoxicological evaluation of tetrahydrofuran by inhalation in rats.Drug Chem. Toxicol.24201–219. 10.1081/DCT-100103720
35
MasudaH.McClayK.SteffanR. J.ZylstraG. J. (2012). Biodegradation of tetrahydrofuran and 14-dioxane by soluble diiron monooxygenase in Pseudonocardia sp. strain ENV478.J. Mol. Microbiol. Biotechnol.22312–316. 10.1159/000343817
36
MikeskovaH.NovotnyC.SvobodovaK. (2012). Interspecific interactions in mixed microbial cultures in a biodegradation perspective.Appl. Microbiol. Biotechnol.95861–870. 10.1007/s00253-012-4234-6
37
MortazaviA.WilliamsB. A.McCueK.SchaefferL.WoldB. (2008). Mapping and quantifying mammalian transcriptomes by RNA-Seq.Nat. Methods5621–628. 10.1038/nmeth.1226
38
NobuM. K.NarihiroT.RinkeC.KamagataY.TringeS. G.WoykeT.et al (2015). Microbial dark matter ecogenomics reveals complex synergistic networks in a methanogenic bioreactor.ISME J.91710–1722. 10.1038/ismej.2014.256
39
OrroA.CappellettiM.D’UrsiP.MilanesiL.Di CanitoA.ZampolliJ.et al (2015). Genome and phenotype microarray analyses of Rhodococcus sp. BCP1 and Rhodococcus opacus R7: genetic determinants and metabolic abilities with environmental relevance.PLOS ONE10:e0139467. 10.1371/journal.pone.0139467
40
RenC. G.BaiY. J.KongC. C.BianB.XieZ. H. (2016). Synergistic interactions between salt-tolerant rhizobia and arbuscular mycorrhizal fungi on salinity tolerance of Sesbania cannabina plants.J. Plant Growth Regul.351098–1107. 10.1007/s00344-016-9607-0
41
RodriguesC. D.Ramírez-GuadianaF. H.MeeskeA. J.WangX.RudnerD. Z. (2016). GerM is required to assemble the basal platform of the SpoIIIA–SpoIIQ transenvelope complex during sporulation in Bacillus subtilis.Mol. Microbiol.102260–273. 10.1111/mmi.13457
42
RosenthalA. Z.MatsonE. G.EldarA.LeadbetterJ. R. (2011). RNA-seq reveals cooperative metabolic interactions between two termite-gut spirochete species in co-culture.ISME J.51133–1142. 10.1038/ismej.2011.3
43
ShtratnikovaV. Y.SchelkunovM. I.FokinaV. V.PekovY. A.IvashinaT.DonovaM. V. (2016). Genome-wide bioinformatics analysis of steroid metabolism-associated genes in Nocardioides simplex VKM Ac-2033D.Curr. Genet.62643–656. 10.1007/s00294-016-0568-4
44
SorensenJ. G.LoeschckeV. (2007). Studying stress responses in the post-genomic era: its ecological and evolutionary role.J. Biosci.32447–456. 10.1007/s12038-007-0044-x
45
TajimaT.HayashidaN.MatsumuraR.OmuraA.NakashimadaY.KatoJ. (2012). Isolation and characterization of tetrahydrofuran-degrading Rhodococcus aetherivorans strain M8.Process. Biochem.471665–1669. 10.1016/j.procbio.2011.08.009
46
UrlacherV. B.GirhardM. (2012). Cytochrome P450 monooxygenases: an update on perspectives for synthetic application.Trends Biotechnol.3026–36. 10.1016/j.tibtech.2011.06.012
47
Vilchez-VargasR.JuncaH.PieperD. H. (2010). Metabolic networks, microbial ecology and ‘omics’ technologies: towards understanding in situ biodegradation processes.Environ. Microbiol.123089–3104. 10.1111/j.1462-2920.2010.02340.x
48
WongY. H.SunJ.HeL. S.ChenL. G.QiuJ. W.QianP. Y. (2015). High-throughput transcriptome sequencing of the cold seep mussel Bathymodiolus platifrons.Sci. Rep.5:16597. 10.1038/srep16597
49
WuX. L.LiangR. X.DaiQ. Y.JinD. C.WangY. Y.ChaoW. L. (2010). Complete degradation of di-n-octyl phthalate by biochemical cooperation between Gordonia sp. strain JDC-2 and Arthrobacter sp. strain JDC-32 isolated from activated sludge.J. Hazard. Mater.176262–268. 10.1016/j.jhazmat.2009.11.022
50
YamK. C.OkamotoS.RobertsJ. N.EltisL. D. (2011). Adventures in Rhodococcus - from steroids to explosives.Can. J. Microbiol.57155–168. 10.1139/W10-115
51
YamasakiS.AndersonP. (2008). Reprogramming mRNA translation during stress.Curr. Opin. Cell. Biol.20222–226. 10.1016/j.ceb.2008.01.013
52
YaoY.LuZ.ZhuF.MinH.BianC. (2013). Successful bioaugmentation of an activated sludge reactor with Rhodococcus sp. YYL for efficient tetrahydrofuran degradation.J. Hazard. Mater.261550–558. 10.1016/j.jhazmat.2013.08.007
53
YaoY. L.LvZ. M.MinH.LvZ. H.JiaoH. P. (2009). Isolation, identification and characterization of a novel Rhodococcus sp. strain in biodegradation of tetrahydrofuran and its medium optimization using sequential statistics-based experimental designs.Bioresour. Technol.1002762–2769. 10.1016/j.biortech.2009.01.006
54
ZhangS. D.GantT. W. (2005). Effect of pooling samples on the efficiency of comparative studies using microarrays.Bioinformatics214378–4383. 10.1093/bioinformatics/bti717
55
ZhaoY.BiK.GaoZ.ChenT.LiuH. Q.XieJ. T.et al (2017). Transcriptome analysis of Arabidopsis thaliana in response to Plasmodiophora brassicae during early infection.Front. Microbiol.2017:673. 10.3389/fmicb.2017.00673
56
ZhaoZ.WongJ. W. (2009). Biosurfactants from Acinetobacter calcoaceticus BU03 enhance the solubility and biodegradation of phenanthrene.Environ. Technol.30291–299. 10.1080/09593330802630801
Summary
Keywords
metatranscriptome, Rhodococcus ruber YYL, Bacillus cereus MLY1, cooperation, tetrahydrofuran degradation, low pH stress
Citation
Liu Z, He Z, Huang H, Ran X, Oluwafunmilayo AO and Lu Z (2017) pH Stress-Induced Cooperation between Rhodococcus ruber YYL and Bacillus cereus MLY1 in Biodegradation of Tetrahydrofuran. Front. Microbiol. 8:2297. doi: 10.3389/fmicb.2017.02297
Received
25 July 2017
Accepted
07 November 2017
Published
21 November 2017
Volume
8 - 2017
Edited by
Florence Abram, National University of Ireland Galway, Ireland
Reviewed by
Bernd Wemheuer, University of New South Wales, Australia; Svetlana Yurgel, Dalhousie University, Canada
Updates
Copyright
© 2017 Liu, He, Huang, Ran, Oluwafunmilayo and Lu.
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) or licensor 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: Zhenmei Lu, lzhenmei@zju.edu.cn
†These authors have contributed equally to this work.
This article was submitted to Microbial Symbioses, a section of the journal Frontiers in Microbiology
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.