Screening and Genomic Characterization of Filamentous Hemagglutinin-Deficient Bordetella pertussis

Pertactin deficiency of Bordetella pertussis: Insights into epidemiology, and perspectives on surveillance and public health impact

Pertussis resurgence has been documented even in countries with high pediatric vaccine coverage. The proportion of Bordetella pertussis isolates not expressing pertactin (PRN) has increased in several countries where acellular pertussis (aP) vaccines are used. We systematically reviewed published literature up to July 2023 on PRN-negative B. pertussis isolates in MEDLINE and Embase with no geographical limitations, complemented with a gray literature search. An increase in the proportion of PRN-negative isolates was observed in countries where aP vaccines were used, while such isolates seem to be absent in countries using whole-cell pertussis vaccination. We reviewed the data supporting aP vaccine-driven evolution of B. pertussis, explored the effects of PRN deficiency on the clinical presentation of pertussis, summarized the evidence for preserved aP vaccine effectiveness, and proposed actions to further improve assessment of the clinical significance of PRN deficiency and its potential impact on pertussis prevention.

The seroepidemiology of immunoglobulin G antibodies against pertussis toxin and filamentous hemagglutinin in the east of China during the COVID-19 pandemic

This study employed sero-epidemiological methods to estimate the incidence of pertussis within a healthy population located in eastern China. The aim was to gain deeper insights into the epidemiological characteristics and burden of pertussis within the country. Blood samples were collected from healthy individuals in Jiangsu Province between June 2019 and December 2022. The levels of IgG antibodies against pertussis toxin (anti-PT) and filamentous hemagglutinin (anti-FHA) in the serum were quantitatively measured using enzyme-linked immunosorbent assay (ELISA). Additionally, pertussis case data reported in Jiangsu Province were collected from the China Information System for Disease Control and Prevention and compared with the results of this study. In 2022, the reported incidence of pertussis stood at 1.0 per 100,000 individuals, marking the highest rate observed in the past two decades. Among 1,909 patients examined, the geometric mean concentration (GMC) of anti-PT IgG antibody was 20.2 (18.5-21.9) IU/ml, while that of anti-FHA IgG antibody was 27.0 (25.4-28.7) IU/ml. The IgG-PT and IgG-FHA seropositivity rate (>20.0 IU/ml) was highest in the 1 ~ 2 y old group and decreased rapidly to the lowest in the 3 ~ 4 y old group and then increased gradually with age. The estimated rate of pertussis infection based on seroprevalence was approximately 25,625-fold higher than the reported notification rate in the ≥15 year age group. Our findings highlight decreased immunity post-vaccination, stressing the importance of additional booster shots for adolescents and adults to maintain immunity and reduce severe illness. Additionally, they offer vital guidance for policymakers to enhance immunization strategies.

Emergence of Erythromycin-Resistant and Pertactin- and Filamentous Hemagglutinin-Deficient Bordetella pertussis Strains - Beijing, China, 2022-2023

What is already known about this topic?

Pertussis has reemerged as a significant public health threat, primarily due to variations in Bordetella pertussis strains, antimicrobial resistance, and vaccine evasion.

What is added by this report?

All isolated strains were identified as ptxA1/ptxC2/ptxP3/prn150/fim2-1/fim3-1/fhaB1/tcfA2 type and exhibited resistance to erythromycin. Two strains showed a deficiency in Fha, thirty in Prn, and one strain exhibited multiple immunogen deficiencies.

What are the implications for public health practice?

The emergence and spread of immunogen-deficient strains likely result from prolonged vaccine selection pressure, posing challenges to the efficacy of pertussis vaccines. Additionally, the ongoing dissemination of ptxP3 strains with high-level macrolide resistance presents a significant obstacle to clinical treatment strategies.

Nanopore DNA sequencing technologies and their applications towards single-molecule proteomics

