The Effects of Plants on Pollutant Removal, Clogging, and Bacterial Community Structure in Palm Mulch-Based Vertical Flow Constructed Wetlands

: In this study, the effects of plants on the performance and bacterial community structure of palm mulch-based vertical ﬂow constructed wetlands was studied. The wetlands were built in August 2013; one of them was planted with Canna indica and Xanthosoma sp., and the other one was not planted and used as a control. The experimental period started in September 2014 and ﬁnished in June 2015. The inﬂuent was domestic wastewater, and the average hydraulic surface loading was 208 L/m 2 d, and those of COD, BOD, and TSS were 77, 57, and 19 g/m 2 d, respectively. Although the bed without plants initially performed better, the ﬁrst symptoms of clogging appeared in December 2014, and then, its performance started to fail. Afterwards, the wetland with plants provided better removals. The terminal restriction fragment length polymorphism (T-RFLP) analysis of Enterococci and Escherichia coli in the efﬂuents suggests that a reduction in their biodiversity was caused by the presence of the plants. Thus, it can be concluded that the plants helped achieve better removals, delay clogging, and reduce Enterococci and E. coli biodiversity in the efﬂuents.


Introduction
Constructed wetlands (CWs) are particularly well suited for wastewater treatment in small communities, i.e., those with a <2.000 p.e. (person equivalent) because of their low cost, easy maintenance, high treatment efficiency, and visual appeal [1]. Similarly, CWs are environmentally friendly and a real alternative for the treatment of different types of wastewaters [2]. CWs are open, shallow reactors composed basically of an impervious layer, a mineral substrate (usually gravel or sand), helophytes, water, and the associated microbes. According to the water flow, CWs can be classified as horizontal flow (HF) or vertical flow (VF) and surface or subsurface flow. One of the main constraints of CWs is their high surface area, usually 3-5 m 2 /p.e. VF CWs demand less surface and are more efficient than HFs because of the higher aeration of the substrates [3]. However, VF CWs usually include smaller-sized substrates, such as sand, to achieve a good surface influent distribution

Location of Constructed Wetland and Sampling
The studied CW was designed to treat raw wastewater from a group of 4 households located in a rural zone of the Island of Gran Canaria, Spain. The households have a permanent population of 5 people, which is increased to 15 people during holidays and weekends.
Wastewater was discharged directly into a partially clogged, non-waterproofed, 36 m 3 seepage pit ( Figure 1) with an approximate hydraulic retention time of 1.5 days. Two submergible, 700-W pumps, which were timer-controlled to function for 1 min each hour, were installed. Influent pulses were of approximately 15-20 L. The pit was supposed to provide partial sedimentation and homogenization. Thus, the influent can be considered to be partially decanted domestic wastewater.
were of approximately 15-20 L. The pit was supposed to provide partial sedimentation and homogenization. Thus, the influent can be considered to be partially decanted domestic wastewater.
The CW was placed at 650 m above sea level. The average temperature is 18 °C, with January being the coldest month with an average temperature of 12 °C and August the warmest month with an average temperature of 23 °C. The average relative air moisture is 80% and the predominant wind direction is NNE. The average pluviosity is 485 mm/year, with the rainy season being between October and March.
The two VF CWs were constructed with 1.080-L plastic, cubic-like recipients. A layer of washed gravel with a height of 10-15 cm and an average diameter of 1 cm was deposited at the bottom of the recipients. The organic substrate was placed directly above the gravel and consisted of approximately 1-m-high mulch made of dry branches of Canarian palm tree (Phoenix canariensis). Additionally, two flow meters and valves were placed at the inlet of each reactor to measure and control the inflow. The influent was distributed on the surface of the VF CWs by means of 5 perforated tubes with a 12-mm diameter.
In order to value the contribution of the plants, one VF CW was planted with Canna indica and Xanthosoma sp. and compared with another VF CW (the control), which was not planted. These species of plant are characterized by large roots and big tubers, respectively. The analyses were started one year after the construction of the VF CWs, because the mulch substrate had been reduced by degradation and compaction.
Chemical and microbiological analyses were completed in unfiltered and homogenized samples of the influent and the effluents of the planted (1p) and unplanted (1np) VF CWs.

