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.

Complete Genome Sequences of Four Macrolide-Resistant Nondiphtheritic Corynebacterium Isolates.

This report describes the complete genome sequences of four isolates of the nondiphtheritic Corynebacterium (NDC) species Corynebacterium pseudodiphtheriticum and Corynebacterium propinquum, recovered during investigation of a large diphtheria outbreak in Bangladesh. These data will assist in better delineating the boundary between these related species and understanding their virulence potential.

Respiratory Illness Caused by Corynebacterium diphtheriae and C. ulcerans, and Use of Diphtheria Antitoxin in the United States, 1996-2018.

BACKGROUND: Respiratory diphtheria is a toxin-mediated disease caused by Corynebacterium diphtheriae. Diphtheria-like illness, clinically indistinguishable from diphtheria, is caused by Corynebacterium ulcerans, a zoonotic bacterium that can also produce diphtheria toxin. In the United States, respiratory diphtheria is nationally notifiable: specimens from suspected cases are submitted to the Centers for Disease Control and Prevention (CDC) for species and toxin confirmation, and diphtheria antitoxin (DAT) is obtained from CDC for treatment. We summarize the epidemiology of respiratory diphtheria and diphtheria-like illness and describe DAT use during 1996-2018 in the United States. METHODS: We described respiratory diphtheria cases reported to the National Notifiable Diseases Surveillance System (NNDSS) and C. ulcerans-related diphtheria-like illness identified through specimen submissions to CDC during 1996-2018. We reviewed DAT requests from 1997 to 2018. RESULTS: From 1996 to 2018, 14 respiratory diphtheria cases were reported to NNDSS. Among these 14 cases, 1 was toxigenic and 3 were nontoxigenic C. diphtheriae by culture and Elek, 6 were culture-negative but polymerase chain reaction (PCR)-positive for diphtheria toxin gene, 1 was culture-positive without further testing, and the remaining 3 were either not tested or tested negative. Five cases of respiratory diphtheria-like illness caused by toxigenic C. ulcerans were identified. DAT was requested by healthcare providers for 151 suspected diphtheria cases between 1997 and 2018, with an average of 11 requests per year from 1997 to 2007, and 3 per year from 2008 to 2018. CONCLUSIONS: Respiratory diphtheria remains rare in the United States, and requests for DAT have declined. Incidental identification of C. ulcerans-related diphtheria-like illness suggests surveillance of this condition might be warranted.

Investigation of a Large Diphtheria Outbreak and Cocirculation of Corynebacterium pseudodiphtheriticum Among Forcibly Displaced Myanmar Nationals, 2017-2019.

BACKGROUND: Diphtheria, a life-threatening respiratory disease, is caused mainly by toxin-producing strains of Corynebacterium diphtheriae, while nontoxigenic corynebacteria (eg, Corynebacterium pseudodiphtheriticum) rarely causes diphtheria-like illness. Recently, global diphtheria outbreaks have resulted from breakdown of health care infrastructures, particularly in countries experiencing political conflict. This report summarizes a laboratory and epidemiological investigation of a diphtheria outbreak among forcibly displaced Myanmar nationals in Bangladesh. METHODS: Specimens and clinical information were collected from patients presenting at diphtheria treatment centers. Swabs were tested for toxin gene (tox)-bearing C. diphtheriae by real-time polymerase chain reaction (RT-PCR) and culture. The isolation of another Corynebacterium species prompted further laboratory investigation. RESULTS: Among 382 patients, 153 (40%) tested tox positive for C. diphtheriae by RT-PCR; 31 (20%) PCR-positive swabs were culture confirmed. RT-PCR revealed 78% (298/382) of patients tested positive for C. pseudodiphtheriticum. Of patients positive for only C. diphtheriae, 63% (17/27) had severe disease compared to 55% (69/126) positive for both Corynebacterium species, and 38% (66/172) for only C. pseudodiphtheriticum. CONCLUSIONS: We report confirmation of a diphtheria outbreak and identification of a cocirculating Corynebacterium species. The high proportion of C. pseudodiphtheriticum codetection may explain why many suspected patients testing negative for C. diphtheriae presented with diphtheria-like symptoms.

Genomic epidemiology of nontoxigenic Corynebacterium diphtheriae from King County, Washington State, USA between July 2018 and May 2019.

