Anaerobic Conversion of Saline Phenol-Containing Wastewater Under Thermophilic Conditions in a Membrane Bioreactor

Closing water loops in chemical industries result in hot and highly saline residual streams, often characterized by high strength and the presence of refractory or toxic compounds. These streams are attractive for anaerobic technologies, provided the chemical compounds are biodegradable. However, under such harsh conditions, effective biomass immobilization is difficult, limiting the use of the commonly applied sludge bed reactors. In this study, we assessed the long-term phenol conversion capacity of a lab-scale anaerobic membrane bioreactor (AnMBR) operated at 55°C, and high salinity (18 gNa+.L–1). Over 388 days, bioreactor performance and microbial community dynamics were monitored using specific methanogenic activity (SMA) assays, phenol conversion rate assays, volatile fatty acids permeate characterization and Illumina MiSeq analysis of 16S rRNA gene sequences. Phenol accumulation to concentrations exceeding 600 mgPh.L–1 in the reactor significantly reduced methanogenesis at different phases of operation, while applying a phenol volumetric loading rate of 0.12 gPh.L–1.d–1. Stable AnMBR reactor performance could be attained by applying a sludge phenol loading rate of about 20 mgPh.gVSS–1.d–1. In situ maximum phenol conversion rates of 21.3 mgPh.gVSS–1.d–1 were achieved, whereas conversion rates of 32.8 mgPh.gVSS–1.d–1 were assessed in ex situ batch tests at the end of the operation. The absence of caproate as intermediate inferred that the phenol conversion pathway likely occurred via carboxylation to benzoate. Strikingly, the hydrogenotrophic SMA of 0.34 gCOD-CH4.gVSS–1.d–1 of the AnMBR biomass significantly exceeded the acetotrophic SMA, which only reached 0.15 gCOD-CH4.gVSS–1.d–1. Our results indicated that during the course of the experiment, acetate conversion gradually changed from acetoclastic methanogenesis to acetate oxidation coupled to hydrogenotrophic methanogenesis. Correspondingly, hydrogenotrophic methanogens of the class Methanomicrobia, together with Synergistia, Thermotogae, and Clostridia classes, dominated the microbial community and were enriched during the three phases of operation, while the aceticlastic Methanosaeta species remarkably decreased. Our findings clearly showed that highly saline phenolic wastewaters could be satisfactorily treated in a thermophilic AnMBR and that the specific phenol conversion capacity was limiting the treatment process. The possibility of efficient chemical wastewater treatment under the challenging studied conditions would represent a major breakthrough for the widespread application of AnMBR technology.


