Skip to main content

ORIGINAL RESEARCH article

Front. Microbiol., 05 August 2020
Sec. Microbiotechnology
This article is part of the Research Topic Microbial Transformations of Pollutants in Constructed Wetlands View all 4 articles

Nutrient Removal Process and Cathodic Microbial Community Composition in Integrated Vertical-Flow Constructed Wetland – Microbial Fuel Cells Filled With Different Substrates

\r\nFei Zhong*Fei Zhong1*Chunmei YuChunmei Yu1Yanhong ChenYanhong Chen1Xue WuXue Wu2Juan Wu,Juan Wu2,3Guoyuan LiuGuoyuan Liu1Jian ZhangJian Zhang1Zifa DengZifa Deng1Shuiping Cheng,*Shuiping Cheng2,3*
  • 1School of Life Sciences, Nantong University, Nantong, China
  • 2Key Laboratory of Yangtze River Water Environment, Ministry of Education, College of Environmental Science and Engineering, Tongji University, Shanghai, China
  • 3Shanghai Institute of Pollution Control and Ecological Security, Shanghai, China

An integrated vertical-flow constructed wetland-microbial fuel cell system (CW-MFC), consisting of an up-flow chamber and a down-flow chamber, was constructed to treat synthetic sewage wastewater. The performance of CW-MFCs filled with different substrates [i.e., ceramsite (CM-A), quartz (CM-B), and zeolite (CM-C) granules] under various hydraulic retention times (HRTs, 7.6, 4.0, and 2.8 d) was evaluated. Efficient and stable nitrogen (N) and phosphorus (P) removals were observed in CM-A under different HRTs, while the voltage outputs of the CW-MFCs was greatly reduced as the HRTs decreased. With an HRT of 2.8 d, the ammonium (NH4+-N) and orthophosphate (PO43–-P) removal efficiencies in CM-A were as high as 93.8 and 99.6%, respectively. Bacterial community analysis indicates that the N removal in the cathode area of CM-A could potentially benefit from the appearance of nitrifying bacteria (e.g., Nitrosomonas and Nitrospira) and relatively high abundance of denitrifiers involved in simultaneous nitrification and denitrification (e.g., Hydrogenophaga, Zoogloea, and Dechloromonas) and denitrifying sulfide removal (e.g., Thauera). Additionally, the difference in N removal efficiency among the CW-MFCs could be partly explained by higher iron (Fe) content in milled ceramsite granules and higher abundance of denitrifiers with nitrate reduction and ferrous ions oxidation capabilities in CM-A compared with that in CM-B and CM-C. Efficient PO43–-P removal in CM-A was mainly ascribed to substrate adsorption and denitrifying phosphorus (P) removal. Concerning the substantial purification performance in CM-A, ceramsite granules could be used to improve the nutrient removal efficiency in integrated vertical-flow CW-MFC.

Introduction

Water eutrophication caused by the enrichment of nitrogen (N) and phosphorus (P) has become a worldwide environmental problem in recent years. To protect surface water bodies from eutrophication, rigorous N and P discharge regulations have been established. This, in turn, encourages the development of novel technologies for wastewater treatment.

Constructed wetlands (CWs) have been widely accepted as a preferable alternative for wastewater treatment because of their low initial investment costs, easy operation and maintenance, and good landscape integration. Meanwhile, CWs are being upgraded continuously to meet increasingly strict water pollutant discharge standards. Recently, the integration of CWs with microbial fuel cells (i.e., constructed wetland-microbial fuel cell system, CW-MFCs) has attracted worldwide attention (Doherty et al., 2015). In general, the efficiency of power generation and wastewater treatment is the focus of CW-MFC research. Concerning research on power generation, great efforts have been made to increase power density (Shen et al., 2018; Xu et al., 2019), though the usefulness of CW-MFCs system as an electricity generating device in real scenario has not yet been well examined. For research on wastewater treatment, attention has been extended from chemical oxygen demand (COD) removal to N removal (Xu et al., 2018).

More specifically for electricity generation, past research paid great attention to up-flow CW-MFCs, because such system could minimize dissolved oxygen (DO) at the anode while ensuring maximum availability in the cathode region, maximizing the redox gradient to generate an electrical current (Corbella et al., 2014). However, under the up-flow pattern, N removal might be affected due to insufficient nitrification and denitrification at the anode and cathode, respectively. Nitrification should be accomplished before the effluent enters into the cathode area to facilitate cathodic denitrification when ammonium (NH4+) is the dominant N species in the influent. Oon et al. (2017) improved the nitrification efficiency by supplying artificial aeration in the cathode region. However, it was usually difficult to balance nitrification and denitrification processes with supplementary aeration, especially given the fact that a high effluent nitrate concentration was often observed.

Wu et al. (2017) successfully improved the TN removal efficiency in an up-flow CW-MFC with closely spaced electrodes by applying central aeration and effluent recirculation. The improvement of the TN removal was ascribed to a higher abundance of denitrifiers and anaerobic ammonium oxidation bacteria in the bottom layer of the system. By comparing open- and closed-circuit CW-MFCs (down-flow), Wang et al. (2016) found that MFCs could improve N removal in CWs when the dominant N species in the influent was nitrate. The analysis of the anodic bacterial community showed that the relative abundance of nitrate-reducing bacteria (Dechloromonas, Desulfobulbus, Flavobacterium, Propionivibrio, and Geobacter) had been significantly increased in the closed-circuit CW-MFCs. Compared with effluent recirculation, treating the effluent from an up-flow CW-MFC by a follow-up down-flow CW-MFC would be a promising alternative, as this combined system decreases energy demand. Liu et al. (2019) developed an integrated vertical-flow (down-flow and up-flow in sequence) CW-MFC in which the anode and cathode were placed in the bottom of the down-flow chamber and the top of the up-flow chamber, respectively. In this configuration, the organic matter was consumed in front of the anode area because of aerobic oxidation. The traditional N removal pathway (i.e., aerobic nitrification and heterotrophic denitrification) might have been mainly responsible for the N removed. The study conducted by Liu et al. (2019) suggests that this type of integrated vertical-flow CW-MFC would be suitable for the treatment of high strength wastewater (e.g., swine wastewater). However, for wastewater with relatively low organic content, it is necessary to evaluate the performance of integrated vertical-flow CW-MFC consisting of an up-flow chamber and a down-flow chamber in sequence.

