Assessing the Impact of the 2020 Council of State and Territorial Epidemiologists Case Definition for Pertussis on Reported Pertussis Cases.

BACKGROUND: In 2020, the Council of State and Territorial Epidemiologists (CSTE) pertussis case definition was modified; the main change was classifying polymerase chain reaction (PCR)-positive cases as confirmed, regardless of cough duration. Pertussis data reported through Enhanced Pertussis Surveillance (EPS) in 7 sites and the National Notifiable Diseases Surveillance System (NNDSS) were used to evaluate the impact of the new case definition. METHODS: We compared the number of EPS cases with cough onset in 2020 to the number that would have been reported based on the prior (2014) CSTE case definition. To assess the impact of the change nationally, the proportion of EPS cases newly reportable under the 2020 CSTE case definition was applied to 2020 NNDSS data to estimate how many additional cases were captured nationally. RESULTS: Among 442 confirmed and probable cases reported to EPS states in 2020, 42 (9.5%) were newly reportable according to the 2020 case definition. Applying this proportion to the 6124 confirmed and probable cases reported nationally in 2020, we estimated that the new definition added 582 cases. Had the case definition not changed, reported cases in 2020 would have decreased by 70% from 2019; the observed decrease was 67%. CONCLUSIONS: Despite a substantial decrease in reported pertussis cases in the setting of coronavirus disease 2019 (COVID-19), our data show that the 2020 pertussis case definition change resulted in additional case reporting compared with the previous case definition, providing greater opportunities for public health interventions such as prophylaxis of close contacts.

Evaluation of Asymptomatic Bordetella Carriage in a Convenience Sample of Children and Adolescents in Atlanta, Georgia, United States.

Few data exist on asymptomatic carriage of Bordetella species among populations receiving acellular pertussis vaccine. We conducted a cross-sectional study among acellular-vaccinated children presenting to an emergency department (ED). Bordetella pertussis carriage prevalence was <1% in this population, a lower prevalence than that found in recent studies among whole-cell pertussis-vaccinated participants.

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.

Genomic characterization of cocirculating Corynebacterium diphtheriae and non-diphtheritic Corynebacterium species among forcibly displaced Myanmar nationals, 2017-2019.

Respiratory diphtheria is a serious infection caused by toxigenic Corynebacterium diphtheriae, and disease transmission mainly occurs through respiratory droplets. Between 2017 and 2019, a large diphtheria outbreak among forcibly displaced Myanmar nationals densely settled in Bangladesh was investigated. Here we utilized whole-genome sequencing (WGS) to characterize recovered isolates of C. diphtheriae and two co-circulating non-diphtheritic Corynebacterium (NDC) species - C. pseudodiphtheriticum and C. propinquum. C. diphtheriae isolates recovered from all 53 positive cases in this study were identified as toxigenic biovar mitis, exhibiting intermediate resistance to penicillin, and formed four phylogenetic clusters circulating among multiple refugee camps. Additional sequenced isolates collected from two patients showed co-colonization with non-toxigenic C. diphtheriae biovar gravis, one of which exhibited decreased susceptibility to the first-line antibiotics and harboured a novel 23-kb multidrug resistance plasmid. Results of phylogenetic reconstruction and virulence-related gene contents of the recovered NDC isolates indicated they were likely commensal organisms, though 80.4 %(45/56) were not susceptible to erythromycin, and most showed high minimum inhibition concentrations against azithromycin. These results demonstrate the high resolution with which WGS can aid molecular investigation of diphtheria outbreaks, through the quantification of bacterial genetic relatedness, as well as the detection of virulence factors and antibiotic resistance markers among case isolates.

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.

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.

Detection and Characterization of Diphtheria Toxin Gene-Bearing Corynebacterium Species through a New Real-Time PCR Assay.

Respiratory diphtheria, characterized by a firmly adherent pseudomembrane, is caused by toxin-producing strains of Corynebacterium diphtheriae, with similar illness produced occasionally by toxigenic Corynebacterium ulcerans or, rarely, Corynebacterium pseudotuberculosis While diphtheria laboratory confirmation requires culture methods to determine toxigenicity, real-time PCR (RT-PCR) provides a faster method to detect the toxin gene (tox). Nontoxigenic tox-bearing (NTTB) Corynebacterium isolates have been described, but impact of these isolates on the accuracy of molecular diagnostics is not well characterized. Here, we describe a new triplex RT-PCR assay to detect tox and distinguish C. diphtheriae from the closely related species C. ulcerans and C. pseudotuberculosis Analytical sensitivity and specificity of the assay were assessed in comparison to culture using 690 previously characterized microbial isolates. The new triplex assay characterized Corynebacterium isolates accurately, with 100% analytical sensitivity for all targets. Analytical specificity with isolates was 94.1%, 100%, and 99.5% for tox, Diph_rpoB, and CUP_rpoB targets, respectively. Twenty-nine NTTB Corynebacterium isolates, representing 5.9% of 494 nontoxigenic isolates tested, were detected by RT-PCR. Whole-genome sequencing of NTTB isolates revealed varied mutations putatively underlying their lack of toxin production, as well as eight isolates with no mutation in tox or the promoter region. This new Corynebacterium RT-PCR method provides a rapid tool to screen isolates and identify probable diphtheria cases directly from specimens. However, the sporadic occurrence of NTTB isolates reinforces the viewpoint that diphtheria culture diagnostics continue to provide the most accurate case confirmation.

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.

