Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Effect of FHA and Prn on Bordetella pertussis colonization of mice is dependent on vaccine type and anatomical site

  • Anne Zeddeman ,

    Contributed equally to this work with: Anne Zeddeman, Evi van Schuppen

    Roles Conceptualization, Formal analysis, Investigation, Writing – original draft, Writing – review & editing

    Affiliations Section Pediatric Infectious Diseases, Laboratory of Medical Immunology, Radboud Institute for Molecular Life Sciences, Radboudumc, Nijmegen, The Netherlands, Radboud Center for Infectious Diseases, Radboudumc, Nijmegen, The Netherlands, Centre for Infectious Diseases Research, Diagnostics and Screening (IDS), National Institute of Public Health and the Environment (RIVM), Bilthoven, The Netherlands

  • Evi van Schuppen ,

    Contributed equally to this work with: Anne Zeddeman, Evi van Schuppen

    Roles Conceptualization, Formal analysis, Investigation, Writing – original draft, Writing – review & editing

    Affiliations Section Pediatric Infectious Diseases, Laboratory of Medical Immunology, Radboud Institute for Molecular Life Sciences, Radboudumc, Nijmegen, The Netherlands, Radboud Center for Infectious Diseases, Radboudumc, Nijmegen, The Netherlands

  • Kristianne E. Kok,

    Roles Investigation, Writing – review & editing

    Affiliation Centre for Infectious Diseases Research, Diagnostics and Screening (IDS), National Institute of Public Health and the Environment (RIVM), Bilthoven, The Netherlands

  • Marjolein van Gent,

    Roles Investigation, Writing – review & editing

    Affiliation Centre for Infectious Diseases Research, Diagnostics and Screening (IDS), National Institute of Public Health and the Environment (RIVM), Bilthoven, The Netherlands

  • Kees J. Heuvelman,

    Roles Formal analysis, Investigation, Writing – review & editing

    Affiliation Centre for Infectious Diseases Research, Diagnostics and Screening (IDS), National Institute of Public Health and the Environment (RIVM), Bilthoven, The Netherlands

  • Marieke J. Bart,

    Roles Formal analysis, Writing – review & editing

    Affiliation Centre for Infectious Diseases Research, Diagnostics and Screening (IDS), National Institute of Public Health and the Environment (RIVM), Bilthoven, The Netherlands

  • Han G. J. van der Heide,

    Roles Formal analysis, Writing – review & editing

    Affiliation Centre for Infectious Diseases Research, Diagnostics and Screening (IDS), National Institute of Public Health and the Environment (RIVM), Bilthoven, The Netherlands

  • Joshua Gillard,

    Roles Formal analysis, Writing – review & editing

    Affiliations Section Pediatric Infectious Diseases, Laboratory of Medical Immunology, Radboud Institute for Molecular Life Sciences, Radboudumc, Nijmegen, The Netherlands, Radboud Center for Infectious Diseases, Radboudumc, Nijmegen, The Netherlands

  • Elles Simonetti,

    Roles Investigation, Writing – review & editing

    Affiliations Section Pediatric Infectious Diseases, Laboratory of Medical Immunology, Radboud Institute for Molecular Life Sciences, Radboudumc, Nijmegen, The Netherlands, Radboud Center for Infectious Diseases, Radboudumc, Nijmegen, The Netherlands

  • Marc J. Eleveld,

    Roles Investigation, Writing – review & editing

    Affiliations Section Pediatric Infectious Diseases, Laboratory of Medical Immunology, Radboud Institute for Molecular Life Sciences, Radboudumc, Nijmegen, The Netherlands, Radboud Center for Infectious Diseases, Radboudumc, Nijmegen, The Netherlands

  • Fred J. H. van Opzeeland,

    Roles Investigation, Writing – review & editing

    Affiliations Section Pediatric Infectious Diseases, Laboratory of Medical Immunology, Radboud Institute for Molecular Life Sciences, Radboudumc, Nijmegen, The Netherlands, Radboud Center for Infectious Diseases, Radboudumc, Nijmegen, The Netherlands

  • Saskia van Selm,

    Roles Investigation, Writing – review & editing

    Affiliations Section Pediatric Infectious Diseases, Laboratory of Medical Immunology, Radboud Institute for Molecular Life Sciences, Radboudumc, Nijmegen, The Netherlands, Radboud Center for Infectious Diseases, Radboudumc, Nijmegen, The Netherlands

  • Ronald de Groot,

    Roles Formal analysis, Writing – review & editing

    Affiliations Section Pediatric Infectious Diseases, Laboratory of Medical Immunology, Radboud Institute for Molecular Life Sciences, Radboudumc, Nijmegen, The Netherlands, Radboud Center for Infectious Diseases, Radboudumc, Nijmegen, The Netherlands

  • Marien I. de Jonge,

    Roles Formal analysis, Writing – review & editing

    Affiliations Section Pediatric Infectious Diseases, Laboratory of Medical Immunology, Radboud Institute for Molecular Life Sciences, Radboudumc, Nijmegen, The Netherlands, Radboud Center for Infectious Diseases, Radboudumc, Nijmegen, The Netherlands

  • Frits R. Mooi,

    Roles Conceptualization, Formal analysis, Writing – original draft, Writing – review & editing

    Affiliations Section Pediatric Infectious Diseases, Laboratory of Medical Immunology, Radboud Institute for Molecular Life Sciences, Radboudumc, Nijmegen, The Netherlands, Radboud Center for Infectious Diseases, Radboudumc, Nijmegen, The Netherlands, Centre for Infectious Diseases Research, Diagnostics and Screening (IDS), National Institute of Public Health and the Environment (RIVM), Bilthoven, The Netherlands, School of Biotechnology and Biomolecular Sciences, University of New South Wales Sydney, Sydney, Australia

  •  [ ... ],
  • Dimitri A. Diavatopoulos

    Roles Conceptualization, Formal analysis, Investigation, Writing – original draft, Writing – review & editing

    Dimitri.Diavatopoulos@Radboudumc.nl

    Affiliations Section Pediatric Infectious Diseases, Laboratory of Medical Immunology, Radboud Institute for Molecular Life Sciences, Radboudumc, Nijmegen, The Netherlands, Radboud Center for Infectious Diseases, Radboudumc, Nijmegen, The Netherlands

  • [ view all ]
  • [ view less ]

