Insight into evolution of Bordetella pertussis from comparative genomic analysis: evidence of vaccine-driven selection

Genome-based prediction of cross-protective, HLA-DR-presented epitopes as putative vaccine antigens for multiple Bordetella species

IMPORTANCE

Pertussis, caused by Bordetella pertussis, can cause debilitating respiratory symptoms, so whole-cell pertussis vaccines (wPVs) were introduced in the 1940s. However, reactogenicity of wPV necessitated the development of acellular pertussis vaccines (aPVs) that were introduced in the 1990s. Since then, until the COVID-19 pandemic began, reported pertussis incidence was increasing, suggesting that aPVs do not induce long-lasting immunity and may not effectively prevent transmission. Additionally, aPVs do not provide protection against other Bordetella species that are observed during outbreaks. The significance of this work is in determining potential new vaccine antigens for multiple Bordetella species that are predicted to elicit long-term immune responses. Genome-based approaches have aided the development of novel vaccines; here, these methods identified Bordetella vaccine candidates that may be cross-protective and predicted to induce strong memory responses. These targets can lead to an improved vaccine with a strong safety profile while also strengthening the longevity of the immune response.

Bordetella bronchiseptica and Bordetella pertussis: Similarities and Differences in Infection, Immuno-Modulation, and Vaccine Considerations

Bordetella pertussis and Bordetella bronchiseptica belong to the genus Bordetella, which comprises 14 other species. B. pertussis is responsible for whooping cough in humans, a severe infection in children and less severe or chronic in adults. These infections are restricted to humans and currently increasing worldwide. B. bronchiseptica is involved in diverse respiratory infections in a wide range of mammals. For instance, the canine infectious respiratory disease complex (CIRDC), characterized by a chronic cough in dogs. At the same time, it is increasingly implicated in human infections, while remaining an important pathogen in the veterinary field. Both Bordetella can evade and modulate host immune responses to support their persistence, although it is more pronounced in B. bronchiseptica infection. The protective immune responses elicited by both pathogens are comparable, while there are important characteristics in the mechanisms that differ. However, B. pertussis pathogenesis is more difficult to decipher in animal models than those of B. bronchiseptica because of its restriction to humans. Nevertheless, the licensed vaccines for each Bordetella are different in terms of formulation, route of administration and immune responses induced, with no known cross-reaction between them. Moreover, the target of the mucosal tissues and the induction of long-lasting cellular and humoral responses are required to control and eliminate Bordetella. In addition, the interaction between both veterinary and human fields are essential for the control of this genus, by preventing the infections in animals and the subsequent zoonotic transmission to humans.

Integrating proteomic data with metabolic modeling provides insight into key pathways of Bordetella pertussis biofilms

Pertussis, commonly known as whooping cough is a severe respiratory disease caused by the bacterium, Bordetella pertussis. Despite widespread vaccination, pertussis resurgence has been observed globally. The development of the current acellular vaccine (ACV) has been based on planktonic studies. However, recent studies have shown that B. pertussis readily forms biofilms. A better understanding of B. pertussis biofilms is important for developing novel vaccines that can target all aspects of B. pertussis infection. This study compared the proteomic expression of biofilm and planktonic B. pertussis cells to identify key changes between the conditions. Major differences were identified in virulence factors including an upregulation of toxins (adenylate cyclase toxin and dermonecrotic toxin) and downregulation of pertactin and type III secretion system proteins in biofilm cells. To further dissect metabolic pathways that are altered during the biofilm lifestyle, the proteomic data was then incorporated into a genome scale metabolic model using the Integrative Metabolic Analysis Tool (iMAT). The generated models predicted that planktonic cells utilised the glyoxylate shunt while biofilm cells completed the full tricarboxylic acid cycle. Differences in processing aspartate, arginine and alanine were identified as well as unique export of valine out of biofilm cells which may have a role in inter-bacterial communication and regulation. Finally, increased polyhydroxybutyrate accumulation and superoxide dismutase activity in biofilm cells may contribute to increased persistence during infection. Taken together, this study modeled major proteomic and metabolic changes that occur in biofilm cells which helps lay the groundwork for further understanding B. pertussis pathogenesis.

One species, different diseases: the unique molecular mechanisms that underlie the pathogenesis of typhoidal Salmonella infections

Salmonella enterica is one of the most widespread bacterial pathogens found worldwide, resulting in approximately 100 million infections and over 200 000 deaths per year. Salmonella isolates, termed 'serovars', can largely be classified as either nontyphoidal or typhoidal Salmonella, which differ in regard to disease manifestation and host tropism. Nontyphoidal Salmonella causes gastroenteritis in many hosts, while typhoidal Salmonella is human-restricted and causes typhoid fever, a systemic disease with a mortality rate of up to 30% without treatment. There has been considerable interest in understanding how different Salmonella serovars cause different diseases, but the molecular details that underlie these infections have not yet been fully characterized, especially in the case of typhoidal Salmonella. In this review, we highlight the current state of research into understanding the pathogenesis of both nontyphoidal and typhoidal Salmonella, with a specific interest in serovar-specific traits that allow human-adapted strains of Salmonella to cause enteric fever. Overall, a more detailed molecular understanding of how different Salmonella isolates infect humans will provide critical insights into how we can eradicate these dangerous enteric pathogens.

Genomic dissection of the microevolution of Australian epidemic Bordetella pertussis

Whooping cough (pertussis) is a highly contagious respiratory disease caused by the bacterium Bordetella pertussis. Despite high vaccine coverage, pertussis has re-emerged in many countries including Australia and caused two large epidemics in Australia since 2007. Here, we undertook a genomic and phylogeographic study of 385 Australian B. pertussis isolates collected from 2008 to 2017. The Australian B. pertussis population was found to be composed of mostly ptxP3 strains carrying different fim3 alleles, with ptxP3-fim3A genotype expanding far more than ptxP3-fim3B. Within the former, there were six co-circulating epidemic lineages (EL1 to EL6). The multiple ELs emerged, expanded, and then declined at different time points over the two epidemics. In population genetics terms, both hard and soft selective sweeps through vaccine selection pressures have determined the population dynamics of Australian B. pertussis. Relative risk estimation suggests that once a new B. pertussis lineage emerged, it was more likely to spread locally within the first 1.5 years. However, after 1.5 years, any new lineage was likely to expand to a wider region. Phylogenetic analysis revealed the expansion of ptxP3 strains was also associated with replacement of the type III secretion system allele bscI1 with bscI3. bscI3 is associated with decreased T3SS secretion and may allow B. pertussis to reduce immune recognition. This study advanced our understanding of the epidemic population structure and spatial and temporal dynamics of B. pertussis in a highly immunized population.

