James M. Schupp

Multiple-Locus Variable-Number Tandem Repeat Analysis Reveals Genetic Relationships within Bacillus anthracis

Bacillus anthracis is one of the most genetically homogeneous pathogens described, making strain discrimination particularly difficult. In this paper, we present a novel molecular typing system based on rapidly evolving variable-number tandem repeat (VNTR) loci. Multiple-locus VNTR analysis (MLVA) uses the combined power of multiple alleles at several marker loci. In our system, fluorescently labeled PCR primers are used to produce PCR amplification products from eight VNTR regions in the B. anthracis genome. These are detected and their sizes are determined using an ABI377 automated DNA sequencer. Five of these eight loci were discovered by sequence characterization of molecular markers (vrrC(1), vrrC(2), vrrB(1), vrrB(2), and CG3), two were discovered by searching complete plasmid nucleotide sequences (pXO1-aat and pXO2-at), and one was known previously (vrrA). MLVA characterization of 426 B. anthracis isolates identified 89 distinct genotypes. VNTR markers frequently identified multiple alleles (from two to nine), with Neis diversity values between 0.3 and 0.8. Unweighted pair-group method arithmetic average cluster analysis identified six genetically distinct groups that appear to be derived from clones. Some of these clones show worldwide distribution, while others are restricted to particular geographic regions. Human commerce doubtlessly has contributed to the dispersal of particular clones in ancient and modern times.

Staphylococcus aureus CC398: Host Adaptation and Emergence of Methicillin Resistance in Livestock

ABSTRACT Since its discovery in the early 2000s, methicillin-resistant Staphylococcus aureus (MRSA) clonal complex 398 (CC398) has become a rapidly emerging cause of human infections, most often associated with livestock exposure. We applied whole-genome sequence typing to characterize a diverse collection of CC398 isolates (n = 89), including MRSA and methicillin-susceptible S. aureus (MSSA) from animals and humans spanning 19 countries and four continents. We identified 4,238 single nucleotide polymorphisms (SNPs) among the 89 core genomes. Minimal homoplasy (consistency index = 0.9591) was detected among parsimony-informative SNPs, allowing for the generation of a highly accurate phylogenetic reconstruction of the CC398 clonal lineage. Phylogenetic analyses revealed that MSSA from humans formed the most ancestral clades. The most derived lineages were composed predominantly of livestock-associated MRSA possessing three different staphylococcal cassette chromosome mec element (SCCmec) types (IV, V, and VII-like) including nine subtypes. The human-associated isolates from the basal clades carried phages encoding human innate immune modulators that were largely missing among the livestock-associated isolates. Our results strongly suggest that livestock-associated MRSA CC398 originated in humans as MSSA. The lineage appears to have undergone a rapid radiation in conjunction with the jump from humans to livestock, where it subsequently acquired tetracycline and methicillin resistance. Further analyses are required to estimate the number of independent genetic events leading to the methicillin-resistant sublineages, but the diversity of SCCmec subtypes is suggestive of strong and diverse antimicrobial selection associated with food animal production. IMPORTANCE Modern food animal production is characterized by densely concentrated animals and routine antibiotic use, which may facilitate the emergence of novel antibiotic-resistant zoonotic pathogens. Our findings strongly support the idea that livestock-associated MRSA CC398 originated as MSSA in humans. The jump of CC398 from humans to livestock was accompanied by the loss of phage-carried human virulence genes, which likely attenuated its zoonotic potential, but it was also accompanied by the acquisition of tetracycline and methicillin resistance. Our findings exemplify a bidirectional zoonotic exchange and underscore the potential public health risks of widespread antibiotic use in food animal production. Modern food animal production is characterized by densely concentrated animals and routine antibiotic use, which may facilitate the emergence of novel antibiotic-resistant zoonotic pathogens. Our findings strongly support the idea that livestock-associated MRSA CC398 originated as MSSA in humans. The jump of CC398 from humans to livestock was accompanied by the loss of phage-carried human virulence genes, which likely attenuated its zoonotic potential, but it was also accompanied by the acquisition of tetracycline and methicillin resistance. Our findings exemplify a bidirectional zoonotic exchange and underscore the potential public health risks of widespread antibiotic use in food animal production.

Population Genetics of Vibrio cholerae from Nepal in 2010: Evidence on the Origin of the Haitian Outbreak

