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

Prevalence and characterization of pertactin deficient Bordetella pertussis strains in Brazil, a whole-cell vaccine country

Many countries have reported antigenic divergence among circulating Bordetella pertussis strains, mainly in those countries which introduced the acellular pertussis (aP) vaccine. This phenomenon can be seen, for example, with the recent rise of pertactin (Prn)-deficient B. pertussis strains, one of the antigens included in aP vaccine formulas. The whole cell pertussis (wP) vaccine has been used in Brazil since 1977 for the primary pertussis, diphtheria and tetanus immunization series. In 2014, the aP vaccine was recommended for women during pregnancy to protect infants in the first months of life. Our objective was to determine the prevalence of Prn-deficiency in 511 isolates of B. pertussis collected in Brazil during 2010-2016. All isolates were characterized, through PFGE and serotyping, and screened for the loss of Prn by ELISA. Prn-deficiency was confirmed by immunoblotting, and identification of the possible genetic markers was performed with PCR and Sanger sequencing. Results indicate that 110 PFGE profiles are currently circulating, with five profiles representing the majority, and the predominant serotype 3, has been gradually replaced by serotype 2 and serotype 2,3. ELISA screening and immunoblotting identified three Prn-deficient isolates. Genotypic characterization by PCR and sequencing indicated that one isolate had a promoter mutation in prn, while the other two did not have an obvious genetic explanation for their deficiency. While the lack of Prn was identified in a few isolates, this study did not detect a relevant occurrence of Prn-deficiency, until 2016, confirming previous observations that Prn-deficiency is likely aP vaccine-driven.

Genomic Surveillance and Improved Molecular Typing of Bordetella pertussis Using wgMLST

Multilocus sequence typing (MLST) provides allele-based characterization of bacterial pathogens in a standardized framework. However, classical MLST schemes for Bordetella pertussis, the causative agent of whooping cough, seldom reveal diversity among the small number of gene targets and thereby fail to delineate population structure. To improve the discriminatory power of allele-based molecular typing of B. pertussis, we have developed a whole-genome MLST (wgMLST) scheme from 225 reference-quality genome assemblies. Iterative refinement and allele curation resulted in a scheme of 3,506 coding sequences and covering 81.4% of the B. pertussis genome. This wgMLST scheme was further evaluated with data from a convenience sample of 2,389 B. pertussis isolates sequenced on Illumina instruments, including isolates from known outbreaks and epidemics previously characterized by existing molecular assays, as well as replicates collected from individual patients. wgMLST demonstrated concordance with whole-genome single nucleotide polymorphism (SNP) profiles, accurately resolved outbreak and sporadic cases in a retrospective comparison, and clustered replicate isolates collected from individual patients during diagnostic confirmation. Additionally, a reanalysis of isolates from two statewide epidemics using wgMLST reconstructed the population structures of circulating strains with increased resolution, revealing new clusters of related cases. Comparison with an existing core genome (cgMLST) scheme highlights the stable gene content of this bacterium and forms the initial foundation for necessary standardization. These results demonstrate the utility of wgMLST for improving B. pertussis characterization and genomic surveillance during the current pertussis disease resurgence.

Rapid and simple SNP genotyping for Bordetella pertussis epidemic strain MT27 based on a multiplexed single-base extension assay

Multilocus variable-number tandem repeat analysis (MLVA) is widely used for genotyping of Bordetella pertussis, the causative bacteria for pertussis. However, MLVA genotyping is losing its discriminate power because prevalence of the epidemic MT27 strain (MLVA-27) is increasing worldwide. To address this, we developed a single nucleotide polymorphism (SNP) genotyping method for MT27 based on multiplexed single-base extension (SBE) assay. A total of 237 MT27 isolates collected in Japan during 1999–2018 were genotyped and classified into ten SNP genotypes (SG1 to SG10) with a Simpson’s diversity index (DI) of 0.79 (95% CI 0.76–0.82). Temporal trends showed a marked increase in the genotypic diversity in the 2010s: Simpson’s DI was zero in 1999–2004, 0.16 in 2005–2009, 0.83 in 2010–2014, and 0.76 in 2015–2018. This indicates that the SNP genotyping is applicable to the recently circulating MT27 strain. Additionally, almost all outbreak-associated MT27 isolates were classified into the same SNP genotypes for each outbreak. Multiplexed SBE assay allows for rapid and simple genotyping, indicating that the SNP genotyping can potentially be a useful tool for subtyping the B. pertussis MT27 strain in routine surveillance and outbreak investigations.

