Interplay Between Virulence and Variability Factors as a Potential Driver of Invasive Meningococcal Disease

Prevalence, genomic features, and antibiotic sensitivities of Neisseria meningitidis isolates from patients with invasive meningococcal disease and healthy carriers in Zhejiang Province, 2015-2023

OBJECTIVE

A comprehensive understanding of invasive meningococcal disease (IMD) and healthy carriers is critical to monitor, control, and prevent the disease. This study investigated the epidemiology of IMD cases and carriage, and compared population-specific genetic variations and antimicrobial susceptibility of Neisseria meningitidis (N. meningitidis) strains isolated from patients with IMD and carriers.

METHODS

Surveillance data from 2015 to 2023 on patients with epidemic meningitis and healthy carriers in Zhejiang Province, China. We successfully collected 21 isolates from meningitis patients and 16 isolates from healthy individuals during this period. Serogroups of a total of 37 N. meningitidis isolates were determined by polymerase chain reaction (PCR) and slide agglutination, as well as whole genome sequencing to assess various genes, single nucleotide polymorphisms (SNPs), and core-pan genome differences. The antibiotic susceptibility of 37 isolates to 12 antibiotics was evaluated using the E-Test on Mueller-Hinton agar supplemented with 5 % sheep blood.

RESULTS

The annual incidence of IMD and carriage rates remained relatively low from 2015 to 2023. IMD cases were primarily observed in infants under 12 months-of-age. Healthy carriers were predominantly 5-9 and 30-59 years-of-age. Population gene analysis revealed no significant difference in genes between the two groups. Strains of patient and carrier groups were both highly resistant to quinolones and sulfonamides.

CONCLUSIONS

The findings enhance the understanding of N. meningitidis carriage in the context of prevalent invasive meningococcal strains. The findings will facilitate the development and updating of the immunization program of meningitis vaccine, and are critical in understanding the spread and drug use strategies of N. meningitidis.

Whole-genome sequencing resolves biochemical misidentification of Neisseria species from urogenital specimens

Neisseria meningitidis (Nm) and Neisseria gonorrhoeae (Ng) are human pathogens that sometimes occupy the same anatomical niche. Ng, the causative agent of gonorrhea, infects 87 million individuals annually worldwide and is an urgent threat due to increasing drug resistance. Ng is a pathogen of the urogenital tract and may infect the oropharyngeal or rectal site, often asymptomatically. Conversely, Nm is an opportunistic pathogen. While often a commensal in the oropharyngeal tract, it is also the leading cause of bacterial meningitis with 1.2 million cases globally, causing significant morbidity and mortality. Horizontal gene transfer (HGT) is likely to occur between Ng and Nm due to their shared anatomical niches and genetic similarity, which poses challenges for accurate detection and treatment. Routine surveillance through the Gonococcal Isolate Surveillance Project and Strengthening the U.S. Response to Resistant Gonorrhea detected six concerning urogenital Neisseria isolates with contradicting species identification in Milwaukee (MIL). While all six isolates were positive for Ng using nucleic acid amplification testing (NAAT) and matrix-assisted laser desorption/ionization time of flight identified the isolates as Ng, two biochemical tests, Gonochek-II and API NH, classified them as Nm. To address this discrepancy, we performed whole-genome sequencing (WGS) using Illumina MiSeq on all isolates and employed various bioinformatics tools. Species detection analysis using BMScan, which uses WGS data, identified all isolates as Ng. Furthermore, Kraken revealed over 98% of WGS reads mapped to the Ng genome and <1% to Nm. Recombination analysis identified putative HGT in all MIL isolates within the γ-glutamyl transpeptidase (ggt) gene, a key component in the biochemical tests used to differentiate between Nm and Ng. Further analysis identified Nm as the source of HGT event. Specifically, the active Nm ggt gene replaced the Ng pseudogenes, ggt1 and ggt2. Together, this study demonstrates that closely related Neisseria species sharing a niche underwent HGT, which led to the misidentification of species following biochemical testing. Importantly, NAAT accurately detected Ng. The misidentification highlights the importance of using WGS to continually evaluate diagnostic or bacterial identification tests.

Classification of Neisseria meningitidis genomes with a bag-of-words approach and machine learning

