Canonical Single Nucleotide Polymorphisms Research Articles

Deep-amplicon sequencing of the complete beta-tubulin gene in Trichuris trichiura before and after albendazole treatment

Concerns about the emergence of benzimidazole resistance in soil-transmitted helminths (STH) infections, particularly against Trichuris trichiura, have arisen. Previous studies of veterinary nematodes have linked benzimidazole resistance to single-nucleotide polymorphisms (SNPs) at three specific codons in the beta-tubulin gene, but similar associations in STH have not been consistently observed. In this work, we screened the complete beta-tubulin gene previously linked to benzimidazole resistance in T. trichiura by deep-amplicon sequencing to identify genetic variants and associate levels of diversity with drug response to albendazole. We used 99 DNA samples extracted from T. trichiura pooled eggs, previously semi-purified from human stool samples collected in Manhiça district, Mozambique. We obtained a set of 39 amplicons of the complete gene by subjecting the pooled eggs to long-read PCR and subsequently sequencing them. Of those amplicons, 22 and 17 were obtained from stool samples collected before, and 21 days after albendazole treatment, respectively. We observed genetic variation across the whole gene sequence, in both exons and introns; however, none were associated with the previously proposed resistance-associated SNPs, and none were predicted to significantly affect protein function. No significant differences in genetic diversity were observed between pre- and post-treatment samples. Using publicly available genome-wide data, we also analysed a second beta-tubulin isotype in the T. trichiura genome. We focused on detecting the canonical SNPs and assessing for signatures of genetic selection around this second isotype gene. This analysis did not reveal evidence supporting this second isotype's role in anthelmintic resistance. Despite the limitations of our study, such as a small sample size, particularly paired pre- and post-treatment samples (n=6), or a restricted geographical area, we found no evidence linking either of the two beta-tubulin genes to benzimidazole resistance in T. trichiura, suggesting that genetic markers of drug resistance likely exist outside the beta-tubulin genes.

New Insights into the Phylogeny of the A.Br.161 ("A.Br.Heroin") Clade of Bacillus anthracis.

Bacillus anthracis is a rare but highly dangerous zoonotic bacterial pathogen. At the beginning of this century, a new manifestation of the disease, injectional anthrax, emerged as a result of recreational heroin consumption involving contaminated drugs. The organisms associated with this 13-year-lasting outbreak event in European drug consumers were all grouped into the canonical single-nucleotide polymorphism (canSNP) clade A-branch (A.Br.) 161 of B. anthracis. Related clade A.Br.161 strains of B. anthracis not associated with heroin consumption have also been identified from different countries, mostly in Asia. Because of inadvertent spread by anthropogenic activities, other strains of this A.Br.161 lineage were, however, isolated from several countries. Thus, without additional isolates from this clade, its origin of evolution or its autochthonous region remains obscure. Here, we genomically characterized six new A.Br.161 group isolates, some of which were from Iran, with others likely historically introduced into Germany. All the chromosomes of these isolates could be grouped into a distinct sub-clade within the A.Br.161 clade. This sub-clade is separated from the main A.Br.161 lineage by a single SNP. We have developed this SNP into a PCR assay facilitating the future attribution of strains to this group.

Whole Genome-Sequencing and Phylogenetic Analysis of Bacillus anthracis from 2019-2023 in Ningxia Hui Autonomous Region, China.

Since 2019, the incidence of anthrax in the Ningxia Hui Autonomous Region has increased significantly compared with previous years, so in this situation the anthrax in the Ningxia region not only had a detrimental impact on public health, but also inflicted significant economic repercussions. Therefore, we conducted a molecular epidemiological study of 20 strains from 2019-2023 isolates. This study investigated the origin of Bacillus anthracis and its genetic diversity. We conducted canonical single-nucleotide polymorphisms (CanSNPs) typing and whole genome sequencing based on the extracted nucleic acid of Bacillus anthracis. Based on the whole genome drafts, we studied the genomic characteristics of 20 isolates. Meanwhile, we performed phylogenetic studies based on genome-wide core single-nucleotide polymorphisms (SNPs) using MEGA's Maximum Likelihood (ML) method and core-genome-based multilocus sequence typing (cgMLST) of the core genomes of these strains using BioNumerics' minimum spanning tree (MST) model. The 20 isolates were categorized into sub-lineages A.Br.001/002, and comparative genomic analyses of these strains with other isolates from other parts of the world showed that the strains from Ningxia were correlated with isolates from Europe, Indonesia, Georgia (USA), and Beijing (China). For the 20 isolates in Ningxia, the genetic relationship of the isolates isolated from the same year or region was relatively close. The A.Br.001/002 subgroup was the dominant endemic strain in Ningxia. The genetic relationship and phylogenesis between isolates from Ningxia and strains from Europe and Indonesia suggest that anthrax spread around the globe through ancient trade routes.

