M. Lucia Tondella

Global Population Structure and Evolution of Bordetella pertussis and Their Relationship with Vaccination

ABSTRACT Bordetella pertussis causes pertussis, a respiratory disease that is most severe for infants. Vaccination was introduced in the 1950s, and in recent years, a resurgence of disease was observed worldwide, with significant mortality in infants. Possible causes for this include the switch from whole-cell vaccines (WCVs) to less effective acellular vaccines (ACVs), waning immunity, and pathogen adaptation. Pathogen adaptation is suggested by antigenic divergence between vaccine strains and circulating strains and by the emergence of strains with increased pertussis toxin production. We applied comparative genomics to a worldwide collection of 343 B. pertussis strains isolated between 1920 and 2010. The global phylogeny showed two deep branches; the largest of these contained 98% of all strains, and its expansion correlated temporally with the first descriptions of pertussis outbreaks in Europe in the 16th century. We found little evidence of recent geographical clustering of the strains within this lineage, suggesting rapid strain flow between countries. We observed that changes in genes encoding proteins implicated in protective immunity that are included in ACVs occurred after the introduction of WCVs but before the switch to ACVs. Furthermore, our analyses consistently suggested that virulence-associated genes and genes coding for surface-exposed proteins were involved in adaptation. However, many of the putative adaptive loci identified have a physiological role, and further studies of these loci may reveal less obvious ways in which B. pertussis and the host interact. This work provides insight into ways in which pathogens may adapt to vaccination and suggests ways to improve pertussis vaccines. IMPORTANCE Whooping cough is mainly caused by Bordetella pertussis, and current vaccines are targeted against this organism. Recently, there have been increasing outbreaks of whooping cough, even where vaccine coverage is high. Analysis of the genomes of 343 B. pertussis isolates from around the world over the last 100 years suggests that the organism has emerged within the last 500 years, consistent with historical records. We show that global transmission of new strains is very rapid and that the worldwide population of B. pertussis is evolving in response to vaccine introduction, potentially enabling vaccine escape. Whooping cough is mainly caused by Bordetella pertussis, and current vaccines are targeted against this organism. Recently, there have been increasing outbreaks of whooping cough, even where vaccine coverage is high. Analysis of the genomes of 343 B. pertussis isolates from around the world over the last 100 years suggests that the organism has emerged within the last 500 years, consistent with historical records. We show that global transmission of new strains is very rapid and that the worldwide population of B. pertussis is evolving in response to vaccine introduction, potentially enabling vaccine escape.

Population diversity among Bordetella pertussis isolates, United States, 1935-2009.

Resurgence of pertussis was not directly correlated with changes in vaccine composition or schedule.

Bordetella pertussis Strain Lacking Pertactin and Pertussis Toxin.

A Bordetella pertussis strain lacking 2 acellular vaccine immunogens, pertussis toxin and pertactin, was isolated from an unvaccinated infant in New York State in 2013. Comparison with a French strain that was pertussis toxin–deficient, pertactin wild-type showed that the strains carry the same 28-kb deletion in similar genomes.

Harmonization of Bordetella pertussis Real-Time PCR Diagnostics in the United States in 2012

ABSTRACT Real-time PCR (rt-PCR) is an important diagnostic tool for the identification of Bordetella pertussis, Bordetella holmesii, and Bordetella parapertussis. Most U.S. public health laboratories (USPHLs) target IS481, present in 218 to 238 copies in the B. pertussis genome and 32 to 65 copies in B. holmesii. The CDC developed a multitarget PCR assay to differentiate B. pertussis, B. holmesii, and B. parapertussis and provided protocols and training to 19 USPHLs. The 2012 performance exercise (PE) assessed the capability of USPHLs to detect these three Bordetella species in clinical samples. Laboratories were recruited by the Wisconsin State Proficiency Testing program through the Association of Public Health Laboratories, in partnership with the CDC. Spring and fall PE panels contained 12 samples each of viable Bordetella and non-Bordetella species in saline. Fifty and 53 USPHLs participated in the spring and fall PEs, respectively, using a variety of nucleic acid extraction methods, PCR platforms, and assays. Ninety-six percent and 94% of laboratories targeted IS481 in spring and fall, respectively, in either singleplex or multiplex assays. In spring and fall, respectively, 72% and 79% of USPHLs differentiated B. pertussis and B. holmesii and 68% and 72% identified B. parapertussis. IS481 cycle threshold (CT ) values for B. pertussis samples had coefficients of variation (CV) ranging from 10% to 28%. Of the USPHLs that differentiated B. pertussis and B. holmesii, sensitivity was 96% and specificity was 95% for the combined panels. The 2012 PE demonstrated increased harmonization of rt-PCR Bordetella diagnostic protocols in USPHLs compared to that of the previous survey.

