Meso-scale ecology of anthrax in southern Africa: a pilot study of diversity and clustering

It has only recently been possible to detect sufficient genetic diversity among anthrax isolates to allow genotype grouping (Keim et al. 1997). Early results of such grouping suggest that the southern African subcontinent may be the geographical origin of Bacillus anthracis. This report describes a pilot investigation of the genetic diversity of a study group of isolates from the Kruger National Park, South Africa, and efforts to detect spatio-temporal clustering within the study group. This study has also served as further validation for the newly developed Multi-Locus VNTR Analysis (MLVA), designed to simplify genotyping of B. anthracis isolates. The results reveal a diverse range of genotypes within the park allied with three genotype reference groups, and show that the MLVA procedure is robust for rapid analysis of B. anthracis genotypes. We also observed multiple genotype groups within epidemics and between geographically and temporally close epidemic episodes. This is in contrast to earlier characterizations of anthrax epidemics. The result of a Mantel test for time-space clustering indicates clustering of the anthrax isolates selected for the study.     

Population structure and evolution of the Bacillus cereus group

Representative strains of the Bacillus cereus group of bacteria, including Bacillus anthracis (11 isolates), B. cereus (38 isolates), Bacillus mycoides (1 isolate), Bacillus thuringiensis (53 isolates from 17 serovars), and Bacillus weihenstephanensis (2 isolates) were assigned to 59 sequence types (STs) derived from the nucleotide sequences of seven alleles, glpF, gmk, ilvD, pta, pur, pycA, and tpi. Comparisons of the maximum likelihood (ML) tree of the concatenated sequences with individual gene trees showed more congruence than expected by chance, indicating a generally clonal structure to the population. The STs followed two major lines of descent. Clade 1 comprised B. anthracis strains, numerous B. cereus strains, and rare B. thuringiensis strains, while clade 2 included the majority of the B. thuringiensis strains together with some B. cereus strains. Other species were allocated to a third, heterogeneous clade. The ML trees and split decomposition analysis were used to assign STs to eight lineages within clades 1 and 2. These lineages were defined by bootstrap analysis and by a preponderance of fixed differences over shared polymorphisms among the STs. Lineages were named with reference to existing designations: Anthracis, Cereus I, Cereus II, Cereus III, Kurstaki, Sotto, Thuringiensis, and Tolworthi. Strains from some B. thuringiensis serovars were wholly or largely assigned to a single ST, for example, serovar aizawai isolates were assigned to ST-15, serovar kenyae isolates were assigned to ST-13, and serovar tolworthi isolates were assigned to ST-23, while other serovars, such as serovar canadensis, were genetically heterogeneous. We suggest a revision of the nomenclature in which the lineage and clone are recognized through name and ST designations in accordance with the clonal structure of the population.     

Population structure and mouse-virulence of Toxoplasma gondii in Brazil

Recent studies found that isolates of Toxoplasma gondii from Brazil were biologically and genetically different from those in North America and Europe. However, to date only a small number of isolates have been analysed from different animal hosts in Brazil. In the present study DNA samples of 46 T. gondii isolates from cats in 11 counties in S?o Paulo state, Brazil were genetically characterised using 10 PCR restriction fragment length polymorphism markers including SAG1, SAG2, SAG3, BTUB, GRA6, c22-8, c29-2, L358, PK1 and Apico. An additional marker, CS3, that locates on chromosome VIIa and has previously been shown to be linked to acute virulence of T. gondii was also used to determine its association to virulence in mice. Genotyping of these 46 isolates revealed a high genetic diversity with 20 genotypes but no clonal Type I, II or III lineage was found. Two of the 46 isolates showed mixed infections. Combining genotyping data in this study with recent reported results from chickens, dogs and cats in Brazil (total 125 isolates) identified 48 genotypes and 26 of these genotypes had single isolates. Four of the 48 genotypes with multiple isolates identified from different hosts and locations are considered the common clonal lineages in Brazil. These lineages are designated as Types BrI, BrII, BrIII and BrIV. These results indicate that the T. gondii population in Brazil is highly diverse with a few successful clonal lineages expanded into wide geographical areas. In contrast to North America and Europe, where the Type II clonal lineage is overwhelmingly predominant, no Type II strain was identified from the 125 Brazil isolates. Analysis of mortality rates in infected mice indicates that Type BrI is highly virulent, Type BrIII is non-virulent, whilst Type BrII and BrIV lineages are intermediately virulent. In addition, allele types at the CS3 locus are strongly linked to mouse-virulence of the parasite. Thus, T. gondii has an epidemic population structure in Brazil and the major lineages have different biological traits.     