Clinical evaluation and validation of laboratory methods for the diagnosis of Bordetella pertussis infection: Culture, polymerase chain reaction (PCR) and anti-pertussis toxin IgG serology (IgG-PT).

INTRODUCTION: The appropriate use of clinically accurate diagnostic tests is essential for the detection of pertussis, a poorly controlled vaccine-preventable disease. The purpose of this study was to estimate the sensitivity and specificity of different diagnostic criteria including culture, multi-target polymerase chain reaction (PCR), anti-pertussis toxin IgG (IgG-PT) serology, and the use of a clinical case definition. An additional objective was to describe the optimal timing of specimen collection for the various tests. METHODS: Clinical specimens were collected from patients with cough illness at seven locations across the United States between 2007 and 2011. Nasopharyngeal and blood specimens were collected from each patient during the enrollment visit. Patients who had been coughing for ≤ 2 weeks were asked to return in 2-4 weeks for collection of a second, convalescent blood specimen. Sensitivity and specificity of each diagnostic test were estimated using three methods-pertussis culture as the "gold standard," composite reference standard analysis (CRS), and latent class analysis (LCA). RESULTS: Overall, 868 patients were enrolled and 13.6% were B. pertussis positive by at least one diagnostic test. In a sample of 545 participants with non-missing data on all four diagnostic criteria, culture was 64.0% sensitive, PCR was 90.6% sensitive, and both were 100% specific by LCA. CRS and LCA methods increased the sensitivity estimates for convalescent serology and the clinical case definition over the culture-based estimates. Culture and PCR were most sensitive when performed during the first two weeks of cough; serology was optimally sensitive after the second week of cough. CONCLUSIONS: Timing of specimen collection in relation to onset of illness should be considered when ordering diagnostic tests for pertussis. Consideration should be given to including IgG-PT serology as a confirmatory test in the Council of State and Territorial Epidemiologists (CSTE) case definition for pertussis.

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.

Development of a qualitative assay for screening of Bordetella pertussis isolates for pertussis toxin production.

Bordetella pertussis infection has been increasing in the US, with reported cases reaching over 50,000 in 2012, a number last observed in the 1950s. Concurrently, B. pertussis lacking the pertactin protein, one of the immunogens included in the acellular vaccine formulations, has rapidly emerged since 2010, and has become the predominant circulating phenotype. Monitoring the production of the remaining acellular vaccine immunogens, such as pertussis toxin (Pt), is a critical next step. To date, methods for screening Pt have been either through genomic sequencing means or by conventional ELISAs. However, sequencing limits detection to the DNA level, missing potential disruptions in transcription or translation. Conventional ELISAs are beneficial for detecting the protein; however, they can often suffer from poor sensitivity and specificity. Here we describe a rapid, highly sensitive and specific electrochemiluminescent capture ELISA that can detect Pt production in prepared inactivated bacterial suspensions. Over 340 isolates were analyzed and analytical validation parameters, such as precision, reproducibility, and stability, were rigorously tested. Intra-plate and inter-plate variability measured at 9.8% and 11.5%, respectively. Refrigerated samples remained stable for two months and variability was unaffected (coefficient of variation was 12%). Interestingly, despite the intention of being a qualitative method, the assay was sensitive enough to detect a small, but statistically significant, difference in protein production between different pertussis promoter allelic groups of strains, ptxP1 and ptxP3. This technology has the ability to perform screening of multiple antigens at one time, thus, improving testing characteristics while minimizing costs, specimen volume, and testing time.

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.

Evaluation of Commercial Assays for Single-Point Diagnosis of Pertussis in the US.