Evolution of Bordetella pertussis in the acellular vaccine era in Norway, 1996 to 2019

We described the population structure of Bordetella pertussis (B. pertussis) in Norway from 1996 to 2019 and determined if there were evolutionary shifts and whether these correlated with changes in the childhood immunization program. We selected 180 B. pertussis isolates, 22 from the whole cell vaccine (WCV) era (1996-1997) and 158 from the acellular vaccine (ACV) era (1998-2019). We conducted whole genome sequencing and determined the distribution and frequency of allelic variants and temporal changes of ACV genes. Norwegian B. pertussis isolates were evenly distributed across a phylogenetic tree that included global strains. We identified seven different allelic profiles of ACV genes (A-F), in which profiles A1, A2, and B dominated (89%), all having pertussis toxin (ptxA) allele 1, pertussis toxin promoter (ptxP) allele 3, and pertactin (prn) allele 2 present. Isolates with ptxP1 and prn1 were not detected after 2007, whereas the prn2 allele likely emerged prior to 1972, and ptxP3 before the early 1980s. Allele conversions of ACV genes all occurred prior to the introduction of ACV. Sixteen percent of our isolates showed mutations within the prn gene. ACV and its booster doses (implemented for children in 2007 and adolescents in 2013) might have contributed to evolvement of a more uniform B. pertussis population, with recent circulating strains having ptxA1, ptxP3, and prn2 present, and an increasing number of prn mutations. These strains clearly deviate from ACV strains (ptxA1, ptxP1, prn1), and this could have implications for vaccine efficiency and, therefore, prevention and control of pertussis.

Pertussis Vaccine Candidate Based on Outer Membrane Vesicles Derived From Biofilm Culture

Outer membrane vesicles (OMV) derived from Bordetella pertussis-the etiologic agent of the resurgent disease called pertussis-are safe and effective in preventing bacterial colonization in the lungs of immunized mice. Vaccine formulations containing those OMV are capable of inducing a mixed Th1/Th2/Th17 profile, but even more interestingly, they may induce a tissue-resident memory immune response. This immune response is recommended for the new generation of pertussis-vaccines that must be developed to overcome the weaknesses of current commercial acellular vaccines (second-generation of pertussis vaccine). The third-generation of pertussis vaccine should also deal with infections caused by bacteria that currently circulate in the population and are phenotypically and genotypically different [in particular those deficient in the expression of pertactin antigen, PRN(-)] from those that circulated in the past. Here we evaluated the protective capacity of OMV derived from bacteria grown in biofilm, since it was observed that, by difference with older culture collection vaccine strains, circulating clinical B. pertussis isolates possess higher capacity for this lifestyle. Therefore, we performed studies with a clinical isolate with good biofilm-forming capacity. Biofilm lifestyle was confirmed by both scanning electron microscopy and proteomics. While scanning electron microscopy revealed typical biofilm structures in these cultures, BipA, fimbria, and other adhesins described as typical of the biofilm lifestyle were overexpressed in the biofilm culture in comparison with planktonic culture. OMV derived from biofilm (OMVbiof) or planktonic lifestyle (OMVplank) were used to formulate vaccines to compare their immunogenicity and protective capacities against infection with PRN(+) or PRN(-) B. pertussis clinical isolates. Using the mouse protection model, we detected that OMVbiof-vaccine was more immunogenic than OMVplank-vaccine in terms of both specific antibody titers and quality, since OMVbiof-vaccine induced antibodies with higher avidity. Moreover, when OMV were administered at suboptimal quantity for protection, OMVbiof-vaccine exhibited a significantly adequate and higher protective capacity against PRN(+) or PRN(-) than OMVplank-vaccine. Our findings indicate that the vaccine based on B. pertussis biofilm-derived OMV induces high protection also against pertactin-deficient strains, with a robust immune response.

Evaluation of Anti-PT Antibody Response after Pertussis Vaccination and Infection: The Importance of Both Quantity and Quality

Pertussis toxin (PT) is considered the main virulence factor causing whooping cough or pertussis. The protein is widely studied and its composition was revealed and sequenced already during the 1980s. The human immune system creates a good response against PT when measured in quantity. However, the serum anti-PT antibodies wane rapidly, and only a small amount of these antibodies are found a few years after vaccination/infection. Therefore, multiple approaches to study the functionality (quality) of these antibodies, e.g., avidity, neutralizing capacity, and epitope specificity, have been investigated. In addition, the long-term B cell memory (Bmem) to PT is crucial for good protection throughout life. In this review, we summarize the findings from functional PT antibody and Bmem studies. These results are discussed in line with the quantity of serum anti-PT antibodies. PT neutralizing antibodies and anti-PT antibodies with proper avidity are crucial for good protection against the disease, and certain epitopes have been identified to have multiple functions in the protection. Although PT-specific Bmem responses are detectable at least five years after vaccination, long-term surveillance is lacking. Variation of the natural boosting of circulating Bordetella pertussis in communities is an important confounding factor in these memory studies.

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.

Live attenuated pertussis vaccine BPZE1 induces a broad antibody response in humans

BACKGROUNDThe live attenuated BPZE1 vaccine candidate induces protection against B. pertussis and prevents nasal colonization in animal models. Here we report on the responses in humans receiving a single intranasal administration of BPZE1.METHODSWe performed multiple assays to dissect the immune responses induced in humans (n = 12) receiving BPZE1, with particular emphasis on the magnitude and characteristics of the antibody responses. Such responses were benchmarked to adolescents (n = 12) receiving the complete vaccination program of the currently used acellular pertussis vaccine (aPV). Using immunoproteomics analysis, potentially novel immunogenic B. pertussis antigens were identified.RESULTSAll BPZE1 vaccinees showed robust B. pertussis-specific antibody responses with regard to significant increase in 1 or more of the following parameters: IgG, IgA, and memory B cells to B. pertussis antigens. BPZE1-specific T cells showed a Th1 phenotype, and the IgG exclusively consisted of IgG1 and IgG3. In contrast, all aPV vaccines showed a Th2-biased response. Immunoproteomics profiling revealed that BPZE1 elicited broader and different antibody specificities to B. pertussis antigens as compared with the aPV that primarily induced antibodies to the vaccine antigens. Moreover, BPZE1 was superior at inducing opsonizing antibodies that stimulated ROS production in neutrophils and enhanced bactericidal function, which was in line with the finding that antibodies against adenylate cyclase toxin were only elicited by BPZE1.CONCLUSIONThe breadth of the antibodies, the Th1-type cellular response, and killing mechanisms elicited by BPZE1 may hold prospects of improving vaccine efficacy and protection against B. pertussis transmission.TRIAL REGISTRATIONClinicalTrials.gov NCT02453048, NCT00870350.FUNDINGILiAD Biotechnologies, Swedish Research Council (Vetenskapsrådet), Swedish Heart-Lung Foundation.