INTRODUCTION
Phenols are major contaminants found in wastewaters of several chemical industries, which are often discharged at high temperatures (Busca et al., 2008;Rosenkranz et al., 2013;Wang et al., 2017). Additionally, closing water loops in the chemical sector often result in concentrated, highly saline wastewaters, which increases the complexity of the produced wastewater . Despite the existing physicochemical processes applied for phenol removal, i.e., membrane distillation, pervaporation, adsorption, extraction, nanofiltration, reverse osmosis, and oxidation processes (wet air, electrochemical, ozonation, UV/H 2 O 2 , Fenton) (Villegas et al., 2016;Raza et al., 2019); biological treatment processes are preferred due to its cost-effectiveness. In this regard, anaerobic treatment offers the advantages of minimal energy requirement, low sludge production, and the conversion of organic pollutants into energy-rich biogas. Under both saline and high-temperature conditions, effective biomass immobilization becomes cumbersome, constraining the application of anaerobic sludge bed systems to treat these wastewaters (Dereli et al., 2012;van Lier et al., 2015). Moreover, the phenol degrading capacity of methanogenic consortia is generally very low and is expected to develop only slowly. Therefore, combining anaerobic treatment with membrane assisted biomass separation is an attractive option when phenol conversion is required at high salinity, and thermophilic conditions (Lin et al., 2013;Muñoz Sierra et al., 2018b).
When chemical wastewaters are at high temperatures, direct thermophilic treatment becomes of interest because it reduces the need for cooling the wastewater. Particularly when process water reclamation is envisaged, maintaining the temperature reduces the overall energy requirement. Despite its potentials (van Lier, 2008), thus far, full-scale anaerobic membrane bioreactors (AnMBRs) are not applied for high-temperature chemical wastewater treatment (Duncan et al., 2017). Only a few previous studies have shown the potential of thermophilic conditions for treating phenolic compounds in continuous reactors (Wang et al., 2011;Ramakrishnan and Surampalli, 2013). Wang et al. (2011) compared UASB reactors under mesophilic and thermophilic conditions and concluded that thermophilic anaerobic digestion improves about 30% the degradation of phenolic compounds. Likewise, Ramakrishnan and Surampalli (2013) suggested that thermophilic conversion of phenolic wastewater in an anaerobic hybrid reactor is superior to mesophilic in terms of methane yield, effluent quality, and stability.
Conversely, other studies have also found drawbacks of treating phenol containing wastewater under thermophilic conditions. Fang et al. (2006) indicated that the phenol conversion rate at 55 • C in a UASB reactor was significantly lower than under mesophilic conditions. Furthermore, Levén and Schnürer (2005) concluded that phenolics are mineralized to methane and carbon dioxide under mesophilic conditions, whereas under thermophilic conditions, only benzoic acid is degraded. Muñoz Sierra et al. (2018b) suggested that under mesophilic and hyper-mesophilic conditions (42-45 • C), the phenol conversion capacity of an AnMBR at high salinity is more stable compared to thermophilic conditions. However, because the operation at 55 • C in that study was carried out only during a short-term, it remains unclear whether or not a stable phenol degrading methanogenic consortium may develop. Therefore, this study aims to determine the maximum conversion capacity of a laboratory-scale AnMBR during a long-term operation at 55 • C, treating phenol-containing wastewater at 18 gNa +. L −1 . Moreover, the microbial community activity and structure in response to increasing phenol loading rates along with three phases of operation were evaluated.

Experimental Set-Up and Operation
The experiments were performed by using a 6.5 L laboratoryscale AnMBR reactor, equipped with an ultra-filtration (UF) sidestream membrane module (Figure 1). A tubular polyvinylidene fluoride membrane (X-flow, Pentair, Netherlands) with 5.2 mm inner diameter, 0.64 m length and 30 nm nominal pore size was used. The reactor was equipped with feed, recycle, and effluent pumps with 4-20 mA variable speed (120U/DV, 220Du, Watson-Marlow, Netherlands), pH, and temperature sensors (Memosens, Endress & Hauser, Germany), and a biogas meter (Milligas Counter MGC-1 PMMA, Ritter, Germany). Transmembrane pressure (TMP) was measured by using three pressure sensors (AE Sensors ATM, Netherlands). The temperature of the jacketed reactor was controlled under thermophilic conditions by a thermostatic water bath (Tamson Instruments, Netherlands). The entire set-up was controlled by a programmable logic controller (PLC) connected to a PC with LabVIEW software (version 15.0.1f1, National Instruments, United States).
The AnMBR was operated at 55.0 ± 0.8 • C for 388 days. During this time, the reactor was fed with synthetic phenolic wastewater with phenol concentration between 50 and 800 mgPh . L −1 and a sodium concentration of 18 gNa +. L −1 . The experiment was divided into three phases (I, II, and III). In all of the phases, the phenol concentration was increased step-wise to determine the maximum conversion capacity and the attainable phenol loading rate of the AnMBR. The applied total organic loading rates (OLRs) were between 2.0 and 4.0 gCOD . L −1. d −1 during all phases ( Table 1). The OLRs were calculated as OLR = influent COD (g.L −1 ) * Flow rate (L.d −1 )/AnMBR volume (L). The applied flow rate was 1.0 L.d −1 . An average solids retention time (SRT) of about 120 ± 13 days was maintained. The average SRT was calculated periodically as SRT av = X (g VSS in the AnMBR) av /X removed net (gVSS.d −1 ) , where X removed net resulted from the biomass removed.day −1 (sampling for tests) and the biomass returned.day −1 (after some tests).
The AnMBR was completely mixed applying a turnover of 170 times.d −1 .The membrane unit was operated at a cross-flow velocity of 0.65 m s −1 . A cyclic membrane filtration operation was carried out, consisting of 500 s filtration and 30 s backwash. Backwash was done by reversing the permeate pump flow. An operational flux of 4.0 L . m −2. h −1 was applied as a result of the experimental settings. The permeate flow was controlled with the variable speed of the effluent and influent and pumps, and it was regularly double-checked manually. An average TMP of 177 ± 92 mbar and a total membrane filtration resistance of 9.77 × 10 12 [1/m] could be maintained during the AnMBR long-term operation.

