Comparative genomics of Bordetella pertussis reveals progressive gene loss in Finnish strains

Genome Characteristic of Bordetella parapertussis Isolated from Iran

Pertussis also known as whooping cough is a respiratory infection in humans particularly with severe symptoms in infants and usually caused by Bordetella pertussis. However, Bordetella parapertussis can also cause a similar clinical syndrome. During 2012 to 2015, from nasal swabs sent from different provinces to the pertussis reference laboratory of Pasture Institute of Iran for pertussis confirmation, seven B. parapertussis isolates were identified by bacterial culture, biochemical tests, and the presence of IS1001 insertion in the genome. The expression of pertactin (Prn) as one the major virulence factor for bacterial adhesion was investigated using western blot. Moreover, the genomic characteristic of one recently collected isolate, IRBP134, from a seven-month infant was investigated using Illumina NextSeq sequencing protocol. The results revealed the genome with G+C content 65% and genome size 4.7 Mbp. A total of 81 single nucleotide polymorphisms and 13 short insertions and deletions were found in the genome compared to the B. parapertussis 12822 as a reference genome showing ongoing evolutionary changes. A phylogeny relationship of IRBP134 was also investigated using global B. parapertussis available genomes.

Towards comprehensive understanding of bacterial genetic diversity: large-scale amplifications in Bordetella pertussis and Mycobacterium tuberculosis

Bacterial genetic diversity is often described solely using base-pair changes despite a wide variety of other mutation types likely being major contributors. Tandem duplication/amplifications are thought to be widespread among bacteria but due to their often-intractable size and instability, comprehensive studies of these mutations are rare. We define a methodology to investigate amplifications in bacterial genomes based on read depth of genome sequence data as a proxy for copy number. We demonstrate the approach with Bordetella pertussis, whose insertion sequence element-rich genome provides extensive scope for amplifications to occur. Analysis of data for 2430 B. pertussis isolates identified 272 putative amplifications, of which 94 % were located at 11 hotspot loci. We demonstrate limited phylogenetic connection for the occurrence of amplifications, suggesting unstable and sporadic characteristics. Genome instability was further described in vitro using long-read sequencing via the Nanopore platform, which revealed that clonally derived laboratory cultures produced heterogenous populations rapidly. We extended this research to analyse a population of 1000 isolates of another important pathogen, Mycobacterium tuberculosis. We found 590 amplifications in M. tuberculosis, and like B. pertussis, these occurred primarily at hotspots. Genes amplified in B. pertussis include those involved in motility and respiration, whilst in M. tuberuclosis, functions included intracellular growth and regulation of virulence. Using publicly available short-read data we predicted previously unrecognized, large amplifications in B. pertussis and M. tuberculosis. This reveals the unrecognized and dynamic genetic diversity of B. pertussis and M. tuberculosis, highlighting the need for a more holistic understanding of bacterial genetics.

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.

Genomic epidemiology of Iranian Bordetella pertussis: 50 years after the implementation of whole cell vaccine

Pertussis caused by Bordetella pertussis, remains a public health problem worldwide, despite high vaccine coverage in infants and children in many countries. Iran has been using whole cell vaccine for the last 50 years with more than 95% vaccination rate since 1988 and has experienced pertussis resurgence in recent years. Here, we sequenced 55 B. pertussis isolates mostly collected from three provinces with the highest number of pertussis cases in Iran, including Tehran, Mazandaran, and Eastern-Azarbayjan from the period of 2008-2016. Most isolates carried ptxP3/prn2 alleles (42/55, 76%), the same genotype as isolates circulating in acellular vaccine-administrating countries. The second most frequent genotype was ptxP3/prn9 (8/55, 14%). Only three isolates (5%) were ptxP1. Phylogenetic analysis showed that Iranian ptxP3 isolates can be divided into eight clades (Clades 1-8) with no temporal association. Most of the isolates from Tehran grouped together as one distinctive clade (Clade 8) with six unique single nucleotide polymorphisms (SNPs). In addition, the prn9 isolates were grouped together as Clade 5 with 12 clade-supporting SNPs. No pertactin deficient isolates were found among the 55 Iranian isolates. Our findings suggest that there is an ongoing adaptation and evolution of B. pertussis regardless of the types of vaccine used.

How Genomics Is Changing What We Know About the Evolution and Genome of Bordetella pertussis

The evolution of Bordetella pertussis from a common ancestor similar to Bordetella bronchiseptica has occurred through large-scale gene loss, inactivation and rearrangements, largely driven by the spread of insertion sequence element repeats throughout the genome. B. pertussis is widely considered to be monomorphic, and recent evolution of the B. pertussis genome appears to, at least in part, be driven by vaccine-based selection. Given the recent global resurgence of whooping cough despite the wide-spread use of vaccination, a more thorough understanding of B. pertussis genomics could be highly informative. In this chapter we discuss the evolution of B. pertussis, including how vaccination is changing the circulating B. pertussis population at the gene-level, and how new sequencing technologies are revealing previously unknown levels of inter- and intra-strain variation at the genome-level.

Resolving the complex Bordetella pertussis genome using barcoded nanopore sequencing