Sequencing of nucleic acids with nanopores has emerged as a powerful tool offering rapid readout, high accuracy, low cost and portability. This label-free method for sequencing at the single-molecule level is an achievement on its own. However, nanopores also show promise for the technologically even more challenging sequencing of polypeptides, something that could considerably benefit biological discovery, clinical diagnostics and homeland security, as current techniques lack portability and speed. Here we survey the biochemical innovations underpinning commercial and academic nanopore DNA/RNA sequencing techniques, and explore how these advances can fuel developments in future protein sequencing with nanopores.

Genome-based prediction of cross-protective, HLA-DR-presented epitopes as putative vaccine antigens for multiple Bordetella species

IMPORTANCE

Pertussis, caused by Bordetella pertussis, can cause debilitating respiratory symptoms, so whole-cell pertussis vaccines (wPVs) were introduced in the 1940s. However, reactogenicity of wPV necessitated the development of acellular pertussis vaccines (aPVs) that were introduced in the 1990s. Since then, until the COVID-19 pandemic began, reported pertussis incidence was increasing, suggesting that aPVs do not induce long-lasting immunity and may not effectively prevent transmission. Additionally, aPVs do not provide protection against other Bordetella species that are observed during outbreaks. The significance of this work is in determining potential new vaccine antigens for multiple Bordetella species that are predicted to elicit long-term immune responses. Genome-based approaches have aided the development of novel vaccines; here, these methods identified Bordetella vaccine candidates that may be cross-protective and predicted to induce strong memory responses. These targets can lead to an improved vaccine with a strong safety profile while also strengthening the longevity of the immune response.

Intranasal challenge with B. pertussis leads to more severe disease manifestations in mice than aerosol challenge

The murine Bordetella pertussis challenge model has been utilized in preclinical research for decades. Currently, inconsistent methodologies are employed by researchers across the globe, making it difficult to compare findings. The objective of this work was to utilize the CD-1 mouse model with two routes of challenge, intranasal and aerosol administration of B. pertussis, to understand the differences in disease manifestation elicited via each route. We observed that both routes of B. pertussis challenge result in dose-dependent colonization of the respiratory tract, but overall, intranasal challenge led to higher bacterial burden in the nasal lavage, trachea, and lung. Furthermore, high dose intranasal challenge results in induction of leukocytosis and pro-inflammatory cytokine responses compared to aerosol challenge. These data highlight crucial differences in B. pertussis challenge routes that should be considered during experimental design.

Cross-reactivity between vaccine antigens from the chitin deacetylase protein family improves survival in a mouse model of cryptococcosis

Meningitis due to the fungal pathogen Cryptococcus neoformans is estimated to cause nearly 200,000 deaths annually, mostly in resource-limited regions. We previously identified cryptococcal protein antigens which, when delivered in glucan particles, afford vaccine-mediated protection against an otherwise lethal infection. Many of these proteins exhibit significant homology to other similar cryptococcal proteins leading us to hypothesize that protection may be augmented by immunologic cross-reactivity to multiple members of a protein family. To examine the significance of protein cross-reactivity in vaccination, we utilized strains of Cryptococcus that are genetically deficient in select antigens, yet are still lethal in mice. Vaccination with a protein without homologs (e.g., Mep1 and Lhc1) protected against challenge with wild-type Cryptococcus, but not against a deletion strain lacking that protein. Contrastingly, vaccination with a single chitin deacetylase (Cda) protein protected against the corresponding deletion strain, presumably due to host recognition of one or more other family members still expressed in this strain. Vaccination with a single Cda protein induced cross-reactive antibody and interferon-gamma (IFNγ) immune responses to other Cda protein family members. Paradoxically, we saw no evidence of cross-protection within the carboxypeptidase family of proteins. Factors such as in vivo protein expression and the degree of homology across the family could inform the extent to which vaccine-mediated immunity is amplified. Together, these data suggest a role for prioritizing protein families in fungal vaccine design: increasing the number of immune targets generated by a single antigen may improve efficacy while diminishing the risk of vaccine-resistant strains arising from gene mutations.

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.

The Fim and FhaB adhesins play a crucial role in nasal cavity infection and Bordetella pertussis transmission in a novel mouse catarrhal infection model