Volatile Fatty Acids
Prior analysis, 10 mL of the AnMBR permeate samples was filtrated over 0.45 µm filter paper. The filtrated liquid was diluted with pentanol (300 mg.L −1 ). 10 µL of formic acid (purity > 99%) was added into the final 1.5 mL vials. Volatile fatty acids (VFAs) were measured by gas chromatography (GC) using an Agilent 19091F-112, 25 m × 320 µm × 0.5 µm column and an FID detector (Agilent 7890A, United States). The sample injection volume was 1 µL. Helium was used as carrier gas with a total flow rate of 67 mL/min and a split ratio of 25:1. The GC oven temperature was programmed to increase from 80 to 180 • C in 10.5 min. The temperatures of the injector and detector were 80 and 240 • C, respectively.

Permeate Characterization
Phenol concentrations were measured using high-performance liquid chromatography HPLC LC-20AT (Shimadzu, Japan) equipped with a 4.6 mm reversed-phase C18 column (Phenomenex, Netherlands) and a UV detector at a wavelength of 280 nm. The mobile phase used was 25% (v/v) acetonitrile at a flow rate of 0.95 mL . min −1 . The column oven was set at 30 • C. Fast phenol measurements were also carried out by Merck -Spectroquant Phenol cell kits by using a spectrophotometer NOVA60 (Merck, Germany). Hach Lange kits were used to measure chemical oxygen demand (COD). The COD was measured using a VIS -spectrophotometer (DR3900, Hach Lange, Germany) making proper dilutions to minimize interference by high chloride concentrations, without compromising the accuracy of the measurement.

Anaerobic Phenol Conversion Rates
Batch tests were conducted in triplicate at the end of phase II to assess the phenol conversion recovery after AnMBR performance perturbation. The volatile suspended solid (VSS) were analyzed according to standard protocols using the lowest possible sample volume (APHA, 2005). Biomass samples were taken with a 150 mL syringe and transferred to 500 mL Schott glass bottles. The bottles were filled to a volume of 400 mL with AnMBR biomass (0.68 g VSS), and a medium containing acetate (3.1-4.6 g . L −1 ), phenol (60-109 mgPh . L −1 ), 6 mL . L −1 macro-and 0.6 mL . L −1 micronutrients solution, 10 mM phosphate buffer solution, and 18 gNa +. L −1 . Three consecutive feedings of the medium were applied. Initial COD and phenol concentrations varied between 3.4-5.1 gCOD . L −1 and 60-109 mgPh . L −1 , respectively. Temperature and mixing were controlled in an orbital incubator shaker (New Brunswick Biological Shakers Innova 44/44R, United States) at 55 • C and 120 rpm respectively. Periodically, liquid samples were taken, and phenol and COD concentrations were measured. Phenol conversion rates [mgPh . gVSS −1. d −1 ] were calculated by using the slope of the phenol concentration vs time curve in each bottle. After the batch tests were finished, the supernatant was removed and biomass was returned to the AnMBR. Similarly, at the end of the AnMBR operation at day 388, biomass samples were taken, and phenol conversion batch tests were carried out with initial phenol concentrations of 40 and 60 mgPh . L −1 .

Specific Acetotrophic and Hydrogenotrophic Methanogenic Activity
Specific acetotrophic methanogenic activity (SMA) tests were performed in triplicate using an automated methane potential test system (AMPTS, Bioprocess Control, Sweden). All the SMA tests were carried out at 55 • C, following the method described by Spanjers and Vanrolleghem (2016). For the hydrogenotrophic methanogenic activity, 250 mL Schott glass bottles were filled with biomass (0.57 g VSS) and medium (6 and 0.6 mL . L −1 macro-and micro-nutrients solution, respectively and 10 mM phosphate buffer solution at pH 7.0) to a liquid volume of 200 mL. The gas-phase of the bottles was exchanged by using a gas exchange board (G.R Instruments B.V, Netherlands) with a gas mixture of 80% CO 2 and 20% H 2 to an end pressure of 0.5 bar during 5 continuous automated cycles to ensure the complete absence of oxygen. Bottles were incubated in an orbital shaker (New Brunswick Biological Shakers Innova 44/44R, United States) at 55 • C and 120 rpm for 10 days, and biogas production was calculated using the exact headspace volume and the drop in headspace pressure versus time. The headspace pressure was measured as described by Coates et al. (1996) using a pressure transducer. The methane content of the biogas was analyzed by using a gas chromatograph 7890A (GC) (Agilent Technologies, United States) equipped with a front thermal conductivity detector (TCD). The temperature of the oven was 45 • C for 6 min, then 25 • C/min to 100 • C. The temperatures of the front inlet, and a front detector were both 200 • C.