COVID-19 Evolution in the Post-Vaccination Phase: Endemic or Extinct?

This piece addresses the urge to assess the evolutionary fate of SARS-CoV-2 in a post-vaccination phase. The possibilities of COVID-19 becoming endemic or extinct are weighed against verifiable properties of extant vaccines and observed genetic trends already apparent under the mild selection pressure exerted almost exclusively by social rules.

Nonsynonymous Polymorphism Counts in Bacterial Genomes: a Comparative Examination

Genomic data reveal single-nucleotide polymorphisms (SNPs) that may carry information about the evolutionary history of bacteria. However, it remains unclear what inferences about selection can be made from genomic SNP data. Bacterial species are often sampled during epidemic outbreaks or within hosts during the course of chronic infections. SNPs obtained from genomic analysis of these data are not necessarily fixed. Treating them as fixed during analysis by using measures such as the ratio of nonsynonymous to synonymous evolutionary changes (dN/dS) may lead to incorrect inferences about the strength and direction of selection. In this study, we consider data from a range of whole-genome sequencing studies of bacterial pathogens and explore patterns of nonsynonymous variation to assess whether evidence of selection can be identified by investigating SNP counts alone across multiple WGS studies. We visualize these SNP data in ways that highlight their relationship to neutral baseline expectations. These neutral expectations are based on a simple model of mutation, from which we simulate SNP accumulation to investigate how SNP counts are distributed under alternative assumptions about positive and negative selection. We compare these patterns with empirical SNP data and illustrate the general difficulty of detecting positive selection from SNP data. Finally, we consider whether SNP counts observed at the between-host population level differ from those observed at the within-host level and find some evidence that suggests that dynamics across these two scales are driven by different underlying processes.IMPORTANCE Identifying selection from SNP data obtained from whole-genome sequencing studies is challenging. Some current measures used to identify and quantify selection acting on genomes rely on fixed differences; thus, these are inappropriate for SNP data where variants are not fixed. With the increase in whole-genome sequencing studies, it is important to consider SNP data in the context of evolutionary processes. How SNPs are counted and analyzed can help in understanding mutation accumulation and trajectories of strains. We developed a tool for identifying possible evidence of selection and for comparative analysis with other SNP data. We propose a model that provides a rule-of-thumb guideline and two new visualization techniques that can be used to interpret and compare SNP data. We quantify the expected proportion of nonsynonymous SNPs in coding regions under neutrality and demonstrate its use in identifying evidence of positive and negative selection from simulations and empirical data.

Monitor for COVID-19 vaccine resistance evolution during clinical trials

Although less common than the evolution of antimicrobial drug resistance, vaccine resistance can and has evolved. How likely is it that COVID-19 vaccines currently in development will be undermined by viral evolution? We argue that this can be determined by repurposing samples that are already being collected as part of clinical trials. Such information would be useful for prioritizing investment among candidate vaccines and maximizing the potential long-term impact of COVID-19 vaccines.

Next-Generation Pertussis Vaccines Based on the Induction of Protective T Cells in the Respiratory Tract

Immunization with current acellular pertussis (aP) vaccines protects against severe pertussis, but immunity wanes rapidly after vaccination and these vaccines do not prevent nasal colonization with Bordetella pertussis. Studies in mouse and baboon models have demonstrated that Th1 and Th17 responses are integral to protective immunity induced by previous infection with B. pertussis and immunization with whole cell pertussis (wP) vaccines. Mucosal Th17 cells, IL-17 and secretory IgA (sIgA) are particularly important in generating sustained sterilizing immunity in the nasal cavity. Current aP vaccines induce potent IgG and Th2-skewed T cell responses but are less effective at generating Th1 and Th17 responses and fail to prime respiratory tissue-resident memory T (TRM) cells, that maintain long-term immunity at mucosal sites. In contrast, a live attenuated pertussis vaccine, pertussis outer membrane vesicle (OMV) vaccines or aP vaccines formulated with novel adjuvants do induce cellular immune responses in the respiratory tract, especially when delivered by the intranasal route. An increased understanding of the mechanisms of sustained protective immunity, especially the role of respiratory TRM cells, will facilitate the development of next generation pertussis vaccines that not only protect against pertussis disease, but prevent nasal colonization and transmission of B. pertussis.

Comparative genomics of whole-cell pertussis vaccine strains from India

BACKGROUND

Despite high vaccination coverage using acellular (ACV) and whole-cell pertussis (WCV) vaccines, the resurgence of pertussis is observed globally. Genetic divergence in circulating strains of Bordetella pertussis has been reported as one of the contributing factors for the resurgence of the disease. Our current knowledge of B. pertussis genetic evolution in circulating strains is mostly based on studies conducted in countries using ACVs targeting only a few antigens used in the production of ACVs. To better understand the adaptation to vaccine-induced selection pressure, it will be essential to study B. pertussis populations in developing countries which are using WCVs. India is a significant user and global supplier of WCVs. We report here comparative genome analyses of vaccine and clinical isolates reported from India. Whole-genome sequences obtained from vaccine strains: WCV (J445, J446, J447 and J448), ACV (BP165) were compared with Tohama-I reference strain and recently reported clinical isolates from India (BPD1, BPD2). Core genome-based phylogenetic analysis was also performed using 166 isolates reported from countries using ACV.

RESULTS