Pulmonary infections caused by Bordetella pertussis used to be the prime cause of infant mortality in the pre-vaccine era and mouse models of pertussis pneumonia served in characterization of B. pertussis virulence mechanisms. However, the biologically most relevant catarrhal disease stage and B. pertussis transmission has not been adequately reproduced in adult mice due to limited proliferation of the human-adapted pathogen on murine nasopharyngeal mucosa. We used immunodeficient C57BL/6J MyD88 KO mice to achieve B. pertussis proliferation to human-like high counts of 108 viable bacteria per nasal cavity to elicit rhinosinusitis accompanied by robust shedding and transmission of B. pertussis bacteria to adult co-housed MyD88 KO mice. Experiments with a comprehensive set of B. pertussis mutants revealed that pertussis toxin, adenylate cyclase toxin-hemolysin, the T3SS effector BteA/BopC and several other known virulence factors were dispensable for nasal cavity infection and B. pertussis transmission in the immunocompromised MyD88 KO mice. In contrast, mutants lacking the filamentous hemagglutinin (FhaB) or fimbriae (Fim) adhesins infected the nasal cavity poorly, shed at low levels and failed to productively infect co-housed MyD88 KO or C57BL/6J mice. FhaB and fimbriae thus appear to play a critical role in B. pertussis transmission. The here-described novel murine model of B. pertussis-induced nasal catarrh opens the way to genetic dissection of host mechanisms involved in B. pertussis shedding and to validation of key bacterial transmission factors that ought to be targeted by future pertussis vaccines.

Towards comprehensive understanding of bacterial genetic diversity: large-scale amplifications in Bordetella pertussis and Mycobacterium tuberculosis

Bacterial genetic diversity is often described solely using base-pair changes despite a wide variety of other mutation types likely being major contributors. Tandem duplication/amplifications are thought to be widespread among bacteria but due to their often-intractable size and instability, comprehensive studies of these mutations are rare. We define a methodology to investigate amplifications in bacterial genomes based on read depth of genome sequence data as a proxy for copy number. We demonstrate the approach with Bordetella pertussis, whose insertion sequence element-rich genome provides extensive scope for amplifications to occur. Analysis of data for 2430 B. pertussis isolates identified 272 putative amplifications, of which 94 % were located at 11 hotspot loci. We demonstrate limited phylogenetic connection for the occurrence of amplifications, suggesting unstable and sporadic characteristics. Genome instability was further described in vitro using long-read sequencing via the Nanopore platform, which revealed that clonally derived laboratory cultures produced heterogenous populations rapidly. We extended this research to analyse a population of 1000 isolates of another important pathogen, Mycobacterium tuberculosis. We found 590 amplifications in M. tuberculosis, and like B. pertussis, these occurred primarily at hotspots. Genes amplified in B. pertussis include those involved in motility and respiration, whilst in M. tuberuclosis, functions included intracellular growth and regulation of virulence. Using publicly available short-read data we predicted previously unrecognized, large amplifications in B. pertussis and M. tuberculosis. This reveals the unrecognized and dynamic genetic diversity of B. pertussis and M. tuberculosis, highlighting the need for a more holistic understanding of bacterial genetics.

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.

Fha Deficient Bordetella pertussis Isolates in Iran with 50 Years Whole Cell Pertussis Vaccination

Background:

Bordetella pertussis, a highly contagious respiratory. Notably, the resurgence of pertussis has recently been associated with the lacking production of vaccine virulence factors. This study aimed to screen pertactin (Prn) and filamentous hemagglutinin (Fha) production in Iran with 50 years’ whole cell vaccine (WCV) immunization program.

Methods:

Overall, 130 B. pertussis isolates collected from Pertussis Reference Laboratory of Iran during 2005–2018. Real-time PCR was performed by targeting IS481, ptxP, IS1001 and IS1002 for species confirmation of B. pertussis. Western-blot was used to evaluate the expression of virulence factors (pertactin and filamentous hemagglutinin).

Results:

All tested B. pertussis isolates expressed Prn and all except two isolates expressed Fha. We have sequenced genomes of these strains and identified differences compared with genome reference B. pertussis Tohama I.

Conclusion:

Many countries reporting Prn and Fha-deficiency due to acellular vaccine (ACV) pressure. Our results demonstrate in a country with WCV history, Fha-deficient isolates may rise independently. However, Prn-deficient isolates are more under the ACV pressure in B. pertussis isolates. Continues surveillance will provide a better understanding of the effect of WCV on the evolution of the pathogen deficiency.

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.

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.

Infant rhesus macaques as a non-human primate model of Bordetella pertussis infection

Background

The prevalent resurgence of pertussis has recently become a critical public health problem worldwide. To understand pertussis pathogenesis and the host response to both the pathogen and vaccines, a suitable pertussis animal model, particularly a non-human primate model, is necessary. Recently, a non-human primate pertussis model was successfully established with baboons. Rhesus macaques have been shown to be ideal animal models for several infectious diseases, but a model of infectious pertussis has not been established in these organisms. Studies on rhesus macaque models of pertussis were performed in the 1920s–1930s, but limited experimental details are available. Recent monkey pertussis models have not been successful because the typical clinical symptoms and transmission have not been achieved.

Methods

In the present study, infant rhesus macaques were challenged with Bordetella pertussis (B.p) using an aerosol method to evaluate the feasibility of this system as an animal model of pertussis.

Results

Upon aerosol infection, monkeys infected with the recently clinically isolated B.p strain 2016-CY-41 developed the typical whooping cough, leukocytosis, bacteria-positive nasopharyngeal wash (NPW), and interanimal transmission of pertussis. Both systemic and mucosal humoral responses were induced by B.p.

Conclusion

These results demonstrate that a model of pertussis was successfully established in infant rhesus macaques. This model provides a valuable platform for research on pertussis pathogenesis and evaluation of vaccine candidates.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-021-06090-y.

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.

Effect of FHA and Prn on Bordetella pertussis colonization of mice is dependent on vaccine type and anatomical site

Bordetella pertussis vaccine escape mutants that lack expression of the pertussis antigen pertactin (Prn) have emerged in vaccinated populations in the last 10–20 years. Additionally, clinical isolates lacking another acellular pertussis (aP) vaccine component, filamentous hemagglutinin (FHA), have been found sporadically. Here, we show that both whole-cell pertussis (wP) and aP vaccines induced protection in the lungs of mice, but that the wP vaccine was more effective in nasal clearance. Importantly, bacterial populations isolated from the lungs shifted to an FHA-negative phenotype due to frameshift mutations in the fhaB gene. Loss of FHA expression was strongly selected for in Prn-deficient strains in the lungs following aP but not wP vaccination. The combined loss of Prn and FHA led to complete abrogation of bacterial surface binding by aP-induced serum antibodies. This study demonstrates vaccine- and anatomical site-dependent adaptation of B. pertussis and has major implications for the design of improved pertussis vaccines.

Optimization of sample preparation for culture-independent sequencing of Bordetella pertussis