Loss of Heterozygosity Drives Clonal Diversity of Phytophthora capsici in China

Phytophthora capsici causes significant loss to pepper (Capsicum annum) in China and our goal was to develop single nucleotide polymorphism (SNP) markers for P. capsici and characterize genetic diversity nationwide. Eighteen isolates of P. capsici from locations worldwide were re-sequenced and candidate nuclear and mitochondrial SNPs identified. From 2006 to 2012, 276 isolates of P. capsici were recovered from 136 locations in 27 provinces and genotyped using 45 nuclear and 2 mitochondrial SNPs. There were two main mitochondrial haplotypes and 95 multi-locus genotypes (MLGs) identified. Genetic diversity was geographically structured with a high level of genotypic diversity in the north and on Hainan Island in the south, suggesting outcrossing contributes to diversity in these areas. The remaining areas of China are dominated by four clonal lineages that share mitochondrial haplotypes, are almost exclusively the A1 or A2 mating type and appear to exhibit extensive diversity based on loss of heterozygosity (LOH). Analysis of SNPs directly from infected peppers confirmed LOH in field populations. One clonal lineage is dominant throughout much of the country. The overall implications for long-lived genetically diverse clonal lineages amidst a widely dispersed sexual population are discussed.     

Molecular evolution and diversity in Bacillus anthracis as detected by amplified fragment length polymorphism markers.

Bacillus anthracis causes anthrax and represents one of the most molecularly monomorphic bacteria known. We have used AFLP (amplified fragment length polymorphism) DNA markers to analyze 78 B. anthracis isolates and six related Bacillus species for molecular variation. AFLP markers are extremely sensitive to even small sequence variation, using PCR and high-resolution electrophoresis to examine restriction fragments. Using this approach, we examined ca. 6.3% of the Bacillus genome for length mutations and ca. 0.36% for point mutations. Extensive variation was observed among taxa, and both cladistic and phenetic analyses were used to construct a phylogeny of B. anthracis and its closest relatives. This genome-wide analysis of 357 AFLP characters (polymorphic fragments) indicates that B. cereus and B. thuringiensis are the closest taxa to B. anthracis, with B. mycoides slightly more distant. B. subtilis, B. polymyxa, and B. stearothermophilus shared few AFLP markers with B. anthracis and were used as outgroups to root the analysis. In contrast to the variation among taxa, only rare AFLP marker variation was observed within B. anthracis, which may be the most genetically uniform bacterial species known. However, AFLP markers did establish the presence or absence of the pXO1 and pXO2 plasmids and detected 31 polymorphic chromosomal regions among the 79 B. anthracis isolates. Cluster analysis identified two very distinct genetic lineages among the B. anthracis isolates. The level of variation and its geographic distribution are consistent with a historically recent African origin for this pathogenic organism. Based on AFLP marker similarity, the ongoing anthrax epidemic in Canada and the northern United States is due to a single strain introduction that has remained stable over at least 30 years and a 1,000-mile distribution.     

Identification of strain specific markers in Bacillus anthracis by random amplification of polymorphic DNA

Classification and differentiation of Bacillus anthracis isolates by genetic markers play an important role in anthrax research. We used a PCR based method--Random Amplification of Polymorphic DNA (RAPD)--to identify genetic markers in B. anthracis strains. Twenty-five differential genetic markers were identified which divided the strains into five different groups. Three selected RAPD-markers were cloned and sequenced. The five RAPD-derived genotypes could be defined by integration of these three markers. This system offers a simple non-expensive method to classify B. anthracis strains in laboratories involved in the research of this bacterium.     

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.     

Identification of an African Bacillus anthracis lineage that lacks expression of the spore surface-associated anthrose-containing oligosaccharide

