Jeannine M. Petersen

Yersinia pestis genome sequencing identifies patterns of global phylogenetic diversity

Pandemic infectious diseases have accompanied humans since their origins1, and have shaped the form of civilizations2. Of these, plague is possibly historically the most dramatic. We reconstructed historical patterns of plague transmission through sequence variation in 17 complete genome sequences and 933 single nucleotide polymorphisms (SNPs) within a global collection of 286 Yersinia pestis isolates. Y. pestis evolved in or near China, and has been transmitted via multiple epidemics that followed various routes, probably including transmissions to West Asia via the Silk Road and to Africa by Chinese marine voyages. In 1894, Y. pestis spread to India and radiated to diverse parts of the globe, leading to country-specific lineages that can be traced by lineage-specific SNPs. All 626 current isolates from the U.S.A. reflect one radiation and 82 isolates from Madagascar represent a second. Subsequent local microevolution of Y. pestis is marked by sequential, geographically-specific SNPs.Plague is a pandemic human invasive disease caused by the bacterial agent Yersinia pestis. We here report a comparison of 17 whole genomes of Y. pestis isolates from global sources. We also screened a global collection of 286 Y. pestis isolates for 933 SNPs using Sequenom MassArray SNP typing. We conducted phylogenetic analyses on this sequence variation dataset, assigned isolates to populations based on maximum parsimony and, from these results, made inferences regarding historical transmission routes. Our phylogenetic analysis suggests that Y. pestis evolved in or near China and spread through multiple radiations to Europe, South America, Africa and Southeast Asia, leading to country-specific lineages that can be traced by lineage-specific SNPs. All 626 current isolates from the United States reflect one radiation, and 82 isolates from Madagascar represent a second radiation. Subsequent local microevolution of Y. pestis is marked by sequential, geographically specific SNPs.

Supplementary information to support : 'Yersinia pestis genome sequencing identifies patterns of global phylogenetic diversity'

Pandemic infectious diseases have accompanied humans since their origins1, and have shaped the form of civilizations2. Of these, plague is possibly historically the most dramatic. We reconstructed historical patterns of plague transmission through sequence variation in 17 complete genome sequences and 933 single nucleotide polymorphisms (SNPs) within a global collection of 286 Yersinia pestis isolates. Y. pestis evolved in or near China, and has been transmitted via multiple epidemics that followed various routes, probably including transmissions to West Asia via the Silk Road and to Africa by Chinese marine voyages. In 1894, Y. pestis spread to India and radiated to diverse parts of the globe, leading to country-specific lineages that can be traced by lineage-specific SNPs. All 626 current isolates from the U.S.A. reflect one radiation and 82 isolates from Madagascar represent a second. Subsequent local microevolution of Y. pestis is marked by sequential, geographically-specific SNPs.Plague is a pandemic human invasive disease caused by the bacterial agent Yersinia pestis. We here report a comparison of 17 whole genomes of Y. pestis isolates from global sources. We also screened a global collection of 286 Y. pestis isolates for 933 SNPs using Sequenom MassArray SNP typing. We conducted phylogenetic analyses on this sequence variation dataset, assigned isolates to populations based on maximum parsimony and, from these results, made inferences regarding historical transmission routes. Our phylogenetic analysis suggests that Y. pestis evolved in or near China and spread through multiple radiations to Europe, South America, Africa and Southeast Asia, leading to country-specific lineages that can be traced by lineage-specific SNPs. All 626 current isolates from the United States reflect one radiation, and 82 isolates from Madagascar represent a second radiation. Subsequent local microevolution of Y. pestis is marked by sequential, geographically specific SNPs.

Development of a Multitarget Real-Time TaqMan PCR Assay for Enhanced Detection of Francisella tularensis in Complex Specimens

