Genotyping and population diversity of Bacillus anthracis in China based on MLVA and canSNP analysis

In pastoral parts of China, anthrax still presents a major risk to livestock and threatens the health of local human populations. Currently, whole-genome-based molecular markers, such as single-nucleotide polymorphisms (SNPs) and variable number tandem repeats (VNTRs), are the most effective tools for genotyping Bacillus anthracis. In this study, 191 isolates were selected to assess the diversity of B. anthracis in China. Five isolates were confirmed not to be B. anthracis by clustered regularly interspaced short palindromic repeat analysis, while the remaining 186 isolates were typed using canonical SNP (canSNP) and VNTR analyses. Five sublineages/subgroups, A.Br.001/002, A.Br.Vollum, A.Br.Aust.94, A.Br.Ames, and A.Br.008/009, were detected based on 13 canSNP sites. The 186 isolates were further assigned 114 sequence types based on 27 VNTR loci, with major branches correlating with the canSNP analysis. We then used a simplified multiple-locus variable number tandem repeat analysis (MLVA) protocol (MLVAmin) based on eight high-resolution VNTR sites to analyze the Chinese isolates, with the resulting phylogeny again agreeing with the canSNP analysis. We also developed two schemes, MLVAc and MLVAp, using various numbers of VNTRs to analyze different canSNP sublineages to increase the typing resolution of the canSNP protocol. The results showed a highly imbalanced geographical distribution of the B. anthracis population, with four different sublineages observed in Xinjiang Province, while only one sublineage, A.Br.001/002, was found in the other six provinces, except for three A.Br.Ames strains isolated from Inner Mongolia. Based on the MLVA and canSNP analysis, the spread of B. anthracis appears to have occurred from west to east via three independent routes.