Whole-genome sequencing reveals the effect of vaccination on the evolution of Bordetella pertussis

Whole-Genome Analysis of Bordetella pertussis MT27 Isolates from School-Associated Outbreaks: Single-Nucleotide Polymorphism Diversity and Threshold of the Outbreak Strains

Bordetella pertussis, the causative agent of whooping cough, can cause pertussis outbreaks in humans, especially in school-aged children. Here, we performed whole-genome sequencing of 51 B. pertussis isolates (epidemic strain MT27) collected from patients infected during 6 school-associated outbreaks lasting less than 4 months. We compared their genetic diversity with that of 28 sporadic isolates (non-outbreak MT27 isolates) based on single-nucleotide polymorphisms (SNPs). Our temporal SNP diversity analysis revealed a mean SNP accumulation rate (time-weighted average) of 0.21 SNPs/genome/year during the outbreaks. The outbreak isolates showed a mean of 0.74 SNP differences (median, 0; range, 0 to 5) between 238 isolate pairs, whereas the sporadic isolates had a mean of 16.12 SNP differences (median, 17; range 0 to 36) between 378 isolate pairs. A low SNP diversity was observed in the outbreak isolates. Receiver operating characteristic analysis demonstrated that the optimal cutoff value to distinguish between the outbreak and sporadic isolates was 3 SNPs (Youden's index of 0.90 with a true-positive rate of 0.97 and a false-positive rate of 0.07). Based on these results, we propose an epidemiological threshold of ≤3 SNPs per genome as a reliable marker of B. pertussis strain identity during pertussis outbreaks that span less than 4 months. IMPORTANCE Bordetella pertussis is a highly infectious bacterium that easily causes pertussis outbreaks in humans, especially in school-aged children. In detection and investigation of outbreaks, excluding non-outbreak isolates is important for understanding the bacterial transmission routes. Currently, whole-genome sequencing is widely used for outbreak investigations, and the genetic relatedness of outbreak isolates is assessed based on differences in the number of single-nucleotide polymorphisms (SNPs) in the genomes of different isolates. The optimal SNP threshold defining strain identity has been proposed for many bacterial pathogens, but not for B. pertussis. In this study, we performed whole-genome sequencing of 51 B. pertussis outbreak isolates and identified a genetic threshold of ≤3 SNPs per genome as a marker defining the strain identity during pertussis outbreaks. This study provides a useful marker for identifying and analyzing pertussis outbreaks and can serve as a basis for future epidemiological studies on pertussis.

Coping Strategies for Pertussis Resurgence

Pertussis (whooping cough) is a respiratory disease caused primarily by Bordetella pertussis, a Gram-negative bacteria. Pertussis is a relatively contagious infectious disease in people of all ages, mainly affecting newborns and infants under 2 months of age. Pertussis is undergoing a resurgence despite decades of high rates of vaccination. To better cope with the challenge of pertussis resurgence, we evaluated its possible causes and potential countermeasures in the narrative review. Expanded vaccination coverage, optimized vaccination strategies, and the development of a new pertussis vaccine may contribute to the control of pertussis.

A Cross-Sectional Study Revealing the Emergence of Erythromycin-Resistant Bordetella pertussis Carrying ptxP3 Alleles in China

Background

Previous limited studies have identified that Bordetella pertussis (B. pertussis) isolates circulating in China possess distinct molecular features and high rates of erythromycin-resistance (ER). Their evolution and potential impact on the prevention and control of global pertussis are worthy of attention.

Methods

The present cross-sectional study involved 311 non-duplicate and unrelated B. pertussis strains isolated from Chinese children from 2017 to 2019. Their antimicrobial susceptibilities were assessed using both E-test strips and Kirby-Bauer (KB) disk diffusion methods. Seven virulence-related genes (ptxA, ptxC, ptxP, prn, fim2, fim3, and tcfA2) and the A2047G mutation in the 23S rRNA gene were detected by PCR. Based on the susceptibilities and genotypes, 50 isolates were selected for multi-locus variable-number tandem-repeat analysis (MLVA) typing and whole-genome sequencing.

Results

A total of 311 B. pertussis strains were isolated from children with a median age of 4 months (interquartile range: 2–9 months). Strains carrying the ptxP1 allele were more frequent (84.9%, 264/311), were always ER (except for one strain), and were mainly related to ptxA1/ptxC1/prn1 alleles (99.6%, 263/264). The remaining 47 (15.1%) strains carried the ptxP3 allele, mainly harboring the ptxA1/ptxC2/prn2 alleles (93.6%, 44/47), and were sensitive to erythromycin (except for two strains). The two ER-ptxP3 isolates were first identified in China, belonged to MT27 and MT28 according to MLVA, and were classified into sub-lineage IVd by phylogenetic analysis of their genome sequences. This sub-lineage also includes many strains carrying the ptxP3 allele spreading in developed countries. For each tested antimicrobial, the susceptibilities judged by KB disks were consistent with those determined by E-test strips.

Conclusion

The present results reveal that B. pertussis strains with the ptxP1-ER profile still dominate in China, and a few strains carrying the ptxP3 allele have acquired the A2047G mutation in the 23S rRNA gene and the ER phenotype. The surveillance of the drug susceptibility of B. pertussis is necessary for all countries, and the KB disk method can be adopted as a screening test.

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.

Macrolide susceptibility and molecular characteristics of Bordetella pertussis

Objective

To analyse macrolide resistance and molecular characteristics of Bordetella pertussis clinical isolates from western China, and to explore the relationship between macrolide-resistance and genotypes.

Methods

Susceptibilities of B. pertussis clinical isolates to erythromycin, azithromycin and clarithromycin were determined by epsilometer test (E-test). Isolated strains were sequenced to ascertain the presence of the 23S rRNA gene A2047G mutation. Strains were typed using multilocus antigen sequence typing, multilocus variable-number tandem-repeat analysis (MLVA) and pulsed-field gel electrophoresis (PFGE).

Results

Of 58 B. pertussis strains isolated in this study, 46 were macrolide-resistant and 12 were macrolide sensitive. All macrolide-resistant strains carried the A2047G mutation and were the prn1/ptxP1/ptxA1/fim3-1/fim2-1 genotype; the MLVA types were MT195 (19/58), MT55 (13/58) and MT104 (14/58), and the PFGE profiles were classified into BpSR23 (17/58) and BpFINR9 (29/58) types. None of the macrolide-sensitive strains carried the A2047G mutation; genotypes were (prn9 or prn2)/ptxP3/ptxA1/fim3-1/fim2-1, and all were MT27. PFGE profiles differed from the macrolide-resistant strains.

Conclusions

B. pertussis clinical isolates from western China were severely resistant to macrolides. Genotypes differed between macrolide-resistant and macrolide-sensitive strains, and there may be a correlation between acquisition of macrolide resistance and changes in specific molecular types.

Comparative genomics of Bordetella pertussis isolates from New Zealand, a country with an uncommonly high incidence of whooping cough

Whooping cough, the respiratory disease caused by Bordetella pertussis, has undergone a wide-spread resurgence over the last several decades. Previously, we developed a pipeline to assemble the repetitive B. pertussis genome into closed sequences using hybrid nanopore and Illumina sequencing. Here, this sequencing pipeline was used to conduct a more high-throughput, longitudinal screen of 66 strains isolated between 1982 and 2018 in New Zealand. New Zealand has a higher incidence of whooping cough than many other countries; usually at least twice as many cases per 100000 people as the USA and UK and often even higher, despite similar rates of vaccine uptake. To the best of our knowledge, these strains are the first New Zealand B. pertussis isolates to be sequenced. The analyses here show that, on the whole, genomic trends in New Zealand B. pertussis isolates, such as changing allelic profile in vaccine-related genes and increasing pertactin deficiency, have paralleled those seen elsewhere in the world. At the same time, phylogenetic comparisons of the New Zealand isolates with global isolates suggest that a number of strains are circulating in New Zealand, which cluster separately from other global strains, but which are closely related to each other. The results of this study add to a growing body of knowledge regarding recent changes to the B. pertussis genome, and are the first genetic investigation into B. pertussis isolates from New Zealand.

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.

Genetic Diversity of Meningococcal Serogroup B Vaccine Antigens among Carriage Isolates Collected from Students at Three Universities in the United States, 2015-2016

Carriage evaluations were conducted during 2015 to 2016 at two U.S. universities in conjunction with the response to disease outbreaks caused by Neisseria meningitidis serogroup B and at a university where outbreak and response activities had not occurred. All eligible students at the two universities received the serogroup B meningococcal factor H binding protein vaccine (MenB-FHbp); 5.2% of students (181/3,509) at one university received MenB-4C. A total of 1,514 meningococcal carriage isolates were obtained from 8,905 oropharyngeal swabs from 7,001 unique participants. Whole-genome sequencing data were analyzed to understand MenB-FHbp’s impact on carriage and antigen genetic diversity and distribution. Of 1,422 isolates from carriers with known vaccination status (726 [51.0%] from MenB-FHbp-vaccinated, 42 [3.0%] from MenB-4C-vaccinated, and 654 [46.0%] from unvaccinated participants), 1,406 (98.9%) had intact fHbp alleles (716 from MenB-FHbp-vaccinated participants). Of 726 isolates from MenB-FHbp-vaccinated participants, 250 (34.4%) harbored FHbp peptides that may be covered by MenB-FHbp. Genogroup B was detected in 122/1,422 (8.6%) and 112/1,422 (7.9%) isolates from MenB-FHbp-vaccinated and unvaccinated participants, respectively. FHbp subfamily and peptide distributions between MenB-FHbp-vaccinated and unvaccinated participants were not statistically different. Eighteen of 161 MenB-FHbp-vaccinated repeat carriers (11.2%) acquired a new strain containing one or more new vaccine antigen peptides during multiple rounds of sample collection, which was not statistically different (P = 0.3176) from the unvaccinated repeat carriers (1/30; 3.3%). Our findings suggest that lack of MenB vaccine impact on carriage was not due to missing the intact fHbp gene; MenB-FHbp did not affect antigen genetic diversity and distribution during the study period.

Genomic Surveillance and Improved Molecular Typing of Bordetella pertussis Using wgMLST

Multilocus sequence typing (MLST) provides allele-based characterization of bacterial pathogens in a standardized framework. However, classical MLST schemes for Bordetella pertussis, the causative agent of whooping cough, seldom reveal diversity among the small number of gene targets and thereby fail to delineate population structure. To improve the discriminatory power of allele-based molecular typing of B. pertussis, we have developed a whole-genome MLST (wgMLST) scheme from 225 reference-quality genome assemblies. Iterative refinement and allele curation resulted in a scheme of 3,506 coding sequences and covering 81.4% of the B. pertussis genome. This wgMLST scheme was further evaluated with data from a convenience sample of 2,389 B. pertussis isolates sequenced on Illumina instruments, including isolates from known outbreaks and epidemics previously characterized by existing molecular assays, as well as replicates collected from individual patients. wgMLST demonstrated concordance with whole-genome single nucleotide polymorphism (SNP) profiles, accurately resolved outbreak and sporadic cases in a retrospective comparison, and clustered replicate isolates collected from individual patients during diagnostic confirmation. Additionally, a reanalysis of isolates from two statewide epidemics using wgMLST reconstructed the population structures of circulating strains with increased resolution, revealing new clusters of related cases. Comparison with an existing core genome (cgMLST) scheme highlights the stable gene content of this bacterium and forms the initial foundation for necessary standardization. These results demonstrate the utility of wgMLST for improving B. pertussis characterization and genomic surveillance during the current pertussis disease resurgence.

Rapid and simple SNP genotyping for Bordetella pertussis epidemic strain MT27 based on a multiplexed single-base extension assay

Multilocus variable-number tandem repeat analysis (MLVA) is widely used for genotyping of Bordetella pertussis, the causative bacteria for pertussis. However, MLVA genotyping is losing its discriminate power because prevalence of the epidemic MT27 strain (MLVA-27) is increasing worldwide. To address this, we developed a single nucleotide polymorphism (SNP) genotyping method for MT27 based on multiplexed single-base extension (SBE) assay. A total of 237 MT27 isolates collected in Japan during 1999–2018 were genotyped and classified into ten SNP genotypes (SG1 to SG10) with a Simpson’s diversity index (DI) of 0.79 (95% CI 0.76–0.82). Temporal trends showed a marked increase in the genotypic diversity in the 2010s: Simpson’s DI was zero in 1999–2004, 0.16 in 2005–2009, 0.83 in 2010–2014, and 0.76 in 2015–2018. This indicates that the SNP genotyping is applicable to the recently circulating MT27 strain. Additionally, almost all outbreak-associated MT27 isolates were classified into the same SNP genotypes for each outbreak. Multiplexed SBE assay allows for rapid and simple genotyping, indicating that the SNP genotyping can potentially be a useful tool for subtyping the B. pertussis MT27 strain in routine surveillance and outbreak investigations.

Nonsynonymous Polymorphism Counts in Bacterial Genomes: a Comparative Examination

Genomic data reveal single-nucleotide polymorphisms (SNPs) that may carry information about the evolutionary history of bacteria. However, it remains unclear what inferences about selection can be made from genomic SNP data. Bacterial species are often sampled during epidemic outbreaks or within hosts during the course of chronic infections. SNPs obtained from genomic analysis of these data are not necessarily fixed. Treating them as fixed during analysis by using measures such as the ratio of nonsynonymous to synonymous evolutionary changes (dN/dS) may lead to incorrect inferences about the strength and direction of selection. In this study, we consider data from a range of whole-genome sequencing studies of bacterial pathogens and explore patterns of nonsynonymous variation to assess whether evidence of selection can be identified by investigating SNP counts alone across multiple WGS studies. We visualize these SNP data in ways that highlight their relationship to neutral baseline expectations. These neutral expectations are based on a simple model of mutation, from which we simulate SNP accumulation to investigate how SNP counts are distributed under alternative assumptions about positive and negative selection. We compare these patterns with empirical SNP data and illustrate the general difficulty of detecting positive selection from SNP data. Finally, we consider whether SNP counts observed at the between-host population level differ from those observed at the within-host level and find some evidence that suggests that dynamics across these two scales are driven by different underlying processes.IMPORTANCE Identifying selection from SNP data obtained from whole-genome sequencing studies is challenging. Some current measures used to identify and quantify selection acting on genomes rely on fixed differences; thus, these are inappropriate for SNP data where variants are not fixed. With the increase in whole-genome sequencing studies, it is important to consider SNP data in the context of evolutionary processes. How SNPs are counted and analyzed can help in understanding mutation accumulation and trajectories of strains. We developed a tool for identifying possible evidence of selection and for comparative analysis with other SNP data. We propose a model that provides a rule-of-thumb guideline and two new visualization techniques that can be used to interpret and compare SNP data. We quantify the expected proportion of nonsynonymous SNPs in coding regions under neutrality and demonstrate its use in identifying evidence of positive and negative selection from simulations and empirical data.

The epidemic of erythromycin-resistant Bordetella pertussis with limited genome variation associated with pertussis resurgence in China

BACKGROUND

The resurgence of Bordetella pertussis infections leading to whooping cough is a concern in many parts of the world. The number of pertussis cases in China has increased significantly since 2013.

RESEARCH DESIGN AND METHODS

In this study, whole-genome sequencing analysis was performed for 388 B. pertussis strains isolated in China from the 1970s to 2018, combining 594 published strains from around the world.

RESULTS

This study revealed that lineage V diverged about 50 years ago in China, while lineage IV is dominant in the other countries. It also revealed that the erythromycin-resistant sub-lineages Va, Vb, and Vc with limited genomic variation emerged 11 ~ 12 years ago. These three sub-lineages were identified after the co-purified acellular vaccines (cp-ACVs) completely replaced the previous whole cell vaccines (WCVs) after the national immunization program of 2012. It suggests that the cp-ACVs cannot induce immunity that is potent enough to restrict the spread of the lineage V, antibiotic abuse further favors the spread of this lineage in China.

CONCLUSIONS

These findings demand a reassessment of the immunization strategy and development of new vaccines in China to stop the resurgence and drug resistance of B. pertussis.

Acellular Pertussis Vaccines Induce Anti-pertactin Bactericidal Antibodies Which Drives the Emergence of Pertactin-Negative Strains

Despite high vaccination coverage, Bordetella pertussis the causative agent of whooping cough is still a health concern worldwide. A resurgence of pertussis cases has been reported, particularly in countries using acellular vaccines with waning immunity and pathogen adaptation thought to be responsible. A better understanding of protective immune responses is needed for the development of improved vaccines. In our study, B. pertussis strain B1917 variants presenting a single gene deletion were generated to analyze the role of vaccine components or candidate vaccine antigens as targets for bactericidal antibodies generated after acellular vaccination or natural infection. Our results show that acellular vaccination generates bactericidal antibodies that are only directed against pertactin. Serum bactericidal assay performed with convalescent samples show that disease induces bactericidal antibodies against Prn but against other antigen(s) as well. Four candidate vaccine antigens (CyaA, Vag8, BrkA, and TcfA) have been studied but were not targets for complement-mediated bactericidal antibodies after natural infection. We confirm that Vag8 and BrkA are involved in complement resistance and would be targeted by blocking antibodies. Our study suggests that the emergence and the widespread circulation of Prn-deficient strains is driven by acellular vaccination and the generation of bactericidal antibodies targeting Prn.

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.

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

Pertussis caused by Bordetella pertussis, remains a public health problem worldwide, despite high vaccine coverage in infants and children in many countries. Iran has been using whole cell vaccine for the last 50 years with more than 95% vaccination rate since 1988 and has experienced pertussis resurgence in recent years. Here, we sequenced 55 B. pertussis isolates mostly collected from three provinces with the highest number of pertussis cases in Iran, including Tehran, Mazandaran, and Eastern-Azarbayjan from the period of 2008-2016. Most isolates carried ptxP3/prn2 alleles (42/55, 76%), the same genotype as isolates circulating in acellular vaccine-administrating countries. The second most frequent genotype was ptxP3/prn9 (8/55, 14%). Only three isolates (5%) were ptxP1. Phylogenetic analysis showed that Iranian ptxP3 isolates can be divided into eight clades (Clades 1-8) with no temporal association. Most of the isolates from Tehran grouped together as one distinctive clade (Clade 8) with six unique single nucleotide polymorphisms (SNPs). In addition, the prn9 isolates were grouped together as Clade 5 with 12 clade-supporting SNPs. No pertactin deficient isolates were found among the 55 Iranian isolates. Our findings suggest that there is an ongoing adaptation and evolution of B. pertussis regardless of the types of vaccine used.

Status and potential of bacterial genomics for public health practice: a scoping review

BACKGROUND

Next-generation sequencing (NGS) is increasingly being translated into routine public health practice, affecting the surveillance and control of many pathogens. The purpose of this scoping review is to identify and characterize the recent literature concerning the application of bacterial pathogen genomics for public health practice and to assess the added value, challenges, and needs related to its implementation from an epidemiologist's perspective.

METHODS

In this scoping review, a systematic PubMed search with forward and backward snowballing was performed to identify manuscripts in English published between January 2015 and September 2018. Included studies had to describe the application of NGS on bacterial isolates within a public health setting. The studied pathogen, year of publication, country, number of isolates, sampling fraction, setting, public health application, study aim, level of implementation, time orientation of the NGS analyses, and key findings were extracted from each study. Due to a large heterogeneity of settings, applications, pathogens, and study measurements, a descriptive narrative synthesis of the eligible studies was performed.

RESULTS

Out of the 275 included articles, 164 were outbreak investigations, 70 focused on strategy-oriented surveillance, and 41 on control-oriented surveillance. Main applications included the use of whole-genome sequencing (WGS) data for (1) source tracing, (2) early outbreak detection, (3) unraveling transmission dynamics, (4) monitoring drug resistance, (5) detecting cross-border transmission events, (6) identifying the emergence of strains with enhanced virulence or zoonotic potential, and (7) assessing the impact of prevention and control programs. The superior resolution over conventional typing methods to infer transmission routes was reported as an added value, as well as the ability to simultaneously characterize the resistome and virulome of the studied pathogen. However, the full potential of pathogen genomics can only be reached through its integration with high-quality contextual data.

CONCLUSIONS

For several pathogens, it is time for a shift from proof-of-concept studies to routine use of WGS during outbreak investigations and surveillance activities. However, some implementation challenges from the epidemiologist's perspective remain, such as data integration, quality of contextual data, sampling strategies, and meaningful interpretations. Interdisciplinary, inter-sectoral, and international collaborations are key for an appropriate genomics-informed surveillance.

Genomic epidemiology of erythromycin-resistant Bordetella pertussis in China

Macrolides such as erythromycin are the empirical treatment of Bordetella pertussis infections. China has experienced an increase in erythromycin-resistant B. pertussis isolates since they were first reported in 2013. Here, we undertook a genomic study on Chinese B. pertussis isolates from 2012 to 2015 to elucidate the origins and phylogenetic relationships of erythromycin-resistant B. pertussis isolates in China. A total of 167 Chinese B. pertussis isolates were used for antibiotic sensitivity testing and multiple locus variable-number tandem repeat (VNTR) analysis (MLVA). All except four isolates were erythromycin-resistant and of the four erythromycin-sensitive isolates, three were non-ptxP1. MLVA types (MT), MT55, MT104 and MT195 were the predominant types. Fifty of those isolates were used for whole genome sequencing and phylogenetic analysis. Genome sequencing and phylogenetic analysis revealed three independent erythromycin-resistant lineages and all resistant isolates carried a mutation in the 23S rRNA gene. A novel fhaB3 allele was found uniquely in Chinese ptxP1 isolates and these Chinese ptxP1-ptxA1-fhaB3 had a 5-fold higher mutation rate than the global ptxP1-ptxA1 B. pertussis population. Our results suggest that the evolution of Chinese B. pertussis is likely to be driven by selection pressure from both vaccination and antibiotics. The emergence of the new non-vaccine fhaB3 allele in Chinese B. pertussis population may be a result of selection from vaccination, whereas the expansion of ptxP1-fhaB3 lineages was most likely to be the result of selection pressure from antibiotics. Further monitoring of B. pertussis in China is required to better understand the evolution of the pathogen.

Genome characteristics of Bordetella pertussis isolates from Tunisia

The recent increase in pertussis cases observed in some countries may have several causes, including the evolution of Bordetella pertussis populations towards escape of vaccine-induced immunity. Most genomic studies of B. pertussis isolates performed so far are from countries that use acellular vaccines. The objective was to analyse genomic sequences of isolates collected during the 2014 whooping cough epidemic in Tunisia, a country where whole-cell vaccines are used. Ten Tunisian isolates and four vaccine strains were sequenced and compared to 169 isolates from countries where acellular vaccines are used. Phylogenetic analysis showed that Tunisian isolates are diverse, demonstrating a multi-strain 2014 epidemic peak, and are intermixed with those circulating in other world regions, showing inter-country transmission. Consistently, Tunisian isolates have antigen variant composition observed in other world regions. No pertactin-deficient strain was observed. The Tunisian B. pertussis population appears to be largely connected with populations from other countries.

Proteomics of Bordetella pertussis whole-cell and acellular vaccines

OBJECTIVES

Bordetella pertussis is the etiological agent of whooping cough, a bacterial infection of especially children, which may be fatal without treatment. In frame of studies to investigate putative effects of vaccination on host-pathogen interaction and clonal distribution of strains, in addition to Corynebacterium diphtheriae and Clostridium tetani toxoid vaccines, also whole-cell and acellular pertussis vaccines were analyzed by mass spectrometry.

DATA DESCRIPTION

LC-MS/MS spectra were generated and analyzed using B. pertussis genome data and proteins present in whole-cell and acellular pertussis vaccines were identified. Subcellular localization of proteins and presence of signal peptides was determined bioinformatically.

Recombinant outer membrane protein Q and putative lipoprotein from Bordetella pertussis inducing strong humoral response were not protective alone in the murine lung colonization model

Despite high vaccination coverage after introduction of whole cell (wP) and acellular pertussis (aP) vaccines, pertussis resurgence has been reported in many countries. aP vaccines are commonly preferred due to side effects of wP vaccines and formulated with aluminum hydroxide (Alum), which is not an effective adjuvant to eliminate Bordetella pertussis. Low efficiency of current aP vaccines is thought to be the main reason for the resurgence for which newer generation aP vaccines are needed. In the present study, immunogenicity and protective efficacy of outer membrane protein Q (OmpQ) and a putative lipoprotein (Lpp) from B. pertussis were investigated in mice by using two diefrent adjuvants, monophosphoryl lipid A (MPLA) or Alum. OmpQ and putative Lpp were cloned, expressed, and purified from Escherichia coli. The proteins were formulated to immunize mice. Both recombinant OmpQ and putative Lpp induced a significant increase in immunoglobulin G1 (IgG1) and immunoglobulin G2a (IgG2a) responses compared to the control group. Moreover, MPLA-adjuvanted formulations resulted in higher IgG2a levels than Alum-adjuvanted ones. However, there were no significant differences between test and control groups regarding interferon-gamma (IFN-γ) levels, and the mice lung colonization experiments indicated that neither rOmpQ nor rLpp could confer protection alone against B. pertussis challenge.

How Genomics Is Changing What We Know About the Evolution and Genome of Bordetella pertussis

The evolution of Bordetella pertussis from a common ancestor similar to Bordetella bronchiseptica has occurred through large-scale gene loss, inactivation and rearrangements, largely driven by the spread of insertion sequence element repeats throughout the genome. B. pertussis is widely considered to be monomorphic, and recent evolution of the B. pertussis genome appears to, at least in part, be driven by vaccine-based selection. Given the recent global resurgence of whooping cough despite the wide-spread use of vaccination, a more thorough understanding of B. pertussis genomics could be highly informative. In this chapter we discuss the evolution of B. pertussis, including how vaccination is changing the circulating B. pertussis population at the gene-level, and how new sequencing technologies are revealing previously unknown levels of inter- and intra-strain variation at the genome-level.

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.

Prevention of pertussis: An unresolved problem

Pertussis is a highly contagious respiratory disease caused by Bordetella pertussis. However, after the introduction of the whole-cell pertussis vaccine (wP), the annual incidence rates of the disease progressively declined. Despite this result, the inclusion of wP in the national immunization schedule of infants and young children was debated regarding its safety. Several efforts to produce vaccines based on B. pertussis components capable of evoking protective immunity with no or limited adverse events were made. Of these others, five pertussis antigens were considered possible components of acellular vaccines (aPs): pertussis toxin (PT), filamentous haemagglutinin (FHA), pertactin (PRN) and fimbria proteins 2 and 3. However, the introduction of aPs was followed by a slight but progressive increase in the incidence of pertussis. This paper discusses the potential reasons for reduced aPs efficacy. Moreover, it attempts to evaluate the real effectiveness of aPs and the potential differences between available preparations. Data analysis showed that several boosters are needed to maintain protection against pertussis and additional studies are needed to confirm the antigens that should be included in aPs to improve the prevention of pertussis.

In Search of Factors Negatively Affecting Vaccine Immunity to Pertussis in Preschool Children Before the Administration of the First Booster

INTRODUCTION

The top priority for active immunoprophylaxis of pertussis is the immunisation of infants as they can sometimes develop severe multiple-organ complications.

OBJECTIVES

The aim of the work is the identification of factors negatively affecting vaccine immunity to pertussis in preschool children prior to the administration of the first booster.

PATIENTS AND METHODS

The research was conducted on 352 children from 4.5 to 5.9 years of age who were hospitalised in the University Children's Hospital in Lublin (Poland) from 1 January 2012 to 31 December 2015. The children taking part in the study had been administered all the mandatory vaccines from their birth to the age of 2 or 2.5 years old according to the Polish Immunisation Program 2008⁻2009. The immunoenzymatic method ELISA (enzyme-linked immunosorbent assay) was applied to assess vaccine immunity to tetanus, diphtheria, pertussis, Haemophilus influenzae type b (Hib), poliomyelitis (IPV), mumps, rubella and measles. The level of vaccine antibodies to hepatitis type B was determined chemilumiscently.

RESULTS

The protective antibody titre was not found in 41 (11.65%) children before the administration of the booster. To verify the collective impact of parameters analysed on antibody titre to pertussis, the Generalized Linear Model (GLZ) was used. Gender, type of vaccine, asthma, Hib and mumps antibody titres have been shown to be predictors of vaccine immunity to pertussis.

CONCLUSIONS

Immunomodulation considered on the example of titre of IgG antibody to pertussis can serve as a useful model of the assessment of development of acquired immunity after mandatory vaccinations.

Primary transcriptome analysis reveals importance of IS elements for the shaping of the transcriptional landscape of Bordetella pertussis

Bordetella pertussis is the causative agent of whooping cough, a respiratory disease still considered as a major public health threat and for which recent re-emergence has been observed. Constant reshuffling of Bordetella pertussis genome organization was observed during evolution. These rearrangements are essentially mediated by Insertion Sequences (IS), a mobile genetic elements present in more than 230 copies in the genome, which are supposed to be one of the driving forces enabling the pathogen to escape from vaccine-induced immunity. Here we use high-throughput sequencing approaches (RNA-seq and differential RNA-seq), to decipher Bordetella pertussis transcriptome characteristics and to evaluate the impact of IS elements on transcriptome architecture. Transcriptional organization was determined by identification of transcription start sites and revealed also a large variety of non-coding RNAs including sRNAs, leaderless mRNAs or long 3' and 5'UTR including seven riboswitches. Unusual topological organizations, such as overlapping 5'- or 3'-extremities between oppositely orientated mRNA were also unveiled. The pivotal role of IS elements in the transcriptome architecture and their effect on the transcription of neighboring genes was examined. This effect is mediated by the introduction of IS harbored promoters or by emergence of hybrid promoters. This study revealed that in addition to their impact on genome rearrangements, most of the IS also impact on the expression of their flanking genes. Furthermore, the transcripts produced by IS are strain-specific due to the strain to strain variation in IS copy number and genomic context.

Bordetella pertussis population dynamics and phylogeny in Japan after adoption of acellular pertussis vaccines

Bordetella pertussis, the causative agent of whooping cough, has experienced a resurgence in the past 15 years, despite the existence of both whole-cell and acellular vaccines. Here, we performed whole genome sequencing analysis of 149 clinical strains, provided by the National Institute of Infectious Diseases (NIID), Japan, isolated in 1982–2014, after Japan became the first country to adopt acellular vaccines against B. pertussis. Additionally, we sequenced 39 strains provided by the Konan Kosei Hospital in Aichi prefecture, Japan, isolated in 2008–2013. The genome sequences afforded insight into B. pertussis genome variability and population dynamics in Japan, and revealed that the B. pertussis population in Japan was characterized by two major clades that divided more than 40 years ago. The pertactin gene was disrupted in about 20 % of the 149 NIID isolates, by either a deletion within the signal sequence (ΔSS) or the insertion of IS element IS481 (prn :: IS481). Phylogeny suggests that the parent clones for these isolates originated in Japan. Divergence dating traced the first generation of the pertactin-deficient mutants in Japan to around 1990, and indicated that strains containing the alternative pertactin allele prn2 may have appeared in Japan around 1974. Molecular clock data suggested that observed fluctuations in B. pertussis population size may have coincided with changes in vaccine usage in the country. The continuing failure to eradicate the disease warrants an exploration of novel vaccine compositions.

Screening and Genomic Characterization of Filamentous Hemagglutinin-Deficient Bordetella pertussis

Despite high vaccine coverage, pertussis cases in the United States have increased over the last decade. Growing evidence suggests that disease resurgence results, in part, from genetic divergence of circulating strain populations away from vaccine references. The United States employs acellular vaccines exclusively, and current Bordetella pertussis isolates are predominantly deficient in at least one immunogen, pertactin (Prn). First detected in the United States retrospectively in a 1994 isolate, the rapid spread of Prn deficiency is likely vaccine driven, raising concerns about whether other acellular vaccine immunogens experience similar pressures, as further antigenic changes could potentially threaten vaccine efficacy. We developed an electrochemiluminescent antibody capture assay to monitor the production of the acellular vaccine immunogen filamentous hemagglutinin (Fha). Screening 722 U.S. surveillance isolates collected from 2010 to 2016 identified two that were both Prn and Fha deficient. Three additional Fha-deficient laboratory strains were also identified from a historic collection of 65 isolates dating back to 1935. Whole-genome sequencing of deficient isolates revealed putative, underlying genetic changes. Only four isolates harbored mutations to known genes involved in Fha production, highlighting the complexity of its regulation. The chromosomes of two Fha-deficient isolates included unexpected structural variation that did not appear to influence Fha production. Furthermore, insertion sequence disruption of fhaB was also detected in a previously identified pertussis toxin-deficient isolate that still produced normal levels of Fha. These results demonstrate the genetic potential for additional vaccine immunogen deficiency and underscore the importance of continued surveillance of circulating B. pertussis evolution in response to vaccine pressure.

Genomic analysis of endemic clones of toxigenic and non-toxigenic Corynebacterium diphtheriae in Belarus during and after the major epidemic in 1990s

BACKGROUND

Diphtheria remains a major public health concern with multiple recent outbreaks around the world. Moreover, invasive non-toxigenic strains have emerged globally causing severe infections. A diphtheria epidemic in the former Soviet Union in the 1990s resulted in ~5000 deaths. In this study, we analysed the genome sequences of a collection of 93 C. diphtheriae strains collected during and after this outbreak (1996 - 2014) in a former Soviet State, Belarus to understand the evolutionary dynamics and virulence capacities of these strains.

RESULTS

C. diphtheriae strains from Belarus belong to ten sequence types (STs). Two major clones, non-toxigenic ST5 and toxigenic ST8, encompassed 76% of the isolates that are associated with sore throat and diphtheria in patients, respectively. Core genomic diversity is limited within outbreak-associated ST8 with relatively higher mutation rates (8.9 × 10-7 substitutions per strain per year) than ST5 (5.6 × 10-7 substitutions per strain per year) where most of the diversity was introduced by recombination. A variation in the virulence gene repertoire including the presence of tox gene is likely responsible for pathogenic differences between different strains. However, strains with similar virulence potential can cause disease in some individuals and remain asymptomatic in others. Eight synonymous single nucleotide polymorphisms were observed between the tox genes of the vaccine strain PW8 and other toxigenic strains of ST8, ST25, ST28, ST41 and non-toxigenic tox gene-bearing (NTTB) ST40 strains. A single nucleotide deletion at position 52 in the tox gene resulted in the frameshift in ST40 isolates, converting them into NTTB strains.

CONCLUSIONS

Non-toxigenic C. diphtheriae ST5 and toxigenic ST8 strains have been endemic in Belarus both during and after the epidemic in 1990s. A high vaccine coverage has effectively controlled diphtheria in Belarus; however, non-toxigenic strains continue to circulate in the population. Recombination is an important evolutionary force in shaping the genomic diversity in C. diphtheriae. However, the relative role of recombination and mutations in diversification varies between different clones.

Bordetella pertussis: new concepts in pathogenesis and treatment

PURPOSE OF REVIEW

The purpose of this review is to summarize and discuss recent findings and selected topics of interest in Bordetella pertussis virulence and pathogenesis and treatment of pertussis. It is not intended to cover issues on immune responses to B. pertussis infection or problems with currently used pertussis vaccines.

RECENT FINDINGS

Studies on the activities of various B. pertussis virulence factors include the immunomodulatory activities of filamentous hemagglutinin, fimbriae, and adenylate cyclase toxin. Recently emerging B. pertussis strains show evidence of genetic selection for vaccine escape mutants, with changes in vaccine antigen-expressing genes, some of which may have increased the virulence of this pathogen. Severe and fatal pertussis in young infants continues to be a problem, with several studies highlighting predictors of fatality, including the extreme leukocytosis associated with this infection. Treatments for pertussis are extremely limited, though early antibiotic intervention may be beneficial. Neutralizing pertussis toxin activity may be an effective strategy, as well as targeting two host proteins, pendrin and sphingosine-1-phosphate receptors, as novel potential therapeutic interventions.

SUMMARY

Pertussis is reemerging as a major public health problem and continued basic research is revealing information on bacterial virulence and disease pathogenesis, as well as potential novel strategies for vaccination and targets for therapeutic intervention.

Multi-locus variable-number tandem repeat analysis of Bordetella pertussis isolates circulating in Poland in the period 1959-2013

PURPOSE

Despite the long history of pertussis vaccination and high vaccination coverage in Poland and many other developed countries, pertussis incidence rates have increased substantially, making whooping cough one of the most prevalent vaccine-preventable diseases. Among the factors potentially involved in pertussis resurgence, the adaptation of the Bordetella pertussis population to country-specific vaccine-induced immunity through selection of non-vaccine-type strains still needs detailed studies.

METHODOLOGY

Multi-locus variable-number tandem repeat analysis (MLVA), also linked to MLST and PFGE profiling, was applied to trace the genetic changes in the B. pertussis population circulating in Poland in the period 1959-2013 versus country-specific vaccine strains.

RESULTS

Generally, among 174 B. pertussis isolates, 31 MLVA types were detected, of which 11 were not described previously. The predominant MLVA types of recent isolates in Poland were different from those of the typical isolates circulating in other European countries. The MT27 type, currently predominant in Europe, was rarely seen and detected in only five isolates among all studied. The features of the vaccine strains used for production of the pertussis component of a national whole-cell diphtheria-tetanus-pertussis (DTP) vaccine, as studied by MLVA and MLST tools, were found to not match those observed in the currently circulating B. pertussis isolates in Poland.

CONCLUSIONS

Differences traced by MLVA in relation to the MLST and PFGE profiling confirmed that the B. pertussis strain types currently observed elsewhere in Europe, even if appearing in Poland, were not able to successfully disseminate within a human population in Poland that has been vaccinated with a whole-cell pertussis vaccine not used in other countries.

The History of Bordetella pertussis Genome Evolution Includes Structural Rearrangement

Despite high pertussis vaccine coverage, reported cases of whooping cough (pertussis) have increased over the last decade in the United States and other developed countries. Although Bordetella pertussis is well known for its limited gene sequence variation, recent advances in long-read sequencing technology have begun to reveal genomic structural heterogeneity among otherwise indistinguishable isolates, even within geographically or temporally defined epidemics. We have compared rearrangements among complete genome assemblies from 257 B. pertussis isolates to examine the potential evolution of the chromosomal structure in a pathogen with minimal gene nucleotide sequence diversity. Discrete changes in gene order were identified that differentiated genomes from vaccine reference strains and clinical isolates of various genotypes, frequently along phylogenetic boundaries defined by single nucleotide polymorphisms. The observed rearrangements were primarily large inversions centered on the replication origin or terminus and flanked by IS481, a mobile genetic element with >240 copies per genome and previously suspected to mediate rearrangements and deletions by homologous recombination. These data illustrate that structural genome evolution in B. pertussis is not limited to reduction but also includes rearrangement. Therefore, although genomes of clinical isolates are structurally diverse, specific changes in gene order are conserved, perhaps due to positive selection, providing novel information for investigating disease resurgence and molecular epidemiology.IMPORTANCE Whooping cough, primarily caused by Bordetella pertussis, has resurged in the United States even though the coverage with pertussis-containing vaccines remains high. The rise in reported cases has included increased disease rates among all vaccinated age groups, provoking questions about the pathogen's evolution. The chromosome of B. pertussis includes a large number of repetitive mobile genetic elements that obstruct genome analysis. However, these mobile elements facilitate large rearrangements that alter the order and orientation of essential protein-encoding genes, which otherwise exhibit little nucleotide sequence diversity. By comparing the complete genome assemblies from 257 isolates, we show that specific rearrangements have been conserved throughout recent evolutionary history, perhaps by eliciting changes in gene expression, which may also provide useful information for molecular epidemiology.

Complete Genome Sequences of Four Bordetella pertussis Vaccine Reference Strains from Serum Institute of India

Serum Institute of India is among the world’s largest vaccine producers. Here, we report the complete genome sequences for four Bordetella pertussis strains used by Serum Institute of India in the production of whole-cell pertussis vaccines.

Complete Genome Sequences of Bordetella pertussis Vaccine Reference Strains 134 and 10536

Vaccine formulations and vaccination programs against whooping cough (pertussis) vary worldwide. Here, we report the complete genome sequences of two divergent Bordetella pertussis reference strains used in the production of pertussis vaccines.

Hajj vaccinations-facts, challenges, and hope

Vaccination is an effective preventive measure that has been used in the unique Hajj pilgrimage setting to control the transmission of infectious diseases. The current vaccination policy applied during Hajj is reviewed herein, highlighting the effectiveness of the approaches applied and identifying research gaps that need to be filled in order to improve the development and dissemination of Hajj vaccination strategies.

Genomic dissection of Australian Bordetella pertussis isolates from the 2008-2012 epidemic

OBJECTIVES

Despite high pertussis vaccination coverage, Australia experienced a prolonged epidemic in 2008-2012. The predominant Bordetella pertussis genotype harboured pertussis toxin promoter allele, ptxP3, and pertactin gene allele, prn2. The emergence and expansion of prn non-expressing isolates (Prn negative), were also observed. We aimed to investigate the microevolution and genomic diversity of epidemic B. pertussis isolates.

METHODS

We sequenced 22 B. pertussis isolates collected in 2008-2012 from two states of Australia which are geographically widely separated. Ten of the 22 were Prn negative isolates with three different modes of silencing of prn (prn::IS481F, prn::IS481R and prn::IS1002). Five pre-epidemic isolates were also sequenced for comparison.

RESULTS

Five single nucleotide polymorphisms were common in the epidemic isolates and differentiated them from pre-epidemic isolates. The Australian epidemic isolates can be divided into five lineages (EL1-EL5) with EL1 containing only Prn negative isolates. Comparison with global isolates showed that three lineages remained geographically and temporally distinct whereas two lineages mixed with isolates from 2012 UK outbreak.

CONCLUSION

Our results suggest significant diversification and the microevolution of B. pertussis within the 2008-2012 Australian epidemic.