As an indispensable part of CW-MFC, substrate plays an essential role in N and P removal. Substrates (e.g., zeolite and ceramsite) with larger surface area and porous structure can improve the N and P adsorption and promote the development of biofilm (Wang et al., 2020). Yakar et al. (2018) investigated the effect of various types of substrates (sand, zeolite and volcanic cinder) on the performance of CW-MFC. Their study recommended that zeolite could be used to increase N and P removal efficiencies and bioelectric production in up-flow CW-MFC. However, it remains unclear whether or not the higher N removal efficiency in CW-MFC filled with zeolite is due to the promotion of active microbial groups involved in N removal. Using a pyrite-based down-flow CW-MFC, Ge et al. (2020) found enhanced nitrate and nitrite (NOX) removal via the intensification of the autotrophic denitrification process in the anode area. MFCs had been successfully operated using nitrate as a terminal electron acceptor in the cathode (Clauwaert et al., 2007). Wang et al. (2019) claimed that the MFC could reduce the dependence of denitrification on organic carbon sources in CWs, and the N removal in the CW-MFC under low influent COD/N ratios could be ascribed to bioelectrochemical denitrification with autotrophic denitrifying bacteria by accepting electrons from the cathode. The potential influence of substrates selection on N removal process in the cathode area of CW-MFCs is to be investigated.

This study compared the performance of integrated vertical-flow CW-MFCs (consisting of an up-flow and a down-flow chamber) filled with different substrates (i.e., ceramsite, quartz, and zeolite granules) treating synthetic sewage wastewater. It aims to examine the potential nutrient removal processes in this type of CW-MFC, and to explore the effects of substrate selection on the nutrient removal processes. Moreover, the composition of bacterial community in the cathode area was analyzed using 16S rRNA gene sequencing to investigate the possible pathways of N and P removal. This study could be beneficial to optimize the integrated vertical-flow CW-MFCs and encourage the development of novel CW-MFCs for use in reality.

Materials and Methods

Construction of Integrated Vertical-Flow CW-MFCs

A schematic diagram of the integrated vertical flow CW-MFCs used in this study is shown in Figure 1. Along the water flow direction, each CW-MFC was separated into an up-flow chamber and a down-flow chamber. Both chambers were made of polyvinylchloride (PVC) tubes with dimensions of 120 mm in diameter and 350 mm in height, and were connected at the top. The total volume of each CW-MFC was 7.92 L, with a liquid volume of around 4.12 L.

FIGURE 1
www.frontiersin.org

Figure 1. Schematic diagram of the integrated vertical-flow constructed wetland – microbial fuel cell systems. U4-U1 and D1-D4 are sampling ports located in the up-flow and down-flow chambers, respectively.

In the up-flow chamber, Canna indica was planted, and an air-cathode made of carbon felt was placed at the water–air interface. Granular graphite (supported by a stainless-steel net) were used as the anode, and the electrodes were connected by an external electrical resistor of 1,000 Ω. In the down-flow chamber, both electrodes were made of granular graphite, with an external electrical resistor of 3,000 Ω used to connect them. The anode and cathode of the CW-MFCs were inoculated with anaerobic digested sludge and activated sludge, respectively. The sludges were collected from a wastewater treatment plant nearby. The voltage outputs of the three CW-MFCs were recorded with an automatic recorder (Wuhan YAV Electronic Technology Co., Ltd., China, Smart V). The power density was calculated by dividing the power with the working volume of the anode.

There were eight sampling ports in the CW-MFCs along the water flow direction [U4-U1 (up-flow chamber) and D1-D4 (down-flow chamber)]. The CW-MFCs were filled with different substrates, i.e., ceramsite (CM-A), quartz (CM-B), and zeolite (CM-C) granules with a diameter of 4–8 mm, which were the same as those reported by Zhong et al. (2015). Time replicates (n = 6) were used to compare the differences in purification performances among the three CW-MFCs filled with different substrates. During the experiment, the CW-MFCs were operated under different HRTs [i.e., 7.6 d (H1), 4.0 d (H2), and 2.8 d (H3)]. The CW-MFCs were operated in intermittent flow mode using self-starting pumps three times each day (i.e., at 6:00, 12:00, and 18:00). The pump running time (within 30 min) was set according to the hydraulic loading rates. The hydraulic loading rates were 0.54 L/d (H1), 1.03 L/d (H2), and 1.47 L/d (H3), respectively. The wetland systems were operated continuously for a period of around 120 days from June to October in 2018. Before the experiment, the CW-MFCs had been stabilized for 3 months. The composition of synthetic wastewater was similar to that reported by Villaseñor et al. (2013) with small changes. It was composed of glucose (160 mg/L), CH3COONa⋅3H2O (160 mg/L), NaHCO3 (111 mg/L), KH2PO4 (12 mg/L), K2HPO4 (40 mg/L), MgCl2⋅6H2O (37.1 mg/L), CaCl2⋅2H2O (30.1 mg/L), and NH4Cl (90 mg/L).

Physicochemical Analysis

Water samples were collected weekly from the influent, effluent and eight sampling ports between 12:00 and 12:30. COD was determined using COD digestion vials (HACH, Loveland CO, 2125825) and a portable spectrophotometer (HACH, Loveland CO, DR2800). The concentrations of NH4+-N, NOX-N, orthophosphate (PO43–-P), sulfide (S2–-S), and sulfate (SO42–-S) were measured using an automatic chemical analyzer (DeChem-Tech, Hamburg, CleverChem anna). DO and pH were measured with the Thermo Orion 5 Star portable meter (Thermo-Orion Inc., Waltham, MA, United States).

Bacterial Community Analysis

The morphologies and elemental microanalysis of the original substrates and the substrates sampled at different layers of CW-MFCs were observed by scanning electron microscopy/energy dispersive X-ray spectrometry (SEM/EDS, Hitachi, Japan, S-3400N II) at the end of the experiment.

In each CW-MFC, three substrate samples were taken from the cathode area and mixed in polyethylene bags. A subsample was taken from each bag for bacterial community analysis. Any visible root or plant material was manually removed prior to homogenization. Samples were stored in a freezer at −80°C until Illumina sequencing analysis was conducted.