The genome of Bordetella pertussis is complex, with high G+C content and many repeats, each longer than 1000 bp. Long-read sequencing offers the opportunity to produce single-contig B. pertussis assemblies using sequencing reads which are longer than the repetitive sections, with the potential to reveal genomic features which were previously unobservable in multi-contig assemblies produced by short-read sequencing alone. We used an R9.4 MinION flow cell and barcoding to sequence five B. pertussis strains in a single sequencing run. We then trialled combinations of the many nanopore user community-built long-read analysis tools to establish the current optimal assembly pipeline for B. pertussis genome sequences. This pipeline produced closed genome sequences for four strains, allowing visualization of inter-strain genomic rearrangement. Read mapping to the Tohama I reference genome suggests that the remaining strain contains an ultra-long duplicated region (almost 200 kbp), which was not resolved by our pipeline; further investigation also revealed that a second strain that was seemingly resolved by our pipeline may contain an even longer duplication, albeit in a small subset of cells. We have therefore demonstrated the ability to resolve the structure of several B. pertussis strains per single barcoded nanopore flow cell, but the genomes with highest complexity (e.g. very large duplicated regions) remain only partially resolved using the standard library preparation and will require an alternative library preparation method. For full strain characterization, we recommend hybrid assembly of long and short reads together; for comparison of genome arrangement, assembly using long reads alone is sufficient.

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.

Complete Genome Sequence of Bordetella pertussis Pelita III, the Production Strain for an Indonesian Whole-Cell Pertussis Vaccine

PT Bio Farma, the sole World Health Organization-approved Indonesian vaccine producer, manufactures a whole-cell whooping cough vaccine (wP) that, as part of a pentavalent diphtheria-tetanus-pertussis/hepatitis B/Haemophilus influenzae b (DTP/HB/Hib) vaccine, is used in Indonesia and many other countries. We report here the whole-genome sequence for Bordetella pertussis Pelita III, PT Bio Farma’s wP production strain.

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.

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.

Bordetella pertussis Isolates from Argentinean Whooping Cough Patients Display Enhanced Biofilm Formation Capacity Compared to Tohama I Reference Strain

Pertussis is a highly contagious disease mainly caused by Bordetella pertussis. Despite the massive use of vaccines, since the 1950s the disease has become re-emergent in 2000 with a shift in incidence from infants to adolescents and adults. Clearly, the efficacy of current cellular or acellular vaccines, formulated from bacteria grown in stirred bioreactors is limited, presenting a challenge for future vaccine development. For gaining insights into the role of B. pertussis biofilm development for host colonization and persistence within the host, we examined the biofilm forming capacity of eight argentinean clinical isolates recovered from 2001 to 2007. All clinical isolates showed an enhanced potential for biofilm formation compared to the reference strain Tohama I. We further selected the clinical isolate B. pertussis 2723, exhibiting the highest biofilm biomass production, for quantitative proteomic profiling by means of two-dimensional fluorescence difference gel electrophoresis (2D-DIGE) coupled with mass spectrometry, which was accompanied by targeted transcriptional analysis. Results revealed an elevated expression of several virulence factors, including adhesins involved in biofilm development. In addition, we observed a higher expression of energy metabolism enzymes in the clinical isolate compared to the Tohama I strain. Furthermore, all clinical isolates carried a polymorphism in the bvgS gene. This mutation was associated to an increased sensitivity to modulation and a faster rate of adhesion to abiotic surfaces. Thus, the phenotypic biofilm characteristics shown by the clinical isolates might represent an important, hitherto underestimated, adaptive strategy for host colonization and long time persistence within the host.

Laboratory Diagnosis of Pertussis

The introduction of vaccination in the 1950s significantly reduced the morbidity and mortality of pertussis. However, since the 1990s, a resurgence of pertussis has been observed in vaccinated populations, and a number of causes have been proposed for this phenomenon, including improved diagnostics, increased awareness, waning immunity, and pathogen adaptation. The resurgence of pertussis highlights the importance of standardized, sensitive, and specific laboratory diagnoses, the lack of which is responsible for the large differences in pertussis notifications between countries. Accurate laboratory diagnosis is also important for distinguishing between the several etiologic agents of pertussis-like diseases, which involve both viruses and bacteria. If pertussis is diagnosed in a timely manner, antibiotic treatment of the patient can mitigate the symptoms and prevent transmission. During an outbreak, timely diagnosis of pertussis allows prophylactic treatment of infants too young to be (fully) vaccinated, for whom pertussis is a severe, sometimes fatal disease. Finally, reliable diagnosis of pertussis is required to reveal trends in the (age-specific) disease incidence, which may point to changes in vaccine efficacy, waning immunity, and the emergence of vaccine-adapted strains. Here we review current approaches to the diagnosis of pertussis and discuss their limitations and strengths. In particular, we emphasize that the optimal diagnostic procedure depends on the stage of the disease, the age of the patient, and the vaccination status of the patient.

Genetic Variation of Bordetella pertussis in Austria

