Reduction, evolutionary pattern and positive selection of genes encoding formate dehydrogenase in Wood-Ljungdahl pathway of gastrointestinal acetogens suggests their adaptation to formate-rich habitats

Acetogens are anaerobes using Wood-Ljungdahl pathway (WLP) as the terminal electron acceptor for both assimilation and dissimilation of CO2 and widely distributed in diverse habitats. However, their habitat adaptation is often unclear. Given that bacterial genome evolution is often the result of environmental selective pressure, hereby we analysed gene copy number, phylogeny and selective pressure of genes involved in WLP within known genomes of 43 species to study the habitat adaption of gastrointestinal acetogens. The gene copy number of formate dehydrogenase (FDH) in gastrointestinal acetogens was much lower than that of non-gastrointestinal acetogens, and in five cases, no FDH genes were found in the genomes of five gastrointestinal acetogens, but that of the other WLP genes showed no difference. The evolutionary pattern of FDH genes was significantly different from that of the other enzymes. Additionally, seven positively selected sites were only identified in the fdhF genes, which means fdhF mutations favoured their adaptation. Collectively, reduction or loss of FDH genes and their evolutionary pattern as well as positive selection in gastrointestinal acetogens indicated their adaptation to formate-rich habitats, implying that FDH genes catalysing CO2 reduction to formate as the first step of methyl branch of WLP may have evolved independently.