ABSTRACT Cholera continues to be an important cause of human infections, and outbreaks are often observed after natural disasters, such as the one following the 2010 earthquake in Haiti. Once the cholera outbreak was confirmed, rumors spread that the disease was brought to Haiti by a battalion of Nepalese soldiers serving as United Nations peacekeepers. This possible connection has never been confirmed. We used whole-genome sequence typing (WGST), pulsed-field gel electrophoresis (PFGE), and antimicrobial susceptibility testing to characterize 24 recent Vibrio cholerae isolates from Nepal and evaluate the suggested epidemiological link with the Haitian outbreak. The isolates were obtained from 30 July to 1 November 2010 from five different districts in Nepal. We compared the 24 genomes to 10 previously sequenced V. cholerae isolates, including 3 from the Haitian outbreak (began July 2010). Antimicrobial susceptibility and PFGE patterns were consistent with an epidemiological link between the isolates from Nepal and Haiti. WGST showed that all 24 V. cholerae isolates from Nepal belonged to a single monophyletic group that also contained isolates from Bangladesh and Haiti. The Nepalese isolates were divided into four closely related clusters. One cluster contained three Nepalese isolates and three Haitian isolates that were almost identical, with only 1- or 2-bp differences. Results in this study are consistent with Nepal as the origin of the Haitian outbreak. This highlights how rapidly infectious diseases might be transmitted globally through international travel and how public health officials need advanced molecular tools along with standard epidemiological analyses to quickly determine the sources of outbreaks. IMPORTANCE Cholera is one of the ancient classical diseases and particularly prone to cause major outbreaks following major natural disasters, such as earthquakes and hurricanes, where the normal separation between sewage and drinking water is destroyed. This was the case following the 2010 earthquake in Haiti. Rumors spread that the disease was brought to Haiti by a battalion of Nepalese soldiers serving as United Nations peacekeepers. This possible connection has never been confirmed. Sequencing the genomes of bacteria can give detailed information on whether isolates from different sites share a common origin. We used this technology to sequence isolates of Vibrio cholerae from Nepal, identify single-nucleotide polymorphisms (SNPs), and compare these high-resolution genotypes to the complete genome sequences of isolates from the Haiti outbreak. We provide support for the hypothesis that the isolates were brought to Haiti from Nepal. Cholera is one of the ancient classical diseases and particularly prone to cause major outbreaks following major natural disasters, such as earthquakes and hurricanes, where the normal separation between sewage and drinking water is destroyed. This was the case following the 2010 earthquake in Haiti. Rumors spread that the disease was brought to Haiti by a battalion of Nepalese soldiers serving as United Nations peacekeepers. This possible connection has never been confirmed. Sequencing the genomes of bacteria can give detailed information on whether isolates from different sites share a common origin. We used this technology to sequence isolates of Vibrio cholerae from Nepal, identify single-nucleotide polymorphisms (SNPs), and compare these high-resolution genotypes to the complete genome sequences of isolates from the Haiti outbreak. We provide support for the hypothesis that the isolates were brought to Haiti from Nepal.

Identification and Characterization of Variable-Number Tandem Repeats in the Yersinia pestis Genome

ABSTRACT Yersinia pestis, the infamous plague-causing pathogen, appears to have emerged in relatively recent history. Evidence of this fact comes from several studies that document a lack of nucleotide diversity in the Y. pestisgenome. In contrast, we report that variable-number tandem repeat (VNTR) sequences are common in the Y. pestis genome and occur frequently in gene coding regions. Larger tandem repeat arrays, most useful for phylogenetic analysis, are present at an average of 2.18 arrays per 10 kbp and are distributed evenly throughout the genome and the two virulence plasmids, pCD1 and pMT1. We examined allelic diversity at 42 chromosomal VNTR loci in 24 selected isolates (12 globally distributed and 12 from Siskiyou County, Calif.). Vast differences in diversity were observed among the 42 VNTR loci, ranging from 2 to 11 alleles. We found that the maximum copy number of repeats in an array was highly correlated with diversity (R = 0.86). VNTR-based phylogenetic analysis of the 24 strains successfully grouped isolates from biovar orientalis and most of the antiqua and mediaevalis strains. Hence, multiple-locus VNTR analysis (MLVA) appears capable of both distinguishing closely related strains and successfully classifying more distant relationships. Harnessing the power of MLVA to establish standardized databases will enable researchers to better understand plague ecology and evolution around the world.

Genetic mapping of soybean cyst nematode race-3 resistance loci in the soybean PI 437.654.

