Genome Structural Diversity among 31 Bordetella pertussis Isolates from Two Recent U.S. Whooping Cough Statewide Epidemics

Strengthening Bordetella pertussis genomic surveillance by direct sequencing of residual positive specimens

Whole-genome sequencing (WGS) of microbial pathogens recovered from patients with infectious disease facilitates high-resolution strain characterization and molecular epidemiology. However, increasing reliance on culture-independent methods to diagnose infectious diseases has resulted in few isolates available for WGS. Here, we report a novel culture-independent approach to genome characterization of Bordetella pertussis, the causative agent of pertussis and a paradigm for insufficient genomic surveillance due to limited culture of clinical isolates. Sequencing libraries constructed directly from residual pertussis-positive diagnostic nasopharyngeal specimens were hybridized with biotinylated RNA "baits" targeting B. pertussis fragments within complex mixtures that contained high concentrations of host and microbial background DNA. Recovery of B. pertussis genome sequence data was evaluated with mock and pooled negative clinical specimens spiked with reducing concentrations of either purified DNA or inactivated cells. Targeted enrichment increased the yield of B. pertussis sequencing reads up to 90% while simultaneously decreasing host reads to less than 10%. Filtered sequencing reads provided sufficient genome coverage to perform characterization via whole-genome single nucleotide polymorphisms and whole-genome multilocus sequencing typing. Moreover, these data were concordant with sequenced isolates recovered from the same specimens such that phylogenetic reconstructions from either consistently clustered the same putatively linked cases. The optimized protocol is suitable for nasopharyngeal specimens with diagnostic IS481 Ct < 35 and >10 ng DNA. Routine implementation of these methods could strengthen surveillance and study of pertussis resurgence by capturing additional cases with genomic characterization.

Whole-genome comparison of two same-genotype macrolide-resistant Bordetella pertussis isolates collected in Japan

The emergence of macrolide-resistant Bordetella pertussis (MRBP) is a significant problem because it reduces treatment options for pertussis and exacerbates the severity and spread of the disease. MRBP has been widely prevalent in mainland China since the 2010s and has been sporadically detected in other Asian countries. In Japan, two MRBP clinical strains were first isolated in Tokyo and Osaka between June and July 2018. The isolates BP616 in Osaka and BP625 in Tokyo harbored the same virulence-associated allelic genes (including ptxP1, ptxA1, prn1, fim3A, and fhaB3) and MT195 genotype and exhibited similar antimicrobial susceptibility profiles. However, despite their simultaneous occurrence, a distinguishable epidemiological link between these isolates could not be established. To gain further insight into the genetic relationship between these isolates in this study, we performed whole-genome analyses. Phylogenetic analysis based on genome-wide single-nucleotide polymorphisms revealed that the isolates belonged to one of the three clades of Chinese MRBP isolates, but there were 11 single-nucleotide polymorphism differences between BP616 and BP625. Genome structure analysis revealed two large inversions (202 and 523 kbp) and one small transposition (3.8 kbp) between the genomes. These findings indicate that the two Japanese MRBP isolates are closely related to Chinese MRBP isolates but are genomically distinct, suggesting that they were introduced into Japan from mainland China through different transmission routes.

Genomic characterization of Bordetella pertussis in South Africa, 2015-2019

Pertussis remains a public health concern in South Africa, with an increase in reported cases and outbreaks in recent years. Whole genome sequencing was performed on 32 Bordetella pertussis isolates sourced from three different surveillance programmes in South Africa between 2015 and 2019. Genome sequences were characterized using multilocus sequence typing, vaccine antigen genes (ptxP, ptxA, ptxB, prn and fimH) and overall genome structure. All isolates were sequence type 2 and harboured the pertussis toxin promoter allele ptxP3. The dominant genotype was ptxP3-ptxA1-ptxB2-prn2-fimH2 (31/32, 96.9 %), with no pertactin-deficient or other mutations in vaccine antigen genes identified. Amongst 21 isolates yielding closed genome assemblies, eight distinct genome structures were detected, with 61.9 % (13/21) of the isolates exhibiting three predominant structures. Increases in case numbers are probably not due to evolutionary changes in the genome but possibly due to other factors such as the cyclical nature of B. pertussis disease, waning immunity due to the use of acellular vaccines and/or population immunity gaps.

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.