Whole-genome analysis of vaccine and clinical isolates reported from India revealed high genetic similarity and conserved genome among strains. Phylogenetic analysis showed that clinical and vaccine strains share genetic closeness with reference strain Tohama-I. The allelic profile of vaccine strains (J445:ptxP1/ptxA2/prn1/fim2-1/fim3-1; J446: ptxP2/ptxA4/prn7/fim2-2/fim3-1; J447 and J448: ptxP1/ptxA1/ prn1/fim2-1/fim3-1), which matched entirely with clinical isolates (BPD1:ptxP1/ptxA1/prn1/fim2-1 and BPD2: ptxP1/ptxA1/prn1/fim2-1) reported from India. Multi-locus sequence typing (MLST) demonstrated the presence of dominant sequence types ST2 and primitive ST1 in vaccine strains which will allow better coverage against circulating strains of B. pertussis.

CONCLUSIONS

The study provides a detailed characterization of vaccine and clinical strains reported from India, which will further facilitate epidemiological studies on genetic shifts in countries which are using WCVs in their immunization programs.

Pertactin-Negative and Filamentous Hemagglutinin-Negative Bordetella pertussis, Australia, 2013-2017

During the 2008–2012 pertussis epidemic in Australia, pertactin (Prn)–negative Bordetella pertussis emerged. We analyzed 78 isolates from the 2013–2017 epidemic and documented continued expansion of Prn-negative ptxP3 B. pertussis strains. We also detected a filamentous hemagglutinin-negative and Prn-negative B. pertussis isolate.

Immunomodulation as a Novel Strategy for Prevention and Treatment of Bordetella spp. Infections

Well-adapted pathogens have evolved to survive the many challenges of a robust immune response. Defending against all host antimicrobials simultaneously would be exceedingly difficult, if not impossible, so many co-evolved organisms utilize immunomodulatory tools to subvert, distract, and/or evade the host immune response. Bordetella spp. present many examples of the diversity of immunomodulators and an exceptional experimental system in which to study them. Recent advances in this experimental system suggest strategies for interventions that tweak immunity to disrupt bacterial immunomodulation, engaging more effective host immunity to better prevent and treat infections. Here we review advances in the understanding of respiratory pathogens, with special focus on Bordetella spp., and prospects for the use of immune-stimulatory interventions in the prevention and treatment of infection.

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.

Surfaceome analysis of Australian epidemic Bordetella pertussis reveals potential vaccine antigens

Since acellular vaccines (ACV) were introduced in Australia, epidemic Bordetella pertussis strains changed from single nucleotide polymorphism (SNP) cluster II to SNP cluster I. Our previous proteomic analysis identified potential proteomic adaptations in the whole cell and secretome of SNP cluster I. Additionally, current ACVs were shown to be less efficacious against cluster I in mice models and there is a pressing need to discover new antigens to improve the ACV. One important source of novel antigens is the surfaceome. Therefore, in this study we established surface shaving in B. pertussis to compare the surfaceome of SNP cluster I (L1423) and II (L1191), and identify novel surface antigens for vaccine development. Surface shaving using 1 μg of trypsin for 5 min identified 126 proteins with the most abundant being virulence-associated and known outer membrane proteins. Cell viability counts showed minimal lysis from shaving. The proportion of immunogenic proteins was higher in the surfaceome than in the whole cell and secretome. Key differences in the surfaceome were identified between SNP cluster I and II, consistent with those identified in the whole cell proteome and secretome. These differences include unique transport proteins and decreased immunogenic proteins in L1423, and provides further evidence of proteomic adaptation in SNP cluster I. Finally, a comparison of proteins in each sub-proteome identified 22 common proteins. These included 11 virulence proteins (Prn, PtxA, FhaB, CyaA, TcfA, SphB1, Vag8, BrkA, BopD, Bsp22 and BipA) and 11 housekeeping proteins (TuF, CtpA, TsF, OmpH, GltA, SucC, SucD, FusA, GroEL, BP3330 and BP3561) which were immunogenic, essential and consistently expressed thus demonstrating their potential as future targets. This study established surface shaving in B. pertussis, confirmed key expression differences and identified unknown surface proteins which may be potential vaccine antigens.

Genomic epidemiology of erythromycin-resistant Bordetella pertussis in China

Macrolides such as erythromycin are the empirical treatment of Bordetella pertussis infections. China has experienced an increase in erythromycin-resistant B. pertussis isolates since they were first reported in 2013. Here, we undertook a genomic study on Chinese B. pertussis isolates from 2012 to 2015 to elucidate the origins and phylogenetic relationships of erythromycin-resistant B. pertussis isolates in China. A total of 167 Chinese B. pertussis isolates were used for antibiotic sensitivity testing and multiple locus variable-number tandem repeat (VNTR) analysis (MLVA). All except four isolates were erythromycin-resistant and of the four erythromycin-sensitive isolates, three were non-ptxP1. MLVA types (MT), MT55, MT104 and MT195 were the predominant types. Fifty of those isolates were used for whole genome sequencing and phylogenetic analysis. Genome sequencing and phylogenetic analysis revealed three independent erythromycin-resistant lineages and all resistant isolates carried a mutation in the 23S rRNA gene. A novel fhaB3 allele was found uniquely in Chinese ptxP1 isolates and these Chinese ptxP1-ptxA1-fhaB3 had a 5-fold higher mutation rate than the global ptxP1-ptxA1 B. pertussis population. Our results suggest that the evolution of Chinese B. pertussis is likely to be driven by selection pressure from both vaccination and antibiotics. The emergence of the new non-vaccine fhaB3 allele in Chinese B. pertussis population may be a result of selection from vaccination, whereas the expansion of ptxP1-fhaB3 lineages was most likely to be the result of selection pressure from antibiotics. Further monitoring of B. pertussis in China is required to better understand the evolution of the pathogen.

Immunogenicity and protective potential of Bordetella pertussis biofilm and its associated antigens in a murine model

The resurgence of whooping cough reflects novel genetic variants of Bordetella pertussis and inadequate protection conferred by current acellular vaccines (aP). Biofilm is a source of novel vaccine candidates, including membrane protein assembly factor (BamB) and lipopolysaccharide assembly protein (LptD). Responses of BALB/c mice to candidate vaccines included IFN-γ and IL-17a production by spleen and lymph node cells, and serum IgG1 and IgG2a reactive with whole bacteria or aP. Protection was determined using bacterial cultured from lungs of vaccinated mice challenged with virulent B. pertussis. Mice vaccinated with biofilm produced efficient IFN-γ responses and more IL-17a and IgG2a than mice vaccinated with planktonic cells, aP or adjuvant alone. Vaccination with aP produced abundant IgG1 with little IgG2a. Mice vaccinated with aP plus BamB and LptD retained lower bacterial loads than mice vaccinated with aP alone. Whooping cough vaccines formulated with biofilm antigens, including BamB and LptD, may have clinical value.

