Abdolreza Advani

Global Population Structure and Evolution of Bordetella pertussis and Their Relationship with Vaccination

ABSTRACT 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. 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.

Shifts of Bordetella pertussis variants in sweden from 1970 to 2003, during three periods marked by different vaccination programs

ABSTRACT The Swedish population of Bordetella pertussis strains was characterized from 1,247 isolates covering a whole-cell vaccine program up to 1979, a 17-year period without vaccination (1979 to 1996), and a period after the introduction of general vaccination among newborns with acellular pertussis vaccines (1997 to 2003). Strains were characterized by serotyping and genotyping of pertactin and ptxA and by means of pulsed-field gel electrophoresis (PFGE). With emphasis on vaccine-related markers, the vast majority of circulating strains were of nonvaccine type. There were shifts of serotype connected with shifts of vaccination program. Serotype Fim3 was most frequent during the periods with general vaccination schedules, whereas serotype Fim2 was predominant during the 17-year vaccine-free period. Pertactin 1 was predominant during the pertussis whole-cell (Pw) vaccine period but was thereafter replaced by prn2 and has not reappeared after the introduction of acellular pertussis (Pa) vaccines. ptxA (1) was predominant over all three decades. There was a significant difference in the distribution of serotypes between vaccinated and unvaccinated individuals, but not for pertactin. A few PFGE profiles were predominant over the years: BpSR25 (serotype Fim3 prn1/7) and BpSR18 (serotype Fim3 prn2) during the Pw period, BpSR1 (serotype Fim2 prn2) during the 17 years without general vaccination, and BpSR11 (serotype Fim3 prn2) after the reintroduction of general vaccination in 1996. Despite differences between the pertactin and toxin types of Pa vaccines and circulating strains, there is no evidence that there is a threat, i.e., the vaccination program so far has been effective against whooping cough, and there seems to be no impact on the effectiveness of the vaccination program from the bacterial polymorphism.

Analysis of Bordetella pertussis Populations in European Countries with Different Vaccination Policies

ABSTRACT Despite the widespread use of pertussis vaccines during the last decades, pertussis has remained an endemic disease with frequent epidemic outbreaks. Currently two types of vaccines are used: whole-cell vaccines (WCVs) and recently developed acellular vaccines (ACVs). The long-term aim of our studies is to assess the effect of different vaccination policies on the population structure of Bordetella pertussis and ultimately on the disease burden in Europe. In the present study, a total of 102 B. pertussis isolates from the period 1998 to 2001 from five European countries (Finland, Sweden, Germany, The Netherlands, and France) were characterized. The isolates were analyzed by typing based on variable number of tandem repeats (VNTR); by sequencing of polymorphic genes encoding the surface proteins pertussis toxin S1 and S3 subunits (ptxA and ptxC), pertactin (prn), and tracheal colonization factor (tcfA); and by fimbrial serotyping. The results reveal a relationship between geographic location and VNTR types, the frequency of the ptxC alleles, and serotypes. We have not observed a relationship between the strain characteristics we studied and vaccination programs. Our results provide a baseline which can be used to reveal changes in the B. pertussis population in Europe in the coming years.

Reference system for characterization of Bordetella pertussis pulsed-field gel electrophoresis profiles.

ABSTRACT Pulsed-field gel electrophoresis (PFGE) has been used as an epidemiological tool for surveillance studies of Bordetella pertussis since the early 1990s. To date there is no standardized procedure for comparison of results, and therefore it has been difficult to directly compare PFGE results between laboratories. We propose a profile-based reference system for PFGE characterization of B. pertussis strain variation and to establish traceability of B. pertussis PFGE results. We initially suggest 35 Swedish reference strains as reference material for PFGE traceability. This reference material is deposited at the Culture Collection of the University of Gothenburg, Gothenburg, Sweden. Altogether, 1,810 Swedish clinical isolates from between 1970 and 2003 were studied, together with the Swedish Pw vaccine strain, six reference strains, and two U.S. isolates. Our system provides evidence that profiles obtained by using only one enzyme, i.e., XbaI, give enough data to analyze the epidemiological relationship between them. Characterization with one enzyme is far less labor intensive, yielding results in half the time than when a two-enzyme procedure is used. Also, we can see that there is a correlation between PFGE profile and pertactin type. One common PFGE profile, BpSR11 (n = 455), showed 100% prn2 and 100% Fim3 when analyzed for pertactin type and serotype. On the other hand, strains with the same profile may express various serotypes when isolated over longer periods of time. Subculturing of the same isolate eight times or lyophilization caused no change in PFGE profile.

Strain Variation among Bordetella pertussis Isolates in Finland, Where the Whole-Cell Pertussis Vaccine Has Been Used for 50 Years