Abstract

Bordetella pertussis vaccine escape mutants that lack expression of the pertussis antigen pertactin (Prn) have emerged in vaccinated populations in the last 10–20 years. Additionally, clinical isolates lacking another acellular pertussis (aP) vaccine component, filamentous hemagglutinin (FHA), have been found sporadically. Here, we show that both whole-cell pertussis (wP) and aP vaccines induced protection in the lungs of mice, but that the wP vaccine was more effective in nasal clearance. Importantly, bacterial populations isolated from the lungs shifted to an FHA-negative phenotype due to frameshift mutations in the fhaB gene. Loss of FHA expression was strongly selected for in Prn-deficient strains in the lungs following aP but not wP vaccination. The combined loss of Prn and FHA led to complete abrogation of bacterial surface binding by aP-induced serum antibodies. This study demonstrates vaccine- and anatomical site-dependent adaptation of B. pertussis and has major implications for the design of improved pertussis vaccines.

Introduction

The first vaccines against pertussis–a highly contagious respiratory disease primarily caused by Bordetella pertussis—were comprised of killed whole bacteria and introduced in the 1940s-1960s [1]. Widespread implementation of diphtheria toxoid, tetanus toxoid and whole-cell pertussis (DTwP) combination vaccines significantly reduced pertussis morbidity and mortality [2]. However, DTwP vaccines also had significant side effects and showed variable vaccine effectiveness, ranging from excellent to almost completely ineffective [3, 4]. These factors led to decreased vaccine acceptance and even temporary cessation of pertussis vaccination in some countries [5]. Consequently, less reactogenic acellular pertussis combination vaccines were developed (DTaP), comprised of 1–5 purified pertussis antigens (DTaP), i.e. pertussis toxin (PT), filamentous hemagglutinin (FHA), pertactin (Prn) and fimbriae (Fim2 and Fim3). Most high-income countries switched to DTaP vaccines from the 1990s onwards, whereas the majority of low- and middle-income countries still use DTwP vaccines [6].

The recent emergence of strains that do not express Prn (Prn- strains) [712] has raised concern. Selective pressure by DTaP vaccines may have driven the expansion of Prn- strains, which reach prevalences of 78–85% in some regions [1315]. Although Prn- strains were not attenuated compared to Prn+ strains in naïve mice, Prn- strains did have a selective advantage compared to Prn+ strains in aP-vaccinated mice [16]. Interestingly, in a mixed infection model, Prn- strains outcompeted Prn+ strains in aP-vaccinated mice, which was completely reversed in naïve mice [17]. In this study, we investigated immune selective pressure induced by DTaP and DTwP vaccines on recent Prn+ and Prn- clinical isolates in mice.

Material and methods

Bacterial strains and growth conditions

Bordetella pertussis strains positive (B1865 and B1917, isolated in the Netherlands in 2000) or negative (B3621 and B3629, isolated in France in 2008 and 2009, respectively) for Prn expression were used for mouse challenge [18]. Inoculation stocks were prepared by growing B. pertussis in chemically defined THIJS medium under non-modulating conditions, as previously described [19, 20]. Bacteria were harvested at mid-low growth phase (OD620 0.5–0.6) and stored at -70°C.

B. pertussis vaccination and infection

Animal experiments were approved by the Radboudumc Committee for Animal Ethics and conducted in accordance with the relevant Dutch legislation. Naïve mice were anesthetized and challenged intranasally with the different B. pertussis strains described above. For vaccination experiments, naïve mice were immunized twice with 3-week intervals by subcutaneous injection with DTaP2, DTaP3, or DTwP. Vaccinated mice were then challenged three weeks after the final dose as described above. Nasal and lung bacterial load were determined on day 3, 7, and 14 after challenge as described previously [21]. The area under the curve (AUC) was calculated for the bacterial load for each vaccine and challenge strain, using the trapezium method [22]. The length of the fhaB G-tract of the bacterial population after infection was determined using fhaB PCR. Prn mutations were verified by PCR before and after passage through the mouse as described previously [15]. Detailed procedures can be found in the Materials and methods section of the Supporting Information.

fhaB phase variation

To screen large numbers of samples, a high throughput ligase detection reaction (LDR) was adapted to the fhaB G-tract [23]. LDR was performed on the fhaB PCR product containing the homopolymeric G-tract. For western blotting, the inoculation stocks of the four tested B. pertussis strains and three post-challenge B3629 bacterial samples were run, blotted, and incubated with polyclonal anti-FHA serum. Antibody binding to bacteria was measured by flow cytometry following incubation with pre-challenge serum from the different treatment groups. Detailed procedures can be found in the Materials and methods section of the Supporting Information.

Statistical analyses

Statistical analyses on the fold-differences in post-challenge CFUs and fhaB-G10 percentages to the inoculum stocks were performed using unpaired, two-tailed t-tests. Statistical analyses on the fold-differences in the binding assays were performed using Mann-Whitney U test. Correlation of FHA production and fhaB-G10 percentage was tested using Pearsons R test. Statistical analyses were performed using SPSS22 software (IBM, New York, United States) and EXCEL. Graphs were made using GraphPad Prism 5.03 (GraphPad Software, La Jolla, United States).

Results

Efficacy of pertussis vaccination against B. pertussis infection

