Molecular Characterization of Korean Bacillus anthracis Isolates by Amplified Fragment Length Polymorphism Analysis and Multilocus Variable-Number Tandem Repeat Analysis

We analyzed the genetic relationships and molecular characteristics of 34 Bacillus anthracis isolates from soil and clinical samples in various regions of Korea and 17 related Bacillus species, using the amplified fragment length polymorphism (AFLP) and multilocus variable-number tandem repeat (MLVA) approaches. Triplicate AFLP profiles of these strains showed high reproducibility and identified 376 polymorphisms. AFLP phylogenetic analysis of B. anthracis isolates showed a high level of similarity, 0.93, and this monomorphic fragment profile proved to be useful to differentiate B. anthracis strains from other Bacillus species. The B. cereus group was separated from other Bacillus species at a level of similarity of 0.68. Among them, some B. cereus strains showed genetic interspersion with B. thuringiensis strains. The evolutionary pattern of nucleotide differences among B. anthracis strains with the eight MLVA markers showed nine MLVA types. Three MLVA types, M1 to M3, were pathogenic B. anthracis isolates and were assigned as new genotypes belonging to the A4 and B3 clusters, compared with 89 genotypes deduced from previous data. This indicates that differences in cluster prevalence and distribution may be influenced more by MLVA markers on two plasmids loci and human activity. Consequently, we suggest that the novel MLVA type may represent significant evidence for historic adaptation to environmental conditions of the Asian continent, particularly Korea. Therefore, MLVA techniques may be available for molecular monitoring on anthrax-release-related bioterrorism and further study is required for the continuous epidemiological study of variable anthrax collections.     

Use of multiple‐locus variable number tandem repeat analysis and phage typing for subtyping ofSalmonella Enteritidis from sporadic human cases in the United States

Aims: To investigate the genetic diversity among S. Enteritidis isolates from different geographic regions to evaluate the relationship between phage types (PTs) and variable number tandem repeat analysis (VNTR) loci. Methods and Results: We performed multiple‐locus variable number tandem repeat analysis (MLVA) and phage typing on 245 S. Enteritidis isolates collected from sporadic human clinical cases in Michigan, Minnesota, New York, and Washington states between 2000 and 2007. Ninety‐four MLVA types and 22 different PTs were identified. Specific PTs were associated with a predominant allele for certain VNTR loci. Cluster analysis using a minimum‐spanning tree demonstrated two major clusters (I, II) and one minor cluster of isolates. PTs 8, 13a, 13 and 34 were significantly associated with MLVA cluster I. Phage types 1, 4, 6a, and 18 were significantly associated with MLVA cluster II. Conclusions: We found significant association between MLVA‐based clusters and PTs. Certain VNTR loci were associated with specific PTs and could serve as useful molecular markers for S. Enteritidis in epidemiological investigations. Significance and Impact of the Study: MLVA genotyping in combination with phage typing can be used for effective characterization of S. Enteritidis isolates. It can also be useful for tracing possible sources during investigations of sporadic and outbreak cases of S. Enteritidis.     

Distribution and molecular evolution of bacillus anthracis genotypes in Namibia

The recent development of genetic markers for Bacillus anthracis has made it possible to monitor the spread and distribution of this pathogen during and between anthrax outbreaks. In Namibia, anthrax outbreaks occur annually in the Etosha National Park (ENP) and on private game and livestock farms. We genotyped 384 B. anthracis isolates collected between 1983-2010 to identify the possible epidemiological correlations of anthrax outbreaks within and outside the ENP and to analyze genetic relationships between isolates from domestic and wild animals. The isolates came from 20 animal species and from the environment and were genotyped using a 31-marker multi-locus-VNTR-analysis (MLVA) and, in part, by twelve single nucleotide polymorphism (SNP) markers and four single nucleotide repeat (SNR) markers. A total of 37 genotypes (GT) were identified by MLVA, belonging to four SNP-groups. All GTs belonged to the A-branch in the cluster- and SNP-analyses. Thirteen GTs were found only outside the ENP, 18 only within the ENP and 6 both inside and outside. Genetic distances between isolates increased with increasing time between isolations. However, genetic distance between isolates at the beginning and end of the study period was relatively small, indicating that while the majority of GTs were only found sporadically, three genetically close GTs, accounting for more than four fifths of all the ENP isolates, appeared dominant throughout the study period. Genetic distances among isolates were significantly greater for isolates from different host species, but this effect was small, suggesting that while species-specific ecological factors may affect exposure processes, transmission cycles in different host species are still highly interrelated. The MLVA data were further used to establish a model of the probable evolution of GTs within the endemic region of the ENP. SNR-analysis was helpful in correlating an isolate with its source but did not elucidate epidemiological relationships.     

