Genomic epidemiology of Iranian Bordetella pertussis: 50 years after the implementation of whole cell vaccine

The effect of in vitro consecutive passages and culture medium on the genetic variations in BCG Pasteur 1173P2 vaccine

Since the introduction of the Bacillus Calmette-Guérin (BCG) vaccine, the genomes of vaccine strains have undergone variations due to repeated passages in different laboratories and vaccine production facilities. Genetic variations have been considered as one of the effective factors in the BCG variable protective efficacy. Consecutive subcultures have been shown to play an essential role in causing genetic variations in several microorganisms, including Mycobacterium bovis BCG. Therefore, the world health organization (WHO) recommendation to limit the passages of master seed lot in the BCG vaccine production should be considered. Besides, the role of other external variables such as quality of the raw ingredients of the culture media, the type of the culture medium and the cultivation methods in the vaccine production has been poorly studied. Here, the effect of passages and culture medium on genetic variations in a BCG seed lot was investigated during a year. The findings of this study revealed a total of 19 variants compared to seed lot while the passages were more than the number recommended by WHO. The first culture of seed lot in the Sauton broth and Middlebrook 7H9 media, and the last subculture in Sauton broth had the least and the most variants, respectively. The observation of the higher number of variants in the last cultures on Sauton broth and Middlebrook 7H9 in comparison to the first and the middle cultures may indicate the effect of passages on the genetic variations in BCG. Additionally, more variants in BCG grown in the Sauton broth do not necessarily represent the greater ability of this medium to cause genetic mutations. For a better conclusion, it is required to examine the medium components as independent variables.

Genomic dissection of the microevolution of Australian epidemic Bordetella pertussis

Whooping cough (pertussis) is a highly contagious respiratory disease caused by the bacterium Bordetella pertussis. Despite high vaccine coverage, pertussis has re-emerged in many countries including Australia and caused two large epidemics in Australia since 2007. Here, we undertook a genomic and phylogeographic study of 385 Australian B. pertussis isolates collected from 2008 to 2017. The Australian B. pertussis population was found to be composed of mostly ptxP3 strains carrying different fim3 alleles, with ptxP3-fim3A genotype expanding far more than ptxP3-fim3B. Within the former, there were six co-circulating epidemic lineages (EL1 to EL6). The multiple ELs emerged, expanded, and then declined at different time points over the two epidemics. In population genetics terms, both hard and soft selective sweeps through vaccine selection pressures have determined the population dynamics of Australian B. pertussis. Relative risk estimation suggests that once a new B. pertussis lineage emerged, it was more likely to spread locally within the first 1.5 years. However, after 1.5 years, any new lineage was likely to expand to a wider region. Phylogenetic analysis revealed the expansion of ptxP3 strains was also associated with replacement of the type III secretion system allele bscI1 with bscI3. bscI3 is associated with decreased T3SS secretion and may allow B. pertussis to reduce immune recognition. This study advanced our understanding of the epidemic population structure and spatial and temporal dynamics of B. pertussis in a highly immunized population.

Genome Characteristic of Bordetella parapertussis Isolated from Iran

Pertussis also known as whooping cough is a respiratory infection in humans particularly with severe symptoms in infants and usually caused by Bordetella pertussis. However, Bordetella parapertussis can also cause a similar clinical syndrome. During 2012 to 2015, from nasal swabs sent from different provinces to the pertussis reference laboratory of Pasture Institute of Iran for pertussis confirmation, seven B. parapertussis isolates were identified by bacterial culture, biochemical tests, and the presence of IS1001 insertion in the genome. The expression of pertactin (Prn) as one the major virulence factor for bacterial adhesion was investigated using western blot. Moreover, the genomic characteristic of one recently collected isolate, IRBP134, from a seven-month infant was investigated using Illumina NextSeq sequencing protocol. The results revealed the genome with G+C content 65% and genome size 4.7 Mbp. A total of 81 single nucleotide polymorphisms and 13 short insertions and deletions were found in the genome compared to the B. parapertussis 12822 as a reference genome showing ongoing evolutionary changes. A phylogeny relationship of IRBP134 was also investigated using global B. parapertussis available genomes.

Genomic characteristics of two most widely used BCG vaccine strains: Danish 1331 and Pasteur 1173P2

Background

