This article describes Bacillus anthracis strains isolated during an outbreak of anthrax on the Yamal Peninsula in the summer of 2016 and independently in Yakutia in 2015. A common feature of these strains is their conservation in permafrost, from which they were extracted either due to the thawing of permafrost (Yamal strains) or as the result of paleontological excavations (Yakut strains). All strains isolated on the Yamal share an identical genotype belonging to lineage B.Br.001/002, pointing to a common source of infection in a territory over 250 km in length. In contrast, during the excavations in Yakutia, three genetically different strains were recovered from a single pit. One strain belongs to B.Br.001/002, and whole genome sequence analysis showed that it is most closely related to the Yamal strains in spite of the remoteness of Yamal from Yakutia. The two other strains contribute to two different branches of A.Br.008/011, one of the remarkable polytomies described so far in the B . anthracis species. The geographic distribution of the strains belonging to A.Br.008/011 is suggesting that the polytomy emerged in the thirteenth century, in combination with the constitution of a unified Mongol empire extending from China to Eastern Europe. We propose an evolutionary model for B . anthracis recent evolution in which the B lineage spread throughout Eurasia and was subsequently replaced by the A lineage except in some geographically isolated areas.

20This article describes Bacillus anthracis strains isolated during an outbreak of anthrax on 21 the Yamal Peninsula in the summer of 2016 and independently in Yakutia in 2015. A common 22 feature of these strains is their conservation in permafrost, from which they were extracted either 23 due to the thawing of permafrost (Yamal strains) or as the result of paleontological excavations 24 (Yakut strains). All strains isolated on the Yamal share an identical genotype belonging to lineage 2 25 B.Br.001/002, pointing to a common source of infection in a territory over 250 km in length. In 26 contrast, during the excavations in Yakutia, three genetically different strains were recovered from 27 a single pit. One strain belongs to B.Br.001/002, as the Yamal strains. Despite the remoteness of 28 Yamal from Yakutia, whole genome sequence analysis showed that the B.Br.001/002 strains are 29 very closely related. The two other strains contribute to two different branches of A.Br.008/011, 30 one of the remarkable polytomies described so far in B. anthracis population. The geographic 31 distribution of the strains belonging to this polytomy is suggesting that this polytomy emerged in 32 the thirteenth century, in combination with the constitution of a unified Mongol empire extending 33 from China to Eastern Europe. We propose an evolutionary model for B. anthracis recent evolution 34 in which the B lineage spread throughout Eurasia and was subsequently replaced by the A lineage 35 except in some geographically isolated areas.36

Tularemia is a severe infectious disease caused by the Gram-negative bacteria Fracisella tularensis. There are four subspecies of F.tularensis: holarctica, tularensis, mediasiatica, and novicida, which differ in their virulence and geographic distribution. One of them, subsp. mediasiatica remains extremely poorly studied, primarily due to the fact that it is found only in the sparsely populated regions of Central Asia and Russia. In particular there is little information in the literature on the virulence and pathogenicity of subsp. mediasiatica. In the present article, we evaluated the comparative virulence of subsp. mediasiatica in vaccinated laboratory animals which we infected with virulent strains: subsp. mediasiatica 678, subsp. holarctica 503, and subsp. tularensis SCHU within 60 to 180 days after vaccination. We found that subsp. mediasiatica is comparable in pathogenicity in mice with subsp. tularensis and in guinea pigs with subsp. holarctica. We also found that the live vaccine does not fully protect mice from subsp. mediasiatica but completely protects guinea pigs for at least six months. In general, our data suggest that subsp. mediasiatica occupies an intermediate position in virulence between spp. tularensis and holarctica.

The main pathogenic factor of Bacillus anthracis is a three-component toxin encoded by the pagA, lef, and cya genes, which are located on the pXO1 plasmid. The atxA gene, which encodes the primary regulator of pathogenicity factor expression, is located on the same plasmid. In this work, we evaluated the polymorphism of the pagA, lef, cya, and atxA genes for 85 B. anthracis strains from different evolutionary lineages and canSNP groups. We have found a strong correlation of 19 genotypes with the main evolutionary lineages, but the correlation with the canSNP group of the strain was not as strong. We have detected several genetic markers indicating the geographical origin of the strains, for example, their source from the steppe zone of the former USSR. We also found that strains of the B.Br.001/002 group caused an anthrax epidemic in Russia in 2016 and strains isolated during paleontological excavations in the Russian Arctic have the same genotype as the strains of the B.Br.CNEVA group circulating in Central Europe. This data could testify in favor of the genetic relationship of these two groups of strains and hypothesize the ways of distribution of their ancestral forms between Europe and the Arctic.

Climate change brings new risks of emergence of especially dangerous diseases. The paper reports the possibility of assessing the pathogenic potential of bacteria as demonstrated by studying the allelic polymorphism of anthrax bacterium pathogenicity factor genes, which is a prerequisite for assessing the associated microbiological risks. The allelic polymorphism of the capBCADE operon (capB, capC, capA, capD, and capE genes) encoding the capsule biosynthesis proteins of Bacillus anthracis, and the acpA and acpB genes encoding the expression regulators of this operon have been studied for the first time. A number of single nucleotide polymorphisms (SNPs) were described in the strains of the studied sample, including 5 SNPs in the capB gene, 3 in capC, 4 in capA, 14 in capD, 2 in capE, and 15 in acpB, as well as 7 SNPs and one insertion in the acpA gene. As a result, the sample has been divided into sequence types for each gene and 17 genotypes, which are combinations of the identified sequence types. In silico translation of the detected alleles of the studied genes revealed three isoforms of the CapB and CapA proteins, two isoforms of the CapC and CapE proteins, six isoforms of the CapD protein, five isoforms of the AcpA protein, and four isoforms in the AcpB protein. It has been demonstrated that the SNP in the 351A → G position of capC is a marker of A.Br.Aust94 group strains. Based on the results, A.Br.Vollum group strains were divided into two subgroups. The strains in the evolutionary lines B and C differed from the line A strains by the presence of an 853G → A SNP in the acpA gene. In addition, a previously unknown variable number tandem repeat (VNTR), has been found in the acpA gene and the possibility of using it for differentiating and genotyping of B. anthracis strains has been demonstrated.