Insights to Genetic Characterization Tools for Epidemiological Tracking of Francisella tularensis in Sweden

Tularaemia, caused by the bacterium Francisella tularensis, is endemic in Sweden and is poorly understood. The aim of this study was to evaluate the effectiveness of three different genetic typing systems to link a genetic type to the source and place of tularemia infection in Sweden. Canonical single nucleotide polymorphisms (canSNPs), MLVA including five variable number of tandem repeat loci and PmeI-PFGE were tested on 127 F. tularensis positive specimens collected from Swedish case-patients. All three typing methods identified two major genetic groups with near-perfect agreement. Higher genetic resolution was obtained with canSNP and MLVA compared to PFGE; F. tularensis samples were first assigned into ten phylogroups based on canSNPs followed by 33 unique MLVA types. Phylogroups were geographically analysed to reveal complex phylogeographic patterns in Sweden. The extensive phylogenetic diversity found within individual counties posed a challenge to linking specific genetic types with specific geographic locations. Despite this, a single phylogroup (B.22), defined by a SNP marker specific to a lone Swedish sequenced strain, did link genetic type with a likely geographic place. This result suggests that SNP markers, highly specific to a particular reference genome, may be found most frequently among samples recovered from the same location where the reference genome originated. This insight compels us to consider whole-genome sequencing (WGS) as the appropriate tool for effectively linking specific genetic type to geography. Comparing the WGS of an unknown sample to WGS databases of archived Swedish strains maximizes the likelihood of revealing those rare geographically informative SNPs.     

Bacillus anthracis Diversity and Geographic Potential across Nigeria,Cameroon and Chad: Further Support of a Novel West African Lineage

Zoonoses, diseases affecting both humans and animals, can exert tremendous pressures on human and veterinary health systems, particularly in resource limited countries. Anthrax is one such zoonosis of concern and is a disease requiring greater public health attention in Nigeria. Here we describe the genetic diversity of Bacillus anthracis in Nigeria and compare it to Chad, Cameroon and a broader global dataset based on the multiple locus variable number tandem repeat (MLVA-25) genetic typing system. Nigerian B. anthracis isolates had identical MLVA genotypes and could only be resolved by measuring highly mutable single nucleotide repeats (SNRs). The Nigerian MLVA genotype was identical or highly genetically similar to those in the neighboring countries, confirming the strains belong to this unique West African lineage. Interestingly, sequence data from a Nigerian isolate shares the anthrose deficient genotypes previously described for strains in this region, which may be associated with vaccine evasion. Strains in this study were isolated over six decades, indicating a high level of temporal strain stability regionally. Ecological niche models were used to predict the geographic distribution of the pathogen for all three countries. We describe a west-east habitat corridor through northern Nigeria extending into Chad and Cameroon. Ecological niche models and genetic results show B. anthracis to be ecologically established in Nigeria. These findings expand our understanding of the global B. anthracis population structure and can guide regional anthrax surveillance and control planning.     

Comparison of Minisatellite Polymorphisms in the Bacillus cereus Complex: a Simple Assay for Large-Scale Screening and Identification of Strains Most Closely Related to Bacillus anthracis

Polymorphism of five tandem repeats that are monomorphic in Bacillus anthracis was investigated in 230 isolates of the B. cereus group and in 5 sequenced B. cereus genomes in search for markers allowing identification of B. cereus and B. thuringiensis strains most closely related to B. anthracis. Using this multiple-locus variable number of tandem repeat analysis (MLVA), a cluster of 30 strains was selected for further characterization. Eventually, six of these were characterized by multilocus sequence type analysis. One of the strains is only six point mutations (of almost 3,000 bp) away from B. anthracis and was also proposed to be closest to B. anthracis by MLVA analysis. However, this strain remains separated from B. anthracis by a number of significant genetic events observed in B. anthracis, including the loss of the hemolysin activity, the presence of four prophages, and the presence of the two virulence plasmids, pXO1 and pXO2. One particular minisatellite marker provides an efficient assay to identify the subset of B. cereus and B. thuringiensis strains closely related to B. anthracis. Based on these results, a very simple assay is proposed that allows the screening of hundreds of strains from the B. cereus complex, with modest equipment and at a low cost, to eventually fill the gap with B. anthracis and better understand the origin and making of this dangerous pathogen.     

