Comparative genomic profiling of Dutch clinical Bordetella pertussis isolates using DNA microarrays: Identification of genes absent from epidemic strains
© King et al; licensee BioMed Central Ltd. 2008
Received: 29 January 2008
Accepted: 30 June 2008
Published: 30 June 2008
The Erratum to this article has been published in BMC Genomics 2010 11:196
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).
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.
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.
Bordetella bronchiseptica, Bordetella parapertussis and Bordetella pertussis are closely related respiratory pathogens that infect mammalian species.B. bronchiseptica causes chronic respiratory tract infections in a wide variety of mammals but has only been sporadically isolated from humans [1, 2]. B. parapertussis consists of two distinct lineages that infect sheep and humans, respectively . B. pertussis is is a human pathogen that causes pertussis, also known as whooping cough, a disease that is particularly severe in infants.
For about 50 years, many countries have been immunizing young children with pertussis vaccines to control disease. Despite vaccination, pertussis remains endemic, with regular epidemic outbreaks. World-wide, whooping cough causes up to 300,000 deaths every year, mostly among unvaccinated infants . Even in countries with a high vaccination coverage, a significant increase in the pertussis incidence has been observed since the 1990s [5–7]. In the Netherlands, such an increase in pertussis has been observed since 1996 .
The resurgence of pertussis in countries with high vaccination coverage has been attributed to a number of factors, including increased awareness with regard to the disease, improved diagnostics, waning vaccine-induced immunity and pathogen adaptation [4, 7]. Consistent with pathogen adaptation, numerous studies have demonstrated that the B. pertussis population has changed in several countries where vaccination has been implemented since the 1950s [9–13]. In particular, antigenic divergence was found between circulating strains and vaccine strains with respect to pertussis toxin (Ptx) and pertactin (Prn).
In the Netherlands this divergence between vaccine and circulating strains has played a role in the reemergence of pertussis [7, 14]. More recent studies have also found polymorphism in other surface proteins, including tracheal colonization factor A (TcfA) [15, 16], the serotype 2 and 3 fimbrial subunits (Fim2 and Fim3) [15, 17–19] and in the promoter region for the ptx operon . The latter mutation was found in the Netherlands since the 1990s and we observed that strains with a particular allele of the ptx promoter (ptxP), i.e. ptxP3, have expanded in the Dutch B. pertussis population. The increase in frequency of ptxP3 strains coincided with the increase in pertussis notifications in the Netherlands. Moreover B. pertussis strains carrying this novel allele for the pertussis toxin promoter were shown to confer increased virulence .
In recent years genomic tools like DNA microarrays have been used for comparison of the genetic composition of different strains within a species. Several microarray-platforms have been developed and applied to specifically address questions related to the genome of B. pertussis [22–25].
In the study presented here we used Multi-Locus Sequence Typing (MLST), Multiple-Locus Variable Number of Tandem repeat Analysis (MLVA) and a microarray-based Comparative Genomic Hybridization (CGH) method to investigate the heterogeneity of the strains dominating the current Dutch pertussis epidemics, in particularly the ptxP3 strains. We have developed an oligonucleotide microarray representing 94% of the gene repertoire of the B. pertussis Tohama I strain, of which the genome has been sequenced . This microarray allowed us to study the gene contents of the strains involved in the Dutch epidemics and thus to identify (other) polymorphic loci associated with epidemic phenotype. Representative Dutch isolates from 1993 to 2004 were selected since pertussis epidemics have been occurring regularly in the Netherlands since 1996. We aimed to identify MLST, MLVA and gene profiles associated with the Dutch epidemics and to determine the characteristics and heterogeneity of those epidemic strains.
Dutch strains were selected from the period 1993–2004. The period 1993–1995 was characterized by a relative low level of notifications. In 1996 a sudden increase in notifications was observed and this high level has been maintained until today .
Three genes known to be polymorphic in the period 1993–2004, ptxP, fim3 and prn, were selected for the MLST analysis of 158 strains. The ptxP gene occurred as two alleles, ptxP1 and ptxP3. The fim3 gene occurred as four alleles (fim3-1, fim3-2, fim3-3, fim3-4), while three prn alleles were observed (prn1,prn2 and prn3).
In contrast to ptxP1, ptxP3 was linked to two fim3 alleles, fim3-1 and fim3-2 (linkage, respectively, 52% and 48%). As was observed for ptxP1, ptxP3 was associated with prn1, prn2 and prn3 (linkage 2%, 96%, 2%). Thus in contrast to ptxP1, ptxP3 was mostly linked to prn2. Two predominant MLSTs associated with ptxP3 were observed, MLST312 and MLST322. MLST312 was found throughout the whole period and predominated in 1998 and 2004 (frequencies 44% and 58%, respectively). MLST322 was first detected in 1997 (frequency 6%) and predominated in the period 2000 to 2003 (frequencies 50% to 85%).
Presence of regions of difference in Dutch clinical isolates from 1993 to 2004
Region of difference Locus:
No. of genes
0% (nr strains)
0% (nr strains)
0% (nr strains)
0% (nr strains)
94% (nr strains)
100% (nr strains)
82% (nr strains)
100% (nr strains)
100% (nr strains)
0% (nr strains)
94% (nr strains)
100% (nr strains)
94% (nr strains)
100% (nr strains)
94% (nr strains)
92% (nr strains)
Classification of Genes in RDs.
Tot in RD
% in RD
Amino acid biosynthesis
Biosynthesis of cofactors, prosthetic groups, and carriers
Central Intermediary metabolism
Mobile and extrachromosomal element functions
Transport and binding
Total nr of genes in RD
The ptxP1 and ptxP3 lineages were distinguished by RD-10, a region comprising 18 genes and 22.7 kb (Table 2, and see Additional file 5). RD-10 was present in all ptxP1 strains analyzed. The ptxP3 allele is linked to two fim3 alleles, fim3-1 and fim3-2, and RD-10 is absent in all strains belonging to these lineages. Hierarchical clustering based on CGH data using Pearson correlation grouped the ptxP1 and the ptxP3 lineages separately (not shown). Of the 18 genes in RD-10 (see Additional file 5), four (21%) genes were classified in energy metabolism, four (21%) genes were involved in transport and binding. Further several single genes (5%) coded for a putative exported protein, a transcriptional regulator and a protein involved in iron transport, respectively (Table 2). Five of the 18 genes (26%) were pseudo genes in Tohama I. RD-10 is flanked by IS481 elements on both sides in the Tohama I strain.
Four genes present in RD-10 (BP1955-BP1958) were duplicated in the Tohama I genome at positions BP0579-BP0582 . In agreement with this, the log ratio values for these genes in strains missing RD-10, were intermediate between those for genes which were absent or present as a single copy, respectively. In contrast, identical log ratios were observed between Tohama I and ptxP1 strains indicating that the genes encoded by BP0579-BP0582 are also present at positions BP1955-BP1958 in ptxP1 strains.
B. pertussis strains carrying the ptxP3 allele were isolated for the first time in the Netherlands in 1988. In order to analyze the homogeneity of the ptxP3 strains further, we extended our microarray based CGH analysis to six strains carrying the ptxP3 allele isolated prior to 1993 (data not shown). As was observed for strains from the period 1993 to 2004, RD-10 was missing in all 6 strains. An additional deletion of two genes (RD-1*, BP0515-BP0516) was observed in four of the early ptxP3 strains. The genes in BP0515 and BP0516 encode for a phage-related hypothetical protein and a hypothetical protein respectively (see Additional file 6).
The presence of RD-10 was investigated with PCR in, respectively, 20 and 35 additional ptxP1 and ptxP3 strains. This set of strains included isolates from Denmark, Finland, France, Sweden and the Netherlands. In all cases the absence of RD-10 was linked to the ptxP3 allele.
Interestingly, in strains missing RD-10, the PCR spanning the deletion resulted in two distinct PCR fragments (sizes approximately 1.9 and 2.9 kb). DNA sequencing showed that in the 1.9 kb PCR fragment the genes flanking RD-10, BP1946 and BP1969, were separated by one copy of IS481, while in the 2.9 kb PCR fragment, two copies of the IS481 separated these genes (data not shown). The second type was found only in strains isolated in 1993–1997.
A number of RDs were found infrequently in 11 strains analyzed with CGH and were not linked to a particular lineage as defined by MLST or MLVA; RD-6 (19 genes, 18.7 kb), RD-27 (1 gene, 0.8 kb), RD-28 (2 genes, 1.7 kb), RD-29 (17 genes, 16.9 kb) and RD-18 (9 genes, 9.4 kb) (See also Additional files 7, 8, 9, 10 and 11). Sequencing analysis showed that only part (0. 2 kb) of the RD-27 gene was deleted. The five RDs contained 48 genes of which four (8%) coded for proteins involved in transport and binding, five (10%) for cell envelope proteins and four (8%) for regulatory proteins. Seven (15%) of these 48 genes are pseudo genes (Table 2). All loci were flanked by IS481 sequences on at least one side in the Tohama I strain. RD-6, RD-28 and RD-29, were each found in a single strain belonging to the ptxP1 lineage. RD-27 was found in 6 strains of the ptxP1 lineage. Finally, RD-18 was found in both the ptxP1 lineage and the ptxP3 lineage in, respectively, 1 and 4 strains. The number of pseudo genes found in the RDs (18%) was higher compared to the whole genome (10%) .
In the Netherlands the upsurge of pertussis was characterized by the expansion of strains carrying a novel allele for the Ptx promoter (ptxP3), which completely replaced the resident ptxP1 strains [20, 21]. It is possible that the ptxP3 allele has increased bacterial fitness and contributed to the upsurge of pertussis. However, the ptxP3 allele may also be linked to other polymorphic loci important for fitness. In this study we have characterized the ptxP1 and ptxP3 lineages using MLVA, MLST and CGH, to establish genetic relationships and to identify loci uniquely associated with ptxP3 strains. Further, this study highlights the microevolution of B. pertussis within closely related lineages.
The genetic diversity, based on MLVA frequencies, of the B. pertussis population showed a gradual decrease in the period 1993 to 2004 from 80% to 40%. This reflects the expansion of ptxP3-fim3-2 strains, the population of which was less diverse than the ptxP1-fim3-1 and ptxP3-fim3-1 populations (genetic diversities, respectively, 44%, 63% and 71%). The low genetic diversity of the ptxP3-fim3-2 population may reflect recent the expansion of a single clone, consistent with its recent appearance in the Dutch population.
The emergence of fim3-2 has also been observed in Canada and in Russia [18, 19]. In view of the many surface antigens expressed by B. pertussis, it is remarkable that a single amino acid change in fim3 may affect the competitive balance between strains. The alleles for the Ptx S1 subunit and Prn code for protein variants which differ in one to seven amino acids [9, 27], yet large changes in allele frequencies have been observed over the years, also suggesting that minor mutations may have a noticeable effect on strain fitness.
Consistent with studies from other groups [23, 24, 24, 25], we observed a number for RDs comprising large deletions in B. pertussis isolates which were flanked by one or more IS481 copies. A comparison of these studies is complicated by the different designations used for the RDs. We propose a common nomenclature for the RDs observed in B. pertussis according to the article by Brinig et al, . In the latter work, 26 RDs were described and we propose that novel RDs should be numbered starting from RD-27. Here RDs were numbered according to Brinig et al .
Most frequent were deletions in RD-3, RD-5 and RD-10. RD-3 and RD-5 were absent from all analyzed strains isolated in the period 1993 to 2004. RD-5 contains several genes that are predicted to be involved in iron uptake, which is essential for survival of most bacteria in the host. However, the B. pertussis genome contains several distinct gene clusters involved in iron uptake not all of which may be essential . The absence of RD-10 was associated with ptxP3 strains. RD-10 is comprised of 18 genes, four of which are duplicated and found elsewhere in the Tohama I chromosome: maleate cis-trans isomerase, a probable hydrolase (a pseudo gene in Tohama I), a conserved hypothetical protein and a putative isochorismatase. CGH and PCR analyses indicated that these four genes are duplicated in all strains analyzed belonging to the ptxP1 lineage. The duplication of part of RD-10 may have allowed the preservation of genes important for fitness. Among the RD-10 genes lost in the ptxP3 lineage are those encoding for an ABC transporter system an iron uptake protein a putative transcriptional regulator and a putative exported protein. The deleted segment did not harbor any known virulence-related genes, nevertheless some of these genes may be involved in virulence. ABC transporter components can be surface associated and may be involved in virulence .
In the ptxP3 strains the region flanking the RD-10 deletion occurred in two forms. In one form the DNA regions flanking the deletion were separated by one IS481 copy and in the second form by two IS481 copies. The form with two IS481 copies was found mostly in strains isolated in 1993–1997 but not in strains isolated after that date. All strains isolated after 1998 were found to have the one copy of IS481 in this region. This suggests that the form with two copies of the IS481 is the precursor form, from which the second form arose after a second recombination step.
A number of RDs were deleted from one to six strains. The RD-18 deletion was observed in two strains characterized by the allele combinations ptxP1-fim3-1 and ptxP3-fim3-2, respectively, and these deletions therefore may reflect independent genetic events, suggesting that deleting RD-18 is beneficial for the bacterium.
As yet it is not clear if and how these deletions affect fitness. As observed by others [23–25], we found the deleted regions to be enriched for pseudo genes compared to the genome, suggesting (part of) the deleted regions were not essential. Further, the loss of (functional) genes may increase fitness by lowering metabolic costs and the number of immune targets. In several bacterial species such as Shigella, Chlamydophila, Mycobacterium tuberculosis, Yersinia pestis and Salmonella enterica, the absence of certain genes has been correlated with some beneficial effect such as an increase in pathogenicity or the onset of virulence .
Brinig et al. , also analyzed Dutch ptxP1 strains with CGH and did not observe deletion of RD-10 in these strains. This study did not include any Dutch ptxP3 strains. Strains in which RD-10 was deleted were found in Australia, Italy and the USA . Similarly, Caro et al.  and Heikkinen et al.  described the loss of RD-10 in B. pertussis isolates in French and Finnish strains (designated RD-4 and locus 3, respectively). In the latter paper the loss of RD-10 was found in all strains of a certain PFGE profile (SR11) that is associated with nationwide epidemics. It would be interesting to determine if in other countries RD-10 loss is linked to ptxP3 as shown in our study.
B. pertussis isolates were selected from the strain collection of the Dutch National Institute of Public Health. For this study we focused on B. pertussis isolates isolated between 1993 and 2004. The selection criteria included the location of isolation and serotype. A total of 222 strains were used for MLVA-typing , 133 for MLST typing [9, 13], and 43 for CGH analysis (See Additional file 12). The strains selected for CGH analysis were all typed with MLVA and MLST and were representative of the most prevalent MLVA- and MLST-types from 1993 to 2004. This work does not require approval of the ethical commission.
MLVA- typing and MLST-typing
MLVA typing was essentially performed as described previously in Schouls et al, 2004). For MLST, gene fragments of the pertactin gene (prn), the pertussis toxin promoter (ptxP) and the fimbriae 3 gene (fim3) were sequenced essentially as described previously [15, 18, 20]. For the MLST designation the different fim3 alleles were numbered as follows: fim3A = fim3-1, fim3B = fim3-2, fim3C = fim3-3 and fim3A* = fim3-4. For each strain, every unique sequence of the ptxP, fim3 and prn loci received a distinct allele number. The compositions of these loci were expressed in an allelic profile based on the allele numbers. E.g. ptxP3-fim3-1-prn2, was abbreviated to 312. The MLVA and MLST profiles were clustered with the Bionumerics software by using a categorical coefficient and a graphing method called minimum spanning tree .
Culture of strains and DNA isolation
The B. pertussis strains used for microarray analysis in this study are listed in Additional file 12. Strains were cultured for 72 hours on Bordet-Gengou agar plates at 35°C. Subsequently, they were grown at 35°C in Verweij medium (NVI, Netherlands) with 200 μg ml-1 heptakis (2,6-di-o-methyl)-β-cyclodextrine for 24 hours while being shaken at 200 rpm. Chromosomal DNA was isolated using the Promega Wizard® Genomic DNA Purification Kit (Promega, Madison, USA) according to the manufacturer's instructions. The precipitated DNA was dissolved in 100 μl of EB (elution buffer, 10 mM Tris, pH 8.0) (Qiagen, Hilden, Germany).
B. pertussis DNA array design and construction
70-mer oligonucleotides were selected from the complete genome sequence of B. pertussis strain Tohama I in collaboration with OPERON by using microarray probe design methodology as described by OPERON . Each oligonucleotide was designed with an optimal specificity and sensitivity to its target gene and has an optimized Tm value. A total of 3751 unique 70-mer oligonucleotides, representing 94% of the genes of a B. pertussis Tohama I strain, were chosen to form the microarray. In addition, 97 oligonucleotides derived from intergenic sequences and 34 control oligonucleotides were included in the microarray. All oligonucleotides were dissolved in 50% dimethyl sulfoxide and spotted on UltraGAPS coated slides in three replicates (Corning, NY, USA) by ServiceXS (Leiden, Netherlands).
Labeling of genomic DNA
For each CGH hybridization we labeled the DNA of a test strain and the reference strain ATCC BAA-589 Tohama I. We mixed 4 μg chromosomal DNA of B. pertussis with 20 μl of the 2.5 × Random Primer Mix (BioPrime DNA labeling kit; Invitrogen) in a total volume of 41 μl of water, boiled the mixture for 5 minutes, and then placed it on ice. The samples were centrifuged for 2 minutes at 13,200 rpm and were mixed with 5 μl 10× dNTP mix [2 mM dATP, dCTP and dGTP, 0.5 mM dTTP (Roche)], 2.5 μl of 1 mM Cy3 dUTP (for reference strain) or Cy5 dUTP (for the test strain) (Amersham Biosciences, UK) and 1 μl Klenow polymerase (40 U μl-1) (BioPrime DNA labeling kit; Invitrogen). The samples were incubated for 3 hours at 37°C. Subsequently, the test and reference samples were purified separately using CyScribe GFX Purification Kit (Amersham Biosciences, UK) and then eluted in 60 μl elution buffer. The incorporation of the Cy dyes in the labeled target sequences was measured with a NanoDrop spectrophotometer (NanoDrop Technologies). For each CGH experiment equal amounts of Cy3 dye (reference) and Cy5 dye (test) were used. The volume of the samples was decreased using a speedvac concentrator (New Brunswick Scientific, Edison, USA). The labeled test and reference samples were mixed together with 79.2 μl of hybridization solution [3.44 × SSC (Invitrogen), 0.32 % SDS (Invitrogen), 1.0 mg yeast tRNA/ml]. Before loading on the microarray the hybridization solution was heated for 3 minutes at 100°C and 5.8 μl of 10 × DIG blocking buffer was added.
Before hybridization of the microarrays, 600 mJ of UV energy was applied to the Corning microarrays using a UV-chamber (Bio-Rad, Richmond, CA, USA) and the microarrays were prehybridized for 45–60 minutes at 42°C in pre-warmed prehybridization solution (5× SSC, 0.1% SDS and 0.1 mgml-1 BSA). The prehybridized microarrays were washed twice for 5 minutes with 0.1× SSC at room temperature. Finally the microarrays were washed with purified water for 30 seconds and were dried by blowing air with a Quick-Dry-filtered Air Gun (Matrix Technologies Corporation, Hudson, NH, USA). The labeled DNA hybridization mixture was applied to the microarray using lifterslip coverslips and hybridization was carried out in a Genemachines Hybchamber (Genomic Solutions, Ann Arbor, MI, USA) for 14–18 hours at 60°C in a water bath. The arrays were then disassembled in 2× SSC, 0.1% SDS at 60°C and washed for 5 minutes with 2× SSC, 0.1% SDS, in a hybridization oven at 60°C, followed by two washes 0.1× SSC, 0.1% SDS for 5 minutes at room temperature and four washes with 0.1× SSC at room temperature for one minute. Finally the slides were washed for 10 seconds in 0.01× SSC and then dried using a Quick-Dry-filtered Air Gun (Matrix Technologies Corporation, Hudson, NH, USA). Slides were scanned using a ScanArray 4000XL microarray scanner (Packard BioChip Technologies/Perkin Elmer) equipped with ScanArray express software.
Eighteen and 25 strains analyzed were hybridized at 42°C and 60°C, respectively. Essentially identical results were obtained, but the higher wash temperature resulted in a lower background and these results are shown in Fig. 5. All subsequent washes were performed according to the manufacturer's instruction (Corning, NY, USA).
Image and Data analysis
The images were analyzed using Genepix Pro 5.1 software (Axon Instruments). Two-color array image data were submitted to an in-house microarray database. Raw microarray signal data (no background subtracted) were collected from the database, and after exclusion of control and blank spots the data were normalized in R . The normalization employed a four-step approach consisting of: (1) log2-transformation, (2) quantile normalization of all scans, (3) calculating log2 (Cy5/Cy3) ratios per spot, and (4) taking the median from replicate spots. Visualizing the normalized log2-ratios as density plots and quantile-quantile plots showed an approximately normal distribution with an average of 0.012 and a standard deviation of 0.315. Based on this normal distribution, setting a cutoff for deletion at -1 would predict deleted genes with a False Discovery Rate (FDR, ) of 7.7 %; setting the cutoff at -1.5 corresponds to a FDR of 0.014%. On the basis of preliminary PCR validations, we determined a cutoff for the normalized intensity ratio of < -1.5, indicating the absence of the gene. These genes were selected for further analysis with PCR or sequence analysis. The combined normalized data were visualized with TIGR MultiExperimentViewer . The logarithm of the hybridization ratio [log2 (Cy5/Cy3)] is indicated in the yellow-black-blue color scale. [log2 (Cy5/Cy3)] = -3 = Yellow, [log2 (Cy5/Cy3)] = 0 = Black, [log2 (Cy5/Cy3)] = +3 = Blue.
PCR confirmation and sequence analysis
Primers used in this study.
Comparative Genomic Hybridization
Multi-Locus Sequence Typing
Multiple-Locus Variable number tandem repeat Analysis.
- Gueirard P, Weber C, Le Coustumier A, Guiso N: Human Bordetella bronchiseptica infection related to contact with infected animals: persistence of bacteria in host. J Clin Microbiol. 1995, 33: 2002-2006.PubMedPubMed CentralGoogle Scholar
- Diavatopoulos DA, Cummings CA, Schouls LM, Brinig MM, Relman DA, Mooi FR: Bordetella pertussis, the causative agent of whooping cough, evolved from a distinct, human-associated lineage of B. bronchiseptica. PLoS Pathog. 2005, 1: e45-10.1371/journal.ppat.0010045.PubMedPubMed CentralView ArticleGoogle Scholar
- van der Zee A, Groenendijk H, Peeters M, Mooi FR: The differentiation of Bordetella parapertussis and Bordetella bronchiseptica from humans and animals as determined by DNA polymorphism mediated by two different insertion sequence elements suggests their phylogenetic relationship. Int J Syst Bacteriol. 1996, 46: 640-647.PubMedView ArticleGoogle Scholar
- Crowcroft NS, Pebody RG: Recent developments in pertussis. Lancet. 2006, 367: 1926-1936. 10.1016/S0140-6736(06)68848-X.PubMedView ArticleGoogle Scholar
- de Melker HE, Schellekens JF, Neppelenbroek SE, Mooi FR, Rumke HC, Conyn-van Spaendonck MA: Reemergence of pertussis in the highly vaccinated population of the Netherlands: observations on surveillance data. Emerg Infect Dis. 2000, 6: 348-357.PubMedPubMed CentralView ArticleGoogle Scholar
- Guris D, Strebel PM, Bardenheier B, Brennan M, Tachdjian R, Finch E, Wharton M, Livengood JR: Changing epidemiology of pertussis in the United States: increasing reported incidence among adolescents and adults, 1990-1996. Clin Infect Dis. 1999, 28: 1230-1237. 10.1086/514776.PubMedView ArticleGoogle Scholar
- Mooi FR, Van Loo IH, King AJ: Adaptation of Bordetella pertussis to vaccination: a cause for its reemergence?. Emerg Infect Dis. 2001, 7: 526-528.PubMedPubMed CentralView ArticleGoogle Scholar
- de Melker HE, Conyn-van Spaendonck MA, Rumke HC, van Wijngaarden JK, Mooi FR, Schellekens JF: Pertussis in The Netherlands: an outbreak despite high levels of immunization with whole-cell vaccine. Emerg Infect Dis. 1997, 3: 175-178.PubMedPubMed CentralView ArticleGoogle Scholar
- Mooi FR, van Oirschot H, Heuvelman K, van der Heide HG, Gaastra W, Willems RJ: Polymorphism in the Bordetella pertussis virulence factors P.69/pertactin and pertussis toxin in The Netherlands: temporal trends and evidence for vaccine-driven evolution. Infect Immun. 1998, 66: 670-675.PubMedPubMed CentralGoogle Scholar
- Mooi FR, He Q, van Oirschot H, Mertsola J: Variation in the Bordetella pertussis virulence factors pertussis toxin and pertactin in vaccine strains and clinical isolates in Finland. Infect Immun. 1999, 67: 3133-3134.PubMedPubMed CentralGoogle Scholar
- Mastrantonio P, Spigaglia P, van Oirschot H, van der Heide HG, Heuvelman K, Stefanelli P, Mooi FR: Antigenic variants in Bordetella pertussis strains isolated from vaccinated and unvaccinated children. Microbiology. 1999, 145 ( Pt 8): 2069-2075.View ArticleGoogle Scholar
- Cassiday P, Sanden G, Heuvelman K, Mooi F, Bisgard KM, Popovic T: Polymorphism in Bordetella pertussis pertactin and pertussis toxin virulence factors in the United States, 1935-1999. J Infect Dis. 2000, 182: 1402-1408. 10.1086/315881.PubMedView ArticleGoogle Scholar
- Van Loo IH, Mooi FR: Changes in the Dutch Bordetella pertussis population in the first 20 years after the introduction of whole-cell vaccines. Microbiology. 2002, 148: 2011-2018.PubMedView ArticleGoogle Scholar
- van Boven M, de Melker HE, Schellekens JF, Kretzschmar M: A model based evaluation of the 1996-7 pertussis epidemic in The Netherlands. Epidemiol Infect. 2001, 127: 73-85. 10.1017/S0950268801005684.PubMedPubMed CentralView ArticleGoogle Scholar
- Van Loo IH, Heuvelman KJ, King AJ, Mooi FR: Multilocus sequence typing of Bordetella pertussis based on surface protein genes. J Clin Microbiol. 2002, 40: 1994-2001. 10.1128/JCM.40.6.1994-2001.2002.PubMedPubMed CentralView ArticleGoogle Scholar
- Van Amersfoorth SC, Schouls LM, van der Heide HG, Advani A, Hallander HO, Bondeson K, Von Konig CH, Riffelmann M, Vahrenholz C, Guiso N, Caro V, Njamkepo E, He Q, Mertsola J, Mooi FR: Analysis of Bordetella pertussis populations in European countries with different vaccination policies. J Clin Microbiol. 2005, 43: 2837-2843. 10.1128/JCM.43.6.2837-2843.2005.PubMedPubMed CentralView ArticleGoogle Scholar
- Packard ER, Parton R, Coote JG, Fry NK: Sequence variation and conservation in virulence-related genes of Bordetella pertussis isolates from the UK. J Med Microbiol. 2004, 53: 355-365. 10.1099/jmm.0.05515-0.PubMedView ArticleGoogle Scholar
- Tsang RS, Lau AK, Sill ML, Halperin SA, Van Caeseele P, Jamieson F, Martin IE: Polymorphisms of the fimbria fim3 gene of Bordetella pertussis strains isolated in Canada. J Clin Microbiol. 2004, 42: 5364-5367. 10.1128/JCM.42.11.5364-5367.2004.PubMedPubMed CentralView ArticleGoogle Scholar
- Borisova O, Kombarova SY, Zakharova NS, van Gent M, Aleshkin VA, Mazurova I, Mooi FR: Antigenic divergence between Bordetella pertussis clinical isolates from Moscow, Russia, and vaccine strains. Clin vaccine Immunol. 2007, 14: 234-238. 10.1128/CVI.00294-06.PubMedPubMed CentralView ArticleGoogle Scholar
- Schouls LM, van der Heide HG, Vauterin L, Vauterin P, Mooi FR: Multiple-locus variable-number tandem repeat analysis of Dutch Bordetella pertussis strains reveals rapid genetic changes with clonal expansion during the late 1990s. J Bacteriol. 2004, 186: 5496-5505. 10.1128/JB.186.16.5496-5505.2004.PubMedPubMed CentralView ArticleGoogle Scholar
- Mooi FR, van Loo IH, He Q, van Gent M, Heuvelman CJ, de Greeff SC, Diavatopoulos DA, Teunis PF, Nagelkerke NJ, Mertsola J: Bordetella Pertussis Strains With a Novel Allele for the Pertussis Toxin Promoter Are Associated With Increased Virulence and the Resurgence of Whooping Cough. 2008,Google Scholar
- Cummings CA, Brinig MM, Lepp PW, van de PS, Relman DA: Bordetella species are distinguished by patterns of substantial gene loss and host adaptation. J Bacteriol. 2004, 186: 1484-1492. 10.1128/JB.186.5.1484-1492.2004.PubMedPubMed CentralView ArticleGoogle Scholar
- Brinig MM, Cummings CA, Sanden GN, Stefanelli P, Lawrence A, Relman DA: Significant gene order and expression differences in Bordetella pertussis despite limited gene content variation. J Bacteriol. 2006, 188: 2375-2382. 10.1128/JB.188.7.2375-2382.2006.PubMedPubMed CentralView ArticleGoogle Scholar
- Caro V, Hot D, Guigon G, Hubans C, Arrive M, Soubigou G, Renauld-Mongenie G, Antoine R, Locht C, Lemoine Y, Guiso N: Temporal analysis of French Bordetella pertussis isolates by comparative whole-genome hybridization. Microbes Infect. 2006, 8: 2228-2235. 10.1016/j.micinf.2006.04.014.PubMedView ArticleGoogle Scholar
- Heikkinen E, Kallonen T, Saarinen L, Sara R, King AJ, Mooi FR, Soini JT, Mertsola J, He Q: Comparative Genomics of Bordetella pertussis Reveals Progressive Gene Loss in Finnish Strains. PLoS ONE. 2007, 2: e904-10.1371/journal.pone.0000904.PubMedPubMed CentralView ArticleGoogle Scholar
- Parkhill J, Sebaihia M, Preston A, Murphy LD, Thomson N, Harris DE, Holden MT, Churcher CM, Bentley SD, Mungall KL, Cerdeno-Tarraga AM, Temple L, James K, Harris B, Quail MA, Achtman M, Atkin R, Baker S, Basham D, Bason N, Cherevach I, Chillingworth T, Collins M, Cronin A, Davis P, Doggett J, Feltwell T, Goble A, Hamlin N, Hauser H, Holroyd S, Jagels K, Leather S, Moule S, Norberczak H, O'Neil S, Ormond D, Price C, Rabbinowitsch E, Rutter S, Sanders M, Saunders D, Seeger K, Sharp S, Simmonds M, Skelton J, Squares R, Squares S, Stevens K, Unwin L, Whitehead S, Barrell BG, Maskell DJ: Comparative analysis of the genome sequences of Bordetella pertussis, Bordetella parapertussis and Bordetella bronchiseptica. Nat Genet. 2003, 35: 32-40. 10.1038/ng1227.PubMedView ArticleGoogle Scholar
- Mooi FR, Hallander H, Wirsing von Konig CH, Hoet B, Guiso N: Epidemiological typing of Bordetella pertussis isolates: recommendations for a standard methodology. Eur J Clin Microbiol Infect Dis. 2000, 19: 174-181. 10.1007/s100960050455.PubMedView ArticleGoogle Scholar
- Garmory HS, Titball RW: ATP-binding cassette transporters are targets for the development of antibacterial vaccines and therapies. Infect Immun. 2004, 72: 6757-6763. 10.1128/IAI.72.12.6757-6763.2004.PubMedPubMed CentralView ArticleGoogle Scholar
- Lawrence JG: Common themes in the genome strategies of pathogens. Curr Opin Genet Dev. 2005, 15: 584-588. 10.1016/j.gde.2005.09.007.PubMedView ArticleGoogle Scholar
- Operon web site. 2008, [http://www.operon.com/arrays]
- r-project. 2008, [http://www.r-project.org]
- Benjamini Y, Hochberg Y: Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society Series B (Methodological). 1995, 57: 289-300.Google Scholar
- Saeed AI, Sharov V, White J, Li J, Liang W, Bhagabati N, Braisted J, Klapa M, Currier T, Thiagarajan M, Sturn A, Snuffin M, Rezantsev A, Popov D, Ryltsov A, Kostukovich E, Borisovsky I, Liu Z, Vinsavich A, Trush V, Quackenbush J: TM4: a free, open-source system for microarray data management and analysis. Biotechniques. 2003, 34: 374-378.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.