Resistance to the soybean cyst nematode (SCN) (Heterodera glycines Ichinohe) is difficult to evaluate in soybean [Glycine max (L.) Merr.] breeding. PI 437.654 has resistance to more SCN race isolates than any other known soybean. We screened 298 F6∶7 recombinant-inbred lines from a cross between PI 437.654 and ‘BSR101’ for SCN race-3 resistance, genetically mapped 355 RFLP markers and the I locus, and tested these markers for association with resistance loci. The Rhg4 resistance locus was within 1 cM of the I locus on linkage group A. Two additional QTLs associated with SCN resistance were located within 3cM of markers on groups G and M. These two loci were not independent because 91 of 96 lines that had a resistant-parent marker type on group G also had a resistant-parent marker type on group M. Rhg4 and the QTL on G showed a significant interaction by together providing complete resistance to SCN race-3. Individually, the QTL on G had greater effect on resistance than did Rhg4, but neither locus alone provided a degree of resistance much different from the susceptible parent. The nearest markers to the mapped QTLs on groups A and G had allele frequencies from the resistant parent indicating 52 resistant lines in this population, a number not significantly different from the 55 resistant lines found. Therefore, no QTLs from PI 437.654 other than those mapped here are expected to be required for resistance to SCN race-3. All 50 lines that had the PI 437.654 marker type at the nearest marker to each of the QTLs on groups A and G were resistant to SCN race-3. We believe markers near to these QTLs can be used effectively to select for SCN race-3 resistance, thereby improving the ability to breed SCN-resistant soybean varieties.

Necrotizing Cutaneous Mucormycosis after a Tornado in Joplin, Missouri, in 2011

BACKGROUND Mucormycosis is a fungal infection caused by environmentally acquired molds. We investigated a cluster of cases of cutaneous mucormycosis among persons injured during the May 22, 2011, tornado in Joplin, Missouri. METHODS We defined a case as a soft-tissue infection in a person injured during the tornado, with evidence of a mucormycete on culture or immunohistochemical testing plus DNA sequencing. We conducted a case-control study by reviewing medical records and conducting interviews with case patients and hospitalized controls. DNA sequencing and whole-genome sequencing were performed on clinical specimens to identify species and assess strain-level differences, respectively. RESULTS A total of 13 case patients were identified, 5 of whom (38%) died. The patients had a median of 5 wounds (range, 1 to 7); 11 patients (85%) had at least one fracture, 9 (69%) had blunt trauma, and 5 (38%) had penetrating trauma. All case patients had been located in the zone that sustained the most severe damage during the tornado. On multivariate analysis, infection was associated with penetrating trauma (adjusted odds ratio for case patients vs. controls, 8.8; 95% confidence interval [CI], 1.1 to 69.2) and an increased number of wounds (adjusted odds ratio, 2.0 for each additional wound; 95% CI, 1.2 to 3.2). Sequencing of the D1-D2 region of the 28S ribosomal DNA yielded Apophysomyces trapeziformis in all 13 case patients. Whole-genome sequencing showed that the apophysomyces isolates were four separate strains. CONCLUSIONS We report a cluster of cases of cutaneous mucormycosis among Joplin tornado survivors that were associated with substantial morbidity and mortality. Increased awareness of fungi as a cause of necrotizing soft-tissue infections after a natural disaster is warranted.

DNA methylation and AFLP marker distribution in the soybean genome

Abstract Amplified fragment length polymorphisms (AFLPs) have become important markers for genetic mapping because of their ability to reliably detect variation at a large number of loci. We report here the dissimilar distribution of two types of AFLP markers generated using restriction enzymes with varying sensitivities to cytosine methylation in the soybean genome. Initially, AFLP markers were placed on a scaffold map of 165 RFLP markers mapped in 42 recombinant inbred (F6:7) lines. These markers were selected from a map of over 500 RFLPs analyzed in 300 recombinant inbred (F6:7) lines generated by crossing BSR101×PI437.654. The randomness of AFLP marker map position was tested using a Poisson-model distribution. We found that AFLP markers generated using EcoRI/MseI deviated significantly from a random distribution, with 34% of the markers displaying dense clustering. In contrast to the EcoRI/MseI AFLP markers, PstI/MseI-generated AFLP markers did not cluster and were under represented in the EcoRI/MseI marker clusters. The restriction enzyme PstI is notably sensitive to cytosine methylation, and these results suggest that this sensitivity affected the distribution of the AFLP markers generated using this enzyme in the soybean genome. The common presence of one EcoRI/MseI AFLP cluster per linkage group and the infrequent presence of markers sensitive to methylation in these clusters are consistent with the low recombination frequency and the high level of cytosine methylation observed in the heterochromatic regions surrounding centromeres. Thus, the dense EcoRI/MseI AFLP marker clusters may be revealing structural features of the soybean genome, including the genetic locations of centromeres.

Molecular diversity in Bacillus anthracis