In Austria, vaccination coverage against Bordetella pertussis infections during infancy is estimated at around 90%. Within the last years, however, the number of pertussis cases has increased steadily, not only in children but also in adolescents and adults, indicating both insufficient herd immunity and vaccine coverage. Waning immunity in the host and/or adaptation of the bacterium to the immunised hosts could contribute to the observed re-emergence of pertussis. In this study we therefore addressed the genetic variability in B. pertussis strains from several Austrian cities. Between the years 2002 and 2008, 110 samples were collected from Vienna (n = 32), Linz (n = 63) and Graz (n = 15) by nasopharyngeal swabs. DNA was extracted from the swabs, and bacterial sequence polymorphisms were examined by MLVA (multiple-locus variable number of tandem repeat analysis) (n = 77), by PCR amplification and conventional Sanger sequencing of the polymorphic regions of the prn (pertactin) gene (n = 110), and by amplification refractory mutation system quantitative PCR (ARMS-qPCR) (n = 110) to directly address polymorphisms in the genes encoding two pertussis toxin subunits (ptxA and ptxB), a fimbrial adhesin (fimD), tracheal colonisation factor (tcfA), and the virulence sensor protein (bvgS). Finally, the ptxP promoter region was screened by ARMS-qPCR for the presence of the ptxP3 allele, which has been associated with elevated production of pertussis toxin. The MLVA analysis revealed the highest level of polymorphisms with an absence of MLVA Type 29, which is found outside Austria. Only Prn subtypes Prn1/7, Prn2 and Prn3 were found with a predominance of the non-vaccine type Prn2. The analysis of the ptxA, ptxB, fimD, tcfA and bvgS polymorphisms showed a genotype mixed between the vaccine strain Tohama I and a clinical isolate from 2006 (L517). The major part of the samples (93%) displayed the ptxP3 allele. The consequences for the vaccination strategy are discussed.

Extensive strain-level copy-number variation across human gut microbiome species

Within each bacterial species, different strains may vary in the set of genes they encode or in the copy number of these genes. Yet, taxonomic characterization of the human microbiota is often limited to the species level or to previously sequenced strains, and accordingly, the prevalence of intra-species variation, its functional role, and its relation to host health remain unclear. Here, we present a comprehensive large-scale analysis of intra-species copy-number variation in the gut microbiome, introducing a rigorous computational pipeline for detecting such variation directly from shotgun metagenomic data. We uncover a large set of variable genes in numerous species and demonstrate that this variation has significant functional and clinically relevant implications. We additionally infer intra-species compositional profiles, identifying population structure shifts and the presence of yet uncharacterized variants. Our results highlight the complex relationship between microbiome composition and functional capacity, linking metagenome-level compositional shifts to strain-level variation.

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.

Genome-wide gene expression analysis of Bordetella pertussis isolates associated with a resurgence in pertussis: elucidation of factors involved in the increased fitness of epidemic strains

Bordetella pertussis (B. pertussis) is the causative agent of whooping cough, which is a highly contagious disease in the human respiratory tract. Despite vaccination since the 1950s, pertussis remains the most prevalent vaccine-preventable disease in developed countries. A recent resurgence pertussis is associated with the expansion of B. pertussis strains with a novel allele for the pertussis toxin (ptx) promoter ptxP3 in place of resident ptxP1 strains. The recent expansion of ptxP3 strains suggests that these strains carry mutations that have increased their fitness. Compared to the ptxP1 strains, ptxP3 strains produce more Ptx, which results in increased virulence and immune suppression. In this study, we investigated the contribution of gene expression changes of various genes on the increased fitness of the ptxP3 strains. Using genome-wide gene expression profiling, we show that several virulence genes had higher expression levels in the ptxP3 strains compared to the ptxP1 strains. We provide the first evidence that wildtype ptxP3 strains are better colonizers in an intranasal mouse infection model. This study shows that the ptxP3 mutation and the genetic background of ptxP3 strains affect fitness by contributing to the ability to colonize in a mouse infection model. These results show that the genetic background of ptxP3 strains with a higher expression of virulence genes contribute to increased fitness.

Seroprevalence of pertussis in Senegal: a prospective study

BACKGROUND

Pertussis, also known as whooping cough, is a vaccine-preventable respiratory disease caused by Bordetella pertussis infection, against which Senegalese children are immunized with the diphtheria-tetanus-whole cell pertussis vaccine (DTwP). Seroepidemiology of pertussis has been widely described in industrialized countries, but rare are the studies referring to it in developing countries.

METHODS

We conducted a longitudinal survey in Northern Senegal to investigate the epidemiology of B. pertussis by evaluating the IgG antibody (Ab) response against pertussis toxin (PT). A cohort of 410 children aged 1 to 9 from five villages in the Middle Senegal River Valley were followed-up for 18 months. During that period, five visits were made to assess the immunological status of the children.

PRINCIPAL FINDINGS

PT-specific IgG responses were significantly different according to age. Until the age of 3, there was a decrease in the Ab response, which then increased in the older groups. Assessment of IgG antibodies to PT (IgG-PT) suggested evidence of recent exposures to the pathogen. Surprisingly, in one of the five villages the average Ab response to PT was very low at all ages during the first 6 months of the study. At the third visit, IgG-PT concentrations peaked to very high levels, to slightly decline at the end of the survey. This indicates an outbreak of B. pertussis, whereas in the other villages a pertussis endemic profile could be observed.