Bacillus Calmette–Guérin (BCG) refers to a group of vaccine strains with unique genetic characteristics. BCG is the only available vaccine for preventing tuberculosis (TB). Genetic and biochemical variations among the BCG vaccine strains have been considered as one of the significant parameters affecting the variable protective efficacy of the vaccine against pulmonary tuberculosis. To track genetic variations, here two vaccine strains (Danish 1331 and Pasteur 1173P2) popularly used according to the BCG World Atlas were subjected to a comparative analysis against the Mycobacterium tuberculosis H37Rv, Mycobacterium bovis AF2122/97, and Mycobacterium tuberculosis variant bovis BCG str. Pasteur 1173P2 reference genomes. Besides, the presence or absence of the experimentally verified human T cell epitopes was examined.

Results

Only two variants were identified in BCG Danish 1331 that have not been reported previously in any BCG strains with the complete submitted genome yet. Furthermore, we identified a DU1-like 14,577 bp region in BCG Danish 1331; The duplication which was previously seemed to be exclusive to the BCG Pasteur. We also found that 35% of the T cell epitopes are absent from both strains, and epitope sequences are more conserved than the rest of the genome.

Conclusions

We provided a comprehensive catalog of single nucleotide polymorphisms (SNPs) and short insertions and deletions (indels) in BCG Danish 1331 and BCG Pasteur 1173P2. These findings may help determine the effect of genetic variations on the variable protective efficacy of BCG vaccine strains.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08826-9.

Characterization of Bordetella pertussis Strains Isolated from India

Despite high level vaccination and the availability of two different types of vaccines, whole cell (wP) and acellular vaccines (aP), the resurgence of pertussis has been reported in many countries. Antigenic variation within circulating and vaccine strains is the most documented reason reported for the resurgence of pertussis. Research on genetic divergence among circulating and vaccine strains has largely been reported in countries using aP vaccines. There are inadequate data available for antigenic variation in B. pertussis from wP-using countries. India has used wP for more than 40 years in their primary immunization program. The present study reports five clinical isolates of B. pertussis from samples of pediatric patients with pertussis symptoms observed in India. Genotypic and phenotypic characterization of clinical isolates were performed by serotyping, genotyping, whole genome analyses and comparative genomics. All clinical isolates showed serotype 1, 2 and 3 based on the presence of fimbriae 2 and 3. Genotyping showed genetic similarities in allele types for five aP genes within vaccine strains and clinical isolates reported from India. The presence of the ptxP3 genotype was observed in two out of five clinical isolates. Whole-genome sequencing was performed for clinical isolates using the hybrid strategy of combining Illumina (short reads) and oxford nanopore (long reads) sequencing strategies. Clinical isolates (n = 5) and vaccine strains (n = 7) genomes of B. pertussis from India were compared with 744 B. pertussis closed genomes available in the public databases. The phylogenomic comparison of B. pertussis genomes reported from India will be advantageous in better understanding pertussis resurgence reported globally with respect to pathogen adaptation.

Evolution of Bordetella pertussis in the acellular vaccine era in Norway, 1996 to 2019

We described the population structure of Bordetella pertussis (B. pertussis) in Norway from 1996 to 2019 and determined if there were evolutionary shifts and whether these correlated with changes in the childhood immunization program. We selected 180 B. pertussis isolates, 22 from the whole cell vaccine (WCV) era (1996–1997) and 158 from the acellular vaccine (ACV) era (1998–2019). We conducted whole genome sequencing and determined the distribution and frequency of allelic variants and temporal changes of ACV genes. Norwegian B. pertussis isolates were evenly distributed across a phylogenetic tree that included global strains. We identified seven different allelic profiles of ACV genes (A–F), in which profiles A1, A2, and B dominated (89%), all having pertussis toxin (ptxA) allele 1, pertussis toxin promoter (ptxP) allele 3, and pertactin (prn) allele 2 present. Isolates with ptxP1 and prn1 were not detected after 2007, whereas the prn2 allele likely emerged prior to 1972, and ptxP3 before the early 1980s. Allele conversions of ACV genes all occurred prior to the introduction of ACV. Sixteen percent of our isolates showed mutations within the prn gene. ACV and its booster doses (implemented for children in 2007 and adolescents in 2013) might have contributed to evolvement of a more uniform B. pertussis population, with recent circulating strains having ptxA1, ptxP3, and prn2 present, and an increasing number of prn mutations. These strains clearly deviate from ACV strains (ptxA1, ptxP1, prn1), and this could have implications for vaccine efficiency and, therefore, prevention and control of pertussis.