Spatial and phylogenetic patterns reveal hidden infection sources of Bacillus anthracis in an anthrax outbreak in Son La province, Vietnam.

Bacillus anthracis, the bacterial cause of anthrax, is a zoonosis affecting livestock and wildlife often spilling over into humans. In Vietnam, anthrax has been nationally reportable since 2015 with cases occurring annually, mostly in the northern provinces. In April 2022, an outbreak was reported in Son La province following the butchering of a water buffalo, Bubalus bubalis. A total of 137 humans from three villages were likely exposed to contaminated meat from the animal. Early epidemiological investigations suggested a single animal was involved in all exposures. Five B. anthracis isolates were recovered from human clinical cases along with one from the buffalo hide, another from associated maggots, and one from soil at the carcass site. The isolates were whole genome sequenced, allowing global, regional, and local molecular epidemiological analyses of the outbreak strains. All recovered B. anthracis belong to the A.Br.001/002 lineage based on canonical single nucleotide polymorphism analysis (canSNP). Although not previously identified in Vietnam, this lineage has been identified in the nearby countries of China, India, Indonesia, Thailand, as well as Australia. A twenty-five marker multi-locus variable number tandem repeat analysis (MLVA-25) was used to investigate the relationship between human, soil, and buffalo strains. Locally, four MLVA-25 genotypes were identified from the eight isolates. This level of genetic diversity is unusual for the limited geography and timing of cases and differs from past literature using MLVA-25. The coupled spatial and phylogenetic data suggest this outbreak originated from multiple, likely undetected, animal sources. These findings were further supported by local news reports that identified at least two additional buffalo deaths beyond the initial animal sampled in response to the human cases. Future outbreak response should include intensive surveillance for additional animal cases and additional molecular epidemiological traceback to identify pathogen sources.

Phenotypic and genotypic discrimination of Francisella tularensis ssp. holarctica clades

Francisella tularensis is the causative agent of tularemia, a zoonotic disease with a wide host range. F. tularensis ssp. holarctica (Fth) is of clinical relevance for European countries, including Germany. Whole genome sequencing methods, including canonical Single Nucleotide Polymorphism (canSNP) typing and whole genome SNP typing, have revealed that European Fth strains belong to a few monophyletic populations. The majority of German Fth isolates belong to two basal phylogenetic clades B.6 (biovar I) and B.12 (biovar II). Strains of B.6 and B.12 seem to differ in their pathogenicity, and it has been shown that strains of biovar II are resistant against erythromycin. In this study, we present data corroborating our previous data demonstrating that basal clade B.12 can be divided into clades B.71 and B.72. By applying phylogenetic whole genome analysis as well as proteome analysis, we could verify that strains of these two clades are distinct from one another. This was confirmed by measuring the intensity of backscatter light on bacteria grown in liquid media. Strains belonging to clades B.6, B.71 or B.72 showed clade-specific backscatter growth curves. Furthermore, we present the whole genome sequence of strain A-1341, as a reference genome of clade B.71, and whole proteomes comparison of Fth strains belonging to clades B.6, B.71 and B.72. Further research is necessary to investigate phenotypes and putative differences in pathogenicity of the investigated different clades of Fth to better understand the relationship between observed phenotypes, pathogenicity and distribution of Fth strains.

Phylogeography of Francisella tularensis subspecies holarctica and epidemiology of tularemia in Switzerland.