CONCLUSIONS

Pertussis is endemic in Northern Senegal despite the introduction of vaccination. The circulation of the bacteria seems to differ between geographic locations and over time. A more complete understanding of the epidemiology of pertussis and its environmental determinants could provide information to adapt vaccination programs.

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.

Insertion sequences shared by Bordetella species and implications for the biological diagnosis of pertussis syndrome

The molecular diagnosis of pertussis and parapertussis syndromes is based on the detection of insertion sequences (IS) 481 and 1001, respectively. However, these IS are also detected in the genomes of various Bordetella species, such that they are not specific for either B. pertussis or B. parapertussis. Therefore, we screened the genome of recently circulating isolates of Bordetella species to compare the prevalence of IS481, IS1001 and, also IS1002 with previously published data and to sequence all IS detected. We also investigated whether the numbers of IS481 and IS1001 copies vary in recently circulating isolates of the different Bordetella species. We used the polymerase chain reaction (PCR) method for screening the genome of circulating isolates and to prepare the fragments for sequencing. We used Southern blotting and quantitative real-time PCR for quantification of the numbers of IS. We found no significant diversity in the sequences of the IS harboured in the genomes of the Bordetella isolates screened, except for a 71-nucleotide deletion from IS1002 in B. bronchiseptica. The IS copy numbers in the genome of recently circulating isolates were similar to those in reference strains. Our results confirm that biological diagnosis targeting the IS481 and IS1001 elements are not specific and detect the species B. pertussis, B. holmesii and B. bronchiseptica (IS481), and B. parapertussis and B. bronchiseptica (IS1001).

Evolution of French Bordetella pertussis and Bordetella parapertussis isolates: increase of Bordetellae not expressing pertactin

Bordetella pertussis and Bordetella parapertussis are closely related bacterial agents of whooping cough. Whole-cell pertussis (wP) vaccine was introduced in France in 1959. Acellular pertussis (aP) vaccine was introduced in 1998 as an adolescent booster and was rapidly generalized to the whole population, changing herd immunity by specifically targeting the virulence of the bacteria. We performed a temporal analysis of all French B. pertussis and B. parapertussis isolates collected since 2000 under aP vaccine pressure, using pulsed-field gel electrophoresis (PFGE), genotyping and detection of expression of virulence factors. Particular isolates were selected according to their different phenotype and PFGE type and their characteristics were analysed using the murine model of respiratory infection and in vitro cell cytotoxic assay. Since the introduction of the aP vaccines there has been a steady increase in the number of B. pertussis and B. parapertussis isolates collected that are lacking expression of pertactin. These isolates seem to be as virulent as those expressing all virulence factors according to animal and cellular models of infection. Whereas wP vaccine-induced immunity led to a monomorphic population of B. pertussis, aP vaccine-induced immunity enabled the number of circulating B. pertussis and B. parapertussis isolates not expressing virulence factors to increase, sustaining our previous hypothesis.

The Bordetella pertussis model of exquisite gene control by the global transcription factor BvgA

Bordetella pertussis causes whooping cough, an infectious disease that is reemerging despite widespread vaccination. A more complete understanding of B. pertussis pathogenic mechanisms will involve unravelling the regulation of its impressive arsenal of virulence factors. Here we review the action of the B. pertussis response regulator BvgA in the context of what is known about bacterial RNA polymerase and various modes of transcription activation. At most virulence gene promoters, multiple dimers of phosphorylated BvgA (BvgA~P) bind upstream of the core promoter sequence, using a combination of high- and low-affinity sites that fill through cooperativity. Activation by BvgA~P is typically mediated by a novel form of class I/II mechanisms, but two virulence genes, fim2 and fim3, which encode serologically distinct fimbrial subunits, are regulated using a previously unrecognized RNA polymerase/activator architecture. In addition, the fim genes undergo phase variation because of an extended cytosine (C) tract within the promoter sequences that is subject to slipped-strand mispairing during replication. These sophisticated systems of regulation demonstrate one aspect whereby B. pertussis, which is highly clonal and lacks the extensive genetic diversity observed in many other bacterial pathogens, has been highly successful as an obligate human pathogen.

Prevalence and genetic characterization of pertactin-deficient Bordetella pertussis in Japan

