Skip to main content

Molecular epidemiology of Bordetella pertussis and analysis of vaccine antigen genes from clinical isolates from Shenzhen, China



Although pertussis cases globally have been controlled through the Expanded Programme on Immunization (EPI), the incidence of pertussis has increased significantly in recent years, with a “resurgence” of pertussis occurring in developed countries with high immunization coverage. Attracted by its fast-developing economy, the population of Shenzhen has reached 14 million and has become one of the top five largest cities by population size in China. The incidence of pertussis here was about 2.02/100,000, far exceeding that of the whole province and the whole country (both < 1/100,000). There are increasing numbers of reports demonstrating variation in Bordetella pertussis antigens and genes, which may be associated with the increased incidence. Fifty strains of Bordetella pertussis isolated from 387 suspected cases were collected in Shenzhen in 2018 for genotypic and molecular epidemiological analysis.


There were 387 suspected cases of pertussis enrolled at surveillance sites in Shenzhen from June to August 2018. Nasopharyngeal swabs from suspected pertussis cases were collected for bacterial culture and the identity of putative Bordetella pertussis isolates was confirmed by real-time PCR. The immunization history of each patient was taken. The acellular pertussis vaccine (APV) antigen genes for pertussis toxin (ptxA, ptxC), pertactin (prn) and fimbriae (fim2 and fim3) together with the pertussis toxin promoter region (ptxP) were analyzed by second-generation sequencing. Genetic and phylogenetic analysis was performed using sequences publicly available from GenBank, National Institutes of Health, Bethesda, MD, USA ( The antimicrobial susceptibility was test by Kirby-Bauer disk diffusion.


Fifty strains of Bordetella pertussis were successfully isolated from nasopharyngeal swabs of 387 suspected cases, with a positivity rate of 16.79%, including 28 males and 22 females, accounting for 56.0% and 44.0% respectively. Thirty-eight of the 50 (76%) patients were found to be positive for B. pertussis by culture. Among the positive cases with a history of vaccination, 30 of 42 (71.4%) cases had an incomplete pertussis vaccination history according to the national recommendation. Three phylogenetic groups (PG1-PG3) were identified each containing a predominant genotype. The two vaccines strains, CS and Tohama I, were distantly related to these three groups. Thirty-one out of fifty (62%) isolates belonged to genotype PG1, with the allelic profile prn2/ptxC2/ptxP3/ptxA1/fim3-1/fim2-1. Eighteen out of fifty (36%) isolates contained the A2047G mutation and were highly resistant to erythromycin, and all belonged to genotype PG3 (prn1/ptxA1/ptxP1/ptxC1/fim3-1/fim2-1), which is closely related to the recent epidemic strains found in northern China.


The positive rate of cases under one-year-old was significantly higher than that of other age groups and should be monitored. The dominant antigen genotypes of 50 Shenzhen isolates are closely related to the epidemic strains in the United States, Australia and many countries in Europe. Despite high rates of immunization with APV, epidemics of pertussis have recently occurred in these countries. Therefore, genomic analysis of circulating isolates of B. pertussis should be continued, for it will benefit the control of whooping cough and development of improved vaccines and therapeutic strategies.


Pertussis is an acute respiratory infectious disease caused by the gram-negative bacterium Bordetella pertussis [1]. The typical symptoms of pertussis include paroxysmal coughing with an inspiratory whoop, post-tussive vomiting, cyanosis, and persistent coryzal symptoms [2]. Pertussis is mainly transmitted by droplets and older children and adults have been increasingly recognized as reservoirs of infection for transmission to unvaccinated or incompletely vaccinated infants who are most at risk of severe disease complications and death [3]. The first pertussis whole-cell vaccine (WPV) made from killed cells was licensed in the USA in 1914. Later chemical inactivation was used and the WPV combined with tetanus and diphtheria toxoids in the 1940s to become widely used diphtheria-tetanus-pertussis (DTP) vaccine, which is about 80% effective in preventing severe illness and death from pertussis. Since then, a successive decline in the incidence of the disease has been observed [4]. However, due to the reactogenicity following WPV immunization (including fever, limb redness and swelling) and the incorrect public opinion linking the WPV with encephalopathies, many parents refused to vaccinate their children, vaccine coverage fell and lawsuits against the vaccine manufacturers forced many of them to stop producing the vaccine [5, 6]. To address the reactogenicity of WPV, vaccine manufacturers switched to producing APV containing a number of highly purified antigens [1]. Since then, there has been an increase in whooping cough cases and a dramatic cycle of epidemics [7, 8]. According to the World Health Organization (WHO), as many as 195,000 children worldwide died of pertussis and its complications in 2008, with 90% of cases occurring in less developed and developing countries [9]. Even in developed countries or countries with high pertussis vaccination rates, the incidence of pertussis has been on the rise in recent years [10]. In 2012, there was an outbreak of pertussis in Washington and other states in the United States, with the highest reported incidence (37.5/100,000) since 1942. It is worth noting that 43% of patients had at least four doses of immunological history of acellular diphtheria, tetanus and pertussis (DTaP) vaccine [11]. According to the different production methods, the common DTaP vaccine used in China can be divided into two types: co-purified vaccines and component vaccines. Domestic manufacturers mainly produced co-purified vaccines. Since there was no separate purification of different antigen components during production procedure, the proportion of antigens in the original solution varied according to the manufacturers and production batches. A study has pointed out that the co-purified pertussis vaccines produced in China contained not only the main antigens purified pertussis toxoid (PT), filamentous hemagglutinin (FHA) and pertactin (PRN), but also fimbriae (FIM) 2 and 3 and other minor protein antigens [12]. For component vaccines, PT, FHA, PRN and other antigens were extracted separately by column chromatography, and then were mixed in a certain proportion, which ensured the consistent quality between batches and had higher purity. Only some imported DTaP vaccines (Pentaxim, Sanofi Pasteur, France and Infanrix, GlaxoSmithKline, Britain) used this technology for production. However, the use of imported vaccines were limited by their high price [13]. Historically, the intracerebral mouse protection test (Kendrick test) has been used to effective determine the potency of whole-cell pertussis vaccines and is the only test that has shown a correlation with protection in children [14].

In recent years, the reported incidence of pertussis in adolescents and adult has increased, due to a combination of better awareness of less severe clinical symptoms and increased availability of diagnostic testing including PCR and serology [15, 16]. Vaccines are the main way to prevent pertussis, which have been widely used for immunization in China since 1960s [17]. Shenzhen adopted whole-cell vaccines before 2008, and gradually replaced it with APV from 2008 to 2010. Since 2010, APV have been used throughout the city. According to the national immunization program, the the recommended age for receipt of the DTaP is 3, 4, 5 and 18 months after birth [18]. The vaccination rate has kept above 99% in recent years. However, according to the statistics of the China Information System for Disease Control and Prevention, the number of reported cases of pertussis in Shenzhen is rising sharply. Studies have shown that pathogens have undergone adaptive changes under the pressure of immune protection pressure induced by the body after vaccination. The alleles of the genes encoding the vaccine antigens including PT, PRN, FIM in the predominant circulating strains in Shenzhen are distinct from the vaccine strains [19, 20].

To further explain the possible impact of the genetic difference on the immune effects, and to better understand the pathogenic characteristics, evolutionary characteristics, and molecular epidemiology of Bordetella pertussis in China, we sequenced 50 strains of Bordetella pertussis isolated from Shenzhen in 2018 to analyze their population structure and sequence characteristics of the major vaccine antigen genes.

Materials and methods

Research object

Three hundred and eight-seven suspected pertussis cases from the outpatient department of Shenzhen Children’s Hospital from June to August 2018.

Main reagents and instruments

Main reagents

Charcoal Agar Base and Bordetella selective supplement (Oxoid company, Canada), Pertussis Bacteria Phase I Standard Serum (REMEL company, Lenexa City, KS, USA) and Whole-genome DNA extraction kit (Qiagen Company, Shanghai, China).

Main instruments

Fluorescence PCR instrument (ABI Company, Oyster Bay, NY, USA) and NanoDrop 1000 Ultra Micro Spectrophotometer (Thermo Fisher Company, Fair Lawn, NJ, USA).


Charcoal agar was prepared according to the instructions of the manufacturer. After autoclaving and cooling to 50 °C, Bordetella selective supplement and 10% defibrinated sheep blood were added into the agar. The samples were inoculated onto the charcoal blood agar plates and incubated at 37 °C for 3 days. Based on colonial appearance, round, moist, protruding, transparent or translucent, presumptive Bordetella pertussis colonies were subcultured to obtain pure cultures and identified by biochemical reactions and agglutination tests. Genomic DNA was extracted from the presumptive B. pertusiss isolates and tested by real-time PCR using the insertion element IS481 and pertussis toxin subunit A gene targets [21].

Genome sequencing and de novo assembly

The genomic DNA of Bordetella pertussis was extracted using the whole-genome DNA extraction kits according to the manufacturer’s instructions. The library was constructed using the Illumina HiSeq 4000 system (Illumina, San Diego, CA, USA) at the Beijing Genomics Institute (Shenzhen, China). DNA sample preparation kit (Illumina, San Diego, CA, USA), and whole gene sequencing and preliminary evaluation filtering was conducted using the Illuminas HiSeq X-Ten platform (BGI, Shenzhen, China). The original data was further filtered by SOAPnuke to obtain valid data via quality control. De novo assembly of the genome was performed by SPAdes gene assembly software (V3.9.1).

Sequence analysis of major antigen genes

The sequence numbers of the reference genes are as follows:

ptxA gene (vaccine CS strain [22] ptxA2 type: WP_010931648.1, Tohama I strain ptxA2 type: NP_882282.1, B592 strain ptxA2 type: AJ245367.1, 287 strain ptxA1 type: AJ006155.1, B6 strain ptxA4 type: AJ506336.1);

ptxC gene (Tohama I strain ptxC1 type:NP_882286.1, 3779 strain ptxC1 type: AAA22985.1, NK strain ptxC2 type: AJ420987.1);

ptxP gene (Tohama I strain ptxP1 type: FN252323.1, B2983 strain ptxP3 type: FN252324.1);

prn gene(vaccine CS strain:WP_010930159 0.1, Tohama I strain: NP_879839.1, B391 strain prn1 type: AJ011091.1, B345 strain prn2 type: AJ011092.1, B343 strain prn3 type: AJ011093.1);

fim2 gene (vaccine CS strain fim2-1 type: WP_010930199.1, Tohama I strain: NP_879898.1, NK strain fim2-2 type: AJ420988.1);

fim3 gene (vaccine CS strain: WP_010930436.1, Tohama I strain fim3-1 type: NP_880302.1, fim3-2 type: AY464179.1, fim3-3 type:AY464180.1, fim3-4 type: AY464181.1). The homologous sequences of virulence factors in the genome sequence of Shenzhen strains were extracted using Blastn, and the sequence alignment and phylogenetic analysis of the homologous genes of ptxA, ptxC, ptxP, prn, fim2 and fim3 sequences was performed using Mega 7.0 software. Genomic relationships of the major antigen was visualize by using GrapeTree software [23].

Antimicrobial susceptibility tests

We thawed fifty B. pertussis isolates on the charcoal agar which containing 10% sheep blood for 4 days. Then subcultured the bacterial suspension with 0.5 McFarland standard on the same kind of culture medium for 72 h. Kirby-Bauer disk diffusion test was used to determine susceptibility to erythromycin. An inhibition diameter > 42 mm suggested complete susceptibility to erythromycin according to some studies [24]. Staphylococcus aureus ATCC25923 and B. pertussis ATCC9797 were used as controls.

Single nucleotide polymorphism (SNP) extraction and phylogenetic typing

Eight hundred and forty-two public genome sequences were downloaded from GenBank (!/prokaryotes/1008/, as of September 2020). Snippy software ( was used to detect the core genome SNP sites of the public sequences and the Shenzhen strains sequences, and the genome sequence of Tohama I strain (GenBank Number: BX470248.1) was used as a reference. TRF software ( and BLASTn were used to identify the recombination region of the reference genome and to remove the SNP sites located in the repeating region. Based on non-repetitive regions of SNP sites, the maximum likelihood phylogenetic tree was built using IQ-Tree software.


Laboratory test results and epidemiological characteristics analysis

Fifty strains of Bordetella pertussis were obtained from nasopharyngeal swabs of 387 children with bedside inoculation, with a positivity rate of 16.79% (65/387). The epidemiological characteristics of patients are shown in Table 1. The vaccination history of eight out of 50 cases was not available. Of the remaining 42 cases, 11 (26.2%) had not received any pertussis-containing vaccine, 5 (11.9%) had completed the primary course (three doses) and 7 (16.7%) had also received a booster (Table 1). Altogether, 30 of 42 (71.4%) cases have an incomplete pertussis vaccination history according to the national recommendation, 14 of 30 (46.7%) cases have not received the recommended number of immunizations according to their age, while the other 16 cases were infected before the recommended age for vaccination.

Table 1 Epidemiological characteristics and vaccination status of children

Genotyping and erythromycin-resistant of B. pertussis isolates

Analysis of the six genes revealed that two had 100% similarity (1242/1242 bases) with those in the vaccine strains CS and Tohama l, namely fim2-1 and fim3-1 respectively. The nucleotide sequences of the remaining four genes were distinct and non-synonymous resulting in a change in antigen type. All ptxA gene sequences in the Shenzen isolates had non-synonymous mutations. All the clinical isolates had the alleleic profile ptxA1 and were distinct from the vaccine strain allele, ptxA2. Non-synonymous mutations also occurred in the pertussis toxin promoter gene (ptxP). Nineteen strains had the same allelic sequence as the vaccine strain (ptxP1), and the remaining 31 were ptxP3. Synonymous mutations were found in the pertussis toxin S3 subunit gene (ptxC). Thirty-one isolates had the allelic profile ptxC2 (Table 2). Three allelic profiles of the pertactin gene (prn) were found in the 50 isolates including, prn1 (n = 18), prn2 (n = 31) and prn3 (n = 1). The major genotypes of the 50 isolates can be divided into three categories, prn1/ptxC1/ptxP1/ptxA1/ fim3-1/fim2-1 (18/50, 36%), prn2/ptxC2/ptxP3/ptxA1/fim3-1/fim2-1(31/50, 62%) and prn3/ptxC1/ ptxP1/ptxA1/fim3-1/fim2-1 (1/50, 2%) (Table 2). Each circle of the minimum spanning tree represents an genotype of the major antigen and discern the genetic relationship among three categories (Fig. 1). All the 18 strains with A2047G mutation site (PG3) were erythromycin resistant while the other 32 strains were sensitive to erythromycin (zone of inhibition > 42 mm).

Table 2 Main virulence genotypes of clinical isolates in Shenzhen
Fig. 1

The minimum spanning tree of 50 Shenzhen stains. The size of the circle indicates the number of strains with this particular genotype. Each color represents one sample. The length of the connecting lines between the circles indicates the genetic distance between different genotypes

Genome evolution analysis

Genetic analysis was performed on the whole genome sequences from the 50 study isolates, the CS and Tohama I strain and 842 publicly available B. pertussis strains (Fig. 2). The Shenzhen isolates in this study were distributed in three evolutionary branches (phylogenetic group, PG1 to 3). The genotype prn2/ptxA1/ptxP3/ptxC2/ fim3-1 /fim2-1 predominated in PG1, genotype prn3/ptxA1/ptxP1/ptxC1/fim3-1/fim2-1 in PG2, and genotype prn1/ptxA1/ptxP1/ptxC1/fim3-1/fim2-1 in PG3. The isolates in PG3 were closely related to the epidemic strains in northern China in recent years, but distant from foreign strains. All isolates were genetically distant from vaccine strains CS and Tohama I (Fig. 2).

Fig. 2

The phylogenetic tree of Shenzhen strains and international strains. A Maximum likelihood phylogenetic tree of 50 Shenzhen strains and 842 international strains. B Maximum likelihood phylogeny tree of Shenzhen strain, phylogenetic branch PG1-PG3 strain and reference strain. Red and blue represented the Shenzhen strain and two reference strains, respectively. For better visualization, the longer phylogenetic branches were artificially shortened (double slashes)


Bordetella pertussis is the only species in the genus Bordetella that can produce PT, which is its major virulence factor and consists of five subunits (S1 to S5). Pertussis toxin belongs to a class of toxins known as A-B toxins. The A domain contains the S1 subunit which has immunoprotective properties and a variety of enzyme activities, and is expressed by ptxA gene [21]. The region between amino acid 65 and 233 is the conserved sequence of pertussis toxin, causing the host cell immune response and is the antigen recognition site of T cells. According to the nucleotide differences of the 68th, 228th, and 232nd amino acids of S1, the ptxA gene is divided into four allelic types, 1, 2, 3 and 4 [22]. The ptxA gene subtypes of pertussis vaccine strains and strains isolated before or early in vaccination were allelic types ptxA2 or ptxA4 in many countries. The reference strains Tohama I and CS in this study belong to the ptxA2 allele, while all of the 50 isolates from Shenzhen had non-synonymous mutations (A to G), that is, from the original isoleucine (I), coded for by ATA to methionine (M), coded for by ATG and belong to the ptxA1 allele.The protective effect of the vaccine produced by the vaccine containing the ptxA2 vaccine strain is weak against the strain containing the ptxA1 gene [25]. The pertussis toxin promoter (ptxP) region is close to 170 bases located upstream of the pertussis toxin gene in the B. pertussis genome. The lack of expression of pertussis toxin in B. parapertussis and B. bronchiseptica is due to differences in the promoter regions of the operons, which have species-specific differences [26]. This region contains binding sites for RNA polymerase and 6 binding sites of BvgA regulatory protein dimers, which can promote the transcriptional expression of PT by interacting with the regulatory protein BvgA. In recent years, researchers have discovered that Bordetella pertussis ptxP presents a certain polymorphism, and according to the mutation bases at specific positions, Bordetella pertussis ptxP is divided into 11 subtypes. Dutch researchers have confirmed that strains containing ptxP3 can highly express PT and have gradually become the dominant epidemic strains in the Netherlands. It is speculated that the “pertussis recurrence” may be related to this. Therefore, the ptxP gene of Bordetella pertussis has attracted more and more attention and has been used as a molecular genetic marker for the evolution of Bordetella pertussis [10]. The ptxP3 allele accounted for 62% of the 50 pertussis strains isolated in this study, which is consistent with the international trend.

The pertactin protein is a non-ciliary outer membrane protein produced by Bordetella pertussis, as an important virulence factor, PRN protein plays an important role in bacterial infection and adhesion to the epithelial cell membrane of the host respiratory system. PRN protein is also an important protective antigen that can induce humoral and cellular immune responses in mice. Pertactin has shown a good protection rate in the respiratory tract attack animal model test of Bordetella pertussis [21]. According to the difference in the structure of this region and amino acids at specific positions, the prn gene is divided into eighteen allelic types prn1-18 [22, 27, 28]. Studies in some countries have confirmed that prn genotypes of strains isolated before or in the early stages of inoculation and pertussis vaccine strain are mainly prn1 genotypes. Since the 1980s, prn2 subtypes have appeared in some European countries. By 2000s, prn2 or prn3 had been the dominant genotypes of the isolates in the Netherlands, Finland, France, the United States, and other countries [29, 30]. Bioinformatics analyses of the amino acid secondary structure deduced from the prn gene sequence of strains isolated in different ages show that the prn genotype structure of the strains isolated after the 2000s changed to a certain extent, and the hydrophilic region appeared, which led to the change of its immunogenicity. However, the previously isolated prn1, prn7, prn10 and prn11 are completely hydrophobic region in this region 1 [31]. In this study, 62% of the isolates belonged to prn2 and 2% belonged to prn3, which was consistent with the evolution trend of prn genes in the world.

Clinical trials in some countries have reported that there is a good correlation between the anti-Fim protein antibody titer in the pertussis vaccine and the immune response in mice and the immune vaccine protection in the children. WHO guidelines for pertussis vaccine recommend that strains expressing of type 2 and 3 Fim protein antigens should be selected as the production strain in the pertussis vaccine production. As reported, these two fimbriae proteins also have certain immunogenicity and have been used as components of APV by some vaccine manufacturers. In this study, the isolates contained the same fim2-1 and fim3-1 genes as the vaccine strain, indicating that these two antigen-related genes are relatively conservative.

In this study, only 16.7% of the children had completed the full vaccination of pertussis vaccine, and about 71.4% of positive patients had not completed the basic process of APV vaccination at the time of onset, mainly because these children were relatively young at the age of onset and had not yet reached the planned immunization. The current pertussis immunization targets in China do not include pregnant women, resulting in infants receiving little pertussis antibody from the mother, and the first dose of vaccination for the child is 3 months old, leading to a longer unprotected window a long-unprotected window for infants of younger months. These may be important reasons for the high incidence of pertussis in China. Another possible reason is that the effectiveness of the pertussis vaccine declines rapidly over time [32, 33]. Studies have shown that the pertussis vaccine is highly effective within three years of vaccination, but then the immunity gradually weakens, with little protection after 7 years [34, 35].

Earlier studies have shown that the genotypes of the pertussis strains were different from those of the vaccine strains, and the hypothesis of pertussis antigenic shifts was put forward. It was hypothesized that changes in the bacterial genes were driven by immune selection pressure leading to the increase of these main epidemic strains. At the same time, the researchers also speculated that due to the antigenic drift of B. pertussis, the antigenic protein subtypes between the vaccine strains and clinical isolates were different, which weaken the immunity effect of pertussis vaccine, and may also be one of the reasons for the reemergence of pertussis in highly vaccinated countries [25]. The genetic diversity of B. pertussis varied greatly prior to mass vaccination with WCVs, and significantly decreased from the 1990s, followed by clonal expansion coinciding with the epidemic periods [36]. Therefore, molecular monitoring of epidemic strains is particularly important. This study recommends that whole genome sequencing should be used to add information on pertussis epidemic strains in different countries or regions to the established database, and to gradually establish a pertussis epidemic surveillance system that can be shared by global public health departments and researchers.


In conclusion, this study presents some explanations for pertussis resurgence, including improved laboratory detection methods and adaptation of Bordetella pertussis. The dominant antigen genotype is prn2/ptxC2/ptxP3/ptxA1/ fim3-1/fim2-1, which is closely related to the epidemic strains in the United States, Australia and many countries in Europe. Despite high rates of immunization with APV, epidemics of pertussis have recently occurred in these countries. Therefore, genomic analysis of circulating isolates of B. pertussis should be continued, for it will benefit the control of whooping cough and development of improved vaccines and therapeutic strategies.

Availability of data and materials

The data that support the findings of this study have been deposited into CNGB Sequence Archive (CNSA) [34] of China National GeneBank DataBase (CNGBdb) [35] with Accession Number CNP0001528. The datasets generated for this study are available on request to the corresponding author.


  1. 1.

    Esposito S, Stefanelli P, Fry NK, Fedele G, He Q, Paterson P, et al. Pertussis prevention: reasons for resurgence, and differences in the current acellular pertussis vaccines. Front Immunol. 2019;10:1344.

    CAS  Article  Google Scholar 

  2. 2.

    Zimmerman LI, Papin JF, Warfel J, Wolf RF, Kosanke SD, Merkel TJ. Histopathology of Bordetella pertussis in the Baboon model. Infect Immun. 2018.

    Article  PubMed  PubMed Central  Google Scholar 

  3. 3.

    Schellekens J, von König CH, Gardner P. Pertussis sources of infection and routes of transmission in the vaccination era. Pediatr Infect Dis J. 2005;24(5 Suppl):S19-24.

    Article  Google Scholar 

  4. 4.

    Kuchar E, Karlikowska-Skwarnik M, Han S, Nitsch-Osuch A. Pertussis: history of the disease and current prevention failure. Adv Exp Med Biol. 2016;934:77–82.

    CAS  Article  Google Scholar 

  5. 5.

    Berkovic SF, Harkin L, McMahon JM, Pelekanos JT, Zuberi SM, Wirrell EC, et al. De-novo mutations of the sodium channel gene SCN1A in alleged vaccine encephalopathy: a retrospective study. Lancet Neurol. 2006;5(6):488–92.

    CAS  Article  Google Scholar 

  6. 6.

    Reyes IS, Hsieh DT, Laux LC, Wilfong AA. Alleged cases of vaccine encephalopathy rediagnosed years later as Dravet syndrome. Pediatrics. 2011;128(3):e699-702.

    PubMed  Google Scholar 

  7. 7.

    de Cellès MD, Magpantay FM, King AA, Rohani P. The pertussis enigma: reconciling epidemiology, immunology and evolution. Proc Biol Sci. 2016.

    Article  Google Scholar 

  8. 8.

    Jackson DW, Rohani P. Perplexities of pertussis: recent global epidemiological trends and their potential causes. Epidemiol Infect. 2014;142(4):672–84.

    CAS  Article  Google Scholar 

  9. 9.

    Robbins JB, Schneerson R, Trollfors B. Pertussis in developed countries. The Lancet. 2002;360(9334):657–8.

    Article  Google Scholar 

  10. 10.

    Mooi FR, van Loo IH, van Gent M, He Q, Bart MJ, Heuvelman KJ, et al. Bordetella pertussis strains with increased toxin production associated with pertussis resurgence. Emerg Infect Dis. 2009;15(8):1206–13.

    CAS  Article  Google Scholar 

  11. 11.

    Pertussis epidemic–Washington. MMWR morbidity and mortality weekly report. 2012; 61(28):517–22.

  12. 12.

    Xu Y, Tan Y, Asokanathan C, Zhang S, Xing D, Wang J. Characterization of co-purified acellular pertussis vaccines. Hum Vaccin Immunother. 2015;11(2):421–7.

    Article  Google Scholar 

  13. 13.

    Hu Y, Duan K, Li X, Pan C, Tian C, Tang R, et al. Development of a novel process for purification of acellular component pertussis vaccine. Chin J New Drugs. 2018;27(21):2498–504.

    Google Scholar 

  14. 14.

    Xing D, Das RG, O’Neill T, Corbel M, Dellepiane N, Milstien J. Laboratory testing of whole cell pertussis vaccine: a WHO proficiency study using the Kendrick test. Vaccine. 2001;20(3–4):342–51.

    CAS  Article  PubMed  Google Scholar 

  15. 15.

    Zhang R, Li Z, Li G, Tie Y, Li X, Gao Y, et al. A highly sensitive one-tube nested quantitative real-time PCR assay for specific detection of Bordetella pertussis using the LNA technique. Int J Infect Dis. 2020;93:224–30.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Tao Y, Tang M, Luo L, Xiang L, Xia Y, Li B, et al. Identification of etiologic agents and clinical characteristics for patients suspected of having pertussis in a large Children’s Hospital in China. Ann Transl Med. 2019;7(18):443.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  17. 17.

    Yu W, Lee LA, Liu Y, Scherpbier RW, Wen N, Zhang G, et al. Vaccine-preventable disease control in the People’s Republic of China: 1949–2016. Vaccine. 2018;36(52):8131–7.

    Article  Google Scholar 

  18. 18.

    Zhang Q, Zheng H, Liu M, Ke H, Shu J, Wu C, et al. The seroepidemiology of immunoglobulin G antibodies against pertussis toxin in China: a cross sectional study. BMC Infect Dis. 2012;12(1):138–44.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Bart MJ, Harris SR, Advani A, Arakawa Y, Bottero D, Bouchez V, et al. Global population structure and evolution of Bordetella pertussis and their relationship with vaccination. mBio. 2014;5(2): e01074.

    Article  Google Scholar 

  20. 20.

    Martin SW, Pawloski L, Williams M, Weening K, DeBolt C, Qin X, et al. Pertactin-negative Bordetella pertussis strains: evidence for a possible selective advantage. Clin Infect Dis. 2015;60(2):223–7.

    CAS  Article  Google Scholar 

  21. 21.

    Gu W, Wang K, Zhang X, Hao C, Lu Y, Wu M, et al. Pathogen analysis of pertussis-like syndrome in children. BMC Infect Dis. 2020;20(1):353–60.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Zhang S, Xu Y, Zhou Z, Wang S, Yang R, Wang J, et al. Complete genome sequence of Bordetella pertussis CS, a Chinese pertussis vaccine strain. J Bacteriol. 2011;193(15):4017–8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Zhou Z, Alikhan NF, Sergeant MJ, Luhmann N, Vaz C, Francisco AP, et al. GrapeTree: visualization of core genomic relationships among 100,000 bacterial pathogens. Genome Res. 2018;28(9):1395–404.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  24. 24.

    Bartkus JM, Juni BA, Ehresmann K, Miller CA, Sanden GN, Cassiday PK, et al. Identification of a mutation associated with erythromycin resistance in Bordetella pertussis: implications for surveillance of antimicrobial resistance. J Clin Microbiol. 2003;41(3):1167–72.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  25. 25.

    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(2):670–5.

    CAS  Article  Google Scholar 

  26. 26.

    Parkhill J, Sebaihia M, Preston A, Murphy LD, Thomson N, Harris DE, et al. Comparative analysis of the genome sequences of Bordetella pertussis, Bordetella parapertussis and Bordetella bronchiseptica. Nat Genet. 2003;35(1):32–40.

    Article  Google Scholar 

  27. 27.

    Gzyl A, Augustynowicz E, Gniadek G, Rabczenko D, Dulny G, Slusarczyk J. Sequence variation in pertussis S1 subunit toxin and pertussis genes in Bordetella pertussis strains used for the whole-cell pertussis vaccine produced in Poland since 1960: efficiency of the DTwP vaccine-induced immunity against currently circulating B. pertussis isolates. Vaccine. 2004;22(17–18):2122–8.

    CAS  Article  Google Scholar 

  28. 28.

    Simmonds K, Fathima S, Chui L, Lovgren M, Shook P, Shuel M, et al. Dominance of two genotypes of Bordetella pertussis during a period of increased pertussis activity in Alberta, Canada: January to August 2012. Int J Infect Dis. 2014;29:223–5.

    Article  PubMed  Google Scholar 

  29. 29.

    Preston NW, Carter EJ. Serotype specificity of vaccine-induced immunity to pertussis. Commun Dis Rep CDR Rev. 1992;2(13):R155–6.

    CAS  PubMed  Google Scholar 

  30. 30.

    Olin P, Rasmussen F, Gustafsson L, Hallander HO, Heijbel H, Ad Hoc Group for the Study of Pertussis Vaccines. Randomised controlled trial of two-component, three-component, and five-component acellular pertussis vaccines compared with whole-cell pertussis vaccine. Lancet. 1997;350(9091):1569–77.

    CAS  Article  Google Scholar 

  31. 31.

    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(Pt 5):355–65.

    CAS  Article  Google Scholar 

  32. 32.

    Warfel JM, Edwards KM. Pertussis vaccines and the challenge of inducing durable immunity. Curr Opin Immunol. 2015;35:48–54.

    CAS  Article  Google Scholar 

  33. 33.

    McGirr A, Fisman DN. Duration of pertussis immunity after DTaP immunization: a meta-analysis. Pediatrics. 2015;135(2):331–43.

    Article  Google Scholar 

  34. 34.

    Schwartz KL, Kwong JC, Deeks SL, Campitelli MA, Jamieson FB, Marchand-Austin A, et al. Effectiveness of pertussis vaccination and duration of immunity. CMAJ. 2016;188(16):E399-406.

    Article  Google Scholar 

  35. 35.

    Pool V, Tomovici A, Johnson DR, Greenberg DP, Decker MD. Humoral immunity 10 years after booster immunization with an adolescent and adult formulation combined tetanus, diphtheria, and 5-component acellular pertussis vaccine in the USA. Vaccine. 2018;36(17):2282–7.

    Article  Google Scholar 

  36. 36.

    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(3):680–8.

    Article  PubMed  PubMed Central  Google Scholar 

  37. 37.

    Smith AM, Guzmán CA, Walker MJ. The virulence factors of Bordetella pertussis: a matter of control. FEMS Microbiol Rev. 2001;25(3):309–33.

    CAS  Article  Google Scholar 

  38. 38.

    Kallonen T, He Q. Bordetella pertussis strain variation and evolution postvaccination. Expert Rev Vaccines. 2009;8(7):863–75.

    CAS  Article  Google Scholar 

Download references


The authors would like to express their deepest gratitude to the Laboratory of Medicine of Shenzhen Children’s Hospital and the National Institute for Communicable Disease Control and Prevention (ICDC) of China CDC.


This research was supported by the National Major Science and Technology Projects of China (No. 2017ZX10303406, No. 2018ZX10714002-003-010), Sanming Project of Medicine in Shenzhen (No. SZSM201811071), National Natural Science Foundation of China (No. 81773436, No. 81673174).

Author information




Conceived and designed the experiments: SW, XS, QH; performed the experiments: SW, HC; contributed analysis: SW, QH, CY, HZ, YZ, MJ; wrote the paper: SW, QH, XS, YH. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Xiaolu Shi.

Ethics declarations

Ethics approval and consent to participate

Bordetella pertussis isolates from the Shenzhen Center for Disease Control and Prevention were de-identified and anonymized to protect patient privacy and confidentiality; therefore, ethical clearance was not required.

Consent for publication

Not applicable.

Competing interests

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Wu, S., Hu, Q., Yang, C. et al. Molecular epidemiology of Bordetella pertussis and analysis of vaccine antigen genes from clinical isolates from Shenzhen, China. Ann Clin Microbiol Antimicrob 20, 53 (2021).

Download citation


  • Bordetella pertussis
  • Genotype
  • Pertussis vaccine
  • Molecular epidemiology
  • Erythromycin-resistant