To investigate vaccine efficacy against B. pertussis, groups of mice were immunized with DTaP vaccines containing either two (DTaP2; PT and FHA) or three (DTaP3; PT, FHA and Prn) pertussis components, or a DTwP vaccine containing inactivated whole B. pertussis. Unvaccinated mice were used as a control. For challenge, we performed single infections with two Prn+ (B1865 and B1917) and two Prn- (B3621 and B3629) strains that were recently isolated from pertussis patients. All challenge strains belong to the ptxP3 lineage, which has been highly prevalent since the 1990s in countries using DTaP vaccines [18, 24]. Apart from the Prn mutations, challenge strains were genetically nearly identical. Bacterial counts were determined in nose and lung at three, seven, and 14 days after intranasal challenge (Fig 1A). To rule out that there were major differences between bacteria recovered by lavage versus the complete bacterial pool in the nose and lungs, we compared lavage to homogenized tissue samples, showing a strong correlation between the two sampling methods (S1 Fig and Supporting Information; Material and methods for a detailed description). Bacterial loads from the challenge strains were plotted separately in Fig 1. Since the recovered bacterial load per strain was typically very low in vaccinated mice, especially in the lungs, for statistical analysis we pooled the two Prn+ strains, and we also pooled the two Prn- stains. No statistically significant differences were observed in colonization dynamics between the Prn+ and the Prn- pools in the nose of unvaccinated mice, suggesting that loss of Prn expression does not significantly attenuate colonization (Fig 1B). DTaP vaccination did not significantly reduce nasal colonization compared to unvaccinated mice, except on day seven when a significant reduction of Prn- strains (8.5-fold, p = 0.04) was observed in DTaP2-vaccinated mice and a significant reduction of Prn+ strains (6.1-fold, p = 0.01) in DTaP3-vaccinated mice (Fig 1B). Minimal bacterial clearance occurred in DTaP2- and DTaP3-vaccinated mice from day seven to 14, with no significant differences in bacterial load at day 14 between vaccinated and naïve mice (Prn+_DTaP2: 1.4-fold, p = 0.6; Prn+_DTaP3: 1.3-fold, p = 0.7; Prn-_DTaP2: 3.7-fold, p = 0.4; Prn-_DTaP3: 5.7-fold, p = 0.2). Although not significant, there was a trend towards reduction in nasal colonization in DTwP-vaccinated mice on day seven (Prn+ strains: 23.5-fold, p = 0.06; Prn- strains: 3.4-fold, p = 0.06), resulting in almost complete clearance of both Prn+ and Prn- strains by day 14. Overall, these data suggest that none of the pertussis vaccines completely prevented infection of the upper respiratory tract, although there was a trend towards enhanced clearance in DTwP-vaccinated mice, particularly towards the later stage of infection.

thumbnail
Fig 1. Colonization in the nose and lungs of naive and vaccinated mice.

(a) Design of the study. Mice were vaccinated at day 0 (d0) and at day 21 (d21) by subcutaneous (s.c.) injection with DTaP2, DTaP3, DTwP, or mock. Three weeks after the final vaccination (d42), mice were infected intranasally (i.n.) with one of the Prn+ or Prn- strains. Three (d45), seven (d49), and 14 (d56) days after challenge, bacterial load was determined in nose (b) and lung lavage (c). For each treatment group and challenge strain, the bacterial load over time is shown as the median log10 CFU ± interquartile range (n = 5–18 mice per group and time point). Dashed line indicates lower limit of detection.

https://doi.org/10.1371/journal.pone.0237394.g001

Since severe pertussis is associated with progression of B. pertussis to the lungs, we determined lung bacterial load. Infection dynamics in the lungs were similar between the Prn+ and Prn- strains in unvaccinated mice, except for day three when Prn- strains reached slightly higher bacterial loads than Prn+ strains (2.9-fold, p = 0.002). This small but significant difference in mean bacterial load was mainly attributable to strain B3629 (Fig 1C). All vaccines used in this study induced significant clearance of the B. pertussis strains in the lungs at day three compared to unvaccinated controls (Fig 1C, Prn+_DTaP2: 22.7-fold, p = 0.00005; Prn+_DTaP3: 898.1-fold, p<0.00005; Prn+_DTwP: 559.4-fold, p<0.00005; Prn-_DTaP2: 115.1-fold, p = 0.0007; Prn-_DTaP3: 74.8-fold, p = <0.00005; Prn-_DTwP: 229.2-fold, p<0.00005), demonstrating that protection against lung infection may be achieved regardless of vaccine type and composition. Although all vaccines induced significant protection against infection in the lungs compared to naïve mice, several differences were observed between vaccines. DTaP3, which contains Prn, protected significantly better against Prn+ strains than DTaP2 early after challenge (39.5-fold, p = 0.004 at day three, Fig 1C). This was also observed for wP vaccination, which was significantly more protective against Prn+ strains than DTaP2 early after challenge (24.6-fold at day three, p = 0.04). Conversely, no significant differences were observed in vaccine efficacy against Prn- strains between vaccines containing (DTwP and DTaP3) or lacking (DTaP2) Prn (Fig 1C). At day seven and 14, all vaccines significantly reduced the bacterial load of the Prn+ strains and Prn- strains compared to unvaccinated controls (Fig 1C, Day 7: Prn+_DTaP2: 523.6-fold, p = 0.008; Prn+_DTaP3: 10703-fold, p = 0.0002; Prn+_DTwP: 500.9-fold, p = 0.01; Prn-_DTaP2: 299.1-fold, p = 0.007; Prn-_DTaP3: 84.7-fold, p = <0.00005; Prn-_DTwP: 118.9-fold, p = 0.001. Day 14: Prn+_DTaP2: 23.0-fold, p = 0.03; Prn+_DTaP3: 93.6-fold, p = 0.02; Prn+_DTwP: 150.2-fold, p = 0.02; Prn-_DTaP2: 163.8-fold, p = 0.05; Prn-_DTaP3: 95.3-fold, p = 0.01; Prn-_DTwP: 4.6-fold, p = 0.04). At day seven, DTaP3 vaccination was still more effective in clearing Prn+ strains compared to DTaP2 vaccination (20.4-fold, p = 0.0004). Of note, even though DTwP-vaccinated mice showed nearly sterilizing immunity in the lungs against all strains, a slightly better survival of one of the Prn- strains (B3621) was observed early after challenge compared to Prn+ strains on day 3 (12.5-fold, p = 0.0002, Fig 1C). We also calculated the area under the curve (AUC) for each vaccine and challenge strain in the nose and lungs, which is an inverse measure of total vaccine efficacy, confirming the results described above (Fig 1B and 1C).