Genomic DNA was extracted from approximately 8 g of mixed substrate samples using an E.Z.N.A.® soil DNA Kit (Omega Bio-tek, Norcross, GA, United States.) according to the manufacturer’s instructions. The DNA extract was checked on a 1% agarose gel, and the DNA concentration and purity were determined with a NanoDrop 2000 UV-vis spectrophotometer (Thermo Scientific, Wilmington, DE, United States). The hypervariable V3–V4 region of the bacterial 16S rRNA gene was amplified with the primer pair 338F (5′-ACTCCTACGGGAGGCAGCAG-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′) by an ABI GeneAmp® 9700 PCR thermocycler (ABI, CA, United States). Purified amplicons were pooled in equimolar amounts and paired-end sequenced (2 × 300) on an Illumina MiSeq platform (Illumina, San Diego, CA, United States) according to the standard protocols by Majorbio Bio-Pharm Technology Co. Ltd. (Shanghai, China).

The raw 16S rRNA gene sequencing reads were demultiplexed, quality-filtered by Trimmomatic and merged by FLASH (Magoč and Salzberg, 2011). Operational taxonomic units (OTUs) with a 97% similarity cutoff were clustered using UPARSE1 (version 7.1), and chimeric sequences were identified and removed. The taxonomy of each OTU representative sequence was analyzed by RDP Classifier2 against the 16S rRNA database (e.g., Silva SSU128) using a confidence threshold of 0.7. Raw sequence data obtained was deposited in NCBI under the BioProject accession PRJNA629887.

Statistical Analysis

All data are presented as the means ± SD (standard deviation); n refers to the number of the samples. Statistical analysis was performed using the SPSS 19.0 software package for Windows. Differences in effluent water quality among the CW-MFCs filled with different substrates and operated under different HRTs were evaluated using one-way ANOVA.

Results

Water Purification and Electricity Generation Performance

The variations in N and P concentrations inside the CW-MFCs, both along the water flow direction and under different HRTs, are shown in Figure 2. It was observed that the NH4+-N concentration decreased sharply from layer U4 to layer D2, especially in CM-A. As the HRT decreased from 7.6 to 2.8 d, increased NH4+-N concentrations were observed in the sampling ports of the up-flow chamber of CM-A, the down-flow chamber of CM-B, and both chambers in CM-C. However, the CM-A effluent NH4+-N concentrations were mostly below 1.50 mg/L when the HRT decreased to 2.8 d. It is noteworthy that the effluent NOx-N concentrations in all sampling ports of the CW-MFCs were mostly below 0.3 mg/L. The variation in PO43–-P concentrations under different HRTs was similar to that in NH4+-N concentrations except for that of CM-B, demonstrating impressive PO43–-P treatment efficiency in CM-A. According to Table 1, both the effluent NH4+-N and PO43–-P concentrations in CM-A were significantly lower than those in CM-B and CM-C under different HRTs (P < 0.05).

FIGURE 2
www.frontiersin.org

Figure 2. Variation in NH4+-N, NOx-N, and PO43–-P concentrations at different layers [influent (in), U4-U1, D1-D4 and effluent (out)] of CW-MFCs filled with ceramsite (CM-A), quartz (CM-B), and zeolite (CM-C) granules under various hydraulic retention times [7.6 d (H1), 4.0 d (H2), and 2.8 d (H3)].

TABLE 1
www.frontiersin.org

Table 1. N, P, and COD concentrations in the influent and effluent of the CW-MFCs under different HRTs.

As the HRT decreased from H1 (7.6 d) to H3 (2.8 d), the COD concentrations in the effluent increased accordingly (Table 1). The lowest effluent COD concentrations in CM-B were observed among the three CW-MFCs under HRTs of 4.0 and 2.8 d. It is interesting to find that the effluent SO42–-S concentrations in CM-A were higher than those in the influent, while the effluent S2–-S concentrations of the CW-MFCs were negligible except for those of CM-B (Figure 3).

FIGURE 3
www.frontiersin.org

Figure 3. Variation in SO42–-S and S2–-S concentrations at different layers [influent (in), U4-U1, D1-D4 and effluent (out)] of CW-MFCs filled with ceramsite (CM-A), quartz (CM-B), and zeolite (CM-C) granules under various hydraulic retention times [7.6 d (H1), 4.0 d (H2), and 2.8 d (H3)].

With an HRT of 7.6 d (H1), the average voltage outputs in both the up-flow and down-flow chambers in CM-A were higher than those in CM-B and CM-C (Supplementary Figure S1). The average power densities in the up-flow chamber were 120.3, 11.3, and 14.2 mW/m3 in CM-A, CM-B, and CM-C, respectively. However, when the HRT was decreased to 4.0 d (H2) and 2.8 d (H3), low voltage outputs (below 0.02 V) with intermittent fluctuations in the three CW-MFCs were observed. Additionally, there were sharp decreases in voltage output in the up-flow chamber after the inflowing of synthetic wastewater.

Bacterial Morphologies and Community Structures in the CW-MFCs

The morphologies of the ceramsite granules at different layers were observed by SEM, and selected images are shown in Figure 4. Structures like mineral crystals were observed on the surface of ceramsite granule at layers U4 and D2. Via EDS analysis, an increase in calcium and a decrease in aluminum and silicon on the surface of the ceramsite granules were observed in layer AD2 compared with the original ones; and the P content on the surface of ceramsite granules sampled from layers AD1 and AD2 increased compared with the original ones (Supplementary Table S2). The crystalline mineral might be associated with Ca and P deposition. Additionally, the SEM images in Supplementary Figure S2 present morphology and structure of quartz and zeolite granules. Except for quartz granules, ceramsite and zeolite granules displayed porous structure and had large specific surface area. Both the SEM images and the results of EDS analysis (the variation in C and N contents) showed that the surface of substrates at layer D1 of the three CW-MFCs was covered by adhesive matrix, which was possibly excreted by the aggregated microorganisms.

FIGURE 4
www.frontiersin.org

Figure 4. SEM images of the surface of ceramsite granules (a: original ceramsite granule, b–d: ceramsite granules at layers U4, D1 and D2 of CM-A).