The adhesin pertactin (Prn) is one of the major virulence factors of Bordetella pertussis, the etiological agent of whooping cough. However, a significant prevalence of Prn-deficient (Prn⁻) B. pertussis was observed in Japan. The Prn⁻ isolate was first discovered in 1997, and 33 (27%) Prn⁻ isolates were identified among 121 B. pertussis isolates collected from 1990 to 2009. Sequence analysis revealed that all the Prn⁻ isolates harbor exclusively the vaccine-type prn1 allele and that loss of Prn expression is caused by 2 different mutations: an 84-bp deletion of the prn signal sequence (prn1ΔSS, n = 24) and an IS481 insertion in prn1 (prn1::IS481, n = 9). The frequency of Prn⁻ isolates, notably those harboring prn1ΔSS, significantly increased since the early 2000s, and Prn⁻ isolates were subsequently found nationwide. Multilocus variable-number tandem repeat analysis (MLVA) revealed that 24 (73%) of 33 Prn⁻ isolates belong to MLVA-186, and 6 and 3 Prn⁻ isolates belong to MLVA-194 and MLVA-226, respectively. The 3 MLVA types are phylogenetically closely related, suggesting that the 2 Prn⁻ clinical strains (harboring prn1ΔSS and prn1::IS481) have clonally expanded in Japan. Growth competition assays in vitro also demonstrated that Prn⁻ isolates have a higher growth potential than the Prn⁺ back-mutants from which they were derived. Our observations suggested that human host factors (genetic factors and immune status) that select for Prn⁻ strains have arisen and that Prn expression is not essential for fitness under these conditions.

Differences in the genomic content of Bordetella pertussis isolates before and after introduction of pertussis vaccines in four European countries

Resurgence of pertussis has been observed in many countries with high vaccination coverage and clonal expansion of certain Bordetella pertussis strains has been associated with recent epidemics in Europe. It is known that vaccinations have selected strains which are different from those used for vaccine production. However, little is known about the differences in genomic content of strains circulating before the vaccination was introduced. In this study, we compared the genomes of 39 vaccine strains and old clinical isolates (isolated 1941-1984) collected from Finland (n = 5), Poland (n = 14), Serbia (n = 10) and the UK (n = 10). The analysis included genotyping, pulsed-field gel electrophoresis (PFGE) and comparative genomic hybridisation (CGH). Compared to the strain Tohama I, the European isolates analyzed have lost three major regions of difference (RD3, 5 and 29). However, difference in frequency of the absent RDs 3 (BP0910A-BP0934), 5 (BP1135-BP1141) or 29 (BP1225) was observed among isolates from the four countries. Of the isolates with absent RD5, half had also a duplicated region in the genome. All four RDs (RD22 (BB0535-BB0541), 23 (BB0916-BB0921), 24 (BB1140-BB1158) and 26 (BB4880-BB4888)) absent in Tohama I were present in majority of the tested isolates. Results obtained from PFGE analysis correlated well with those of CGH. Recently a novel pertussis toxin promoter allele (ptxP3) was described. Isolates with ptxP3 have replaced resident ptxP1 isolates in the countries where this was investigated. When the recent isolates, collected in 2000-2004, selected from the four countries were examined, the ptxP3 allele was found in all countries except Poland. Our result indicates that at least three clusters of B. pertussis circulated in Europe in pre- and early vaccine era and their genomes were distinct from that of the reference strain Tohama I. Although progressive gene loss occurs in B. pertussis population with time, difference in frequency of the lost genes were observed among isolates from the four countries. The observed differences in genomic content might be vaccine-driven.

Differential expression of type III effector BteA protein due to IS481 insertion in Bordetella pertussis

BACKGROUND

Bordetella pertussis is the primary etiologic agent of the disease pertussis. Universal immunization programs have contributed to a significant reduction in morbidity and mortality of pertussis; however, incidence of the disease, especially in adolescents and adults, has increased in several countries despite high vaccination coverage. During the last three decades, strains of Bordetella pertussis in circulation have shifted from the vaccine-type to the nonvaccine-type in many countries. A comparative proteomic analysis of the strains was performed to identify protein(s) involved in the type shift.

METHODOLOGY/PRINCIPAL FINDING

Proteomic analysis identified one differentially expressed protein in the B. pertussis strains: the type III cytotoxic effector protein BteA, which is responsible for host cell death in Bordetella bronchiseptica infections. Immunoblot analysis confirmed the prominent expression of BteA protein in the nonvaccine-type strains but not in the vaccine-type strains. Sequence analysis of the vaccine-type strains revealed an IS481 insertion in the 5' untranslated region of bteA, -136 bp upstream of the bteA start codon. A high level of bteA transcripts from the IS481 promoter was detected in the vaccine-type strains, indicating that the transcript might be an untranslatable form. Furthermore, BteA mutant studies demonstrated that BteA expression in the vaccine-type strains is down-regulated by the IS481 insertion.

CONCLUSION/SIGNIFICANCE

The cytotoxic effector BteA protein is expressed at higher levels in B. pertussis nonvaccine-type strains than in vaccine-type strains. This type-dependent expression is due to an insertion of IS481 in B. pertussis clinical strains, suggesting that augmented expression of BteA protein might play a key role in the type shift of B. pertussis.

Interactions between Bordetella pertussis and the complement inhibitor factor H