ABSTRACT Tularemia is the zoonotic disease caused by the gram-negative coccobacillus Francisella tularensis. Its wide distribution in the environment poses a challenge for understanding the transmission, ecology, and epidemiology of the disease. F. tularensis is also considered a potential biological weapon due to its extreme infectivity. We have developed a multitarget real-time TaqMan PCR assay capable of rapidly and accurately detecting F. tularensis in complex specimens. Targeted regions included the ISFtu2 element and the 23kDa, fopA, and tul4 genes. Analysis of the four TaqMan assays demonstrated that three (ISFtu2, 23kDa, and tul4) performed within our established criterion of a detection limit of one organism. The combined use of the three assays was highly specific, displaying no cross-reactivity with the non-Francisella bacteria tested and capable of differentially diagnosing both F. tularensis and Francisella philomiragia. When the multitarget TaqMan assay (ISFtu2, 23kDa, and tul4) was compared to culturing, using environmentally contaminated specimens, the TaqMan PCR assay was significantly more sensitive than culturing (P ≤ 0.05). The sensitive and specific nature of this rapid multitarget TaqMan assay provides a valuable new tool that with future evaluations can be used for analyzing clinical specimens, field samples during bioterrorism threat assessment, and samples from outbreaks and for improving our understanding of the ecology and environmental prevalence of F. tularensis.

Epidemiologic and molecular analysis of human tularemia, United States, 1964-2004.

Distinct subpopulations of F. tularensis differ in their clinical manifestations, geographic distribution, and likely modes of transmission.

Francisella tularensis in the United States.

Subpopulations A.I and A.II. of Francisella tularensis subsp. tularensis are associated with unique biotic and abiotic factors that maintain disease foci.

Francisella tularensis: an arthropod-borne pathogen.

Arthropod transmission of tularemia occurs throughout the northern hemisphere. Few pathogens show the adaptability of Francisella tularensis to such a wide array of arthropod vectors. Nonetheless, arthropod transmission of F. tularensis was last actively investigated in the first half of the 20th century. This review will focus on arthropod transmission to humans with respect to vector species, modes of transmission, geographic differences and F. tularensis subspecies and clades.

Phylogeography of Francisella tularensis: global expansion of a highly fit clone.

Francisella tularensis contains several highly pathogenic subspecies, including Francisella tularensis subsp. holarctica, whose distribution is circumpolar in the northern hemisphere. The phylogeography of these subspecies and their subclades was examined using whole-genome single nucleotide polymorphism (SNP) analysis, high-density microarray SNP genotyping, and real-time-PCR-based canonical SNP (canSNP) assays. Almost 30,000 SNPs were identified among 13 whole genomes for phylogenetic analysis. We selected 1,655 SNPs to genotype 95 isolates on a high-density microarray platform. Finally, 23 clade- and subclade-specific canSNPs were identified and used to genotype 496 isolates to establish global geographic genetic patterns. We confirm previous findings concerning the four subspecies and two Francisella tularensis subsp. tularensis subpopulations and identify additional structure within these groups. We identify 11 subclades within F. tularensis subsp. holarctica, including a new, genetically distinct subclade that appears intermediate between Japanese F. tularensis subsp. holarctica isolates and the common F. tularensis subsp. holarctica isolates associated with the radiation event (the B radiation) wherein this subspecies spread throughout the northern hemisphere. Phylogenetic analyses suggest a North American origin for this B-radiation clade and multiple dispersal events between North America and Eurasia. These findings indicate a complex transmission history for F. tularensis subsp. holarctica.

Identification of a novel pathogenic Borrelia species causing Lyme borreliosis with unusually high spirochaetaemia: a descriptive study

