| Home | E-Submission | Sitemap | Contact Us |  
Environ Eng Res > Volume 28(6); 2023 > Article
Zhou, Chen, Zhang, Huang, Tong, Lu, Zhang, and Ma: Succession of fungal communities and fungal-bacterial interactions in biofilm samples within a multistage bio-contact oxidation reactor during the treatment of low-COD and high-salinity produced water


The multistage bio-contact oxidation reactor (MBCOR) has been extensively applicated in the oily wastewater treatment. Nonetheless, the underlying mechanisms, especially the fungal community information during the bioremediation of low-COD and high-salinity oilfield produced water is still poorly understood. In this study, the fungal community succession, and fungi-bacteria interactions of a MBCOR were investigated. The result of Illumina high-throughput sequencing showed that the chemical oxygen demand was the main variable in influencing the community structure within the first tank. Besides, the linear discriminant analysis effect size analysis disclosed the feature fungal taxa at different stages, revealing that the fungal community shift was driven by its functional needs. Additionally, network analysis revealed that fungal and bacterial communities were highly interconnected and possessed cooperative relationships. Furthermore, the equilibrium of fungal-bacterial communities was primarily influenced by species richness and diversity. This study provides a deeper insight into the role of fungi within MBCOR in treating oily wastewater.

1. Introduction

The large volume of oilfield wastewater is one of the major environmental concerns due to its inclusion of various toxic compounds such as heavy metals, alkanes, BTEX (i.e., Benzene, Toluene, Ethylbenzene, and Xylene), polycyclic aromatic hydrocarbons (PAHs), ammonia nitrogen, and sulfide, posing potential hazards to human health [1]. Conventional physical and chemical methods used to address this problem could be cost-intensive and energy-consuming [2]. The use of microorganisms for bioremediation and detoxification of wastewater is considered as an efficient and lucrative alternative [3]. By using microbial systems, ecologically safe and environmentally friendly technologies can be designed to alleviate pollution.
Among the biological strategies for wastewater treatment, bio-contact oxidation (BCO) is a commonly used technology and is characterized as energy-efficient and cost-effective [4]. Of this system, contaminants were removed and biodegraded through the thorough contact of the biofilm, which was attached on the carriers and formed when the aerated wastewater flowed through the carriers [5]. Successful applications of the BCO reactor in oily wastewater treatment have been shown in previous literatures. For instance, Lu et al. [6] reported that the application of a hydrolysis acidification/BCO system effectively reduced the chemical oxygen demand (COD), ammonia nitrogen, total suspended solid, and total petroleum hydrocarbons in high-salinity oilfield-produced water. Similarly, a study conducted by He et al. [7] showed that the aerobic BCO reactor efficiently removed most of the investigated pollutants from oily wastewater. These studies proved BCO as a promising approach in the bio-treatment of oilfield wastewater.
To date, most studies related to wastewater bioremediation high-lighted the role of bacteria in pollutants degradation. Nevertheless, fungi are also potentially important contributors in various processes during the remediation and have been reported to contain a series of enzyme systems involved in the degradation, especially the well-known cytochrome P450 (CYP450) enzyme systems [8, 9]. They can participate in multiple pathways in the wastewater treatment systems, such as organic compounds degradation [10] and sludge flocs construction [9]. Some white rot fungi can produce extracellular polymeric substances, which could play a role in the biofilm formation [11]. Besides, fungi generally exhibited stronger tolerance against harsh environment than bacteria and could convert toxic compounds into low toxic and nontoxic end products, providing suitable conditions for bacterial growth [12]. Moreover, the filament of fungi could serve as a carrier for the immobilization of bacteria [13]. Overall, these observations suggest the necessities to comprehensively characterize the role of fungal communities in the wastewater treatment systems.
Currently, the pieces of information on fungal communities in wastewater bio-treatment and pollutants degradation remain vastly unexplored and the functions of fungi are still underestimated. Besides, recent studies have shown the strong ecological linkage and co-occurrence patterns between fungi and bacteria in different environments such as PAH-contaminated soil [14], chicken manure compost [15], and tillage [16]. A deep understanding of their interactions will provide critical information and instructive guidance for the modification and improvement of the performance of pollutants removal to meet the increasingly strict wastewater discharge standard. Nevertheless, few studies have been conducted to investigate the interconnection of fungal-bacterial communities in a multistage BCO reactor (MBCOR) during the treatment of oily wastewater, especially for the wastewater characterized by both low COD and high salinity, which could greatly inhibit the degradation process.
Since biofilms play major roles in decomposing pollutants in MBCOR, the analysis of microbial communities within biofilm samples is essential to uncover the bioremediation mechanisms. The rapid development of sequencing technologies provides convenience for the exploration of microbial community in wastewater treatment systems [17]. Our previous study has revealed the significant role of bacterial community within a MBCOR, which was located in Karamay, Xinjiang Uyghur Autonomous Region (45.68° N, 85.05° E), in treating low-COD and high-salinity oilfield wastewater [18]. In this study, we investigated the shift of fungal community in biofilm samples during the wastewater treatment through Illumina high-throughput sequencing. Meanwhile, the feature fungal taxa at different stages were analyzed via linear discriminant analysis (LDA) effect size (LEfSe) analysis. Furthermore, the interactions of fungal and bacterial communities were explored through network analysis. To the best of our knowledge, this study is the first to evaluate the fungal community succession in a MBCOR during the treatment of low-COD and high-salinity oilfield wastewater, providing useful information for deepening the understanding of pollutants removal in such wastewater treatment systems.