To observe the differences in vaccine efficacy in the lungs between Prn+ and Prn- strains after DTaP2 and DTaP3 vaccination, a re-organized version of Fig 1B and 1C is included as a Supporting file where each panel represents the bacterial load over time per strain for the different treatment groups (S2 Fig).

DTaP vaccination is associated with phase-variable loss of FHA expression in the lungs in Prn- strains

Whole genome sequencing of fhaB in 49 clinical isolates identified phase variation in its homopolymeric G-tract. The wild type fhaB allele contains 10 Gs (fhaB-G10), while the two major mutant alleles contain 9 Gs (fhaB-G9) and 11 Gs (fhaB-G11) (Fig 2A). Notably, all 49 strains contained a mixture of fhaB-G10 with either fhaB-G9 or fhaB-G11. We speculate that variation in the fhaB G-tract would affect FHA expression. Indeed, we found a positive relationship between FHA levels and the wild type to mutant fhaB allele ratio in the 49 clinical strains (S3 Fig).

thumbnail
Fig 2. Changes in FHA expression in B. pertussis after passage through nose and lungs of naive and vaccinated mice.

(a) DNA sequence of the three homopolymeric fhaB-G tract variants observed in B. pertussis. The wild type fhaB gene is 10,773 nucleotides (nt) long, with a homopolymeric tract of 10 Gs at nucleotide position 1,078–1,088. Frameshift mutations in this G-tract result in the introduction of stop codons at nucleotides 1,541 and 1,090 for fhaB-G11 and fhaB-G9, respectively. (b-c) Bacteria from the nose (b) or lungs (c) were pooled per mouse. Each symbol represents the log-transformed fhaB-G10 percentage of the recovered bacteria normalized to the fhaB-G10 percentage of its respective inoculum. Data is shown for all time points combined. Horizontal lines represent median ± interquartile range. Dashed line indicates no change in fhaB-G10 percentage of the recovered bacteria compared to its inoculum. (d) Expression of FHA in inoculum batches and in post-challenge bacterial samples was determined by western blot. Polyclonal antibodies to FHA were used; purified FHA and a fhaB knockout strain (ΔfhaB) were included as a positive and negative control, respectively. B3629-m089, B3629-m086 and B3629-m093 contain predominantly wild type, mutant, and mutant fhaB alleles, respectively. The position of molecular weight marker (in kDa) is shown.

https://doi.org/10.1371/journal.pone.0237394.g002

To assess the role of Prn and FHA in immunity against Prn+ and Prn- strains, we studied the impact of vaccination on bacterial populations infecting the nose or lungs of mice by quantifying shifts in fhaB allele frequencies and bacterial loads.

For each strain we normalized the fhaB-G10 percentage of the B. pertussis population recovered from the challenged mice to that of its respective inoculum, where a value of 0 indicates a similar composition of the recovered bacteria and the inoculum. Conversely, a value higher than 0 indicates a proportional shift towards mutant fhaB alleles. The effect of vaccination was then assessed. Since the recovered bacterial load in the lungs was typically very low in vaccinated mice, results from the two Prn+ strains were pooled and results from the two Prn- strains were also pooled for statistical analyses. In unvaccinated mice, no shifts in fhaB-G10 percentage were detected in the nose. A small but significant shift towards mutant fhaB alleles was observed in the nose of DTaP3-vaccinated mice for both Prn+ and Prn- strains, compared to unvaccinated controls (Fig 2B, 1.6-fold, p = 0.046 and 1.6-fold, p = 0.019, respectively). This difference was mainly attributable to strain B3629, which showed a significant increase in mutant fhaB alleles compared to strain B3621. However, overall the fhaB allele composition of the recovered nasal bacterial populations was very similar to the inocula. Conversely, a strikingly different pattern emerged in the lungs. All B. pertussis populations isolated from the lungs of mice vaccinated with DTaP2 or DTaP3 showed a significant shift towards mutant fhaB alleles compared to unvaccinated mice (Fig 2C, Prn+ strains: 3.9-fold, p = 0.001 and 4.4-fold, p = 0.00005 for DTaP2 and DTaP3, respectively; Prn- strains: 2.9-fold, p = <0.00001 and 3.1-fold, p = <0.00001 for DTaP2 and DTaP3, respectively). A small but statistically significant shift towards mutant fhaB alleles was also observed for DTwP compared to unvaccinated controls, but only for Prn+ strains (2.9-fold, p = 0.03). Intriguingly, overall increases in mutant allele percentages in the recovered bacterial populations were much more pronounced in strains lacking Prn than in Prn+ strains. Nonetheless, shifts towards mutant fhaB alleles were also observed in Prn+ strains in mice vaccinated with DTaP vaccines. When comparing the pooled Prn- strains with the pooled Prn+ strains, the former had shifted significantly more towards mutant fhaB alleles in naïve, DTaP2-, and DTaP3-vaccinated mice (3.1-fold, p = 0.00001; 2.3-fold, p = 0.004 and 2.2-fold, p = 0.002, respectively), suggesting that Prn+ and Prn- strains may differ with regards to the fitness of mutant fhaB alleles. Together, these results suggest that DTaP but not DTwP vaccines induce strong immune selective pressure on FHA expression in the lungs but not the nose, which is most evident in a Prn- background.

To determine whether the observed shift in fhaB genotype correlated with FHA protein expression, western blot was performed on the recovered bacterial samples that contained predominately wild type (B3629-m089 containing 63% fhaB-G10) or mutant (B3629-m086 containing 37% fhaB-G10 and B3629-m093 containing 16% fhaB-G10) fhaB alleles. FHA expression in these samples was compared to the inoculum stocks. Bacterial lysates were incubated with polyclonal heat-inactivated serum from mice vaccinated with purified FHA (Fig 2D and S4 Fig). All inoculum stocks expressed FHA at the time of infection. FHA expression was detected in B3629-m089, but not in samples predominantly containing mutant fhaB alleles.