Tularemia, an endemic disease that mainly affects wild animals and humans, is caused by Francisella tularensis subsp. holarctica (Fth) in Switzerland. The Swiss Fth population consist of multiple different subclades which are distributed throughout the country. The aim of this study is to characterize the genetic diversity of Fth in Switzerland and to describe the phylogeographic relationship of isolates by single nucleotide polymorphism (SNP) analysis. This analysis is combined with human surveillance data from reported cases over the last 10 years and in vitro and in silico antibiotic resistance tests to provide insight into the epidemiology of tularemia in Switzerland. We sequenced the whole genomes of 52 Fth strains of human or tick origin collected in Switzerland between 2009 and 2022 and analyzed together with all publicly available sequencing data of Swiss and European Fth. Next, we performed a preliminary classification with the established canonical single nucleotide polymorphism nomenclature. Furthermore, we tested 20 isolates from all main Swiss clades for antimicrobial susceptibility against a panel of antimicrobial agents. All 52 sequenced isolates from Switzerland belong to major clade B.6, specifically subclades B.45 and B.46, previously described in Western Europe. We were able to accurately reconstruct the population structure according to the global phylogenetic framework. No resistance to clinically recommended antibiotics could be identified in vitro or in silico in the western B.6 strains.

Identification of the molecular characteristics of Bacillus anthracis (1982-2020) isolates in East Indonesia using multilocus variable-number tandem repeat analysis.

Background and Aim:Anthrax is one of the endemic strategic diseases in East Indonesia, particularly in the provinces of South Sulawesi, West Sulawesi, Gorontalo, East Nusa Tenggara, and West Nusa Tenggara. Anthrax is an important disease due to its zoonotic and economic impact on the livestock industry. This study aimed to identify the molecular characteristics of Bacillus anthracis in East Indonesia using multilocus variable-number tandem repeat (VNTR) analysis (MLVA).Materials and Methods:Isolates were obtained from an investigation of anthrax outbreaks in five provinces of East Indonesia from 1982 to 2020. Conventional polymerase chain reaction for B. anthracis was used to identify MLVA-8. Deoxyribonucleic acid sequencing analysis was based on MLVA-8 primers for VNTR identification of the phylogenetic relationship among 24 isolates of B. anthracis obtained from 17 distinct districts/cities in East Indonesia. Tandem Repeats Finder was used for VNTR identification, and Molecular Evolutionary Genetics Analysis X was used to construct phylogenetic analysis.Results:In this study, 24 isolates were classified as genotype or lineage A. There were four subgroups of B. anthracis circulating in East Indonesia based on eight molecular marker loci sequence results.Conclusion:The findings of this study show that MLVA-8 typing might be useful as a subtyping tool for the epidemiological investigation of identical genotypes and low genetic diversity of B. anthracis. No other lineage of B. anthracis was circulating in East Indonesia. Other molecular methods are needed, such as extended MLVA, whole-genome sequencing, and canonical single-nucleotide polymorphism, for a more precise study of B. anthracis genetic diversity.

Investigation on Anthrax in Bangladesh during the Outbreaks of 2011 and Definition of the Epidemiological Correlations

In 2011, in Bangladesh, 11 anthrax outbreaks occurred in six districts of the country. Different types of samples were collected from May to September in the six districts where anthrax had occurred in order to detect and type Bacillus anthracis (B. anthracis) strains. Anthrax was detected in 46.6% of the samples analysed, in particular in soils, but also in bone samples, water, animal feed, and rumen ingesta of dead animals. Canonical single nucleotide polymorphisms (CanSNPs) analysis showed that all the isolates belonged to the major lineage A, sublineage A.Br.001/002 of China and Southeast Asia while the multi-locus variable number of tandem repeats (VNTRs) analysis (MLVA) with 15 VNTRs demonstrated the presence of five genotypes, of which two resulted to be new genotypes. The single nucleotide repeats (SNRs) analysis showed 13 SNR types; nevertheless, due to its higher discriminatory power, the presence of two isolates with different SNR-type polymorphisms was detected within two MLVA genotypes. This study assumes that soil is not the only reason for the spread of the disease in Bangladesh; contaminated feed and water can also play an important role in the epidemiology of anthrax. Possible explanations for these epidemiological relationships are discussed.

High-resolution melting PCR assay as a powerful tool for the epidemiological surveillance of tularemia in Western Europe.

In Europe, tularemia is caused by Francisella tularensis subsp. holarctica and is a sporadic disease affecting mainly wildlife animals and humans. Classification of this species relies on canonical single nucleotide polymorphisms (canSNPs). Four main clades have been described for F. tularensis subsp. holarctica: B.4, B.6, B.12 and B.16. Phylogeographic studies have shown that clade B.6 is predominant in Western Europe and B.12 in Eastern and Central Europe. Based on this global phylogeny, we aimed to design a molecular typing assay for all genetic subclades of subclade B.11, which is the predominant subclade in clade B.6. We designed high-resolution melting (HRM) primers for the screening of 109 canSNPs divided in seven orders of discrimination for the molecular epidemiology analysis and tracking of Francisella tularensis subsp. holarctica in Western Europe.

