Effect of Cathode Material and Its Size on the Abundance of Nitrogen Removal Functional Genes in Microcosms of Integrated Bioelectrochemical-Wetland Systems

Constructed wetland-microbial electrochemical snorkel (CW-MES) systems, which are short-circuited microbial fuel cells (MFC), have emerged as a novel tool for wastewater management, although the system mechanisms are insufficiently studied in process-based or environmental contexts. Based on quantitative polymerase chain reaction assays, we assessed the prevalence of different nitrogen removal processes for treating nitrate-rich waters with varying cathode materials (stainless steel, graphite felt, and copper) and sizes in the CW-MES systems and correlated them to the changes of N2O emissions. The nitrate and nitrite removal efficiencies were in range of 40% to 75% and over 98%, respectively. In response to the electrochemical manipulation, the abundances of most of the nitrogen-transforming microbial groups decreased in general. Graphite felt cathodes supported nitrifiers, but nirK-type denitrifiers were inhibited. Anaerobic ammonium oxidation (ANAMMOX) bacteria were less abundant in the electrochemically manipulated treatments compared to the controls. ANAMMOX and denitrification are the main nitrogen reducers in CW-MES systems. The treatments with 1:1 graphite felt, copper, plastic, and stainless-steel cathodes showed higher N2O emissions. nirSand nosZI-type denitrifiers are mainly responsible for producing and reducing N2O emissions, respectively. Hence, electrochemical manipulation supported dissimilatory nitrate reduction to ammonium (DNRA) microbes may play a crucial role in producing N2O in CW-MES systems.


Introduction
A fundamental change in the nitrogen (N) cycle has been happening since the beginning of 20th century as anthropogenic activities have considerably increased the availability of N in the environment through the legume cultivation, fossil fuel combustion, and the chemical fixation process [1]. These processes are extremely important for humans because of transport and agriculture,

Experimental Set-Up
The water and sediments used in the experiment were collected from the constructed wetland treating N rich agricultural diffuse pollution [28]. CW-MES system reactors were plastic vessels with the volume of 2 L with the height of 13 cm, diameter of the bottom was 10.5 cm and the diameter of top was 12 cm (Figure 1). The level of the sediment in the reactors was approximately 7 cm. The anodes were embedded in the sediment 3.5 cm above the bottom of the reactor. The material of anodes was graphite felt and their diameter was 10 cm and thickness 1.5 cm. The anode was connected to cathode via stainless steel wires, which were connected with copper tape. The cathode was resting on the sediment. For most of the time, the electrodes were connected in shortcut, in CW-MES configuration, and were disconnected only for the electrochemical studies. Two parameters of the cathodes were varied and studied, namely their size and the material of cathodes. We express the varying size of graphite felt cathodes (GFD-4,6-PA) as a ratio of geometrical surfaces between the electrode in the sediment and the one in water. Therefore, the ratio described as 1:1 had the same geometrical surface of anode and cathode, while in the ratios 1:2 and 1:4, the surface of cathode was respectively 2 and 4 times smaller than the surface of anode. Diameters were 9.4, 6.24, and 3.68 cm respectively, with the graphite felt thickness of 0.5 cm. In a series of experiments when the changing parameter was the cathode material, the following materials were used: graphite felt (GFD-4,6-PA) copper, plastic, and stainless steel (specifications of materials: in the Supplementary Information). In this case, the diameter of the cathode was always 8 cm. The graphite felt (GFD-4,6-PA) was pre-treated in an ultrasonic bath in ethanol/water solutions of decreasing concentrations.