A Review on T Cell Epitopes Identified Using Prediction and Cell-Mediated Immune Models for Mycobacterium tuberculosis and Bordetella pertussis

In the present review, we summarize work from our as well as other groups related to the characterization of bacterial T cell epitopes, with a specific focus on two important pathogens, namely, Mycobacterium tuberculosis (Mtb), the bacterium that causes tuberculosis (TB), and Bordetella pertussis (BP), the bacterium that causes whooping cough. Both bacteria and their associated diseases are of large societal significance. Although vaccines exist for both pathogens, their efficacy is incomplete. It is widely thought that defects and/or alteration in T cell compartments are associated with limited vaccine effectiveness. As discussed below, a full genome-wide map was performed in the case of Mtb. For BP, our focus has thus far been on the antigens contained in the acellular vaccine; a full genome-wide screen is in the planning stage. Nevertheless, the sum-total of the results in the two different bacterial systems allows us to exemplify approaches and techniques that we believe are generally applicable to the mapping and characterization of human immune responses to bacterial pathogens. Finally, we add, as a disclaimer, that this review by design is focused on the work produced by our laboratory as an illustration of approaches to the study of T cell responses to Mtb and BP, and is not meant to be comprehensive, nor to detract from the excellent work performed by many other groups.

Comparison of the Whole Cell Proteome and Secretome of Epidemic Bordetella pertussis Strains From the 2008-2012 Australian Epidemic Under Sulfate-Modulating Conditions

Sulfate is an important modulator for virulence factor expression in Bordetella pertussis, the causative organism for whooping cough. During infection, sulfate is released when respiratory epithelial cells are damaged which can affect gene expression. The current predominant strains in Australia are found in single nucleotide polymorphism (SNP) cluster I (ptxP3/prn2). It has been reported that ptxP3 strains have higher mRNA expression of virulence genes than ptxP1 strains under intermediate sulfate-modulating conditions (5 mM MgSO₄). Our previous proteomic study compared L1423 (cluster I, ptxP3) and L1191 (cluster II, ptxP1) in Thalen-IJssel (THIJS) media without sulfate modulation and identified an upregulation of transport proteins and a downregulation of immunogenic proteins. To determine whether proteomic differences exist between cluster I and cluster II strains in intermediate modulating conditions, this study compared the whole cell proteome and secretome between L1423 and L1191 grown in THIJS media with 5 mM MgSO₄ using iTRAQ and high-resolution multiple reaction monitoring (MRM-hr). Two proteins (BP0200 and BP1175) in the whole cell were upregulated in L1423 [fold change (FC) >1.2, false discovery rate (FDR) <0.05]. In the secretome, four proteins from the type III secretion system (T3SS) effectors were downregulated (FC < 0.8, FDR < 0.05) while six proteins, including two adhesins, pertactin (Prn) and tracheal colonization factor A (TcfA), were upregulated which were consistent with our previous proteomic study. The upregulation of Prn and TcfA in SNP cluster I may result in improved adhesion while the downregulation of the T3SS and other immunogenic proteins may reduce immune recognition, which may contribute to the increased fitness of cluster I B. pertussis strains.

Proteomic Adaptation of Australian Epidemic Bordetella pertussis

Bordetella pertussis causes whooping cough. The predominant strains in Australia changed to single nucleotide polymorphism (SNP) cluster I (pertussis toxin promoter allele ptxP3/pertactin gene allele prn2) from cluster II (non-ptxP3/non-prn2). Cluster I was mostly responsible for the 2008-2012 Australian epidemic and was found to have higher fitness compared to cluster II using an in vivo mouse competition assay, regardless of host's immunization status. This study aimed to identify proteomic differences that explain higher fitness in cluster I using isobaric tags for relative and absolute quantification (iTRAQ), and high-resolution multiple reaction monitoring (MRM-hr). A few key differences in the whole cell and secretome were identified between the cluster I and II strains tested. In the whole cell, nine proteins were upregulated (>1.2 fold change, q < 0.05) and three were downregulated (<0.8 fold change, q < 0.05) in cluster I. One downregulated protein was BP1569, a TLR2 agonist for Th1 immunity. In the secretome, 12 proteins were upregulated and 1 was downregulated which was Bsp22, a type III secretion system (T3SS) protein. Furthermore, there was a trend of downregulation in three T3SS effectors and other virulence factors. Three proteins were upregulated in both whole cell and supernatant: BP0200, molybdate ABC transporter (ModB), and tracheal colonization factor A (TcfA). Important expression differences in lipoprotein, T3SS, and transport proteins between the cluster I and II strains were identified. These differences may affect immune evasion, virulence and metabolism, and play a role in increased fitness of cluster I.

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.

Multi-locus variable-number tandem repeat analysis of Bordetella pertussis isolates circulating in Poland in the period 1959-2013

PURPOSE

Despite the long history of pertussis vaccination and high vaccination coverage in Poland and many other developed countries, pertussis incidence rates have increased substantially, making whooping cough one of the most prevalent vaccine-preventable diseases. Among the factors potentially involved in pertussis resurgence, the adaptation of the Bordetella pertussis population to country-specific vaccine-induced immunity through selection of non-vaccine-type strains still needs detailed studies.

METHODOLOGY

Multi-locus variable-number tandem repeat analysis (MLVA), also linked to MLST and PFGE profiling, was applied to trace the genetic changes in the B. pertussis population circulating in Poland in the period 1959-2013 versus country-specific vaccine strains.

RESULTS

Generally, among 174 B. pertussis isolates, 31 MLVA types were detected, of which 11 were not described previously. The predominant MLVA types of recent isolates in Poland were different from those of the typical isolates circulating in other European countries. The MT27 type, currently predominant in Europe, was rarely seen and detected in only five isolates among all studied. The features of the vaccine strains used for production of the pertussis component of a national whole-cell diphtheria-tetanus-pertussis (DTP) vaccine, as studied by MLVA and MLST tools, were found to not match those observed in the currently circulating B. pertussis isolates in Poland.

CONCLUSIONS

Differences traced by MLVA in relation to the MLST and PFGE profiling confirmed that the B. pertussis strain types currently observed elsewhere in Europe, even if appearing in Poland, were not able to successfully disseminate within a human population in Poland that has been vaccinated with a whole-cell pertussis vaccine not used in other countries.

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.