To evaluate the potential influences of substrate selection on the bacterial groups involved in N and P removal in the cathode region, the bacterial communities close to the cathodes of the CW-MFCs were determined by 16S rRNA gene sequencing. The results of the richness indices, including OTU numbers and the ACE, Simpson, Sobs and Chao1 indices, are summarized in Supplementary Table S3. The Chao1 index values, which are used to estimate the richness of the total bacterial community, were 559.60 for CM-A, 663.17 for CM-B and 628.40 for CM-C. A similar order was observed in the values of indicators such as OTU numbers, ACE, and Sobs. These findings show that the samples from CM-B had the highest richness, followed by those from CM-C and CM-A. According to the values of Shannon index and Simpson index, the bacterial diversities were in the order of CM-B > CM-A > CM-C. Good’s coverage values of all the three samples were above 0.99.

According to the community barplot analysis at the phylum level (Figure 5A), the majority of the sequences belong to the phyla Proteobacteria and Bacteroidetes, comprising 80.6, 65.1, and 80.7% of the total sequences in CM-A, CM-B, and CM-C, respectively. Within Proteobacteria, the β-subdivision was the predominant group, followed by γ-Proteobacteria, α-Proteobacteria, δ-Proteobacteria, and ε-Proteobacteria. At the family level, Comamonadaceae, Lentimicrobiaceae, and Rhodocyclaceae were dominant families in the CW-MFCs. The dominant families accounted for 50.9, 33.0, and 42.8% of the classified sequences in CM-A, CM-B, and CM-C, respectively.

FIGURE 5
www.frontiersin.org

Figure 5. Bacterial community composition at the phylum level (A) and relative abundance of potential functional genera responsible for nitrogen and sulfur cycling and electricity generation (B) in the CW-MFCs. The color intensity in each panel indicates the similarity or difference in characteristics among CM-A, CM-B, and CM-C; red and blue colors represent low or high enrichment of a genus, respectively (B).

Relative abundances of bacterial taxa that are potentially related to nitrification, denitrification, SO42– reduction, and electricity generation are presented in Figure 5B. The genera Nitrosomonas and Nitrospira accounted for 0.3, 0.8, and 0.3% of the detected OTUs in CM-A, CM-B and CM-C, respectively. High abundance of potential denitrifying bacteria (e.g., Azospira, Thauera, Hydrogenophaga, Zoogloea, Arenimonas, Rhodobacter, Flavobacterium, and Dechloromonas) at the genus level were observed, accounting for 35.8, 14.3 and 15.2% of the detected OTUs in CM-A, CM-B, and CM-C, respectively. The majority of the above genera belong to the β-subclass of the Proteobacteria. The abundance of the genera Desulfovibrio and Desulfomonile was higher in CM-B (2.2%) than in CM-A (0.5%) and CM-C (0.6%). A slightly higher abundance of the genera Klebsiella, Bacteroides, Rhodopseudomonas, and Arcobacter was observed in CM-B (3.1%) than in CM-A (2.6%) and CM-C (0.4%). Notably, the highest abundance of the genus Sphaerotilus (a potential iron-oxidizing bacterium) was observed in CM-A (2.2%) compared with CM-B (0.1%) and CM-C (0.1%).

Discussion

Nitrogen Removal Pathways

From layer AU4 to AU2, it was unexpected that efficient NH4+-N removal was observed. In contrast, the effluent NH4+-N concentrations remained almost the same from layer BU4 to BU2 under all the three HRTs, while the effluent NH4+-N concentrations decreased from layer CU4 to CU2 with an HRT of 7.6 d. The effluent pH values in the up-flow chamber of CM-A (above 8.00) were higher than those in CM-B and CM-C under all the three HRTs (Supplementary Table S1). The relatively higher pH values in the up-flow chamber of CM-A could contribute to the ammonia loss via volatilization. In CM-C, the NH4+-N could be adsorbed by the zeolite granules under relatively long HRT (7.6 d). The results indicated that the NH4+-N removal from layer U4 to U2 might be associated with the substrates selected.

Although the average effluent DO concentrations (in the range of 2.03 to 3.60 mg/L) in the up-flow chamber of the three CW-MFCs under three HRTs were not so low as to inhibit nitrification (Supplementary Table S1), it was unexpected to notice low abundance of ammonium-oxidizing bacteria (i.e., Nitrosomonas) and nitrite-oxidizing bacteria (i.e., Nitrospira) in the cathode area of CW-MFCs. However, promising N removal efficiency was observed in CM-A without the assistance of extra energy-demanding measures (i.e., aeration and effluent recirculation). Especially in the cathode area (from layer U2 to D2) of CM-A, efficient NH4+-N removal and low effluent NOx-N concentrations were observed as expected. The abundance of denitrifying bacteria in the cathode area in CM-A was the highest among the three CW-MFCs. Some of the denitrifying bacteria might serve versatile functions in the N removal process. Feng et al. (2015) indicated that certain groups of bacteria (e.g., Hydrogenophaga) were capable of heterotrophic nitrification and aerobic denitrification (HNAD) in the cathode chamber of an MFC. In a tidal flow CW, the genera Hydrogenophaga, Dechloromonas, and Zoogloea were recognized as HNAD bacteria, which were essential to achieve good N removal performance via simultaneous nitrification and denitrification (SND) process (Tan et al., 2020). In accordance with these studies, relatively higher abundance of Hydrogenophaga, Dechloromonas, and Zoogloea was observed in CM-A compared with that in CM-B and CM-C, which might contribute to the efficient N removal from layer AU2 to AD2 in this study. Moreover, Virdis et al. (2011) indicated that the SND process in the cathode chamber could be promoted in relatively high DO concentrations (5.73 ± 0.03 mg/L) if biofilm stratification (i.e., the nitrifying bacteria appear in the outer layer of the biofilm and the putative denitrifying organisms occupy the inner layer) was achieved. According to the SEM images (Figure 4C) and the composition of bacterial community, the SND process in the cathode area of CM-A could possibly benefit from the formation of biofilm on the surface of the ceramsite granules.