Genome Structural Diversity among 31 Bordetella pertussis Isolates from Two Recent U.S. Whooping Cough Statewide Epidemics.

Pertussis, or whooping cough, is the most poorly controlled vaccine-preventable bacterial disease in the United States, which has experienced a resurgence for more than a decade. Once viewed as a monomorphic pathogen, B. pertussis strains circulating during epidemics exhibit diversity visible on a genome structural level, previously undetectable by traditional sequence analysis using short-read technologies. For the first time, we combine short- and long-read sequencing platforms with restriction optical mapping for single-contig, de novo assembly of 31 isolates to investigate two geographically and temporally independent U.S. pertussis epidemics. These complete genomes reshape our understanding of B. pertussis evolution and strengthen molecular epidemiology toward one day understanding the resurgence of pertussis. ABSTRACT During 2010 and 2012, California and Vermont, respectively, experienced statewide epidemics of pertussis with differences seen in the demographic affected, case clinical presentation, and molecular epidemiology of the circulating strains. To overcome limitations of the current molecular typing methods for pertussis, we utilized whole-genome sequencing to gain a broader understanding of how current circulating strains are causing large epidemics. Through the use of combined next-generation sequencing technologies, this study compared de novo, single-contig genome assemblies from 31 out of 33 Bordetella pertussis isolates collected during two separate pertussis statewide epidemics and 2 resequenced vaccine strains. Final genome architecture assemblies were verified with whole-genome optical mapping. Sixteen distinct genome rearrangement profiles were observed in epidemic isolate genomes, all of which were distinct from the genome structures of the two resequenced vaccine strains. These rearrangements appear to be mediated by repetitive sequence elements, such as high-copy-number mobile genetic elements and rRNA operons. Additionally, novel and previously identified single nucleotide polymorphisms were detected in 10 virulence-related genes in the epidemic isolates. Whole-genome variation analysis identified state-specific variants, and coding regions bearing nonsynonymous mutations were classified into functional annotated orthologous groups. Comprehensive studies on whole genomes are needed to understand the resurgence of pertussis and develop novel tools to better characterize the molecular epidemiology of evolving B. pertussis populations. IMPORTANCE Pertussis, or whooping cough, is the most poorly controlled vaccine-preventable bacterial disease in the United States, which has experienced a resurgence for more than a decade. Once viewed as a monomorphic pathogen, B. pertussis strains circulating during epidemics exhibit diversity visible on a genome structural level, previously undetectable by traditional sequence analysis using short-read technologies. For the first time, we combine short- and long-read sequencing platforms with restriction optical mapping for single-contig, de novo assembly of 31 isolates to investigate two geographically and temporally independent U.S. pertussis epidemics. These complete genomes reshape our understanding of B. pertussis evolution and strengthen molecular epidemiology toward one day understanding the resurgence of pertussis.

The History of Bordetella pertussis Genome Evolution Includes Structural Rearrangement.

Despite high pertussis vaccine coverage, reported cases of whooping cough (pertussis) have increased over the last decade in the United States and other developed countries. Although Bordetella pertussis is well known for its limited gene sequence variation, recent advances in long-read sequencing technology have begun to reveal genomic structural heterogeneity among otherwise indistinguishable isolates, even within geographically or temporally defined epidemics. We have compared rearrangements among complete genome assemblies from 257 B. pertussis isolates to examine the potential evolution of the chromosomal structure in a pathogen with minimal gene nucleotide sequence diversity. Discrete changes in gene order were identified that differentiated genomes from vaccine reference strains and clinical isolates of various genotypes, frequently along phylogenetic boundaries defined by single nucleotide polymorphisms. The observed rearrangements were primarily large inversions centered on the replication origin or terminus and flanked by IS481, a mobile genetic element with >240 copies per genome and previously suspected to mediate rearrangements and deletions by homologous recombination. These data illustrate that structural genome evolution in B. pertussis is not limited to reduction but also includes rearrangement. Therefore, although genomes of clinical isolates are structurally diverse, specific changes in gene order are conserved, perhaps due to positive selection, providing novel information for investigating disease resurgence and molecular epidemiology.IMPORTANCE Whooping cough, primarily caused by Bordetella pertussis, has resurged in the United States even though the coverage with pertussis-containing vaccines remains high. The rise in reported cases has included increased disease rates among all vaccinated age groups, provoking questions about the pathogens evolution. The chromosome of B. pertussis includes a large number of repetitive mobile genetic elements that obstruct genome analysis. However, these mobile elements facilitate large rearrangements that alter the order and orientation of essential protein-encoding genes, which otherwise exhibit little nucleotide sequence diversity. By comparing the complete genome assemblies from 257 isolates, we show that specific rearrangements have been conserved throughout recent evolutionary history, perhaps by eliciting changes in gene expression, which may also provide useful information for molecular epidemiology.

