Phylogenetics of Bacillus anthracis isolates from Russia and bordering countries

Abstract

Anthrax is a concern for public health and veterinary medicine in Russia. The available phylogenetic data on isolates from Russia and neighboring CIS countries are clearly not enough to gain a better understanding of their position in the global phylogenetic population structure of this pathogen. In this study, we analyzed the genomes of 66 Bacillus anthracis strains, which were isolated between 1935 and 2019 from different sources in Russia, as well as in Ukraine, Azerbaijan, Georgia, Armenia and Moldova. Whole genome SNP analysis of genomes of 66 strains obtained in this study along with 222 B. anthracis genomes available in the GenBank database revealed 7242 SNPs used to construct a phylogenetic reconstruction with the method of Maximum Likelihood. Studied strains belong to 6 different genetic groups: A.Br.008(A.Br.008/009), A.Br.081(Ames), A.Br.014(A.Br.Aust94), A.Br.082(A.Br.001/002), A.Br.034(A.Br.005/006, Ancient A) and B.Br.002 (B.Br.001/002). Within the group A.Br.014(A.Br.Aust94) a subcluster A.Br.029 of strains isolated in Georgia, Armenia, Azerbaijan, Russia (Republic of Dagestan) and Turkey, named Caucasus-East Anatolia (CEA), was identified. In the subgroup A.Br.105(Tsiankovskii) the cluster A.Br.117 of strains from Russia, Ukraine and Slovakia are assigned, in the subgroup A.Br118 (STI) - cluster A.Br.123 with strains from Russia and Georgia and cluster A.Br.125 with strains from Republic of Dagestan. New subclusters B.Br.017(“EUROPE”) were identified in the B.Br.002(B.Br.001/002) cluster, represented by strains from the European part of Russia, as well as from South Korea and Finland. For 8 clusters and subclusters, the SNP markers were identified. The study confirmed a significant genetic diversity of the strains isolated in Russia and border countries and clarified their position in the phylogenetic structure of the global B. anthracis population. New genetic clusters A.Br.029 (CEA), A.Br.117, A.Br.123, A.Br.125, and B.Br.017 («EUROPE») were defined. 96 marker SNPs specific for these clusters were identified