Population dynamics and antigenic drift of Bordetella pertussis following whole cell vaccine replacement, Barcelona, Spain, 1986-2015

Among the factors associated with the resurgence of whooping cough, special emphasis has been given to pathogen adaptation after the introduction of the acellular vaccine (ACV). To assess the impact of the vaccine transition strategy from whole-cell vaccine (WCV) to ACV on population dynamics of Bordetella pertussis in Barcelona (Spain), we studied 339 isolates collected from 1986 to 2015 by PFGE and multi-locus variable-number tandem repeat analysis (MLVA). Additionally, allelic variants for the pertussis toxin and its promoter, pertactin, type 3 fimbriae and fimbrial serotyping were assessed to determine its antigenic drift. A shift was observed in the B. pertussis population as well as in its antigenic profile concurrently with the introduction of ACV in Barcelona. Four out of the five most prevalent PFGE profiles were replaced by new profiles following the ACV introduction. MLVA type 27 was the dominant genotype, and its frequency increased from 25% to 79.3% after WCV replacement. Antigen typing demonstrated the emergence of prn2, ptxP3, fim3-2 and a shift from the fimbriae 3 to the fimbriae 2 serotypes after the ACV introduction. Our findings support the presence of population and antigenic dynamic changes in B. pertussis likely driven by the introduction of ACV.

Genomic Sequencing of Bordetella pertussis for Epidemiology and Global Surveillance of Whooping Cough

Bordetella pertussis causes whooping cough, a highly contagious respiratory disease that is reemerging in many world regions. The spread of antigen-deficient strains may threaten acellular vaccine efficacy. Dynamics of strain transmission are poorly defined because of shortcomings in current strain genotyping methods. Our objective was to develop a whole-genome genotyping strategy with sufficient resolution for local epidemiologic questions and sufficient reproducibility to enable international comparisons of clinical isolates. We defined a core genome multilocus sequence typing scheme comprising 2,038 loci and demonstrated its congruence with whole-genome single-nucleotide polymorphism variation. Most cases of intrafamilial groups of isolates or of multiple isolates recovered from the same patient were distinguished from temporally and geographically cocirculating isolates. However, epidemiologically unrelated isolates were sometimes nearly undistinguishable. We set up a publicly accessible core genome multilocus sequence typing database to enable global comparisons of B. pertussis isolates, opening the way for internationally coordinated surveillance.

Molecular Epidemiology of Bordetella pertussis

Although vaccination has been effective, Bordetella pertussis is increasingly causing epidemics, especially in industrialized countries using acellular vaccines (aPs). One factor behind the increased circulation is the molecular changes on the pathogen level. After pertussis vaccinations were introduced, changes in the fimbrial (Fim) serotype of the circulating strains was observed. When bacterial typing methods improved, further changes between the vaccine and circulating strains, especially among the common virulence genes including pertussis toxin (PT) and pertactin (PRN) were noticed. Moreover, development of genome based techniques including pulsed-field gel electrophoresis (PFGE), multiple-locus variable number tandem repeat analysis (MLVA) and whole-genome sequencing (WGS) have offered a better resolution to monitor B. pertussis strains. After the introduction of aP vaccines, B. pertussis strains that are deficient to vaccine antigens, especially PRN, have appeared widely. On the other hand, antimicrobial resistance to first line drugs (macrolides) against B. pertussis is still low in many countries and therefore no globally evaluated antimicrobial susceptibility test values have been recommended. In this review, we focus on the molecular changes in the bacteria, which have or may have affected the past and current epidemiology of pertussis.

Pertactin-deficient Bordetella pertussis isolates in Poland-a country with whole-cell pertussis primary vaccination