Between July 2018 and May 2019, Corynebacterium diphtheriae was isolated from eight patients with non-respiratory infections, seven of whom experienced homelessness and had stayed at shelters in King County, WA, USA. All isolates were microbiologically identified as nontoxigenic C. diphtheriae biovar mitis. Whole-genome sequencing confirmed that all case isolates were genetically related, associated with sequence type 445 and differing by fewer than 24 single-nucleotide polymorphisms (SNPs). Compared to publicly available C. diphtheriae genomic data, these WA isolates formed a discrete cluster with SNP variation consistent with previously reported outbreaks. Virulence-related gene content variation within the highly related WA cluster isolates was also observed. These results indicated that genome characterization can readily support epidemiology of nontoxigenic C. diphtheriae.

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.

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.

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.

Risk Factors Associated With Bordetella pertussis Among Infants ≤4 Months of Age in the Pre-Tdap Era: United States, 2002-2005.

BACKGROUND: In the United States, infants have the highest reported pertussis incidence and death rates. Improved understanding of infant risk factors is needed to optimize prevention strategies. METHODS: We prospectively enrolled infants ≤4 months of age with incident-confirmed pertussis from 4 sites during 2002-2005 (preceding pertussis antigen-containing vaccination recommendations for adolescents/adults); each case-patient was age and site matched with 2 control subjects. Caregivers completed structured interviews. Infants and their contacts ≥11 years of age were offered serologic testing for IgG; being seropositive was defined as ≥94 antipertussis toxin IgG enzyme-linked immunosorbent assay units per milliliter. RESULTS: Enrolled subjects (115 case-patients; 230 control subjects) had 4396 contacts during incubation periods; 83 (72%) case-patients had ≥1 contact with prolonged (≥5 days) new cough in primary or secondary households. In multivariable analysis, the odds for pertussis were higher for infants with primary/secondary household contacts who had a prolonged new cough, compared with infants who did not. These contacts included mother [adjusted matched odds ratio (aMOR), 43.8; 95% confidence interval (CI), 6.45-298.0] and ≥1 nonmother contact (aMOR, 20.1; 95% CI, 6.48-62.7). Infants receiving breast milk with 0-1 formula feedings daily had decreased pertussis odds (aMOR, 0.27; 95% CI, 0.08-0.89), compared with those receiving more formula. Of 41 tested case-patients, 37 (90%) were seropositive. CONCLUSIONS: Pertussis in infants was associated with prolonged new cough (≥5 days) in infants' household contacts. Findings suggest that breastfeeding protects against pertussis and warrants recommendation with pertussis prevention strategies, which currently include pertussis vaccination of pregnant mothers and infants' close contacts.

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.

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.

Pertactin-negative Bordetella pertussis strains: evidence for a possible selective advantage.

BACKGROUND: A recent increase in Bordetella pertussis without the pertactin protein, an acellular vaccine immunogen, has been reported in the United States. Determining whether pertactin-deficient (PRN(-)) B. pertussis is evading vaccine-induced immunity or altering the severity of illness is needed. METHODS: We retrospectively assessed for associations between pertactin production and both clinical presentation and vaccine history. Cases with isolates collected between May 2011 and February 2013 from 8 states were included. We calculated unadjusted and adjusted odds ratios (ORs) using multivariable logistic regression analysis. RESULTS: Among 753 isolates, 640 (85%) were PRN(-). The age distribution differed between cases caused by PRN(-) B. pertussis and cases caused by B. pertussis producing pertactin (PRN(+)) (P = .01). The proportion reporting individual pertussis symptoms was similar between the 2 groups, except a higher proportion of PRN(+) case-patients reported apnea (P = .005). Twenty-two case-patients were hospitalized; 6% in the PRN(+) group compared to 3% in the PRN(-) group (P = .11). Case-patients having received at least 1 pertussis vaccine dose had a higher odds of having PRN(-) B. pertussis compared with unvaccinated case-patients (adjusted OR = 2.2; 95% confidence interval [CI], 1.3-4.0). When restricted to case-patients at least 1 year of age and those age-appropriately vaccinated, the adjusted OR increased to 2.7 (95% CI, 1.2-6.1). CONCLUSIONS: The significant association between vaccination and isolate pertactin production suggests that the likelihood of having reported disease caused by PRN(-) compared with PRN(+) strains is greater in vaccinated persons. Additional studies are needed to assess whether vaccine effectiveness is diminished against PRN(-) strains.