Complete Genome Sequences of Four Bordetella pertussis Vaccine Reference Strains from Serum Institute of India

Serum Institute of India is among the world’s largest vaccine producers. Here, we report the complete genome sequences for four Bordetella pertussis strains used by Serum Institute of India in the production of whole-cell pertussis vaccines.

Pertussis: Biology, epidemiology and prevention

Despite long-standing vaccination programs, substantial increases in reported cases of pertussis have been described in several countries during the last 5years. Cases among very young infants who are at greatest risk of pertussis-related hospitalizations and mortality are the most alarming. Multiple hypotheses including but not limited to the availability of more sensitive diagnostic tests, greater awareness, and waning vaccine-induced immunity over time have been posited for the current challenges with pertussis. The conference "Pertussis: biology, epidemiology and prevention" held in Annecy-France (November 11-13, 2015) brought together experts and interested individuals to examine these issues and to formulate recommendations for optimal use of current vaccines, with a particular focus on strategies to minimize severe morbidity and mortality among infants during the first months of life. The expert panel concluded that improving vaccination strategies with current vaccines and development of new highly immunogenic and efficacious pertussis vaccines that have acceptable adverse event profiles are currently the two main areas of investigation for the control of pertussis. Some possible pathways forward to address these main challenges are discussed in this report.

Biofilm forming potential and antimicrobial susceptibility of newly emerged Western Australian Bordetella pertussis clinical isolates

Whooping cough caused by Bordetella pertussis is increasing in several countries despite high vaccine coverage. One potential reason for the resurgence is the emergence of genetic variants of the bacterium. Biofilm formation has recently been associated with the pathogenesis of B. pertussis. Biofilm formation of 21 Western Australian B. pertussis clinical isolates was investigated. All isolates formed thicker biofilms than the reference vaccine strain Tohama I while retaining susceptibility to ampicillin, erythromycin, azithromycin and streptomycin. When two biofilm-forming clinical isolates were compared with Tohama I, minimum bactericidal concentrations of antimicrobial agents increased. Isobaric tags for relative and absolute quantitation (iTRAQ)-based proteomic analysis revealed significant differences in protein expression in B. pertussis biofilms, providing an opportunity for identification of novel biofilm-associated antigens for incorporation in current pertussis vaccines to improve their protective efficacy. The study also highlights the importance of determining antibiograms for biofilms to formulate improved antimicrobial therapeutic regimens.

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

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.

Better colonisation of newly emerged Bordetella pertussis in the co-infection mouse model study

Molecular epidemiological data indicates that the resurgence of pertussis (whooping cough) in populations with high vaccine coverage is associated with genomic adaptation of Bordetella pertussis, the causative agent of the disease, to vaccine selection pressure. We have previously shown that in the period after the introduction of acellular pertussis vaccine (ACV), the majority of circulating strains in Australia switched to single nucleotide polymorphism (SNP) cluster I (carrying ptxP3/prn2), replacing SNP cluster II (carrying ptxP1/prn3). In this study, we carried out an in vivo competition assay using a mouse model infected with SNP cluster I and II B. pertussis strains from Australia. We found that the SNP cluster I strain colonised better than the SNP cluster II strain, in both naïve and immunised mice, suggesting that SNP cluster I strains had better fitness regardless of immunisation status of the host, consistent with SNP cluster I strains replacing SNP cluster II. Nevertheless, we found that ACV enhanced clearance of both SNP cluster I and II strains from the mouse respiratory tract.

Genomic dissection of Australian Bordetella pertussis isolates from the 2008-2012 epidemic

OBJECTIVES

Despite high pertussis vaccination coverage, Australia experienced a prolonged epidemic in 2008-2012. The predominant Bordetella pertussis genotype harboured pertussis toxin promoter allele, ptxP3, and pertactin gene allele, prn2. The emergence and expansion of prn non-expressing isolates (Prn negative), were also observed. We aimed to investigate the microevolution and genomic diversity of epidemic B. pertussis isolates.

METHODS

We sequenced 22 B. pertussis isolates collected in 2008-2012 from two states of Australia which are geographically widely separated. Ten of the 22 were Prn negative isolates with three different modes of silencing of prn (prn::IS481F, prn::IS481R and prn::IS1002). Five pre-epidemic isolates were also sequenced for comparison.

RESULTS

Five single nucleotide polymorphisms were common in the epidemic isolates and differentiated them from pre-epidemic isolates. The Australian epidemic isolates can be divided into five lineages (EL1-EL5) with EL1 containing only Prn negative isolates. Comparison with global isolates showed that three lineages remained geographically and temporally distinct whereas two lineages mixed with isolates from 2012 UK outbreak.

CONCLUSION

Our results suggest significant diversification and the microevolution of B. pertussis within the 2008-2012 Australian epidemic.

Transcriptional profiling of Bordetella pertussis reveals requirement of RNA chaperone Hfq for Type III secretion system functionality

Bordetella pertussis, the causative agent of human whooping cough (pertussis) produces a complex array of virulence factors in order to establish efficient infection in the host. The RNA chaperone Hfq and small regulatory RNAs are key players in posttranscriptional regulation in bacteria and have been shown to play an essential role in virulence of a broad spectrum of bacterial pathogens. This study represents the first attempt to characterize the Hfq regulon of the human pathogen B. pertussis under laboratory conditions as well as upon passage in the host and indicates that loss of Hfq has a profound effect on gene expression in B. pertussis. Comparative transcriptional profiling revealed that Hfq is required for expression of several virulence factors in B. pertussis cells including the Type III secretion system (T3SS). In striking contrast to the wt strain, T3SS did not become operational in the hfq mutant passaged either through mice or macrophages thereby proving that Hfq is required for the functionality of the B. pertussis T3SS. Likewise, expression of virulence factors vag8 and tcfA encoding autotransporter and tracheal colonization factor, respectively, was strongly reduced in the hfq mutant. Importantly, for the first time we demonstrate that B. pertussis T3SS can be activated upon contact with macrophage cells in vitro.

Pertactin negative Bordetella pertussis demonstrates higher fitness under vaccine selection pressure in a mixed infection model