BACKGROUND: Pertussis serodiagnosis is increasingly being used in the United States despite the lack of a US Food and Drug Administration-approved, commercially available assay. To better understand the utility of these assays in diagnosing pertussis, serology assays were evaluated for analytical parameters and clinical accuracy. METHODS: Forty-three antigen-antibody combinations were evaluated for single-point diagnosis of pertussis. Serum panels included sera from laboratory-confirmed cases, an international reference standard, and healthy donors. Phase I panel (n = 20) of sera was used to assess precision, linearity, and accuracy; Phase II panel (n = 226) followed with positive percent agreement (PPA) and negative percent agreement (NPA) estimates. Analytical analyses included coefficients of variation (CV) and concordance correlation coefficients (rc). RESULTS: Intra-analyst variability was found to be relatively low among samples per assay, with only 6% (78 of 1240) having CV >20%, primarily with the highly concentrated immunoglobulin (Ig)G anti-pertussis toxin (PT) specimens and IgM assays. The rc measurements to assess linearity ranged between 0.282 and 0.994, 0.332 and 0.999, and -0.056 and 0.482 for IgA, IgG, and IgM, respectively. Analytical accuracy for calibrated IgG anti-PT assays was 86%-115%. The PPA and NPA varied greatly for all assays; PPA/NPA ranges for IgA, IgG, and IgM assays, with culture and/or polymerase chain reaction positivity as control, were 29-90/13-100, 26-96/27-100, and 0-73/42-100, respectively. In IgG assays, mixing filamentous hemagglutinin antigen with PT increased PPA but decreased NPA. CONCLUSIONS: Seroassays varied substantially under both analytical and clinical parameters; however, those that were calibrated to a reference standard were highly accurate. Our findings support incorporation of calibrated pertussis seroassays to the pertussis case definition for improved diagnosis and surveillance.

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

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

Complete Genome Sequences of 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.

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

During 2010 and 2012, California and Vermont, respectively, experienced statewide epidemics of pertussis with differences seen in the demographic affected, case clinical presentation, and molecular epidemiology of the circulating strains. To overcome limitations of the current molecular typing methods for pertussis, we utilized whole-genome sequencing to gain a broader understanding of how current circulating strains are causing large epidemics. Through the use of combined next-generation sequencing technologies, this study compared de novo, single-contig genome assemblies from 31 out of 33 Bordetella pertussis isolates collected during two separate pertussis statewide epidemics and 2 resequenced vaccine strains. Final genome architecture assemblies were verified with whole-genome optical mapping. Sixteen distinct genome rearrangement profiles were observed in epidemic isolate genomes, all of which were distinct from the genome structures of the two resequenced vaccine strains. These rearrangements appear to be mediated by repetitive sequence elements, such as high-copy-number mobile genetic elements and rRNA operons. Additionally, novel and previously identified single nucleotide polymorphisms were detected in 10 virulence-related genes in the epidemic isolates. Whole-genome variation analysis identified state-specific variants, and coding regions bearing nonsynonymous mutations were classified into functional annotated orthologous groups. Comprehensive studies on whole genomes are needed to understand the resurgence of pertussis and develop novel tools to better characterize the molecular epidemiology of evolving B. pertussis populations. IMPORTANCE Pertussis, or whooping cough, is the most poorly controlled vaccine-preventable bacterial disease in the United States, which has experienced a resurgence for more than a decade. Once viewed as a monomorphic pathogen, B. pertussis strains circulating during epidemics exhibit diversity visible on a genome structural level, previously undetectable by traditional sequence analysis using short-read technologies. For the first time, we combine short- and long-read sequencing platforms with restriction optical mapping for single-contig, de novo assembly of 31 isolates to investigate two geographically and temporally independent U.S. pertussis epidemics. These complete genomes reshape our understanding of B. pertussis evolution and strengthen molecular epidemiology toward one day understanding the resurgence of pertussis.

Changes in Predominance of Pulsed-Field Gel Electrophoresis Profiles of Bordetella pertussis Isolates, United States, 2000-2012.

To clarify the characteristics of circulating Bordetella pertussis isolates, we used pulsed-field gel electrophoresis (PFGE) to analyze 5,262 isolates collected in the United States during 2000-2012. We found 199 PFGE profiles; 5 profiles accounted for 72% of isolates. The most common profile, CDC013, accounted for 35%-46% of isolates tested from 2000-2009; however, the proportion of isolates of this profile rapidly decreased in 2010. Profile CDC237, first seen in 2009, increased rapidly and accounted for 29% of 2012 isolates. No location bias was observed among profiles during 2000-2010, but differences were observed among isolates from different states during 2012. Predominant profiles match those observed in recent European PFGE studies. PFGE profile changes are concurrent with other recent molecular changes in B. pertussis and may be contributing to the reemergence of pertussis in the United States. Continued PFGE monitoring is critical for understanding the changing epidemiology of pertussis.

Bordetella pertussis Strain Lacking Pertactin and Pertussis Toxin.

A Bordetella pertussis strain lacking 2 acellular vaccine immunogens, pertussis toxin and pertactin, was isolated from an unvaccinated infant in New York State in 2013. Comparison with a French strain that was pertussis toxin-deficient, pertactin wild-type showed that the strains carry the same 28-kb deletion in similar genomes.