Bps polysaccharide of Bordetella pertussis resists antimicrobial peptides by functioning as a dual surface shield and decoy and converts Escherichia coli into a respiratory pathogen

Infections and disease caused by the obligate human pathogen Bordetella pertussis (Bp) are increasing, despite widespread vaccinations. The current acellular pertussis vaccines remain ineffective against nasopharyngeal colonization, carriage, and transmission. In this work, we tested the hypothesis that Bordetella polysaccharide (Bps), a member of the poly-β-1,6-N-acetyl-D-glucosamine (PNAG/PGA) family of polysaccharides promotes respiratory tract colonization of Bp by resisting killing by antimicrobial peptides (AMPs). Genetic deletion of the bpsA-D locus, as well as treatment with the specific glycoside hydrolase Dispersin B, increased susceptibility to AMP-mediated killing. Bps was found to be both cell surface-associated and released during laboratory growth and mouse infections. Addition of bacterial supernatants containing Bps and purified Bps increased B. pertussis resistance to AMPs. By utilizing ELISA, immunoblot and flow cytometry assays, we show that Bps functions as a dual surface shield and decoy. Co-inoculation of C57BL/6J mice with a Bps-proficient strain enhanced respiratory tract survival of the Bps-deficient strain. In combination, the presented results highlight the critical role of Bps as a central driver of B. pertussis pathogenesis. Heterologous production of Bps in a non-pathogenic E. coli K12 strain increased AMP resistance in vitro, and augmented bacterial survival and pathology in the mouse respiratory tract. These studies can serve as a foundation for other PNAG/PGA polysaccharides and for the development of an effective Bp vaccine that includes Bps.

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.

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.

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.

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

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

Comparative Omics Analysis of Historic and Recent Isolates of Bordetella pertussis and Effects of Genome Rearrangements on Evolution

Despite high vaccination coverage, pertussis is increasing in many industrialized countries, including the Czech Republic. To better understand Bordetella pertussis resurgence, we analyzed historic strains and recent clinical isolates by using a comparative omics approach. Whole-genome sequencing showed that historic and recent isolates of B. pertussis have substantial variation in genome organization and form separate phylogenetic clusters. Subsequent RNA sequence analysis and liquid chromatography with mass tandem spectrometry analyses showed that these variations translated into discretely separated transcriptomic and proteomic profiles. When compared with historic strains, recent isolates showed increased expression of flagellar genes and genes involved in lipopolysaccharide biosynthesis and decreased expression of polysaccharide capsule genes. Compared with reference strain Tohama I, all strains had increased expression and production of the type III secretion system apparatus. We detected the potential link between observed effects and insertion sequence element–induced changes in gene context only for a few genes.

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.

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.

Chemical Synthesis and Immunological Evaluation of a Pentasaccharide Bearing Multiple Rare Sugars as a Potential Anti-pertussis Vaccine

With the infection rate of Bordetella pertussis at a 60-year high, there is an urgent need for new anti-pertussis vaccines. The lipopolysaccharide (LPS) of B. pertussis is an attractive antigen for vaccine development. With the presence of multiple rare sugars and unusual glycosyl linkages, the B. pertussis LPS is a highly challenging synthetic target. In this work, aided by molecular dynamics simulation and modeling, a pertussis-LPS-like pentasaccharide was chemically synthesized for the first time. The pentasaccharide was conjugated with a powerful carrier, bacteriophage Qβ, as a vaccine candidate. Immunization of mice with the conjugate induced robust anti-glycan IgG responses with IgG titers reaching several million enzyme-linked immunosorbent assay (ELISA) units. The antibodies generated were long lasting and boostable and could recognize multiple clinical strains of B. pertussis, highlighting the potential of Qβ-glycan as a new anti-pertussis vaccine.

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.

