Localized Hypermutation is the Major Driver of Meningococcal Genetic Variability during Persistent Asymptomatic Carriage

High-throughput phenotype-to-genotype testing of meningococcal carriage and disease isolates detects genetic determinants of disease-relevant phenotypic traits

Genome-wide association studies (GWAS) with binary or single phenotype data have successfully identified disease-associated genotypes and determinants of antimicrobial resistance. We describe a novel phenotype-to-genotype approach for a major bacterial pathogen that involves simultaneously testing for associations among multiple disease-related phenotypes and linkages between phenotypic variation and genetic determinants. High-throughput assays quantified variation among 163 Neisseria meningitidis serogroup W ST-11 clonal complex isolates for 11 phenotypic traits. A comparison of carriage and two disease subgroups detected significant differences between groups for eight phenotypic traits. Candidate genotypic testing indicated that indels in csw, a capsular biosynthesis gene, were associated with reduced survival in antibody-depleted heat-inactivated serum. GWAS testing detected 341 significant genetic variants (3 single-nucleotide polymorphisms and 338 unitigs) across all traits except serum bactericidal antibody-depleted assays. Growth traits were associated with variants of capsular biosynthesis genes, carbonic anhydrase, and an iron-uptake system while adhesion-linked variation was in pilC2, marR, and mutS. Multiple phase variation states or combinatorial phasotypes were associated with significant differences in multiple phenotypes. Controlling for group effects through regression and recursive random forest approaches detected group-independent effects for nalP with biofilm formation and fetA with a growth trait. Through random forest testing, nine phenotypes were weakly predictive of MenW:cc11 sub-lineage, original or 2013, for disease isolates while three characteristics separated carriage and disease isolates with >80% accuracy. This study demonstrates the power of combining high-throughput phenotypic testing of pathogenically relevant isolate collections with genomics for identifying genetic determinants of specific disease-relevant phenotypes and the pathobiology of microbial pathogens.IMPORTANCENext-generation sequencing technologies have led to the creation of extensive microbial genome sequence databases for several bacterial pathogens. Mining of these databases is now imperative for unlocking the maximum benefits of these resources. We describe a high-throughput methodology for detecting associations between phenotypic variation in multiple disease-relevant traits and a range of genetic determinants for Neisseria meningitidis, a major causative agent of meningitis and septicemia. Phenotypic variation in 11 disease-related traits was determined for 163 isolates of the hypervirulent ST-11 lineage and linked to specific single-nucleotide polymorphisms, short sequence variants, and phase variation states. Application of machine learning algorithms to our data outputs identified combinatorial phenotypic traits and genetic variants predictive of a disease association. This approach overcomes the limitations of generic meta-data, such as disease versus carriage, and provides an avenue to explore the multi-faceted nature of bacterial disease, carriage, and transmissibility traits.

Microevolution and Its Impact on Hypervirulence, Antimicrobial Resistance, and Vaccine Escape in Neisseria meningitidis

Neisseria meningitidis is commensal of the human pharynx and occasionally invades the host, causing the life-threatening illness invasive meningococcal disease. The meningococcus is a highly diverse and adaptable organism thanks to natural competence, a propensity for recombination, and a highly repetitive genome. These mechanisms together result in a high level of antigenic variation to invade diverse human hosts and evade their innate and adaptive immune responses. This review explores the ways in which this diversity contributes to the evolutionary history and population structure of the meningococcus, with a particular focus on microevolution. It examines studies on meningococcal microevolution in the context of within-host evolution and persistent carriage; microevolution in the context of meningococcal outbreaks and epidemics; and the potential of microevolution to contribute to antimicrobial resistance and vaccine escape. A persistent theme is the idea that the process of microevolution contributes to the development of new hyperinvasive meningococcal variants. As such, microevolution in this species has significant potential to drive future public health threats in the form of hypervirulent, antibiotic-resistant, vaccine-escape variants. The implications of this on current vaccination strategies are explored.

Targeted Hypermutation as a Survival Strategy: A Theoretical Approach

Targeted hypermutation has proven to be a useful survival strategy for bacteria under severe stress and is also used by multicellular organisms in specific instances such as the mammalian immune system. This might appear surprising, given the generally observed deleterious effects of poor replication fidelity/high mutation rate. A previous theoretical model designed to explore the role of replication fidelity in the origin of life was applied to a simulated hypermutation scenario. The results confirmed that the same model is useful for analyzing hypermutation and can predict the effects of the same parameters (survival probability, replication fidelity, mutation effect, and others) on the survival of cellular populations undergoing hypermutation as a result of severe stress.