Bordetella pertussis causes whooping cough in humans, a highly contagious disease of the upper respiratory tract. An increase in cases of whooping cough in adolescents and adults in many countries has been reported, despite high immunization rates in children. To efficiently colonize the host the bacteria have to resist complement, the first defence line of innate immunity. B. pertussis has previously been shown to bind the classical pathway inhibitors C4b-binding protein and C1-inhibitor being thereby able to escape the classical pathway of complement. In this study recent clinical isolates of B. pertussis and B. parapertussis were found to survive alternative pathway attack in fresh non-immune serum better than the reference B. pertussis strain, Tohama I. By using adsorption assays, flow cytometry and a radioligand binding assay we observed that both B. pertussis and B. parapertussis bound the alternative pathway inhibitor factor H (FH) from normal human serum. The surface attached FH maintained its complement regulatory activity and promoted factor I-mediated cleavage of C3b. The main binding region was located to the C-terminal part of FH, into short consensus repeat domains 19-20. In contrast, the avian pathogen B. avium did not bind FH and was sensitive to the alternative pathway of human complement. In conclusion, the human pathogens B. pertussis and B. parapertussis are able to evade the alternative complement pathway by surface acquisition of the host complement regulator FH.

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.

Dose response of attenuated Bordetella pertussis BPZE1-induced protection in mice

Despite the availability of efficacious vaccines, the incidence of whooping cough is still high in many countries and is even increasing in countries with high vaccine coverage. Most severe and life-threatening pertussis cases occur in infants who are too young to be sufficiently protected by current vaccine regimens. As a potential solution to this problem, we have developed an attenuated live Bordetella pertussis vaccine strain, named BPZE1. Here, we show that after a single administration, BPZE1 induces dose-dependent protection against challenge with virulent B. pertussis in low-dose and in high-dose intranasal mouse lung colonization models. In addition, we observed BPZE1 dose-dependent antibody titers to B. pertussis antigens, as well as cell-mediated immunity, evidenced by the amounts of gamma interferon (IFN-gamma) released from spleen cells upon stimulation with B. pertussis antigens. These two parameters may perhaps be used as readouts in clinical trials in humans that are currently being planned.

Changes in the genomic content of circulating Bordetella pertussis strains isolated from the Netherlands, Sweden, Japan and Australia: adaptive evolution or drift?

Background

Bordetella pertussis is the causative agent of human whooping cough (pertussis) and is particularly severe in infants. Despite worldwide vaccinations, whooping cough remains a public health problem. A significant increase in the incidence of whooping cough has been observed in many countries since the 1990s. Several reasons for the re-emergence of this highly contagious disease have been suggested. A particularly intriguing possibility is based on evidence indicating that pathogen adaptation may play a role in this process. In an attempt to gain insight into the genomic make-up of B. pertussis over the last 60 years, we used an oligonucleotide DNA microarray to compare the genomic contents of a collection of 171 strains of B. pertussis isolates from different countries.

Results

The CGH microarray analysis estimated the core genome of B. pertussis, to consist of 3,281 CDSs that are conserved among all B. pertussis strains, and represent 84.8% of all CDSs found in the 171 B. pertussis strains. A total of 64 regions of difference consisting of one or more contiguous CDSs were identified among the variable genes. CGH data also revealed that the genome size of B. pertussis strains is decreasing progressively over the past 60 years. Phylogenetic analysis of microarray data generated a minimum spanning tree that depicted the phylogenetic structure of the strains. B. pertussis strains with the same gene content were found in several different countries. However, geographic specificity of the B. pertussis strains was not observed. The gene content was determined to highly correlate with the ptxP-type of the strains.

Conclusions

An overview of genomic contents of a large collection of isolates from different countries allowed us to derive a core genome and a phylogenetic structure of B. pertussis. Our results show that B. pertussis is a dynamic organism that continues to evolve.

Bordetella pertussis and vaccination: the persistence of a genetically monomorphic pathogen

Before childhood vaccination was introduced in the 1950s, pertussis or whooping cough was a major cause of infant death worldwide. Widespread vaccination of children was successful in significantly reducing morbidity and mortality. However, despite vaccination, pertussis has persisted and, in the 1990s, resurged in a number of countries with highly vaccinated populations. Indeed, pertussis has become the most prevalent vaccine-preventable disease in developed countries with estimated infection frequencies of 1-6%. Recently vaccinated children are well protected against pertussis disease and its increase is mainly seen in adolescents and adults in which disease symptoms are often mild. The etiologic agent of pertussis, Bordetella pertussis, is extremely monomorphic and its ability to persist in the face of intensive vaccination is intriguing. Numerous studies have shown that B. pertussis populations changed after the introduction of vaccination suggesting adaptation. These adaptations did not involve the acquisition of novel genes but small genetic changes, mainly SNPs, and occurred in successive steps in a period of 40 years. The earliest adaptations resulted in antigenic divergence with vaccine strains. More recently, strains emerged with increased pertussis toxin (Ptx) production. Here I argue that the resurgence of pertussis is the compound effect of pathogen adaptation and waning immunity. I propose that the removal by vaccination of naïve infants as the major source for transmission was the crucial event which has driven the changes in B. pertussis populations. This has selected for strains which are more efficiently transmitted by primed hosts in which immunity has waned. The adaptation of B. pertussis to primed hosts involved delaying an effective immune response by antigenic divergence with vaccine strains and by increasing immune suppression through higher levels of Ptx production. Higher levels of Ptx may also benefit transmission by enhancing clinical symptoms. The study of B. pertussis populations has not only increased our understanding of pathogen evolution, but also suggests way to improve pertussis vaccines, underlining the public health significance of population-based studies of pathogens.