In the up-flow chamber of the CW-MFCs, the SO42–-S concentrations decreased from layers U4 to U2 at most times, while the S2–-S concentrations increased simultaneously compared with those in the influent (Figure 3). In the down-flow chamber of the CW-MFCs, the S2–-S was consumed from layers U2 to D1, suggesting the S2–-S might be involved in the denitrifying sulfide removal (DSR) process. Huang et al. (2017) stated that not only autotrophic bacteria can perform the DSR process, but also some heterotrophic bacteria are capable of DSR, achieving the simultaneous removal of NOx-N, S2–-S and organic carbon. Thauera was recognized as a sulfide-tolerant and facultative denitrifier (Mao et al., 2013; Liang et al., 2020), and it became the predominant bacterial genus when both sulfide and organic carbon were provided (Zhang et al., 2019). Liu et al. (2015) reported Thauera as one of the key strains in a denitrifying ammonium oxidation reactor treating organics-deficient wastewater with excess S0 production. In this study, relatively high abundance of genus Thauera in CM-A (7.72%) might also play a vital role in efficient N removal. In CM-B, the dominance of sulfate-reducing bacteria (e.g., Desulfovibrio and Desulfomonile) in the up-flow chamber might contribute to the increased COD removal efficiency and S2– accumulation when the HRT decreased to 2.8 d. Due to the low NH4+-N removal rate and nitrate production rate, the S2– consumption by DSR process would be decreased, which might be the reason why the highest effluent S2– concentration was observed in CM-B.

Moreover, the results of bacterial community analysis indicate that autotrophic denitrification using ferrous iron as electron donor might be involved in N removal process in the cathode area of CM-A. Sphaerotilus strains have the ability to deposit insoluble oxides or hydroxides of iron on their sheaths, providing resistance to growth inhibition caused by high ferric concentrations (Rogers and Anderson, 1976). In this study, relatively higher abundance of Sphaerotilus was observed in CM-A (2.21%) compared with that in CM-B (0.09%) and in CM-C (0.04%). The substrates filled in the CW-MFCs were identical to those used in Zhong et al. (2015), in which iron (Fe) content in the milled ceramsite, quartz, and zeolite granules were 41,791 ± 452 mg/kg, 8,096 ± 111 mg/kg, and 219 ± 19 mg/kg, respectively. The high Fe content in ceramsite granules might encourage the growth of Sphaerotilus in CM-A. The oxidation of ferrous ions and reduction of ferric ions, which are carried out by the Fe transformation functional bacterial groups, may have synergetic relationships with denitrifiers and contribute to denitrification (Zhang et al., 2018). Li et al. (2016) reported that Azospira, Zoogloea, and Dechloromonas played important roles in both the nitrate reduction and ferrous ions oxidation in the soil suspension. In CM-A, the relative abundances of Azospira, Zoogloea, and Dechloromonas were 11.17, 2.96, and 1.29%, respectively. Their dominance in CM-A would be beneficial to the denitrification process by donating electrons via oxidizing ferrous ions to ferric ions. Although the variations in ferric and ferrous concentrations were not measured in this study, it could be inferred that the efficient N removal in CM-A might be partly ascribed to the high Fe content in ceramsite granules and the internal transformation and cycling of ferric and ferrous ions.

The co-existence of autotrophic and heterotrophic denitrification processes had been reported in wastewater treatment systems under different influent COD/N ratios (Xu D. et al., 2017). Some genera (e.g., Zoogloea) are capable of both autotrophic and heterotrophic denitrification (Su et al., 2020), which might play a vital role under mixotrophic conditions. Additionally, the synthetic wastewater with relatively high COD/N ratio was first pumped into the anode area of the up-flow chamber in this study. The organic compounds would be partly consumed before arriving at the cathode area as reported (Oon et al., 2017; Yakar et al., 2018). The impact of organic compounds on the activity of the autotrophic denitrifiers would be limited. Xu D. et al. (2017) investigated the effects of COD/N on the simultaneous heterotrophic and hydrogen-based autotrophic denitrification performance in a bioelectrochemically assisted CW, and their study found that the N removal efficiency were significantly increased with increasing COD/N ratio. The combined heterotrophic and autotrophic denitrification processes could contribute to the relatively high and stable N removal performance in the cathode area of CM-A in this study.

Phosphorous Removal Pathways

Compared with quartz and zeolite granules, Zhong et al. (2015) recommended ceramsite granules as the best substrate for efficient P removal in a horizontal subsurface flow CW. The substrates filled in the CW-MFC were the same as those reported by Zhong et al. (2015), and the high calcium and Fe content and the heating process used to produce porous ceramsite granules could also be the reasons for the highest PO43–-P removal efficiencies in CM-A. The EDS analysis confirmed that the P content on the surface of ceramsite granules at layers AD1 and AD2 increased compared with the original ones. It has been reported that PO43–-P can be removed by metal oxides adsorbents, and the occurrence of ion exchange between SO42–-S and PO43–-P was observed (Wu et al., 2014). The increased effluent SO42–-S concentrations in CM-A could be partly related to the efficient removal of PO43–-P by the ceramsite granules.

Biological P removal might also contribute to the efficient PO43–-P removal in CM-A. Wu et al. (2013) reported that a synergistic relationship may exist between sulfur cycle and biological P removal. In the sulfur cycle-associated enhanced biological P removal (EBPR) process, poly-S served as one of the energy sources for P uptake and poly-P synthesis and ended up as sulfate; the genus Thauera was regarded as a key organism in the sulfur cycle-associated EBPR process and were assumed to participate in the denitrifying P and S2–-S removal from sulfate-rich wastewater (Guo et al., 2019). Additionally, the genus Dechloromonas has been reported as denitrifying poly-P accumulating organism that is capable of denitrifying P removal (Yuan et al., 2016). In this study, relatively high abundance of Thauera and Dechloromonas and SO42–-S concentrations were observed in CM-A, suggesting denitrifying P removal might take place in this type of CW-MFC filled with ceramsite granules.

Electricity Generation

The genera Klebsiella, Arcobacter, Bacteroides, and Rhodopseudomonas have been widely accepted as potential electrochemically active bacteria in the anode (Xia et al., 2010; Huang et al., 2014; Xiao et al., 2015). In this study, their appearance in the cathode area was observed, and only the abundance of Klebsiella was higher in CM-A compared with that in CM-B and CM-C. Whether or not they could be still electrochemically active in the cathode area deserve further investigation.