Global spatial dynamics and vaccine-induced fitness changes of Bordetella pertussis

As with other pathogens, competitive interactions between Bordetella pertussis strains drive infection risk. Vaccines are thought to perturb strain diversity through shifts in immune pressures; however, this has rarely been measured because of inadequate data and analytical tools. We used 3344 sequences from 23 countries to show that, on average, there are 28.1 transmission chains circulating within a subnational region, with the number of chains strongly associated with host population size. It took 5 to 10 years for B. pertussis to be homogeneously distributed throughout Europe, with the same time frame required for the United States. Increased fitness of pertactin-deficient strains after implementation of acellular vaccines, but reduced fitness otherwise, can explain long-term genotype dynamics. These findings highlight the role of vaccine policy in shifting local diversity of a pathogen that is responsible for 160,000 deaths annually.

Pertactin-Deficient Bordetella pertussis with Unusual Mechanism of Pertactin Disruption, Spain, 1986-2018

Bordetella pertussis not expressing pertactin has increased in countries using acellular pertussis vaccines (ACV). The deficiency is mostly caused by pertactin gene disruption by IS481. To assess the effect of the transition from whole-cell vaccine to ACV on the emergence of B. pertussis not expressing pertactin in Spain, we studied 342 isolates collected during 1986–2018. We identified 93 pertactin-deficient isolates. All were detected after introduction of ACV and represented 38% of isolates collected during the ACV period; 58.1% belonged to a genetic cluster of isolates carrying the unusual prn::del(–292, 1340) mutation. Pertactin inactivation by IS481 insertion was identified in 23.7% of pertactin-deficient isolates, arising independently multiple times and in different phylogenetic branches. Our findings support the emergence and dissemination of a cluster of B. pertussis with an infrequent mechanism of pertactin disruption in Spain, probably resulting from introduction of ACV.

Genomic evidence for stability of the Bacillus Calmette-Guérin (BCG) vaccine strain (Pasteur 1173P2) from different batches in Iran

AIMS

Investigate the genetic stability of the BCG vaccine produced in Iran from different batches compared to the reference strain.

METHODS AND RESULTS

We comparatively analyzed the whole genome sequences of the vaccine batches from different years. Eleven vials of different batches from 2010, 2018, and 2019 were included. Complete genome analyses revealed no difference between the old (2010) and new (2018 and 2019) vaccine batches. Additionally, minor genetic changes include five single nucleotide polymorphisms (SNPs) and insertions/deletions (indels) were observed compared to the BCG Pasteur 1173P2 reference strain, which were shared among all batches. Besides, the batches were identical to the reference strain in terms of antibiotic resistance genes, prophage sequences, and clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated (Cas) systems.

CONCLUSIONS

High genetic stability of the BCG vaccine used in the national immunization program was confirmed, which indicates the optimal conditions in the vaccine production process.

SIGNIFICANCE AND IMPACT OF THE STUDY

Genetic differences within and between vaccine strains have been declared as one of the main parameters related to the BCG vaccine variable protective efficacy. No study has been done to investigate the genetic variations of the vaccine batches at the single-base level.

Evaluation of Outer Membrane Vesicles Obtained from Predominant Local Isolate of Boredetella pertussis as a Vaccine Candidate

Background:

Pertussis is a current contagious bacterial disease caused by Bp. Given the prevalence of pertussis, development of new vaccines is important. This study was attempted to evaluate the expression of main virulence factors (PTX, PRN, and FHA) from Bp predominant strains and also compare the expression of these factors in the OMVs obtained from predominant circulating Bp isolate.

Methods:

The physicochemical features of the prepared OMVs were analyzed by electron microscopy and SDS-PAGE. The presence of the mentioned virulence factors was confirmed by Western blotting. BALB/c mice (n = 21) immunized with characterized OMVs were challenged intranasally with sublethal doses of Bp, to examine their protective capacity.

Results:

Electron microscopic examination of the OMVs indicated vesicles within the range of 40 to 200 nm. SDS-PAGE and Western blotting demonstrated the expression of all three main protective immunogens (PTX, PRN, and FHA), prevalent in the predominant, challenge, and vaccine strains, and OMVs of the predominant IR37 strain and BP134 vaccine strain. Significant differences were observed in lung bacterial counts between the immunized mice with OMV (30 CFU/lung) compared to the negative control group ((6 10⁴ CFU/lung; p < 0.001). In mice immunized with OMVs (3 µg), the number of lungs recovered colonies after five days dropped at least five orders of magnitude compared to the control group.

Conclusion:

OMVs obtained from circulating isolates with the predominant profile may constitute a highly promising vaccine quality. They also can be proposed as a potential basic material for the development of new pertussis vaccine candidate.

Evolution of Bordetella pertussis over a 23-year period in France, 1996 to 2018

BackgroundBordetella pertussis is the main agent of whooping cough. Vaccination with acellular pertussis vaccines has been largely implemented in high-income countries. These vaccines contain 1 to 5 antigens: pertussis toxin (PT), filamentous haemagglutinin (FHA), pertactin (PRN) and/or fimbrial proteins (FIM2 and FIM3). Monitoring the emergence of B. pertussis isolates that might partially escape vaccine-induced immunity is an essential component of public health strategies to control whooping cough.AimWe aimed to investigate temporal trends of fimbriae serotypes and vaccine antigen-expression in B. pertussis over a 23-year period in France (1996-2018).MethodsIsolates (n = 2,280) were collected through hospital surveillance, capturing one third of hospitalised paediatric pertussis cases. We assayed PT, FHA and PRN production by Western blot (n = 1,428) and fimbriae production by serotyping (n = 1,058). Molecular events underlying antigen deficiency were investigated by genomic sequencing.ResultsThe proportion of PRN-deficient B. pertussis isolates has increased steadily from 0% (0/38) in 2003 to 48.4% (31/64) in 2018 (chi-squared test for trend, p < 0.0001), whereas only 5 PT-, 5 FHA- and 9 FIM-deficient isolates were found. Impairment of PRN production was predominantly due to IS481 insertion within the prn gene or a 22 kb genomic inversion involving the prn promoter sequence, indicative of convergent evolution. FIM2-expressing isolates have emerged since 2011 at the expense of FIM3.ConclusionsB. pertussis is evolving through the rapid increase of PRN-deficient isolates and a recent shift from FIM3 to FIM2 expression. Excluding PRN, the loss of vaccine antigen expression by circulating B. pertussis isolates is epidemiologically insignificant.

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.

Isolate-Based Surveillance of Bordetella pertussis, Austria, 2018-2020

Pertussis is a vaccine-preventable disease, and its recent resurgence might be attributable to the emergence of strains that differ genetically from the vaccine strain. We describe a novel pertussis isolate-based surveillance system and a core genome multilocus sequence typing scheme to assess Bordetella pertussis genetic variability and investigate the increased incidence of pertussis in Austria. During 2018–2020, we obtained 123 B. pertussis isolates and typed them with the new scheme (2,983 targets and preliminary cluster threshold of <6 alleles). B. pertussis isolates in Austria differed genetically from the vaccine strain, both in their core genomes and in their vaccine antigen genes; 31.7% of the isolates were pertactin-deficient. We detected 8 clusters, 1 of them with pertactin-deficient isolates and possibly part of a local outbreak. National expansion of the isolate-based surveillance system is needed to implement pertussis-control strategies.

Evolutionary genomics of recent clinical Bordetella pertussis isolates from Iran: wide circulation of multiple ptxP3 lineages and report of the first ptxP3 filamentous hemagglutinin-negative B. pertussis

Here we investigated nationwide clinical Bordetella pertussis isolated during 2005-2017 from different provinces of Iran, a country with more than 50 years whole cell vaccine immunisation history. Our results revealed the ongoing increase in the population of ptxP3/fim3-2 B. pertussis isolates in different provinces which were differentiated into nine clades. The largest clade (clade 8) which was previously found to be prevalent in Tehran was also prevalent across the country and clade 5 with ptxP3/prn9 genotype has also increased in frequency (14% of all ptxP3 isolates) in recent years. Furthermore, we detected the first ptxP3 B. pertussis isolates that does not express filamentous hemagglutinin (FhaB) as one of the major antigens of the pathogen and a key component of the acellular pertussis vaccine.

Hospital-based prospective study of pertussis in infants and close contacts in Tehran, Iran

Background

Pertussis remain a global health concern, especially in infants too young to initiate their vaccination. Effective vaccination and high coverage limit the circulation of the pathogen, yet duration of protection is limited and boosters are recommended during a lifetime. In Iran, boosters are given at 18 months and 6 years old using whole pertussis vaccines for which efficacy is not known, and pertussis surveillance is scant with only sporadic biological diagnosis. Burden of pertussis is not well understood and local data are needed.

Methods

Hospital-based prospective study implementing molecular laboratory testing in infants aged ≤6 months and presenting ≥5 days of cough associated to one pertussis-like symptom in Tehran. Household and non-household contact cases of positive infants were evaluated by comprehensive pertussis diagnosis (molecular testing and serology) regardless of clinical signs. Clinical evaluation and source of infection were described.

Results

A total of 247 infants and 130 contact cases were enrolled. Pertussis diagnosis result was obtained for 199 infants and 104 contact cases. Infant population was mostly < 3 months old (79.9%; 157/199) and unvaccinated (62.3%; 124/199), 20.1% (40/199) of them were confirmed having B. pertussis infection. Greater cough duration and lymphocyte counts were the only symptoms associated to positivity. Half of the contact cases (51.0%; 53/104) had a B. pertussis infection, median age was 31 years old. A proportion of 28.3% (15/53) positive contacts did not report any symptom. However, 67.9% (36/53) and 3.8% (2/53) of them reported cough at inclusion or during the study, including 20.8% (11/53) who started coughing ≥7 days before infant cough onset. Overall, only five samples were successfully cultured.

Conclusion

These data evidenced the significant prevalence of pertussis infection among paucy or poorly symptomatic contacts of infants with pertussis infection. Widespread usage of molecular testing should be implemented to identify B. pertussis infections.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-021-06266-6.

Prevalence and characterization of pertactin deficient Bordetella pertussis strains in Brazil, a whole-cell vaccine country

Many countries have reported antigenic divergence among circulating Bordetella pertussis strains, mainly in those countries which introduced the acellular pertussis (aP) vaccine. This phenomenon can be seen, for example, with the recent rise of pertactin (Prn)-deficient B. pertussis strains, one of the antigens included in aP vaccine formulas. The whole cell pertussis (wP) vaccine has been used in Brazil since 1977 for the primary pertussis, diphtheria and tetanus immunization series. In 2014, the aP vaccine was recommended for women during pregnancy to protect infants in the first months of life. Our objective was to determine the prevalence of Prn-deficiency in 511 isolates of B. pertussis collected in Brazil during 2010-2016. All isolates were characterized, through PFGE and serotyping, and screened for the loss of Prn by ELISA. Prn-deficiency was confirmed by immunoblotting, and identification of the possible genetic markers was performed with PCR and Sanger sequencing. Results indicate that 110 PFGE profiles are currently circulating, with five profiles representing the majority, and the predominant serotype 3, has been gradually replaced by serotype 2 and serotype 2,3. ELISA screening and immunoblotting identified three Prn-deficient isolates. Genotypic characterization by PCR and sequencing indicated that one isolate had a promoter mutation in prn, while the other two did not have an obvious genetic explanation for their deficiency. While the lack of Prn was identified in a few isolates, this study did not detect a relevant occurrence of Prn-deficiency, until 2016, confirming previous observations that Prn-deficiency is likely aP vaccine-driven.

Bordetella pertussis in School-Age Children, Adolescents, and Adults: A Systematic Review of Epidemiology, Burden, and Mortality in the Middle East

Despite modern diphtheria-tetanus-pertussis (DTP) vaccines and high vaccine coverage, a resurgence of pertussis (whooping cough) has been observed globally. In North America and Europe, high vaccine coverage in children has led to a shift in the age-specific peak incidence of infection away from infants and towards older children and adolescents. However, much less is known about the prevalence of pertussis in older children and adults in the Middle East. A systematic search of MEDLINE, EMBASE, and BIOSIS was undertaken to identify studies published between 1 January 1990 and 17 June 2019, with information on pertussis epidemiology, burden of illness, and mortality in school-aged children, adolescents, and adults in the Middle East. Studies identified for inclusion were reviewed narratively because a statistical comparison was not possible because of the mix of methodologies used. The results showed that surveillance data are weak or missing in most Middle Eastern countries, and among 24 epidemiological studies identified, most were from Iran (14), Israel (4), and Turkey (3), with single studies from the United Arab Emirates and Iraq. Despite various surveillance periods, clinical definitions, and antibody cut-off values used across the studies, the reported seroprevalence of pertussis antibodies suggested that adolescents and adults are commonly exposed to pertussis in the community and that vaccine-acquired immunity from childhood wanes. Few countries in the Middle East include a diphtheria-tetanus-acellular pertussis (Tdap) booster for adolescents on the national schedule. Israel was the only country with epidemiological data in a population that received Tdap, and the study showed that after the introduction of the adolescent booster dose, there was decrease in pertussis among children aged 5–14 years. To conclude, results from the Middle East suggest that in common with other regions, pertussis is widely circulating and that it might be shifting towards older age groups.

The global prevalence ptxP3 lineage of Bordetella pertussis was rare in young children with the co-purified aPV vaccination: a 5 years retrospective study

Background

The global prevalent ptxP3 strains varies from about 10% to about 50% of circulating B. pertussis population in different areas of China.

Methods

To investigate the difference of vaccination status between different genotypes in the circulating B. pertussis after 10 years of acellular pertussis vaccine (aPV) used in China. The nasopharyngeal swabs and isolates of B. pertussis from these patients were used to perform genotyping of antigen genes. We use antibiotic susceptibility test against erythromycin and sequencing methods for site 2047 of 23S rRNA to determine the resistance status.

Results

The ptxP1 allele with erythromycin resistant (ER) B. pertussis infection (total of 449 subjects) consisted of 84.70 to 96.70% from 2012 to 2016 in this study. Vaccinated with co-purified aPV was found in 133(133/403,33.0%), 1(1/9,11.1%) and 2(2/21,9.5%) in ptxP1/fhaB3-ER, ptxP1/fhaB2-ES and ptxP3/fhaB2-ES B. pertussis infected children each, which showed a significant difference (χ² = 6.87, P = 0.032).

Conclusions

The ptxP3-ES B. pertussis was rare while the ptxP1-ER B. pertussis was steadily increased in Xi’an, China from 2012 to 2016, where co-purified aPV was prevalent used. This pose a hypothesis that the co-purified aPV might protect against ptxP3 strains more efficient, which generated a rare chance for ptxP3 strains to be under the antibiotic pressure and further developed to be erythromycin resistance. A further cohort study and the mechanisms of the additional antigen proteins of co-purified aPV protected against B. pertussis should be consideration.

Comparative genomics of whole-cell pertussis vaccine strains from India

BACKGROUND

Despite high vaccination coverage using acellular (ACV) and whole-cell pertussis (WCV) vaccines, the resurgence of pertussis is observed globally. Genetic divergence in circulating strains of Bordetella pertussis has been reported as one of the contributing factors for the resurgence of the disease. Our current knowledge of B. pertussis genetic evolution in circulating strains is mostly based on studies conducted in countries using ACVs targeting only a few antigens used in the production of ACVs. To better understand the adaptation to vaccine-induced selection pressure, it will be essential to study B. pertussis populations in developing countries which are using WCVs. India is a significant user and global supplier of WCVs. We report here comparative genome analyses of vaccine and clinical isolates reported from India. Whole-genome sequences obtained from vaccine strains: WCV (J445, J446, J447 and J448), ACV (BP165) were compared with Tohama-I reference strain and recently reported clinical isolates from India (BPD1, BPD2). Core genome-based phylogenetic analysis was also performed using 166 isolates reported from countries using ACV.

RESULTS

Whole-genome analysis of vaccine and clinical isolates reported from India revealed high genetic similarity and conserved genome among strains. Phylogenetic analysis showed that clinical and vaccine strains share genetic closeness with reference strain Tohama-I. The allelic profile of vaccine strains (J445:ptxP1/ptxA2/prn1/fim2-1/fim3-1; J446: ptxP2/ptxA4/prn7/fim2-2/fim3-1; J447 and J448: ptxP1/ptxA1/ prn1/fim2-1/fim3-1), which matched entirely with clinical isolates (BPD1:ptxP1/ptxA1/prn1/fim2-1 and BPD2: ptxP1/ptxA1/prn1/fim2-1) reported from India. Multi-locus sequence typing (MLST) demonstrated the presence of dominant sequence types ST2 and primitive ST1 in vaccine strains which will allow better coverage against circulating strains of B. pertussis.

CONCLUSIONS

The study provides a detailed characterization of vaccine and clinical strains reported from India, which will further facilitate epidemiological studies on genetic shifts in countries which are using WCVs in their immunization programs.