Bordetella pertussis, the aetiological agent of whooping cough, is re-emerging globally despite widespread vaccination. B. pertussis is highly infectious and, prior to vaccination programmes, was the leading cause of infant mortality. The WHO estimated that over 600 000 deaths are prevented annually by pertussis vaccination, but B. pertussis infection was still responsible for over 63 000 deaths globally in 2013. The re-emergence of B. pertussis has been linked to strains with inactive or absent major virulence factors included in vaccines such as pertactin, pertussis toxin and filamentous haemagglutinin. Thus, the molecular surveillance of currently circulating strains is critical in understanding and controlling B. pertussis. Such information provides data on strains to inform control measures and the identification of future vaccine antigens. Current surveillance and typing methods for B. pertussis rely on the availability of clinical isolates. However, since the 1990s, the majority of pertussis cases have been diagnosed by PCR, where an isolate is not needed. The rapid decline in the availability of B. pertussis isolates impacts our ability to monitor this infection. The growing uptake of next-generation sequencing (NGS) has offered the opportunity for culture-independent genome sequencing and typing of this fastidious pathogen. Therefore, the objective of the study was to optimize respiratory sample preparation, independent of culture, in order to type B. pertussis using NGS. The study compared commercial depletion kits and specimen-processing methods using selective lysis detergents. The goal was to deplete human DNA, the major obstacle for sequencing a pathogen directly from a clinical sample. Samples spiked with a clinically relevant amount of B. pertussis were used to provide comparison between the different methods. Commercial depletion kits including the MolYsis, Qiagen Microbiome and NEBNext Kits were tested. Previously published methods, for Saponin and TritonX-100, were also trialled as a depletion. The ratio of B. pertussis to human DNA was determined by real-time PCR for ERV3 and IS481 (as markers of human and B. pertussis DNA, respectively), then samples were sequenced using the Illumina NextSeq 500 platform. The number of human and B. pertussis sequenced reads were then compared between treatments. The results showed that commercial kits reduced the human DNA present, but also reduced the concentration of target B. pertussis. However, selective lysis with Saponin treatment resulted in almost undetectable levels of human DNA, with minimal loss of target B. pertussis DNA. Sequencing read depth improved five-fold in reads to B. pertussis. Our investigation delivered a potential protocol that will enable the public health laboratory surveillance of B. pertussis in the era of culture-independent testing.

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.

Resolving the complex Bordetella pertussis genome using barcoded nanopore sequencing

The genome of Bordetella pertussis is complex, with high G+C content and many repeats, each longer than 1000 bp. Long-read sequencing offers the opportunity to produce single-contig B. pertussis assemblies using sequencing reads which are longer than the repetitive sections, with the potential to reveal genomic features which were previously unobservable in multi-contig assemblies produced by short-read sequencing alone. We used an R9.4 MinION flow cell and barcoding to sequence five B. pertussis strains in a single sequencing run. We then trialled combinations of the many nanopore user community-built long-read analysis tools to establish the current optimal assembly pipeline for B. pertussis genome sequences. This pipeline produced closed genome sequences for four strains, allowing visualization of inter-strain genomic rearrangement. Read mapping to the Tohama I reference genome suggests that the remaining strain contains an ultra-long duplicated region (almost 200 kbp), which was not resolved by our pipeline; further investigation also revealed that a second strain that was seemingly resolved by our pipeline may contain an even longer duplication, albeit in a small subset of cells. We have therefore demonstrated the ability to resolve the structure of several B. pertussis strains per single barcoded nanopore flow cell, but the genomes with highest complexity (e.g. very large duplicated regions) remain only partially resolved using the standard library preparation and will require an alternative library preparation method. For full strain characterization, we recommend hybrid assembly of long and short reads together; for comparison of genome arrangement, assembly using long reads alone is sufficient.

Comparative genomics of Czech vaccine strains of Bordetella pertussis

Bordetella pertussis is a strictly human pathogen causing the respiratory infectious disease called whooping cough or pertussis. B. pertussis adaptation to acellular pertussis vaccine pressure has been repeatedly highlighted, but recent data indicate that adaptation of circulating strains started already in the era of the whole cell pertussis vaccine (wP) use. We sequenced the genomes of five B. pertussis wP vaccine strains isolated in the former Czechoslovakia in the pre-wP (1954-1957) and early wP (1958-1965) eras, when only limited population travel into and out of the country was possible. Four isolates exhibit a similar genome organization and form a distinct phylogenetic cluster with a geographic signature. The fifth strain is rather distinct, both in genome organization and SNP-based phylogeny. Surprisingly, despite isolation of this strain before 1966, its closest sequenced relative appears to be a recent isolate from the US. On the genome content level, the five vaccine strains contained both new and already described regions of difference. One of the new regions contains duplicated genes potentially associated with transport across the membrane. The prevalence of this region in recent isolates indicates that its spread might be associated with selective advantage leading to increased strain fitness.