Water Analysis
Sampling was performed in the morning at 08:00 h, coinciding with when the pumps started a working period. The influent sample was taken during pumping from the influent outlet. The volume sampled was 1 L for each sampling point. The sampling recipients were 1-L plastic bottles for standard parameters, sterilized 50-mL glass bottles for E. coli and total E. coli enumeration, and sterilized 1-L amber glass bottles for DNA studies. The CW was placed at 650 m above sea level. The average temperature is 18 • C, with January being the coldest month with an average temperature of 12 • C and August the warmest month with an average temperature of 23 • C. The average relative air moisture is 80% and the predominant wind direction is NNE. The average pluviosity is 485 mm/year, with the rainy season being between October and March.
The two VF CWs were constructed with 1.080-L plastic, cubic-like recipients. A layer of washed gravel with a height of 10-15 cm and an average diameter of 1 cm was deposited at the bottom of the recipients. The organic substrate was placed directly above the gravel and consisted of approximately 1-m-high mulch made of dry branches of Canarian palm tree (Phoenix canariensis). Additionally, two flow meters and valves were placed at the inlet of each reactor to measure and control the inflow. The influent was distributed on the surface of the VF CWs by means of 5 perforated tubes with a 12-mm diameter.
In order to value the contribution of the plants, one VF CW was planted with Canna indica and Xanthosoma sp. and compared with another VF CW (the control), which was not planted. These species of plant are characterized by large roots and big tubers, respectively. The analyses were started one year after the construction of the VF CWs, because the mulch substrate had been reduced by degradation and compaction.
Chemical and microbiological analyses were completed in unfiltered and homogenized samples of the influent and the effluents of the planted (1p) and unplanted (1np) VF CWs.

Water Analysis
Sampling was performed in the morning at 08:00 h, coinciding with when the pumps started a working period. The influent sample was taken during pumping from the influent outlet. The volume sampled was 1 L for each sampling point. The sampling recipients were 1-L plastic bottles for standard parameters, sterilized 50-mL glass bottles for E. coli and total E. coli enumeration, and sterilized 1-L amber glass bottles for DNA studies.
The water quality parameters were measured according to standard methods (APHA, 2005). BOD 5  The concentrations of Escherichia coli and total coliforms (TC) were determined by the membrane filter method with incubation at 37 • C for 24 h with a chromogenic agar (Panreac, Spain). E. coli displayed dark blue to violet colonies, while TC were assumed to be the sum of the E. coli colonies plus the salmon-red colonies.

Molecular Method
In order to analyze molecularly the bacterial community structure, a terminal restriction fragment length polymorphism (T-RFLP) analysis was carried out in the effluents of both CWs. The analyses were performed in December 2014, when the initial clogging symptoms were observed in 1np (light clogging), and in May 2015, when the clogging in 1 np was severe (severe clogging).

