Skip directly to site content Skip directly to page options Skip directly to A-Z link Skip directly to A-Z link Skip directly to A-Z link
Volume 16, Number 2—February 2010

Bordetella pertussis Clones Identified by Multilocus Variable-Number Tandem-Repeat Analysis

Jacob Kurniawan, Ram P. Maharjan, Wai-Fong Chan, Peter R. Reeves, Vitali Sintchenko, Gwendolyn L. Gilbert, Frits R. Mooi, and Ruiting LanComments to Author 
Author affiliations: University of New South Wales, Sydney, New South Wales, Australia (J. Kurniawan, R.P. Maharjan, W.-F. Chan, R. Lan); University of Sydney, Sydney (P.R. Reeves, V. Sintchenko, G.L. Gilbert); Westmead Hospital, Westmead, New South Wales, Australia (V. Sintchenko, G.L. Gilbert); National Institute for Public Health and the Environment, Bilthoven, the Netherlands (F.R. Mooi).

Cite This Article


Multilocus variable-number tandem-repeat analysis (MLVA) of 316 Bordetella pertussis isolates collected over 40 years from Australia and 3 other continents identified 66 MLVA types (MTs), including 6 predominant MTs. Typing of genes encoding acellular vaccine antigens showed changes that may be vaccine driven in 2 MTs prevalent in Australia.

Despite longstanding vaccination programs, pertussis remains endemic to many industrialized countries, including Australia, Canada, Italy, Japan, the Netherlands, Switzerland, and the United States, all of which have reported recent increases in incidence (1). Although pertussis is classically a disease of infants and children, this increase has been mainly among adults and adolescents (2,3). Factors contributing to pertussis resurgence remain unclear, but possible causes are waning immunity, suboptimal vaccine coverage, improved surveillance and diagnosis, the switch from whole cell vaccine (WCV) to acellular vaccine (ACV), and adaptation of circulating Bordetella pertussis strains (49). To determine the global epidemiology of pertussis, we analyzed an international collection of B. pertussis isolates collected mainly over the past 40 years.

The Study

We used 8 variable-number tandem-repeats, including 6 from Schouls et al. (10), to develop a multiplex PCR multilocus variable-number tandem-repeats analysis (MLVA) assay (Table 1; Technical Appendix) and used it to characterize 316 B. pertussis isolates from 12 countries on 4 continents, including 208 isolates from Australia and 87 isolates representative of common pulsed field gel electrophoresis types from Canada, Japan, Finland, and the United States (complete list available from authors). The Simpson index of diversity (D) ranged from 0.02 to 0.73 per locus with a combined D of 0.911 (Table 2). The isolates were resolved into 66 MLVA types (MTs) (complete list available from authors). Thirty-five MTs were represented by single isolates, including 15 of 208 isolates from Australia and 10 of 49 isolates from Japan. Thirty-seven MTs were previously found in Europe (10,11) and 27 were novel. Fourteen MTs were found in >2 countries or regions.

Figure 1

Thumbnail of Temporal trends of predominant multilocus variable-number tandem-repeat analysis (MLVA) types in Australia. Isolates of 4 major MLVA types (MT70, MT27, MT29, and MT64) obtained in Australia were divided into 3 periods: whole cell vaccine (WCV) (before 1997), transition from WCV to acellular vaccine (ACV) (1997–1999), and ACV (2000 onward).

Figure 1. Temporal trends of predominant multilocus variable-number tandem-repeat analysis (MLVA) types in Australia. Isolates of 4 major MLVA types (MT70, MT27, MT29, and MT64) obtained in Australia were divided into 3 periods:...

The 208 isolates from Australia were grouped into 37 MTs, of which the 4 most prevalent represented 65.4% of the isolates: MT27, 13.5%, including 1 isolate from 1973 and the others from the 1990s to 2008; MT29, 21.6%, observed since 1972; MT70, 21.2%, 1996–2005, mostly since ACV introduction in 1997; and MT64, 9.1%, during 1989–2002. Prevalence trends of the 4 most common MTs were analyzed for 3 periods determined by vaccine type(s) in use: WCV (prior to 1997), the transition period of both WCV and ACV (1997–1999), and ACV only (2000 onward) (Figure 1). MT64 prevalence was steady over time. MT29 decreased while MT27 and MT70 increased. Trends in Australia for MT27 and MT29 were similar to those observed in the United Kingdom (11) and the Netherlands (10). D values were 0.86, 0.83, and 0.83 for WCV, transition, and ACV periods, respectively. This slight decrease in genetic diversity might indicate expansion of clones that are better adapted to ACV-induced immunity.