The introduction of pertussis vaccination in the 1950s resulted in a significant decrease in the incidence of disease. However, since the 1990s many highly vaccinated countries have observed the re-emergence of the disease. One of the causes of this phenomenon might be related to the adaptation of Bordetella pertussis to vaccination. The purpose of the presented study was an investigation of the emergence and spread of vaccine antigen-deficient B. pertussis isolates in Poland and genomic characterization of the currently circulating pathogen population using PFGE, MLVA and MAST. The results revealed that all tested isolates expressed Ptx, FHA and ACT antigens but 15.4% (4/26) of isolates from 2010 to 2016 were Prn-deficient. Moreover, one TcfA-deficient isolate was collected in 2015. The genotyping showed a genetic distinction between the isolates circulating in 2010-2016 and isolates from previous periods. The majority of currently circulating isolates belonged to PFGE group IV (96.2%), type MT27 (73.1%), and carried ptxA1-ptxC2-ptxP3-prn2-tcfA2-fim2-1-fim3-1 alleles (61.5%). The unique genetic structure of the B. pertussis population in Poland has changed since 2010 and became similar to that observed in countries with aP vaccination. This could be a result of increasing use of aP vaccines (60% of primary vaccination in 2013) over wP vaccines, which have been broadly used for primary vaccination in Poland for decades.

Screening and Genomic Characterization of Filamentous Hemagglutinin-Deficient Bordetella pertussis

Despite high vaccine coverage, pertussis cases in the United States have increased over the last decade. Growing evidence suggests that disease resurgence results, in part, from genetic divergence of circulating strain populations away from vaccine references. The United States employs acellular vaccines exclusively, and current Bordetella pertussis isolates are predominantly deficient in at least one immunogen, pertactin (Prn). First detected in the United States retrospectively in a 1994 isolate, the rapid spread of Prn deficiency is likely vaccine driven, raising concerns about whether other acellular vaccine immunogens experience similar pressures, as further antigenic changes could potentially threaten vaccine efficacy. We developed an electrochemiluminescent antibody capture assay to monitor the production of the acellular vaccine immunogen filamentous hemagglutinin (Fha). Screening 722 U.S. surveillance isolates collected from 2010 to 2016 identified two that were both Prn and Fha deficient. Three additional Fha-deficient laboratory strains were also identified from a historic collection of 65 isolates dating back to 1935. Whole-genome sequencing of deficient isolates revealed putative, underlying genetic changes. Only four isolates harbored mutations to known genes involved in Fha production, highlighting the complexity of its regulation. The chromosomes of two Fha-deficient isolates included unexpected structural variation that did not appear to influence Fha production. Furthermore, insertion sequence disruption of fhaB was also detected in a previously identified pertussis toxin-deficient isolate that still produced normal levels of Fha. These results demonstrate the genetic potential for additional vaccine immunogen deficiency and underscore the importance of continued surveillance of circulating B. pertussis evolution in response to vaccine pressure.

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 pathogen's 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.

The characterization of Bordetella pertussis strains isolated in the Central-Western region of Brazil suggests the selection of a specific genetic profile during 2012-2014 outbreaks

Pertussis is a worldwide acute respiratory disease caused by the bacterium Bordetella pertussis. Despite high vaccine coverage, the bacterium continues to circulate in populations and is still one of the most common vaccine-preventable diseases. In Brazil, pertussis incidence has presented a significant decrease since 1990 but since 2011 a sudden increase in incidence has been observed. Thus, the aim of this study was to perform a molecular epidemiological characterization of B. pertussis strains isolated in the Central-Western region (specifically in Distrito Federal) of Brazil from August 2012 to August 2014. During this period, 92 B. pertussis strains were isolated from the outbreaks. All strains were characterized by serotyping and XbaI pulsed-field gel electrophoresis profiles. From August to December 2012, the most prevalent serotype observed was 1,3 (13/17). During 2013 the prevalence of serotype 1,3 decreased (13/30) and from January 2014 to August 2014 the most prevalent serotype was 1,2 (33/45). Fourteen PFGE profiles were identified. Of these, BP-XbaI0039 prevalence increased from 3/17 in 2012 to 10/30 in 2013, and 35/45 in 2014. These results evidence the selection of a specific genetic profile during this period, suggesting the occurrence of a bacterial genomic profile with high circulation potential.