Whole genome sequencing of bacteria is important to enable strain classification. Using entire genomes as an input to machine learning (ML) models would allow rapid classification of strains while using information from multiple genetic elements. We developed a "bag-of-words" approach to encode, using SentencePiece or k-mer tokenization, entire bacterial genomes and analyze these with ML. Initial model selection identified SentencePiece with 8,000 and 32,000 words as the best approach for genome tokenization. We then classified in Neisseria meningitidis genomes the capsule B group genotype with 99.6% accuracy and the multifactor invasive phenotype with 90.2% accuracy, in an independent test set. Subsequently, in silico knockouts of 2,808 genes confirmed that the ML model predictions aligned with our current understanding of the underlying biology. To our knowledge, this is the first ML method using entire bacterial genomes to classify strains and identify genes considered relevant by the classifier.

Whole-Genome Sequencing Reveals Differences among Kingella kingae Strains from Carriers and Patients with Invasive Infections

As a result of the increasing use of sensitive nucleic acid amplification tests, Kingella kingae is being recognized as a common pathogen of early childhood, causing medical conditions ranging from asymptomatic oropharyngeal colonization to bacteremia, osteoarthritis, and life-threatening endocarditis. However, the genomic determinants associated with the different clinical outcomes are unknown. Employing whole-genome sequencing, we studied 125 international K. kingae isolates derived from 23 healthy carriers and 102 patients with invasive infections, including bacteremia (n = 23), osteoarthritis (n = 61), and endocarditis (n = 18). We compared their genomic structures and contents to identify genomic determinants associated with the different clinical conditions. The mean genome size of the strains was 2,024,228 bp, and the pangenome comprised 4,026 predicted genes, of which 1,460 (36.3%) were core genes shared by >99% of the isolates. No single gene discriminated between carried and invasive strains; however, 43 genes were significantly more frequent in invasive isolates, compared to asymptomatically carried organisms, and a few showed a significant differential distribution among isolates from skeletal system infections, bacteremia, and endocarditis. The gene encoding the iron-regulated protein FrpC was uniformly absent in all 18 endocarditis-associated strains but was present in one-third of other invasive isolates. Similar to other members of the Neisseriaceae family, the K. kingae differences in invasiveness and tropism for specific body tissues appear to depend on combinations of multiple virulence-associated determinants that are widely distributed throughout the genome. The potential role of the absence of the FrpC protein in the pathogenesis of endocardial invasion deserves further investigation. IMPORTANCE The wide range of clinical severities exhibited by invasive Kingella kingae infections strongly suggests that isolates differ in their genomic contents, and strains associated with life-threatening endocarditis may harbor distinct genomic determinants that result in cardiac tropism and severe tissue damage. The results of the present study show that no single gene discriminated between asymptomatically carried isolates and invasive strains. However, 43 putative genes were significantly more frequent among invasive isolates than among pharyngeal colonizers. In addition, several genes displayed a significant differential distribution among isolates from bacteremia, skeletal system infections, and endocarditis, suggesting that the virulence and tissue tropism of K. kingae are multifactorial and polygenic, depending on changes in the allele content and genomic organization. Further analysis of these putative genes may identify genomic determinants of the invasiveness of K. kingae and its affinity for specific body tissues and potential targets for a future protective vaccine.

Comprehensive Genome Analysis of Neisseria meningitidis from South America Reveals a Distinctive Pathogenicity-Related Prophage Repertoire

Neisseria meningitidis, a bacterium that colonizes in the human nasopharynx, occasionally causes invasive meningococcal disease leading to meningitis or septicemia. Different serogroups and lineages (clonal complexes) are related to the occurrence and epidemiology of N. meningitidis. Despite vaccines for most serogroups, N. meningitidis lineages causing unusual clinical manifestations and a higher fatality rate compared to other lineages have been reported in South America. The present study focused on exploring the diversity of N. meningitidis prophages from South America and their relationship with the epidemiological variables of these strains. We found a high diversity of prophages among the different clonal complexes. By comparing them with previously described N. meningitidis phages and prophages, we revealed groups of prophages sharing similar compositions, which could be useful for prophage comparison in N. meningitidis. Furthermore, we observed a high correlation between the prophage content and epidemiological features, e.g., pathogenicity or clonal complex. Additionally, a distinctive filamentous prophage named here as IMSAR-11 (Invasive Meningococci from South America Related to cc11) was identified. Interestingly, two versions of IMSAR-11, circular and chromosomally integrated, were found. Overall, this study reinforces the importance of the genomic characterization of circulating N. meningitidis lineages to generate new targets for lineage monitoring, diagnosis, or appropriateness of vaccine development. Further studies are necessary to understand the role of these prophages in the persistence, dispersal, and virulence of N. meningitidis in the world.