Experimental Set-up
The water and sediments used in the experiment were collected from the constructed wetland treating N rich agricultural diffuse pollution [28]. CW-MES system reactors were plastic vessels with the volume of 2 L with the height of 13 cm, diameter of the bottom was 10.5 cm and the diameter of top was 12 cm (Figure 1). The level of the sediment in the reactors was approximately 7 cm. The anodes were embedded in the sediment 3.5 cm above the bottom of the reactor. The material of anodes was graphite felt and their diameter was 10 cm and thickness 1.5 cm. The anode was connected to cathode via stainless steel wires, which were connected with copper tape. The cathode was resting on the sediment. For most of the time, the electrodes were connected in shortcut, in CW-MES configuration, and were disconnected only for the electrochemical studies. Two parameters of the cathodes were varied and studied, namely their size and the material of cathodes. We express the varying size of graphite felt cathodes (GFD-4,6-PA) as a ratio of geometrical surfaces between the electrode in the sediment and the one in water. Therefore, the ratio described as 1:1 had the same geometrical surface of anode and cathode, while in the ratios 1:2 and 1:4, the surface of cathode was respectively 2 and 4 times smaller than the surface of anode. Diameters were 9.4, 6.24, and 3.68 cm respectively, with the graphite felt thickness of 0.5 cm. In a series of experiments when the changing parameter was the cathode material, the following materials were used: graphite felt (GFD-4,6-PA) copper, plastic, and stainless steel (specifications of materials: in the Supplementary Information). In this case, the diameter of the cathode was always 8 cm. The graphite felt (GFD-4,6-PA) was pretreated in an ultrasonic bath in ethanol/water solutions of decreasing concentrations. The different materials, i.e., graphite felt, plastic, stainless steel, and copper and varying ratios, i.e., 1:1 cathode, 1:2 cathode, and 1:4 cathode were used in order to study the effect of cathode material and its size on the abundance of the N transformation genes present in the bioelectrochemical systems. For the control, we used two systems where one system was without any electrodes and the second without cathode but had an anode embedded in the sediment. Every treatment was carried The different materials, i.e., graphite felt, plastic, stainless steel, and copper and varying ratios, i.e., 1:1 cathode, 1:2 cathode, and 1:4 cathode were used in order to study the effect of cathode material and its size on the abundance of the N transformation genes present in the bioelectrochemical systems. For the control, we used two systems where one system was without any electrodes and the second without cathode but had an anode embedded in the sediment. Every treatment was carried out in triplicate, and in all, 27 separate CW-MES systems were created. All other connections were made using copper wires.
The systems were operated for 45 days in a dark room at room temperature. To imitate the nitrate-rich waters, the reactors were added with 50 mg N-NO 3 -/l on days 9, 13, and 24. No external polarization was applied to the electrodes that were left connected in short circuit.

Sediment and Water Sampling
Sediment samples (ca 5-20 g each) from CW-MES reactors were collected in the beginning of the experiment (day 1), in the middle of the experiment (day 18) and at the end of the experiment (day 45) and stored at −20 • C. Dissolved nitrogen (DN), dissolved carbon (DC), and dissolved organic carbon (DOC) were determined in the sediment samples using a Vario TOC cube (Elementar Analysensysteme GmbG, Langenselbold, Germany), following EVS-EN 12260, EVS-EN 1484, and EVS-EN 1484 standards, respectively.
Water samples from the CW-MES reactors were collected every three days. Nitrate, nitrite, and sulfate concentrations were determined from the water samples using Ion chromatograph Metrohm 930 Compact IC Flex, the column used was Metrosep A Supp5 100/4.0 with the effluent composition of 1.0 mM NaHCO 3 and 3.2 mM Na 2 CO 3 .

N 2 O Sampling and Analysis
The CW-MES reactors were sampled once a week for N 2 O. For gas sampling, the vessel of the CW-MES reactor was closed tightly with the rubber seal and plastic cap, where the latter had connections for carrier gas inlet and gas sampling ( Figure 1). Helium 6.0 was used to purged to the headspace of the vessel. The silicon tubes were evacuated with syringe before sampling. Gas samples from the headspace of the vessel were collected into a previously evacuated (0.3 mbar) 30 mL glass vial in every 40 min during a 2-h session (sample time-points: 0, 40, 80, and 120 min) using a polypropylene syringe. The N 2 O concentration from the samples was analyzed using Shimadzu GC-2014 gas chromatography equipped with an electron capture detector (ECD) and Loftfield's autosampler [29].