Molecular epidemiology of the pertussis epidemic in Washington State in 2012.

Although pertussis disease is vaccine preventable, Washington State experienced a substantial rise in pertussis incidence beginning in 2011. By June 2012, the reported cases reached 2,520 (37.5 cases per 100,000 residents), a 1,300% increase compared with the same period in 2011. We assessed the molecular epidemiology of this statewide epidemic using 240 isolates collected from case patients reported from 19 of 39 Washington counties during 2012 to 2013. The typing methods included pulsed-field gel electrophoresis (PFGE), multilocus variable number tandem repeat analysis (MLVA), multilocus sequence typing (MLST), and pertactin gene (prn) mutational analysis. Using the scheme PFGE-MLVA-MLST-prn mutations-Prn deficiency, the 240 isolates comprised 65 distinct typing profiles. Thirty-one PFGE types were found, with the most common types, CDC013 (n = 51), CDC237 (n = 44), and CDC002 (n = 42), accounting for 57% of them. Eleven MLVA types were observed, mainly comprising type 27 (n = 183, 76%). Seven MLST types were identified, with the majority of the isolates typing as prn2-ptxP3-ptxA1-fim3-1 (n = 157, 65%). Four different prn mutations accounted for the 76% of isolates exhibiting pertactin deficiency. PFGE provided the highest discriminatory power (D = 0.87) and was found to be a more powerful typing method than MLVA and MLST combined (D = 0.67). This study provides evidence for the continued predominance of MLVA 27 and prn2-ptxP3-ptxA1 alleles, along with the reemergence of the fim3-1 allele. Our results indicate that the Bordetella pertussis population causing this epidemic was diverse, with a few molecular types predominating. The PFGE, MLVA, and MLST profiles were consistent with the predominate types circulating in the United States and other countries. For prn, several mutations were present in multiple molecular types.

Genome Sequences of Nine Bordetella holmesii Strains Isolated in the United States.

An increasing number of pertussis-like cases are attributed to the emergent pathogen Bordetella holmesii. The genomes of 9 clinical isolates show that they are clonal, lack the virulence factors encoded by B. pertussis, and are more similar to nonpertussis bordetellae. New markers for B. holmesii can be developed using these sequences.

Global population structure and evolution of Bordetella pertussis and their relationship with vaccination.

Bordetella pertussis causes pertussis, a respiratory disease that is most severe for infants. Vaccination was introduced in the 1950s, and in recent years, a resurgence of disease was observed worldwide, with significant mortality in infants. Possible causes for this include the switch from whole-cell vaccines (WCVs) to less effective acellular vaccines (ACVs), waning immunity, and pathogen adaptation. Pathogen adaptation is suggested by antigenic divergence between vaccine strains and circulating strains and by the emergence of strains with increased pertussis toxin production. We applied comparative genomics to a worldwide collection of 343 B. pertussis strains isolated between 1920 and 2010. The global phylogeny showed two deep branches; the largest of these contained 98% of all strains, and its expansion correlated temporally with the first descriptions of pertussis outbreaks in Europe in the 16th century. We found little evidence of recent geographical clustering of the strains within this lineage, suggesting rapid strain flow between countries. We observed that changes in genes encoding proteins implicated in protective immunity that are included in ACVs occurred after the introduction of WCVs but before the switch to ACVs. Furthermore, our analyses consistently suggested that virulence-associated genes and genes coding for surface-exposed proteins were involved in adaptation. However, many of the putative adaptive loci identified have a physiological role, and further studies of these loci may reveal less obvious ways in which B. pertussis and the host interact. This work provides insight into ways in which pathogens may adapt to vaccination and suggests ways to improve pertussis vaccines. IMPORTANCE Whooping cough is mainly caused by Bordetella pertussis, and current vaccines are targeted against this organism. Recently, there have been increasing outbreaks of whooping cough, even where vaccine coverage is high. Analysis of the genomes of 343 B. pertussis isolates from around the world over the last 100 years suggests that the organism has emerged within the last 500 years, consistent with historical records. We show that global transmission of new strains is very rapid and that the worldwide population of B. pertussis is evolving in response to vaccine introduction, potentially enabling vaccine escape.