Genomic Survey of Bordetella pertussis Diversity, United States, 2000-2013

We characterized 170 complete genome assemblies from clinical Bordetella pertussis isolates representing geographic and temporal diversity in the United States. These data capture genotypic shifts, including increased pertactin deficiency, occurring amid the current pertussis disease resurgence and provide a foundation for needed research to direct future public health control strategies.

Association of CMV genomic mutations with symptomatic infection and hearing loss in congenital CMV infection

Background

Congenital cytomegalovirus (cCMV) infection is the most common congenital infection and a leading cause of long-term neurological and sensory sequelae, the most common being sensorineural hearing loss (SNHL). Despite extensive research, clinical or laboratory markers to identify CMV infected children with increased risk for disease have not been identified. This study utilizes viral whole-genome next generation-sequencing (NGS) of specimens from congenitally infected infants to explore viral diversity and specific viral variants that may be associated with symptomatic infection and SNHL.

Methods

CMV DNA from urine specimens of 30 infants (17 asymptomatic, 13 symptomatic) was target enriched and next generation sequenced resulting in 93% coverage of the CMV genome allowing analysis of viral diversity.

Results

Variant frequency distribution was compared between children with symptomatic and asymptomatic cCMV and those with (n = 13) and without (n = 17) hearing loss. The CMV genes UL48A, UL88, US19 and US22 were found to have an increase in nucleotide diversity in symptomatic children; while UL57, UL20, UL104, US14, UL115, and UL35 had an increase in diversity in children with hearing loss. An analysis of single variant differences between symptomatic and asymptomatic children found UL55 to have the highest number, while the most variants associated with SNHL were in the RL11 gene family. In asymptomatic infants with SNHL, mutations were observed more frequently in UL33 and UL20.

Conclusion

CMV genomes from infected newborns can be mapped to 93% of the genome at a depth allowing accurate and reproducible analysis of polymorphisms for variant and gene discovery that may be linked to symptomatic and hearing loss outcomes.

Conserved Patterns of Symmetric Inversion in the Genome Evolution of Bordetella Respiratory Pathogens

Whooping cough (pertussis), primarily caused by Bordetella pertussis, has resurged in the United States, and circulating strains exhibit considerable chromosome structural fluidity in the form of rearrangement and deletion. The genus Bordetella includes additional pathogenic species infecting various animals, some even causing pertussis-like respiratory disease in humans; however, investigation of their genome evolution has been limited. We studied chromosome structure in complete genome sequences from 167 Bordetella species isolates, as well as 469 B. pertussis isolates, to gain a generalized understanding of rearrangement patterns among these related pathogens. Observed changes in gene order primarily resulted from large inversions and were only detected in species with genomes harboring multicopy insertion sequence (IS) elements, most notably B. holmesii and B. parapertussis While genomes of B. pertussis contain >240 copies of IS481, IS elements appear less numerous in other species and yield less chromosome structural diversity through rearrangement. These data were further used to predict all possible rearrangements between IS element copies present in Bordetella genomes, revealing that only a subset is observed among circulating strains. Therefore, while it appears that rearrangement occurs less frequently in other species than in B. pertussis, these clinically relevant respiratory pathogens likely experience similar mutation of gene order. The resulting chromosome structural fluidity presents both challenges and opportunity for the study of Bordetella respiratory pathogens.IMPORTANCE Bordetella pertussis is the primary agent of whooping cough (pertussis). The Bordetella genus includes additional pathogens of animals and humans, including some that cause pertussis-like respiratory illness. The chromosome of B. pertussis has previously been shown to exhibit considerable structural rearrangement, but insufficient data have prevented comparable investigation in related species. In this study, we analyze chromosome structure variation in several Bordetella species to gain a generalized understanding of rearrangement patterns in this genus. Just as in B. pertussis, we observed inversions in other species that likely result from common mutational processes. We used these data to further predict additional, unobserved inversions, suggesting that specific genome structures may be preferred in each species.