DNA Extraction
The DNA was extracted from 0.25 g of the sediment samples using the PowerSoil DNA Isolation kit (MO BIO Laboratories Inc, Carlsbad, CA, USA) according to the described manufacturer's protocol. Homogenizer Precellys 24 (Berlin Technologies, Montigny-le-Bretonneux, France) was used for the homogenization of the samples at 5000 rpm for 20 s. The concentrations and the quality of the extracted DNA was analyzed using the spectrophotometer Infinite M200 (Tecan AG, Grodig, Austria). The extracted DNA samples were stored in −20 • C.
The quantification data were analyzed with RotorGene Series Software v. 2.0.2 (Qiagen, Hilden, Germany) and LinRegPCR program v. 2020.0 [31]. The abundance of each target gene was presented per gram of dry soil (copies/g dw). The abundance of total prokaryotic organisms was calculated by summing the bacterial and archaeal 16S rRNA gene abundances. The proportion of target genes and microbial groups in soil microbial community was estimated by normalizations against prokaryotic community. The used qPCR method is described in more detail by Espenberg et al. [6].

Electrochemical Study
Small pieces of the cathodes of each pilot were collected at day 20 and were then studied by cyclic voltammetry as the working electrode using a potentiostat (910 PSTAT mini Metrohm) in a three electrodes configuration cell. The reference electrode was Ag/AgCl (+0.197 V vs. SHE) and the scan rate was 1 mV s −1 . A stainless-steel rod was used as counter electrode. Two cycles were recorded and only the second one is shown. The cathode reactor samples were studied in a 0.1 M sodium chloride electrolyte prepared in pure water. The electrolyte was degassed with helium bubbling prior to the experiment in order to remove dissolved oxygen. The cyclic voltammograms were recorded before and after addition of 10 mM sodium nitrate (NaNO 3 ) to test for potential nitrate bio-electrochemical reduction by microbes at the surface of the cathode samples.
The cathode-and anode electrode polarizations were identified every other day in each pilot. For that a reference electrode was submerged into the sediment and the potential between the reference electrode and the cathode-or anode electrode was measured with a multimeter. During the experiment two different reference electrodes were used Ag/AgCl (+0.197 V vs. SHE) and saturated calomel electrode (SCE +0.244 V vs. SHE). To make data comparison easier, the polarizations were converted against SHE.

Statistical Analysis
We used R version 3.6.3 (R Core Team, 2020) to explore, harmonize, analyze, and visualize the data. To evaluate the significance of the differences between studied treatments with respect to physicochemical variables, gene parameters, and emission values, multivariate linear models were constructed, and models were tested using the ANOVA function in R package mvabund v. 4.1.3 [32]. Spearman's rank correlation coefficients were used to assess the relationships between target gene parameters and environmental factors. The p-values were also adjusted for the false discovery rate by the Benjamini-Hochberg method with significance at p < 0.05 and p < 0.01. Spearman correlation networks were visualized with Cytoscape v. 3.7.2 [33].