Detection of virulent Escherichia coli O157 strains using multiplex PCR and single base sequencing for SNP characterization

Aims: To compare 167 Norwegian human and nonhuman Escherichia coli O157:H7/NM (nonmotile) isolates with respect to an A/T single nucleotide polymorphism (SNP) in the tir gene and to detect specific SNPs that differentiate STEC O157 into distinct virulence clades (1–3 and 8). Methods and Results: We developed a multiplex PCR followed by single base sequencing for detection of the SNPs, and examined the association among SNP genotype, virulence profile (stx and eae status), multilocus variable number of tandem repeats analysis (MLVA) profile and clinical outcome. We found an over‐representation of the T allele among human strains compared to nonhuman strains, including 5/6 haemolytic‐uraemic syndrome cases. Fourteen strains belonged to clade 8, followed by two clade 2 strains. No clade 1 nor 3 isolates were observed. stx1 in combination with either stx2_EDL933 or stx2c were frequently observed among human strains, whereas stx2c was dominating in nonhuman strains. MLVA indicated that only single cases or small outbreaks with E. coli O157 have been observed in Norway through the years 1993–2008. Conclusion: We observed that the tir‐255 A/T SNP and the stx status were different between human and nonhuman O157 strains. No major outbreaks were observed, and only a few strains were differentiated into the virulence clades 2 and 8. Significance and Impact of the Study: The detection of virulence clade‐specific SNPs enables the rapid designation of virulent E. coli O157 strains, especially in outbreak situations.     

Development of polymorphic microsatellite and single nucleotide polymorphism markers for Cercospora beticola (Mycosphaerellaceae)

The plant pathogenic fungus, Cercospora beticola, causes the most important foliage disease of sugar beet. A previous study has shown that isolates of opposite mating types are present in equal proportions in natural populations; therefore, the aim of this study was to develop highly reproducible polymorphic markers for analysing populations of C. beticola. Five microsatellite and four single nucleotide polymorphism (SNP) markers were developed that allow rapid screening of genetic diversity in C. beticola. Six populations were screened with these markers and all were found to be in gametic equilibrium, indicating random mating in C. beticola.     

Pathogenic potential and horizontal gene transfer in ovine gastrointestinal Escherichia coli

Aims: To demonstrate that a thorough characterization and virulotyping of Escherichia coli strains isolated from sheep over time leads to new insights into ovine E. coli potentially becoming human pathogens through horizontal gene transfer. Methods and Results: One hundred and fifty E. coli isolates from two sheep, sampled over 3 weeks, were characterized by serotyping, virulotyping, genotyping using multiple locus variable number tandem repeats analysis (MLVA) and susceptibility to phage infection in vitro. The 35 MLVA profiles and the serotype and virulotypes of the strains were closely associated. Many MLVA profiles differed in one locus independent of serotypes. Escherichia coli isolates of the same serotype or virulotype had identical or very similar MLVA profiles. No transductants that incorporated the bacteriophages were found in vivo, but six E. coli isolates were susceptible to the phage infection in vitro. Changes in MLVA profiles were seen after acquisition of Stx phages in vitro only. Conclusions: The sheep carried Stx phage susceptible E. coli that possessed virulence markers associated with human pathogenicity. Changes in bacterial genomes by phage transfer may complicate outbreak source investigations. Serotype has to be taken into account when evaluating strain relationships by MLVA. Significance and Impact of the Study: Sheep carry E. coli that encode for virulence markers and belong to serogroups known to be human pathogens. In addition, a selection of isolates was found to be susceptible to horizontal transfer of Shiga toxin genes by means of bacteriophages in vitro, and the transfer resulted in a discernible change of the MLVA patterns of E. coli.     

SNR analysis: molecular investigation of an anthrax epidemic

Background  