A Pertussis Outer Membrane Vesicle-Based Vaccine Induces Lung-Resident Memory CD4 T Cells and Protection Against Bordetella pertussis, Including Pertactin Deficient Strains

Pertussis is a respiratory infectious disease that has been resurged during the last decades. The change from the traditional multi-antigen whole-cell pertussis (wP) vaccines to acellular pertussis (aP) vaccines that consist of a few antigens formulated with alum, appears to be a key factor in the resurgence of pertussis in many countries. Though current aP vaccines have helped to reduce the morbidity and mortality associated with pertussis, they do not provide durable immunity or adequate protection against the disease caused by the current circulating strains of Bordetella pertussis, which have evolved in the face of the selection pressure induced by the vaccines. Based on the hypothesis that a new vaccine containing multiple antigens could overcome deficiencies in the current aP vaccines, we have designed and characterized a vaccine candidate based on outer membrane vesicle (OMVs). Here we show that the OMVs vaccine, but not an aP vaccine, protected mice against lung infection with a circulating pertactin (PRN)-deficient isolate. Using isogenic bacteria that in principle only differ in PRN expression, we found that deficiency in PRN appears to be largely responsible for the failure of the aP vaccine to protect against this circulating clinical isolates. Regarding the durability of induced immunity, we have already reported that the OMV vaccine is able to induce long-lasting immune responses that effectively prevent infection with B. pertussis. Consistent with this, here we found that CD4 T cells with a tissue-resident memory (T_RM) cell phenotype (CD44⁺CD62L^lowCD69⁺ and/or CD103⁺) accumulated in the lungs of mice 14 days after immunization with 2 doses of the OMVs vaccine. CD4 T_RM cells, which have previously been shown to play a critical role sustained protective immunity against B. pertussis, were also detected in mice immunized with wP vaccine, but not in the animals immunized with a commercial aP vaccine. The CD4 T_RM cells secreted IFN-γ and IL-17 and were significantly expanded through local proliferation following respiratory challenge of mice with B. pertussis. Our findings that the OMVs vaccine induce respiratory CD4 T_RM cells may explain the ability of this vaccine to induce long-term protection and is therefore an ideal candidate for a third generation vaccine against B. pertussis.

Genomic Sequencing of Bordetella pertussis for Epidemiology and Global Surveillance of Whooping Cough

Bordetella pertussis causes whooping cough, a highly contagious respiratory disease that is reemerging in many world regions. The spread of antigen-deficient strains may threaten acellular vaccine efficacy. Dynamics of strain transmission are poorly defined because of shortcomings in current strain genotyping methods. Our objective was to develop a whole-genome genotyping strategy with sufficient resolution for local epidemiologic questions and sufficient reproducibility to enable international comparisons of clinical isolates. We defined a core genome multilocus sequence typing scheme comprising 2,038 loci and demonstrated its congruence with whole-genome single-nucleotide polymorphism variation. Most cases of intrafamilial groups of isolates or of multiple isolates recovered from the same patient were distinguished from temporally and geographically cocirculating isolates. However, epidemiologically unrelated isolates were sometimes nearly undistinguishable. We set up a publicly accessible core genome multilocus sequence typing database to enable global comparisons of B. pertussis isolates, opening the way for internationally coordinated surveillance.

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.

New Pertussis Vaccines: A Need and a Challenge