Physicochemical Properties, Nitrate Removal and N 2 O Flux
Sediment DC and DOC values of CW-MES showed significant difference between treatments (p < 0.05 and p < 0.01, respectively), although sediment DN concentration remained similar (p > 0.05) ( Figure 2). Throughout the study period DC, DOC, and DN values showed completely different temporal trend in Controls compared to other treatments as their concentration rather increased and was the highest in the last day of experiment (Supplementary Figure S1). The CW-MES were significantly different regarding the nitrate and sulfate concentration in the water (p < 0.05 and p < 0.01, respectively). The nitrite concentrations were very low and showed no significant difference between treatments (p > 0.05). Treatment 1:1 graphite felt increases nitrite for a while but not in a long-term perspective (Supplementary Figure S1). The water parameters for inorganic substrates, i.e., nitrate, nitrite and sulfate were also analyzed over the course of the experiment, nitrate and sulfate showed significant differences between treatments (p < 0.05 in both cases) ( Figure 3). The trend in nitrate, nitrite, and sulfate concentrations remains similar during the sediment sampling days and throughout the experiment. No statistical differences were discernible in N2O emission values directly before soil sampling between treatments (p > 0.05) ( Figure 3). The N2O emissions throughout the experimental duration displayed highest emissions for the treatments with cathodes 1:1 graphite felt, copper, plastic, and stainless steel.  The water parameters for inorganic substrates, i.e., nitrate, nitrite and sulfate were also analyzed over the course of the experiment, nitrate and sulfate showed significant differences between treatments (p < 0.05 in both cases) ( Figure 3). The trend in nitrate, nitrite, and sulfate concentrations remains similar during the sediment sampling days and throughout the experiment. No statistical differences were discernible in N 2 O emission values directly before soil sampling between treatments (p > 0.05) ( Figure 3). The N 2 O emissions throughout the experimental duration displayed highest emissions for the treatments with cathodes 1:1 graphite felt, copper, plastic, and stainless steel. The water parameters for inorganic substrates, i.e., nitrate, nitrite and sulfate were also analyzed over the course of the experiment, nitrate and sulfate showed significant differences between treatments (p < 0.05 in both cases) ( Figure 3). The trend in nitrate, nitrite, and sulfate concentrations remains similar during the sediment sampling days and throughout the experiment. No statistical differences were discernible in N2O emission values directly before soil sampling between treatments (p > 0.05) ( Figure 3). The N2O emissions throughout the experimental duration displayed highest emissions for the treatments with cathodes 1:1 graphite felt, copper, plastic, and stainless steel.  Nitrate removal efficiency was highly variable in different treatments (Table 1). Treatments with plastic cathodes showed the highest nitrate removal efficiency followed by control and no cathode systems. All other systems showed slightly lower nitrate removal efficiency, which was less than 50%. Nitrite removal efficiency on the other hand was very high (over 98%) in all the studied systems. qPCR results showed that the total abundance of archaea differed significantly between the sediments of treatments (p < 0.05), although that was not the case for the total abundance of bacteria (p > 0.05) (Figure 4). Latter still showed somewhat higher abundance in the two controls (control and no cathode). No significant differences were not recognized in proportions of bacteria and archaea in prokaryotic communities over the CW-MES sediments (p > 0.05 in both cases) (Supplementary Figure S2). The bacterial 16S rRNA genes were found in more abundance than the archaeal 16SrRNA genes in the sediment samples. In addition, it seems that 1:1 graphite felt decreases the total gene abundances of bacteria and archaea as well as functional genes compared to other treatments, although the same trend was not noticed in the target gene proportions of the prokaryotic communities. Throughout the whole experiment period in most of the treatments, the temporal dynamics of the bacterial 16S rRNA gene parameters showed somewhat decreasing trend and the archaeal 16S rRNA gene parameters showed somewhat increasing trend, except for the control in both cases ( Figure 5; Supplementary Figure S3).
nirS, nirK, nosZI, nosZII, bacterial and archaeal amoA, ANAMMOX-specific 16S rRNA genes were detected in all sediment samples, whereas nrfA genes were not detected in sediments of one parallel of No Cathode, Graphite Felt, and 1:2 Graphite Felt treatments (Figure 4; Supplementary Figure S2). Based on average values of the gene parameters, the abundances of the functional gene copy numbers and the functional gene proportions in the prokaryotic communities were declining in order as follows: nirS, nirK, nosZII, nosZI, ANAMMOX-specific 16S rRNA, archaeal amoA, bacterial amoA, nrfA. In the all sediments, the archaeal amoA genes abundance exceeded the bacterial amoA abundance by ca one order of magnitude and nirS and nosZII genes were found to be around two-fold more abundant than nirK and nosZI genes, respectively.
The total abundance of archaeal amoA, nirK, nosZII, and ANAMMOX-specific 16S rRNA gene copy numbers varied significantly between the sediments of the treatments (p < 0.05, p < 0.01, p < 0.05 and p < 0.01, respectively) ( Figure 4). In addition, nirS proportions in the prokaryotic communities showed significant difference in the treatments' sediments (p < 0.05) (Supplementary Figure S2). The total abundance of almost all studied functional genes was roughly higher in the two controls (control and no cathode) than other treatments. The only exception was determined for nrfA genes (both the total abundances and the proportions in the prokaryotic community), which were in the same level or showed even higher values in the treatments with electrochemical manipulation compared to the controls. All denitrification genes (nirS, nirK, nosZII, nosZI) showed lower abundances in treatments of 1:1 Graphite Felt and Graphite Felt compared to others. A decrease in the ANAMMOX-specific 16S rRNA gene total abundances and proportions in the prokaryotic community was reported regarding the treatments with electrochemical manipulation.     The temporal trend followed by all treatments were similar in terms for the different functional genes in most cases ( Figure 5; Supplementary Figure S3). In general, nir, nosZI, nrfA and archaeal amoA genes have shown a rise in their copy numbers and also in their proportions in prokaryotic community, when to compare the beginning and end of the experiment. The nosZII genes showed a fall in both abundances and proportions (in prokaryotic community) throughout the experiment. The substantial decline in described gene parameters was observed of the ANAMMOX-specific 16S rRNA gene during the full experiment period. Some opposite trends were noticed between treatments in temporal scale regarding bacterial amoA gene parameters. In case of treatments with the graphite felt cathode, clear opposite trends are apparent for bacterial amoA gene and archaeal amoA genes, while nirK genes seemed to be inhibited.