Retrospective Analysis of the Relationship between Two Anthrax Outbreaks in Kazakhstan Based on Genomic Data.

We present a retrospective analysis of strains from two anthrax outbreaks in western Kazakhstan in 2009. The outbreaks occurred during the same period and in the same area located close to main roads, favoring a single source of infection. However, multilocus variable-number tandem-repeat analysis (MLVA), canonical single-nucleotide polymorphism (CanSNP) analysis, and genome-wide analysis demonstrated that the outbreaks were not connected.

Genotyping of Francisella tularensis subsp. holarctica from Hares in Germany.

Francisella tularensis is the causative agent of the zoonotic disease tularemia. In Germany, most human infections are caused by contact with infected hares. The aim of this study was to characterize Francisella tularensis subsp. holarctica strains isolated from hares in Germany and to develop bioinformatics tools to analyze their genetic relatedness. In total, 257 German isolates—obtained mainly from hares (n = 233), other vertebrate animals, and ticks, but also from humans (n = 3)—were analyzed within this study. Publically available sequence data from 49 isolates were used to put our isolates into an epidemiological context and to compare isolates from natural foci and humans. Whole-genome sequences were analyzed using core-genome Multi-Locus-Sequence-Typing, canonical Single Nucleotide Polymorphism (SNP) typing and whole-genome SNP typing. An overall conformity of genotype clustering between the typing methods was found, albeit with a lower resolution for canonical single SNP typing. The subclade distribution, both on local and national levels, among strains from humans and hares was similar, suggesting circulation of the same genotypes both in animals and humans. Whilst close to identical isolates of the same subclade were found distributed over large areas, small geographical foci often harbored members of different subclades. In conclusion, although genomic high-resolution typing was shown to be robust, reproducible and allowed the identification of highly closely related strains, genetic profiling alone is not always conclusive for epidemiological linkage of F. tularensis strains.

Changed epidemiology of anthrax and molecular characteristics of Bacillus anthracis in Inner Mongolia Autonomous Region, China.

Anthrax is a natural foci disease in Inner Mongolia, which poses a severe threat to public health. In this study, the incidence number, rate and constituent ratio were used to describe the epidemiological characteristics of anthrax in the region from 1956-2018. The molecular correlation and genetic characteristics of the strains were investigated using canonical single nucleotide polymorphisms (CanSNP), multiple-locus variable-number tandem repeat analysis (MLVA-15) and whole genome sequencing (WGS). The epidemiological characteristics of anthrax in Inner Mongolia have altered significantly. The incidence of anthrax has decreased annually without vaccination, and the regional distribution of anthrax gradually transferred from central and western regions to the eastern. Moreover, the occupation distribution evolved from multiple early occupations to predominated by farmers and herdsmen. This change is closely related to policy factors and to changes in the means of production and the living habits of the local population. This indicates that reformulating the control and prevention strategies is essential. Both A. Br. Ames and A. Br. 001/002 subgroups were the predominant CanSNP genotypes of Bacillus anthracis in Inner Mongolia. A total of 36 strains constituted six shared MLVA-15 genotypes, suggesting an epidemiological link between the strains of each shared genotype. The six shared genotypes ([GT1, 9, 11 and 15] and [GT8 and 12]) consisting of 2-7 strains confirmed the occurrence of multiple point outbreaks and cross-regional transmission caused by multiple common sources of infection. Phylogenetic analysis based on the WGS core genome showed that strains from this study formed an independent clade (C.V.), and they were positioned close to each other, suggesting a common origin. Further comparison analysis should be performed to ascertain the geographic origin of these strains.

Molecular Typing of Invasive Staphylococcus aureus from the Emerging Infections Program (EIP) Using Whole-Genome Sequencing

