Comparative genomic analysis reveals distinct population structure in Legionella anisa

Legionella anisa has been frequently isolated from engineered water systems; however, its population structure remains understudied compared to Legionella pneumophila. Here, we generated complete genome sequences for four L. anisa isolates recovered from a healthcare facility in Rimouski, Canada. Further the population structure of this species was investigated by performing comparative genomic analyses of the genomes generated in this study together with publicly available L. anisa genomes. Genome-wide phylogenetic analysis revealed the presence of three distinct clades separated by substantial genetic divergence (∼500 SNP), with the Rimouski isolates forming a tightly clustered group, suggesting a clonal lineage. Comparative pangenome analysis indicated moderate core genome conservation accompanied by a highly variable accessory genome (∼50%). The isolates characterized in this study harbored multiple plasmids encoding genes associated with conjugation, heavy metal resistance, and other stress-related functions, suggesting potential roles in environmental persistence. Previous studies have shown that L. anisa can proliferate within protozoan host cells, although outcomes vary depending on the host species. Our isolates showed efficient proliferation within Acanthamoeba castellanii, but not within Vermamoeba vermiformis, under the conditions tested. Together, these findings underscore the genomic diversity of this understudied Legionella species and provide a framework for future investigations regarding environmental persistence and potential pathogenicity.