In Italy, anthrax is endemic but occurs sporadically. During the summer of 2004, in the Pollino National Park, Basilicata, Southern Italy, an anthrax epidemic consisting of 41 outbreaks occurred; it claimed the lives of 124 animals belonging to different mammal species. This study is a retrospective molecular epidemiological investigation carried out on 53 isolates collected during the epidemic. A 25-loci Multiple Locus VNTR Analysis (MLVA) MLVA was initially performed to define genetic relationships, followed by an investigation of genetic diversity between epidemic strains through Single Nucleotide Repeat (SNR) analysis.     

Use of high-resolution melting and melting temperature-shift assays for specific detection and identification of Bacillus anthracis based on single nucleotide discrimination

Single nucleotide polymorphisms (SNPs) are important diagnostic markers for the detection and differentiation of Bacillus anthracis. High-Resolution Melting (HRM) and Melting Temperature (Tm)-shift methods are two approaches that enable SNP detection without the need for expensive labeled probes. We evaluated the potential diagnostic capability of those methods to discriminate B. anthracis from the other members of the B. cereus group. Two assays targeting B. anthracis-specific SNPs in the plcR and gyrA genes were designed for each method and used to genotype a panel of 155 Bacilli strains. All B. anthracis isolates (n = 65) were correctly and unambiguously identified. Assays also proved to be appropriate for the direct genotyping of biological samples. They could reliably detect B. anthracis in contaminated organs containing as little as 10³ CFU/ml, corresponding to a few genome equivalents per reaction. The HRM and Tm-shift applications described here represent valuable tools for specific identification of B. anthracis at reduced cost.     

Molecular characterization of Korean Bacillus anthracis isolates by amplified fragment length polymorphism analysis and multilocus variable-number tandem repeat analysis

We analyzed the genetic relationships and molecular characteristics of 34 Bacillus anthracis isolates from soil and clinical samples in various regions of Korea and 17 related Bacillus species, using the amplified fragment length polymorphism (AFLP) and multilocus variable-number tandem repeat (MLVA) approaches. Triplicate AFLP profiles of these strains showed high reproducibility and identified 376 polymorphisms. AFLP phylogenetic analysis of B. anthracis isolates showed a high level of similarity, 0.93, and this monomorphic fragment profile proved to be useful to differentiate B. anthracis strains from other Bacillus species. The B. cereus group was separated from other Bacillus species at a level of similarity of 0.68. Among them, some B. cereus strains showed genetic interspersion with B. thuringiensis strains. The evolutionary pattern of nucleotide differences among B. anthracis strains with the eight MLVA markers showed nine MLVA types. Three MLVA types, M1 to M3, were pathogenic B. anthracis isolates and were assigned as new genotypes belonging to the A4 and B3 clusters, compared with 89 genotypes deduced from previous data. This indicates that differences in cluster prevalence and distribution may be influenced more by MLVA markers on two plasmids loci and human activity. Consequently, we suggest that the novel MLVA type may represent significant evidence for historic adaptation to environmental conditions of the Asian continent, particularly Korea. Therefore, MLVA techniques may be available for molecular monitoring on anthrax-release-related bioterrorism and further study is required for the continuous epidemiological study of variable anthrax collections.     

MLVA genotyping of Chinese human Brucella melitensis biovar 1, 2 and 3 isolates

Background

Since 1950, Brucella melitensis has been the predominant strain associated with human brucellosis in China. In this study we investigated the genotypic characteristics of B. melitensis isolates from China using a multiple-locus variable-number tandem-repeat analysis (MLVA) and evaluated the utility of MLVA with regards to epidemiological trace-back investigation.

Results

A total of 105 B. melitensis strains isolated from throughout China were divided into 69 MLVA types using MLVA-16. Nei's genetic diversity indices for the various loci ranged between 0.00 - 0.84. 12 out 16 loci were the low diversity with values < 0.2 and the most discriminatory markers were bruce16 and bruce30 with a diversity index of > 0.75 and containing 8 and 7 alleles, respectively. Many isolates were single-locus or double-locus variants of closely related B. melitensis isolates from different regions, including the north and south of China. Using panel 1, the majority of strains (84/105) were genotype 42 clustering to the 'East Mediterranean' B. melitensis group. Chinese B. melitensis are classified in limited number of closely related genotypes showing variation mainly at the panel 2B loci.