Molecular typing of Bacillus anthracis has been extremely difficult due to the lack of polymorphic DNA markers. We have identified nine novel variable number tandemly repeated loci from previously known amplified fragment length polymorphism markers or from the DNA sequence. In combination with the previously known vrrA locus, these markers provide discrimination power to genetically characterize B. anthracis isolates. The variable number tandem repeat (VNTR) loci are found in both gene coding (genic) and non‐coding (non‐genic) regions. The genic differences are ‘in frame’ and result in additions or deletion of amino acids to the predicted proteins. Due the rarity of molecular differences, the VNTR changes represent a significant portion of the genetic variation found within B. anthracis. This variation could represent an important adaptive mechanism. Marker similarity and differences among diverse isolates have identified seven major diversity groups that may represent the only world‐wide B. anthracis clones. The lineages reconstructed using these data may reflect the dispersal and evolution of this pathogen.

In Vitro Selection and Characterization of Bacillus anthracis Mutants with High-Level Resistance to Ciprofloxacin

ABSTRACT Mutants of attenuated Bacillus anthracis with high-level ciprofloxacin resistance were isolated using a three-step in vitro selection. Ciprofloxacin MICs were 0.5 μg/ml for first-step mutants, which had one of two gyrA quinolone resistance-determining region (QRDR) mutations. Ciprofloxacin MICs were 8 and 16 μg/ml for second-step mutants, which had one of three parC QRDR mutations. Ciprofloxacin MICs for third-step mutants were 32 and 64 μg/ml. Mutants for which MICs were 64 μg/ml had one of two additional mutations within the gyrA QRDR or one of two mutations within the gyrB QRDR. Mutants for which MICs were 32 μg/ml had no additional target modifications but showed evidence of enhanced ciprofloxacin efflux.

Cryptococcus gattii in North American Pacific Northwest: Whole-Population Genome Analysis Provides Insights into Species Evolution and Dispersal

ABSTRACT The emergence of distinct populations of Cryptococcus gattii in the temperate North American Pacific Northwest (PNW) was surprising, as this species was previously thought to be confined to tropical and semitropical regions. Beyond a new habitat niche, the dominant emergent population displayed increased virulence and caused primary pulmonary disease, as opposed to the predominantly neurologic disease seen previously elsewhere. Whole-genome sequencing was performed on 118 C. gattii isolates, including the PNW subtypes and the global diversity of molecular type VGII, to better ascertain the natural source and genomic adaptations leading to the emergence of infection in the PNW. Overall, the VGII population was highly diverse, demonstrating large numbers of mutational and recombinational events; however, the three dominant subtypes from the PNW were of low diversity and were completely clonal. Although strains of VGII were found on at least five continents, all genetic subpopulations were represented or were most closely related to strains from South America. The phylogenetic data are consistent with multiple dispersal events from South America to North America and elsewhere. Numerous gene content differences were identified between the emergent clones and other VGII lineages, including genes potentially related to habitat adaptation, virulence, and pathology. Evidence was also found for possible gene introgression from Cryptococcus neoformans var. grubii that is rarely seen in global C. gattii but that was present in all PNW populations. These findings provide greater understanding of C. gattii evolution in North America and support extensive evolution in, and dispersal from, South America. IMPORTANCE Cryptococcus gattii emerged in the temperate North American Pacific Northwest (PNW) in the late 1990s. Beyond a new environmental niche, these emergent populations displayed increased virulence and resulted in a different pattern of clinical disease. In particular, severe pulmonary infections predominated in contrast to presentation with neurologic disease as seen previously elsewhere. We employed population-level whole-genome sequencing and analysis to explore the genetic relationships and gene content of the PNW C. gattii populations. We provide evidence that the PNW strains originated from South America and identified numerous genes potentially related to habitat adaptation, virulence expression, and clinical presentation. Characterization of these genetic features may lead to improved diagnostics and therapies for such fungal infections. The data indicate that there were multiple recent introductions of C. gattii into the PNW. Public health vigilance is warranted for emergence in regions where C. gattii is not thought to be endemic. Cryptococcus gattii emerged in the temperate North American Pacific Northwest (PNW) in the late 1990s. Beyond a new environmental niche, these emergent populations displayed increased virulence and resulted in a different pattern of clinical disease. In particular, severe pulmonary infections predominated in contrast to presentation with neurologic disease as seen previously elsewhere. We employed population-level whole-genome sequencing and analysis to explore the genetic relationships and gene content of the PNW C. gattii populations. We provide evidence that the PNW strains originated from South America and identified numerous genes potentially related to habitat adaptation, virulence expression, and clinical presentation. Characterization of these genetic features may lead to improved diagnostics and therapies for such fungal infections. The data indicate that there were multiple recent introductions of C. gattii into the PNW. Public health vigilance is warranted for emergence in regions where C. gattii is not thought to be endemic.