Microbial Community Analysis
Biomass samples were taken from the AnMBR on days 88, 241, and 376 to evaluate the microbial community. The DNA extraction was performed from 0.5 g of biomass by using the DNeasy UltraClean Microbial Kit (Qiagen, Hilden, Germany). Agarose gel electrophoresis and Qubit3.0 DNA detection (Qubit dsDNA HS Assay Kit, Life Technologies, United States) were used for quality and quantity control of the DNA. The amplification of the 16S rRNA gene (V3-V4 region) was performed and followed by high throughput sequencing using the MiSeq Illumina platform (BaseClear, Leiden, Netherlands) using the primers 341F (5 -CCTACGGGNGGCWGCAG-3 ) and 785R (5 -GACTACHVGGGTATCTAATCC-3 ). The Illumina fastq reads (2 × 250) were processed in the QIIME2 pipeline (2018.7) (Bolyen et al., 2019). Reads were quality filtered, denoised, and the amplicon sequences variants (ASVs) were resolved with the DADA2 plugin (Callahan et al., 2016), removing chimeras with the "consensus" method. The taxonomic classification of the representative sequences of ASVs was performed with the "classify-consensus-vsearch" plugin (Rognes et al., 2016) using the SILVA (132) database as a reference. The representative sequences were aligned with the MAFFT algorithm (Katoh and Standley, 2013), and a phylogenetic tree was constructed with FastTree (Price et al., 2010). The feature table and tree were exported to the R environment. Differential abundance analyses between reactor operation phases were performed with the DESeq2 library (Love et al., 2014). The abundance and the tree were visualized with the phyloseq library (McMurdie and Holmes, 2013). ASVs with differential abundances within the operational phases were analyzed with BLAST against the refseq RNA database to identify the closest related species. The sequences reported in this paper have been deposited at ENA under the study accession number PRJEB38467.

Thermophilic AnMBR Performance
Influent and effluent phenol and COD concentrations during the long-term thermophilic AnMBR operation are shown in Figure 2. During days 0-96 in phase I, the effluent COD concentrations were in the range of 2.0-10.0 gCOD . L −1 , and the corresponding removal efficiencies were between 59.0 and 92.3% ( Figure 2B) at an average OLR of 3.8 gCOD . L −1. d −1 . The increase in the phenol loading rate ( Table 1) from 0.01 to 0.09 gPh L −1 d −1 concomitantly increased the phenol removal efficiency from 54 to 95% (Figure 2A). At the end of phase I, the phenol concentration in the reactor reached about 738 mgPh . L −1 (Figure 2A).

Volatile Fatty Acids (VFAs) Spectrum
Throughout the entire thermophilic operation, VFAs were detected in the AnMBR permeate (Figure 3), indicating a limiting methanogenic conversion capacity. The effluent COD mainly consisted of acetate (0.02-9.6 g . L −1 ), which peaked at almost 9.6 g . L −1 between 133 and 144 days when phenol accumulation occurred. Concomitantly, the butyrate concentration increased to 616 mg . L −1 in this period. In phase II, high concentrations of acetate (5.2 g . L −1 ) and to a lesser extent butyrate (295 mg . L −1 ) were found at day 260 when phenol again accumulated after an increase in the phenol loading rate. Propionate was most of the time present in the reactor effluent with an average concentration of 129 ± 57 mg . L −1 . The valerate concentration increased to 254 and 156 mg . L −1 when the reactor performance deteriorated in phases I and II, respectively. In phase III, on day 350, an increase in all VFAs was observed when reactor phenol concentration increased to 124 mgPh . L −1 .