Background: The CDC has performed surveillance for invasive Staphylococcus aureus (iSA) infections through the Emerging Infections Program (EIP) since 2004. SCCmec and spa typing for clonal complex (CC) assignment and genomic markers have been used to characterize isolates. In 2019, whole-genome sequencing (WGS) of isolates began, allowing for high-resolution assessment of genomic diversity. Here, we evaluate the reliability of SCCmec typing, spa typing, and CC assignment using WGS data compared to traditional methods to ensure that backwards compatibility is maintained. Methods:S. aureus isolates were obtained from a convenience sample of iSA cases reported through the EIP surveillance system. Overall, 78 iSA isolates with diverse spa repeat patterns, CCs, SCCmec types, and antimicrobial susceptibility profiles were sequenced (MiSeq, Illumina). Real-time PCR and Sanger sequencing were used as the SCCmec and spa typing reference methods, respectively. spa-MLST mapping (Ridom SpaServer) served as the reference method for CC assignment. WGS assembly and multilocus sequence typing (MLST) were performed using the CDC QuAISAR-H pipeline. WGS-based MLST CCs were assigned using eBURST and SCCmec types using SCCmecFinder. spa types were assigned from WGS assemblies using BioNumerics. For isolate subtyping, previously published and validated canonical single-nucleotide polymorphisms (canSNPs) as well as the presence of the Panton-Valentine leukocidin (PVL) toxin and arginine catabolic mobile element (ACME) virulence factor were assessed for all genome assemblies. Results: All isolates were assigned WGS-based spa types, which were 100% concordant (78 of 78) with Sanger-based spa typing. SCCmecFinder assigned 91% of isolates (71 of 78) SCCmec types, which were 100% concordant with reference method results. Also, 7 isolates had multiple cassettes predicted or an incomplete SCCmec region assembly. Using WGS data, 96% (75 of 78) of isolates were assigned CCs; 3 isolates had unknown sequence types that were single-locus variants of established sequence types. Overall, 70 isolates had CCs assigned by the reference method; 100% (70 of 70) concordance was observed with WGS-based CCs. Analysis of canSNPs placed 42% (33 of 78) of isolates into CC8, with 17 (52%) of these isolates classified as USA300. PVL and ACME were not accurate markers for inferring the USA300 subtype as 24% (4 of 17) of isolates did not contain these markers. Conclusions:S. aureus CCs, SCCmec, and spa types can be reliably determined using WGS. Incorporation of canSNP analysis represents a more efficient method for CC8 assignment than the use of genomic markers alone. WGS allows for the replacement of multiple typing methods for increased laboratory efficiency, while maintaining backward compatibility with historical typing nomenclature.Funding: NoneDisclosures: None

Molecular characterization of B. anthracis isolates from the anthrax outbreak among cattle in Karnataka, India

BackgroundAnthrax, a zoonotic disease is caused by the Gram positive bacterium Bacillus anthracis. During January 2013, an anthrax outbreak among cattle was reported in Gundlupet Taluk, neighboring Bandipur National Park and tiger reserve, India. The present study aims at the molecular identification and characterization of 12 B. anthracis isolates from this outbreak by 16S rRNA gene sequencing, screening B. anthracis specific prophages and chromosomal markers, protective antigen (pag) gene and canonical single nucleotide polymorphism (canSNP) analysis to subtype the isolates into one of the twelve globally identified clonal sub-lineages of B. anthracis.ResultsThese isolates had identical 16S rDNA nucleotide sequences with B. anthracis specific dual peaks showing mixed base pair R (G/A) at position 1139 with visual inspection while the automated basecaller software indicated a G. Alternatively the nucleotide A at 1146 position was indicative of the 16S rDNA type 7. Multiple sequence alignment with additional 170 (16S rDNA) sequences of B. cereus sensu lato group from GenBank database revealed 28 new 16S types in addition to eleven 16S types reported earlier. The twelve B. anthracis isolates were found to harbor the four B. anthracis specific prophages (lambdaBa01, lambdaBa02, lambdaBa03, and lambdaBa04) along with its four specific loci markers (dhp 61.183, dhp 77.002, dhp 73.019, and dhp 73.017). The pag gene sequencing identified the isolates as protective antigen (PA) genotype I with phenylalanine-proline-alanine phenotype (FPA phenotype). However, sequence clustering with additional 34 pag sequences from GenBank revealed two additional missense mutations at nucleotide positions 196 bp and 869 bp of the 2294 bp pag sequence among the 5 B. cereus strains with pXO1 like plasmids. The canSNP analysis showed that the isolates belong to A.Br.Aust94 sub-lineage that is distributed geographically in countries of Asia, Africa, Europe and Australia.ConclusionsThe analysis of 16S rDNA sequences reiterated the earlier findings that visual inspection of electropherogram for position 1139 having nucleotide R could be used for B. anthracis identification and not the consensus sequence from base caller. The canSNP results indicated that the anthrax outbreak among cattle was caused by B. anthracis of A.Br.Aust94 sub-lineage.

