BMScan: using whole genome similarity to rapidly and accurately identify bacterial meningitis causing species

Abstract Background Bacterial meningitis is a life-threatening infection that remains a public health concern. Bacterial meningitis is commonly caused by the following species: Neisseria meningitidis, Streptococcus pneumoniae, Listeria monocytogenes, Haemophilus influenzae and Escherichia coli. Here, we describe BMScan (Bacterial Meningitis Scan), a whole-genome analysis tool for the species identification of bacterial meningitis-causing and closely-related pathogens, an essential step for case management and disease surveillance. BMScan relies on a reference collection that contains genomes for 17 focal species to scan against to identify a given species. We established this reference collection by supplementing publically available genomes from RefSeq with genomes from the isolate collections of the Centers for Disease Control Bacterial Meningitis Laboratory and the Minnesota Department of Health Public Health Laboratory, and then filtered them down to a representative set of genomes which capture the diversity for each species. Using this reference collection, we evaluated two genomic comparison algorithms, Mash and Average Nucleotide Identity, for their ability to accurately and rapidly identify our focal species. Results We found that the results of Mash were strongly correlated with the results of ANI for species identification, while providing a significant reduction in run-time. This drastic difference in run-time enabled the rapid scanning of large reference genome collections, which, when combined with species-specific threshold values, facilitated the development of BMScan. Using a validation set of 15,503 genomes of our species of interest, BMScan accurately identified 99.97% of the species within 16 min 47 s. Conclusions Identification of the bacterial meningitis pathogenic species is a critical step for case confirmation and further strain characterization. BMScan employs species-specific thresholds for previously-validated, genome-wide similarity statistics compiled from a curated reference genome collection to rapidly and accurately identify the species of uncharacterized bacterial meningitis pathogens and closely related pathogens. BMScan will facilitate the transition in public health laboratories from traditional phenotypic detection methods to whole genome sequencing based methods for species identification

Genetic Similarity of Gonococcal Homologs to Meningococcal Outer Membrane Proteins of Serogroup B Vaccine