Variable disruption of epithelial monolayers by Neisseria meningitidis carriage isolates of the hypervirulent MenW cc11 and MenY cc23 lineages

Colonization of mucosal tissues by Neisseria meningitidis requires adhesion mediated by the type IV pilus and multiple outer-membrane proteins. Penetration of the mucosa and invasion of epithelial cells are thought to contribute to host persistence and invasive disease. Using Calu-3 cell monolayers grown at an air-liquid interface, we examined adhesion, invasion and monolayer disruption by carriage isolates of two clonal complexes of N. meningitidis. Carriage isolates of both the serogroup Y cc23 and the hypervirulent serogroup W cc11 lineages exhibited high levels of cellular adhesion, and a variable disruption phenotype across independent isolates. Inactivation of the gene encoding the main pilus sub-unit in multiple cc11 isolates abrogated both adhesive capacity and ability to disrupt epithelial monolayers. Contrastingly, inactivation of the phase-variable opa or nadA genes reduced adhesion and invasion, but not disruption of monolayer integrity. Adherence of tissue-disruptive meningococci correlated with loss of staining for the tight junction protein, occludin. Intriguingly, in a pilus-negative strain background, we observed compensatory ON switching of opa genes, which facilitated continued adhesion. We conclude that disruption of epithelial monolayers occurs in multiple meningococcal lineages but can vary during carriage and is intimately linked to pilus-mediated adhesion.

Whole genome sequence analysis of Neisseria meningitidis strains circulating in Kazakhstan, 2017-2018

Neisseria meningitidis (meningococcus) is a cosmopolitan bacterium that is often found in the upper respiratory tract of asymptomatic humans. However, N. meningitidis also causes meningeal inflammation and/or sepsis in humans with a periodic resurgence in incidence and high mortality rates. The pathogen is highly diverse genetically and antigenically, so that genotyping is considered important for vaccine matching to circulating strains. Annual incidence of meningococcal disease in Kazakhstan ranges between 0.2 and 2.5 cases per 100 thousand population. In total, 78 strains of N. meningitidis were isolated from clinical patients and contact persons during the years 2017-2018 in Kazakhstan. Of these, 41 strains including four from the patients and 37 from contacts, were sequenced using Illumina MiSeq. In silico typing was completed using the Neisseria pipeline 1.2 on the Galaxy Workflow Management System and PubMLST. Whole genome SNP (single nucleotide polymorphisms) trees were built using BioNumerics 8. Seven-gene multilocus sequence typing (MLST) identified ten sequence types (ST), two of which have not been previously described (ST-16025; ST-16027). ST-16025 was detected in two patients with invasive meningococcal disease in 2017 and 2018 in Akmola region and 16 contacts in 2017 in Turkistan region. This prevalent type ST-16025 demonstrates considerable intertypic diversity as it consists of three subcomplexes with a distance of more than 2000 SNPs. Invasive and carrier strains belong to different serogroups (MenB and MenC), PorA and FetA_VR. Two invasive strains were MenB, one MenC and one MenW (Hajj lineage). The strains from the contact persons were: MenC (n = 18), cnl (n = 9), MenY (n = 7), MenW (n = 1), MenB (n = 1) and one unidentifiable. Different numbers of alleles were present: 12, 11, 7, and 7 alleles for PorA, FetA, fHbp, and NHBA, respectively. This study is the first report of the genetic diversity of N. meningitidis strains in Kazakhstan. Despite limitations with the studied sample size, important conclusions can be drawn based on data produced. This study provides evidence for regulatory authorities with regard to changing routine diagnostic protocols to increase the collecting of samples for WGS.

Within-Host Genetic Variation in Neisseria gonorrhoeae over the Course of Infection

