Genus Thermus is the main focus of researcher among the thermophiles. Members of this genus are the inhabitants of both natural and artificial thermal environments. We performed phylogenomic analyses and comparative genomic studies to unravel the genomic diversity among the strains belonging to the genus Thermus in geographically different thermal springs. Sixteen Thermus strains were isolated and sequenced from hot springs, Qucai hot springs in Tibet and Tengchong hot springs in Yunnan, China. 16S rRNA gene based phylogeny and phylogenomic analyses based on concatenated set of 971 Orthologous Protein Families (supermatrix and gene content methods) revealed a mixed distribution of the Thermus strains. Whole genome based phylogenetic analysis showed, all 16 Thermus strains belong to five species; Thermus oshimai (YIM QC-2-109, YIM 1640, YIM 1627, 77359, 77923, 77838), Thermus antranikianii (YIM 73052, 77412, 77311, 71206), Thermus brokianus (YIM 73518, 71318, 72351), Thermus hydrothermalis (YIM 730264 and 77927) and one potential novel species 77420 forming clade with Thermus thalpophilus SYSU G00506T. Although the genomes of different strains of Thermus of same species were highly similar in their metabolic pathways, but subtle differences were found. CRISPR loci were detected through genome-wide screening, which showed that Thermus isolates from two different thermal locations had well developed defense system against viruses and adopt similar strategy for survival. Additionally, comparative genome analysis screened competence loci across all the Thermus genomes which could be helpful to acquire DNA from environment. In the present study it was found that Thermus isolates use two mechanism of incomplete denitrification pathway, some Thermus strains produces nitric oxide while others nitrious oxide (dinitrogen oxide), which show the heterotrophic lifestyle of Thermus genus. All isolated organisms encoded complete pathways for glycolysis, tricarboxylic acid and pentose phosphate. Calvin Benson Bassham cycle genes were identified in genomes of T. oshimai and T. antranikianii strains, while genomes of all T. brokianus strains and organism 77420 were lacking. Arsenic, cadmium and cobalt-zinc-cadmium resistant genes were detected in genomes of all sequenced Thermus strains. Strains 77,420, 77,311, 73,518, 77,412 and 72,351 genomes were found harboring genes for siderophores production. Sox gene clusters were identified in all sequenced genomes, except strain YIM 730264, suggesting a mode of chemolithotrophy. Through the comparative genomic analysis, we also identified 77420 as the genome type species and its validity as novel organism was confirmed by whole genome sequences comparison. Although isolate 77420 had 99.0% 16S rRNA gene sequence similarity with T. thalpophilus SYSU G00506T but based on ANI 95.86% (Jspecies) and digital DDH 68.80% (GGDC) values differentiate it as a potential novel species. Similarly, in the phylogenomic tree, the novel isolate 77,420 forming a separate branch with their closest reference type strain T. thalpophilus SYSU G00506T.