Epidemiological Investigation and Etiological Analysis of a Cutaneous Anthrax Epidemic Caused by Butchering Sick Cattle in Guizhou, China.

A suspected human cutaneous anthrax epidemic caused by butchering sick cattle occurred in Zhijin County of Guizhou Province, Southwest of China, in 2016. Epidemiological investigation and etiological analysis were performed to provide a scientific basis for the source tracking of the epidemic. The epidemic was epidemiologically investigated, and skin blister samples collected from patients and soil samples collected from the butchering spots were used for Bacillus anthracis isolation. The suspicious B. anthracis isolates were identified using conventional methods and PCR, followed by genotyping using multiple-locus variable-number tandem repeats (VNTRs) analysis (MLVA-15) and canonical single-nucleotide polymorphism (canSNP). The genetic relationship of epidemic strains and isolates collected from other regions was analyzed. Epidemiological investigation results showed that the patients may be infected by B. anthracis during butchering sick cattle. Two suspected B. anthracis strains were isolated from blood samples and blister fluids, respectively. Conventional methods identified the two suspected isolates as B. anthracis, while PCR results showed that anti-protective antigen (PA) and capsule (CAP) gene were positive in the two isolates. MLVA-15 showed that the MLVA profiles of the two isolates were 9-20-12-53-16-2-8-8-8-4-4-4-4-10-4, which is different from the MLVA profiles of representative strains from other regions. CanSNP analysis showed that the two strains belonged to cluster A.Br.001/002. Clustering analysis and minimum spanning tree (MST) demonstrated that the two isolates were clustered with strains previously isolated from Guizhou Province. The results indicated that B. anthracis was the pathogen for this epidemic, and the patients were infected during butchering the sick. The genetic characteristics and the relationship of the B. anthracis isolates to strains from other regions indicated that the epidemic was a local occurrence.

Phylogeography and Genetic Diversity of Francisella tularensis subsp. holarctica in France (1947-2018).

In France, tularemia is caused by Francisella tularensis subsp. holarctica and is a sporadic disease affecting mainly wildlife animals and humans. F. tularensis species presents low genetic diversity that remains poorly described in France, as only a few genomes of isolates from the country are available so far. The objective of this study was to characterize the genetic diversity of F. tularensis in France and describe the phylogenetic distribution of isolates through whole-genome sequencing and molecular typing. Whole genomes of 350 strains of human or animal origin, collected from 1947 to 2018 in France and neighboring countries, were sequenced. A preliminary classification using the established canonical single nucleotide polymorphism (canSNP) nomenclature was performed. All isolates from France (except four) belonged to clade B.44, previously described in Western Europe. To increase the resolution power, a whole-genome SNP analysis was carried out. We were able to accurately reconstruct the population structure according to the global phylogenetic framework, and highlight numerous novel subclades. Whole-genome SNP analysis identified 87 new canSNPs specific to these subclades, among which 82 belonged to clade B.44. Identifying genomic features that are specific to sublineages is highly relevant in epidemiology and public health. We highlighted a large number of clusters among a single clade (B.44), which shows for the first time some genetic diversity among F. tularensis isolates from France, and the star phylogeny observed in clade B.44-subclades revealed that F. tularensis biodiversity in the country is relatively recent and resulted from clonal expansion of a single population. No association between clades and hosts or clinical forms of the disease was detected, but spatiotemporal clusters were identified for the first time in France. This is consistent with the hypothesis of persistence of F. tularensis strains found in Western Europe in the environment, associated with slow replication rates. Moreover, the presence of identical genotypes across long periods of time, and across long distances, supports this hypothesis but also suggests long-distance dispersal of the bacterium.

Genetic diversity of Bacillus anthracis Ames lineage strains in China

