Divergent dual C-H isotopic fractionation pattern during anaerobic biodegradation of toluene within Aromatoleum species under nitrate-reducing conditions

Toluene is a pollutant frequently detected in contaminated groundwater, mostly due to leakage from underground gasoline storage tanks and pipeline ruptures. Multi-element compound-specific isotope analysis provides a framework to understand transformation processes and design efficient remediation strategies. In this study, we enriched an anaerobic bacterial culture derived from a BTEX-contaminated aquifer that couples toluene and phenol oxidation with nitrate reduction and the concomitant production of carbon dioxide and biomass. The 16S rRNA gene amplicon data indicated that the toluene-degrading consortium was dominated by an Aromatoleum population (87 ± 2 % relative abundance), and metagenome sequencing confirmed that the genome of this Aromatoleum sp. encoded glycyl-radical enzyme benzylsuccinate synthase (BssABC) and phenylphospate synthase (PpsA1BC) homologous genes involved in the first step of toluene and phenol transformation, respectively. Carbon and hydrogen isotopic fractionation were ε