Knowledge of within-host genetic variation informs studies on transmission dynamics. We studied within-host genetic variation in Neisseria gonorrhoeae over the course of infection and across different anatomical locations. Isolates were obtained during a clinical trial, and isolates from consecutive time points reflected persistent infections after treatment failure. We compared sequence types (STs) and recombination unfiltered- and filtered core genome single nucleotide polymorphism (SNP) distances in 65 within-host isolate pairs from the same anatomical location over time-obtained with a median interval of 7 days-and 65 isolate pairs across different anatomical locations at one time point. Isolates with different Multi-Locus Sequence Types (MLST), NG-Sequence Types for Antimicrobial Resistance (NG-STAR) and NG-Multi Antigen Sequence Types (NG-MAST) had a median of 1466 recombination filtered SNPs, whereas a median of 1 SNP was found between isolates with identical STs or a different NG-MAST only. The threshold for differentiating between strains was set at 10 recombination filtered SNPs, showing that isolates from persistent infections could have different NG-MASTs. Antibiotic pressure applied through treatment did not lead to an increase in genetic variation in specific genes or in overall extent of variation, compared to variation across anatomical locations. Instead, within-host genetic variation was proposedly driven by the host immune response, as it was concentrated in genomic regions encoding surface exposed proteins involved in host-microbe interaction. Ultimately, 15/228 (6.5%) between-host pairs contained a single strain, suggesting between-host transmission. However, patient reported data are needed to differentiate within-host persistence from between-host transmission. IMPORTANCE Understanding transmission dynamics of Neisseria gonorrhoeae (Ng) is based on the identification of transmission events. These can be identified by assessing genetic relatedness between Ng isolates, expressed as core genome SNP distances. However, a SNP threshold to differentiate between strains needs to be defined, using knowledge on within- and between-host genetic variation. Here, we assessed within-host genetic variation, using a unique set of within-host Ng isolates from the same anatomical location over time or across different anatomical locations at one time point. The insights in genetic variation that occurred during the infection period contribute to the understanding of infection dynamics. In addition, the obtained knowledge can be used for future research on transmission dynamics and development of public health interventions based on bacterial genomic data.

The Host-Pathogen Interactions and Epicellular Lifestyle of Neisseria meningitidis

Neisseria meningitidis is a gram-negative diplococcus and a transient commensal of the human nasopharynx. It shares and competes for this niche with a number of other Neisseria species including N. lactamica, N. cinerea and N. mucosa. Unlike these other members of the genus, N. meningitidis may become invasive, crossing the epithelium of the nasopharynx and entering the bloodstream, where it rapidly proliferates causing a syndrome known as Invasive Meningococcal Disease (IMD). IMD progresses rapidly to cause septic shock and meningitis and is often fatal despite aggressive antibiotic therapy. While many of the ways in which meningococci survive in the host environment have been well studied, recent insights into the interactions between N. meningitidis and the epithelial, serum, and endothelial environments have expanded our understanding of how IMD develops. This review seeks to incorporate recent work into the established model of pathogenesis. In particular, we focus on the competition that N. meningitidis faces in the nasopharynx from other Neisseria species, and how the genetic diversity of the meningococcus contributes to the wide range of inflammatory and pathogenic potentials observed among different lineages.

Neisseria gonorrhoeae physiology and pathogenesis

Neisseria gonorrhoeae is an obligate human pathogen that is the cause of the sexually transmitted disease gonorrhoea. Recently, there has been a surge in gonorrhoea cases that has been exacerbated by the rapid rise in gonococcal multidrug resistance to all useful antimicrobials resulting in this organism becoming a significant public health burden. Therefore, there is a clear and present need to understand the organism's biology through its physiology and pathogenesis to help develop new intervention strategies. The gonococcus initially colonises and adheres to host mucosal surfaces utilising a type IV pilus that helps with microcolony formation. Other adhesion strategies include the porin, PorB, and the phase variable outer membrane protein Opa. The gonococcus is able to subvert complement mediated killing and opsonisation by sialylation of its lipooligosaccharide and deploys a series of anti-phagocytic mechanisms. N. gonorrhoeae is a fastidious organism that is able to grow on a limited number of primary carbon sources such as glucose and lactate. The utilization of lactate by the gonococcus has been implicated in a number of pathogenicity mechanisms. The bacterium lives mainly in microaerobic environments and can grow both aerobically and anaerobically with the aid of nitrite. The gonococcus does not produce siderophores for scavenging iron but can utilize some produced by other bacteria, and it is able to successful chelate iron from host haem, transferrin and lactoferrin. The gonococcus is an incredibly versatile human pathogen; in the following chapter, we detail the intricate mechanisms used by the bacterium to invade and survive within the host.

Modelling evolutionary pathways for commensalism and hypervirulence in Neisseria meningitidis

Neisseria meningitidis, the meningococcus, resides exclusively in humans and causes invasive meningococcal disease (IMD). The population of N. meningitidis is structured into stable clonal complexes by limited horizontal recombination in this naturally transformable species. N. meningitidis is an opportunistic pathogen, with some clonal complexes, such as cc53, effectively acting as commensal colonizers, while other genetic lineages, such as cc11, are rarely colonizers but are over-represented in IMD and are termed hypervirulent. This study examined theoretical evolutionary pathways for pathogenic and commensal lineages by examining the prevalence of horizontally acquired genomic islands (GIs) and loss-of-function (LOF) mutations. Using a collection of 4850 genomes from the BIGSdb database, we identified 82 GIs in the pan-genome of 11 lineages (10 hypervirulent and one commensal lineage). A new computational tool, Phaser, was used to identify frameshift mutations, which were examined for statistically significant association with genetic lineage. Phaser identified a total of 144 frameshift loci of which 105 were shown to have a statistically significant non-random distribution in phase status. The 82 GIs, but not the LOF loci, were associated with genetic lineage and invasiveness using the disease carriage ratio metric. These observations have been integrated into a new model that infers the early events of the evolution of the human adapted meningococcus. These pathways are enriched for GIs that are involved in modulating attachment to the host, growth rate, iron uptake and toxin expression which are proposed to increase competition within the meningococcal population for the limited environmental niche of the human nasopharynx. We surmise that competition for the host mucosal surface with the nasopharyngeal microbiome has led to the selection of isolates with traits that enable access to cell types (non-phagocytic and phagocytic) in the submucosal tissues leading to an increased risk for IMD.