Loss of recognition by DTaP-induced antibodies

To determine the impact of FHA and Prn non-expression on immunological recognition, antibody opsonization experiments were performed. Bacteria were opsonized with heat-inactivated post-vaccination, pre-challenge sera from vaccinated mice, after which IgG and IgM binding was determined by flow cytometry (Fig 3A and 3B). DTaP2 vaccination induced a moderate increase in IgG binding to B1917 (Prn+ FHA+) compared to unvaccinated controls (10.6-fold increase, p = 0.0022). Vaccination with DTaP3 and DTwP further increased IgG binding (53.1-fold, p = 0.0022 and 55.2-fold, p = 0.0022, respectively). Prn non-expression reduced binding by DTaP3-induced IgG to a level equivalent to DTaP2 (11.5-fold compared to the Prn+ strain, p = 0.0087). Prn-deficiency also affected DTwP-induced IgG opsonization levels (3.6-fold lower binding than to the Prn+ strain, p = 0.0022). The subsequent loss of FHA expression resulted in a small 1.5-fold reduction in IgG opsonization by DTwP compared to the Prn-deficient strain (p = 0.065). Strikingly, the combined loss of Prn and FHA expression almost completely abrogated IgG binding by either DTaP2 or DTaP3, with opsonization levels comparable to unvaccinated mice. IgM binding patterns were essentially similar to IgG (Fig 3B). These results imply a selective advantage for Prn- FHA- variants in vaccinated mice and possibly also humans.

thumbnail
Fig 3. Antibody opsonization of B. pertussis.

Binding of IgG (a) and IgM (b) antibodies to B1917 (Prn+ FHA+), B3629-m089 (Prn- FHA+), and B3629-m093 (Prn- FHA-) using 10% mouse sera from the different vaccination groups. Binding was measured using flow cytometry, MFI indicates the mean fluorescence intensity. Horizontal lines represent median. * p≤0.05, ** p≤0.01, ** p≤0.001.

https://doi.org/10.1371/journal.pone.0237394.g003

Discussion

The causes of the resurgence of pertussis have been much debated, including the relationship between pertussis vaccination and temporal changes in the B. pertussis population [25, 26]. A recent phenomenon has been the emergence and expansion of strains that have lost expression of Prn [79, 14, 15]. These strains have increased in frequency in countries that have implemented DTaP vaccines but not in countries using DTwP vaccines. This difference is very likely caused by the higher Prn antibody levels induced by DTaP vaccines compared to DTwP vaccines [27].Our findings demonstrate that different pertussis vaccines may exert distinct immune selective pressure on B. pertussis in the nose and lungs. We showed that B. pertussis is able to colonize the upper and lower respiratory tract. Both biotic (microbiome, receptors, immune responses) and abiotic (temperature, O2 pressure) factors differ at these anatomical sites. We propose that phase variation of FHA (and possibly other genes) plays a role in the adaptation of B. pertussis to the different environments it encounters in the host.

Using a mouse model, which shares many similarities to humans with regards to immunological mechanisms of B. pertussis [28], we found essentially no attenuation of Prn- strains in the absence of vaccine-induced immunity, consistent with previous literature [16]. In the nose, we observed a trend towards enhanced clearance in DTwP-vaccinated mice compared to mice vaccinated with DTaP, particularly towards the late stage of infection, consistent with recent findings [29, 30].

In the lungs, all pertussis vaccines induced early clearance of B. pertussis. As expected, DTaP vaccines containing Prn were more efficacious against Prn+ strains than against Prn- strains in the lungs early after challenge. These findings are in line with previous literature and support the hypothesis that Prn- strains have emerged because they are more fit in individuals vaccinated with Prn-containing DTaP vaccines [8, 10, 15, 16, 3032]. Indeed, pertussis patients receiving at least one pertussis vaccine dose were more likely to be infected by a Prn- strain than unvaccinated patients [14]. Of note, a recent study using pertussis disease as the clinical endpoint showed that vaccine effectiveness of DTaP vaccines against Prn- strains remains high [33]. An important difference with this study is that we do not have a clinical endpoint of disease due to the use of mice.

We provide evidence that FHA phase variation allows B. pertussis populations to adapt to different anatomical sites in the mouse as reflected by changes in the fhaB allele frequencies. Specific selection of mutant fhaB alleles was most prominent in the lungs of DTaP-vaccinated mice, as evident from the stronger proportional shift of Prn- than Prn+ bacteria towards mutant fhaB alleles. The primary loss of Prn expression provides a selective advantage in the lungs of DTaP-vaccinated hosts, reflected in the 115-fold and 75-fold difference in bacterial load early after infection in the lungs after DTaP2 and DTaP3 vaccination, respectively. Consequently, for these strains a larger bacterial pool may be available from which bacteria bearing mutant fhaB alleles can be selected. FHA phase variants were only rarely detected in the nose, suggesting that FHA is essential during early colonization of the nose but not the lungs, in line with previous findings [3437]. An alternative, not mutually exclusive explanation is that DTaP-induced immunity is simply less effective in the nose than in the lungs and therefore exerts less selective pressure in the nose. Previous human vaccination studies have also raised questions about the contribution of FHA to protection [38, 39].

Besides differences in anatomy there are also other differences in the immune system between mice and humans and consequently care should be taken to extrapolate our findings directly to humans. However, the general immunological mechanisms by which vaccines protect against pertussis have been shown to be quite similar between humans and mice [40] and studies in mice have provided important insights into immunity to B. pertussis infection. We show here that none of the vaccines completely prevent infection in the nose, but all the vaccines induced significant protection against infection in the lungs compared to naïve mice. We need further studies to confirm these findings in humans, but from our experience with the controlled human pertussis infection model [41] human volunteers who were vaccinated during infancy and who were challenged intranasally with B. pertussis can be readily colonized. At present it is unclear what the role is of FHA phase variation in humans. This study may give incentive to analyze fhaB phase variation in human nasopharynx and lung samples.

