Volume 18, Number 5—May 2012
Temporal Trends in Bordetella pertussis Populations, Denmark, 1949–2010
We used multilocus variable-number tandem repeat analysis and multiple antigen sequence typing to characterize isolates of Bordetella pertussis strains circulating in Denmark during periods with and without pertussis vaccination coverage. Our results show substantial shifts in the B. pertussis population over time and a reduction in genetic diversity. These changes might have resulted from the introduction of pertussis vaccines in Denmark and other parts of Europe. The predominant strains currently circulating in Denmark resemble those in other European countries.
Whooping cough is a vaccine-preventable disease caused by the bacterium Bordetella pertussis; however, the protection conferred by vaccination does not last throughout life. In Denmark, as in other industrialized countries with high vaccination coverage, the disease is still endemic despite ≈50 years of vaccination. In 1961, a whole-cell B. pertussis (wP) vaccine was introduced into the childhood vaccination schedule in Denmark. In January 1997, it was replaced by a monocomponent acellular B. pertussis (aP) vaccine (DiTeKiPol or DiTeKiPol/Act-Hib; Statens Serum Institut [SSI], Copenhagen, Denmark) containing hydrogen peroxide–inactivated B. pertussis toxoid as the sole pertussis antigen. This vaccine was originally described and patented by the National Institutes of Health (Bethesda, MD, USA) (1,2). In 2003, a booster dose for preschool-age children (5 years of age) was introduced (diTekiBooster or DiTeKiPol Booster; SSI). Coverage among infants in Denmark for the third dose of pertussis vaccine has traditionally been high (≈87%–91%) (3). Acellular pertussis vaccines can contain up to 5 different antigens from B. pertussis, and Denmark is the only country using a monocomponent vaccine for both primary and booster vaccination (4).
The incidences of whooping cough and related deaths in Denmark have decreased dramatically since the introduction of pertussis vaccines. During a whooping cough epidemic in the early 1950s, before vaccine was introduced, the incidence of infection was ≈11,000 per 100,000 infants (5). In 2010, however, the incidence of infection among infants 0–1 year old had dropped to 110 per 100,000 infants and the incidence among the whole population had dropped to 7 per 100,000 persons (6). In addition, the last 2 pertussis–related deaths among infants in Denmark were notified in 2010 and 2005 (7). In general, the population-wide incidence of whooping cough in Denmark has been low since implementation of the pertussis vaccination program. However, occasional epidemic peaks have occurred, and the latest 2 were in 2002 (incidence, 36 cases/100,000 persons) and 2004 (incidence, 24 cases/100,000 persons) (8).
To determine the predominant strains of B. pertussis circulating in Denmark, we characterized clinical isolates obtained during 3 periods. Period 1, comprised 1949–1961, the year vaccine was introduced; period 2 comprised 1962–1996, during which wP vaccine was used; and period 3 comprised the years after 1996, during which aP vaccine has been used. We characterized the isolates by using multilocus variable-number tandem repeat analysis (MLVA) and multiple antigen sequence typing (MAST) to partially sequence the genes encoding pertactin (prn), B. pertussis toxin S1 subunit (ptxA, also designated ptxS1), B. pertussis toxin promoter (ptxP), and tracheal colonization factor A (tcfA). Because there is no consensus about which genes should be included in MAST of B. pertussis, we designated the results obtained in Denmark as MASTdk.
B. pertussis isolates were available from the strain collection at SSI. Strains in the collection for 1950–1974 were lyophilized; for 1974–1994, isolates were either lyophilized or stored in liquid nitrogen, and after 1995, they were mostly stored at −80°C. The isolates were cultured on B. pertussis charcoal agar plates containing cephalexin and 10% horse blood (SSI). Plates were incubated at 36°C for up to 5 days.
Strains were selected to cover whooping cough patients from all over Denmark during 1 prevaccine and 2 postvaccine periods. A total of 24 isolates were selected from period 1 (1961 and earlier); 51 were selected from period 2 (1962–1996); and 185 were selected from period 3 (1997 forward). Strains were selected regardless of the patient’s sex, age, clinical complications, and vaccine history; limited or no information was known about the patients.
DNA and Primers
Approximately 1 μL (1 inoculation loop) from a pure 4-day culture of B. pertussis on agar was suspended in 200 μL of sample preparation reagent (PrepMan; Applied Biosystems, Foster City, CA, USA), boiled for 10 min, and centrifuged for 10 min at maximum speed. For PCR, 2 μL of a 1:400 dilution of lysate was used. Primers used in the study are shown in Table 1.
Six loci, originally identified by Schouls et al. (9), were included in the typing of isolates: variable-number tandem repeats (VNTRs) 1, 3a, 3b, 4, 5, and 6. VNTR3b is a duplication of VNTR3a and is present only in a subset of isolates.VNTR2 was excluded from the analysis because it did not exhibit substantial variation between isolates (9). Amplification was performed in a total volume of 20 μL; 2 μL of a 1:400 dilution of lysate was added to a mix of 1 μL of each primer (10 pmol/μL), 10 μL of HotStarTaq Master Mix Kit (QIAGEN, Hilden, Germany), and 4 μL of 5M betaine (Sigma-Aldrich Chemie, Zwijndrecht, the Netherlands) for VNTR1, 3, 4, and 5 or 6 μL of 5 M betaine for VNTR 6.
The PCR protocol we used was a modified version of protocols published by Schouls et al. (9) and Kurniawan et al. (14). PCRs were set up as monoplex reactions and amplified in similar PCR programs, except that the annealing temperature was 68°C for primers 1, 5, and 6 and 60°C for primers 3–4. Amplification was initiated by denaturation at 96°C for 15 min and followed by 40 cycles at 95°C for 20 s, 68°C or 60°C for 30 s, 72°C for 90 s, and a final extension step at 72°C for 20 min. Final PCR products of VNTR 1, 5, and 6 and VNTR 3 and 4 were diluted 1:100 and 1:200, respectively, before the fragments were separated on an ABI 3130 DNA analyzer (Applied Biosystems).
We converted the DNA analyzer–determined size of each VNTR into the number of repeat units by using a custom-made conversion table. Data were imported to BioNumerics version 6.1 software (Applied Maths, Sint-Martens-Latem, Belgium) and analyzed. Each isolate was defined by an MLVA profile containing a string of numbers representing the number of repeats at each allele in the following order: VNTR1, 3a, 3b, 4, 5, and 6. Each unique MLVA profile was assigned an MLVA type (MT) and named according to the Dutch scheme (15).
To verify the correct conversion from sequence sizes to number of repeats, we sequenced 1–3 representatives of each VNTR size. Sequencing results showed that there is a difference between the size determined by fragment analysis and the size obtained by sequencing in our setup of the ABI 3130 DNA analyzer. VNTR1 (15-bp unit) had a difference of 26 bp, VNTR3 (5-bp unit) 10 bp, VNTR4 (12 bp-unit) 3–4 bp, VNTR5 (6 bp-unit) 9–11 bp, and VNTR6 (9-bp unit) 9–12 bp (data not shown). In all cases, the sequencing size was larger than the DNA analyzer–determined size, and the difference for each VNTR was stable for all investigated sizes. Schouls et al. suggested that the inaccurate sizing obtained by the DNA analyzer resulted from the secondary structure of the PCR product (9).
Four genetic loci known to be polymorphic for B. pertussis were selected for the MASTdk analysis: ptxP (9–11), ptxA (11,12), and 2 genes encoding surface proteins, tracheal colonization factor A (tcfA) (9,11) and pertactin (prn) (11–13). Sequence typing was performed as described for the individual genes (Table 1).
MLVA profiles were clustered in the BioNumerics version 6.1 software by using a categorical coefficient and visualized by using the minimum spanning tree method. Temporal tendencies of MLVA and MAST were determined by using the statistics tool in the BioNumerics program. The genetic diversity (Simpson’s index of diversity) was calculated by using Comparing Partitions, an online tool for quantitative assessment of partition congruence (http://www.comparingpartitions.info).
We used MLVA to type 260 clinical B. pertussis isolates collected in Denmark during 1949–2010. The isolates were resolved into 40 MTs, 27 of which have been published (15). These new types derived from the entire study period. Two MTs, MT27 and MT29, were predominant among all isolates investigated, representing 47% and 19% of the isolates, respectively. We found 7 MTs in >4 isolates, and 31 were in 1–3 isolates only and were compiled into a single group, denoted minor types, that represented ≈20% of all types (Table 2). The allelic profiles of all MTs found in this study are summarized in the Table A1).
The clonal relationship between different MTs was investigated and visualized by constructing a minimum spanning tree based on the categorical clustering of MLVA profiles (Figure 1). The tree showed that the majority of isolates belonged to 1 of the dominant MTs, MT27 or MT29, or were single-locus variants to 1 of these. MT27 and MT29 differed from each other by a single locus. A smaller proportion of isolates were more divergent, showing multiple locus differences between the MTs. The more divergent MTs primarily derived from period 1, and most isolates from period 2 were MT29 or single-locus variants of this type; period 3 was dominated by MT27 or single-locus variants of this type (Figure 2; Table 2). MT29 was first detected in Denmark in 1951, and during period 2, the proportion of MT29 increased to constitute 53% of isolates during that period. During period 3, MT27 gradually replaced MT29 to become the predominant type, constituting 65% of the isolates in period 3.
We determined the genetic diversity (Simpson’s index of diversity) of MLVA for periods 1, 2, and 3 to be 0.84, 0.72, and 0.57, respectively (Table 2). These numbers reflect a tendency of decreasing genetic diversity from period 1 to 3, i.e., a change from a more even distribution of isolates among the detected MTs toward single dominant types (Figure 2). Throughout the 3 study periods, we detected a proportion of minor and new MTs, which indicated the continuous appearance of new genetic B. pertussis types, regardless of the use or type of vaccine. Our results are in line with observations from the Netherlands (9), the United Kingdom (16), and Australia (14).
We analyzed the following genes of a selection of B. pertussis isolates by using previously published sequencing methods: prn (153 isolates), ptxA (151 isolates), ptxP (148 isolates), and tcfA (153 isolates). These genes were selected because they are the most polymorphic genes found in B. pertussis, and we were particularly interested in variation in PtxA and the Ptx promoter because Denmark has used a monocomponent vaccine containing only pertussis toxoid since 1997. Before this study, 13 prn, 8 ptxA (17), 4 tcfA (11), and 14 ptxP (10,18) alleles had been identified. In our study, we detected 4 prn alleles (prn1, prn2, prn3, and prn7), 3 ptxA alleles (ptxA1, ptxA2, and ptxA4), 5 ptxP alleles (ptxP1, ptxP2, ptxP3, ptxP15, and ptxP17), and 2 tcfA alleles (tcfA2 and tcfA3).
We analyzed the frequencies of the individual alleles and their temporal trends (Table 3; Figure 3). In the case of prn (Figure 3, panel A), there was a replacement of predominant alleles from period 1 to periods 2 and 3. The prn1 allele dominated in period 1 and 2 but was rarely observed in period 3. The prn2 allele appeared in period 2 and dominated in period 3 (86%). The prn3 allele was present at a low level in period 2 and 3, and the prn7 allele was present in period 1 but was not found later.
We observed a similar replacement of dominant alleles for ptxA. Three alleles were detected in period 1: ptxA1 (13%), pxtA2 (35%), and ptxA4 (52%). In period 2, ptxA4 was not detected, ptxA2 decreased to 11%, and ptxA1 increased to 89%. In period 3, the ptxA1 allele was the only allele detected. For tcfA, the tcfA2 allele dominated throughout all 3 periods (Figure 3, panel D). The tcfA3 allele was detected in few isolates during 1969–2001. In the case of ptxP, the ptxP1 allele gradually replaced ptxP2 during period 1, and ptxP2 was not found in isolates collected after 1952. The ptxP1 allele dominated in period 2 (96%); however, in 1995 the ptxP3 allele was detected and increased to constitute 63% of all isolates in period 3. Sporadic isolates of ptxP15 and ptxP17 were detected in period 3.
Similar observations of the temporal changes of prn and ptxA alleles have been reported from other European countries, such as United Kingdom (16), Finland (19,20), France (20), and the Netherlands (11,12). The current population of B. pertussis in Europe is now dominated by the prn2 allele and almost exclusively the ptxA1 allele (16,20,21).
The recent emergence and current dominance of the ptxP3 allele has also been reported for other European countries (9,10,22,23). In the Netherlands, ptxP3 strains produced more B. pertussis toxin than did ptxP1 strains; in addition, the emergence of ptxP3 was associated with increased whooping cough notifications and a shift in disease prevalence toward older age groups (10). In Denmark, however, the emergence of ptxP3 has not been associated with increased whooping cough notification, although a shift in age distribution has been observed (6). The emergence of ptxP3 strains is unrelated to the introduction of acellular vaccines because in Denmark and the Netherlands, these strains circulated when wP vaccines were used. In the Netherlands, ptxP3 strains reached frequencies of >90% when the wP vaccine was used for the primary series (10). The changes in the B. pertussis population in Denmark appear to be independent of epidemics because isolates from the 2002 and 2004 epidemics follow the ongoing trends in the period.
We identified 14 combinations of the sequence types of prn, ptxA, ptxP, and tcfA (MASTdk); some MASTdk types appeared to be dominant and were detected in multiple isolates, whereas others were seen in few or a single isolate (Table 4). Similar to our MLVA findings, we found that the genetic diversity determined by using Simpson’s index of diversity based on the MASTdk results decreased from 0.78 in period 1 to 0.77 in period 3 to 0.56 in period 3 (Figure 4; Table 4).
The trends for B. pertussis MTs and sequence types (or MAST) are similar in Denmark and the Netherlands (9,10,23). Different vaccines have been used in these countries, which may suggest that the most important factor driving these changes is not the type of vaccine used, but the removal, by vaccination, of immunologically naive infants as a major source of B. pertussis transmission. This might, as suggested (10), have selected for strains that are more efficiently transmitted by adolescents and adults, in whom immunity has waned.
By using MLVA and sequence typing of B. pertussis isolates collected in Denmark during 1949–2010, we showed that the population of this pathogen has changed over time. In general, the B. pertussis population in Denmark has changed from having a higher genetic diversity, as measured by MLVA and sequence typing, toward the dominance of single types. The predominant MLVA- and MASTdk-identified strains currently circulating in Denmark resemble the types observed in other European countries.
The genetic diversity of the B. pertussis isolates in Denmark was highest during period 1 (1949–1961), i.e., before introduction of pertussis vaccine. Period 2, when wP vaccine was used, was dominated by few or single types, as is the case in period 3, when aP vaccine has been used. The observed genetic changes of B. pertussis could therefore be related to the introduction of vaccines. However, there is no evidence that the wP or aP vaccines used in Denmark have selected for other dominant MT or sequence types than those observed in other European countries. Also, since the aP vaccine used in Denmark contains only B. pertussis toxoid, changes in the prn alleles in period 3 must have occurred independently of that vaccine. Travel within Europe has increased substantially since the prevaccine era; thus, B. pertussis imports from neighboring countries might explain the shift in the B. pertussis populations in Denmark. This explanation is supported by the appearance of similar B. pertussis types around Europe and by the occurrence of such types in Denmark before the introduction of change in pertussis vaccines in Denmark. The true explanation for the changes in genetic diversity among B. pertussis isolates in Denmark is probably a combination of those 2 theories.
Mrs Petersen is a researcher in the Department of Microbiological Diagnostics, Statens Serum Institute. Her research interests include molecular epidemiology and evolution of B. pertussis.
We thank Han van der Heide for assistance with the setup for sequence typing and Anette Petersen and Merete Holmberg for technical assistance.
This research was conducted as part of the national surveillance program of respiratory tract infections in Denmark.
- Sekura RD, Zhang YL, Roberson R, Acton B, Trollfors B, Tolson N, Clinical, metabolic, and antibody responses of adult volunteers to an investigational vaccine composed of pertussis toxin inactivated by hydrogen peroxide. J Pediatr. 1988;113:806–13.
- Sekura RD, Fish F, Manclark CR, Meade B, Zhang YL. Pertussis toxin. Affinity purification of a new ADP-ribosyltransferase. J Biol Chem. 1983;258:14647–51.
- Statens Serum Institut. Vaccinationsdaekning [in Danish] [cited 2011 May 5]. http://www.ssi.dk/Vaccination/Vaccinationsdaekning.aspx
- Hviid A, Stellfeld M, Andersen PH, Wohlfahrt J, Melbye M. Impact of routine vaccination with a pertussis toxoid vaccine in Denmark. Vaccine. 2004;22:3530–4.
- Nielsen A, Larsen SO. Whooping cough epidemiology in Denmark prior to and after the introduction of whooping cough vaccination. Protective effect of the vaccine and herd immunity [in Danish]. Ugeskr Laeger. 1990;152:597–604.
- Dalby T, Christensen JJ. Whooping cough 2008. EPI-NEWS. 2009; 44 [cited 2011 May 5]. http://www.ssi.dk/English/News/EPI-NEWS/~/media/Indhold/EN%20-%20engelsk/EPI-NEWS/2009/pdf/EPI-NEWS%20-%202009%20-%20No%2044.ashx
- Knudsen LK, Andersen PH. Whooping cough in children <2 years. EPI-NEWS. 2011; 42–43 [cited 2011 May 5]. http://www.ssi.dk/English/News/EPI-NEWS/2011/No%2042-43%20-%202011.aspx
- Kaltoft MS, Madsen J, Jarløv JO, Jensen TG, Prag J. Laboratory diagnosed whooping cough 2002–2004. EPI-NEWS. 2005;46 [cited 2011 May 5]. http://www.ssi.dk/English/News/EPI-NEWS/~/media/Indhold/EN%20-%20engelsk/EPI-NEWS/2005/PDF/EPI-NEWS%20-%202005%20-%20No%2046.ashx
- 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–505.
- Mooi FR, van Loo IHM, van Gent M, He Q, Bart MJ, Heuvelman KJ, Bordetella pertussis strains with increased toxin production associated with pertussis resurgence. Emerg Infect Dis. 2009;15:1206–13.
- van Loo IHM, Heuvelman KJ, King AJ, Mooi FR. Multilocus sequence typing of Bordetella pertussis based on surface protein genes. J Clin Microbiol. 2002;40:1994–2001.
- 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–5.
- King AJ, Berbers G, van Oirschot HF, Hoogerhout P, Knipping K, Mooi FR. Role of the polymorphic region 1 of the Bordetella pertussis protein pertactin in immunity. Microbiology. 2001;147:2885–95.
- Kurniawan J, Maharjan RP, Chan WF, Reeves PR, Sintchenko V, Gilbert GL, Bordetella pertussis clones identified by multilocus variable-number tandem-repeat analysis. Emerg Infect Dis. 2010;16:297–300.
- National Institute for Public Health and the Environment. MLVA: Bordetella pertussis [cited 2011 Feb 2]. http://www.mlva.net/bpertussis/default.asp
- Litt DJ, Neal SE, Fry NK. Changes in genetic diversity of the Bordetella pertussis population in the United Kingdom between 1920 and 2006 reflect vaccination coverage and emergence of a single dominant clonal type. J Clin Microbiol. 2009;47:680–8.
- Mooi FR. Bordetella pertussis and vaccination: the persistence of a genetically monomorphic pathogen. Infect Genet Evol. 2010;10:36–49.
- Advani A, Gustafsson L, Ahren C, Mooi FR, Hallander HO. Appearance of Fim3 and ptxP3-Bordetella pertussis strains, in two regions of Sweden with different vaccination programs. Vaccine. 2011;29:3438–42.
- 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–4.
- Caro V, Elomaa A, Brun D, Mertsola J, He Q, Guiso N. Bordetella pertussis, Finland and France. Emerg Infect Dis. 2006;12:987–9.
- van Amersfoorth SC, Schouls LM, van der Heide HG, Advani A, Hallander HO, Bondeson K, Analysis of Bordetella pertussis populations in European countries with different vaccination policies. J Clin Microbiol. 2005;43:2837–43.
- Advani A, van der Heide HG, Hallander HO, Mooi FR. Analysis of Swedish Bordetella pertussis isolates with three typing methods: characterization of an epidemic lineage. J Microbiol Methods. 2009;78:297–301.
- King AJ, van Gorkom T, Pennings JL, van der Heide HG, He Q, Diavatopoulos D, Comparative genomic profiling of Dutch clinical Bordetella pertussis isolates using DNA microarrays: identification of genes absent from epidemic strains. BMC Genomics. 2008;9:311.