Neisseria genes required for persistence identified via in vivo screening of a transposon mutant library

The mechanisms used by human adapted commensal Neisseria to shape and maintain a niche in their host are poorly defined. These organisms are common members of the mucosal microbiota and share many putative host interaction factors with Neisseria meningitidis and Neisseria gonorrhoeae. Evaluating the role of these shared factors during host carriage may provide insight into bacterial mechanisms driving both commensalism and asymptomatic infection across the genus. We identified host interaction factors required for niche development and maintenance through in vivo screening of a transposon mutant library of Neisseria musculi, a commensal of wild-caught mice which persistently and asymptomatically colonizes the oral cavity and gut of CAST/EiJ and A/J mice. Approximately 500 candidate genes involved in long-term host interaction were identified. These included homologs of putative N. meningitidis and N. gonorrhoeae virulence factors which have been shown to modulate host interactions in vitro. Importantly, many candidate genes have no assigned function, illustrating how much remains to be learned about Neisseria persistence. Many genes of unknown function are conserved in human adapted Neisseria species; they are likely to provide a gateway for understanding the mechanisms allowing pathogenic and commensal Neisseria to establish and maintain a niche in their natural hosts. Validation of a subset of candidate genes confirmed a role for a polysaccharide capsule in N. musculi persistence but not colonization. Our findings highlight the potential utility of the Neisseria musculi-mouse model as a tool for studying the pathogenic Neisseria; our work represents a first step towards the identification of novel host interaction factors conserved across the genus.

The Neisseria genus contains many genetically related commensals of animals and humans, and two human pathogens, Neisseria gonorrhoeae and Neisseria meningitidis. The mechanisms allowing commensal Neisseria to maintain a niche in their host is little understood. To identify genes required for persistence, we screened a library of transposon mutants of Neisseria musculi, a commensal of wild-caught mice, in CAST/EiJ mice, which persistently and asymptomatically colonizes. Approximately 500 candidate host interaction genes were identified. A subset of these are homologs of N. meningitidis and N. gonorrhoeae genes known to modulate pathogen-host interactions in vitro. Many candidate genes have no known function, demonstrating how much remains to be learned about N. musculi niche maintenance. As many genes of unknown function are conserved in human adapted Neisseria, they provide a gateway for understanding Neisseria persistence mechanisms in general.

Conserved G-Quadruplex Motifs Regulate Gene Expression in Neisseria meningitidis

The noncanonical structures, G-quadruplexes (GQs), formed in the guanine-rich region of nucleic acids regulate various biological and molecular functions in prokaryotes and eukaryotes. Neisseria meningitidis is a commensal residing in a human's upper respiratory tract but occasionally becomes virulent, causing life-threatening septicemia and meningitis. The factors causing these changes in phenotypes are not fully understood. At the molecular level, regulatory components help in a clearer understanding of the pathogen's virulence and pathogenesis. Herein, genome analysis followed by biophysical assays and cell-based experiments revealed the presence of conserved GQ motifs in N. meningitidis. These GQs are linked to the essential genes involved in cell adhesion, pathogenesis, virulence, transport, DNA repair, and recombination. Primer extension stop assay, reporter assays, and quantitative real-time polymerase chain reaction (qRT-PCR) further affirmed the formation of stable GQs in vitro and in vivo. These results support the existence of evolutionarily conserved GQ motifs in N. meningitidis and uphold the usage of GQ-specific ligands as novel antimeningococcal therapeutics.

The minor pilin PilV provides a conserved adhesion site throughout the antigenically variable meningococcal type IV pilus

