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.

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

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.

Bordetella pertussis strains with increased toxin production associated with pertussis resurgence.

Before childhood vaccination was introduced in the 1940s, pertussis was a major cause of infant death worldwide. Widespread vaccination of children succeeded in reducing illness and death. In the 1990s, a resurgence of pertussis was observed in a number of countries with highly vaccinated populations, and pertussis has become the most prevalent vaccine-preventable disease in industrialized countries. We present evidence that in the Netherlands the dramatic increase in pertussis is temporally associated with the emergence of Bordetella pertussis strains carrying a novel allele for the pertussis toxin promoter, which confers increased pertussis toxin (Ptx) production. Epidemiologic data suggest that these strains are more virulent in humans. We discuss changes in the ecology of B. pertussis that may have driven this adaptation. Our results underline the importance of Ptx in transmission, suggest that vaccination may select for increased virulence, and indicate ways to control pertussis more effectively.

Small mutations in Bordetella pertussis are associated with selective sweeps.

Bordetella pertussis is the causative agent of pertussis, a highly contagious disease of the human respiratory tract. Despite high vaccination coverage, pertussis has resurged and has become one of the most prevalent vaccine-preventable diseases in developed countries. We have proposed that both waning immunity and pathogen adaptation have contributed to the persistence and resurgence of pertussis. Allelic variation has been found in virulence-associated genes coding for the pertussis toxin A subunit (ptxA), pertactin (prn), serotype 2 fimbriae (fim2), serotype 3 fimbriae (fim3) and the promoter for pertussis toxin (ptxP). In this study, we investigated how more than 60 years of vaccination has affected the Dutch B. pertussis population by combining data from phylogeny, genomics and temporal trends in strain frequencies. Our main focus was on the ptxA, prn, fim3 and ptxP genes. However, we also compared the genomes of 11 Dutch strains belonging to successful lineages. Our results showed that, between 1949 and 2010, the Dutch B. pertussis population has undergone as least four selective sweeps that were associated with small mutations in ptxA, prn, fim3 and ptxP. Phylogenetic analysis revealed a stepwise adaptation in which mutations accumulated clonally. Genomic analysis revealed a number of additional mutations which may have a contributed to the selective sweeps. Five large deletions were identified which were fixed in the pathogen population. However, only one was linked to a selective sweep. No evidence was found for a role of gene acquisition in pathogen adaptation. Our results suggest that the B. pertussis gene repertoire is already well adapted to its current niche and required only fine tuning to persist in the face of vaccination. Further, this work shows that small mutations, even single SNPs, can drive large changes in the populations of bacterial pathogens within a time span of six to 19 years.

Studies on Prn variation in the mouse model and comparison with epidemiological data.

The virulence factor pertactin (Prn) is a component of pertussis vaccines and one of the most polymorphic Bordetella pertussis antigens. After the introduction of vaccination shifts in predominant Prn types were observed and strains with the Prn vaccine type (Prn1) were replaced by strains carrying non-vaccine types (Prn2 and Prn3), suggesting vaccine-driven selection. The aim of this study was to elucidate the shifts observed in Prn variants. We show that, although Prn2 and Prn3 circulated in similar frequencies in the 1970s and 1980s, in the 1990s Prn2 strains expanded and Prn3 strains disappeared, suggesting that in vaccinated populations Prn2 strains are fitter than Prn3 strains. We established a role for Prn in the mouse model by showing that a Prn knock-out (Prn-ko) mutation reduced colonization in trachea and lungs. Restoration of the mutation resulted in a significant increase in colonization compared to the knock-out mutant. The ability of clinical isolates with different Prn variants to colonize the mouse lung was compared. Although these isolates were also polymorphic at other loci, only variation in the promoter for pertussis toxin (ptxP) and Prn were found to contribute significantly to differences in colonization. Analysis of a subset of strains with the same ptxP allele revealed that the ability to colonize mice decreased in the order Prn1>Prn2 and Prn3. Our results are consistent with the predominance of Prn1 strains in unvaccinated populations. Our results show that ability to colonize mice is practically the same for Prn2 and Prn3. Therefore other factors may have contributed to the predominance of Prn2 in vaccinated populations. The mouse model may be useful to assess and predict changes in the B. pertussis population due to vaccination.

SNP-based typing: a useful tool to study Bordetella pertussis populations.

To monitor changes in Bordetella pertussis populations, mainly two typing methods are used; Pulsed-Field Gel Electrophoresis (PFGE) and Multiple-Locus Variable-Number Tandem Repeat Analysis (MLVA). In this study, a single nucleotide polymorphism (SNP) typing method, based on 87 SNPs, was developed and compared with PFGE and MLVA. The discriminatory indices of SNP typing, PFGE and MLVA were found to be 0.85, 0.95 and 0.83, respectively. Phylogenetic analysis, using SNP typing as Gold Standard, revealed false homoplasies in the PFGE and MLVA trees. Further, in contrast to the SNP-based tree, the PFGE- and MLVA-based trees did not reveal a positive correlation between root-to-tip distance and the isolation year of strains. Thus PFGE and MLVA do not allow an estimation of the relative age of the selected strains. In conclusion, SNP typing was found to be phylogenetically more informative than PFGE and more discriminative than MLVA. Further, in contrast to PFGE, it is readily standardized allowing interlaboratory comparisons. We applied SNP typing to study strains with a novel allele for the pertussis toxin promoter, ptxP3, which have a worldwide distribution and which have replaced the resident ptxP1 strains in the last 20 years. Previously, we showed that ptxP3 strains showed increased pertussis toxin expression and that their emergence was associated with increased notification in The Netherlands. SNP typing showed that the ptxP3 strains isolated in the Americas, Asia, Australia and Europe formed a monophyletic branch which recently diverged from ptxP1 strains. Two predominant ptxP3 SNP types were identified which spread worldwide. The widespread use of SNP typing will enhance our understanding of the evolution and global epidemiology of B. pertussis.

Multilocus sequence typing of Bordetella pertussis based on surface protein genes.

Despite more than 50 years of vaccination, Bordetella pertussis has remained endemic in The Netherlands, causing epidemic outbreaks every 3 to 5 years. Strain variation may play a role in the persistence of B. pertussis and was studied by sequencing 15 genes coding for surface proteins, including genes for all five components of acellular pertussis vaccines: pertussis toxin (Ptx), pertactin (Prn), filamentous hemagglutinin, and fimbriae (Fim2 and Fim3). A low level of allelic variation was observed, confirming a recent evolutionary origin of B. pertussis. In modern isolates, polymorphism was observed only in prn, ptxS1, ptxS3, and tcfA. Polymorphism in ptxS1, ptxS3, and tcfA was used to categorize isolates in multilocus sequence types (MLSTs). Analysis of Dutch isolates from 1949 to 1999 revealed five MLSTs, which showed a highly dynamic temporal behavior. We observed significant changes in the MLSTs after the introduction of pertussis vaccination in The Netherlands. Epidemic years were found to be associated with the expansion of MLST-4 or MLST-5. MLST-5 showed a remarkable expansion from 10% in 1997 to 80% in 1999. The MLST analysis was extended to a number of widely separated geographic regions: Finland, Italy, Japan, and the United States. MLST-4 and MLST-5 were found to dominate in Italy and the United States. In Finland and Japan, MLST-3 and MLST-2, respectively, were predominant. Thus, although each region showed distinctive MLST frequencies, in three of the five regions MLST-4 and MLST-5 were predominant. These types may represent newly emerged, successful clones. The identification of highly successful clones may shed light on the question of how B. pertussis is able to maintain itself in vaccinated populations.