Effective diphtheria, tetanus toxoids, whole-cell pertussis (wP) vaccines were used for massive immunization in the 1950s. The broad use of these vaccines significantly reduced the morbidity and mortality associated with pertussis. Because of reports on the induction of adverse reactions, less-reactogenic acellular vaccines (aP) were later developed and in many countries, especially the industrialized ones, the use of wP was changed to aP. For many years, the situation of pertussis seemed to be controlled with the use of these vaccines, however in the last decades the number of pertussis cases increased in several countries. The loss of the immunity conferred by the vaccines, which is faster in the individuals vaccinated with the acellular vaccines, and the evolution of the pathogen towards geno/phenotypes that escape more easily the immunity conferred by the vaccines were proposed as the main causes of the disease resurgence. According to their composition of few immunogens, the aP vaccines seem to be exerting a greater selection pressure on the circulating bacterial population causing the prevalence of bacterial isolates defective in the expression of vaccine antigens. Under this context, it is clear that new vaccines against pertussis should be developed. Several vaccine candidates are in preclinical development and few others have recently completed phaseI/phaseII trials. Vaccine candidate based on OMVs is a promising candidate since appeared overcoming the major weaknesses of current aP-vaccines. The most advanced development is the live attenuated-vaccine BPZE1 which has successfully completed a first-in-man clinical trial.

Pertactin-deficient Bordetella pertussis isolates in Poland-a country with whole-cell pertussis primary vaccination

The introduction of pertussis vaccination in the 1950s resulted in a significant decrease in the incidence of disease. However, since the 1990s many highly vaccinated countries have observed the re-emergence of the disease. One of the causes of this phenomenon might be related to the adaptation of Bordetella pertussis to vaccination. The purpose of the presented study was an investigation of the emergence and spread of vaccine antigen-deficient B. pertussis isolates in Poland and genomic characterization of the currently circulating pathogen population using PFGE, MLVA and MAST. The results revealed that all tested isolates expressed Ptx, FHA and ACT antigens but 15.4% (4/26) of isolates from 2010 to 2016 were Prn-deficient. Moreover, one TcfA-deficient isolate was collected in 2015. The genotyping showed a genetic distinction between the isolates circulating in 2010-2016 and isolates from previous periods. The majority of currently circulating isolates belonged to PFGE group IV (96.2%), type MT27 (73.1%), and carried ptxA1-ptxC2-ptxP3-prn2-tcfA2-fim2-1-fim3-1 alleles (61.5%). The unique genetic structure of the B. pertussis population in Poland has changed since 2010 and became similar to that observed in countries with aP vaccination. This could be a result of increasing use of aP vaccines (60% of primary vaccination in 2013) over wP vaccines, which have been broadly used for primary vaccination in Poland for decades.

Complete Genome Sequences of Bordetella pertussis Clinical Isolate FR5810 and Reference Strain Tohama from Combined Oxford Nanopore and Illumina Sequencing

Here, we describe the complete genome sequences of two Bordetella pertussis strains, FR5810, a clinical isolate recovered from the respiratory tract of an infant, and Tohama, a key reference strain for the species. Sequences were obtained using a hybrid approach combining Oxford Nanopore Technologies MinION and Illumina NextSeq 500 sequence data.

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.

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.

Surveillance of Circulating Bordetella pertussis Strains in Europe during 1998 to 2015

One reason for increased pertussis incidence is the adaptation of Bordetella pertussis to vaccine-induced immunity by modulating its genomic structure. This study, EUpert IV, includes 265 isolates collected from nine European countries during 2012 to 2015 (n = 265) and compares the results to previous EUpert I to III studies (1998 to 2009). The analyses included genotyping, serotyping, pulsed-field gel electrophoresis (PFGE), and multilocus variable-number tandem-repeat analysis (MLVA). Genotyping results showed only small variations among the common virulence genes of B. pertussis The frequencies of serotypes Fim2 and Fim3 varied among the four collections. Genomic analyses showed that MLVA type 27 increased to 80% between the periods of 1998 to 2001 and 2012 to 2015. Two PFGE profiles, BpSR3 (29.4%) and BpSR10 (27.2%), constituted more than 50% of the circulating isolates in the present collection. Our study indicates that the European B. pertussis population is changing and became more homogenous after the introduction of acellular pertussis vaccines.

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.