Neisseria meningitidis utilizes type IV pili (T4P) to adhere to and colonize host endothelial cells, a process at the heart of meningococcal invasive diseases leading to meningitis and sepsis. T4P are polymers of an antigenically variable major pilin building block, PilE, plus several core minor pilins that initiate pilus assembly and are thought to be located at the pilus tip. Adhesion of N. meningitidis to human endothelial cells requires both PilE and a conserved noncore minor pilin PilV, but the localization of PilV and its precise role in this process remains to be clarified. Here, we show that both PilE and PilV promote adhesion to endothelial vessels in vivo. The substantial adhesion defect observed for pilV mutants suggests it is the main adhesin. Consistent with this observation, superresolution microscopy showed the abundant distribution of PilV throughout the pilus. We determined the crystal structure of PilV and modeled it within the pilus filament. The small size of PilV causes it to be recessed relative to adjacent PilE subunits, which are dominated by a prominent hypervariable loop. Nonetheless, we identified a conserved surface-exposed adhesive loop on PilV by alanine scanning mutagenesis. Critically, antibodies directed against PilV inhibit N. meningitidis colonization of human skin grafts. These findings explain how N. meningitidis T4P undergo antigenic variation to evade the humoral immune response while maintaining their adhesive function and establish the potential of this highly conserved minor pilin as a vaccine and therapeutic target for the prevention and treatment of N. meningitidis infections.

Forest and Trees: Exploring Bacterial Virulence with Genome-wide Association Studies and Machine Learning

The advent of inexpensive and rapid sequencing technologies has allowed bacterial whole-genome sequences to be generated at an unprecedented pace. This wealth of information has revealed an unanticipated degree of strain-to-strain genetic diversity within many bacterial species. Awareness of this genetic heterogeneity has corresponded with a greater appreciation of intraspecies variation in virulence. A number of comparative genomic strategies have been developed to link these genotypic and pathogenic differences with the aim of discovering novel virulence factors. Here, we review recent advances in comparative genomic approaches to identify bacterial virulence determinants, with a focus on genome-wide association studies and machine learning.

High coverage of diverse invasive meningococcal serogroup B strains by the 4-component vaccine 4CMenB in Australia, 2007-2011: Concordant predictions between MATS and genetic MATS

Meningococcal serogroup B (MenB) accounts for an important proportion of invasive meningococcal disease (IMD). The 4-component vaccine against MenB (4CMenB) is composed of factor H binding protein (fHbp), neisserial heparin-binding antigen (NHBA), Neisseria adhesin A (NadA), and outer membrane vesicles of the New Zealand strain with Porin 1.4. A meningococcal antigen typing system (MATS) and a fully genomic approach, genetic MATS (gMATS), were developed to predict coverage of MenB strains by 4CMenB. We characterized 520 MenB invasive disease isolates collected over a 5-year period (January 2007-December 2011) from all Australian states/territories by multilocus sequence typing and estimated strain coverage by 4CMenB. The clonal complexes most frequently identified were ST-41/44 CC/Lineage 3 (39.4%) and ST-32 CC/ET-5 CC (23.7%). The overall MATS predicted coverage was 74.6% (95% coverage interval: 61.1%-85.6%). The overall gMATS prediction was 81.0% (lower-upper limit: 75.0-86.9%), showing 91.5% accuracy compared with MATS. Overall, 23.7% and 13.1% (MATS) and 26.0% and 14.0% (gMATS) of isolates were covered by at least 2 and 3 vaccine antigens, respectively, with fHbp and NHBA contributing the most to coverage. When stratified by year of isolate collection, state/territory and age group, MATS and gMATS strain coverage predictions were consistent across all strata. The high coverage predicted by MATS and gMATS indicates that 4CMenB vaccination may have an impact on the burden of MenB-caused IMD in Australia. gMATS can be used in the future to monitor variations in 4CMenB strain coverage over time and geographical areas even for non-culture confirmed IMD cases.

Deconvolution of intergenic polymorphisms determining high expression of Factor H binding protein in meningococcus and their association with invasive disease