DNA Isolation
First, 1 liter of water of each sample was centrifuged, and pellets pooled to extract genomic DNA. The DNA isolation was performed according to the CTAB method described by Murray and Thompson [28]. In detail, CTAB was made of CTAB 0.1% (w/v), PVPP 0.1% (w/v), TRIS-HCl pH 8.6, SDS 10, EDTA 0.5M pH 8, NaCl 4M, and β-Mercaptoethanol 2% (v/v). To obtain bacterial samples, 800-µg pellets were introduced into microtubes. Then, an extraction buffer (800 µL) was added to a microtube containing the source of the DNA template. The sample was held for 1 h in a bath at 65 • C, and mixed gently by inversion approximately every 20 min. Later, 800 µL of CIA (chloroform: isoamyl alcohol, 24:1 v/v) was added and centrifuged for 20 min at 3000 rpm in a Beckman centrifuge.
Successive washes in CIA and centrifugations were carried out until the supernatant became whitish. To continue, 2/3 of isopropyl alcohol at −20 • C was added, and a centrifugation for 30 min at 15,000× g was carried out. Afterwards, the isopropyl alcohol was removed, and 20 µL of ethanol (80%) was added. Finally, the samples were centrifuged for 5 min at 15,000× g. The supernatant was then discarded, and the pellets were resuspended in 15 µL of H 2 O and stored at 4 • C until used.
The yield and purity of genomic DNA were calculated from the A260/A280 ratio measured using a NanoDrop spectrophotometer (ThermoFisher Scientific, Waltham, MA, USA). For each VF CW type and date, an assay was repeated three times with two replicates of each one.
The amplifications were performed in a thermocycler MyCyclerTM (Biorad, Hercules, CA, USA), with an initial denaturation at 96 • C for 4 min followed by 40 cycles at 94 • C for 1 min, 45 • C for 1 min, and 72 • C for 5 min. The final elongation was carried out at 72 • C for 10 min. The PCR-amplification products were checked by agarose (1%) electrophoresis at 75 V, resulting in gene fragments of 733 pb and 876 pb for Enterococcus sp. and E. coli, respectively. The amplifications were carried out by triplicate.
PCR-amplification products were separately purified using a Wizard ® SV Gel and PCR Clean-Up System (Promega, Madison WI, USA). Purified products from E. coli and Enterococcus sp. were checked by electrophoresis on a 2% agarose gel. The concentrations were estimated as previously mentioned.
T-RFLP analysis was carried out based on the previous enzyme digestion assay. The NEBcutter v2.0 (New England BioLabs) software program was used to estimate a range of terminal restriction fragments (T-RFs) of varying sizes in silico. The restriction enzyme MseI (New England BioLabs) was selected in silico according to the yielded T-RFs of the respective PCR amplicons for E. coli and Enterococcus sp.
Enzyme digestion was performed with 100 ng of purified PCR amplicon for 3 h at 37 • C and inactivated at 65 • C for 20 min. Amplicon was then cut with 0.5 U of the restriction enzyme MseI (recognition site: T'TAA) according to the manufacturer's instructions (BioLabs, New England BioLabs). Restriction digestion was carried out 3 times for the E. coli and Enterococcus sp., respectively.
Sequencing was carried out on an ABI 3730 capillary sequencer (Life Technologies) by Secugen (Madrid, Spain), and the resulting sequences were analyzed using GeneMarker v1.85 (SoftGenetics). Optimizing was performed with different dilutions according to Secugen.

Chemical and Microbiological Analysis
Average values of concentrations, surface loadings, and removals were compared by means of ANOVA if the data were homoscedastic (Bartlett test) and normally distributed (Shapiro-Wilk test). If these conditions were not met, the Kruskal-Wallis non-parametric test was used. In all cases, a significance level of 95% (p-value > 0.05) was utilized.

T-RFLP Analysis
The T-RFLP fragment sequences were obtained in .fsa format and GeneMarker software (SoftGenetics) was used for data scoring. The analyses of the polymorphic bands were only scored when they were highly reproducible and according to the peak height. Thus, small T-RFs ranging from 0 to 100 bp were excluded from the analysis, as they were considered to be PCR artifacts. Moreover, those T-RFs with intensities lower than 1% florescence, regarded as background interference, were also excluded from the matrices. Next, relative peak areas were recalculated according to the removal of artifactual peaks. The relative abundance of T-RFs was determined by calculating the ratio between the peak height of each peak and the total peak height of all the peaks within 1 sample and assuming that unambiguous T-RFs ranged in size from 100 to 850 bp.
For a comparative analysis of T-RFLP profiles, alignment was carried out through T-Align software with a default value threshold of 0.5 [31]. This software allowed for the compiling of peak profiles of different electropherograms, which corresponded to separate sampling days and sampling sites. In addition, the resulting consensus profile (henceforth overall consensus electropherogram) valued the presence or absence of T-RFs through a binary matrix by comparing with the other consensus profile of each of the other sites. For the purposes of qualitative valuation, the number of peaks in the consensus electropherograms was interpreted to be the operational taxonomic units (OTUs) [27], although non-sequencing was carried out. Statistical comparisons of a number of T-RFs were performed using R software [32]. A one-way ANOVA followed by post hoc Tukey HSD and Dunnett T3 tests was used to detect significant differences (p ≤ 0.1) between 1p and 1np under light and severe clogging. Additionally, Venn diagrams were performed with Venn Diagram Plotter (Integrative Omics, Pacific Northwest National Laboratory) [33].

Features of the Influent
As mentioned above, wastewater from the households was discharged directly into the seepage pit. The pit acted as a primary settler, making the influent more stable and reducing the concentration of the total suspended solids (TSS) as compared with other systems without a primary settler [16]. Table 1 shows the chemical and microbiological parameters of the influent. The number of data points changes for each parameter because only analytically valid data and non-outliers are included. Other parameters such as Na, electrical conductivity, and Sodium Absorption Ratio (SAR)were analyzed less frequently with the goal to obtain the agronomic quality of the effluents. As can be observed, the biodegradability (BOD/COD) ratio was 71% and turbidity was 521 NTU. With these traits, the influent can be considered to range from medium to high strength for organic matter (COD and BOD), low strength for TSS, and medium strength for nutrients N and P [34].

Average Surface Loading Rates and Removals
The hydraulic loading rate (HLR) and organic loading rates (OLR) are represented in Table 2. The average BOD-LRs applied in the present study were 51 and 70 g/m 2 d, for 1np and 1p, respectively. These values fall in the upper range of those found in the literature. For instance, the OLR applied to the French System's first stage ranged between 40 and 50 g BOD/m 2 d [5].
The average HLR of 1p (233 L/m 2 d) was significantly higher (p: 0.026) than that of 1np (189 L/m 2 d). The same result was obtained for the BOD-LR (1p: 70 g/m 2 d, 1np: 51 g/m 2 d, p: 0.0233) and ammonia-N (1p: 9 g/m 2 d, 1np: 7 g/m 2 d, p: 0.0212). Non-significant differences were reported for the COD-LR (1p: 90 g/m 2 d and 1np: 70 g/m 2 d, p: 0.0965) and TSS-LR (1p: 22 g/m 2 d and 1np: 18 g/m 2 d, p: 0.129). Considering that in the cases of dissimilarity, the p-values were close to 0.05 Sustainability 2019, 11, 632 7 of 18 (HLR, COD, and ammonia-N) and that there was no difference for COD and TSS, it can be concluded that both VF CWs received almost similar surface loadings and their differences in removal can be mainly ascribed to the presence of plants. Although for some parameters, 1p received higher surface loadings, its performance was equal to or better than 1np ( Table 3). The removals were slightly higher for 1p, although not significantly (p > 0.1) for COD (1p: 58%, 1np: 50%), BOD (1p: 62%, 1np: 59%), TSS (1p: 77%, 1np: 72%), E. coli (1p: 79%, 1np: 73%), and TC (1p: 82%, 1np: 74%). The exception was turbidity removal (1p: 90%, 1np: 73%, p: 0.005) for which the presence of plants improved clearly. This result is of interest, because according to Spanish legislation, turbidity is a basic parameter for the reuse of treated wastewater [35]. The removals of ammonia-N, phosphate-P, and sulfate were negative in all cases. The increased concentrations of the ions can only come from the substrate, because it is an organic material and from the accumulated sludge. As mentioned above, high surface loadings were applied, and no rest periods were allowed. Additionally, the reactors were very efficient in removing TSS, which were retained as sludge, mainly on the surface of the CWs. The negative removals of ammonia-N, phosphate, and sulfate ions are more likely to come from the leaching or degradation of the accumulated sludge, rather than from the mulch, which had been in operation for 1 year. In other works (e.g., [16]) with palm mulch, we have observed good physical and chemical stability of the substrate, both in vertical and horizontal CWs, once it became stabilized, after 2-3 months in operation. However, under the stringent conditions imposed to the reactors, a further degradation of the substrate cannot be rejected.
It can be thought that because the substrate was mainly of an organic nature, the release of these nutrients can be caused by its decomposition. In fact, Saeed and Sun [36] observed that HF CWs, which employed wood mulch and gravel-mulch media, released organic matter, phosphorus, and TSS. This behavior was not observed in the VF CWs tested in our study. Additionally, in the present research, the substrate height was reduced by 40% by compaction and degradation during the first few months of operation. However, the analyses were initiated 1 year after construction. Thus, it seems more probable that the released ions came from the fast mineralization of the particulate matter retained on the surface of the VF CWs under the appropriate conditions of temperature, moisture, and aeration.

Performance Evolution of the Wetlands
The evolution of parameters such as COD, BOD, turbidity, TSS, ammonia-N, E. coli, and TC are represented for the influent and the effluents of 1p and 1np ( Figure 2). Missing points have been removed after the statistical analysis of outliers. research, the substrate height was reduced by 40% by compaction and degradation during the first few months of operation. However, the analyses were initiated 1 year after construction. Thus, it seems more probable that the released ions came from the fast mineralization of the particulate matter retained on the surface of the VF CWs under the appropriate conditions of temperature, moisture, and aeration.
(a)   Figure 3 shows the timeline of removals. During December 2014, the first clogging symptoms were observed in 1np. After that, the relative performance of the wetlands changed dramatically. Table 4 shows the average removals of chemical and microbiological parameters before and after clogging. Before clogging, 1np performed better than 1p with respect to turbidity (1np: 95% and 1p: 90%; p: 0.028), E. coli (1p: 59%, 1np: 97%, p: 0.007), and TC (1p: 82%, 1np:98 %, p: 0.001). Nevertheless, the removals of BOD (1p: 83%, 1np: 85%), COD (1p: 64%, 1np: 62%), and TSS (1p: 84%, 1np: 80%) were similar (p > 0.288). After clogging, 1np performance was significantly damaged (p < 0.006) for COD, BOD, turbidity, E. coli, and TC, while that of 1p was improved or remained unaltered.   The better performance of 1np before clogging can be explained by the stronger media compaction in the absence of plants and the consequent increase in the hydraulic retention time (HRT). In 1p, plant roots would be able to keep the media unclogged and, as a consequence, reduce the HRT and removal efficiency. In fact, it was observed that the effluent of 1np exited more slowly than that of 1p. Torrens et al. [37] observed a similar behavior in planted VF CWs, because the effluent from the first batches flowed through the beds more quickly than through unplanted ones. The authors claimed that the plant rhizomes could create preferential pathways. The preventing effect of plant roots regarding clogging has been discussed for a long time. Some studies have suggested that plants can make the substrate more porous [38], while others, such as Teixeira et al., [39] found a reduction in the porosity of the medium of HF CWs caused by the development of roots of Vetiver and Tifton 85 grasses. Yet, they considered that the low root volume, only 3.07% and 4.11% of the total pore space, could not have a great influence on clogging. However, Hua et al. [40] claimed that the role of plants varied throughout the clogging process. In the early stages, plant roots restricted water flow, while in the later stages, growing roots opened new pore spaces in the substrate. Our results are in agreement with those of Hua et al. [40] for the later stage but are the opposite for the early stages. This illustrates the complexity of the clogging process and the role of plant rhizomes. The better performance of 1np before clogging can be explained by the stronger media compaction in the absence of plants and the consequent increase in the hydraulic retention time (HRT). In 1p, plant roots would be able to keep the media unclogged and, as a consequence, reduce the HRT and removal efficiency. In fact, it was observed that the effluent of 1np exited more slowly than that of 1p. Torrens et al. [37] observed a similar behavior in planted VF CWs, because the effluent from the first batches flowed through the beds more quickly than through unplanted ones. The authors claimed that the plant rhizomes could create preferential pathways. The preventing effect of plant roots regarding clogging has been discussed for a long time. Some studies have suggested that plants can make the substrate more porous [38], while others, such as Teixeira et al., [39] found a reduction in the porosity of the medium of HF CWs caused by the development of roots of Vetiver and Tifton 85 grasses. Yet, they considered that the low root volume, only 3.07% and 4.11% of the total pore space, could not have a great influence on clogging. However, Hua et al. [40] claimed that the role of plants varied throughout the clogging process. In the early stages, plant roots restricted water flow, while in the later stages, growing roots opened new pore spaces in the substrate. Our results are in agreement with those of Hua et al. [40] for the later stage but are the opposite for the early stages. This illustrates the complexity of the clogging process and the role of plant rhizomes. The fact that 1np became clogged in such a short time, after about two years in operation, was unexpected. Additionally, in the case of 1p, a negative trend in performance along time can be observed, particularly for BOD (Figure 3). To understand these results, the following must be considered: (i) The mulch employed was not fully stable. Although the analyses were made 1 year after the construction of the beds, the results indicate that the substrate was not fully stabilized. The palm branch is basically composed of a central rachis on which leaflets and spines are inserted. Once the branch is dead, the leaflets degrade relatively fast and release color, organic matter, and ions. Hence, their presence in palm mulch for CW substrate should be minimized. (ii) The average surface loadings were relatively high (BOD: 59-62 g/m 2 d, TSS: 18-22 g/m 2 d, Table 2). (iii) The reactors were in continuous operation without any rest period. In a survey of 169 full-scale VF CWs treating raw domestic wastewater, Paing et al. [5] claimed that the nominal BOD load calculated on the first stage was generally 40-50 g/m 2 d, which is a similar value to that applied in this study. However, the success of the French system is partially justified by the good aeration of the substrate. This was achieved by different measures that include the application of rest periods, allowing the sludge accumulated on the surface to become mineralized.
The reduction in the performance of 1np caused by clogging underlines the importance of the substrate aeration. Regarding pathogen removal, Headley et al. [41] observed that aeration improved E. coli removal in horizontal flow wetlands. The better performance of 1p could be caused by the accumulation of sludge on the bed surface and the consequent HRT increase. At the same time, the presence of the roots would keep the substrate unclogged.
These results show that the effect of plants in mulch-based VF CWs was not univocal. In the early stages of the CWs, the presence of the plants reduced the efficiency of the removal of pollutants by augmenting the substrate porosity with their roots. However, in the long term, they helped to retard clogging and, as a consequence, to achieve better removals.

Analysis of T-RFLP
Terminal restriction fragment length polymorphism (T-RFLP) is a rapid, highly reproducible, and robust molecular tool for the study of bacterial community structure [27,42]. It is known that bacteria communities are affected by several factors, such as physicochemical properties, climate conditions, and the presence of plants [43,44]. Hence, bacterial communities of 1p and 1np were analyzed via T-RFLP in order to ascertain the effects of the presence of plants.
Bacterial diversity, as estimated by the number of T-RFs, showed reproducible patterns across the MseI enzyme. In addition, the size distribution of the T-RFs generated by the restriction enzyme was consistent with that of the T-RFs derived from in silico digestion.
Interestingly, shifts in the sizes of the T-RFs occurred as a result of two different factors: the presence of plants and the progressive clogging of the substrate, as can be inferred from the analysis of the operational taxonomic units (OTUs). The two VF CWs shared 31 T-RFs for E. coli and 19 for Enterococcus sp. (Figure 4). The shared T-RFs indicate the presence of a common bacteria population regardless of the presence of plants ( Figure 4). Furthermore, the size of the T-RFs ranged from 117 to 810 bp for E. coli and from 115 to 714 bp for Enterococcus sp.
In addition, 1np showed a significant increase (p < 0.1) in T-RFs compared with that of 1p for both E. coli and Enterococcus sp. (Figure 5).
was consistent with that of the T-RFs derived from in silico digestion.
Interestingly, shifts in the sizes of the T-RFs occurred as a result of two different factors: the presence of plants and the progressive clogging of the substrate, as can be inferred from the analysis of the operational taxonomic units (OTUs). The two VF CWs shared 31 T-RFs for E. coli and 19 for Enterococcus sp. (Figure 4). The shared T-RFs indicate the presence of a common bacteria population regardless of the presence of plants (Figure 4). Furthermore, the size of the T-RFs ranged from 117 to 810 bp for E. coli and from 115 to 714 bp for Enterococcus sp.  In addition, the results from the light and severe clogging periods showed that 1p reported a significant diminution (p < 0.1) in the T-RFs of E. coli over time (22.22% in light clogging versus 17.68% severe clogging). This behavior was also shown for 1np (38.89% light clogging versus 21.21% severe clogging; Figure 5a). This difference was more significant (p < 0.01) when relative abundance of T-RFs in 1np (38.89%) were compared to that in 1p (22.22%) at light clogging (Figure 5a). This result may infer that 1np was slightly clogged in December 2014 and continued to clog over time. In this sense, the chemical parameters sustain the T-RFs results (Table 4).
Additionally, the same trend in relative abundance of T-RFs was reported for Enterococcus sp., although the lower biodiversity of 1np under severe clogging was unexpected (Figure 5b). Hua et al. [26] claimed that clogging reduced bacterial diversity. Hence, an approach to explain low diversity in 1np was attempted by comparing shared T-RFs in 1p and 1np over time ( Figure 6). Venn diagrams showed that 1np shared a major number of T-RFs between two sampled periods (1.8% light clogging versus 7.8% severe clogging). The presence of low shared-T-RFs, namely OTUs, in 1p would highlight the role of plant roots in combination with palm mulch substrate at retaining FIB bacteria. Kadam et al. [27] argued that the diminution of T-RFs in the planted VF CWs was committed to the diminution of genetic features, which could be correlated with the presence of root-associated microbiota. In addition, we assumed that root presence and microfauna might affect 1p performance.
The results obtained in this work indicate that plants increase the substrate porosity, but the results can be quite different depending on the level of substrate clogging. During the early stage of the experimental period, when the substrate was far from being clogged, the effect was negative as the CW with plants demonstrated worse pollutant removal. In this case, an increase in porosity led to shorter HRT. The sludge accumulated on the substrate would help to increase the HRT and thus the removals. However, later, when the sludge accumulation becomes excessive and clogging starts to develop, the consequences were very positive, because the plants were able to notably delay clogging by increasing the substrate porosity. The strong reduction of the palm mulch during the first year of operation indicates the importance of the substrate selection. The application of rest periods is highly recommended to delay clogging and improve the pollutant removals and the accumulated sludge mineralization. Similarly, the constructed wetland performance was optimized by the presence of plants as the reduction in the number of OTUs demonstrated. light clogging versus 7.8% severe clogging). The presence of low shared-T-RFs, namely OTUs, in 1p would highlight the role of plant roots in combination with palm mulch substrate at retaining FIB bacteria. Kadam et al. [27] argued that the diminution of T-RFs in the planted VF CWs was committed to the diminution of genetic features, which could be correlated with the presence of root-associated microbiota. In addition, we assumed that root presence and microfauna might affect 1p performance. The results obtained in this work indicate that plants increase the substrate porosity, but the results can be quite different depending on the level of substrate clogging. During the early stage of the experimental period, when the substrate was far from being clogged, the effect was negative as the CW with plants demonstrated worse pollutant removal. In this case, an increase in porosity led to shorter HRT. The sludge accumulated on the substrate would help to increase the HRT and thus the removals. However, later, when the sludge accumulation becomes excessive and clogging starts to develop, the consequences were very positive, because the plants were able to notably delay clogging by increasing the substrate porosity.
The strong reduction of the palm mulch during the first year of operation indicates the importance of the substrate selection. The application of rest periods is highly recommended to delay clogging and improve the pollutant removals and the accumulated sludge mineralization. Similarly, the constructed wetland performance was optimized by the presence of plants as the reduction in the number of OTUs demonstrated.

Conflicts of Interest:
The authors declare no conflicts of interest.