We typed 5 genes, the products of which are used in ACVs (prn, ptxA, fim2, fim3, and fhaB), using the method of Chan et al. (12) to assess the effect of the switch from WCV to ACV on prevalence of the 4 predominant MTs (MT27, MT29, MT64, and MT70) in Australia. Isolates from these MTs have the same ptxA1 and fhaB1 alleles but vary in the other 3 genes investigated (complete list available from authors). The predominant ACV used in Australia is from GlaxoSmithKline (GSK) (Research Triangle Park, NC, USA), which contains pertussis toxoid, filamentous hemagglutinin, and pertactin but no fimbriae (FIMs). The strain used for GSK ACV contains the alleles prn1, ptxA2, and fhaB1 (10,13). However, ACV from Sanofi-Aventis (Pasteur, Lyon, France), which contains FIM2+3 in addition to pertussis toxoid, pertactin, and filamentous hemagglutinin with unknown allele types, is also licensed in Australia, complicating interpretation of variation in fim genes. On the basis of their frequencies and late appearance, fim2-2 and fim3-B are not likely to be the vaccine alleles. A significant increase (p<0.005) of prn2 (36% vs. 3%), fim2-2 (34% vs. 8%), and fim3-B (24% vs. 0%) was observed in the ACV period in comparison to the WCV period.

This increase of allelic frequency is better reflected in changes in antigenic profiles. MT27 has 3 profiles (prn1, fim2-1, fim3-A; prn2, fim2-1, fim3-A; and prn2, fim2-1, fim3-B). The first profile was seen once in the WCV period, whereas the other 2 first appeared in the WCV/ACV transition period and increased in frequency in the ACV period; the third profile, which differed by 2 alleles from the first, was more frequent. The appearance of prn2 in the second profile and additional change from fim3-A to fim3-B in the third represent increases in prevalence of alleles absent from ACV. MT29 also has 3 profiles ([prn1, prn2, or prn3], fim2-1, fim3-A), which differ in prn only. Most MT29 isolates carry prn3, and the profile is prevalent in both WCV and ACV periods, with no obvious increase in non-ACV alleles.

MT70 and MT64 both have uniform allelic profiles (prn1, fim2-2, fim3-A and prn1, fim2-1, fim3-A, respectively). However, MT70 (with fim2-2, not likely to be in ACV) increased significantly over the study period while MT64 with all alleles likely to be in ACV remained steady. Overall, the frequency of MT27 and MT70, with non-ACV alleles, increased significantly (p<0.0001) and correlated with the introduction of ACV, suggesting that antigenic changes could be driven by selection pressure.

The 2 MTs predominant in Australia were also prevalent in other countries and possibly have a global distribution. MT27 (18% of isolates) was found in 8 countries and MT29 (17% isolates) in 5. However, absence of an MT in a country might result from the small samples used. MT27 and MT29 were the most common types in the Netherlands (10) and the United Kingdom (11). MT10, MT64, MT70, MT84, and MT186 were also relatively common. MT10 and MT186 were found predominantly in Japan, although each had been found elsewhere, in China (1957) and Hong Kong (2002), respectively. MT64 was predominantly from Australia with 1 isolate from Japan, and MT70 was only found in Australia. However, all of these frequent MTs (except MT186) have been observed before. MT10 was frequent in the United Kingdom in the prepertussis vaccine era, while MT70 was common during 1998–2001 (11).

Nine isolates, including Tohama I, indentified in samples of pertussis strains collected during the 1920s–1950s from 5 countries (China, France, Japan, United Kingdom, and United States) were distributed among 7 MTs: MT10, MT12, MT75, MT83, MT127, MT205, and MT206, 2 of which were also represented among recent strains: MT10, 6 isolates from Japan 1989–2007; and MT75, 1 isolate from France in 1993. The remaining 5 MTs were either unique or shared only among the 9 early isolates.

Figure 2