Rapid Transmission of a Hyper-Virulent Meningococcal Clone Due to High Effective Contact Numbers and Super Spreaders

Rapid transmission, a critical contributory factor in outbreaks of invasive meningococcal disease, requires naïve populations of sufficient size and intermingling. We examined genomic variability and transmission dynamics in a student population subject to an 11-fold increase in carriage of a hypervirulent Neisseria meningitidis serogroup W ST-11 clone. Phylogenetic clusters, mutation and recombination rates were derived by bioinformatic analyses of whole-genome sequencing data. Transmission dynamics were determined by combining observed carriage rates, cluster sizes and distributions with simple SIS models. Between 9 and 15 genetically-distinct clusters were detected and associated with seven residential halls. Clusters had low mutation accumulation rates and infrequent recombination events. Modeling indicated that effective contacts decreased from 10 to 2 per day between the start and mid-point of the university term. Transmission rates fluctuated between 1 and 4% while the R(t) for carriage decreased from an initial rate of 47 to 1. Decreases in transmission values correlated with a rise in vaccine-induced immunity. Observed carriage dynamics could be mimicked by populations containing 20% of super spreaders with 2.3-fold higher effective contact rates. We conclude that spread of this hypervirulent ST-11 meningococcal clone depends on the levels of effective contacts and immunity rather than genomic variability. Additionally, we propose that super-spreaders enhance meningococcal transmission and that a 70% MenACWY immunization level is sufficient to retard, but not fully prevent, meningococcal spread in close-contact populations.

Diversification in immunogenicity genes caused by selective pressures in invasive meningococci

We studied population genomics of 486 Neisseria meningitidis isolates causing meningitis in the Netherlands during the period 1979-2003 and 2006-2013 using whole-genome sequencing to evaluate the impact of a hyperendemic period of serogroup B invasive disease. The majority of serogroup B isolates belonged to ST-41/44 (41 %) and ST-32 complex (16 %). Comparing the time periods, before and after the decline of serogroup B invasive disease, there was a decrease of ST-41/44 complex sequences (P=0.002). We observed the expansion of a sub-lineage within ST-41/44 complex sequences being associated with isolation from the 1979-2003 time period (P=0.014). Isolates belonging to this sub-lineage expansion within ST-41/44 complex were marked by four antigen allele variants. Presence of these allele variants was associated with isolation from the 1979-2003 time period after correction for multiple testing (Wald test, P=0.0043 for FetA 1-5; P=0.0035 for FHbp 14; P=0.012 for PorA 7-2.4 and P=0.0031 for NHBA two peptide allele). These sequences were associated with 4CMenB vaccine coverage (Fisher's exact test, P<0.001). Outside of the sub-lineage expansion, isolates with markedly lower levels of predicted vaccine coverage clustered in phylogenetic groups showing a trend towards isolation in the 2006-2013 time period (P=0.08). In conclusion, we show the emergence and decline of a sub-lineage expansion within ST-41/44 complex isolates concurrent with a hyperendemic period in meningococcal meningitis. The expansion was marked by specific antigen peptide allele combinations. We observed preliminary evidence for decreasing 4CMenB vaccine coverage in the post-hyperendemic period.

Current Challenges in Vaccinology

The development of vaccines, which prime the immune system to respond to future infections, has led to global declines in morbidity and mortality from dreadful infectious communicable diseases. However, many pathogens of public health importance are highly complex and/or rapidly evolving, posing unique challenges to vaccine development. Several of these challenges include an incomplete understanding of how immunity develops, host and pathogen genetic variability, and an increased societal skepticism regarding vaccine safety. In particular, new high-dimensional omics technologies, aided by bioinformatics, are driving new vaccine development (vaccinomics). Informed by recent insights into pathogen biology, host genetic diversity, and immunology, the increasing use of genomic approaches is leading to new models and understanding of host immune system responses that may provide solutions in the rapid development of novel vaccine candidates.