2. Materials and Methods

2.1. Descriptions of MBCOR and Analytical Methods

The investigated MBCOR began to work on September 7th, 2018. The first two months were domestication stage which allowed the microorganisms to adapt to the environment of oilfield wastewater. Some basic characteristics of the wastewater were shown in Table S1. It was clearly observable that the wastewater possessed the properties of low-COD (249.3±45.6 mg/L), high-salinity (12969.3±1058.6 mg/L), and high-temperature (47.9±1.6°C). The schematic diagram of the MBCOR was shown in Fig. S1. Detailed operational processes of the MBCOR have been described in our previous study [18]. After acclimation, the MBCOR was continuously operated for 181 days under aerobic conditions at ambient temperature. Wastewater samples were collected throughout the operation process, transferred to sterile plastic tubes immediately after sampling, and preserved in a refrigerator (4 °C). COD, biochemical oxygen demand (BOD), oil content, and ammonia nitrogen in influent and effluent were analyzed by chlorine emendation method, dilution method, infrared spectrophotometry, and Nessler’s reagent colorimetry, respectively, as described previously [18].

2.2. Biofilm Sampling and DNA Extraction

A total of 18 biofilm samples were obtained, which were collected from different tanks (i.e., tank I, II, and III), depths (i.e., 1.5 m, 3.0 m, and 4.5 m), and stages (i.e., domestication stage and operation stage). The sampling of biofilm for domestication stage and operation stage was performed on day 10 and day 75, respectively. The collected biofilms could representatively manifest the feature of fungal community at these stages, as reflected by the discrepancy of pollutants removal performance. The biofilm samples were transferred to sterile glass bottles immediately after sampling and stored in a −80 °C freezer before further analysis. Total fungal genomic DNA was extracted using E.Z.N.A.® Water DNA Kit (Omega Bio-Tek, Norcross, USA) following manufacturer’s instructions. The DNA quality was checked by the ratio of absorbance at 260 nm and 280 nm using a NanoDrop 2000 spectrophotometer (NanoDrop Inc., Wilmington, DE) and 0.8% agarose gel electrophoresis.

2.3. Illumina High-throughput Sequencing

The extracted DNA was then used as the template to amplify fungal ITS gene with the primer sets of ITS1F (GGAAGTAAAA GTCGTAACAAGG) and ITS1R (GCTGCGTTCTTCATCGATGC). PCR products were purified using Cycle Pure Kit PCR (Omega Bio-Tek, Norcross, USA) and paired-end sequenced on an Illumina Miseq sequencing platform (Personalbio Co., Ltd., Shanghai, China) according to the standard protocol. Low-quality Reads and chimera sequences were removed using Trimmomatic and UCHIME, respectively. According to reported literatures, operational taxonomic units (OTUs) were clustered based on a 97% sequence identity threshold using UCLUST algorithm [1921]. Alpha diversity indices including Good’s coverage, ACE, Chao1, Shannon, Simpson, and observed species were calculated using Quantitative Insights Into Microbial Ecology (QIIME) software. The UNITE (release 5.0) database was used as the reference database for fungal OTU classification. The obtained raw reads were deposited in the NCBI Sequence Read Archive database (accession number SRP266535).

2.4. Data Analyses