Neisseria meningitidis is a strictly human pathogen and is the major cause of septicemia and meningitis worldwide. Factor H binding protein (fHbp) is a meningococcal surface-exposed lipoprotein that binds the human Complement factor H allowing the bacterium to evade the host innate immune response. FHbp is also a key antigen in two vaccines against N. meningitidis serogroup B. Although the fHbp gene is present in most circulating meningococcal strains, level of fHbp expression varies among isolates and has been correlated to differences in promoter sequences upstream of the gene. Here we elucidated the sequence determinants that control fHbp expression in globally circulating strains. We analyzed the upstream fHbpintergenic region (fIR) of more than 5800 strains representative of the UK circulating isolates and we identified eleven fIR sequence alleles which represent 88% of meningococcal strains. By engineering isogenic recombinant strains where fHbp expression was under the control of each of the eleven fIR alleles, we confirmed that the fIR sequence determines a specific and distinct level of expression. Moreover, we identified the molecular basis for variation in expression through polymorphisms within key regulatory regions that are known to affect fHbp expression. We experimentally established three expression groups, high–medium–low, that correlated directly with the susceptibility to killing mediated by anti-fHbp antibodies and the ability of the meningococcal strain to survive within human serum. By using this sequence classification and information about the variant, we predicted fHbp expression in the panel of UK strains and we observed that strains with higher expressing fIR alleles are more likely associated with invasive disease. Overall, our findings can contribute to understand and predict vaccine coverage mediated by fHbp as well as to shed light on the role of this virulence factor in determining an invasive phenotype.

Complement plays a key role in the immunity against Neisseria meningitidis. The meningococcus uses the Factor H binding protein (fHbp), to bind a negative regulator of the alternative complement pathway, factor H, to its surface thus preventing complement deposition and lysis. The use of fHbp as an antigen in two licensed vaccines highlights its public health relevance. Therefore the levels of this antigen produced by the bacterium are pivotal on the one hand for the survival of N. meningitidis in blood and on the other hand for the susceptibility to vaccine-induced killing antibodies. Here, we identify the predominant nucleotide sequences that drive distinct levels of the fHbp antigen in circulating meningococcal strains. We cluster them into distinct groups with increasing levels and observe that strains expressing higher fHbp amounts are associated with invasive disease. Our findings show that the nucleotide sequence of the fHbp promoter can be used for the prediction of antigen levels of any given strain and consequently for both the assessment of its sensitivity to killing by fHbp antibodies and its likelihood to cause invasive disease.

Porphyromonas gingivalis laboratory strains and clinical isolates exhibit different distribution of cell surface and secreted gingipains

Background: The cell-surface cysteine proteinases RgpA, RgpB (Arg-gingipain), and Kgp (Lys-gingipain) are major virulence factors of P. gingivalis, a keystone pathogen in the development of destructive periodontal disease. The gingipains function as proteinases and transpeptidases utilising small peptides such as glycylglycine as acceptor molecules. However, the characteristics of the gingipains from most P. gingivalis strains have not been determined.

Methods: We determined the phenotypes of a panel of P. gingivalis laboratory strains and global clinical isolates with respect to growth on blood agar plus whole-cell and vesicle-free culture supernatant (VFSN) Arg- and Lys-specific proteinase activities.

Results: The P. gingivalis isolates exhibited different growth characteristics and hydrolysis of haemoglobin in solid media. Whole-cell Arg-gingipain Vmax varied 5.8-fold and the whole cell Lys-gingipain V_max varied 2.1-fold across the strains. Furthermore, the P. gingivalis strains showed more than 107-fold variance in soluble Arg-gingipain activity in VFSN and more than 371-fold variance in soluble Lys-gingipain activity in VFSN. Glycylglycine and cysteine stimulated Arg- and Lys-specific cleavage activities of all strains. The stimulation by cysteine was in addition to its redox effect consistent with both glycylglycine and cysteine promoting transpeptidation.

Conclusion: The global P. gingivalis clinical isolates exhibit different Arg- and Lys‑gingipain activities with substantial variability in the level of soluble proteinases released into the environment.

Genetic incorporation of non-canonical amino acid photocrosslinkers in Neisseria meningitidis: New method provides insights into the physiological function of the function-unknown NMB1345 protein

