Correlations between microbial population dynamics, bamA gene abundance and performance of anaerobic sequencing batch reactor (ASBR) treating increasing concentrations of phenol

•Microbial community composition in the ASBR changed upon phenol concentration increments during the process.•BamA gene abundance was used as a proxy to estimate the dynamics of phenol degrading bacteria in the ASBR system.•The abundance of Syntrophorhabdus, Chloroflexus, Smithella, Methanolinea and Methanosaeta correlated with degradation rate. The relevant microorganims driving efficiency changes in anaerobic digestion of phenol remains uncertain. In this study correlations were established between microbial population and the process performance in an anaerobic sequencing batch reactor (ASBR) treating increasing concentrations of phenol (from 120 to 1200 mg L−1). Sludge samples were taken at different operational stages and microbial community dynamics was analyzed by 16S rRNA sequencing. In addition, bamA gene was quantified in order to evaluate the dynamics of anaerobic aromatic degraders. The microbial community was dominated by Anaerolineae, Bacteroidia, Clostridia, and Methanobacteria classes. Correlation analysis between bamA gene copy number and phenol concentration were highly significant, suggesting that the increase of aromatic degraders targeted by bamA assay was due to an increase in the amount of phenol degraded over time. The incremental phenol concentration affected hydrogenotrophic archaea triggering a linear decrease of Methanobacterium and the growth of Methanobrevibacter. The best performance in the reactor was at 800 mg L−1 of phenol. At this stage, the highest relative abundances of Syntrophorhabdus, Chloroflexus, Smithella, Methanolinea and Methanosaeta were observed and correlated positively with initial degradation rate, suggesting that these microorganisms are relevant players to maintain a good performance in the ASBR.