Crucial role for insertion sequence elements in Lactobacillus helveticus evolution as revealed by interstrain genomic comparison

Lactobacillus helveticus is a versatile dairy bacterium found to possess heterogeneous genotypes depending on the ecosystem from which it was isolated. The recently published genome sequence showed the remarkable flexibility of its structure, demonstrated by a substantial level of insertion sequence (IS) element expansion in association with massive gene decay. To assess this diversity and examine the level of genome plasticity within the L. helveticus species, an array-based comparative genome hybridization (aCGH) experiment was designed in which 10 strains were analyzed. The aCGH experiment revealed 16 clusters of open reading frames (ORFs) flanked by IS elements. Four of these ORFs are associated with restriction/modification which may have played a role in accelerated evolution of strains in a commercially intensive ecosystem undoubtedly challenged through successive phage attack. Furthermore, analysis of the IS-flanked clusters demonstrated that the most frequently encountered ISs were also those most abundant in the genome (IS1201, ISL2, ISLhe1, ISLhe2, ISLhe65, and ISLhe63). These findings contribute to the overall viewpoint of the versatile character of IS elements and the role they may play in bacterial genome plasticity.

First report and detailed characterization of B. pertussis isolates not expressing Pertussis Toxin or Pertactin

Bordetella pertussis isolates not expressing Pertussis Toxin (PT) or Pertactin (PRN) have been collected, for the first time in 2007, in France, a highly vaccinated country with acellular vaccines. Non-expression was due to deletion of the entire ptx locus, to IS481 insertion in the prn gene or deletion of a part of this gene. Genome sequencing does not indicate any regions of differences when compared to other circulating isolates. It nevertheless shows some sequence differences and an increased number of repeated sequences. The infant infected by the isolate not expressing pertussis toxin, did not present hyperlymphocytosis. All isolates were found less pathogen in animal or cellular models; their circulation raises the problem of clinical and biological diagnoses.

Bordetella pertussis isolates in Finland: serotype and fimbrial expression

BACKGROUND

Bordetella pertussis causes whooping cough or pertussis in humans. It produces several virulence factors, of which the fimbriae are considered adhesins and elicit immune responses in the host. B. pertussis has three distinct serotypes Fim2, Fim3 or Fim2,3. Generally, B. pertussis Fim2 strains predominate in unvaccinated populations, whereas Fim3 strains are often isolated in vaccinated populations. In Finland, pertussis vaccination was introduced in 1952. The whole-cell vaccine contained two strains, 18530 (Fim3) since 1962 and strain 1772 (Fim2,3) added in 1976. After that the vaccine has remained the same until 2005 when the whole-cell vaccine was replaced by the acellular vaccine containing pertussis toxin and filamentous hemagglutinin. Our aims were to study serotypes of Finnish B. pertussis isolates from 1974 to 2006 in a population with > 90% vaccination coverage and fimbrial expression of the isolates during infection. Serotyping was done by agglutination and serotype-specific antibody responses were determined by blocking ELISA.

RESULTS

Altogether, 1,109 isolates were serotyped. Before 1976, serotype distributions of Fim2, Fim3 and Fim2,3 were 67%, 19% and 10%, respectively. From 1976 to 1998, 94% of the isolates were Fim2 serotype. Since 1999, the frequency of Fim3 strains started to increase and reached 83% during a nationwide epidemic in 2003. A significant increase in level of serum IgG antibodies against purified fimbriae was observed between paired sera of 37 patients. The patients infected by Fim3 strains had antibodies which blocked the binding of monoclonal antibodies to Fim3 but not to Fim2. Moreover, about one third of the Fim2 strain infected patients developed antibodies capable of blocking of binding of both anti-Fim2 and Fim3 monoclonal antibodies.

CONCLUSION

Despite extensive vaccinations in Finland, B. pertussis Fim2 strains were the most common serotype. Emergence of Fim3 strains started in 1999 and coincided with nationwide epidemics. Results of serotype-specific antibody responses suggest that Fim2 strains could express Fim3 during infection, showing a difference in fimbrial expression between in vivo and in vitro.

Pertussis: review of epidemiology, diagnosis, management and prevention

Bordetella pertussis--the cause of pertussis or whooping cough--is an exclusively human pathogen. Disease elimination by vaccination should, therefore, be possible, but has proved elusive. Many industrialised countries with long established immunisation programs are currently seeing a resurgence of pertussis, despite universal vaccination with high uptake, with the highest burden in the least immunised age groups (infants under 6 months of age and persons over 10 years old). However, low recognition and reporting and insensitive diagnostic tests mean that the true burden of pertussis is still underestimated. Recently, efforts to improve diagnostic yield include the expanded use of polymerase chain reaction and serological tests but both have significant limitations. The range of antibiotics available for treatment and prophylaxis has expanded to include the newer macrolides, azithromycin and clarithromycin, and a range of universal and targeted vaccination strategies have been implemented or proposed. This paper reviews the current epidemiology of pertussis in developed countries, including modes of clinical presentation, diagnosis, management and potential vaccination strategies.