BACKGROUND Lyme borreliosis is the most common tick-borne disease in the northern hemisphere. It is a multisystem disease caused by Borrelia burgdorferi sensu lato genospecies and characterised by tissue localisation and low spirochaetaemia. In this study we aimed to describe a novel Borrelia species causing Lyme borreliosis in the USA. METHODS At the Mayo clinic, from 2003 to 2014, we tested routine clinical diagnostic specimens from patients in the USA with PCR targeting the oppA1 gene of B burgdorferi sensu lato. We identified positive specimens with an atypical PCR result (melting temperature outside of the expected range) by sequencing, microscopy, or culture. We collected Ixodes scapularis ticks from regions of suspected patient tick exposure and tested them by oppA1 PCR. FINDINGS 100 545 specimens were submitted by physicians for routine PCR from Jan 1, 2003 to Sept 30, 2014. From these samples, six clinical specimens (five blood, one synovial fluid) yielded an atypical oppA1 PCR product, but no atypical results were detected before 2012. Five of the six patients with atypical PCR results had presented with fever, four had diffuse or focal rash, three had symptoms suggestive of neurological inclusion, and two were admitted to hospital. The sixth patient presented with knee pain and swelling. Motile spirochaetes were seen in blood samples from one patient and cultured from blood samples from two patients. Among the five blood specimens, the median oppA1 copy number was 180 times higher than that in 13 specimens that tested positive for B burgdorferi sensu stricto during the same time period. Multigene sequencing identified the spirochaete as a novel B burgdorferi sensu lato genospecies. This same genospecies was detected in ticks collected at a probable patient exposure site. INTERPRETATION We describe a new pathogenic Borrelia burgdorferi sensu lato genospecies (candidatus Borrelia mayonii) in the upper midwestern USA, which causes Lyme borreliosis with unusually high spirochaetaemia. Clinicians should be aware of this new B burgdorferi sensu lato genospecies, its distinct clinical features, and the usefulness of oppA1 PCR for diagnosis. FUNDING US Centers for Disease Control and Prevention Epidemiology and Laboratory Capacity for Infectious Diseases (ELC) Cooperative Agreement and Mayo Clinic Small Grant programme.

Molecular Epidemiology of Francisella tularensis in the United States

BACKGROUND In the United States, tularemia is caused by Francisella tularensis subsps. tularensis (type A) and holarctica (type B). Molecular subtyping has further divided type A into 2 subpopulations, A1 and A2. Significant mortality differences were previously identified between human infections caused by A1 (14%), A2 (0%) and type B (7%). To verify these findings and to further define differences among genotypes, we performed a large-scale molecular epidemiologic analysis of F. tularensis isolates from humans and animals. METHODS Pulsed-field gel electrophoresis with PmeI was performed on 302 type A and 61 type B isolates. Pulsed-field gel electrophoresis pattern and epidemiologic analyses were performed. Logistic regression was used to assess factors associated with human mortality. RESULTS Pulsed-field gel electrophoresis typing identified 4 distinct type A genotypes, A1a, A1b, A2a, and A2b, as well as type B. Genotypic and geographic divisions observed among isolates from humans were mirrored among isolates from animals, specifically among animal species that are linked to human infection and to enzootic maintenance of tularemia. Significant differences between human infections caused by different genotypes were identified with respect to patient age, site of organism recovery, and mortality. Human infections due to A1b resulted in significantly higher mortality (24%) than those caused by A1a (4%), A2 (0%), and type B (7%). CONCLUSIONS Three type A genotypes, A1a, A1b, and A2, were shown to be epidemiologically important. Our analysis suggests that A1b strains may be significantly more virulent in humans than A1a, A2, or type B strains. These findings have important implications for disease progression, disease prevention, and basic research programs.

From protein microarrays to diagnostic antigen discovery

MOTIVATION An important application of protein microarray data analysis is identifying a serodiagnostic antigen set that can reliably detect patterns and classify antigen expression profiles. This work addresses this problem using antibody responses to protein markers measured by a novel high-throughput microarray technology. The findings from this study have direct relevance to rapid, broad-based diagnostic and vaccine development. RESULTS Protein microarray chips are probed with sera from individuals infected with the bacteria Francisella tularensis, a category A biodefense pathogen. A two-step approach to the diagnostic process is presented (1) feature (antigen) selection and (2) classification using antigen response measurements obtained from F.tularensis microarrays (244 antigens, 46 infected and 54 healthy human sera measurements). To select antigens, a ranking scheme based on the identification of significant immune responses and differential expression analysis is described. Classification methods including k-nearest neighbors, support vector machines (SVM) and k-Means clustering are applied to training data using selected antigen sets of various sizes. SVM based models yield prediction accuracy rates in the range of approximately 90% on validation data, when antigen set sizes are between 25 and 50. These results strongly indicate that the top-ranked antigens can be considered high-priority candidates for diagnostic development. AVAILABILITY All software programs are written in R and available at http://www.igb.uci.edu/index.php?page=tools and at http://www.r-project.org. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.