Conclusion

The MLVA-16 assay can be useful to reveal the predominant genotypes and strain relatedness in endemic or non-endemic regions of brucellosis. However it is not suitable for biovar differentiation of B. melitensis. Genotype 42 is widely distributed throughout China during a long time. Bruce 16 and bruce 30 in panel 2B markers are most useful for typing Chinese isolates.     

Genotypic and PFGE/MLVA analyses of Vibrio cholerae O1: geographical spread and temporal changes during the 2007-2010 cholera outbreaks in Thailand

Background

Vibrio cholerae O1 El Tor dominated the seventh cholera pandemic which occurred in the 1960s. For two decades, variants of V. cholerae O1 El Tor that produce classical cholera toxin have emerged and spread globally, replacing the prototypic El Tor biotype. This study aims to characterize V. cholerae O1 isolates from outbreaks in Thailand with special reference to genotypic variations over time.

Methods/Findings

A total of 343 isolates of V. cholerae O1 from cholera outbreaks from 2007 to 2010 were investigated, and 99.4% were found to carry the classical cholera toxin B subunit (ctxB) and El Tor rstR genes. Pulsed-field gel electrophoresis (PFGE) differentiated the isolates into 10 distinct pulsotypes, clustered into two major groups, A and B, with an overall similarity of 88%. Ribotyping, multiple-locus variable-number tandem-repeat analysis (MLVA), and PCR to detect Vibrio seventh pandemic island II (VSP-II) related genes of randomly selected isolates from each pulsotype corresponded to the results obtained by PFGE. Epidemiological investigations revealed that MLVA type 2 was strongly associated with a cholera outbreak in northeastern Thailand in 2007, while MLVA type 7 dominated the outbreaks of the southern Gulf areas in 2009 and MLVA type 4 dominated the outbreaks of the central Gulf areas during 2009–2010. Only MLVA type 16 isolates were found in a Thai-Myanmar border area in 2010, whereas those of MLVA types 26, 39, and 41 predominated this border area in 2008. Type 39 then disappeared 1–2 years later as MLVA type 41 became prevalent. Type 41 was also found to infect an outbreak area.

Conclusions

MLVA provided a high-throughput genetic typing tool for understanding the in-depth epidemiology of cholera outbreaks. Our epidemiological surveys suggest that some clones of V. cholerae O1 with similar but distinctive genetic traits circulate in outbreak sites, while others disappear over time.     

Streptococcus suis outbreak investigation using multiple‐locus variable tandem repeat number analysis

Two outbreaks of Streptococcus suis ST7 occurred in humans in 1998 and 2005 in China. PFGE of chromosome restriction fragments found all ST7 isolates to be indistinguishable. Due to the genetic homogeneity of ST7 isolates, development of a rapid sub‐typing method with high discriminatory power for ST7 isolates is required. In this study, a novel method, MLVA, was developed to type S. suis serotype 2 strains. Further, this method was used to analyze outbreak‐associated ST7 strains in China. A total of 144 ST7 S. suis isolates were sub‐typed into 34 MLVA types. Among these, eight isolates from the 1998 outbreak were sub‐typed into five MLVA types, of which four MLVA types were also detected in Sichuan in 2005. These data indicate that the pathogens responsible for the two outbreaks had the same origin. In addition, some observations also provided molecular evidence for the transmission route, possibly indicating that the MLVA method has usefulness in epidemiology. The developed MLVA scheme for S. suis has greater discriminative power than PFGE. The method described here may be useful for identifying the source of S. suis infection and monitoring its spread.     

MLVA Genotyping of Brucella melitensis and Brucella abortus Isolates from Different Animal Species and Humans and Identification of Brucella suis Vaccine Strain S2 from Cattle in China