Using western blot, we only detected FHA expression in recovered bacterial populations that contained predominantly wild type fhaB alleles (B3629-m089) and not in the samples containing predominantly mutant fhaB alleles (B3629-m086 and B3629-m093). Although bacteria bearing mutant alleles were predicted to produce small truncated FHA molecules due to premature termination of translation, immunoblotting with anti-FHA polyclonal serum did not detect such protein fragments. Possibly, the small fragments are not produced, are degraded, are not recognized by the serum used, or the detection method is not sensitive enough to detect FHA fragments in mixed populations containing lower proportions of wild type fhaB alleles. It is conceivable that the FHA fragments are present, but not recognized by the antibodies raised against the native protein. Interestingly, both predicted fragments contain a two-partner secretion (TPS) domain which is essential for secretion [42]. Indeed, a 304 residue N-terminal FHA fragment, nearly identical to the 292 residue fragment predicted to be produced by fhaB-G9, was shown to be secreted by B. pertussis [42]. Thus, it is possible that the fhaB-G9 and/or fhaB-G11 fragments are excreted and have some function. However, apart from the TPS domain, known functional domains of FHA have been located outside the predicted fhaB-G9 and/or fhaB-G11 fragments [43]. Strikingly, B. pertussis bacteria that lost expression of both Prn and FHA were not recognized anymore by DTaP-induced antibodies. With the loss of both Prn and FHA as antigenic targets for opsonization, DTaP-induced immunity against B. pertussis will therefore be mostly dependent on the recognition of PT. Although PT is predominantly a secreted protein and anti-PT antibodies do not bind to the bacterium, monocomponent pertussis vaccines based only on PT have been shown to be effective in Denmark [44]. Because Prn- FHA- variants were ultimately still cleared in the lungs of DTaP-vaccinated animals, the question is thus how protection is achieved. It seems likely that PT antibodies confer passive protection in mice, presumably through neutralization of the biological activity of PT [45]. In addition, cellular immunity may play a role.

Although naturally occurring FHA mutants have been described before [46] and we also identified clinical strains with mutations in the fhaB-G10 region, they have not been identified very frequently. One possible explanation for this is that the standard diagnostic sampling procedure for pertussis patients is to collect a nasopharyngeal swab. Because such samples are not obtained from the deeper airways, it is entirely possible that FHA phase variation in human pertussis patients is a frequent event that remains largely undetected.

An important question that we have not yet addressed in this study is the relation between bacterial load, disease induction and transmission. Since the Prn- FHA- variants were not recognized by DTaP-induced immunity in the lungs, it is likely that non-expression of both Prn and FHA as antigenic targets has a direct effect on survival of the bacteria and hence the duration, disease severity, and symptoms. The selection of FHA-negative phase variants in the lungs and reduced immune recognition of these variants in DTaP-vaccinated hosts could increase the duration of infection, the bacterial load and the severity of clinical symptoms, thereby potentially enhancing transmission of B. pertussis. Although Warfel et al. showed that transmission may occur even in the absence of clinical symptoms [30], it is not unlikely that infection of the lungs and the induction of coughing enhances transmission rates. This would imply that there is positive selective pressure on B. pertussis for infection of the lower respiratory tract of vaccinated individuals and the induction of clinical symptoms. Transmission is thought to occur during the catarrhal phase until up to three weeks after the start of the paroxysmal phase. Since FHA phase variation was virtually undetected in the upper respiratory tract but was commonly detected in the lower respiratory tract, it would be interesting to investigate the transmission and disease capacity of strains deficient in both Prn and FHA and the role of phase-variable antigen expression in the epidemiology of pertussis, for instance in the baboon model.

In conclusion, our study highlights the importance of studying the effects of vaccination on pathogen populations in vivo. Importantly, it suggests that the efficacy of DTaP vaccines against Prn- strains may be further reduced due to phase variable expression of FHA. Immunity to two acellular pertussis vaccine components may therefore be compromised.

Supporting information

S1 Fig. Correlation between bacterial populations recovered from the airways via lavage versus homogenized tissue samples from the same mouse.

(a) Correlation between %fhaB-G10 of lung lavage samples and homogenized lung tissue. (b) Correlation between %fhaB-G10 of nose lavage samples and homogenized nasopharynx tissue. Dashed lines indicate 95% confidence intervals.

https://doi.org/10.1371/journal.pone.0237394.s002

(EPS)

S2 Fig. Colonization in the nose and lungs of naive and vaccinated mice.

Mice were vaccinated at day 0 (d0) and at day 21 (d21) by subcutaneous (s.c.) injection with DTaP2, DTaP3, DTwP, or mock. Three weeks after the final vaccination (d42), mice were infected intranasally (i.n.) with one of the Prn+ or Prn- strains. Three (d45), seven (d49), and 14 (d56) days after challenge, bacterial load was determined in nose (a) and lung lavage (b). Each panel represents the bacterial load over time per strain for the different treatment groups, shown as the median log10 CFU ± interquartile range (n = 5–18 mice per group and time point). Dashed line indicates lower limit of detection.

https://doi.org/10.1371/journal.pone.0237394.s003

(EPS)

S3 Fig. Correlation of FHA production with the percentage of fhaB-G10.

The amount of FHA produced and the length of the homopolymeric G-tract was determined using 49 clinical isolates. The percentage G10 of each isolate was calculated by dividing the G10 (wild type) LDR signal by the total G10, G9, and G11 signal. A significant positive correlation was found (Pearson r 0.49, p<0.05).

https://doi.org/10.1371/journal.pone.0237394.s004

(EPS)

S4 Fig. Unadjusted western blot image.

Original blot containing molecular weight marker and loading order of experimental samples.

https://doi.org/10.1371/journal.pone.0237394.s005

(EPS)

Acknowledgments

We thank E. Gogol and C. Cummings for their valuable advice on the execution of the ligase detection reaction, A. Lutyńska from the Polish National Institute of Public Health–National Institute of Hygiene, for providing the Polish DTwP vaccine and A. Gorringe from Public Health England for providing the anti-FHA polyclonal antibodies.