Relationships between Target Genes, Physicochemical Parameters and N 2 O Flux
The data analysis detected several statistically significant relationships between target gene parameters and environmental factors, while all treatments, except two controls (control and no cathode), were analyzed together ( Figure 6). The sediment DC and DN concentrations showed positive correlations with the abundances of bacterial 16S rRNA, nirS, nosZI, nosZII, and nrfA genes (p < 0.05 in all cases), while negative correlations were observed between DC and DN values and nirK and ANAMMOX-specific 16S rRNA genes proportions in the prokaryotic community (p < 0.05 in all cases). Similarly, the sediment DOC concentration was related to the abundances of bacterial 16S rRNA, nirS, and nrfA genes and ANAMMOX-specific 16S rRNA genes proportions in the prokaryotic community (p < 0.05 in all cases). Positive correlation between water nitrate concentration and archaeal amoA proportion was revealed (p < 0.05), while this chemical factor had negative relationship with archaeal 16S rRNA abundance and nirK proportion (p < 0.05 in both cases). The water nitrite content was related positively to nirS proportion (p < 0.05) and negatively to nosZII and ANAMMOX-specific 16S rRNA genes proportions (p < 0.05 in both cases). Only negative relationships between water sulfate concentration and nosZII and nrfA gene parameters were revealed by a correlation analysis (p < 0.01 and p < 0.05, respectively). The proportion of nirK and ANAMMOX-specific 16S rRNA genes correlated negatively with N2O emission from the CW-MES systems (p < 0.05 in both cases). In addition, the abundance of nirS, nosZI, and nosZII genes and the proportion of nirS and nosZI genes were correlated positively with N2O emission from the CW-MES systems (p < 0.05 in all cases). Several significant relationships more appeared, but they were not confirmed by the Benjamini-Hochberg correction ( Figure S4). Nevertheless, these relationships may have also indications to The proportion of nirK and ANAMMOX-specific 16S rRNA genes correlated negatively with N 2 O emission from the CW-MES systems (p < 0.05 in both cases). In addition, the abundance of nirS, nosZI, and nosZII genes and the proportion of nirS and nosZI genes were correlated positively with N 2 O emission from the CW-MES systems (p < 0.05 in all cases). Several significant relationships more appeared, but they were not confirmed by the Benjamini-Hochberg correction ( Figure S4). Nevertheless, these relationships may have also indications to important processes in the CW-MES systems, although they should be taken into account with reservations. Sediment DC, DOC and DN contents were all positively related to abundances of bacterial and archaeal amoA genes (p < 0.05 in all cases). Water nitrate concentration showed different relationships with nirS, nirK, nrfA and bacterial 16S rRNA genes (p < 0.05 in all cases). Water nitrite concentration was related to different genes (nirS, nosZI, nosZII, bacterial amoA, and archaeal 16S rRNA genes) quite differently (p < 0.05 in all cases). Water sulfate content was still negatively correlated with most of the studied different gene parameters (p < 0.05 in all cases). Functional gene nrfA abundance showed positive correlation with N 2 O from the CW-MES systems (p < 0.05). The proportion of nosZI in the prokaryotic community was positively correlated with one-time measurements of N 2 O (p < 0.05), but negatively with the average N 2 O emission values throughout the experiment (p < 0.05).