Complete Genome Sequences of Bordetella pertussis Isolates with Novel Pertactin-Deficient Deletions

Clinical isolates of the respiratory pathogen Bordetella pertussis in the United States have become predominantly deficient for the acellular vaccine immunogen pertactin through various independent mutations. Here, we report the complete genome sequences for four B. pertussis isolates that harbor novel deletions responsible for pertactin deficiency.

PacBio But Not Illumina Technology Can Achieve Fast, Accurate and Complete Closure of the High GC, Complex Burkholderia pseudomallei Two-Chromosome Genome

Although PacBio third-generation sequencers have improved the read lengths of genome sequencing which facilitates the assembly of complete genomes, no study has reported success in using PacBio data alone to completely sequence a two-chromosome bacterial genome from a single library in a single run. Previous studies using earlier versions of sequencing chemistries have at most been able to finish bacterial genomes containing only one chromosome with de novo assembly. In this study, we compared the robustness of PacBio RS II, using one SMRT cell and the latest P6-C4 chemistry, with Illumina HiSeq 1500 in sequencing the genome of Burkholderia pseudomallei, a bacterium which contains two large circular chromosomes, very high G+C content of 68–69%, highly repetitive regions and substantial genomic diversity, and represents one of the largest and most complex bacterial genomes sequenced, using a reference genome generated by hybrid assembly using PacBio and Illumina datasets with subsequent manual validation. Results showed that PacBio data with de novo assembly, but not Illumina, was able to completely sequence the B. pseudomallei genome without any gaps or mis-assemblies. The two large contigs of the PacBio assembly aligned unambiguously to the reference genome, sharing >99.9% nucleotide identities. Conversely, Illumina data assembled using three different assemblers resulted in fragmented assemblies (201–366 contigs), sharing only 92.2–100% and 92.0–100% nucleotide identities to chromosomes I and II reference sequences, respectively, with no indication that the B. pseudomallei genome consisted of two chromosomes with four copies of ribosomal operons. Among all assemblies, the PacBio assembly recovered the highest number of core and virulence proteins, and housekeeping genes based on whole-genome multilocus sequence typing (wgMLST). Most notably, assembly solely based on PacBio outperformed even hybrid assembly using both PacBio and Illumina datasets. Hybrid approach generated only 74 contigs, while the PacBio data alone with de novo assembly achieved complete closure of the two-chromosome B. pseudomallei genome without additional costly bench work and further sequencing. PacBio RS II using P6-C4 chemistry is highly robust and cost-effective and should be the platform of choice in sequencing bacterial genomes, particularly for those that are well-known to be difficult-to-sequence.

Complete Genome Sequence of Bordetella pertussis Pelita III, the Production Strain for an Indonesian Whole-Cell Pertussis Vaccine

PT Bio Farma, the sole World Health Organization-approved Indonesian vaccine producer, manufactures a whole-cell whooping cough vaccine (wP) that, as part of a pentavalent diphtheria-tetanus-pertussis/hepatitis B/Haemophilus influenzae b (DTP/HB/Hib) vaccine, is used in Indonesia and many other countries. We report here the whole-genome sequence for Bordetella pertussis Pelita III, PT Bio Farma’s wP production strain.

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.

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

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

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 Sequence of Bordetella pertussis Strain VA-190 Isolated from a Vaccinated 10-Year-Old Patient with Whooping Cough

The number of cases of pertussis has increased in the United States despite vaccination. We present the genome of an isolate of Bordetella pertussis from a vaccinated patient from Virginia. The genome was sequenced by long-read methodology and compared to that of a clinical isolate used for laboratory studies, D420.

Complete Genome Sequences of Four Different Bordetella sp. Isolates Causing Human Respiratory Infections

Species of the genus Bordetella associate with various animal hosts, frequently causing respiratory disease. Bordetella pertussis is the primary agent of whooping cough and other Bordetella species can cause similar cough illness. Here, we report four complete genome sequences from isolates of different Bordetella species recovered from human respiratory infections.

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.