Although whole-genome sequencing has provided novel insights into Neisseria meningitidis, many open reading frames have only been annotated as hypothetical proteins with unknown biological functions. Our previous genetic analyses revealed that the hypothetical protein, NMB1345, plays a crucial role in meningococcal infection in human brain microvascular endothelial cells; however, NMB1345 has no homology to any identified protein in databases and its physiological function could not be elucidated using pre-existing methods. Among the many biological technologies to examine transient protein-protein interaction in vivo, one of the developed methods is genetic code expansion with non-canonical amino acids (ncAAs) utilizing a pyrrolysyl-tRNA synthetase/tRNAPyl pair from Methanosarcina species: However, this method has never been applied to assign function-unknown proteins in pathogenic bacteria. In the present study, we developed a new method to genetically incorporate ncAAs-encoded photocrosslinking probes into N. meningitidis by utilizing a pyrrolysyl-tRNA synthetase/tRNAPyl pair and elucidated the biological function(s) of the NMB1345 protein. The results revealed that the NMB1345 protein directly interacts with PilE, a major component of meningococcal pili, and further physicochemical and genetic analyses showed that the interaction between the NMB1345 protein and PilE was important for both functional pilus formation and meningococcal infectious ability in N. meningitidis. The present study using this new methodology for N. meningitidis provides novel insights into meningococcal pathogenesis by assigning the function of a hypothetical protein.

Insights on Population Structure and Within-Host Genetic Changes among Meningococcal Carriage Isolates from U.S. Universities

In 2015 and 2016, meningococcal carriage evaluations were conducted at two universities in the United States following mass vaccination campaigns in response to Neisseria meningitidis serogroup B (NmB) disease outbreaks. A simultaneous carriage evaluation was also conducted at a university near one of the outbreaks, where no NmB cases were reported and no mass vaccination occurred. A total of ten cross-sectional carriage evaluation rounds were conducted, resulting in 1,514 meningococcal carriage isolates collected from 7,001 unique participants; 1,587 individuals were swabbed at multiple time points (repeat participants). All isolates underwent whole-genome sequencing. The most frequently observed clonal complexes (CC) were CC198 (27.3%), followed by CC1157 (17.4%), CC41/44 (9.8%), CC35 (7.4%), and CC32 (5.6%). Phylogenetic analysis identified carriage isolates that were highly similar to the NmB outbreak strains; comparative genomics between these outbreak and carriage isolates revealed genetic changes in virulence genes. Among repeat participants, 348 individuals carried meningococcal bacteria during at least one carriage evaluation round; 50.3% retained N. meningitidis carriage of a strain with the same sequence type (ST) and CC across rounds, 44.3% only carried N. meningitidis in one round, and 5.4% acquired a new N. meningitidis strain between rounds. Recombination, point mutations, deletions, and simple sequence repeats were the most frequent genetic mechanisms found in isolates collected from hosts carrying a strain of the same ST and CC across rounds. Our findings provide insight on the dynamics of meningococcal carriage among a population that is at higher risk for invasive meningococcal disease than the general population.IMPORTANCE U.S. university students are at a higher risk of invasive meningococcal disease than the general population. The responsible pathogen, Neisseria meningitidis, can be carried asymptomatically in the oropharynx; the dynamics of meningococcal carriage and the genetic features that distinguish carriage versus disease states are not completely understood. Through our analyses, we aimed to provide data to address these topics. We whole-genome sequenced 1,514 meningococcal carriage isolates from individuals at three U.S. universities, two of which underwent mass vaccination campaigns following recent meningococcal outbreaks. We describe the within-host genetic changes among individuals carrying a strain with the same molecular type over time, the primary strains being carried in this population, and the genetic differences between closely related outbreak and carriage strains. Our results provide detailed information on the dynamics of meningococcal carriage and the genetic differences in carriage and outbreak strains, which can inform future efforts to reduce the incidence of invasive meningococcal disease.

Neisseria meningitidis: using genomics to understand diversity, evolution and pathogenesis

Meningococcal disease remains an important cause of morbidity and death worldwide despite the development and increasing implementation of effective vaccines. Elimination of the disease is hampered by the enormous diversity and antigenic variability of the causative agent, Neisseria meningitidis, one of the most variable bacteria in nature. These features are attained mainly through high rates of horizontal gene transfer and alteration of protein expression through phase variation. The recent availability of whole-genome sequencing (WGS) of large-scale collections of N. meningitidis isolates from various origins, databases to facilitate storage and sharing of WGS data and the concomitant development of effective bioinformatics tools have led to a much more thorough understanding of the diversity of the species, its evolution and population structure and how virulent traits may emerge. Implementation of WGS is already contributing to enhanced epidemiological surveillance and is essential to ascertain the impact of vaccination strategies. This Review summarizes the recent advances provided by WGS studies in our understanding of the biology of N. meningitidis and the epidemiology of meningococcal disease.