Electrochemical Measurements
Cyclic voltammetry study on cathode reactor samples showed a reduction system for two out of three replicates equipped with graphite felt cathodes (Figure 7 and Supplementary Figure S5). No similar reduction peaks were obtained with the other cathode materials. For the first pilot (graphite felt 1), a reduction peak is visible at −0.3 V vs. Ag/AgCl under helium saturation before nitrate addition and the peak current of this system increases after nitrate addition. The second pilot (graphite felt 2) shows an even clearer difference with no reduction at all in the case of helium saturation and no nitrate addition, with a clear reduction peak appearing after nitrate supplementation at the same potential of  Average polarizations of the cathode electrodes were +0.136, −0.02, and −0.166 V vs. SHE for carbon felt, copper and stainless-steel electrodes, respectively. And the average polarizations of the anode electrodes were −0.03, −0.08, and −228 V vs. SHE for carbon felt, copper, and stainless-steel electrodes, respectively.

Discussion
The capability of microorganisms to produce electrical current or to accept and use electrons continues to be a high-potential approach and strategy for developing sustainable water treatment technologies [11,34]. However, bioelectrochemical systems have present a challenge to face for several reasons, e.g., no conserved functional genes exist to specifically target electrochemically active microorganisms [11]. Thus, one of the best options is to examine mechanisms underlying microbialmediated in CW-MES and other electrochemical systems through process-based genes from the microbiological perspective. Denitrification has received the greatest attention regarding Average polarizations of the cathode electrodes were +0.136, −0.02, and −0.166 V vs. SHE for carbon felt, copper and stainless-steel electrodes, respectively. And the average polarizations of the anode electrodes were −0.03, −0.08, and −228 V vs. SHE for carbon felt, copper, and stainless-steel electrodes, respectively.

