Dominance of two genotypes of Bordetella pertussis during a period of increased pertussis activity in Alberta, Canada: January to August 2012

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

Background

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.

Methods

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 (https://www.ncbi.nlm.nih.gov/genbank/). The antimicrobial susceptibility was test by Kirby-Bauer disk diffusion.

Results

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.

Conclusions

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.

A Call for Caution in Use of Pertussis Vaccine Effectiveness Studies to Estimate Waning Immunity: A Canadian Immunization Research Network Study

BACKGROUND

Vaccine effectiveness (VE) studies provide essential evidence on waning vaccine-derived immunity, a major threat to pertussis control. We evaluated how study design affects estimates by comparing 2 case-control studies conducted in Ontario, Canada.

METHODS

We compared results from a test-negative design (TND) with a frequency-matched design (FMD) case-control study using pertussis cases from 2005-2015. In the first study, we identified test-negative controls from the public health laboratory that diagnosed cases and, in the second, randomly selected controls from patients attending the same physicians that reported cases, frequency matched on age and year. We compared characteristics of cases and controls using standardized differences.

RESULTS

In both designs, VE estimates for the early years postimmunization were consistent with clinical trials (TND, 84%; FMD, 89% at 1-3 years postvaccination) but diverged as time since last vaccination increased (TND, 41%; FMD, 74% by 8 years postvaccination). Overall, we observed lower VE and faster waning in the TND than the FMD. In the TND but not FMD, controls differed from cases in important confounders, being younger, having more comorbidities, and higher healthcare use. Differences between the controls of each design were greater than differences between cases. TND controls were more likely to be unvaccinated or incompletely vaccinated than FMD controls (P < .001).

CONCLUSIONS

The FMD adjusted better for healthcare-seeking behavior than the TND. Duration of protection from pertussis vaccines is unclear because estimates vary by study design. Caution should be exercised by experts, researchers, and decision makers when evaluating evidence on optimal timing of boosters.

Clinical characteristics, molecular epidemiology and antimicrobial susceptibility of pertussis among children in southern China

BACKGROUND

Increasing numbers of pertussis cases have been reported in recent years. The reported cases from Shenzhen Children's Hospital were close to one tenth of all cases in China. The epidemiology of antigenic genotype and antibiotic resistance of circulating strains in children have been unknown in Shenzhen, southern China. The aim of this study was to describe the clinical features and explore the genotypes and antimicrobial susceptibility of circulating Bordetella pertussis among children in Shenzhen.

METHODS

Data of hospitalized children with pertussis in Shenzhen Children's Hospital from August 2015 to April 2017 were collected. The genetic variability of isolates was investigated and Etest was performed for phenotypic susceptibility to erythromycin, azithromycin, clarithromycin, clindamycin, and trimethoprim/sulfamethoxazole.

RESULTS

469 children with pertussis confirmed by real-time quantitative polymerase chain reaction were hospitalized and strains were isolated from 105 patients. White blood cell count ≥ 20 × 10⁹/L and lymphocyte proportion ≥ 60% were observed in 39.29% of infants younger than 3 months. The two predominant profiles of virulence-associated allelic genes were ptxA1/ptxC1/ptxP1/prn1 (48.6%) and ptxA1/ptxC2/ptxP3/prn2 (44.8%). Among the isolates, 48.6% (51/105) were found resistant to macrolides.

CONCLUSIONS

These findings indicate that leukocytosis is not a sensitive indicator of pertussis. Isolates with the gene profile ptxP3/prn2 were highly circulating in Shenzhen and less resistant to macrolides, different from patterns observed in other parts of China.

Molecular Epidemiology of Bordetella pertussis

Although vaccination has been effective, Bordetella pertussis is increasingly causing epidemics, especially in industrialized countries using acellular vaccines (aPs). One factor behind the increased circulation is the molecular changes on the pathogen level. After pertussis vaccinations were introduced, changes in the fimbrial (Fim) serotype of the circulating strains was observed. When bacterial typing methods improved, further changes between the vaccine and circulating strains, especially among the common virulence genes including pertussis toxin (PT) and pertactin (PRN) were noticed. Moreover, development of genome based techniques including pulsed-field gel electrophoresis (PFGE), multiple-locus variable number tandem repeat analysis (MLVA) and whole-genome sequencing (WGS) have offered a better resolution to monitor B. pertussis strains. After the introduction of aP vaccines, B. pertussis strains that are deficient to vaccine antigens, especially PRN, have appeared widely. On the other hand, antimicrobial resistance to first line drugs (macrolides) against B. pertussis is still low in many countries and therefore no globally evaluated antimicrobial susceptibility test values have been recommended. In this review, we focus on the molecular changes in the bacteria, which have or may have affected the past and current epidemiology of pertussis.

Short-Read Whole-Genome Sequencing for Laboratory-Based Surveillance of Bordetella pertussis

Bordetella pertussis is a Gram-negative bacterium that causes respiratory infections in humans. Ongoing molecular surveillance of B. pertussis acellular vaccine (aP) antigens is critical for understanding the interaction between evolutionary pressures, disease pathogenesis, and vaccine effectiveness. Methods currently used to characterize aP components are relatively labor-intensive and low throughput. To address this challenge, we sought to derive aP antigen genotypes from minimally processed short-read whole-genome sequencing data generated from 40 clinical B. pertussis isolates and analyzed using the SRST2 bioinformatic package. SRST2 was able to identify aP antigen genotypes for all antigens with the exception of pertactin, possibly due to low read coverage in GC-rich low-complexity regions of variation. Two main genotypes were observed in addition to a singular third genotype that contained an 84-bp deletion that was identified by SRST2 despite the issues in allele calling. This method has the potential to generate large pools of B. pertussis molecular data that can be linked to clinical and epidemiological information to facilitate research of vaccine effectiveness and disease severity in the context of emerging vaccine antigen-deficient strains.

The characterization of Bordetella pertussis strains isolated in the Central-Western region of Brazil suggests the selection of a specific genetic profile during 2012-2014 outbreaks

Pertussis is a worldwide acute respiratory disease caused by the bacterium Bordetella pertussis. Despite high vaccine coverage, the bacterium continues to circulate in populations and is still one of the most common vaccine-preventable diseases. In Brazil, pertussis incidence has presented a significant decrease since 1990 but since 2011 a sudden increase in incidence has been observed. Thus, the aim of this study was to perform a molecular epidemiological characterization of B. pertussis strains isolated in the Central-Western region (specifically in Distrito Federal) of Brazil from August 2012 to August 2014. During this period, 92 B. pertussis strains were isolated from the outbreaks. All strains were characterized by serotyping and XbaI pulsed-field gel electrophoresis profiles. From August to December 2012, the most prevalent serotype observed was 1,3 (13/17). During 2013 the prevalence of serotype 1,3 decreased (13/30) and from January 2014 to August 2014 the most prevalent serotype was 1,2 (33/45). Fourteen PFGE profiles were identified. Of these, BP-XbaI0039 prevalence increased from 3/17 in 2012 to 10/30 in 2013, and 35/45 in 2014. These results evidence the selection of a specific genetic profile during this period, suggesting the occurrence of a bacterial genomic profile with high circulation potential.

Time series analysis of temporal trends in the pertussis incidence in Mainland China from 2005 to 2016

Short-term forecast of pertussis incidence is helpful for advanced warning and planning resource needs for future epidemics. By utilizing the Auto-Regressive Integrated Moving Average (ARIMA) model and Exponential Smoothing (ETS) model as alterative models with R software, this paper analyzed data from Chinese Center for Disease Control and Prevention (China CDC) between January 2005 and June 2016. The ARIMA (0,1,0)(1,1,1)₁₂ model (AICc = 1342.2 BIC = 1350.3) was selected as the best performing ARIMA model and the ETS (M,N,M) model (AICc = 1678.6, BIC = 1715.4) was selected as the best performing ETS model, and the ETS (M,N,M) model with the minimum RMSE was finally selected for in-sample-simulation and out-of-sample forecasting. Descriptive statistics showed that the reported number of pertussis cases by China CDC increased by 66.20% from 2005 (4058 cases) to 2015 (6744 cases). According to Hodrick-Prescott filter, there was an apparent cyclicity and seasonality in the pertussis reports. In out of sample forecasting, the model forecasted a relatively high incidence cases in 2016, which predicates an increasing risk of ongoing pertussis resurgence in the near future. In this regard, the ETS model would be a useful tool in simulating and forecasting the incidence of pertussis, and helping decision makers to take efficient decisions based on the advanced warning of disease incidence.

Antigenic and genetic characterization of Bordetella pertussis recovered from Quebec, Canada, 2002–2014: detection of a genetic shift

Despite vaccination, cyclical peaks of Bordetella pertussis incidence rates are still observed in Canada and other developed countries, making pertussis one of the most prevalent vaccine preventable bacterial diseases. In the postacellular vaccine era, evolution of bacterial strains has resulted in strains with altered vaccine antigens. Previous Canadian studies have focused on isolates mainly from the provinces of Ontario and Alberta, with only small numbers of isolates from other provinces. Therefore, in this study, we examined a larger sample (n = 52) of isolates from Quebec, Canada, between 2002 and 2014. Isolates were characterized by serotype, sequence type, and prevalence of pertactin deficiency. The Quebec isolates shared characteristics similar to other Canadian isolates and to isolates circulating globally. Although pertactin-deficient isolates were not present, a significant shift in sequence type was observed in more recent years. This study highlights the importance of continually monitoring disease-causing isolates to track evolutionary trends and gain a better understanding of the molecular epidemiology of pertussis in Canada.

Variation in Bordetella pertussis Susceptibility to Erythromycin and Virulence-Related Genotype Changes in China (1970-2014)

OBJECTIVES

To investigate changes in virulence-related genotypes and in the antimicrobial susceptibility of Bordetella pertussis isolates collected from the 1970s to 2014 in the northern part of China.

METHODS

A total of 124 B. pertussis isolates from three periods, the 1970s, 2000-2008, and May 2013-Sept 2014, were typed by multilocus sequence typing (MLST) and tested for antimicrobial susceptibility and virulence-related genes. A fragment of the 23S rRNA gene from each of the 99 isolates from 2013-2014 was amplified and sequenced.

RESULTS

All isolates from 2000-2008 and 2013-2014 were identified as ST2, whereas isolates from the 1970s were ST1. PtxA2/ptxC1/ptxP1/prn1/fim2-1/fim3-1/tcfA2, which was the same as the vaccine strain, was the only type in the 1970s. During the 2000s and 2013-2014, the virulence type ptxA1/ptxC1/ptxP1/prn1/fim2-1/fim3-1/tcfA2 was dominant, with frequencies of 68.4% and 91.9%, respectively. Nine ptxP3 strains, which were more virulent, were detected after 2000. All 124 isolates were susceptible to levofloxacin, sulphamethoxazole/trimethoprim and tetracycline. The isolates from the 1970s and 2000-2008 were susceptible to all tested macrolides, whereas 91.9% of the 2013-2014 isolates were highly resistant (minimal inhibitory concentration, MIC >256 μg/ml). No ptxP3 strain was resistant to macrolides. All erythromycin-resistant strains except for one had the A2047G mutation in the 23S rRNA gene.

CONCLUSIONS

Macrolide resistance of the B. pertussis population has been a serious problem in the northern part of China. Because most of the epidemic clone of the pathogen expresses the same antigen profiles as the vaccine strain, except ptxA, improvements in immunization strategies may prevent the spread of infection and drug resistance.