Principal coordinates analysis (PCoA) based on bray-curtis distance was employed to reveal the fungal community similarities among different biofilm samples. Canonical correlation analysis (CCA) was performed by the genescloud tools (https://www.gene-scloud.cn) to evaluate the relationship between environmental variables and fungal community structure. Spearman correlation analysis was conducted by R (version 3.6.1) using “psych” package. The correlation is significant when P value is below 0.05, while it is extremely significant when P value is below 0.01. LEfSe analysis was applied to detect the potential biomarkers in fungal communities using OmicStudio tools (https://www.omicstudio.cn/tool). Phylogenetic investigation of communities by reconstruction of unobserved states (PICRUSt2) algorithm was employed to predict ecological functions of fungi on Bioincloud platform (https://www.bioincloud.tech/). The co-occurrence networks were constructed by Gephi software (version 0.9.2). The topological properties of Erdös–Réyni random networks were calculated by R using igraph, openxlsx, and tidyverse packages. The bacterial community data used for network construction were derived from the same MBCOR as the fungal community and have been reported in our previous study [18].

3. Results and Discussion

3.1. Fungal Community Characteristics

3.1.1. Fungal richness and diversity

The monthly average contents of COD, BOD, NH4+-N, and oil in influent and effluent were shown in Fig. S2. The average removal efficiencies for COD, BOD, NH4+-N, and oil during the operation stage were 65.5%, 64.8%, 81.0%, and 70.2%, respectively, suggesting that biofilms in MBCOR could effectively decompose pollutants comprised in the high-salinity and low-COD produced water [18].
To explore the roles of fungi in pollutants degradation, fungal community structures of different biofilm samples were investigated through high-throughput sequencing. Meanwhile, the reasonability of fungal abundance for wastewater treatment was verified by real-time quantitative PCR. The mean abundance of ITS gene for domestication and operation stage was 2.38±0.51×107 and 7.64±1.62×107 copies/g wet weight biofilm, respectively, which was close to the result reported by Wei et al. [22] obtained from activated sludge. With respect to the result of high-throughput sequencing, a total of 721,454 (ranging from 27655 to 82424) effective sequences for all samples were obtained. The obtained OTUs ranged from 11 to 48 (Table S2). The coverage indices were in range of 0.992–0.998, demonstrating the effectiveness of constructed sequence libraries in covering the fungal communities (Fig. S3). Alpha diversity indices showed that the species richness was significantly (P<0.05) higher at domestication stage than operation stage (Fig. 1a), demonstrating the stabilization of fungal community after acclimatization. Similarly, Shannon and Simpson indices decreased slightly after the domestication stage, indicating the reduction of fungal diversity of the investigated biofilm samples. The result suggested that certain components in produced water played significant roles in shaping the fungal community structure, eliminating maladaptive microorganisms, and simplifying community compositions. A similar pattern was shown in bacterial community from the same MBCOR [18]. Besides, fungal diversity in tank I was significantly higher than other tanks (Table S2), which is probably due to the high pollutants content in this tank, serving as substrates and supporting the proliferation of diverse fungal genera. Moreover, it was found that the depth did not significantly (P>0.05) affect α-diversity indices of microbial community as tested by one-way ANOVA.

3.1.2. Fungal community structure

The fungal taxonomic distributions of different biofilm samples were investigated to explore the shift of community structure. The multistage structure of MBCOR resulted in a differentiation in fungal community (Fig. 1b), especially the tank I and other tanks, which primarily resulted from the discrepancy of environmental conditions in different tanks [18]. Seven major fungal phyla were detected in all samples, among which Zygomycota (7.37%–78.36%) dominated the fungal community. This phylum was reported to be significantly associated with primary and secondary metabolites such as sugar, amino acids, and organic acids [23]. The sharp increase in relative abundance of this phylum in tank II and tank III is probably related to the high concentration of metabolites in these tanks received from tank I. In addition, the proportion of Ascomycota raised from 0.95%–5.11% at domestication stage to 7.04%–23.76% at operation stage within tank I. This phylum is the largest group in aquatic fungal kingdom, contributing for more than 65% of the currently known fungi and key fungal communities for the decomposition of organic pollutants in nature [24]. Similarly, the abundance of Basidiomycota elevated in in tank I at operation stage (Fig. 1b). Members belonged to phylum are also known to mineralize organic contaminants [24, 25]. The enrichment of these phyla in tank I implies the significant role of this tank in transforming the contaminants, which is in accordance with the result obtained in our previous study [18]. Meanwhile, this phylum possesses the feature of heat tolerant [26], which also enabled its proliferation within the wastewater of high temperature (>45°C) in tank I. Moreover, considerable proportions of fungal populations remained unclassified or uncultured, especially for biofilm samples in tank I, suggesting the fungal community was rather complex.
To obtain a better resolution of fungal community, taxonomic distributions at genus level was further investigated. Nineteen major genera were identified, and the result was presented in Fig. 1c. Since the MBCOR was operated under aerobic conditions, most of these fungal populations were aerobes or facultative aerobes. Among the identified genera, Basidiobolus is predominant in the majority of biofilm communities with the relative abundance ranging from 4.34% to 78.36%. This genus is a typical saprobe and usually exists in soil or plant detritus, playing important roles in the decomposition of organic matters [27]. In this study, the high abundance of this genus implied that it might participate in pollutants transformation. Besides, Bullera (8.03%–23.61%) was found to occupy high proportions in tank I at domestication stage. Bullera was reported to be capable of secreting extracellular polysaccharides which are important substrates for biofilm formation [28]. At the operation stage, when the biofilm was mature and strong enough to resist environmental disturbance, the relative abundance of this genus decreased dramatically. In addition, Aspergillus is a halophilic fungus and was affiliated with PAHs degraders, possessing strong abilities in utilizing PAHs and other xenobiotic compounds in saline conditions [29]. Meanwhile, this genus is capable of producing biosurfactant which can assist the hydrophobic pollutants biodegradation [30]. The increase in relative abundance of this genus in tank I at operation stage demonstrated its potential role in reducing oil content (Fig. 1c). Other fungal genera such as Penicillium, Simplicillium, Candida, Rhizopus, and Scleroderma were all previously reported to be associated with hydrocarbons degradation, biosurfactant production, biofilm formation, and ammonia nitrogen assimilation [3134].
PCoA analysis (OTU level) can reveal the community similarities among different biofilm samples. As shown in Fig. 2a, the axis1 and axis2 of PCoA explained 57.62% and 21.53% of the data distribution, respectively. An apparent fungal community differentiation was observed in biofilm samples between tank I and other tanks. This could be ascribed to the high hydraulic loading and pollutants concentration within the first tank, selectively enriching certain microbial populations that can tolerate harsh environments. An analogous result was observed in bacterial communities [18]. With respect to tank II and tank III, the fungal communities exhibited high similarities, which is probably due to their reception of wastewater with similar constituents. The CCA results further illustrated the importance of environmental factors in shaping the fungal community structure (Fig. 2b). The axis1 and axis2 of CCA explained 49.52% and 18.60% of the microbial community structure, respectively. Apparently, BOD, COD, and oil content were major variables in influencing the fungal communities within tank I, among which COD was the most important factor. A similar result has been reported by Cortés-Lorenzo et al. [35] that COD is one of the most important factor in influencing the fungal community structure within a submerged fixed bed bioreactor during the treatment of urban saline wastewater.

3.2. Biomarker Identification and Functional Profiles Prediction

To identify the feature fungal taxa in biofilm samples at different stages, LEfSe analysis was conducted, and the result was presented in Fig. 3. Forty fungal clades presented statistically significant difference (LDA value > 2.0, P<0.05) at two stages. The representative fungal genera with variations in abundance at domestication stage included Cantharellus, Roccella, and Metschnikowia. These genera have been shown to carry genes related to the synthesis of extracellular lipid, protein, and polysaccharides, such as Acc gene encoding acetyl-CoA carboxylase and Fas gene encoding fatty acid synthase, which enables them to play certain roles in biofilm construction and formation [36, 37]. With respect to the operation stage, members including Phanerochaete and Pichia are typical PAH degraders and have been reported to contain CYP450 genes, implying that they might involve in PAHs biodegradation during the wastewater decontamination [3840]. Besides, Geastrum is capable of producing a variety of exoenzymes such as laccase, manganese peroxidase, and lipase, indicating its potential role in catalyzing the decomposition of organic pollutants [41]. Thus, the feature fungal taxa at different stages suggested the fungal community shift is driven by functional requirements.
The newly updated tool PICRUSt2 was applied to explore the functional profiles of fungal community in different biofilm samples. Enzymes related to carbon metabolism, nitrogen transformation, and biofilm formation such as peroxidase, laccase, monooxygenase, dioxygenase, and synthase, as reported in previous studies, were selected from the results of functional predictions [42, 43]. An obvious division in fungal function was observed at different stages (Fig. 4). Where, enzymes associated with biofilm construction were greatly enriched at domestication stage, while pollutants degradation related enzymes rose at operation stage. The result was in accordance with the pollutants treatment performance at different stages and confirmed that fungal community succession was function-driven. In addition, the result also implied that microorganisms could adapt to harsh environments by regulating the synthesis the required enzymes.

3.3. Co-occurrence Patterns of Fungal and Bacterial Communities

A co-occurrence network was constructed to elucidate the interactions between bacterial and fungal communities at different stages. As shown in Fig. 5, fungal and bacterial genera were highly interconnected and exhibited evident non-random co-occurrence patterns, indicating that fungi and bacteria have a synergistical cooperation during the wastewater treatment. Topological properties of constructed networks were calculated and compared with those of Erdös–Réyni random networks. Properties including clustering coefficient (CC), graph density (GD), modularity (MD), average path length (APL), and network diameter (ND) were all higher than those of the random networks (Table 1). The results suggested that the microbial networks possessed “small world” properties, that is nodes in constructed networks are more connected than in identically sized random networks [44]. Meanwhile, the MD values of networks were both higher than 0.4, suggesting the networks had modular structures [45]. The network at domestication stage (named as “D-network”) consisted of 156 nodes (genera) and 545 edges (367 positive edges), while the network at operation stage (named as “O-network”) comprised 141 nodes (genera) and 496 edges (294 positive edges). In addition, the AD value was slightly higher in O-network than in D-network, indicating a higher network complexity and microbial interactions after biofilm acclimation [46]. A similar pattern was shown in the fungal only networks (Fig. S4). The enhanced interconnections of microbes within and across microbial types could greatly elevate the microbial adaptions and resistance to environmental stress, thus leading to an improved pollutants degradation performance at operation stage [47].
Additionally, both networks contained five major modules, of which modules I–V accounted for 24.4%, 19.9%, 17.3%, 14.7%, and 14.7% of the whole network, respectively, at domestication stage (Fig. 5a). Meanwhile, the five modules were mainly occupied by Pseudomonadota (54.2%), Actinomycetota (7.7%), Chloroflexota (7.0%), and Bacteroidota (6.3%) (Fig. S5a). At operation stage, the microbial linkages were centralized within module III and V, demonstrating that the microbial interactions were enhanced during this process (Fig. 5b). About 18.4%, 17.0%, 17.0%, 14.2%, and 15.6% nodes were distributed in modules I–V, respectively. Pseudomonadota (46.6%), Chloroflexota (10.3%), Actinomycetota (8.6%), Ascomycota (7.8%), and Planctomycetota (6.0%) were the major phyla within these modules (Fig. S5b).
The potential keystone taxa or important nodes are usually identified by betweenness centrality, closeness centrality, degree centrality, and strongly connected ID [48]. In this study, between-ness centrality and degree centrality were used to screen the keystone genera in the fungal-bacterial networks at different stages. As shown in Figs. 5c–d, uncultured_bacterium_f_ Acidimicrobiaceae and Erythrobacter were identified as keystone taxa in microbial communities at domestication stage and operation stage, respectively. Generally, keystone taxa play vital roles in keeping the structural and functional stability of the microbial community under varied environmental conditions [48]. Nevertheless, in the present study, keystone taxa occupied low proportions in the microcosms with relative abundances of 0.07%–0.39% and 0.006%–0.50%, indicating that the proportions of genera were not well correlated with their status in the fungal-bacterial networks.

3.4. Equalization Analysis Between Fungal and Bacterial Community Structures

The community equalization is essential for microcosms to maintain the stability [49]. Therefore, an equalization analysis between fungal and bacterial communities was performed. As shown in Fig. 6a, the ITS/16S ratio decreased significantly (P<0.05) at operation stage, indicating a reduction in microbial equalization during this process [50]. The result suggested toxic compounds in the produced water were more likely to exterminate the sensitive fungal populations in the present study. The correlation between fungal and bacterial diversities was examined by Pearson’s correlation analysis, and negative correlations between the alpha diversity indices including Shannon (r=−0.13, P=0.22), Ace (r=−0.37, P=0.13), and Chao1 (r=−0.31, P=0.21) of the microbial communities were found (Figs. 6b–d). This showed that the microbial equilibrium of the bacterial and fungal communities was regulated by the species richness and diversity [50]. Moreover, it could be speculated that a weak competition between fungal and bacterial populations existed in the community, which was detrimental for pollutants decontamination. To improve the wastewater treatment performance, further research is needed to equalize the fungi-bacteria system in the MBCOR.

4. Conclusion

This study investigated the fungal community succession and microbial interactions in the MBCOR during the treatment of low-COD and high-salinity oilfield produced water. Results showed that COD was the main environmental factor in shaping the fungal community within the first tank. Besides, LEfSe analysis suggested that the community structure succession was driven by functional requirements. Additionally, network analysis revealed a strong synergistical cooperation between fungal and bacterial communities, which is in favor of the wastewater treatment. Moreover, variations in species richness and diversity were the principal cause for the change in fungal-bacterial equilibrium in the MBCOR. The findings obtained in this study give deeper insights into the microbial mechanisms during the oily wastewater treatment, providing an instructive guidance for improving the treatment performance of the MBCOR.

Supplementary Information


This study was financially supported by Shaanxi Province Academy of Sciences Project (2021K-33), Zhoushan Science and Technology Department Project (2019C81056), the Fundamental Research Funds for the Central Universities (2020QNA4045), the Scientific and Technological Aid Project for Xinjiang (2021E02041), Western Young Scholars Project of Chinese Academy of Science (XAB2020YW01), and the Guangxi Key Laboratory of Theory and Technology for Environmental Pollution Control (No. 2001 K004).


Author contributions

H. Z. (Associate professor) conducted the conceptualization, methodology, investigation, formal analysis, data curation, and original draft writing.

C. C. (PhD student) conducted the methodology, software, investigation, formal analysis, and data curation.

N. Z. (Lecturer) conducted the investigation and data curation.

X. H. (Engineer) conducted the investigation and data curation.

Z. T. (PhD student) conducted the investigation and data curation.

M. L. (PhD student) conducted the investigation and data curation.

C. Z. (Associate professor) conducted the investigation and data curation.

Y. M. (Associate professor) conducted the validation, review, and project administration.

Conflict-of-Interest Statement

The authors declare that they have no conflict of interest.


1. Neto AG, Costa CSB. Survival and growth of the dominant salt marsh grass spartina alterniflora in an oil industry saline wastewater. Int. J. Phytoremediat. 2009;11(7)640–650. https://doi.org/10.1080/15226510902861727
crossref pmid

2. Sowani H, Deshpande A, Gupta V, Kulkarni M, Zinjarde S. Biodegradation of squalene and n-hexadecane by Gordonia amicalis HS-11 with concomitant formation of biosurfactant and carotenoids. Int. Biodeterior. Biodegrad. 2019;142:172–181. https://doi.org/10.1016/j.ibiod.2019.05.005

3. Zhou HH, Huang XM, Liang YP, et al. Enhanced bioremediation of hydraulic fracturing flowback and produced water using an indigenous biosurfactant-producing bacteria Acinetobacter sp. Y2. Chem. Eng. J. 2020;397:125348 https://doi.org/10.1016/j.cej.2020.125348

4. Yuan CL, Lu CH, Ma YG, et al. A novel method to treat old landfill leachate combining multi-stage biological contact oxidation (MBCO) and single-stage autotrophic nitrogen removal using anammox and partial nitrification (SNAP). Chem. Eng. J. 2019;359:1635–1643. https://doi.org/10.1016/j.cej.2018.10.233

5. Zheng TL, Li P, Ma XY, et al. Pilot-scale multi-level biological contact oxidation system on the treatment of high concentration poultry manure wastewater. Process Saf. Environ. Prot. 2018;120:187–194. https://doi.org/10.1016/j.psep.2018.09.013

6. Lu M, Zhang ZZ, Yu WY, Zhu W. Biological treatment of oilfield-produced water: A field pilot study. Int. Biodeterior. Biodegrad. 2009;63(3)316–321. https://doi.org/10.1016/j.ibiod.2008.09.009

7. He F, Fu PC, Xu CM. Using biological treatment of Henan oilfield-produced water: Pilot plant test study. Adv. Mater. Res. 2012;(361–363)593–597. https://doi.org/10.4028/www.scientific.net/AMR.361-363.593

8. Ning DL, Wang H, Zhuang Y. Induction of functional cytochrome P450 and its involvement in degradation of benzoic acid by Phanerochaete chrysosporium . Biodegradation. 2010;21(2)297–308. https://doi.org/10.1007/s10532-009-9301-z
crossref pmid

9. Niu LH, Li Y, Xu LL, et al. Ignored fungal community in activated sludge wastewater treatment plants: diversity and altitudinal characteristics. Environ. Sci. Pollut. Res. 2017;24(4)4185–4193. https://doi.org/10.1007/s11356-016-8137-4
crossref pmid

10. Xu LH, Tan ZH, Zhang CF, et al. Performance and microbial diversity of a full-scale oilfield wastewater treatment plant. Desalin. Water Treat. 2017;99:239–247. https://doi.org/10.5004/dwt.2017.21694

11. Li NJ, Liu J, Yang R, Wu L. Distribution, characteristics of extracellular polymeric substances of Phanerochaete chrysosporium under lead ion stress and the influence on Pb removal. Sci. Rep. 2020;10(1)17633 https://doi.org/10.1038/s41598-020-74983-0
crossref pmid pmc

12. Zhou HH, Huang XM, Bu KY, Wen F, Zhang DD, Zhang CF. Fungal proliferation and hydrocarbon removal during bio-stimulation of oily sludge with high total petroleum hydrocarbon. Environ. Sci. Pollut. Res. 2019;26(32)33192–33201. https://doi.org/10.1007/s11356-019-06432-z
crossref pmid

13. Yuan XY, Zhang XY, Chen XP, Kong DW, Liu XY, Shen SY. Synergistic degradation of crude oil by indigenous bacterial consortium and exogenous fungus Scedosporium boydii . Bioresour. Technol. 2018;264:190–197. https://doi.org/10.1016/j.biortech.2018.05.072
crossref pmid

14. Ma XD, Li X, Liu JX, et al. Soil microbial community succession and interactions during combined plant/white-rot fungus remediation of polycyclic aromatic hydrocarbons. Sci. Total Environ. 2021;752:142224 https://doi.org/10.1016/j.scitotenv.2020.142224
crossref pmid

15. Mao HL, Wang K, Wang Z, Peng J, Ren NQ. Metabolic function, trophic mode, organics degradation ability and influence factor of bacterial and fungal communities in chicken manure composting. Bioresour. Technol. 2020;302:122883 https://doi.org/10.1016/j.biortech.2020.122883
crossref pmid

16. Wang YJ, Liu L, Luo Y, et al. Mulching practices alter the bacterial-fungal community and network in favor of soil quality in a semiarid orchard system. Sci. Total Environ. 2020;725:138527 https://doi.org/10.1016/j.scitotenv.2020.138527
crossref pmid

17. Cao JS, Zhang T, Wu Y, et al. Correlations of nitrogen removal and core functional genera in full-scale wastewater treatment plants: Influences of different treatment processes and influent characteristics. Bioresour. Technol. 2020;297:122455 https://doi.org/10.1016/j.biortech.2019.122455
crossref pmid

18. Zhou HH, Chen CL, Zhou SX, et al. Performance and microbial community analysis of a bio-contact oxidation reactor during the treatment of low-COD and high-salinity oilfield produced water. Bioresour. Technol. 2021;335:125267 https://doi.org/10.1016/j.biortech.2021.125267
crossref pmid

19. Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26(19)2460–2461. https://doi.org/10.1093/bioinformatics/btq461
crossref pmid

20. Li XG, Ding CF, Zhang TL, Wang XX. Fungal pathogen accumulation at the expense of plant-beneficial fungi as a consequence of consecutive peanut monoculturing. Soil Biol. Biochem. 2014;72:11–18. https://doi.org/10.1016/j.soilbio.2014.01.019

21. Geml J, Gravendeel B, van der Gaag KJ, et al. The contribution of DNA metabarcoding to fungal conservation: diversity assessment, habitat partitioning and mapping red-listed fungi in protected coastal salix repens communities in the Netherlands. PLoS One. 2014;9:1–17. https://doi.org/10.1371/journal.pone.0099852
crossref pmid pmc

22. Wei ZY, Liu YY, Feng K, et al. The divergence between fungal and bacterial communities in seasonal and spatial variations of wastewater treatment plants. Sci. Total Environ. 2018;(628–629)969–978. https://doi.org/10.1016/j.scitotenv.2018.02.003
crossref pmid

23. Xu MP, Gao DX, Fu SY, et al. Long-term effects of vegetation and soil on the microbial communities following afforestation of farmland with Robinia pseudoacacia plantations. Geoderma. 2020;367:114263 https://doi.org/10.1016/j.geoderma.2020.114263

24. Assress HA, Selvarajan R, Nyoni H, Ntushelo K, Mamba BB, Msagati TAM. Diversity, co-occurrence and implications of fungal communities in wastewater treatment plants. Sci. Rep. 2019;9(1)14056 https://doi.org/10.1038/s41598-019-50624-z
crossref pmid pmc

25. Guest R, Smith DW. A potential new role for fungi in a wastewater MBR biological nitrogen reduction system. J. Environ. Eng. Sci. 2002;1(6)433–437. https://doi.org/10.1139/s02-037

26. Fang C, Yuan XR, Liao KK, et al. Micro-aerobic conditions based on membrane-covered improves the quality of compost products: Insights into fungal community evolution and dissolved organic matter characteristics. Bioresour. Technol. 2022;362:127849 https://doi.org/10.1016/j.biortech.2022.127849
crossref pmid

27. Al-Hatmi AMS, Balkhair A, Al-Busaidi I, et al. Basidiobolus omanensis sp. nov. Causing angioinvasive abdominal basidiobolomycosis. J. Fungi. 2021;7(8)1–14. https://doi.org/10.3390/jof7080653
crossref pmid pmc

28. Elinov NP, Anan’eva EP, Vitovskaya GA, Trushina OA. Extracellular polysaccharides of Bullera alba VKM Y-2141. Chem. Nat. Compd. 1990;26(2)139–142. https://doi.org/10.1007/BF00607528

29. Peidro-Guzmán H, Pérez-Llano Y, González-Abradelo D, et al. Transcriptomic analysis of polyaromatic hydrocarbon degradation by the halophilic fungus Aspergillus sydowii at hypersaline conditions. Environ. Microbiol. 2021;23(7)3435–3459. https://doi.org/10.1111/1462-2920.15166
crossref pmid

30. Kiran GS, Hema TA, Gandhimathi R, et al. Optimization and production of a biosurfactant from the sponge-associated marine fungus Aspergillus ustus MSF3. Colloid Surf. B-Biointerfaces. 2009;73(2)250–256. https://doi.org/10.1016/j.colsurfb.2009.05.025
crossref pmid

31. Luna JM, Rufino RD, Sarubbo LA, Campos-Takaki GM. Characterisation, surface properties and biological activity of a biosurfactant produced from industrial waste by Candida sphaerica UCP0995 for application in the petroleum industry. Colloid Surf. B-Biointerfaces. 2013;102:202–209. https://doi.org/10.1016/j.colsurfb.2012.08.008
crossref pmid

32. Aranda E, Godoy P, Reina R, et al. Isolation of Ascomycota fungi with capability to transform PAHs: Insights into the bio-degradation mechanisms of Penicillium oxalicum . Int. Biodeterior. Biodegrad. 2017;122:141–150. https://doi.org/10.1016/j.ibiod.2017.05.015

33. Liu XH, Xing XY, Dong QL, Liu W, Li WN. Efficient removal of nitrogen/phosphorous by mix-cultivation of Haematococcus pluvialis and Simplicillium lanosoniveum in wastewater supplemented with NaHCO3 . Biochem. Eng. J. 2022;182:108433 https://doi.org/10.1016/j.bej.2022.108433

34. Yu ZD, Song G, Liu J, Wang JY, Zhang PY, Chen KS. Beneficial effects of extracellular polysaccharide from Rhizopus nigricans on the intestinal immunity of colorectal cancer mice. Int. J. Biol. Macromol. 2018;115:718–726. https://doi.org/10.1016/j.ijbiomac.2018.04.128
crossref pmid

35. Cortés-Lorenzo C, González-Martínez A, Smidt H, González-López J, Rodelas B. Influence of salinity on fungal communities in a submerged fixed bed bioreactor for wastewater treatment. Chem. Eng. J. 2016;285:562–572. https://doi.org/10.1016/j.cej.2015.10.009

36. Li Q, Wang DY, Liu XL, Li AJ, Chandran K. Enhanced lipid accumulation in Metschnikowia pulcherrima using volatile fatty acids under non-sterile repeated batch cultivation. Int. Biodeterior. Biodegrad. 2021;163:105256 https://doi.org/10.1016/j.ibiod.2021.105256

37. Li Q, Liao M, Yang M, et al. Characterization of the mitochondrial genomes of three species in the ectomycorrhizal genus Cantharellus and phylogeny of Agaricomycetes . Int. J. Biol. Macromol. 2018;118:756–769. https://doi.org/10.1016/j.ijbiomac.2018.06.129
crossref pmid

38. Wang CP, Yu L, Zhang ZY, Wang BL, Sun HW. Tourmaline combined with Phanerochaete chrysosporium to remediate agricultural soil contaminated with PAHs and OCPs. J. Hazard. Mater. 2014;264:439–448. https://doi.org/10.1016/j.jhazmat.2013.10.073
crossref pmid

39. Pan F, Yang QX, Zhang Y, Zhang SJ, Yang M. Biodegradation of polycyclic aromatic hydrocarbons by Pichia anomala . Biotechnol. Lett. 2004;26(10)803–806. https://doi.org/10.1023/B:BILE.0000025882.33234.91
crossref pmid

40. Loice CN, Sinji H, Chie N, Hiroshi T, Hirofumi I, Hiroyuki W. Gene expression profiling and enzymatic function of Phanerochaete Chrysosporium cytochrome P450s involved in the metabolism of Benzo(a)pyrene. Int. J. Environ. Bioremediat. Biodegrad. 2016;4(3)85–92. https://doi.org/10.12691/ijebb-4-3-3

41. Kuhar F, Castiglia VC, Zamora JC. Detection of manganese peroxidase and other exoenzymes in four isolates of Geastrum (Geastrales) in pure culture. Rev. Argent. Microbiol. 2016;48(4)274–278. https://doi.org/10.1016/j.ram.2016.09.002
crossref pmid

42. Yang XS, Zhang ZJ, Li SZ, et al. Fungal dynamics and potential functions during anaerobic digestion of food waste. Environ. Res. 2022;212:113298 https://doi.org/10.1016/j.envres.2022.113298
crossref pmid

43. Li SZ, Deng Y, Lian SY, Dai CX, Ma Q, Qu YY. Succession of diversity, functions, and interactions of the fungal community in activated sludge under aromatic hydrocarbon stress. Environ. Res. 2022;204:112143 https://doi.org/10.1016/j.envres.2021.112143
crossref pmid

44. Ju F, Xia Y, Guo F, Wang ZP, Zhang T. Taxonomic relatedness shapes bacterial assembly in activated sludge of globally distributed wastewater treatment plants. Environ. Microbiol. 2014;16(8)2421–2432. https://doi.org/10.1111/1462-2920.12355
crossref pmid

45. Jiao S, Liu ZS, Lin YB, Yang J, Chen WM, Wei GH. Bacterial communities in oil contaminated soils: Biogeography and co-occurrence patterns. Soil Biol. Biochem. 2016;98:64–73. https://doi.org/10.1016/j.soilbio.2016.04.005

46. Zhao YN, Chen WD, Wen DH. The effects of crude oil on microbial nitrogen cycling in coastal sediments. Environ. Int. 2020;139:105724 https://doi.org/10.1016/j.envint.2020.105724
crossref pmid

47. Yin HQ, Niu JJ, Ren YH, et al. An integrated insight into the response of sedimentary microbial communities to heavy metal contamination. Sci. Rep. 2015;5(1)14266 https://doi.org/10.1038/srep14266
crossref pmid pmc

48. Li J, Cai MH, Miao Y, et al. Bacterial community structure and predicted function in an acidogenic sulfate-reducing reactor: Effect of organic carbon to sulfate ratios. Bioresour. Technol. 2019;293:122020 https://doi.org/10.1016/j.biortech.2019.122020
crossref pmid

49. Lax S, Cardona C, Zhao D, et al. Microbial and metabolic succession on common building materials under high humidity conditions. Nat. Commun. 2019;10(1)1–12. https://doi.org/10.1038/s41467-019-09764-z
crossref pmid pmc

50. Shi XL, Zhao XH, Ren JY, et al. Influence of peanut, sorghum, and soil salinity on microbial community composition in interspecific interaction zone. Front. Microbiol. 2021;12:1–13. https://doi.org/10.3389/fmicb.2021.678250
crossref pmid pmc

Fig. 1
(a) Alpha diversity indices including ACE, Chao1, Simpson, and Shannon in biofilm samples collected from domestication stage and operation stage; Taxonomic distribution of fungal phyla (b) and identified genera (c) in different biofilm samples. In Fig. 1c, taxa with an abundance of > 1% in at least one sample were shown. a indicates the month in which the biofilm samples were collected. b and c represent the stage (i.e., tank I, tank II, and tank III) and depth (i.e., 1.5 m, 3 m, and 4.5 m) of the oxidation reactor, respectively.
Fig. 2
(a) PCoA to illustrate the similarity of fungal communities in different biofilm samples; (b) CCA to evaluate the environmental factors on fungal community structures. The length of arrow is proportional to the importance of environmental factors in shaping the fungal community structure.
Fig. 3
LEfSe analysis showing the fungal biomarkers at different stages with LDA scores > 2.0. (a) Cladogram; (b) bar plot.
Fig. 4
Heatmap showing the functional profiles of fungal communities predicted by PICRUSt2. Values represent the transformed relative abundance.
Fig. 5
Co-occurring network of microbial communities in biofilm samples during do-mestication stage (a) and operation stage (b). Betweenness centrality and degree of each OTU in the networks of domestication stage (c) and operation stage (d). Co-occurring networks were colored by modularity class. A connection stands for a strong (Spearman’s ρ>0.9) and significant (P<0.01) correlation. Circular nodes represent bacterial genera, while triangle nodes represent fungal genera. The size of each node is proportional to its degree; the thickness of each connection between two nodes (edge) is proportional to the value of Spearman’s correlation coefficients.
Fig. 6
Equalization analysis of the microbial communities at different stages. Ratio of fungal and bacterial OTUs (a). The correlation between bacterial and fungal Shannon indices (b), Ace indices (c), and Chao1 indices (d).
Table 1
Topological Properties of Co-occurring Networks and Erdös Réyni Random Networks
Parameters Node Edge Average degree (AD) Clustering coefficient (CC) Graph Density (GD) Modularity (MD) Average path length (APL) Network diameter (ND)
Sept. network 156 545 6.987 0.221 0.045 0.576 3.668 9.00
Random network 156 545 6.987 0.045 0.045 0.304 2.797 5.15
Nov. network 141 496 7.035 0.296 0.050 0.447 4.224 10.00
Random network 141 496 7.035 0.050 0.050 0.302 2.736 5.08
PDF Links  PDF Links
PubReader  PubReader
Full text via DOI  Full text via DOI
Download Citation  Download Citation
Supplement  Supplement
Web of Science
Editorial Office
464 Cheongpa-ro, #726, Jung-gu, Seoul 04510, Republic of Korea
TEL : +82-2-383-9697   FAX : +82-2-383-9654   E-mail : eer@kosenv.or.kr

Copyright© Korean Society of Environmental Engineers.        Developed in M2PI
About |  Browse Articles |  Current Issue |  For Authors and Reviewers