Gonorrhea, a sexually transmitted disease, causes substantial global morbidity and economic burden. New prevention and control measures for this disease are urgently needed, as strains resistant to almost all classes of antibiotics available for treatment have emerged. Previous reports demonstrate that cross-protection from gonococcal infections may be conferred by meningococcal serogroup B (MenB) outer membrane vesicle (OMV)-based vaccines. Among 1,525 common proteins shared across the genomes of both N. gonorrhoeae and N. meningitidis, 57 proteins were predicted to be surface expressed (outer membrane proteins [OMPs]) and thus preferred targets for vaccine development. The majority of these OMPs showed high sequence identity between the 2 bacterial species. Our results provide valuable insight into the meningococcal antigens present in the current OMV-containing MenB-4C vaccine that may contribute to cross-protection against gonorrhea and may inform next steps in gonorrhea vaccine development.The human pathogens Neisseria gonorrhoeae and Neisseria meningitidis share high genome identity. Retrospective analysis of surveillance data from New Zealand indicates the potential cross-protective effect of outer membrane vesicle (OMV) meningococcal serogroup B vaccine (MeNZB) against N. gonorrhoeae. A licensed OMV-based MenB vaccine, MenB-4C, consists of a recombinant FHbp, NhbA, NadA, and the MeNZB OMV. Previous work has identified several abundantly expressed outer membrane proteins (OMPs) as major components of the MenB-4C OMV with high sequence similarity between N. gonorrhoeae and N. meningitidis, suggesting a mechanism for cross-protection. To build off these findings, we performed comparative genomic analysis on 970 recent N. gonorrhoeae isolates collected through a U.S surveillance system against N. meningitidis serogroup B (NmB) reference sequences. We identified 1,525 proteins that were common to both Neisseria species, of which 57 proteins were predicted to be OMPs using in silico methods. Among the MenB-4C antigens, NhbA showed moderate sequence identity (73%) to the respective gonococcal homolog, was highly conserved within N. gonorrhoeae, and was predicted to be surface expressed. In contrast, the gonococcal FHbp was predicted not to be surface expressed, while NadA was absent in all N. gonorrhoeae isolates. Our work confirmed recent observations (E. A. Semchenko, A. Tan, R. Borrow, and K. L. Seib, Clin Infect Dis, 2018, https://doi.org/10.1093/cid/ciy1061) and describes homologous OMPs from a large panel of epidemiologically relevant N. gonorrhoeae strains in the United States against NmB reference strains. Based on our results, we report a set of OMPs that may contribute to the previously observed cross-protection and provide potential antigen targets to guide the next steps in gonorrhea vaccine development

Phylogenetic relationships and regional spread of meningococcal strains in the meningitis belt, 2011–2016

Background: Historically, the major cause of meningococcal epidemics in the meningitis belt of sub-Saharan Africa has been Neisseria meningitidis serogroup A (NmA), but the incidence has been substantially reduced since the introduction of a serogroup A conjugate vaccine starting in 2010. We performed whole-genome sequencing on isolates collected post-2010 to assess their phylogenetic relationships and inter-country transmission. Methods: A total of 716 invasive meningococcal isolates collected between 2011 and 2016 from 11 meningitis belt countries were whole-genome sequenced for molecular characterization by the three WHO Collaborating Centers for Meningitis. Findings: We identified three previously-reported clonal complexes (CC): CC11 (n = 434), CC181 (n = 62) and CC5 (n = 90) primarily associated with NmW, NmX, and NmA, respectively, and an emerging CC10217 (n = 126) associated with NmC. CC11 expanded throughout the meningitis belt independent of the 2000 Hajj outbreak strain, with isolates from Central African countries forming a distinct sub-lineage within this expansion. Two major sub-lineages were identified for CC181 isolates, one mainly expanding in West African countries and the other found in Chad. CC10217 isolates from the large outbreaks in Nigeria and Niger were more closely related than those from the few cases in Mali and Burkina Faso. Interpretations: Whole-genome based phylogenies revealed geographically distinct strain circulation as well as inter-country transmission events. Our results stress the importance of continued meningococcal molecular surveillance in the region, as well as the development of an affordable vaccine targeting these strains. Fund: Meningitis Research Foundation; CDC's Office of Advanced Molecular Detection; GAVI, the Vaccine Alliance. Keywords: Meningitis, Meningitis belt, Invasive meningococcal disease, Phylogenetics, Neisseria meningitidi

Phylogenetic Structure and Comparative Genomics of Multi-National Invasive Haemophilus influenzae Serotype a Isolates

Recent reports have indicated a rise of invasive disease caused by Haemophilus influenzae serotype a (Hia) in North America and some European countries. The whole-genome sequences for a total of 410 invasive Hia isolates were obtained from 12 countries spanning the years of 1998 to 2019 and underwent phylogenetic and comparative genomic analysis in order to characterize the major strains causing disease and the genetic variation present among factors contributing to virulence and antimicrobial resistance. Among 410 isolate sequences received, 408 passed our quality control and underwent genomic analysis. Phylogenetic analysis revealed that the Hia isolates formed four genetically distinct clades: clade 1 (n = 336), clade 2 (n = 13), clade 3 (n = 3) and clade 4 (n = 56). A low diversity subclade 1.1 was found in clade 1 and contained almost exclusively North American isolates. The predominant sequence types in the Hia collection were ST-56 (n = 125), ST-23 (n = 98) and ST-576 (n = 51), which belonged to clade 1, and ST-62 (n = 54), which belonged to clade 4. Clades 1 and 4 contained predominantly North American isolates, and clades 2 and 3 predominantly contained European isolates. Evidence of the presence of capsule duplication was detected in clade 1 and 2 isolates. Seven of the virulence genes involved in endotoxin biosynthesis were absent from all Hia isolates. In general, the presence of known factors contributing to β-lactam antibiotic resistance was low among Hia isolates. Further tests for virulence and antibiotic susceptibility would be required to determine the impact of these variations among the isolates

Phylogenetic Structure and Comparative Genomics of Multi-National Invasive Haemophilus influenzae Serotype a Isolates

Recent reports have indicated a rise of invasive disease caused by Haemophilus influenzae serotype a (Hia) in North America and some European countries. The whole-genome sequences for a total of 410 invasive Hia isolates were obtained from 12 countries spanning the years of 1998 to 2019 and underwent phylogenetic and comparative genomic analysis in order to characterize the major strains causing disease and the genetic variation present among factors contributing to virulence and antimicrobial resistance. Among 410 isolate sequences received, 408 passed our quality control and underwent genomic analysis. Phylogenetic analysis revealed that the Hia isolates formed four genetically distinct clades: clade 1 (n = 336), clade 2 (n = 13), clade 3 (n = 3) and clade 4 (n = 56). A low diversity subclade 1.1 was found in clade 1 and contained almost exclusively North American isolates. The predominant sequence types in the Hia collection were ST-56 (n = 125), ST-23 (n = 98) and ST-576 (n = 51), which belonged to clade 1, and ST-62 (n = 54), which belonged to clade 4. Clades 1 and 4 contained predominantly North American isolates, and clades 2 and 3 predominantly contained European isolates. Evidence of the presence of capsule duplication was detected in clade 1 and 2 isolates. Seven of the virulence genes involved in endotoxin biosynthesis were absent from all Hia isolates. In general, the presence of known factors contributing to β-lactam antibiotic resistance was low among Hia isolates. Further tests for virulence and antibiotic susceptibility would be required to determine the impact of these variations among the isolates.This work was funded by Centers for Disease Control and Prevention.info:eu-repo/semantics/publishedVersio

Phylogenetic Structure and Comparative Genomics of Multi-National Invasive Haemophilus influenzae Serotype a Isolates

Recent reports have indicated a rise of invasive disease caused by Haemophilus influenzae serotype a (Hia) in North America and some European countries. The whole-genome sequences for a total of 410 invasive Hia isolates were obtained from 12 countries spanning the years of 1998 to 2019 and underwent phylogenetic and comparative genomic analysis in order to characterize the major strains causing disease and the genetic variation present among factors contributing to virulence and antimicrobial resistance. Among 410 isolate sequences received, 408 passed our quality control and underwent genomic analysis. Phylogenetic analysis revealed that the Hia isolates formed four genetically distinct clades: clade 1 (n = 336), clade 2 (n = 13), clade 3 (n = 3) and clade 4 (n = 56). A low diversity subclade 1.1 was found in clade 1 and contained almost exclusively North American isolates. The predominant sequence types in the Hia collection were ST-56 (n = 125), ST-23 (n = 98) and ST-576 (n = 51), which belonged to clade 1, and ST-62 (n = 54), which belonged to clade 4. Clades 1 and 4 contained predominantly North American isolates, and clades 2 and 3 predominantly contained European isolates. Evidence of the presence of capsule duplication was detected in clade 1 and 2 isolates. Seven of the virulence genes involved in endotoxin biosynthesis were absent from all Hia isolates. In general, the presence of known factors contributing to β-lactam antibiotic resistance was low among Hia isolates. Further tests for virulence and antibiotic susceptibility would be required to determine the impact of these variations among the isolates