The surfaces of Bacillus anthracis endospores expose a pentasaccharide containing the monosaccharide anthrose, which has been considered for use as a vaccine or target for specific detection of the spores. In this study B. anthracis strains isolated from cattle carcasses in African countries where anthrax is endemic were tested for their cross-reactivity with monoclonal antibodies (MAbs) specific for anthrose-containing oligosaccharides. Unexpectedly, none of the isolates collected in Chad, Cameroon, and Mali were recognized by the MAbs. Sequencing of the four-gene operon encoding anthrose biosynthetic enzymes revealed the presence of premature stop codons in the aminotransferase and glycosyltransferase genes in all isolates from Chad, Cameroon, and Mali. Both immunological and genetic findings suggest that the West African isolates are unable to produce anthrose. The anthrose-deficient strains from West Africa belong to a particular genetic lineage. Immunization of cattle in Chad with a locally produced vaccine based on anthrose-positive spores of the B. anthracis strain Sterne elicited an anti-carbohydrate IgG response specific for a synthetic anthrose-containing tetrasaccharide as demonstrated by glycan microarray analysis. Competition immunoblots with synthetic pentasaccharide derivatives suggested an immunodominant role of the anthrose-containing carbohydrate in cattle. In West Africa anthrax is highly endemic. Massive vaccination of livestock in this area has taken place over long periods of time using spores of the anthrose-positive vaccine strain Sterne. The spread of anthrose-deficient strains in this region may represent an escape strategy of B. anthracis.     

Genotype diversity of the Bacillus anthracis strains isolated from the Caucasus region

The study of the genotypes of Bacillus anthracis strains isolated from the Caucasus region was performed using MLVA. Among 149 strains, 32 distinctive MLVA-8 genotypes belonging to Ala, A3a, A4 and B1 molecular diversity groups were identified. 9 genotypes were not described previously; 6 genotypes were not found in other geographic regions and could be considered as endemic for Caucasus. The majority of the identified genotypes are widespread not only in Caucasus, but also throughout Eurasia, Africa, and America. Molecular diversity of Caucasian isolates is comparable to the worldwide diversity. This represents historical relations of this region, proximity to ancient trade routes and intensity of the anthrax epizootic in Caucasus. The obtained results are of interest from the theoretical point of view, as well as for the application in epidemiological research of the anthrax outbreaks.     

Evolution of pathogenicity in the Bacillus cereus group

The Bacillus cereus group of bacteria comprises soil-dwelling saprophytes but on occasion these bacteria can cause a wide range of diseases in humans, including food poisoning, systemic infections and highly lethal forms of anthrax. While anthrax is almost invariably caused by strains from a single evolutionary lineage, Bacillus anthracis, variation in the virulence properties of strains from other lineages has not been fully addressed. Using multi-locus sequence data from 667 strains, we reconstructed the evolutionary history of the B. cereus group in terms of both clonal inheritance and recombination. The strains included 155 clinical isolates representing B. anthracis, and isolates from emetic and diarrhoeal food poisoning, septicaemia and related infections, wound, and lung infections. We confirmed the existence of three major clades and found that clinical isolates of B. cereus (with the exception of emetic toxin-producing strains) are evenly distributed between and within clades 1 and 2. B. anthracis in particular and emetic toxin-producing B. cereus show more clonal structure and are restricted to clade 1. Our characterization of the patterns of genetic exchange showed that there exist partial barriers to gene flow between the three clades. The pathogenic strains do not exhibit atypically high or low rates of recombination, consistent with the opportunistic nature of most pathogenic infections. However, there have been a large number of recent imports in clade 1 of strains from external origins, which is indicative of an on-going shift in gene-flow boundaries for this clade.     

Putative origin of clonal lineages of Amylostereum areolatum, the fungal symbiont associated with Sirex noctilio, retrieved from Pinus sylvestris, in eastern Canada

The Eurasian Sirex noctilio-Amylostereum areolatum complex was discovered and has become established close to the North American Great Lakes in the 2000s. This invasive forest insect pest represents a very high risk to native and exotic pines in North America. We investigated the geographical origin of clonal lineages of the fungal symbiont A. areolatum in the recently pest-colonized eastern Canadian region by analyzing mitochondrial and nuclear sequence variations and comparing the genetic diversity of a worldwide collection of fungal symbionts among six countries where the Sirex complex is native and four countries from which the insect-fungal complex has been introduced. In total, 102 isolates were analyzed. While 12 multilocus genotypes (MLGs) are observed in the areas where S. noctilio is native, only two MLGs are retrieved from areas where S. noctilio is not native, indicating the wide spread of clonal lineages in the introduced fungal symbiont of S. noctilio. MLG2 comprises 26% of the Canadian isolates and is also observed in Chile and South Africa, where the insect-fungal complex has also been introduced. MLG3 comprises 74% of the Canadian isolates and is also observed in the USA, but nowhere else in our worldwide collection. Thus, at least one of the Canadian clonal lineages shares a common origin with A.?areolatum isolates from the Southern Hemisphere. The source of the second clonal lineage is still unknown, but phylogenetic analyses show that MLG3 is isolated. More extended sampling is necessary to determine the origin of this fungal clonal lineage and investigate its probable symbiotic association with native North American Sirex.     