Changes in Predominance of Pulsed-Field Gel Electrophoresis Profiles of Bordetella pertussis Isolates, United States, 2000-2012.

These changes are concurrrent with other recent molecular changes and may be contributing to US pertussis reemergence.

Development of a qualitative assay for screening of Bordetella pertussis isolates for pertussis toxin production

Bordetella pertussis infection has been increasing in the US, with reported cases reaching over 50,000 in 2012, a number last observed in the 1950s. Concurrently, B. pertussis lacking the pertactin protein, one of the immunogens included in the acellular vaccine formulations, has rapidly emerged since 2010, and has become the predominant circulating phenotype. Monitoring the production of the remaining acellular vaccine immunogens, such as pertussis toxin (Pt), is a critical next step. To date, methods for screening Pt have been either through genomic sequencing means or by conventional ELISAs. However, sequencing limits detection to the DNA level, missing potential disruptions in transcription or translation. Conventional ELISAs are beneficial for detecting the protein; however, they can often suffer from poor sensitivity and specificity. Here we describe a rapid, highly sensitive and specific electrochemiluminescent capture ELISA that can detect Pt production in prepared inactivated bacterial suspensions. Over 340 isolates were analyzed and analytical validation parameters, such as precision, reproducibility, and stability, were rigorously tested. Intra-plate and inter-plate variability measured at 9.8% and 11.5%, respectively. Refrigerated samples remained stable for two months and variability was unaffected (coefficient of variation was 12%). Interestingly, despite the intention of being a qualitative method, the assay was sensitive enough to detect a small, but statistically significant, difference in protein production between different pertussis promoter allelic groups of strains, ptxP1 and ptxP3. This technology has the ability to perform screening of multiple antigens at one time, thus, improving testing characteristics while minimizing costs, specimen volume, and testing time.

Evaluation of Commercial Assays for Single-Point Diagnosis of Pertussis in the US

Background Pertussis serodiagnosis is increasingly being used in the United States despite the lack of a US Food and Drug Administration-approved, commercially available assay. To better understand the utility of these assays in diagnosing pertussis, serology assays were evaluated for analytical parameters and clinical accuracy. Methods Forty-three antigen-antibody combinations were evaluated for single-point diagnosis of pertussis. Serum panels included sera from laboratory-confirmed cases, an international reference standard, and healthy donors. Phase I panel (n = 20) of sera was used to assess precision, linearity, and accuracy; Phase II panel (n = 226) followed with positive percent agreement (PPA) and negative percent agreement (NPA) estimates. Analytical analyses included coefficients of variation (CV) and concordance correlation coefficients (rc). Results Intra-analyst variability was found to be relatively low among samples per assay, with only 6% (78 of 1240) having CV >20%, primarily with the highly concentrated immunoglobulin (Ig)G anti-pertussis toxin (PT) specimens and IgM assays. The rc measurements to assess linearity ranged between 0.282 and 0.994, 0.332 and 0.999, and -0.056 and 0.482 for IgA, IgG, and IgM, respectively. Analytical accuracy for calibrated IgG anti-PT assays was 86%-115%. The PPA and NPA varied greatly for all assays; PPA/NPA ranges for IgA, IgG, and IgM assays, with culture and/or polymerase chain reaction positivity as control, were 29-90/13-100, 26-96/27-100, and 0-73/42-100, respectively. In IgG assays, mixing filamentous hemagglutinin antigen with PT increased PPA but decreased NPA. Conclusions Seroassays varied substantially under both analytical and clinical parameters; however, those that were calibrated to a reference standard were highly accurate. Our findings support incorporation of calibrated pertussis seroassays to the pertussis case definition for improved diagnosis and surveillance.

Complete Genome Sequences of Bordetella pertussis Vaccine Reference Strains 134 and 10536

ABSTRACT Vaccine formulations and vaccination programs against whooping cough (pertussis) vary worldwide. Here, we report the complete genome sequences of two divergent Bordetella pertussis reference strains used in the production of pertussis vaccines.