References

  1. 1. Mills KH, Ross PJ, Allen AC, Wilk MM. Do we need a new vaccine to control the re-emergence of pertussis? Trends Microbiol. 2014;22(2):49–52. pmid:24485284
  2. 2. Plotkin SA. The effectiveness of whole-cell pertussis vaccines. Dev Biol Stand. 1997;89:171–4. pmid:9272348
  3. 3. Cody CL, Baraff LJ, Cherry JD, Marcy SM, Manclark CR. Nature and rates of adverse reactions associated with DTP and DT immunizations in infants and children. Pediatrics. 1981;68(5):650–60. pmid:7031583
  4. 4. Jefferson T, Rudin M, DiPietrantonj C. Systematic review of the effects of pertussis vaccines in children. Vaccine. 2003;21(17–18):2003–14. pmid:12706690
  5. 5. Edwards KM. Unraveling the challenges of pertussis. Proc Natl Acad Sci U S A. 2014;111(2):575–6. pmid:24381159
  6. 6. Lambert LC. Pertussis vaccine trials in the 1990s. J Infect Dis. 2014;209 Suppl 1:S4–9.
  7. 7. Barkoff AM, Mertsola J, Guillot S, Guiso N, Berbers G, He Q. Appearance of Bordetella pertussis strains not expressing the vaccine antigen pertactin in Finland. Clin Vaccine Immunol. 2012;19(10):1703–4. pmid:22914363
  8. 8. Bouchez V, Brun D, Cantinelli T, Dore G, Njamkepo E, Guiso N. First report and detailed characterization of B. pertussis isolates not expressing Pertussis Toxin or Pertactin. Vaccine. 2009;27(43):6034–41. pmid:19666155
  9. 9. Hegerle N, Paris AS, Brun D, Dore G, Njamkepo E, Guillot S, et al. Evolution of French Bordetella pertussis and Bordetella parapertussis isolates: increase of Bordetellae not expressing pertactin. Clin Microbiol Infect. 2012;18(9):E340–6. pmid:22717007
  10. 10. Otsuka N, Han HJ, Toyoizumi-Ajisaka H, Nakamura Y, Arakawa Y, Shibayama K, et al. Prevalence and genetic characterization of pertactin-deficient Bordetella pertussis in Japan. PLoS One. 2012;7(2):e31985. pmid:22348138
  11. 11. Pawloski LC, Queenan AM, Cassiday PK, Lynch AS, Harrison MJ, Shang W, et al. Prevalence and molecular characterization of pertactin-deficient Bordetella pertussis in the United States. Clin Vaccine Immunol. 2014;21(2):119–25. pmid:24256623
  12. 12. Tsang RS, Shuel M, Jamieson FB, Drews S, Hoang L, Horsman G, et al. Pertactin-negative Bordetella pertussis strains in Canada: characterization of a dozen isolates based on a survey of 224 samples collected in different parts of the country over the last 20 years. Int J Infect Dis. 2014;28:65–9. pmid:25244999
  13. 13. Lam C, Octavia S, Ricafort L, Sintchenko V, Gilbert GL, Wood N, et al. Rapid increase in pertactin-deficient Bordetella pertussis isolates, Australia. Emerg Infect Dis. 2014;20(4):626–33. pmid:24655754
  14. 14. Martin SW, Pawloski L, Williams M, Weening K, DeBolt C, Qin X, et al. Pertactin-negative Bordetella pertussis strains: evidence for a possible selective advantage. Clin Infect Dis. 2015;60(2):223–7. pmid:25301209
  15. 15. Zeddeman A, van Gent M, Heuvelman CJ, van der Heide HG, Bart MJ, Advani A, et al. Investigations into the emergence of pertactin-deficient Bordetella pertussis isolates in six European countries, 1996 to 2012. Euro Surveill. 2014;19(33).
  16. 16. Hegerle N, Dore G, Guiso N. Pertactin deficient Bordetella pertussis present a better fitness in mice immunized with an acellular pertussis vaccine. Vaccine. 2014;32(49):6597–600. pmid:25312274
  17. 17. Safarchi A, Octavia S, Luu LD, Tay CY, Sintchenko V, Wood N, et al. Pertactin negative Bordetella pertussis demonstrates higher fitness under vaccine selection pressure in a mixed infection model. Vaccine. 2015;33(46):6277–81. pmid:26432908
  18. 18. Bart MJ, van der Heide HG, Zeddeman A, Heuvelman K, van Gent M, Mooi FR. Complete Genome Sequences of 11 Bordetella pertussis Strains Representing the Pandemic ptxP3 Lineage. Genome Announc. 2015;3(6).
  19. 19. Thalen M, van den IJ, Jiskoot W, Zomer B, Roholl P, de Gooijer C, et al. Rational medium design for Bordetella pertussis: basic metabolism. J Biotechnol. 1999;75(2–3):147–59. pmid:10553654
  20. 20. de Gouw D, Hermans PW, Bootsma HJ, Zomer A, Heuvelman K, Diavatopoulos DA, et al. Differentially expressed genes in Bordetella pertussis strains belonging to a lineage which recently spread globally. PLoS One. 2014;9(1):e84523. pmid:24416242
  21. 21. Cron LE, Stol K, Burghout P, van Selm S, Simonetti ER, Bootsma HJ, et al. Two DHH subfamily 1 proteins contribute to pneumococcal virulence and confer protection against pneumococcal disease. Infect Immun. 2011;79(9):3697–710. pmid:21768284
  22. 22. Baggot JD. Some aspects of clinical pharmacokinetics in veterinary medicine II. Journal of Veterinary Pharmacology and Therapeutics. 1978;1(2):111–8.
  23. 23. Gogol EB, Cummings CA, Burns RC, Relman DA. Phase variation and microevolution at homopolymeric tracts in Bordetella pertussis. BMC Genomics. 2007;8:122. pmid:17509142
  24. 24. Mooi FR, van Loo IH, van Gent M, He Q, Bart MJ, Heuvelman KJ, et al. Bordetella pertussis strains with increased toxin production associated with pertussis resurgence. Emerg Infect Dis. 2009;15(8):1206–13. pmid:19751581
  25. 25. Mooi FR, van Loo IH, King AJ. Adaptation of Bordetella pertussis to vaccination: a cause for its reemergence? Emerg Infect Dis. 2001;7(3 Suppl):526–8.
  26. 26. Mooi FR, Van Der Maas NA, De Melker HE. Pertussis resurgence: waning immunity and pathogen adaptation—two sides of the same coin. Epidemiol Infect. 2014;142(4):685–94. pmid:23406868
  27. 27. Edwards KM, Meade BD, Decker MD, Reed GF, Rennels MB, Steinhoff MC, et al. Comparison of 13 acellular pertussis vaccines: overview and serologic response. Pediatrics. 1995;96(3 Pt 2):548–57.
  28. 28. Mills KH. Immunity to Bordetella pertussis. Microbes Infect. 2001;3(8):655–77. pmid:11445452
  29. 29. Wilk MM, Borkner L, Misiak A, Curham L, Allen AC, Mills KHG. Immunization with whole cell but not acellular pertussis vaccines primes CD4 TRM cells that sustain protective immunity against nasal colonization with Bordetella pertussis. Emerg Microbes Infect. 2019;8(1):169–85. pmid:30866771
  30. 30. Warfel JM, Zimmerman LI, Merkel TJ. Acellular pertussis vaccines protect against disease but fail to prevent infection and transmission in a nonhuman primate model. Proc Natl Acad Sci U S A. 2014;111(2):787–92. pmid:24277828
  31. 31. van Gent M, van Loo IH, Heuvelman KJ, de Neeling AJ, Teunis P, Mooi FR. Studies on prn variation in the mouse model and comparison with epidemiological data. PLoS One. 2011;6(3):e18014. pmid:21464955
  32. 32. de Gouw D, de Jonge MI, Hermans PW, Wessels HJ, Zomer A, Berends A, et al. Proteomics-identified Bvg-activated autotransporters protect against bordetella pertussis in a mouse model. PLoS One. 2014;9(8):e105011. pmid:25133400
  33. 33. Breakwell L, Kelso P, Finley C, Schoenfeld S, Goode B, Misegades LK, et al. Pertussis Vaccine Effectiveness in the Setting of Pertactin-Deficient Pertussis. Pediatrics. 2016;137(5).
  34. 34. Kimura A, Mountzouros KT, Relman DA, Falkow S, Cowell JL. Bordetella pertussis filamentous hemagglutinin: evaluation as a protective antigen and colonization factor in a mouse respiratory infection model. Infect Immun. 1990;58(1):7–16. pmid:2294058
  35. 35. Geuijen CA, Willems RJ, Bongaerts M, Top J, Gielen H, Mooi FR. Role of the Bordetella pertussis minor fimbrial subunit, FimD, in colonization of the mouse respiratory tract. Infect Immun. 1997;65(10):4222–8. pmid:9317030
  36. 36. Mooi FR, Jansen WH, Brunings H, Gielen H, van der Heide HG, Walvoort HC, et al. Construction and analysis of Bordetella pertussis mutants defective in the production of fimbriae. Microb Pathog. 1992;12(2):127–35. pmid:1350044
  37. 37. Julio SM, Inatsuka CS, Mazar J, Dieterich C, Relman DA, Cotter PA. Natural-host animal models indicate functional interchangeability between the filamentous haemagglutinins of Bordetella pertussis and Bordetella bronchiseptica and reveal a role for the mature C-terminal domain, but not the RGD motif, during infection. Mol Microbiol. 2009;71(6):1574–90. pmid:19220744
  38. 38. Cherry JD, Gornbein J, Heininger U, Stehr K. A search for serologic correlates of immunity to Bordetella pertussis cough illnesses. Vaccine. 1998;16(20):1901–6. pmid:9796041
  39. 39. Olin P, Hallander HO, Gustafsson L, Reizenstein E, Storsaeter J. How to make sense of pertussis immunogenicity data. Clin Infect Dis. 2001;33 Suppl 4:S288–91.
  40. 40. Mills KH, Ryan M, Ryan E, Mahon BP. A murine model in which protection correlates with pertussis vaccine efficacy in children reveals complementary roles for humoral and cell-mediated immunity in protection against Bordetella pertussis. Infect Immun. 1998;66(2):594–602. pmid:9453614
  41. 41. de Graaf H, Ibrahim M, Hill AR, Gbesemete D, Vaughan AT, Gorringe A, et al. Controlled human infection with Bordetella pertussis induces asymptomatic, immunising colonisation. Clin Infect Dis. 2019.
  42. 42. Clantin B, Hodak H, Willery E, Locht C, Jacob-Dubuisson F, Villeret V. The crystal structure of filamentous hemagglutinin secretion domain and its implications for the two-partner secretion pathway. Proc Natl Acad Sci U S A. 2004;101(16):6194–9. pmid:15079085
  43. 43. Villarino Romero R, Osicka R, Sebo P. Filamentous hemagglutinin of Bordetella pertussis: a key adhesin with immunomodulatory properties? Future Microbiol. 2014;9(12):1339–60. pmid:25517899
  44. 44. Thierry-Carstensen B, Dalby T, Stevner MA, Robbins JB, Schneerson R, Trollfors B. Experience with monocomponent acellular pertussis combination vaccines for infants, children, adolescents and adults—a review of safety, immunogenicity, efficacy and effectiveness studies and 15 years of field experience. Vaccine. 2013;31(45):5178–91. pmid:23994021
  45. 45. Sato H, Ito A, Chiba J, Sato Y. Monoclonal antibody against pertussis toxin: effect on toxin activity and pertussis infections. Infect Immun. 1984;46(2):422–8. pmid:6094351
  46. 46. Weigand MR, Pawloski LC, Peng Y, Ju H, Burroughs M, Cassiday PK, et al. Screening and Genomic Characterization of Filamentous Hemagglutinin-Deficient Bordetella pertussis. Infect Immun. 2018;86(4).