Genome-wide analysis of single nucleotide polymorphisms in Bordetella pertussis using comparative genomic sequencing

Bordetella pertussis is known to be a genotypically homogeneous pathogen but the extent of homogeneity at the genomic level is unknown. A currently circulating B. pertussis isolate from Australia was compared with the genome-sequenced Tohama I strain isolated in Japan in the 1950s from a distantly related lineage. Microarray-based comparative genome sequencing (CGS) was used to detect single nucleotide polymorphisms (SNPs) in a total of 1.4 Mb of the 4.09 Mb genome, including 1012 coding-regions, 217 pseudogenes and 268 intergenic regions. The CGS analysis, followed by validation using real-time PCR and DNA sequencing, identified 70 SNPs and five 1-3 bp indels, giving an overall frequency of base changes of 1 per 20 kb. Thirty-two of the 56 SNPs in coding regions were non-synonymous, including five located in virulence-associated genes. The data also allowed us to compare genomic diversity with other "clonal" human pathogens such as Mycobacterium tuberculosis and Yersinia pestis, showing that B. pertussis may be one of the least variable pathogenic bacterial species.

Comparative genomic profiling of Dutch clinical Bordetella pertussis isolates using DNA microarrays: identification of genes absent from epidemic strains

Background

Whooping cough caused by Bordetella pertussis in humans, is re-emerging in many countries despite vaccination. Several studies have shown that significant shifts have occurred in the B. pertussis population resulting in antigenic divergence between vaccine strains and circulating strains and suggesting pathogen adaptation. In the Netherlands, the resurgence of pertussis is associated with the rise of B. pertussis strains with an altered promoter region for pertussis toxin (ptxP3).

Results

We used Multi-Locus Sequence Typing (MLST), Multiple-Locus Variable Number of Tandem Repeat Analysis (MLVA) and microarray-based comparative genomic hybridization (CGH) to characterize the ptxP3 strains associated with the Dutch epidemic. For CGH analysis, we developed an oligonucleotide (70-mers) microarray consisting of 3,581 oligonucleotides representing 94% of the gene repertoire of the B. pertussis strain Tohama I. Nine different MLST profiles and 38 different MLVA types were found in the period 1993 to 2004. Forty-three Dutch clinical isolates were analyzed with CGH, 98 genes were found to be absent in at least one of the B. pertussis strains tested, these genes were clustered in 8 distinct regions of difference.

Conclusion

The presented MLST, MLVA and CGH-analysis identified distinctive characteristics of ptxP3 B. pertussis strains -the most prominent of which was a genomic deletion removing about 23,000 bp. We propose a model for the emergence of ptxP3 strains.

Genomic content of Bordetella pertussis clinical isolates circulating in areas of intensive children vaccination

Background

The objective of the study was to analyse the evolution of Bordetella pertussis population and the influence of herd immunity in different areas of the world where newborns and infants are highly vaccinated.

Methodology

The analysis was performed using DNA microarray on 15 isolates, PCR on 111 isolates as well as GS-FLX sequencing technology on 3 isolates and the B. pertussis reference strain, Tohama I.

Principal Findings

Our analyses demonstrate that the current circulating isolates are continuing to lose genetic material as compared to isolates circulating during the pre-vaccine era whatever the area of the world considered. The lost genetic material does not seem to be important for virulence. Our study confirms that the use of whole cell vaccines has led to the control of isolates that were similar to vaccine strains. GS-FLX sequencing technology shows that current isolates did not acquire any additional material when compared with vaccine strains or with isolates of the pre-vaccine era and that the sequenced strain Tohama I is not representative of the isolates. Furthermore, this technology allowed us to observe that the number of Insertion Sequence elements contained in the genome of the isolates is temporally increasing or varying between isolates.

Conclusions

B. pertussis adaptation to humans is still in progress by losing genetic material via Insertion Sequence elements. Furthermore, recent isolates did not acquire any additional material when compared with vaccine strains or with isolates of the pre-vaccine era. Herd immunity, following intensive vaccination of infants and children with whole cell vaccines, has controlled isolates similar to the vaccine strains without modifying significantly the virulence of the isolates. With the replacement of whole cell vaccines by subunit vaccines, containing only few bacterial antigens targeting the virulence of the bacterium, one could hypothesize the circulation of isolates expressing less or modified vaccine antigens.

Factors contributing to pertussis resurgence

Despite extensive immunization, the disease pertussis remains one of the world’s leading causes of vaccine-preventable deaths. An estimated 50 million cases and 300,000 deaths occur every year. A resurgence of pertussis is observed in highly immunized populations. Increasing numbers of pertussis are reported in adolescents and adults who transmit bacteria to newborns and infants to whom pertussis may be a life-threatening disease. Many studies have shown that the causes for the resurgence are multiple, such as increased awareness of disease, use of better diagnostic tools, improved surveillance methods and waning vaccine-induced immunity. Recently, antigenic divergence has been found between vaccine strains and clinical isolates in many countries with high vaccination coverage. Here, we summarize these findings and discuss the factors contributing to pertussis resurgence in immunized populations.