Conversion Rates and Methanogenic Activity
Specific Methanogenic Activity and Phenol Conversion Rates SMA tests using acetate as the substrate showed a drop in the methanogenic activity of the phenol-degrading biomass in phase I, meanwhile, an increase in the specific phenol conversion rates was observed ( Table 2). During phase II the lowest observed SMA was 0.04 ± 0.02 gCOD-CH 4 . gVSS −1. d −1 , following phenol accumulation at the end of phase I. In phase III when the influent phenol concentration was decreased to 0.2 gPh . L −1 , the SMA of the phenol-degrading biomass increased to 0.13 gCOD-CH 4 . gVSS −1. d −1 , which was similar to the value observed at the beginning of phase I. Table 2 also shows the calculated maximum in-reactor phenol conversion rates in the different phases. At the end of phase I, the phenol conversion rate had increased from an initial value of 1.3 to 21.3 ± 0.2 mgPh . gVSS −1. d −1 . During the recovery periods of perturbation, i.e., days 156-198, the phenol conversion rate was 5.9 ± 0.3 mgPh . gVSS −1. d −1 . From days 220 to 262, the phenol loading rate increased until an average of 16.9 ± 0.6 mgPh . gVSS −1. d −1 . In phase III, the phenol conversion rate decreased to the range 2.6 ± 0.1-7.6 ± 1.6 mgPh . gVSS −1. d −1 .

Ex situ Phenol Conversion Rate After Reactor Perturbation
After reactor perturbation at the end of phase II, the phenol conversion rate was assessed in a batch test. Three different feedings were applied with different initial phenol concentrations (see Table 3). After the first feed, the phenol conversion rate was calculated as 4.0 ± 1.4 mgPh . gVSS −1. d −1 after 7 days of incubation. The phenol conversion rate increased to 9.6 ± 2.6 and 10.5 ± 3.3 mgPh . gVSS −1. d −1 after the second and third consecutive feeding respectively.