Molecular epidemiology of Bacillus anthracis: determining the correct origin

We analyzed and compared strains of Bacillus anthracis isolated from husbandry and industrial anthrax cases in Switzerland between 1952 and 1981 with published data using multiple-locus variable-number tandem repeat analysis. Strains isolated from autochthonous cases of anthrax in cattle belong to genotype B2, together with strains from continental Europe, while human B. anthracis strains clustered with genotype A4. These strains could be traced back to outbreaks of human anthrax that occurred between 1978 and 1981 in a factory processing cashmere wool from the Indian subcontinent. We interpret the worldwide occurrence of B. anthracis strains of cluster A4 to be due to the extensive global trade of untreated cashmere wool during the last century.     

Genetic characterization of Toxoplasma gondii isolates from China

Toxoplasma gondii infections are prevalent in humans and animals worldwide. In North America and Europe, T. gondii is highly clonal, consisting of three distinct lineages (Types I, II and III), whereas in South America, T. gondii is highly diverse with a few lineages expanded in the population. However, there is limited data on the diversity of T. gondii in Asia. Here we report the genetic characterization of T. gondii isolates from different hosts and geographical locations in China using the multilocus PCR–RFLP. A total of 17 T. gondii isolates from humans (3 strains), sheep (1 strain), pigs (5 strains) and cats (8 strains) were typed at 10 genetic markers including 9 nuclear loci SAG1, SAG2, SAG3, BTUB, GRA6, L358, PK1, c22-8, c29-2 and an apicoplast locus Apico. Four genotypes were revealed, including three previously reported and one new genotype. Three isolates belong to the clonal Type I lineage, one isolate belongs to the clonal Type II lineage, and the rest 13 isolates are grouped into two genotypes. This is the first report of genetic typing of T. gondii isolates from different hosts and geographical locations in China using a number of genetic markers, which has implications for the studies of population genetic structures of T. gondii, as well as for the prevention and control of T. gondii infections in humans and animals in China.     

Inferring a population structure for Staphylococcus epidermidis from multilocus sequence typing data

Despite its importance as a human pathogen, information on population structure and global epidemiology of Staphylococcus epidermidis is scarce and the relative importance of the mechanisms contributing to clonal diversification is unknown. In this study, we addressed these issues by analyzing a representative collection of S. epidermidis isolates from diverse geographic and clinical origins using multilocus sequence typing (MLST). Additionally, we characterized the mobile element (SCCmec) carrying the genetic determinant of methicillin resistance. The 217 S. epidermidis isolates from our collection were split by MLST into 74 types, suggesting a high level of genetic diversity. Analysis of MLST data using the eBURST algorithm revealed the existence of nine epidemic clonal lineages that were disseminated worldwide. One single clonal lineage (clonal complex 2) comprised 74% of the isolates, whereas the remaining isolates were clustered into 8 minor clonal lineages and 13 singletons. According to our evolutionary model, SCCmec was acquired at least 56 times by S. epidermidis. Although geographic dissemination of S. epidermidis strains and the value of the index of association between the alleles, 0.2898 (P < 0.05), support the clonality of S. epidermidis species, examination of the sequence changes at MLST loci during clonal diversification showed that recombination gives rise to new alleles approximately twice as frequently as point mutations. We suggest that S. epidermidis has a population with an epidemic structure, in which nine clones have emerged upon a recombining background and evolved quickly through frequent transfer of genetic mobile elements, including SCCmec.     

Quantitative immunofluorescence studies of the serology of Bacillus anthracis spores.

A fluorescein-conjugated antibody against formalin-inactivated spores of Bacillus anthracis Vollum reacted only weakly with a variety of Bacillus species in microfluorometric immunofluorescence assays. A conjugated antibody against spores of B. anthracis Sterne showed little affinity for spores of several B. anthracis isolates including B. anthracis Vollum, indicating that more than one anthrax spore serotype exists.     

Microsatellite evidence for high clonality and limited genetic diversity in Ziziphus celata (Rhamnaceae), an endangered, self-incompatible shrub endemic to the Lake Wales Ridge, Florida, USA