Additionally, Thauera was recognized as the major genus on an autotrophic denitrifying biocathode (Jiang et al., 2017). Yang et al. (2018) reported that selected Thauera-dominated cultures could achieve high performance in electricity generation in parallel with N removal. Besides involved in N removal, Rhodobacter and Hydrogenophaga were also thought to play a role in electron transfer from the cathode surface to the terminal electron acceptor (Huang et al., 2014). The predominance of Thauera, Hydrogenophaga, and Rhodobacter in the cathode area of CM-A might contribute to the enhancement of voltage output in CM-A compared with CM-B and CM-C.

Due to the low voltage output of the CW-MFCs under HRTs of 4.0 and 2.8 d, power density was not presented. During the experiment, the influent was pumped into the CW-MFCs in a short time. When the HRT was decreased, the influent wastewater with low DO concentrations might flow to the cathode area directly, resulting in the reduction in power output of CW-MFCs. While taking efforts to increase voltage output of this system, future studies should probe other approaches that can make full use of the generated electricity.

Output bioelectrical signals (e.g., voltage) of the CW-MFCs have been reported to correspond to changes in COD concentration (Xu L. et al., 2017; Corbella et al., 2019). In this study, it is noteworthy that the voltage output of the CW-MFCs was intermittently fluctuated. As the daily variation in voltage output was highly related to the time of pumping in this study (Supplementary Figure S1), the signals could be used to sense wastewater inflow. If the wastewater inflow was blocked due to clogging problems, it could be expected that the regular fluctuation in output voltage will change. Further studies should be conducted to examine whether bioelectrical signals could be regarded as an indicator of system clogging problems. If so, the decrease in nutrient removal efficiency induced by clogging problems could be timely detected via the abnormal fluctuation of bioelectrical signals.

Conclusion

CW-MFCs consisting of up-flow and down-flow chambers could remove NH4+-N and PO43–-P efficiently when filled with ceramsite granules. The N and P removal efficiencies in the CW-MFCs filled with ceramsite granules were relatively stable under different HRTs, while the voltage outputs were greatly reduced when the HRTs decreased from 7.6 to 2.8 d. The efficient N removal in CM-A could be ascribed to the appearance of nitrifying bacteria (e.g., Nitrosomonas and Nitrospira) and relatively high abundances of denitrifying bacteria in the cathode area that were potentially involved in SND (e.g., Hydrogenophaga, Zoogloea, and Dechloromonas), DSR (e.g., Thauera) and ferrous ion-dependent nitrate removal (e.g., Azospira, Zoogloea, and Dechloromonas). Substrate adsorption and denitrifying P removal via the genera Thauera and Dechloromonas might be the dominant pathways of P removal. Concerning the impressive purification performance, ceramsite granules are recommended rather than quartz and zeolite granules as substrates for the integrated vertical-flow CW-MFCs.

Data Availability Statement

We have deposited the data in NCBI with accession number PRJNA629887.

Author Contributions

FZ: idea. FZ, CY, YC, JW, JZ, and ZD: data analysis. CY, YC, XW, and GL: data collection. FZ and SC: writing and reviewing the manuscript. All authors contributed to the article and approved the submitted version.

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.

Acknowledgments

This research was supported by the Austrian-Chinese Scientific & Technological Cooperation project (CN 18/2016), China National Major Programme of Science and Technology (2017ZX07603-003), the Public Service Platform Project for Jiangsu Ocean High-end Equipment Testing and Application, and Large Instruments Open Foundation of Nantong University (KFJN2044).

Supplementary Material

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

Footnotes

  1. ^ http://drive5.com/uparse/
  2. ^ http://rdp.cme.msu.edu/

References

Clauwaert, P., Rabaey, K., Aelterman, P., De Schamphelaire, L., Pham, T. H., Boeckx, P., et al. (2007). Biological denitrification in microbial fuel cells. Environ. Sci. Technol. 41, 3354–3360. doi: 10.1021/es062580r

PubMed Abstract | CrossRef Full Text | Google Scholar

Corbella, C., Garfi, M., and Puigagut, J. (2014). Vertical redox profiles in treatment wetlands as function of hydraulic regime and macrophytes presence: surveying the optimal scenario for microbial fuel cell implementation. Sci. Total Environ. 470-471, 754–758. doi: 10.1016/j.scitotenv.2013.09.068

PubMed Abstract | CrossRef Full Text | Google Scholar

Corbella, C., Hartl, M., Fernandez-Gatell, M., and Puigagut, J. (2019). MFC-based biosensor for domestic wastewater COD assessment in constructed wetlands. Sci. Total Environ. 660, 218–226. doi: 10.1016/j.scitotenv.2018.12.347

PubMed Abstract | CrossRef Full Text | Google Scholar

Doherty, L., Zhao, Y., Zhao, X., Hu, Y., Hao, X., Xu, L., et al. (2015). A review of a recently emerged technology: constructed wetland-Microbial fuel cells. Water Res. 85, 38–45. doi: 10.1016/j.watres.2015.08.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Feng, C., Huang, L., Yu, H., Yi, X., and Wei, C. (2015). Simultaneous phenol removal, nitrification and denitrification using microbial fuel cell technology. Water Res. 76, 160–170. doi: 10.1016/j.watres.2015.03.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Ge, X., Cao, X., Song, X., Wang, Y., Si, Z., Zhao, Y., et al. (2020). Bioenergy generation and simultaneous nitrate and phosphorus removal in a pyrite-based constructed wetland-microbial fuel cell. Bioresour. Technol. 296:122350. doi: 10.1016/j.biortech.2019.122350

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, G., Ekama, G. A., Wang, Y., Dai, J., Biswal, B. K., Chen, G., et al. (2019). Advances in sulfur conversion-associated enhanced biological phosphorus removal in sulfate-rich wastewater treatment: a review. Bioresour. Technol. 285:121303. doi: 10.1016/j.biortech.2019.03.142

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, C., Liu, Q., Chen, C., Chen, F., Zhao, Y. K., Gao, L. F., et al. (2017). Elemental sulfur recovery and spatial distribution of functional bacteria and expressed genes under different carbon/nitrate/sulfide loadings in up-flow anaerobic sludge blanket reactors. J. Hazard. Mater. 324, 48–53. doi: 10.1016/j.jhazmat.2016.03.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, J., Wang, Z., Zhu, C., Ma, J., Zhang, X., and Wu, Z. (2014). Identification of microbial communities in open and closed circuit bioelectrochemical MBRs by high-throughput 454 pyrosequencing. PLoS One 9:e93842. doi: 10.1371/journal.pone.0093842

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, C., Yang, Q., Wang, D., Zhong, Y., Chen, F., Li, X., et al. (2017). Simultaneous perchlorate and nitrate removal coupled with electricity generation in autotrophic denitrifying biocathode microbial fuel cell. Chem. Eng. J. 308, 783–790. doi: 10.1016/j.cej.2016.09.121

