Volume 22, Number 3—March 2016
Changes in Predominance of Pulsed-Field Gel Electrophoresis Profiles of Bordetella pertussis Isolates, United States, 2000–2012
Despite high coverage with Bordetella pertussis component–containing vaccines, the incidence of reported pertussis has been increasing in the United States, and notable pertussis outbreaks have occurred in recent years (1). More than 48,000 pertussis cases were reported in the United States in 2012, the highest number reported since 1955 (2). Multiple factors have likely contributed to this increase, including increased recognition of pertussis among the general population, increased diagnosis of healthcare providers, improved diagnostic testing and reporting, and waning immunity from pertussis vaccines (3–5).
Concerns over adverse reactions after receipt of vaccines containing whole-cell preparations (WCVs) of B. pertussis led to development of vaccines with less reactogenicity (6). Starting in the 1990s, vaccines with acellular pertussis components (ACVs) began replacing the use of WCVs in the United States, and in 2005 ACVs with a lower concentration of pertussis components (known as Tdap vaccines for their tetanus and diphtheria toxoids and acellular pertussis components) were recommended for the first time as a booster dose among adolescents and adults (7,8).
Researchers in several countries have used methods such as pulsed-field gel electrophoresis (PFGE), multilocus sequence typing (MLST), and multilocus variable-number tandem-repeat analysis (MLVA) to study the evolution of molecular changes in B. pertussis populations over time (9–25). Recent studies have documented genetic changes in currently circulating B. pertussis strains, including shifts in virulence-associated protein phenotypes and predominant molecular types between prevaccine and postvaccine eras. These shifts may also be contributing to the resurgence of pertussis as a result of pathogen adaptation to current pertussis vaccines (26–28). In addition, recent changes in the epidemiology of pertussis have highlighted an increasing number of cases among older children and adolescents who have been fully vaccinated with ACVs; this increase further supports a role for vaccine pressure in the reemergence of pertussis in the United States (1,3–5,7).
PFGE, which can differentiate between individual isolates, has been used to characterize US B. pertussis isolates for >15 years (19,29). Because this method uses the entire genome, PFGE typing is more discriminatory for B. pertussis than PCR-based typing methods, such as MLST and MLVA, which analyze the sequences of a few select loci (22,30). MLST and MLVA typing must be used in tandem to obtain discriminatory power similar to that of PFGE typing (13,30).
Our objective was to increase understanding of the reemergence of pertussis and characteristics of circulating B. pertussis strains in the United States. To determine the current distribution of circulating PFGE profiles and identify changes in profile distributions over time, we analyzed PFGE profiles of B. pertussis isolates collected in the United States during 2000–2012.
A total of 5,262 isolates were available for testing, collected from 32 states during 2000–2012 (Table). Isolates were identified and collected through Enhanced Pertussis Surveillance (EPS) conducted as part of the Emerging Infections Program Network, through routine state and local health department pertussis surveillance, or during localized outbreaks. Isolates from Massachusetts and Minnesota were characterized by PFGE in their respective state public health laboratories, and tagged image file format images of PFGE gels were sent to the Centers for Disease Control and Prevention (CDC) for analysis; isolates from all other states were sent directly to CDC for PFGE characterization (Table). Because running conditions differed from those of the PFGE method commonly used in Europe, images produced by each method could not be directly compared (9,11,12,16–18,23). An additional set of 5 B. pertussis isolates, representing European profiles BpSR3, BpSR5, BpSR10, BpSR11, and BpSR12, was also included in this study with which to compare US isolates with the most commonly circulating PFGE profiles in several European countries during 1998–2009 (12).
PFGE was performed by using restriction enzyme XbaI (19), based on the method developed by Gautom et al. (31) and similar to that currently used by US state health departments participating in CDC’s PulseNet (http://www.cdc.gov/pulsenet/pathogens/index.html) for the typing of foodborne pathogens (32). PFGE patterns were compared with those in a database of B. pertussis isolate profiles maintained at CDC, and profiles were assigned on the basis of bands in the 125- to 450-kb range by using BioNumerics software version 5.01 (Applied Maths, Austin, TX, USA).
To identify predominantly circulating profiles, we analyzed overall circulating profiles for each individual year of the study period. Distributions were also compared during 2 shorter periods (2000–2009 and 2010–2012) to assess differences between recent peak years in disease and the remainder of the study period. To assess how states that submitted a large proportion of isolates affected the overall study findings, we performed 2 subanalyses. First, since Massachusetts and Minnesota together contributed >50% of isolates annually during 2000–2010 (Table), we conducted a subanalysis to check for profile local bias by comparing isolates collected from these states during 2000–2010 (Massachusetts, n = 1,758; Minnesota, n = 897), with isolates collected from all other states combined (n = 1,639). A second subanalysis was conducted in which isolates obtained during 2012 from 2 states that experienced substantial statewide epidemics of pertussis, Vermont (n = 333) and Washington (n = 222), were compared with isolates obtained from the other 10 states (n = 275) in 2012 to assess the distribution of circulating PFGE profiles during a period of widespread disease. The χ2 test was used for comparison of proportions; p<0.005 was considered significant.
A dendrogram of predominant profiles was created using BioNumerics software to visualize the degree of similarity between identified PFGE patterns. Dendrogram analysis also included the 5 predominant profiles found recently circulating in Europe. Clustering was determined by using unweighted pair group method with arithmetic mean with 1% band tolerance and optimization settings.
Overall, we identified 199 distinct PFGE profiles among the study isolates. CDC013, CDC010, CDC082, CDC002, and CDC046 were the predominant 5 profiles among our study population and accounted for 72% of all isolates tested (Figure 1). One additional profile, CDC237, accounted for 5% of isolates overall but was only observed in isolates collected in 2009 and subsequently (Figure 2). None of the remaining 193 profiles accounted for >5% of isolates overall.
When we assessed the most predominant profiles by year, we noted differences in the order of predominance beginning in 2010 (Figure 2). CDC013 was consistently predominant during 2000–2009, but CDC082 was predominant among isolates collected in 2010. CDC002 appeared to be fading from circulation during 2000–2009, but it emerged as the most common profile among 2011 isolates. CDC237 and CDC002 predominated in 2012; each accounted for ≈29% of circulating profiles.
When the 5 most predominant profiles were compared in each period, we found CDC013 in 41% of isolates collected during 2000–2009 (35%–46% annually) but only in 9% of isolates analyzed during 2010–2012 (p<0.0001). Prevalence also declined significantly between periods for CDC010 (from 15% to 4.5%, p<0.0001) and CDC082 (from 12% to 3%, p<0.0001). Conversely, profile CDC002, which comprised 4% of isolates in the earlier period, increased to 25% among isolates collected during 2010–2012 (p<0.0001). No changes were observed for profile CDC046, which accounted for a consistent 6% of isolates in each period. In addition, frequency of 3 profiles (CDC217, CDC237, CDC253), each accounting for <1% of isolates during 2000–2009, increased significantly to 4%, 24%, and 6%, respectively, of isolates collected during 2010–2012 (p<0.0001 for all).
When we compared profiles of Massachusetts and Minnesota isolates to isolates of all other states combined during 2000–2010, our subanalysis revealed that CDC013 was the most common PFGE profile among all 3 groups, ranging from 37% to 44% of isolates (Figure 3, panel A). In addition, CDC010 accounted for a similar proportion of isolates across all groups (12%–17%). Whereas a similar proportion of isolates were of profile CDC082 in Minnesota and in all other states combined (9% and 10%, respectively), this profile was found in significantly higher proportion in Massachusetts (17%, p<0.0001), where it was the second most predominant profile. CDC046 accounted for a consistent 5%–7% of isolates from all 3 groups, and CDC002 frequency ranged between from 2% to 5%.
Our subanalysis of the predominant profiles in 2012 isolates from 2 statewide epidemics (Vermont and Washington) and 10 other states combined revealed marked differences between the collections (Figure 3, panel B). CDC013 predominated among Washington isolates at 22%, which was significantly more than the 1% and 4% of isolates with this profile found among isolates from Vermont and combined states, respectively (p<0.0001). In Vermont, a higher percentage of isolates was CDC237 (41%), compared with 18% and 22% of isolates from Washington and the other states combined, respectively (p<0.0001). CDC046 was also found in significantly higher proportion among Vermont isolates (11%) than among those from Washington (1%; p<0.0001), but the proportion among isolates from the other states (5%, p = 0.0098) was not significantly higher. The most predominant profile from the other states combined group was CDC002 (35%), which accounted for a similar proportion of isolates from Vermont (29%, p = 0.0966), but a significantly higher proportion of isolates from Washington (18%, p<0.0001). CDC010 comprised only 1% of isolates from Vermont, significantly fewer than those from Washington (7%, p = 0.0003) and the other states combined (6%, p = 0.0021). CDC082 accounted for <1% of isolates in each group.
A dendrogram analysis comparing the 5 most predominant profiles overall in the United States during the study period (CDC013, CDC010, CDC082, CDC002, CDC046) with 3 increasingly circulating profiles (CDC217, CDC237, CDC253) revealed 2 clusters with 75% similarity. Cluster 1 contains profiles CDC013, CDC046, and CDC082, first seen among US isolates collected during 1993–1997. Cluster 2 contains profiles CDC002, CDC010, CDC217, CDC237, and CDC253. CDC002 and CDC010 were first seen among US isolates collected in the 1980s, whereas the other 3 profiles were not seen until 2007. In the area of band analysis used in this study (125 kb–450 kb), the profiles of the 5 predominant strains from Europe, (BpSR3, BpSR10, BpSR11, BpSR5, BpSR12) were indistinguishable from the 5 predominantly circulating US profiles (Figure 4).
Overall genetic diversity among the study isolates was 84%. Genetic diversity ranged from 71% to 82% annually for isolates collected during 2000–2009 before a high of 91% was reached in 2010 (Figure 2). By 2012, diversity decreased slightly to 82%.
Of the 199 B. pertussis PFGE profiles observed in this study, a single profile, CDC013, predominated between 2000 and 2009, and 5 profiles comprised most isolates tested from 2000 to 2012. These finding are similar to those previously reported by Hardwick et al., who also found 3 predominant profiles circulating in the United States during 1986–1999 (19). Hardwick et al. reported the 3 predominant profiles as CYXXI-010 (now CDC010) (37%), CYXXI-002 (now CDC002) (11%), and CYXXI-013 (now CDC013) (10%), which are still among the currently circulating predominant profiles in the United States. Additionally, 42% of the profiles observed during 1986–1999 were still circulating in the United States during our study. Other studies have also shown that a small number of PFGE profiles usually make up most of the circulating B. pertussis strains during outbreak and nonoutbreak periods (12,29).
In contrast to the isolates collected during 2000–2009, profile predominance changed rapidly among isolates collected during 2010–2012, and a different profile predominated each year. Rare profiles such as CDC237 quickly became more common than the previously predominant CDC013 profile. Two additional rarely seen profiles, CDC217 and CDC253, which were closely related to the CDC237 profile, also increased significantly in recent years. These findings prompted us to explore possible associations between PFGE profile and other molecular changes that have recently occurred in the organism. Specifically, in recent years, pertactin-deficient B. pertussis isolates have rapidly emerged. Pertactin is a key immunogen included in all ACVs currently used in the United States (21,30), and data have suggested a possible selective advantage of pertactin-deficient mutants among ACV-vaccinated populations (27). Using supplemental data, we observed that 87% of 2012 isolates included in our analysis were pertactin-deficient. Of interest, isolates with PFGE profiles CDC002 and CDC237, the most common profiles among 2012 isolates, were significantly more likely to be pertactin-deficient than isolates with profile CDC013, the most common profile during 2000–2009. That these 2 profiles started becoming more predominant in 2010, coinciding with the rapid increase in pertactin-deficient isolates in the United States (21), suggests a linkage between the 2 changes. To date, no associations have been reported between pertactin-deficiency and specific PCR-based MLVA or MLST types, which reinforces the discrimatory power of PFGE typing methods and highlights the value of this method for elucidating B. pertussis evolution on a more granular level (24,30).
Although Massachusetts and Minnesota contributed a large proportion of isolates during 2000–2010, we observed similar distributions among 4 of the 5 predominant PFGE profiles when we compared isolates collected in these 2 states to isolates in all other states combined. These results suggest that location bias did not affect the findings. This observation contrasts with results of our subanalysis of 2012 isolates, in which significant differences were noted in the predominant strains found in Vermont, Washington, and all other states combined. This finding could be caused by many factors, including clonal expansion within epidemic areas and, as previously noted, selective advantage of pertactin-deficient B. pertussis among vaccinated populations. Exploration of epidemiologic data associated with these isolates may further explain the observed geographic differences in PFGE profiles.
Our analysis showed that the 5 most predominant profiles circulating in several countries in Europe (Denmark, Finland, France, Germany, the Netherlands, Norway, Sweden, United Kingdom) during 1998–2009 (12) were also the 5 most predominant profiles circulating in the United States during 2000–2012. The proportion of the most predominant profile in Europe, BpSR11, changed from 25%–30% of 1998–2005 isolates to 13% of 2007–2009 isolates, which was similar to changes observed in the United States equivalent CDC013 profile (41%, 2000–2009, to 9%, 2010–2012). Comparable increases were also noted for CDC002 (4%, 2000–2009, to 25%, 2010–2012) and its European equivalent, BpSR3 (0%, 1998–2005, to 21%, 2007–2009).
We did not detect a large increase in CDC010 during 2010–2012, such as was reported for BpSR10 during 2007–2009 in Europe, but instead detected a large increase in CDC237, first seen among US isolates in 2009. As shown, CDC010 and CDC0237 are extremely similar and cluster together in our dendrogram (Figure 4), which indicates that they are closely related. Possible explanations for differences in these profiles include a point mutation with the gain of a restriction site in CDC237 isolates or a rearrangement at the chromosomal level (23,35,36).
Similar to the findings for US isolates, the study in Europe also reported increased diversity (12). However, even though circulating PFGE profiles appear to be similar in the United States and Europe, percentages of pertactin-deficient isolates collected during 2010–2014 in European countries are much lower than in the United States, which suggests that even though PFGE changes and pertactin deficiency are occurring simultaneously in the United States, these changes are not necessarily linked and that virulence and transmission of pertussis depend on many factors (37,38).
A limitation of our study is that the differences in switch times and running time of our PFGE method prevents direct comparison with isolates tested by the European PFGE method (23). Our method results in low resolution of 6–8 small (45- to 125-kb) DNA bands, which are not included as part of our analysis, leading to fewer PFGE profiles overall when compared to the European method. However, by using the same PFGE method as in our previous study, we were able to directly compare 6,595 US isolates collected over 77 years. The profile changes observed in the course of the 2 studies have provided useful information for the development of current molecular typing systems and whole-genome sequencing strategies.
Our work reveals that B. pertussis PFGE profiles are changing rapidly in the United States, as they are in Europe (12). The previously held assumption thatB. pertussis was a highly clonal organism has been challenged because recent isolates of B. pertussis have been shown to be quite variable at the genome level, particularly in those genes that code for ACV components (21,27,28). Although more emphasis is now being placed on molecular typing methods that do not require the availability of isolates, PFGE, in tandem with other typing methods, continues to remain a valuable tool for understanding B. pertussis because it can detect a greater level of diversity among circulating strains.
As we move more toward whole-genome sequencing and rearrangement analysis, PFGE can be useful for strain selection when genetic diversity is a desired parameter for analysis. Having a more complete picture of the molecular evolution of B. pertussis through whole-genome methods such as PFGE and sequencing may provide key data to help guide development of the next generation of pertussis vaccines or, in the near term, inform us on ways to more effectively use existing vaccines. Continued monitoring of the molecular epidemiology of circulating B. pertussis continues to be critical for understanding the changing epidemiology of the disease both in the United States and abroad to optimize current prevention and control strategies.
Ms. Cassiday is a research microbiologist in the Meningitis and Vaccine Preventable Diseases Branch, Division of Bacterial Diseases, National Center for Immunization and Respiratory Diseases, CDC. Her research focuses on Bordetella pertussis and Corynebacterium diphtheriae.
We gratefully acknowledge the following state health laboratories, who are former or current participants in CDC’s Emerging Infections Program Enhanced Pertussis Surveillance, for contributing most of the isolates used in this study: Arizona, California, Connecticut, Georgia, Massachusetts, Minnesota, New Mexico, New York, and Oregon. We also thank Hans Hallander of the Swedish Institute for Communicable Disease Control for the donation of the European B. pertussis type strains used in this study and Marina Chesnokova for technical assistance.
- Clark TA, Messonnier NE, Hadler SC. Pertussis control: time for something new? Trends Microbiol. 2012;20:211–3.
- Centers for Disease Control and Prevention. Pertussis (whooping cough)—fast facts. 2013 [cited 2014 Jan 24]. http://www.cdc.gov/pertussis/fast-facts.html
- Centers for Disease Control and Prevention. Pertussis epidemic—Washington, 2012. MMWR Morb Mortal Wkly Rep. 2012;61:517–22 .
- Misegades LK, Winter K, Harriman K, Talarico J, Messonnier NE, Clark TA, Association of childhood pertussis with receipt of 5 doses of pertussis vaccine by time since last vaccine dose, California, 2010. JAMA. 2012;308:2126–32.
- Winter K, Harriman K, Zipprich J, Schechter R, Talarico J, Watt J, California pertussis epidemic, 2010. J Pediatr. 2012;161:1091–6.
- Decker MD, Edwards KM, Steinhoff MC, Rennels MB, Pichichero ME, Englund JA, Comparison of 13 acellular pertussis vaccines: adverse reactions.Pediatrics. 1995;96:557–66 .
- Skoff TH, Cohn AC, Clark TA, Messonnier NE, Martin SW. Early impact of the US Tdap vaccination program on pertussis trends. Arch Pediatr Adolesc Med. 2012;166:344–9 .
- Centers for Disease Control and Prevention. Pertussis vaccination: use of acellular pertussis vaccines among infants and young children. Recommendations of the Advisory Committee on Immunization Practices (ACIP). MMWR Recomm Rep. 1997;46(RR-7):1–25 .
- Advani A, Donnelly D, Gustafsson L, Hallander HO. Changes of the Swedish Bordetella pertussis population in incidence peaks during an acellular pertussis vaccine period between 1997 and 2004. APMIS. 2007;115:299–310.
- 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.
- Advani A, Gustafsson L, Carlsson RM, Donnelly D, Hallander HO. Clinical outcome of pertussis in Sweden: association with pulsed-field gel electrophoresis profiles and serotype. APMIS. 2007;115:736–42.
- Advani A, Hallander HO, Dalby T, Krogfelt KA, Guiso N, Njamkepo E, Pulsed-field gel electrophoresis analysis of Bordetella pertussis isolates circulating in Europe from 1998 to 2009. J Clin Microbiol. 2013;51:422–8.
- 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.
- 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–8.
- Elomaa A, Advani A, Donnelly D, Antila M, Mertsola J, Hallander H, Strain variation among Bordetella pertussis isolates in Finland, where the whole-cell pertussis vaccine has been used for 50 years. J Clin Microbiol. 2005;43:3681–7.
- Elomaa A, He QS, Minh NNT, Mertsola J. Pertussis before and after the introduction of acellular pertussis vaccines in Finland. Vaccine.2009;27:5443–9.
- Hallander H, Advani A, Riffelmann M, von Konig CH, Caro V, Guiso N, Bordetella pertussis strains circulating in Europe in 1999 to 2004 as determined by pulsed-field gel electrophoresis. J Clin Microbiol. 2007;45:3257–62.
- Hallander HO, Advani A, Donnelly D, Gustafsson L, Carlsson RM. Shifts of Bordetella pertussis variants in Sweden from 1970 to 2003, during three periods marked by different vaccination programs. J Clin Microbiol. 2005;43:2856–65.
- Hardwick TH, Cassiday P, Weyant RS, Bisgard KM, Sanden GN. Changes in predominance and diversity of genomic subtypes of Bordetella pertussisisolated in the United States, 1935 to 1999. Emerg Infect Dis. 2002;8:44–9.
- Schmidtke AJ, Boney KO, Martin SW, Skoff TH, Tondella ML, Tatti KM. Population diversity among Bordetella pertussis isolates, United States, 1935–2009. Emerg Infect Dis. 2012;18:1248–55.
- Pawloski LC, Queenan AM, Cassiday PK, Lynch AS, Harrison M, Shang W, Prevalence and molecular characterization of pertactin-deficientBordetella pertussis in the US. Clin Vaccine Immunol. 2014;21:119–25.
- Caro V, Njamkepo E, Van Amersfoorth SC, Mooi FR, Advani A, Hallander HO, Pulsed-field gel electrophoresis analysis of Bordetella pertussispopulations in various European countries with different vaccine policies. Microbes Infect. 2005;7:976–82.
- Advani A, Donnelly D, Hallander H. Reference system for characterization of Bordetella pertussis pulsed-field gel electrophoresis profiles. J Clin Microbiol. 2004;42:2890–7.
- Tsang RS, Shuel M, Jamieson FB, Drews S, Hoang L, Horsman G, Pertactin-negative Bordetella pertussis strains in Canada: characterization of a dozen isolates based on a survey of 224 samples collected in different parts of the country over the last 20 years. Int J Infect Dis. 2014;28:65–9.
- van Gent M, Heuvelman CJ, van der Heide HG, Hallander HO, Advani A, Guiso N, Analysis of Bordetella pertussis clinical isolates circulating in European countries during the period 1998–2012. Eur J Clin Microbiol Infect Dis. 2015;34:821–30.
- Mooi FR, Van Der Mas NAT, De Melker H. Pertussis resurgence: waning immunity and pathogen adaptation—two sides of the same coin. Epidemiol Infect. 2014;142:685–94.
- Martin SW, Pawloski L, Williams M, Weening K, DeBolt C, Qin X, Pertactin-negative Bordetella pertussis strains: evidence for a possible selective advantage. Clin Infect Dis. 2015;60:223–7.
- Bart MJ, Harris SR, Advani A, Arakawa Y, Bottero D, Bouchez V, Global population structure and evolution of Bordetella pertussis and their relationship with vaccination. MBio. 2014;5:e01074–14.
- Bisgard KM, Christie CD, Reising SF, Sanden GN, Cassiday PK, Gomersall C, Molecular epidemiology of Bordetella pertussis by pulsed-field gel electrophoresis profile: Cincinnati, 1989–1996. J Infect Dis. 2001;183:1360–7.
- Bowden KE, Williams MM, Cassiday PK, Milton A, Pawloski L, Harrison M, Molecular epidemiology of the pertussis epidemic in Washington State in 2012. J Clin Microbiol. 2014;52:3549–57.
- Gautom RK. Rapid pulsed-field gel electrophoresis protocol for typing of Escherichia coli O157:H7 and other gram-negative organisms in 1 day. J Clin Microbiol. 1997;35:2977–80 .
- Swaminathan B, Barrett TJ, Hunter SB, Tauxe RV, Force CDCPT. PulseNet: the molecular subtyping network for foodborne bacterial disease surveillance, United States. Emerg Infect Dis. 2001;7:382–9.
- Hunter PR, Gaston MA. Numerical index of the discriminatory ability of typing systems: an application of Simpson’s index of diversity. J Clin Microbiol. 1988;26:2465–6 .
- Grundmann H, Hori S, Tanner G. Determining confidence intervals when measuring genetic diversity and the discriminatory abilities of typing methods for microorganisms. J Clin Microbiol. 2001;39:4190–2.
- Tenover FC, Arbeit RD, Goering RV, Mickelsen PA, Murray BE, Persing DH, Interpreting chromosomal DNA restriction patterns produced by pulsed-field gel electrophoresis: criteria for bacterial strain typing. J Clin Microbiol. 1995;33:2233–9 .
- Stibitz S, Yang MS. Genomic plasticity in natural populations of Bordetella pertussis. J Bacteriol. 1999;181:5512–5 .
- Zeddeman A, van Gent M, Heuvelman CJ, van der Heide HG, Bart MJ, Advani A, Investigations into the emergence of pertactin-deficient Bordetella pertussis isolates in six European counries, 1996–2012. Euro Surveill. 2014;19:20881.
- Hegerle N, Guiso N. Bordetella pertussis and pertactin-deficient clinical isolates: lessons for pertussis vaccines. Expert Rev Vaccines.2014;13:1135–46.
Suggested citation for this article: Cassiday PK, Skoff TH, Jawahir S, Tondella ML. Changes in predominance of pulsed-field gel electrophoresis profiles ofBordetella pertussis isolates, United States, 2000–2012. Emerg Infect Dis. 2016 Mar [date cited]. http://dx.doi.org/10.3201/eid2203.151136