BackgroundAnthrax is an endemic disease that persists in the rural regions of China. The global genetic population structure of B.anthracis has also been defined by the canonical single-nucleotide polymorphisms (canSNP) and multiple-locus variable-number tandem repeat analysis (MLVA). Five canSNP lineages were found in China, and the A.Br.Ames lineage has been the second predominant group in recent years in China. The objective of this study was to reveal genetic diversity of the Ames lineage strains by MLVA.MethodsTwo molecular typing methods, canSNP and MLVA with 15markers were used to study the genetic relationship among the Ames lineage strains. The outbreak information associated with these strains was also collected and investigated.ResultsFrom 2007 to 2018, a total of 21 human anthrax infection outbreaks (68 patients) associated with B. anthracis Ames lineage strains were reported in China. Ames lineage strain-associated human anthrax is mainly distributed in the northern part of China, including the provinces of Inner Mongolia, Liaoning, Gansu, and Xinjiang. In the study, a total of 30 Ames lineage strains were included and 10 MLVA15 genotypes were identified. These strains were mainly found in northeast China, Inner Mongolia and Liaoning. In recent years, the Ames lineage strains were isolated in the two provinces every year. The 18 Ames lineage strains isolated from Inner Mongolia were divided into eight MLVA15 genotypes. From 2010 to 2015, there were continuous reports of outbreaks in Keyouzhongqi County, Inner Mongolia, and the strains that were isolated annually in succession belonged to the MLVA15–30 genotype.ConclusionsThe Ames lineage strains are widely distributed in northern China. Their genetic diversity can be illustrated by the results of the MLVA. The genetic characteristics of the Ames lineage strains from outbreaks in different provinces varied. In some areas, human anthrax outbreaks occurred annually in succession, and these related strains grouped together. These observations indicate that the local environment was persistently contaminated with B. anthracis spores, vaccination of livestock should become the fundamental control measure in the areas.

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.

Genetic characterization of Bacillus anthracis strains circulating in Italy from 1972 to 2018.

In Italy anthrax is an endemic disease, with a few outbreaks occurring almost every year. We surveyed 234 B. anthracis strains from animals (n = 196), humans (n = 3) and the environment (n = 35) isolated during Italian outbreaks in the years 1972-2018. Despite the considerable genetic homogeneity of B. anthracis, the strains were effectively differentiated using canonical single nucleotide polymorphisms (CanSNPs) assay and multiple-locus variable-number tandem repeat analysis (MLVA). The phylogenetic identity was determined through the characterization of 14 CanSNPs. In addition, a subsequent 31-loci MLVA assay was also used to further discriminate B. anthracis genotypes into subgroups. The analysis of 14 CanSNPs allowed for the identification of four main lineages: A.Br.011/009, A.Br.008/011 (respectively belonging to A.Br.008/009 sublineage, also known Trans-Eurasian or TEA group), A.Br.005/006 and B.Br.CNEVA. A.Br.011/009, the most common subgroup of lineage A, is the major genotype of B. anthracis in Italy. The MLVA analysis revealed the presence of 55 different genotypes in Italy. Most of the genotypes are genetically very similar, supporting the hypothesis that all strains evolved from a local common ancestral strain, except for two genotypes representing the branch A.Br.005/006 and B.Br.CNEVA. The genotyping analysis applied in this study remains a very valuable tool for studying the diversity, evolution, and molecular epidemiology of B. anthracis.

Phylogenetic Placement of Isolates Within the Trans-Eurasian Clade A.Br.008/009 of Bacillus anthracis.

The largest phylogenetic lineage known to date of the anthrax pathogen Bacillus anthracis is the wide-spread, so-called Trans-Eurasian clade systematically categorized as the A.Br.008/009 group sharing two defining canonical single-nucleotide polymorphisms (canSNP). In this study, we genome-sequenced a collection of 35 B. anthracis strains of this clade, derived from human infections, animal outbreaks or soil, mostly from European countries isolated between 1936 and 2008. The new data were subjected to comparative chromosomal analysis, together with 75 B. anthracis genomes available in public databases, and the relative placements of these isolates were determined within the global phylogeny of the A.Br.008/009 canSNP group. From this analysis, we have detected 3754 chromosomal SNPs, allowing the assignation of the new chromosomal sequences to established sub-clades, to define new sub-clades, such as two new Spanish, one Bulgarian or one German group(s), or to introduce orphan lineages. SNP-based results were compared with that of a multilocus variable number of tandem repeat analysis (MLVA). This analysis indicated that MLVA typing might provide additional information in cases when genomics yields identical genotypes or shows only minor differences. Introducing the delayed mismatch amplification assay (DMAA) PCR-analysis, we developed a cost-effective method to interrogate for a set of ten phylogenetically informative SNPs within genomes of A.Br.008/009 canSNP clade strains of B. anthracis. By this approach, additional 32 strains could be assigned to five of ten defined clades.