Discussion
The capability of microorganisms to produce electrical current or to accept and use electrons continues to be a high-potential approach and strategy for developing sustainable water treatment technologies [11,34]. However, bioelectrochemical systems have present a challenge to face for several reasons, e.g., no conserved functional genes exist to specifically target electrochemically active microorganisms [11]. Thus, one of the best options is to examine mechanisms underlying microbial-mediated in CW-MES and other electrochemical systems through process-based genes from the microbiological perspective. Denitrification has received the greatest attention regarding electrochemical systems [23][24][25], although there is a rising interest in other and newly discovered processes (e.g., ANAMMOX and DNRA), because they can be crucial in the N dynamics [11].
The results of this study showed nitrate removal efficiencies in range of 40-75% and the highest removal efficiency was detected in treatments with plastic cathodes. All the treatments containing graphite felt cathodes removed nitrate from the CW-MES systems in a quite similar level. In treatments with graphite felt cathodes, nitrate reduction was also demonstrated by cyclic voltammetry. The cathode materials were tested independently after the experiments for their potential nitrate reduction ability: only two graphite felt cathodes showed an electrochemical response to nitrate addition at −0.3 V, which suggests that the biofilm on these cathodes may harbor electroactive denitrifying microorganisms [35,36]. Nitrite removal efficiencies were over 98% in all studied treatments of CW-MES systems.
Our results showed that electrochemically manipulated treatments all in all rather decreased the abundance of bacteria compared to controls, which may indicate the community shift because of the effect of the treatment. The electrotrophic activity among bacteria and archaea seems to be not very widely spread. However, the diversity of electroactive bacterial, archaeal, and eukaryotic microorganisms is constantly increasing due to the discovery of new electroactive microorganisms [11]. However, MFCs are not definable only by electroactive organisms [37][38][39] the microbial cooperation is ubiquitous in nature and helps microorganisms to overcome biological and thermodynamic limitations in different environments [40]. Based on our results, archaea account for ca 5-20% of all prokaryotes in CW-MES systems and change in archaeal abundance between the controls and other treatments seemed to be smaller compared to bacteria. Archaea are widely present in many environments and they possess unique features [41], and maybe, they are also quite resilient to electrochemical manipulation or even this manipulation enhance the abundance of archaea.
In addition, we can notice that the genes in the controls, i.e., control with no electrodes, no cathode control with just the anode electrode, and plastic as the cathode material, i.e., non-conductive electrode, displayed the highest number of abundances for many of the target N transformation genes. The fall in the target gene copy numbers in some other treatments can be attributed to the electrons passing through the system hampering the count of the microbes, thereby reducing the observed copy numbers in the conductive systems. Furthermore, the larger the size of the electrodes more is the conductivity and the reduction rate but when bio-electrochemical systems come into picture the major catalysts are the microbes and therefore, they should not burn out due to higher current through the systems. For example, some versions of graphite felt cathodes treatment (1:1 graphite felt and the systems in which graphite felt was used as a cathode material) lessen the abundances of bacteria, archaea as well as functional genes (especially denitrification genes) compared to other treatments, but this tendency was not seen in the target gene proportions of the prokaryotic community.
We found that the balance between bacterial amoA and archaeal amoA genes was in favor of latter genes in CW-MES systems. In both engineered and natural environments, the domination and activity of bacterial or archaeal ammonia oxidizers depend on mainly on factors like substrate concentration, temperature, light, pH, and oxygen [42]. Nitrification seemed to be suppressed to some extent in different treatments with cathodes as controls showed higher abundance of archaeal amoA possessing microbes. Treatment with graphite felt supported amoA-harbouring archaea or bacteria depending on physicochemical conditions, whereas nirK-type denitrifiers were somewhat inhibited.
The abundance of nir was significantly higher than nosZ in CW-MES systems, although the nirSand nosZII-type denitrifiers were more abundant in prokaryotic community, compared to denitrifiers possessing nirK and nosZI genes, respectively. In general, temporal dynamics of nirSand nosZI-type denitrifiers were very similar and of nirK and nosZII-type denitrifiers showed almost opposite trend during the experiment. This can be attributed to the facts that there is higher frequency of co-occurrence of nosZ with nirS than with nirK (only 10-30% of the nirK-type denitrifiers also carry the nosZ gene) and most of the nosZI-types possess nirK or nirS (ca half of nosZII-harbouring microbes lack nir genes) [43]. Our results indicated that the abundance of all denitrification genes was smaller in the treatments with cathodes than controls, which may mean a similar trend of the community shift, as discussed above about bacteria. Different treatments altered significantly the abundances of nirS, nirK, and nosZII-type denitrifiers.
Electrochemical activity has been recently confirmed in highly enriched cultures of the ANAMMOX bacteria ('Candidatus Brocadia sinica' and 'Candidatus Scalindua sp.'), in spite of the fact that their electroactive potential and characteristics have been really scarcely studied [11]. In our CW-MES systems, the abundance of ANAMMOX bacteria was different between treatments and they were less abundant in the treatments with cathodes than controls. Accordingly, it may indicate their reluctance to electrochemical actions, despite some of their possible electroactive features or it would be also explained due to the presence of low amounts of nitrite in the water, which must have curbed the microbes involved in carrying out ANAMMOX processes.
In our CW-MES systems, microorganisms conducting DNRA seemed not to be affected too much of the electrochemical manipulation, and some of the treatments showed even higher abundances and proportions of DNRA microbes compared to the controls. The rise in the nrfA gene copy numbers may be attributed to their diverse metabolic pathways of DNRA microbes, which allow them to thrive at varying environments [44]. The DNRA retain the amount of nitrogen in the system by converting the nitrite into ammonium and thus, the rise in nrfA possessing organisms could have possibly led to ammonium generation, thereby influencing the rise in the amoA-harbouring microbes.
Our results indicated that gene parameters relationships with sediment DC, DOC and DN and water nitrite and nitrate concentrations revealed a competition regarding those essential nutrients. Rivalry was mainly between denitrification and ANAMMOX and less with nitrification. Water sulfate content inhibited the abundance of all the major N-cycle genes.
From an environmental sustainability point of view, N 2 O emissions should be also estimated and considered in designs of CW-MFC systems [27]. In our CW-MES systems, the treatments with 1:1 graphite felt, copper, plastic, and stainless-steel cathodes showed higher N 2 O emissions compared to others throughout the study period. N 2 O emissions could be mainly attributed to nirS-type denitrifiers, although there were some indications that DNRA microbes may also play a crucial role in producing N 2 O in CW-MES systems. nosZI-type denitrifiers were for the most part responsible for reducing N 2 O emissions.