Whooping cough or pertussis is a highly infectious respiratory disease in humans caused by Bordetella pertussis. The use of acellular vaccines (ACV) has been associated with the recent resurgence of pertussis in developed countries including Australia despite high vaccination coverage where B. pertussis strains that do not express pertactin (Prn), a key antigenic component of the ACV, have emerged and become prevalent. In this study, we used an in vivo competition assay in mice immunised with ACV and in naïve (control) mice to compare the proportion of colonisation with recent clinical Prn positive and Prn negative B. pertussis strains from Australia. The Prn negative strain colonised the respiratory tract more effectively than the Prn positive strain in immunised mice, out-competing the Prn positive strain by day 3 of infection. However, in control mice, the Prn positive strain out-competed the Prn negative strain. Our findings of greater ability of Prn negative strains to colonise ACV-immunised mice are consistent with reports of selective advantage for these strains in ACV-immunised humans.

The adaptive evolution of virulence: a review of theoretical predictions and empirical tests

Why is it that some parasites cause high levels of host damage (i.e. virulence) whereas others are relatively benign? There are now numerous reviews of virulence evolution in the literature but it is nevertheless still difficult to find a comprehensive treatment of the theory and data on the subject that is easily accessible to non-specialists. Here we attempt to do so by distilling the vast theoretical literature on the topic into a set of relatively few robust predictions. We then provide a comprehensive assessment of the available empirical literature that tests these predictions. Our results show that there have been some notable successes in integrating theory and data but also that theory and empiricism in this field do not ‘speak’ to each other very well. We offer a few suggestions for how the connection between the two might be improved.

Investigating pertussis toxin and its impact on vaccination

Whooping cough, caused by Bordetella pertussis, remains a major global health problem. Each year around 40 million of pertussis cases resulting in 200,000-400,000 annual deaths occur worldwide. Pertussis toxin is a major virulence factor of B. pertussis. Murine studies have shown its importance in bacterial colonization and in immunomodulation to evade innate or adaptive immunity. The toxin is composed of an A protomer expressing ADP-ribosyltransferase activity and a B oligomer, responsible for toxin binding to target cells. The toxin is also a major protective antigen in all currently available vaccines. However, vaccine escape mutants with altered toxin expression have recently been isolated in countries with high vaccination coverage illustrating the need for improved pertussis vaccines.

Dynamics of Pertussis Transmission in the United States

Past patterns of infectious disease transmission set the stage on which modern epidemiologic dynamics are played out. Here, we present a comprehensive account of pertussis (whooping cough) transmission in the United States during the early vaccine era. We analyzed recently digitized weekly incidence records from Morbidity and Mortality Weekly Reports from 1938 to 1955, when the whole-cell pertussis vaccine was rolled out, and related them to contemporary patterns of transmission and resurgence documented in monthly incidence data from the National Notifiable Diseases Surveillance System. We found that, during the early vaccine era, pertussis epidemics in US states could be categorized as 1) annual, 2) initially annual and later multiennial, or 3) multiennial. States with predominantly annual cycles tended to have higher per capita birth rates, more household crowding, more children per family, and lower rates of school attendance than the states with multiennial cycles. Additionally, states that exhibited annual epidemics during 1938-1955 have had the highest recent (2001-2010) incidence, while those states that transitioned from annual cycles to multiennial cycles have had relatively low recent incidence. Our study provides an extensive picture of pertussis epidemiology in the United States dating back to the onset of vaccination, a back-story that could aid epidemiologists in understanding contemporary transmission patterns.

Studying Bordetella pertussis populations by use of SNPeX, a simple high-throughput single nucleotide polymorphism typing method

Large outbreaks of pertussis occur despite vaccination. A first step in the analyses of outbreaks is strain typing. However, the typing of Bordetella pertussis, the causative agent of pertussis, is problematic because the available assays are insufficiently discriminatory, not unequivocal, time-consuming, and/or costly. Here, we describe a single nucleotide primer extension assay for the study of B. pertussis populations, SNPeX (single nucleotide primer extension), which addresses these problems. The assay is based on the incorporation of fluorescently labeled dideoxynucleotides (ddNTPs) at the 3' end of allele-specific poly(A)-tailed primers and subsequent analysis with a capillary DNA analyzer. Each single nucleotide polymorphism (SNP) primer has a specific length, and as a result, up to 20 SNPs can be determined in one SNPeX reaction. Importantly, PCR amplification of target DNA is not required. We selected 38 SNPeX targets from the whole-genome sequencing data of 74 B. pertussis strains collected from across the world. The SNPeX-based phylogenetic trees preserved the general tree topology of B. pertussis populations based on whole-genome sequencing, with a minor loss of details. We envisage a strategy whereby SNP types (SnpTs) are quickly identified with the SNPeX assay during an outbreak, followed by whole-genome sequencing (WGS) of a limited number of isolates representing predominant SnpTs and the incorporation of novel SNPs in the SNPeX assay. The flexibility of the SNPeX assay allows the method to evolve along with the pathogen, making it a promising method for studying outbreaks of B. pertussis and other pathogens.

Analysis of Bordetella pertussis clinical isolates circulating in European countries during the period 1998-2012

Despite more than 50 years of vaccination, pertussis is still an endemic disease, with regular epidemic outbreaks. With the exception of Poland, European countries have replaced whole-cell vaccines (WCVs) by acellular vaccines (ACVs) in the 1990s. Worldwide, antigenic divergence in vaccine antigens has been found between vaccine strains and circulating strains. In this work, 466 Bordetella pertussis isolates collected in the period 1998–2012 from 13 European countries were characterised by multi-locus antigen sequence typing (MAST) of the pertussis toxin promoter (ptxP) and of the genes coding for proteins used in the ACVs: pertussis toxin (Ptx), pertactin (Prn), type 2 fimbriae (Fim2) and type 3 fimbriae (Fim3). Isolates were further characterised by fimbrial serotyping, multi-locus variable-number tandem repeat analysis (MLVA) and pulsed-field gel electrophoresis (PFGE). The results showed a very similar B. pertussis population for 12 countries using ACVs, while Poland, which uses a WCV, was quite distinct, suggesting that ACVs and WCVs select for different B. pertussis populations. This study forms a baseline for future studies on the effect of vaccination programmes on B. pertussis populations.

Rapid increase in pertactin-deficient Bordetella pertussis isolates, Australia