Genetic data are often crucial for designing management strategies for rare and endangered species. Ziziphus celata is an endangered sandhill shrub endemic to the Lake Wales Ridge of central Florida. This self-incompatible clonal species is known from only 14 wild populations, most of which are small (under 100 plants). Focusing on the five populations discovered in 2007, we evaluate the level of genetic diversity and identify clonal lineages within the wild populations of the species with a set of microsatellite loci. To account for somatic mutations and genotyping errors, we identified clonal lineages using a threshold cutoff for pair-wise genetic distances among samples. The microsatellites had up to 18 alleles/locus, and, consistent with outcrossing, samples were highly heterozygous (average population level H _o  = 0.69). Most populations of Z. celata consist of a single clone, and the most diverse population has only 10 clones. Overall Z. celata comprises 41 multi-locus genotypes, and 30 clonal lineages. With nearly 1,000 recorded plants (595 genotyped) and only 30 clonal lineages, Ziziphus celata is highly clonal: clonal richness, R = 0.049. The pair-wise distance method facilitates identification of clonal lineages, avoiding overestimation of clonal diversity. In most cases, the samples that grouped into a lineage were one to four plants differing from a surrounding genotype by a single microsatellite repeat insertion/deletion mutation, consistent with these having arisen via somatic mutations. Our data will enable managers to incorporate extant diversity from wild populations into ex situ collections. Additionally, our research demonstrates the utility of microsatellites for conservation of imperiled species, identifying genotypes of high priority for preservation.     

Globally invading populations of the fungal plant pathogen Verticillium dahliae are dominated by multiple divergent lineages

The spread of aggressive fungal pathogens into previously non‐endemic regions is a major threat to plant health and food security. Analyses of the spatial and genetic structure of plant pathogens offer valuable insights into their origin, dispersal mechanisms and evolution, and have been useful to develop successful disease management strategies. Here, we elucidated the genetic diversity, population structure and demographic history of worldwide invasion of the ascomycete Verticillium dahliae, a soil‐borne pathogen, using a global collection of 1100 isolates from multiple plant hosts and countries. Seven well‐differentiated genetic clusters were revealed through discriminant analysis of principal components (DAPC), but no strong associations between these clusters and host/geographic origin of isolates were found. Analyses of clonal evolutionary relationships among multilocus genotypes with the eBURST algorithm and analyses of genetic distances revealed that genetic clusters represented several ancient evolutionary lineages with broad geographic distribution and wide host range. Comparison of different scenarios of demographic history using approximate Bayesian computations revealed the branching order among the different genetic clusters and lineages. The different lineages may represent incipient species, and this raises questions with respect to their evolutionary origin and the factors allowing their maintenance in the same areas and same hosts without evidence of admixture between them. Based on the above findings and the biology of V. dahliae, we conclude that anthropogenic movement has played an important role in spreading V. dahliae lineages. Our findings have implications for the development of management strategies such as quarantine measures and crop resistance breeding.     

Divergence, hybridization, and recombination in the mitochondrial genome of the human pathogenic yeast Cryptococcus gattii

The inheritance of mitochondrial genes and genomes are uniparental in most sexual eukaryotes. This pattern of inheritance makes mitochondrial genomes in natural populations effectively clonal. Here, we examined the mitochondrial population genetics of the emerging human pathogenic fungus Cryptococcus gattii . The DNA sequences for five mitochondrial DNA fragments were obtained from each of 50 isolates belonging to two evolutionary divergent lineages, VGI and VGII. Our analyses revealed a greater sequence diversity within VGI than that within VGII, consistent with observations of the nuclear genes. The combined analyses of all five gene fragments indicated significant divergence between VGI and VGII. However, the five individual genealogies showed different relationships among the isolates, consistent with recent hybridization and mitochondrial gene transfer between the two lineages. Population genetic analyses of the multilocus data identified evidence for predominantly clonal mitochondrial population structures within both lineages. Interestingly, there were clear signatures of recombination among mitochondrial genes within the VGII lineage. Our analyses suggest historical mitochondrial genome divergence within C. gattii , but there is evidence for recent hybridization and recombination in the mitochondrial genome of this important human yeast pathogen.     

Generation and characterization of monoclonal antibodies to protective antigen of Bacillus anthracis

Monoclonal antibodies (MoAbs) were generated following immunization of BALB/c mice with protective antigen (PA) of B. anthracis. Five clones reactive to this protein were stabilized and preserved. These MoAbs could detect nanogram levels of PA when tested in ELISA. In Western blotting, they reacted with all PA preparations tested and no cross-reactivity was observed with lethal factor, edema factor of B. anthracis and with other organisms. These MoAbs could detect PA from 22 confirmed clinical isolates of B. anthracis on Western blotting and hold promise for direct detection of PA in clinical samples for diagnosing anthrax.