Spatial Structure of <i>Yersinia pestis</i> Population Belonging to Altai Biovar, Subspecies central asiatica Acording to Genome-Wide Sequencing Data

Abstract

The aim of the work was to conduct phylogenetic analysis of Y. pestis strains, biovar altaica ssp. central asiatica, isolated in Gorno-Altai high-mountain and Sailyugem natural plague foci on the territory of Russia and Mongolia in 1965–2020, by full-genome sequencing data.Materials and methods. To determine the population structure of the biovar altaica ssp. central asiatica, 34 whole genome sequences were used (including 20 Y. pestis strains of the biovar altaica, 18 of which were sequenced by us). To isolate DNA from the Y. pestis strains, a PureLink Genomic DNA Mini Kit (Invitrogen, USA) was applied. Sequencing of the nucleotide sequences of Y. pestis strains was carried out in Ion PGM system Lifetechnologies. The analysis and processing of the obtained data were performed with the help of Newblergs Assembler 2.6 and IonTorrent Suite software package, 3.4.2. The search for SNPs was performed using the Wombac 2.0 program. The Maximum Likelihood dendrogram was built in the PhyML 3.1. The dendrogram was visualized using the FigTree 1.4.3 software.Results and discussion. Based on the data of whole genome analysis, taking into account the 1871 revealed polymorphic nucleotides, the spatial structure of the biovar altaica ssp. central asiatica has been determined. It includes several phylogeographic branches: the Kurai-Tarkhatinskaya (cluster 0.PE4a-1) and the Ulandryk-Mongolian (0.PE4a-2), which is in agreement with the geographical regions of the isolation of strains forming these branches in the Altai Mountains. The Kurai-Tarkhatinskaya branch is further divided into the Kurai (sub-cluster 0.PE4a-1-1, formed by the strains of 2009–2018) and Tarkhatinskaya (subcluster 0.PE4a-1-2, formed by the strains of 2012–2020) sub-branches, while the Ulandryk-Mongolian branch of evolution is split into sub-branches represented by strains from the Ulandryk meso focus (sub-cluster 0.PE4a-2-2, strains 1965–2010) and the Sailyugem focus of Mongolia (sub-cluster 0.PE4a-2-1, strains 1964–1990).