Specific Methanogenic Activity and Phenol Conversion Rates at the End of the Long-Term Operation
Since low acetate-fed SMA values at the end of the long-term thermophilic operation period were observed of 0.13 ± 0.10 gCOD-CH 4 . gVSS −1. d −1 , SMA tests were performed with both acetate and hydrogen as electron donor. Likewise, the phenol conversion rate was measured in batch test with an initial phenol concentration of 40 and 60 mgPh . L −1 (  phenol conversion rate of 32.8 ± 0.5 mgPh . gVSS −1. d −1 was found, applying an initial phenol concentration of 60 mgPh . L −1 .

Microbial Community Structure Analysis
The microbial community analysis revealed a total of 141 ASVs with differential abundance across samples. Figure 4 shows the genera from both bacteria and archaea domains with main population changes in relative abundance at the different phases of the thermophilic AnMBR at 18 gNa + L −1 . Petrotoga (Thermotogae class) was enriched along with the long-term operation up to 21.1% in phase III. Thermovirga (Synergisitia class) increased from 8.0 to 14.1% from phase I to II and then decreased to 7.95% relative abundance in phase III. Acetomicrobium also belonging to phylum Synergistetes FIGURE 4 | Heatmap of the genera from bacteria and archaea domains in the thermophilic AnMBR that were positive after differential abundance analysis (DESeq2) among the three phases of reactor operation. The color scale ranges from 0 to 22% relative abundance.
(see Supplementary Figure 1) increased from 0.35% during phase I to 3.93% during phase III. A 100% similarity was found with the specie Acetomicrobium hydrogeniformans sp. (see Supplementary Table 1). Syntrophobacter affiliated to the class Deltaproteobacteria increased from 0.1 to 0.5% in relative abundance from phase I to phase II; however it was not detected in phase III. A hit of 100% similarity was found for the halophilic bacteria Syntrophobacter sulfatireducens sp. Similarly, the relative abundance of WS6 bacterium decreased during the long-term operation. In the case of the Clostridia class, the major abundance decrease was observed with the genera Caldicoprobacter and Tepidimicrobium, whereas an increase was observed for Syntrophaceticus (3.3%), Pelotomaculum (1.2%), and Proteiniclasticum (2.4%). Also, Corynebacterium and Enterococcus genera belonging to Actinobacteria and Bacilli class increased significantly to 6.6 and 1.8%, respectively, in phase III.
Methanosaeta was about 19.7% in relative abundance in phase I decreasing to 11.4% in phase II, while it almost disappeared in phase III. The halotolerant Methanosaeta harundinacea sp. obtained a 100% similarity on this genus. Concomitantly, Methanobacterium increased from 3.5% in phase I to 11.5% in phase II, but significantly dropped in phase III. In phase III, the hydrogenotrophic methanogens Methanoculleus increased to 6.1% in relative abundance becoming the dominant archaea while Methanolinea decreased significantly. Figure 4 also illustrates that the most prominent amplicon sequence variants belong to 15 and 2 different classes in the bacterial and archaeal domain, respectively.
For the bacteria microbial community, the abundance of microorganisms of the class Actinobacteria, Bacilli, Synergistia, Gammaproteobacteria, and Thermotogae increased correspondingly to 6.8, 5.9, 22.9, 12.6, and 26.6% at the end of the long-term thermophilic operation.

Thermophilic AnMBR Performance and Volatile Fatty Acids Spectrum
The COD removal efficiency started to decrease at the end of the applied phenol loading rate of 0.09 gPh . L −1. d −1 (Figure 2B), i.e., before day 96 in phase I, and subsequently acetate concentration increased (Figure 3). The observed deterioration possibly can be ascribed to an accumulation of a phenol conversion intermediate, such as benzoate (Wang and Barlaz, 1998;Tay et al., 2001), impacting acetoclastic methanogens or syntrophic acetate oxidizers. After the phenol loading rate was further increased to 0.12 gPh . L −1. d −1 , equivalent to a sludge phenol loading rate of 22.6 mgPh . gVSS −1. d −1 , the OLR was decreased to 2.3 gCOD . L −1. d −1 to avoid high volatile fatty acids concentrations in the reactor. The increase in phenol loading rate at the end of phase I, severely impacted COD and phenol conversion, resulting in removal efficiencies below 10%. The phenol concentration in the reactor broth reached 738 mgPh . L −1 (Figure 2A), which apparently inhibited both phenol conversion and methanogenesis. Similarly, Madigou et al. (2016) reported inhibition of phenol conversion at a concentration of 600 mgPh . L −1 , and biogas production nearly stopped when phenol reached 895 mgPh . L −1 . In order to reduce the phenol concentration in the AnMBR, the phenol loading rate was decreased to 0.02 gPh . L −1. d −1 (3.8 mgPh . gVSS −1. d −1 ) on day 133 to prevent further intoxication, while the OLR was reduced to 2.0 gCOD . L −1. d −1 . In phase II, the reactor performance again deteriorated. Consequently, biogas production almost ceased on day 260 when the reactor phenol concentration reached 616 mgPh . L −1 , inhibiting the methanogenic consortium. Surprisingly in phase III, by applying a phenol loading rate of 0.04 gPh . L −1. d −1 (9.3 mgPh . gVSS −1. d −1 ) on day 352, phenol and COD removal efficiencies decreased to 50 and 80%, respectively. The latter indicates an increased sensitivity of the biomass to phenol compared to phase I.
Volatile fatty acid concentrations in the permeate, mainly acetate, indicated a limiting methanogenic conversion capacity for the OLR applied. The high observed acetate concentration reaching 9.6 g . L −1 (Figure 3) when phenol accumulation occurred, could have caused a secondary inhibitory effect to the anaerobic microorganisms such as propionate-oxidizing bacteria (Van Lier et al., 1993). Despite the peaks of 620 and 300 mg . L −1 of butyrate during the phenol accumulation events, all VFA concentrations (propionate, butyrate, and valerate) remained at a relatively low level, commonly observed in anaerobic reactors under anaerobic thermophilic conditions (Van Lier et al., 1993). Caproate always remained below detection level throughout the entire thermophilic AnMBR operation. Previous studies inferred that at 55 • C, phenol-degrading biomass might degrade phenol via n-caproate (Evans and Fuchs, 1988;Fang et al., 2006). On the other hand, Hoyos-Hernandez et al. (2014) have demonstrated that thermophilic phenol degradation is possible via the benzoate conversion route, similar to mesophilic conditions, contrasting these previous studies. Even though benzoate was not determined analytically in the AnMBR, the absence of caproate in any of the VFA analyses strongly suggests that the prevailing phenol conversion pathway was likely via benzoate carboxylation at 55 • C. Following the proposed pathway, benzoate is subsequently de-aromatized to form cyclohexane carboxylic acid, which is cleaved to heptanoate, degraded further through β-oxidation to form valerate, propionate, and acetate, or directly to propionate and butyrate, which are further degraded to acetate (Liang and Fang, 2010).

Conversion Rates and Methanogenic Activity
The assessed SMA values found are similar to those reported elsewhere for phenol degrading biomass under thermophilic conditions, i.e., 0.1 gCOD-CH 4 . gVSS −1. d −1 (Fang et al., 2006). The phenol conversion rate assessed in ex situ batch tests increased to 10.5 ± 3.3 mgPh . gVSS −1. d −1 , inferring a 62% recovery of the conversion capacity after about 17 days of incubation. Interestingly, also Chen et al. (2008) found a maximum phenol conversion rate after 10 to 14 days of incubation. The ex situ assessed phenol conversion rate at the end of the batch tests was about 61% of the AnMBR conversion rate before phenol accumulation occurred (see Table 3). Apparently, the phenol-degrading biomass could recover from the phenol shocks after a relatively short recovery period, while being exposed to low phenol concentrations in the reactor (<100 mgPh . L −1 ). The observed maximum phenol conversion rate of 32.8 ± 0.5 mgPh . gVSS −1. d −1 in batch test at the end of the experiment is higher than the maximum observed phenol conversion rate in the AnMBR during phase I, which likely can be attributed to longterm adaptation of the biomass to phenol after three phases. Comparing our present results with the different studies summarized in Table 5, which are performed in a broad range of temperatures, the observed maximum phenol conversion rates are comparable to those obtained under both mesophilic (Rosenkranz et al., 2013;Franchi et al., 2020) and psychrophilic (Collins et al., 2005;Scully et al., 2006) conditions, using other anaerobic high-rate reactors configurations. Nonetheless, it should be recalled that our present results were obtained under extreme salinity conditions applying sodium concentrations of 18 gNa +. L −1 .
The acetate-fed SMA obtained at the end of the longterm operation was 0.15 ± 0.04 gCOD-CH 4 . gVSS −1. d −1 similar to what was observed at the beginning of phase I and at the end of phase III. Remarkably, the hydrogenotrophic methanogenic activity was a factor 2.3 higher (0.34 ± 0.08 gCOD-CH 4 . gVSS −1. d −1 ), which made us hypothesize that methanogenesis of acetate proceeded via syntrophic acetate oxidation coupled to hydrogenotrophic methanogenesis, rather than aceticlastic methanogenesis. Of interest is the relatively high hydrogenotrophic SMA, which indicates that acetate is possibly syntrophically methanised via oxidation to hydrogen and carbon dioxide. Note that syntrophic acetate oxidation is energetically more favorable at elevated temperature and high acetate concentration and is more often observed as the dominant pathway in a large number of thermophilic anaerobic reactors (van Lier, 1996;Westerholm et al., 2016;Li et al., 2020). Moreover, hydrogenotrophic methanogenesis is commonly observed at elevated salt concentrations (De Vrieze et al., 2016). Based on our present results, follow-up research in a thermophilic AnMBR should reveal the minimum acetate concentration that is required to enhance phenol conversion and to avoid high VFA concentrations in the permeate. In such research, concentrations in the range of 0.3-1.0 gCOD . L −1 acetate, and 500-800 mgPh . L −1 phenol are recommended. Applying similar salinity conditions, we propose an OLR of 2.0-2.2 gCOD . L −1. d −1 and phenol loading rates of 15-20 mgPh . gVSS −1. d −1 in order to maximize the phenol conversion capacity, without compromising the methanogenic activity.

Microbial Community Structure
Petrotoga was enriched to 21.1% during the AnMBR long-term operation. This bacteria has been enriched under anaerobic thermophilic conditions from an oil reservoir containing mostly halophilic species (Dellagnezze et al., 2016). Similarly, high salinity conditions will enrich for salt-tolerant and halophilic Thermovirga and Clostridium species (Muñoz Sierra et al., 2018a). Acetomicrobium hydrogeniformans sp., which also increased during phase III, is found in oil production water and is capable of producing hydrogen. Some species required NaCl for growth (Cook et al., 2018). In the case of the Clostridia class, the genera found Syntrophaceticus are known for their capability of syntrophic acetate oxidation, and Pelotomaculum play an important role in the conversion of phenol and benzoate under methanogenic conditions (Chen et al., 2008).
The abundance of bacteria of the classes Actinobacteria, Bacilli, Synergistia, Gammaproteobacteria, and Thermotogae increased in phase III. Microorganisms corresponding to these classes have been found in thermophilic saline environments. Clostridia class remained in a comparable relative abundance along with the phases, but with changes at genus level. Especially, microorganisms that belong to this class have been reported as the most essential fermentative bacteria and syntrophic phenol degraders as Pelotomaculum (Chen et al., 2009;Muñoz Sierra et al., 2019).
Our results showed a high relative abundance of uncultured microorganisms in the AnMBR. Bacteria belonging to the class JS1, and the candidate phylum Atribacteria (see Supplementary  Figure 1) were dominant during phases I (9.15%) and II (9.41%). Microorganisms belonging to Atribacteria, mostly have been found in deep-sea methane-rich sediments (Carr et al., 2015). Lee et al. (2018) suggested a fermentative role of these microorganisms, capable of using various substrates, and syntrophic acetate oxidation coupled with hydrogen scavenging methanogens. Recently, the first culturable representative strain of this phylum was isolated and it was confirmed that it plays a role in hydrogenogenic fermentative metabolism (Katayama et al., 2019). In our case, the highest species similarity found was Bacillus alkalitolerans strain T3-209, with only 83%.
In the archaea domain, the microbial dynamics indicated a clear switch from acetotrophic to hydrogenotrophic methanogens in the class Methanomicrobia. Both the microbial community structure and the observation that at the end of the experiment the biomass hydrogenotrophic methanogenic activity was substantially higher than the acetate-fed methanogenic activity, support our hypothesis that acetate conversion switched from aceticlastic methanogenesis to acetate oxidation coupled to hydrogenotrophic methanogenesis.

Perspectives and Further Research
The observed AnMBR performance perturbation following phenol accumulation indicated that the cultivated thermophilic phenol-degrading biomass was dependent on the presence of active methanogens. A drop in the hydrogenotrophic methanogenic activity may have led to the observed reduced phenol conversion capacity. It should be noted that the entire experiment was performed under high salinity (18 gNa +. L −1 ) conditions. Further research will focus on the role of syntrophic acetate oxidation and phenol degradation intermediates (e.g., benzoate). Our present study showed that highly saline phenolic wastewaters indeed could be treated in a thermophilic AnMBR. However, the achievable phenol conversion capacity was restricted to 20 mgPh . gVSS −1. d −1 , determining an applicable phenol loading rate of about 0.12 gPh . L −1. d −1 . Although thermophilic operation will bring operational energy benefits when treating high-temperature industrial wastewaters with the target of process water reuse, the phenol conversion capacity of the reactor will be lower than when opting for mesophilic operation.

CONCLUSION
Maximum COD and phenol removal efficiencies of about 95% were achieved during the long-term thermophilic AnMBR operation at 18 gNa +. L −1 . However, severe perturbations occurred following relatively small increments in the phenol loading rate from 0.01 to 0.12 gPh . L −1. d −1 . Moreover, by exceeding a sludge phenol loading rate of 20 mgPh . gVSS −1. d −1 the system immediately responded in phenol build-up to concentrations higher than 600 mgPh . L −1 leading to significant deterioration of methanogenesis. The observed maximum phenol conversion rates were 21.3 ± 0.2 and 32.8 ± 0.5 mgPh . gVSS −1. d −1 in the AnMBR, and in ex situ batch test at the end of the reactor operation, respectively. The absence of caproate in the VFAs spectrum inferred that the phenol conversion pathway was likely via benzoate carboxylation. The assessed hydrogenotrophic SMA was a factor 2.3 higher than the acetate-fed SMA. Correspondingly, microbial population dynamics indicated that hydrogenotrophic methanogens were enriched during the long-term operation and Clostridia class was dominant. Overall, thermophilic AnMBR operation under high salinity seemed to be susceptible to sudden increase in phenol loading rate or phenol shocks, indicating that the specific phenol conversion capacity under the studied conditions was limiting the treatment process.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ebi.ac.uk/ ena, PRJEB38467.