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.

Comparative genomics of prevaccination and modern Bordetella pertussis strains.

BACKGROUND: Despite vaccination since the 1950s, pertussis has persisted and resurged. It remains a major cause of infant death worldwide and is the most prevalent vaccine-preventable disease in developed countries. The resurgence of pertussis has been associated with the expansion of Bordetella pertussis strains with a novel allele for the pertussis toxin (Ptx) promoter, ptxP3, which have replaced resident ptxP1 strains. Compared to ptxP1 strains, ptxP3 produce more Ptx resulting in increased virulence and immune suppression. To elucidate how B. pertussis has adapted to vaccination, we compared genome sequences of two ptxP3 strains with four strains isolated before and after the introduction vaccination. RESULTS: The distribution of SNPs in regions involved in transcription and translation suggested that changes in gene regulation play an important role in adaptation. No evidence was found for acquisition of novel genes. Modern strains differed significantly from prevaccination strains, both phylogenetically and with respect to particular alleles. The ptxP3 strains were found to have diverged recently from modern ptxP1 strains. Differences between ptxP3 and modern ptxP1 strains included SNPs in a number of pathogenicity-associated genes. Further, both gene inactivation and reactivation was observed in ptxP3 strains relative to modern ptxP1 strains. CONCLUSIONS: Our work suggests that B. pertussis adapted by successive accumulation of SNPs and by gene (in)activation. In particular changes in gene regulation may have played a role in adaptation.

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.

The identification of response regulator-specific binding sites reveals new roles of two-component systems in Bacillus cereus and closely related low-GC Gram-positives.

In bacteria, environmental challenges are often translated into a transcriptional response via the cognate response regulators (RRs) of specialized two-component systems (TCSs). A phylogenetic footprinting/shadowing approach was designed and used to identify many novel RR-specific operators for species of the Bacillus cereus group and related Gram-positives. Analysis of the operator sequences revealed characteristic traits for each RR subfamily. For instance, operators related to the largest subfamily (OmpR) typically consisted of direct repeats (e.g. TTAAGA-N5-TTAAGA), whereas operators related to the second largest family (NarL) consisted of inverted repeats (e.g. ATGACA-N2-TGTCAT). This difference indicates a fundamentally different organization of the bound RR dimers between the two subfamilies. Moreover, the identification of the specific operator motifs allowed relating several RRs to a minimal regulon and thereby to a characteristic transcriptional response. Mostly, these regulons comprised genes encoding transport systems, suggesting a direct coupling of stimulus perception to the transport of target compounds. New biological roles could be attributed to various TCSs, including roles in cytochrome c biogenesis (HssRS), transport of carbohydrates, peptides and/or amino acids (YkoGH, LytSR), and resistance to toxic ions (LiaSR), antimicrobial peptides (BceRS) and beta-lactam antibiotics (BacRS, YcbLM). As more and more bacterial genome sequences are becoming available, the use of comparative analyses such as the approach applied in this study will further increase our knowledge of bacterial signal transduction mechanisms and provide directions for the assessment of their role in bacterial performance and survival strategies.

Complete Genome Sequences of 11 Bordetella pertussis Strains Representing the Pandemic ptxP3 Lineage.

Pathogen adaptation has contributed to the resurgence of pertussis. To facilitate our understanding of this adaptation we report here 11 completely closed and annotated Bordetella pertussis genomes representing the pandemic ptxP3 lineage. Our analyses included six strains which do not produce the vaccine components pertactin and/or filamentous hemagglutinin.

Complete Genome Sequences of Bordetella pertussis Isolates B1917 and B1920, Representing Two Predominant Global Lineages.

Bordetella pertussis is the causative agent of pertussis, a disease which has resurged despite vaccination. We report the complete, annotated genomes of isolates B1917 and B1920, representing two lineages predominating globally in the last 50 years. The B1917 lineage has been associated with the resurgence of pertussis in the 1990s.

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.

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.

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.