CrossRef Full Text | Google Scholar

Li, X., Zhang, W., Liu, T., Chen, L., Chen, P., and Li, F. (2016). Changes in the composition and diversity of microbial communities during anaerobic nitrate reduction and Fe(II) oxidation at circumneutral pH in paddy soil. Soil Biol. Biochem. 94, 70–79. doi: 10.1016/j.soilbio.2015.11.013

CrossRef Full Text | Google Scholar

Liang, Z., Sun, J., Zhan, C., Wu, S., Zhang, L., and Jiang, F. (2020). Effects of sulfide on mixotrophic denitrification by Thauera-dominated denitrifying sludge. Environ. Sci. Water Res. Technol. 6, 1186–1195. doi: 10.1039/C9EW01014A

CrossRef Full Text | Google Scholar

Liu, C., Zhao, D., Yan, L., Wang, A., Gu, Y., and Lee, D. J. (2015). Elemental sulfur formation and nitrogen removal from wastewaters by autotrophic denitrifiers and anammox bacteria. Bioresour. Technol. 191, 332–336. doi: 10.1016/j.biortech.2015.05.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, F., Sun, L., Wan, J., Tang, A., Deng, M., and Wu, R. (2019). Organic matter and ammonia removal by a novel integrated process of constructed wetland and microbial fuel cells. RSC Adv. 9, 5384–5393. doi: 10.1039/C8RA10625H

CrossRef Full Text | Google Scholar

Magoč, T., and Salzberg, S. L. (2011). FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27, 2957–2963. doi: 10.1093/bioinformatics/btr507

PubMed Abstract | CrossRef Full Text | Google Scholar

Mao, Y., Xia, Y., and Zhang, T. (2013). Characterization of Thauera-dominated hydrogen-oxidizing autotrophic denitrifying microbial communities by using high-throughput sequencing. Bioresour. Technol. 128, 703–710. doi: 10.1016/j.biortech.2012.10.106

PubMed Abstract | CrossRef Full Text | Google Scholar

Oon, Y. L., Ong, S. A., Ho, L. N., Wong, Y. S., Dahalan, F. A., Oon, Y. S., et al. (2017). Role of macrophyte and effect of supplementary aeration in up-flow constructed wetland-microbial fuel cell for simultaneous wastewater treatment and energy recovery. Bioresour. Technol. 224, 265–275. doi: 10.1016/j.biortech.2016.10.079

PubMed Abstract | CrossRef Full Text | Google Scholar

Rogers, S. R., and Anderson, J. J. (1976). Role of iron deposition in Sphaerotilus discophorus. J. Bacteriol. 126, 264–271. doi: 10.1128/jb.126.1.264-271.1976

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, X., Zhang, J., Liu, D., Hu, Z., and Liu, H. (2018). Enhance performance of microbial fuel cell coupled surface flow constructed wetland by using submerged plants and enclosed anodes. Chem. Eng. J. 351, 312–318. doi: 10.1016/j.cej.2018.06.117

CrossRef Full Text | Google Scholar

Su, J. F., Li, G. Q., Huang, T. L., and Xue, L. (2020). The mixotrophic denitrification characteristics of Zoogloea sp. L2 accelerated by the redox mediator of 2-hydroxy-1,4-naphthoquinone. Bioresour. Technol. 311:123533. doi: 10.1016/j.biortech.2020.123533

PubMed Abstract | CrossRef Full Text | Google Scholar

Tan, X., Yang, Y. L., Li, X., Zhou, Z. W., Liu, C. J., Liu, Y. W., et al. (2020). Intensified nitrogen removal by heterotrophic nitrification aerobic denitrification bacteria in two pilot-scale tidal flow constructed wetlands: influence of influent C/N ratios and tidal strategies. Bioresour. Technol. 302:122803. doi: 10.1016/j.biortech.2020.122803

PubMed Abstract | CrossRef Full Text | Google Scholar

Villaseñor, J., Capilla, P., Rodrigo, M. A., Cañizares, P., and Fernández, F. J. (2013). Operation of a horizontal subsurface flow constructed wetland-microbial fuel cell treating wastewater under different organic loading rates. Water Res. 47, 6731–6738. doi: 10.1016/j.watres.2013.09.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Virdis, B., Read, S. T., Rabaey, K., Rozendal, R. A., Yuan, Z., and Keller, J. (2011). Biofilm stratification during simultaneous nitrification and denitrification (SND) at a biocathode. Bioresour. Technol. 102, 334–341. doi: 10.1016/j.biortech.2010.06.155

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Song, X., Wang, Y., Abayneh, B., Li, Y., Yan, D., et al. (2016). Nitrate removal and bioenergy production in constructed wetland coupled with microbial fuel cell: establishment of electrochemically active bacteria community on anode. Bioresour. Technol. 221, 358–365. doi: 10.1016/j.biortech.2016.09.054

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Tian, Y., Liu, H., Zhao, X., and Wu, Q. (2019). Effects of influent COD/TN ratio on nitrogen removal in integrated constructed wetland-microbial fuel cell systems. Bioresour. Technol. 271, 492–495. doi: 10.1016/j.biortech.2018.09.039

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Cai, Z., Sheng, S., Pan, F., Chen, F., and Fu, J. (2020). Comprehensive evaluation of substrate materials for contaminants removal in constructed wetlands. Sci. Total Environ. 701:134736. doi: 10.1016/j.scitotenv.2019.134736

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, D., George, A. E., Lu, H., Chui, H. K., Liu, W. T., Brdjanovic, D., et al. (2013). A new biological phosphorus removal process in association with sulfur cycle. Water Res. 47, 3057–3069. doi: 10.1016/j.watres.2013.03.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, K., Liu, T., Ma, C., Chang, B., Chen, R., and Wang, X. (2014). The role of Mn oxide doping in phosphate removal by Al-based bimetal oxides: adsorption behaviors and mechanisms. Environ. Sci. Pollut. R. Int. 21, 620–630. doi: 10.1007/s11356-013-1937-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, S., Lv, T., Lu, Q., Ajmal, Z., and Dong, R. (2017). Treatment of anaerobic digestate supernatant in microbial fuel cell coupled constructed wetlands: evaluation of nitrogen removal, electricity generation, and bacterial community response. Sci. Total Environ. 15, 339–346. doi: 10.1016/j.scitotenv.2016.11.138