ABSTRACT Pertussis is an infectious disease of the respiratory tract caused by Bordetella pertussis. Despite the introduction of mass vaccination against pertussis in Finland in 1952, pertussis has remained an endemic disease with regular epidemics. To monitor changes in the Finnish B. pertussis population, 101 isolates selected from 1991 to 2003 and 21 isolates selected from 1953 to 1982 were studied together with two Finnish vaccine strains. The analyses included serotyping of fimbriae (Fim), genotyping of the pertussis toxin S1 subunit (ptxA) and pertactin (prn), and pulsed-field gel electrophoresis (PFGE) after digestion of B. pertussis genomic DNA with XbaI restriction enzyme. Strains isolated before 1977 were found to harbor the same ptxA as the strains used in the Finnish whole-cell pertussis vaccine, and strains isolated before 1982 harbored the same prn as the strains used in the Finnish whole-cell pertussis vaccine. All recent isolates, however, represented genotypes distinct from those of the two vaccine strains. A marked shift of predominant serotype from Fim serotype 2 (Fim2) to Fim3 has been observed since the late 1990s. Temporal changes were seen in the genome of B. pertussis by PFGE analysis. Three PFGE profiles (BpSR1, BpSR11, and BpSR147) were distinguished by their prevalence between 1991 and 2003. The yearly emergence of the three profiles was distributed periodically. Our study stresses the importance of the continuous monitoring of emerging strains of B. pertussis and the need to obtain a better understanding of the relationship of the evolution of B. pertussis in vaccinated populations.

Appearance of Fim3 and ptxP3-Bordetella pertussis strains, in two regions of Sweden with different vaccination programs

After introduction of a mono-component vaccine, containing only pertussis toxoid (PT), the incidence of pertussis was significantly higher in the Gothenburg area among children during the period October 1, 1997 until end of 2006 compared to the Rest of Sweden where a vaccine containing PT and two other pertussis antigens was used. To investigate a possible cause of this difference, the Bordetella pertussis populations in both regions were compared by determining the fimbrial serotype (Fim), the PFGE-type and the pertussis toxin promoter allele type (ptxP). Strains with the ptxP1 allele were successively replaced by ptxP3 strains producing more pertussis toxin. In Gothenburg compared to the Rest of Sweden, Fim3 and ptxP3 strains were observed earlier and reached higher frequencies in the studied period. Since ptxP3 strains have been shown to be more virulent, their higher prevalence may have contributed to the higher incidence of pertussis in the Gothenburg area. In addition we found a high degree of linkage between PFGE-profile and ptxP3. Our results highlight the importance of strain typing to gain insight into the mechanisms of immunity-associated selection of microbial subtypes and the causes of changes in incidences of infectious diseases.

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

BackgroundBordetella 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.ResultsThe 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.ConclusionsAn 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 Strains Circulating in Europe in 1999 to 2004 as Determined by Pulsed-Field Gel Electrophoresis

ABSTRACT Clinical isolates of Bordetella pertussis collected during the year 2004 (n = 153) in eight European countries, Denmark, Finland, France, Germany, The Netherlands, Poland, Sweden, and United Kingdom, were analyzed by pulsed-field gel electrophoresis (PFGE), and their PFGE profiles were compared with those of isolates collected in 1999 (n = 102). The 255 isolates produced 59 distinct PFGE profiles. Among the 153 isolates from 2004, 36 profiles were found, while within the 102 isolates from 1999, 33 profiles were detected. One PFGE profile, BpSR11, was dominant (30% to 50%) in all countries except Denmark (10%) and Poland (0%). In comparison with 1999, there was an increase in BpSR11 prevalence in Finland in 2004 from 5% to 40%, coinciding with a major incidence peak. Some other PFGE profiles seemed to be associated with limited dissemination. Poland was the only country in which the most common actual European PFGE profiles were not found. In a dendrogram analysis, all common PFGE profiles were identified within PFGE group IV, and BpSR11 clustered together with PFGE subgroup IVβ. Compared to the 1999 isolates, PFGE group V representative for pertactin variant prn3 strains had disappeared, and a new cluster was seen. In conclusion, some PFGE profiles, such as BpSR11, evidently have a higher capacity to spread, suggesting increased fitness to the present immunological environment. It is therefore of major interest to continue with surveillance programs of B. pertussis isolates, as both waning vaccine-derived immunity and strain variation may play a role in the persistence of pertussis.

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.

Titration-free 454 sequencing using Y adapters

We describe a protocol for construction and quantification of libraries for emulsion PCR (emPCR)-based sequencing platforms such as Roche 454 or Ion Torrent PGM. The protocol involves library construction using customized Y adapters, quantification using TaqMan-MGB (minor groove binder) probe–based quantitative PCR (qPCR) and calculation of an optimal template-to-bead ratio based on Poisson statistics, thereby avoiding the need for a laborious titration assay. Unlike other qPCR methods, the TaqMan-MGB probe specifically quantifies effective libraries in molar concentration and does not require specialized equipment. A single quality control step prior to emulsion PCR ensures that libraries contain no adapter dimers and have an optimal length distribution. The presented protocol takes ∼7 h to prepare eight barcoded libraries from genomic DNA into libraries that are ready to use for full-scale emPCR. It will be useful, for example, to allow analyses of precious clinical samples and amplification-free metatranscriptomics.