In China, brucellosis is an endemic disease and the main sources of brucellosis in animals and humans are infected sheep, cattle and swine. Brucella melitensis (biovars 1 and 3) is the predominant species, associated with sporadic cases and outbreak in humans. Isolates of B. abortus, primarily biovars 1 and 3, and B. suis biovars 1 and 3 are also associated with sporadic human brucellosis. In this study, the genetic profiles of B. melitensis and B. abortus isolates from humans and animals were analyzed and compared by multi-locus variable-number tandem-repeat analysis (MLVA). Among the B. melitensis isolates, the majority (74/82) belonged to MLVA8 genotype 42, clustering in the ‘East Mediterranean’ group. Two B. melitensis biovar 1 genotype 47 isolates, belonging to the ‘Americas’ group, were recovered; both were from the Himalayan blue sheep (Pseudois nayaur, a wild animal). The majority of B. abortus isolates (51/70) were biovar 3, genotype 36. Ten B. suis biovar 1 field isolates, including seven outbreak isolates recovered from a cattle farm in Inner Mongolia, were genetically indistinguishable from the vaccine strain S2, based on MLVA cluster analysis. MLVA analysis provided important information for epidemiological trace-back. To the best of our knowledge, this is the first report to associate Brucella cross-infection with the vaccine strain S2 based on molecular comparison of recovered isolates to the vaccine strain. MLVA typing could be an essential assay to improve brucellosis surveillance and control programs.     

Geographical distribution of genotypic and phenotypic markers among Bacillus anthracis isolates and related species by historical movement and horizontal transfer

The geographical distribution of Bacillus anthracis strains and isolates bearing some of the same genetic markers as the Amerithrax Ames isolate was examined and evaluated. At least one mechanism for the horizontal movement of genetic markers was shown amongst isolates and closely related species and the effect of such mixing was demonstrated on phenotype. The results provided potential mechanisms by which attempts to attribute isolates of Bacillus anthracis to certain geographical and isolate sources may be disrupted.     

MLVA16 Typing of Portuguese Human and Animal Brucella melitensis and Brucella abortus Isolates

To investigate the epidemiological relationship of isolates from different Portuguese geographical regions and to assess the diversity among isolates, the MLVA16(Orsay) assay (panels 1, 2A and 2B) was performed with a collection of 126 Brucella melitensis (46 human and 80 animal isolates) and 157 B. abortus field isolates, seven vaccine strains and the representative reference strains of each species. The MLVA16(Orsay) showed a similar high discriminatory power (HGDI 0.972 and 0.902) for both species but panel 1 and 2A markers displayed higher diversity (HGDI 0.693) in B. abortus compared to B. melitensis isolates (HGDI 0.342). The B. melitensis population belong to the "Americas" (17%) and "East Mediterranean" (83%) groups. No isolate belonged to the "West Mediterranean" group. Eighty-five percent of the human isolates (39 in 46) fit in the "East-Mediterranean" group where a single lineage known as MLVA11 genotype 116 is responsible for the vast majority of Brucella infections in humans. B. abortus isolates formed a consistent group with bv1 and bv3 isolates in different clusters. Four MLVA11 genotypes were observed for the first time in isolates from S. Jorge and Terceira islands from Azores. From the collection of isolates analysed in this study we conclude that MLVA16(Orsay) provided a clear view of Brucella spp. population, confirming epidemiological linkage in outbreak investigations. In particular, it suggests recent and ongoing colonisation of Portugal with one B. melitensis lineage usually associated with East Mediterranean countries.     

Molecular typing and antimicrobial resistance of <Emphasis Type="Italic">Salmonella</Emphasis> Enteritidis isolated from poultry,food, and humans in Serbia

Molecular typing and resistotyping coupled with gyrA single nucleotide polymorphism (SNP) of 60 Salmonella Enteritidis (SE) isolates originated from poultry, food, and humans in Serbia is described. Molecular fingerprinting was performed by randomly amplified polymorphic DNA (RAPD) using four primers, and the diversity index (D) was 0.688. In combination with resistotyping and gyrA SNP, D increased to 0.828. A total of 23 genetic groups were obtained. When four RAPD primers were combined, epidemic isolates from a fast-food restaurant outbreak were clustered in a distinctive genetic group. Among 60 SE strains, three had multiple resistances to three or more antibiotics. Nine strains were resistant to nalidixic acid (NAL; a non-fluorinated quinolone). The mutations in quinolone resistance-determining region (QRDR) found in NAL-resistant strains were attributed to Asp⁸⁷ → Asn in six strains, Asp⁸⁷ → Gly in one strain, and Ser⁸³ → Phe in one strain. One NAL-resistant strain had no mutations in QRDR, suggesting another mechanism of resistance.