PubMed Abstract | CrossRef Full Text | Google Scholar

Xia, X., Cao, X., Liang, P., Huang, X., Yang, S., and Zhao, G. (2010). Electricity generation from glucose by a Klebsiella sp. in microbial fuel cells. Appl. Microbiol. Biot. 87, 383–390. doi: 10.1007/s00253-010-2604-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiao, Y., Zheng, Y., Wu, S., Zhang, E. H., Chen, Z., Liang, P., et al. (2015). Pyrosequencing reveals a core community of anodic bacterial biofilms in bioelectrochemical systems from China. Front. Microbiol. 6:1410. doi: 10.3389/fmicb.2015.01410

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, D., Xiao, E., Xu, P., Lin, L., Zhou, Q., Xu, D., et al. (2017). Bacterial community and nitrate removal by simultaneous heterotrophic and autotrophic denitrification in a bioelectrochemically-assisted constructed wetland. Bioresour. Technol. 245, 939–999. doi: 10.1016/j.biortech.2017.09.045

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, F., Cao, F., Kong, Q., Zhou, L., Yuan, Q., Zhu, Y., et al. (2018). Electricity production and evolution of microbial community in the constructed wetland-microbial fuel cell. Chem. Eng. J. 339, 479–486. doi: 10.1016/j.cej.2018.02.003

CrossRef Full Text | Google Scholar

Xu, F., Ouyang, D. L., Rene, E. R., Ng, H. Y., Guo, L. L., Zhu, Y. J., et al. (2019). Electricity production enhancement in a constructed wetland-microbial fuel cell system for treating saline wastewater. Bioresour. Technol. 288:121462. doi: 10.1016/j.biortech.2019.121462

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, L., Zhao, Y., Fan, C., Fan, Z., and Zhao, F. (2017). First study to explore the feasibility of applying microbial fuel cells into constructed wetlands for COD monitoring. Bioresour. Technol. 243, 846–854. doi: 10.1016/j.biortech.2017.06.179

PubMed Abstract | CrossRef Full Text | Google Scholar

Yakar, A., Türe, C., Türker, O. C., Vymazal, J., and Saz, Ç (2018). Impacts of various filtration media on wastewater treatment and bioelectric production in up-flow constructed wetland combined with microbial fuel cell (UCW-MFC). Ecol. Eng. 117, 120–132. doi: 10.1016/j.ecoleng.2018.03.016

CrossRef Full Text | Google Scholar

Yang, N., Zhan, G., Wu, T., Zhang, Y., Jiang, Q., Li, D., et al. (2018). Effect of air-exposed biocathode on the performance of a Thauera-dominated membraneless single-chamber microbial fuel cell (SCMFC). J. Environ. Sci. 66, 216–224. doi: 10.1016/j.jes.2017.05.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, Y., Liu, J., Ma, B., Liu, Y., Wang, B., and Peng, Y. (2016). Improving municipal wastewater nitrogen and phosphorous removal by feeding sludge fermentation products to sequencing batch reactor (SBR). Bioresour. Technol. 222, 326–334. doi: 10.1016/j.biortech.2016.09.103

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, P., Peng, Y., Lu, J., Li, J., Chen, H., and Xiao, L. (2018). Microbial communities and functional genes of nitrogen cycling in an electrolysis augmented constructed wetland treating wastewater treatment plant effluent. Chemosphere 211, 25–33. doi: 10.1016/j.chemosphere.2018.07.067

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, R. C., Xu, X. J., Chen, C., Shao, B., Zhou, X., Yuan, Y., et al. (2019). Bioreactor performance and microbial community analysis of autotrophic denitrification under micro-aerobic condition. Sci. Total Environ. 647, 914–922. doi: 10.1016/j.scitotenv.2018.07.389

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhong, F., Wu, J., Dai, Y., Xiang, D., Cheng, S., and Ji, H. (2015). Performance evaluation of wastewater treatment using horizontal subsurface flow constructed wetlands optimized by micro-aeration and substrate selection. Water Sci. Technol. 71, 1317–1324. doi: 10.2166/wst.2015.090

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: substrate selection, hydraulic retention time, nitrogen transformation, functional bacterial groups, bioelectrical signals

Citation: Zhong F, Yu C, Chen Y, Wu X, Wu J, Liu G, Zhang J, Deng Z and Cheng S (2020) Nutrient Removal Process and Cathodic Microbial Community Composition in Integrated Vertical-Flow Constructed Wetland – Microbial Fuel Cells Filled With Different Substrates. Front. Microbiol. 11:1896. doi: 10.3389/fmicb.2020.01896

Received: 30 April 2020; Accepted: 20 July 2020;
Published: 05 August 2020.

Edited by:

Jochen A. Mueller, Helmholtz Centre for Environmental Research (UFZ), Germany

Reviewed by:

Aqiang Ding, Chongqing University, China
Sunil A. Patil, Indian Institute of Science Education and Research Mohali, India

Copyright © 2020 Zhong, Yu, Chen, Wu, Wu, Liu, Zhang, Deng and Cheng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Fei Zhong, Znpob25nQG50dS5lZHUuY24=; emloZW9AMTYzLmNvbQ==; Shuiping Cheng, c2hwY2hlbmdAdG9uZ2ppLmVkdS5jbg==

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.