Acellular vaccines against Bordetella pertussis were introduced in Australia in 1997. By 2000, these vaccines had replaced whole-cell vaccines. During 2008–2012, a large outbreak of pertussis occurred. During this period, 30% (96/320) of B. pertussis isolates did not express the vaccine antigen pertactin (prn). Multiple mechanisms of prn inactivation were documented, including IS481 and IS1002 disruptions, a variation within a homopolymeric tract, and deletion of the prn gene. The mechanism of lack of expression of prn in 16 (17%) isolates could not be determined at the sequence level. These findings suggest that B. pertussis not expressing prn arose independently multiple times since 2008, rather than by expansion of a single prn-negative clone. All but 1 isolate had ptxA1, prn2, and ptxP3, the alleles representative of currently circulating strains in Australia. This pattern is consistent with continuing evolution of B. pertussis in response to vaccine selection pressure.

Investigating genome reduction of Bordetella pertussis using a multiplex PCR-based reverse line blot assay (mPCR/RLB)

Background

The genetic composition of the bacterium causing whooping cough, Bordetella pertussis, has been investigated using microarray studies in order to examine potential genetic contributors to the disease re-emergence in the past decade. Regions of difference (RDs) have been previously identified as clusters of genes flanked by insertion sequences which are variably present in different sets of isolates, and have also been shown to be potential markers of B. pertussis evolution.

This study used microarray data to identify and select a panel of RDs; primers and probes for these RDs were then designed to test for the presence or absence of these regions in a novel and less expensive multiplex PCR-based reverse line blot (mPCR/RLB) assay. By comparing the presence or absence of RDs, we aimed to determine the genomic variability of a diverse collection of B. pertussis strains and how they have changed over time.

Results

A B. pertussis specific mPCR/RLB using 43 genes representing 30 RDs, was developed and used to characterise a set of 42 B. pertussis isolates. When mapped against the previously identified evolutionary relationships of the strains, the losses of two RDs - BP0910A - BP00930 and BP1948-BP1962 - were found to be associated with significant events in B. pertussis history: the loss of BP0910A - BP00930 coincided with introduction of whole cell vaccines in the 1950s while that of BP1948-BP1962 occurred after the introduction of acellular vaccines. The loss of BP1948-BP1962 also coincided with expansion of the most recent B. pertussis strains.

Conclusions

The mPCR/RLB assay offers an inexpensive and fast method of determining the gene content of B. pertussis strains and also confirms that gene losses are an ongoing feature of B. pertussis evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/1756-0500-7-727) contains supplementary material, which is available to authorized users.

Global population structure and evolution of Bordetella pertussis and their relationship with vaccination

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.

Antibody-mediated inhibition of Bordetella pertussis adenylate cyclase-haemolysin-induced macrophage cytotoxicity is influenced by variations in the bacterial population

Whooping cough is a vaccine-preventable disease presenting with epidemic cycles linked to natural and/or vaccine-driven evolution of the aetiological agent of the disease, Bordetella pertussis. Adenylate cyclase-haemolysin (AC-Hly) is a major toxin produced by this pathogen, which mediates macrophage apoptosis in vitro and in vivo. While current acellular pertussis vaccine (APV) formulations do not include AC-Hly, they all contain pertussis toxin and can comprise filamentous haemagglutinin (FHA), which interacts with AC-Hly, and pertactin (PRN), which has been hypothesized also to interact with AC-Hly. We aimed to study the capacity of specific antibodies to inhibit the in vitro B. pertussis AC-Hly-mediated cytotoxicity of J774A.1 murine macrophages in a background of a changing bacterial population. We demonstrate that: (i) clinical isolates of different types or PRN phenotype are all cytotoxic and lethal in the mouse model of respiratory infection; (ii) lack of PRN production does not impact AC-Hly-related phenotypes; (iii) anti-AC-Hly antibodies inhibit cell lysis whatever the phenotype of the isolate, while anti-PRN antibodies significantly inhibit cell lysis provided the isolate produces this antigen, which might be relevant in vivo for APV-induced immunity; and (iv) anti-FHA antibodies only inhibit lysis induced by isolates collected in 2012, maybe indicating specific characteristics of epidemic lineages of B. pertussis.

Substantial gaps in knowledge of Bordetella pertussis antibody and T cell epitopes relevant for natural immunity and vaccine efficacy

The recent increase in whooping cough in vaccinated populations has been attributed to waning immunity associated with the acellular vaccine. The Immune Epitope Database (IEDB) is a repository of immune epitope data from the published literature and includes T cell and antibody epitopes for human pathogens. The IEDB conducted a review of the epitope literature, which revealed 300 Bordetella pertussis-related epitopes from 39 references. Epitope data are currently available for six virulence factors of B. pertussis: pertussis toxin, pertactin, fimbrial 2, fimbrial 3, adenylate cyclase and filamentous hemagglutinin. The majority of epitopes were defined for antibody reactivity; fewer T cell determinants were reported. Analysis of available protective correlates data revealed a number of candidate epitopes; however few are defined in humans and few have been shown to be protective. Moreover, there are a limited number of studies defining epitopes from natural infection versus whole cell or acellular/subunit vaccines. The relationship between epitope location and structural features, as well as antigenic drift (SNP analysis) was also investigated. We conclude that the cumulative data is yet insufficient to address many fundamental questions related to vaccine failure and this underscores the need for further investigation of B. pertussis immunity at the molecular level.

Pertussis resurgence: waning immunity and pathogen adaptation - two sides of the same coin

Pertussis or whooping cough has persisted and resurged in the face of vaccination and has become one of the most prevalent vaccine-preventable diseases in Western countries. The high circulation rate of Bordetella pertussis poses a threat to infants that have not been (completely) vaccinated and for whom pertussis is a severe, life-threatening, disease. The increase in pertussis is mainly found in age groups in which immunity has waned and this has resulted in the perception that waning immunity is the main or exclusive cause for the resurgence of pertussis. However, significant changes in B. pertussis populations have been observed after the introduction of vaccinations, suggesting a role for pathogen adaptation in the persistence and resurgence of pertussis. These changes include antigenic divergence with vaccine strains and increased production of pertussis toxin. Antigenic divergence will affect both memory recall and the efficacy of antibodies, while higher levels of pertussis toxin may increase suppression of the innate and acquired immune system. We propose these adaptations of B. pertussis have decreased the period in which pertussis vaccines are effective and thus enhanced the waning of immunity. We plead for a more integrated approach to the pertussis problem which includes the characteristics of the vaccines, the B. pertussis populations and the interaction between the two.