Comparative genomic analysis of Thermus provides insights into the evolutionary history of an incomplete denitrification pathway.

Biological denitrification is a crucial process in the nitrogen biogeochemical cycle, and Thermus has been reported to be a significant heterotrophic denitrifier in terrestrial geothermal environments. However, neither the denitrification potential nor the evolutionary history of denitrification genes in the genus Thermus or phylum Deinococcota is well understood. Here, we performed a comparative analysis of 23 Thermus genomes and identified denitrification genes in 15 Thermus strains. We confirmed that Thermus harbors an incomplete denitrification pathway as none of the strains contain the nosZ gene. Ancestral character state reconstructions and phylogenetic analyses showed that narG, nirS, and norB genes were acquired by the last common ancestor of Thermales and were inherited vertically. In contrast, nirK of Thermales was acquired via two distinct horizontal gene transfers from Proteobacteria to the genus Caldithermus and from an unknown donor to the common ancestor of all known Thermus species except Thermus filiformis. This study expands our understanding of the genomic potential for incomplete denitrification in Thermus, revealing a largely vertical evolutionary history of the denitrification pathway in the Thermaceae, and supporting the important role for Thermus as an important heterotrophic denitrifier in geothermal environments.