Conclusions
The results of this study show high nitrate and nitrite removal efficiencies in CW-MES systems. In response to the electrochemical manipulation, the abundances of bacterial, archaeal, and most of the N-transforming microbial groups decreased. Treatment with graphite felt cathodes supported nitrifying archaea or bacteria depending on the physicochemical conditions, but nirK-type denitrifiers were inhibited. ANAMMOX bacteria were less abundant in the treatments with cathodes than controls and their electroactive potential is in ambiguity. However, ANAMMOX bacteria and denitrifiers are the main N reducers in CW-MES systems. The treatments with 1:1 graphite felt, copper, plastic, and stainless-steel cathodes showed higher N 2 O emissions compared to others throughout the study period. nirSand nosZI-type denitrifiers are mainly responsible for producing and reducing N 2 O emissions, respectively. There is some evidence that DNRA microbes are assisted by electrochemical manipulation and may also play a crucial role in producing N 2 O in CW-MES systems.
Supplementary Materials: The following are available online at http://www.mdpi.com/2571-8789/4/3/47/s1, Figure S1. Temporal dynamics of different sediment and water parameters during the experiment (each n = 3). Abbreviations: DN-dissolved nitrogen, DC-dissolved carbon, DOC-dissolved organic carbon, Figure S2. Box plots of target gene copy numbers proportion in the prokaryotic community in nine studied treatments (each n = 6). Statistically significant differences between groups considering different genes are shown with F and p values determined by multivariate analysis. The central line is the median, the edges of the box are the 25th and 75th percentiles, the whiskers represent the 95% confidence interval, and grey dots indicate outliers, Figure S3. Temporal dynamics of different gene proportions in the prokaryotic community during the whole experiment (for each n = 3), Figure S4. Spearman correlation network (p < 0.05) for each target gene parameter and environmental factor pair (n = 42) which were left out after the Benjamini-Hochberg correction. Green is used to represent positive relationships, and red indicates negative relationships (exact R values are shown on the graph); yellow stands for gene abundances, grey for soil physicochemical variables and pink for gaseous parameters. Abbreviations: DN-dissolved nitrogen, DC-dissolved carbon, DOC-dissolved organic carbon, period-average values throughout the experiment, bact-bacterial 16S rRNA gene abundance, arch-archaeal 16S rRNA gene abundance, prok-prokaryotic community (total bacterial and archaeal 16S rRNA gene proportion), amoA_B-bacterial amoA gene abundance, amoA_A-archaeal amoA gene abundance, N 2 O-N 2 O emission, Figure S5. Different cyclic voltammograms of carbon felt i.e., graphite felt before and after Sodium Nitrate addition in the electrolyte under helium saturation and ambient air. Possible nitrate reduction is showed by blue arrows. Funding: This study was supported by the Estonian Research Council (PUTJD715, IUT2-16, PRG352 and PRG676); the EU through the European Regional Development Fund (Centre of Excellence EcolChange, Estonia) and by the European Structural and Investment Funds. The financial support from The French National Research Agency (ANR-17-CE04-0004) is gratefully acknowledged. The financial support from the International mobility support for PhD students-DrEAM (University of Lorraine) is gratefully acknowledged.