Thumbnail of Minimum spanning tree (MST) of multilocus variable-number tandem-repeat analysis (MLVA types for global Bordetella pertussis isolates. The MST produced in Bionumerics (Applied Maths, Kortrijk, Belgium) used categorical coefficient and the eBURST priority rule of the highest number of single-locus changes for the clustering. Each circle represents an MLVA type with the type number in the circle. Thick lines, types differing by a single MLVA locus; thin lines, double-locus variants; d

Figure 2. Minimum spanning tree (MST) of multilocus variable-number tandem-repeat analysis (MLVA types for global Bordetella pertussis isolates. The MST produced in Bionumerics (Applied Maths, Kortrijk, Belgium) used categorical coefficient and the eBURST...

MLVA data were used to construct a minimum spanning tree (MST) (Figure 2). The 66 MTs were grouped into 2 clonal complexes and 9 singletons. Most MTs (54 of 66) belong to 1 clonal complex and 3 (MT186, MT187, and MT194) belong to another. Relationships between singletons with multiple allelic differences are not robust because they can be connected to other nodes equally. Thus, the MST cannot be rooted to infer the direction of change. Two internationally predominant MTs (27 and 29) are closely related with 1 allele difference. MT10, prevalent in Japan, is also closely related to MT29, with 1 allele difference. MT29, first isolated in the prevaccine era in the United Kingdom (11), has the highest number of SLVs and was found over 4 continents, which suggests that it arose early. Because MT10 and MT27 have a high frequency of SLVs, both likely emerged quite early. MT10 was isolated as early as 1957 in China and MT27 in 1950 in the Netherlands (10). Two high-frequency MTs (MT64 and MT70) were found in Australia only recently, with few SLVs, and may have contributed to the resurgence of pertussis in Australia.


Analysis of 208 isolates from Australia and representative isolates of common pulsed-field gel electrophoresis types from Canada, Japan, Finland, and the United States identified 6 predominant MTs (clones). Two (MT27 and MT29) were distributed worldwide, while 4 (MT10, MT64, MT70 and MT186) predominated in specific countries. Several MTs have persisted over long periods, including 3 that have circulated for at least half a century. Typing of genes encoding ACV antigens showed that use of ACV may have driven antigenic changes of 2 MTs now predominant in Australia.

Mr Kurniawan is a PhD student in medical microbiology at the University of New South Wales, Sydney, Australia. His research interests include molecular epidemiology and evolution of B. pertussis.



We appreciate the generous donations of isolates by Kazunari Kamachi, Shane Byrne, Sullivan Nicolaides, John Tapsall, Ian Carter, Margret Ip, Qiushui He, Nicole Guiso, Raymond Tsang, and Lucia Tondella. We thank Sophie Octavia for technical assistance and the reviewers for helpful suggestions.

This research was supported by the National Health and Medical Research Council of Australia. J.K. was supported by an Australian Postgraduate Award.



  1. Tan  T, Trindade  E, Skowronski  D. Epidemiology of pertussis. Pediatr Infect Dis J. 2005;24:S108. DOIPubMedGoogle Scholar
  2. Guris  D, Strebel  PM, Bardenheier  B, Brennan  M, Tachdjian  R, Finch  E, Changing epidemiology of pertussis in the United States: increasing reported incidence among adolescents and adults, 1990–1996. Clin Infect Dis. 1999;28:12307. DOIPubMedGoogle Scholar
  3. Celentano  LPMDM, Massari  MD, Paramatti  DD, Salmaso  SD, Tozzi  AE. EUVAC-NET Group. Resurgence of pertussis in Europe. Pediatr Infect Dis J. 2005;24:7615. DOIPubMedGoogle Scholar
  4. Poynten  M, McIntyre  PB, Mooi  FR, Heuvelman  KJ, Gilbert  GL. Temporal trends in circulating Bordetella pertussis strains in Australia. Epidemiol Infect. 2004;132:18593. DOIPubMedGoogle Scholar
  5. Bamberger  ES, Srugo  I. What is new in pertussis? Eur J Pediatr. 2008;167:1339. DOIPubMedGoogle Scholar
  6. Mooi  FR, van Loo  IH, King  AJ. Adaptation of Bordetella pertussis to vaccination: a cause for its reemergence? Emerg Infect Dis. 2001;7:5268. DOIPubMedGoogle Scholar
  7. Mooi  FR, He  Q, Guiso  N. Phylogeny, evolution, and epidemiology of Bordetellae. In: Locht C, editor. Bordetella molecular microbiology. Wymondham (UK): Horizon Scientific Press; 2007. p. 17–45.
  8. Kodama  A, Kamachi  K, Horiuchi  Y, Konda  T, Arakawa  Y. Antigenic divergence suggested by correlation between antigenic variation and pulsed-field gel electrophoresis profiles of Bordetella pertussis isolates in Japan. J Clin Microbiol. 2004;42:54537. DOIPubMedGoogle Scholar
  9. Byrne  S, Slack  AT. Analysis of Bordetella pertussis pertactin and pertussis toxin types from Queensland, Australia, 1999–2003. BMC Infect Dis. 2006;6:53. DOIPubMedGoogle Scholar
  10. 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:5496505. DOIPubMedGoogle Scholar
  11. 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:6808. DOIPubMedGoogle Scholar
  12. Chan  WF, Maharjan  RP, Reeves  PR, Sintchenko  V, Gilbert  GL, Lan  R. Rapid and accurate typing of Bordetella pertussis targeting genes encoding acellular vaccine antigens using real time PCR and high resolution melt analysis. J Microbiol Methods. 2009;77:3269. DOIPubMedGoogle Scholar
  13. 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:283743. DOIPubMedGoogle Scholar




Cite This Article

DOI: 10.3201/eid1602.081707

Table of Contents – Volume 16, Number 2—February 2010

EID Search Options
presentation_01 Advanced Article Search – Search articles by author and/or keyword.
presentation_01 Articles by Country Search – Search articles by the topic country.
presentation_01 Article Type Search – Search articles by article type and issue.



Please use the form below to submit correspondence to the authors or contact them at the following address:

Ruiting Lan, School of Biotechnology and Biomolecular Sciences, University of New South Wales, Sydney, NSW 2052, Australia

Send To

10000 character(s) remaining.


Page created: December 10, 2010
Page updated: December 10, 2010
Page reviewed: December 10, 2010
The conclusions, findings, and opinions expressed by authors contributing to this journal do not necessarily reflect the official position of the U.S. Department of Health and Human Services, the Public Health Service, the Centers for Disease Control and Prevention, or the authors' affiliated institutions. Use of trade names is for identification only and does not imply endorsement by any of the groups named above.