Genome flexibility in Neisseria meningitidis

Performance Comparison of Two In-House PCR Methods for Detecting Neisseria meningitidis in Asymptomatic Carriers and Antimicrobial Resistance Profiling

Background/Objective: Bacteriological culture has been a widely used method for the detection of meningococcus, but it has low sensitivity and a long turnaround time. Molecular detection targeting capsule transport A (ctrA) gene has been used, but over 16% of meningococcal carriage isolates lack ctrA, resulting in false-negative reports. The Cu-Zn superoxide dismutase gene (sodC) is specific to N. meningitidis, and is not found in other Neisseria species, making it a useful target for improved detection of non-groupable meningococci without intact ctrA. The primary objective of this study was to evaluate the performance compassion of two in-house PCR methods, sodC gene- and ctrA gene-based PCR assays, for detecting N. meningitidis in asymptomatic carriers. The secondary objective was to assess antimicrobial resistance profiling of N. meningitidis isolates. Methods: The performance of sodC gene-based PCR assay compared to ctrA gene-based PCR for detection of N. meningitidis was evaluated using clinical samples (pharyngeal swabs; n = 137) collected from suspected asymptomatic carriers and culture-confirmed meningococci isolates (n = 49). Additionally, the antimicrobial sensitivity of the 49 isolates against antimicrobial drugs was determined using a disk diffusion test. Result: Of 49 DNA samples from culture-positive N. meningitidis isolates, the sodC gene-based PCR accurately identified all 49, whereas the ctrA gene-based PCR identified only 33 out of 49. Of 137 pharyngeal swabs, the sodC gene-based assay detected N. meningitidis DNA in 105 (76.6%), while the ctrA-based assay detected N. meningitidis DNA in 64 (46.7%). Out of the 49 N. meningitidis isolates, 43 (87.8%) were resistant to amoxicillin, 42 (83.7%) to ampicillin, 32 (65.3%) to trimethoprim-sulfamethoxazole, 22 (44.9%) to ceftazidime, 18 (36.7%) to ceftriaxone, and 7 (15.2%) to meropenem. Additionally, the majority of the isolates, 36 (73.5%), were sensitive to cefepime, 31 (63.3%) to ceftriaxone and meropenem, and 26 (53.1%) to ceftazidime. Conclusions: The findings of this study highlight the necessity of adopting non-capsular sodC-based PCR to replace ctrA in resource-constrained laboratories to improve N. meningitidis detection, and underscore the importance of periodic antimicrobial resistance surveillance to inform and adapt treatment strategies.

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.

Genetic characterization of Neisseria meningitidis isolates recovered from patients with invasive meningococcal disease in Lithuania

Introduction

Neisseria meningitidis is a gram-negative bacterium responsible for life-threatening invasive infections known as invasive meningococcal disease and is associated with high fatality rates and serious lifelong disabilities among survivors.

Methods

This study aimed to characterize N. meningitidis isolates cultured from blood and cerebrospinal fluid collected between 2009 and 2021 in Lithuania, assess their genomic relationships with European strains, and evaluate the possibility of using a cost-effective method for strain characterization, thus improving the national molecular surveillance of invasive meningococcal disease. In total, 321 N. meningitidis isolates were collected and analyzed using multilocus restriction typing (MLRT). Amplification of the penA gene and restriction fragment length polymorphism analysis were performed to identify the modified penA genes. Based on the MLRT genotyping results, we selected 10 strains for additional analysis using whole-genome sequencing. The sequenced genomes were incorporated into a dataset of publicly available N. meningitidis genomes to evaluate genomic diversity and establish phylogenetic relationships within the Lithuanian and European circulating strains.

Results

We identified 83 different strains using MLRT genotyping. Genomic diversity of N. meningitidis genomes analysed revealed 21 different sequence types (STs) circulating in Lithuania. Among these, ST34 was the most prevalent. Notably, three isolates displayed unique combinations of seven housekeeping genes and were identified as novel STs: ST16969, ST16901, and ST16959. The analyzed strains were found to possess virulence factors not commonly found in N. meningitidis. Six distinct penA profiles were identified, each with different frequencies. In the present study, we also identified N. meningitidis strains with new penA, NEIS0123, NEIS1320, NEIS1525, NEIS1600, and NEIS1753 loci variants. In our study, using the cgMLST scheme, Minimum Spanning Tree (MST) analysis did not identify significant geographic relationships between Lithuanian N. meningitidis isolates and strains from Europe.

Discussion

Discussion: To our knowledge, this is the first study to employ whole genome sequencing (WGS) method for a comprehensive genetic characterization of invasive N. meningitidis isolates from Lithuania. This approach provides a more detailed and precise analysis of genomic relationships and diversity compared to prior studies relying on traditional molecular typing methods and antigen analysis.

Loss to gain: pseudogenes in microorganisms, focusing on eubacteria, and their biological significance

Pseudogenes are defined as "non-functional" copies of corresponding parent genes. The cognition of pseudogenes continues to be refreshed through accumulating and updating research findings. Previous studies have predominantly focused on mammals, but pseudogenes have received relatively less attention in the field of microbiology. Given the increasing recognition on the importance of pseudogenes, in this review, we focus on several aspects of microorganism pseudogenes, including their classification and characteristics, their generation and fate, their identification, their abundance and distribution, their impact on virulence, their ability to recombine with functional genes, the extent to which some pseudogenes are transcribed and translated, and the relationship between pseudogenes and viruses. By summarizing and organizing the latest research progress, this review will provide a comprehensive perspective and improved understanding on pseudogenes in microorganisms. KEY POINTS: • Concept, classification and characteristics, identification and databases, content, and distribution of microbial pseudogenes are presented. • How pseudogenization contribute to pathogen virulence is highlighted. • Pseudogenes with potential functions in microorganisms are discussed.

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.

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.

Pangenome graphs in infectious disease: a comprehensive genetic variation analysis of Neisseria meningitidis leveraging Oxford Nanopore long reads

Whole genome sequencing has revolutionized infectious disease surveillance for tracking and monitoring the spread and evolution of pathogens. However, using a linear reference genome for genomic analyses may introduce biases, especially when studies are conducted on highly variable bacterial genomes of the same species. Pangenome graphs provide an efficient model for representing and analyzing multiple genomes and their variants as a graph structure that includes all types of variations. In this study, we present a practical bioinformatics pipeline that employs the PanGenome Graph Builder and the Variation Graph toolkit to build pangenomes from assembled genomes, align whole genome sequencing data and call variants against a graph reference. The pangenome graph enables the identification of structural variants, rearrangements, and small variants (e.g., single nucleotide polymorphisms and insertions/deletions) simultaneously. We demonstrate that using a pangenome graph, instead of a single linear reference genome, improves mapping rates and variant calling for both simulated and real datasets of the pathogen Neisseria meningitidis. Overall, pangenome graphs offer a promising approach for comparative genomics and comprehensive genetic variation analysis in infectious disease. Moreover, this innovative pipeline, leveraging pangenome graphs, can bridge variant analysis, genome assembly, population genetics, and evolutionary biology, expanding the reach of genomic understanding and applications.

Genomic Diversity and Chromosomal Rearrangements in Neisseria gonorrhoeae and Neisseria meningitidis

Chromosomal rearrangements in N. gonorrhoeae and N. meningitidis were studied with the determination of mobile elements and their role in rearrangements. The results of whole-genome sequencing and de novo genome assembly for 50 N. gonorrhoeae isolates collected in Russia were compared with 96 genomes of N. gonorrhoeae and 138 genomes of N. meningitidis from the databases. Rearrangement events with the determination of the coordinates of syntenic blocks were analyzed using the SibeliaZ software v.1.2.5, the minimum number of events that allow one genome to pass into another was calculated using the DCJ-indel model using the UniMoG program v.1.0. Population-level analysis revealed a stronger correlation between changes in the gene order and phylogenetic proximity for N. meningitidis in contrast to N. gonorrhoeae. Mobile elements were identified, including Correa elements; Spencer-Smith elements (in N. gonorrhoeae); Neisserial intergenic mosaic elements; IS elements of IS5, IS30, IS110, IS1595 groups; Nf1-Nf3 prophages; NgoФ1-NgoФ9 prophages; and Mu-like prophages Pnm1, Pnm2, MuMenB (in N. meningitidis). More than 44% of the observed rearrangements most likely occurred with the participation of mobile elements, including prophages. No differences were found between the Russian and global N. gonorrhoeae population both in terms of rearrangement events and in the number of transposable elements in genomes.

The effect of recombination on the evolution of a population of Neisseria meningitidis

Neisseria meningitidis (the meningococcus) is a major human pathogen with a history of high invasive disease burden, particularly in sub-Saharan Africa. Our current understanding of the evolution of meningococcal genomes is limited by the rarity of large-scale genomic population studies and lack of in-depth investigation of the genomic events associated with routine pathogen transmission. Here, we fill this knowledge gap by a detailed analysis of 2839 meningococcal genomes obtained through a carriage study of over 50,000 samples collected systematically in Burkina Faso, West Africa, before, during, and after the serogroup A vaccine rollout, 2009-2012. Our findings indicate that the meningococcal genome is highly dynamic, with highly recombinant loci and frequent gene sharing across deeply separated lineages in a structured population. Furthermore, our findings illustrate how population structure can correlate with genome flexibility, as some lineages in Burkina Faso are orders of magnitude more recombinant than others. We also examine the effect of selection on the population, in particular how it is correlated with recombination. We find that recombination principally acts to prevent the accumulation of deleterious mutations, although we do also find an example of recombination acting to speed the adaptation of a gene. In general, we show the importance of recombination in the evolution of a geographically expansive population with deep population structure in a short timescale. This has important consequences for our ability to both foresee the outcomes of vaccination programs and, using surveillance data, predict when lineages of the meningococcus are likely to become a public health concern.

A Narrative Review of the Molecular Epidemiology and Laboratory Surveillance of Vaccine Preventable Bacterial Meningitis Agents: Streptococcus pneumoniae, Neisseria meningitidis, Haemophilus influenzae and Streptococcus agalactiae

This narrative review describes the public health importance of four most common bacterial meningitis agents, Streptococcus pneumoniae, Neisseria meningitidis, Haemophilus influenzae, and S. agalactiae (group B Streptococcus). Three of them are strict human pathogens that normally colonize the nasopharynx and may invade the blood stream to cause systemic infections and meningitis. S. agalactiae colonizes the genito-gastrointestinal tract and is an important meningitis agent in newborns, but also causes invasive infections in infants or adults. These four bacteria have polysaccharide capsules that protect them against the host complement defense. Currently licensed conjugate vaccines (against S. pneumoniae, H. influenza, and N. meningitidis only but not S. agalactiae) can induce protective serum antibodies in infants as young as two months old offering protection to the most vulnerable groups, and the ability to eliminate carriage of homologous serotype strains in vaccinated subjects lending further protection to those not vaccinated through herd immunity. However, the serotype-specific nature of these vaccines have driven the bacteria to adapt by mechanisms that affect the capsule antigens through either capsule switching or capsule replacement in addition to the possibility of unmasking of strains or serotypes not covered by the vaccines. The post-vaccine molecular epidemiology of vaccine-preventable bacterial meningitis is discussed based on findings obtained with newer genomic laboratory surveillance methods.

Identification of integrative and conjugative elements in pathogenic and commensal Neisseriaceae species via genomic distributions of DNA uptake sequence dialects

Mobile genetic elements (MGEs) are key factors responsible for dissemination of virulence determinants and antimicrobial-resistance genes amongst pathogenic bacteria. Conjugative MGEs are notable for their high gene loads donated per transfer event, broad host ranges and phylogenetic ubiquity amongst prokaryotes, with the subclass of chromosomally inserted integrative and conjugative elements (ICEs) being particularly abundant. The focus on a small number of model systems has biased the study of ICEs towards those conferring readily selectable phenotypes to host cells, whereas the identification and characterization of integrated cryptic elements remains challenging. Even though antimicrobial resistance and horizontally acquired virulence genes are major factors aggravating neisserial infection, conjugative MGEs of Neisseria gonorrhoeae and Neisseria meningitidis remain poorly characterized. Using a phenotype-independent approach based on atypical distributions of DNA uptake sequences (DUSs) in MGEs relative to the chromosomal background, we have identified two groups of chromosomally integrated conjugative elements in Neisseria: one found almost exclusively in pathogenic species possibly deriving from the genus Kingella, the other belonging to a group of Neisseria mucosa-like commensals. The former element appears to enable transfer of traditionally gonococcal-specific loci such as the virulence-associated toxin-antitoxin system fitAB to N. meningitidis chromosomes, whilst the circular form of the latter possesses a unique attachment site (attP) sequence seemingly adapted to exploit DUS motifs as chromosomal integration sites. In addition to validating the use of DUS distributions in Neisseriaceae MGE identification, the >170 identified ICE sequences provide a valuable resource for future studies of ICE evolution and host adaptation.

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

Host persistence of bacteria is facilitated by mutational and recombinatorial processes that counteract loss of genetic variation during transmission and selection from evolving host responses. Genetic variation was investigated during persistent asymptomatic carriage of Neisseria meningitidis Interrogation of whole-genome sequences for paired isolates from 25 carriers showed that de novo mutations were infrequent, while horizontal gene transfer occurred in 16% of carriers. Examination of multiple isolates per time point enabled separation of sporadic and transient allelic variation from directional variation. A comprehensive comparative analysis of directional allelic variation with hypermutation of simple sequence repeats and hyperrecombination of class 1 type IV pilus genes detected an average of seven events per carrier and 2:1 bias for changes due to localized hypermutation. Directional genetic variation was focused on the outer membrane with 69% of events occurring in genes encoding enzymatic modifiers of surface structures or outer membrane proteins. Multiple carriers exhibited directional and opposed switching of allelic variants of the surface-located Opa proteins that enables continuous expression of these adhesins alongside antigenic variation. A trend for switching from PilC1 to PilC2 expression was detected, indicating selection for specific alterations in the activities of the type IV pilus, whereas phase variation of restriction modification (RM) systems, as well as associated phasevarions, was infrequent. We conclude that asymptomatic meningococcal carriage on mucosal surfaces is facilitated by frequent localized hypermutation and horizontal gene transfer affecting genes encoding surface modifiers such that optimization of adhesive functions occurs alongside escape of immune responses by antigenic variation.IMPORTANCE Many bacterial pathogens coexist with host organisms, rarely causing disease while adapting to host responses. Neisseria meningitidis, a major cause of meningitis and septicemia, is a frequent persistent colonizer of asymptomatic teenagers/young adults. To assess how genetic variation contributes to host persistence, whole-genome sequencing and hypermutable sequence analyses were performed on multiple isolates obtained from students naturally colonized with meningococci. High frequencies of gene transfer were observed, occurring in 16% of carriers and affecting 51% of all nonhypermutable variable genes. Comparative analyses showed that hypermutable sequences were the major mechanism of variation, causing 2-fold more changes in gene function than other mechanisms. Genetic variation was focused on genes affecting the outer membrane, with directional changes in proteins responsible for bacterial adhesion to host surfaces. This comprehensive examination of genetic plasticity in individual hosts provides a significant new platform for rationale design of approaches to prevent the spread of this pathogen.

Continuous production of Neisseria meningitidis outer membrane vesicles

Outer membrane vesicles (OMVs) are nanoparticles secreted by Gram-negative bacteria that can be used for diverse biotechnological applications. Interesting applications have been developed, where OMVs are the basis of drug delivery, enzyme carriers, adjuvants, and vaccines. Historically, OMV research has mainly focused on vaccines. Therefore, current OMV production processes have been based on batch processes. The production of OMVs in batch mode is characterized by relatively low yields and high costs. Transition of OMV production processes from batch to continuous processes could increase the volumetric productivity, reduce the production and capital costs, and result in a higher quality product. Here, we study the continuous production of Neisseria meningitidis OMVs to improve volumetric productivity. Continuous cultivation of N. meningitidis resulted in a steady state with similar high OMV concentrations as are reached in current batch processes. The steady state was reproducible and could be maintained for at least 600 h. The volumetric productivity of a continuous culture reached 4.0 × 1014 OMVs per liter culture per day, based on a dilution rate of 1/day. The tested characteristics of the OMVs did not change during the experiments showing feasibility of a continuous production process for the production of OMVs for any application.

Whole-genome sequencing and characterization of an antibiotic resistant Neisseria meningitidis B isolate from a military unit in Vietnam

Background

Invasive meningococcal disease (IMD) persists in military units in Vietnam despite the availability of antibiotics and vaccines. A hindrance to reducing the incidence of IMD in Vietnam is a lack of molecular data from isolates of the causative agent, Neisseria meningitidis from this country. Here, we characterized key genetic and epidemiological features of an invasive N. meningitidis isolate from a military unit in Vietnam using whole-genome sequencing.

Methods

Neisseria meningitidis was isolated from a conscript admitted for meningitis and tested against seven antibiotics. DNA from the isolate was extracted and sequenced using the Illumina HiSeq platform. Denovo assembly and scaffolding were performed to construct a draft genome assembly, from which genes were predicted and functionally annotated. Genome analysis included epidemiological characterization, genomic composition and identification of antibiotic resistance genes.

Results

Susceptibility testing of the isolate showed high levels of resistance to chloramphenicol and diminished susceptibility to ampicillin and rifampicin. A draft genome of ~ 2.1 Mb was assembled, containing 2451 protein coding sequences, 49 tRNAs and 3 rRNAs. Fifteen coding sequences sharing ≥ 84% identity with known antibiotic resistance genes were identified. Genome analysis revealed abundant repetitive DNAs and two prophages. Epidemiological typing revealed newly described sequence type, antigenic finetype and Bexsero^® Antigen Sequence Typing (BAST). The BAST profile showed no coverage by either Bexsero^® or Trumenba^®.

Conclusions

Our results present the first genome assembly of an invasive N. meningitidis isolate from a military unit in Vietnam. This study illustrates the usefulness of whole genome sequencing (WGS) analysis for epidemiological and antibiotic resistance studies and surveillance of IMD in Vietnam.

The Meningococcal Cysteine Transport System Plays a Crucial Role in Neisseria meningitidis Survival in Human Brain Microvascular Endothelial Cells

Neisseria meningitidis colonizes at a nasopharynx of human as a unique host and has many strains that are auxotrophs for amino acids for their growth. To cause invasive meningococcal diseases (IMD) such as sepsis and meningitis, N. meningitidis passes through epithelial and endothelial barriers and infiltrates into blood and cerebrospinal fluid as well as epithelial and endothelial cells. However, meningococcal nutrients, including cysteine, become less abundant when it more deeply infiltrates the human body even during inflammation, such that N. meningitidis has to acquire nutrients in order to survive/persist, disseminate, and proliferate in humans. This was the first study to examine the relationship between meningococcal cysteine acquisition and the pathogenesis of meningococcal infections. The results of the present study provide insights into the mechanisms by which pathogens with auxotrophs acquire nutrients in hosts and may also contribute to the development of treatments and prevention strategies for IMD.

While Neisseria meningitidis typically exists in an asymptomatic nasopharyngeal carriage state, it may cause potentially lethal diseases in humans, such as septicemia or meningitis, by invading deeper sites in the body. Since the nutrient compositions of human cells are not always conducive to meningococci, N. meningitidis needs to exploit nutrients from host environments. In the present study, the utilization of cysteine by the meningococcal cysteine transport system (CTS) was analyzed for the pathogenesis of meningococcal infections. A N. meningitidis strain deficient in one of the three cts genes annotated as encoding cysteine-binding protein (cbp) exhibited approximately 100-fold less internalization into human brain microvascular endothelial cells (HBMEC) than the wild-type strain. This deficiency was restored by complementation with the three cts genes together, and the infectious phenotype of HBMEC internalization correlated with cysteine uptake activity. However, efficient accumulation of ezrin was observed beneath the cbp mutant. The intracellular survival of the cbp mutant in HBMEC was markedly reduced, whereas equivalent reductions of glutathione concentrations and of resistance to reactive oxygens species in the cbp mutant were not found. The cbp mutant grew well in complete medium but not in synthetic medium supplemented with less than 300 μM cysteine. Taking cysteine concentrations in human cells and other body fluids, including blood and cerebrospinal fluid, into consideration, the present results collectively suggest that the meningococcal CTS is crucial for the acquisition of cysteine from human cells and participates in meningococcal nutrient virulence.

Tracking a serial killer: Integrating phylogenetic relationships, epidemiology, and geography for two invasive meningococcal disease outbreaks

BACKGROUND

While overall rates of meningococcal disease have been declining in the United States for the past several decades, New York City (NYC) has experienced two serogroup C meningococcal disease outbreaks in 2005-2006 and in 2010-2013. The outbreaks were centered within drug use and sexual networks, were difficult to control, and required vaccine campaigns.

METHODS

Whole Genome Sequencing (WGS) was used to analyze preserved meningococcal isolates collected before and during the two outbreaks. We integrated and analyzed epidemiologic, geographic, and genomic data to better understand transmission networks among patients. Betweenness centrality was used as a metric to understand the most important geographic nodes in the transmission networks. Comparative genomics was used to identify genes associated with the outbreaks.

RESULTS

Neisseria meningitidis serogroup C (ST11/ET-37) was responsible for both outbreaks with each outbreak having distinct phylogenetic clusters. WGS did identify some misclassifications of isolates that were more distant from the outbreak strains, as well as those that should have been included based on high genomic similarity. Genomes for the second outbreak were more similar than the first and no polymorphism was found to either be unique or specific to either outbreak lineage. Betweenness centrality as applied to transmission networks based on phylogenetic analysis demonstrated that the outbreaks were transmitted within focal communities in NYC with few transmission events to other locations.

CONCLUSIONS

Neisseria meningitidis is an ever changing pathogen and comparative genomic analyses can help elucidate how it spreads geographically to facilitate targeted interventions to interrupt transmission.

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

Neisseria meningitidis (Nm) is frequently found in the upper respiratory tract of the human population. Despite its prevalence as a commensal organism, Nm can occasionally invade the pharyngeal mucosal epithelium causing septicemia and life-threatening disease. A number of studies have tried to identify factors that are responsible for the onset of a virulent phenotype. Despite this however, we still miss clear causative elements. Several factors have been identified to be associated to an increased susceptibility to meningococcal disease in humans. None of them, however, could unambiguously discriminate healthy carrier from infected individuals. Similarly, comparative studies of virulent and apathogenic strains failed to identify virulence factors that could explain the emergence of the pathogenic phenotype. In line with this, a recent study of within host evolution found that Nm accumulates genomic changes during the asymptomatic carriage phase and that these are likely to contribute to the shift to a pathogenic phenotype. These results suggest that the presence of virulence factors in the meningococcal genome is not a sufficient condition for developing virulent traits, but is rather the ability to promote phenotypic variation, through the stochastic assortment of the repertoire of such factors, which could explain the occasional and unpredictable onset of IMD. Here, we present a series of argumentations supporting the hypothesis that invasive meningococcal disease comes as a result of the coexistence of bacterial virulence and variability factors in a plot that can be further complicated by additional latent factors, like host pre-existing immune status and genetic predisposition.

Survival of Neisseria meningitidis outside of the host: environmental effects and differences among strains

Neisseria meningitidis is a gram-negative bacterium that lives as a commensal in the human nasopharynx. Meningococci are generally non-invasive, but can invade the nasopharyngeal epithelia and enter the bloodstream causing life-threatening illnesses. It is generally thought that meningococci do not survive for long outside the host, and that transmission requires relatively close contact between hosts. There are some reports, however, that meningococci can survive drying on surfaces, including glass, plastic and cloth. Our examination of N. meningitidis strains dried on glass showed differences in survival of isolates belonging to serogroups B, C and W135, including persistence of Cuban, New Zealand, and Norwegian epidemic strains up to 8 days, depending on temperature and humidity. Survival of a New Zealand epidemic strain isolate NZ98/254 under ambient conditions in the laboratory was greatest in winter suggesting that environmental factors impacted survival. For most isolates, including NZ98/254, survival under controlled conditions at 30 °C was greater at 22% than 30% relative humidity. There were also some differences in survival between carriage and invasive strains. The results suggest that N. meningitidis could be transmitted through contact with surfaces outside the host, potentially including contact through shared drinking vessels.

The primary transcriptome of Neisseria meningitidis and its interaction with the RNA chaperone Hfq

Neisseria meningitidis is a human commensal that can also cause life-threatening meningitis and septicemia. Despite growing evidence for RNA-based regulation in meningococci, their transcriptome structure and output of regulatory small RNAs (sRNAs) are incompletely understood. Using dRNA-seq, we have mapped at single-nucleotide resolution the primary transcriptome of N. meningitidis strain 8013. Annotation of 1625 transcriptional start sites defines transcription units for most protein-coding genes but also reveals a paucity of classical σ70-type promoters, suggesting the existence of activators that compensate for the lack of −35 consensus sequences in N. meningitidis. The transcriptome maps also reveal 65 candidate sRNAs, a third of which were validated by northern blot analysis. Immunoprecipitation with the RNA chaperone Hfq drafts an unexpectedly large post-transcriptional regulatory network in this organism, comprising 23 sRNAs and hundreds of potential mRNA targets. Based on this data, using a newly developed gfp reporter system we validate an Hfq-dependent mRNA repression of the putative colonization factor PrpB by the two trans-acting sRNAs RcoF1/2. Our genome-wide RNA compendium will allow for a better understanding of meningococcal transcriptome organization and riboregulation with implications for colonization of the human nasopharynx.

Investigation of potential targets of Porphyromonas CRISPRs among the genomes of Porphyromonas species

The oral bacterial species Porphyromonas gingivalis, a periodontal pathogen, has plastic genomes that may be driven by homologous recombination with exogenous deoxyribonucleic acid (DNA) that is incorporated by natural transformation and conjugation. However, bacteriophages and plasmids, both of which are main resources of exogenous DNA, do not exist in the known P. gingivalis genomes. This could be associated with an adaptive immunity system conferred by clustered regularly interspaced short palindromic repeat (CRISPR) and CRISPR-associated (cas) genes in P. gingivalis as well as innate immune systems such as a restriction-modification system. In a previous study, few immune targets were predicted for P. gingivalis CRISPR/Cas. In this paper, we analyzed 51 P. gingivalis genomes, which were newly sequenced, and publicly available genomes of 13 P. gingivalis and 46 other Porphyromonas species. We detected 6 CRISPR/Cas types (classified by sequence similarity of repeat) in P. gingivalis and 12 other types in the remaining species. The Porphyromonas CRISPR spacers with potential targets in the genus Porphyromonas were approximately 23 times more abundant than those with potential targets in other genus taxa (1,720/6,896 spacers vs. 74/6,896 spacers). Porphyromonas CRISPR/Cas may be involved in genome plasticity by exhibiting selective interference against intra- and interspecies nucleic acids.

Transcriptomic buffering of cryptic genetic variation contributes to meningococcal virulence

Background

Commensal bacteria like Neisseria meningitidis sometimes cause serious disease. However, genomic comparison of hyperinvasive and apathogenic lineages did not reveal unambiguous hints towards indispensable virulence factors. Here, in a systems biological approach we compared gene expression of the invasive strain MC58 and the carriage strain α522 under different ex vivo conditions mimicking commensal and virulence compartments to assess the strain-specific impact of gene regulation on meningococcal virulence.

Results

Despite indistinguishable ex vivo phenotypes, both strains differed in the expression of over 500 genes under infection mimicking conditions. These differences comprised in particular metabolic and information processing genes as well as genes known to be involved in host-damage such as the nitrite reductase and numerous LOS biosynthesis genes. A model based analysis of the transcriptomic differences in human blood suggested ensuing metabolic flux differences in energy, glutamine and cysteine metabolic pathways along with differences in the activation of the stringent response in both strains. In support of the computational findings, experimental analyses revealed differences in cysteine and glutamine auxotrophy in both strains as well as a strain and condition dependent essentiality of the (p)ppGpp synthetase gene relA and of a short non-coding AT-rich repeat element in its promoter region.

Conclusions

Our data suggest that meningococcal virulence is linked to transcriptional buffering of cryptic genetic variation in metabolic genes including global stress responses. They further highlight the role of regulatory elements for bacterial virulence and the limitations of model strain approaches when studying such genetically diverse species as N. meningitidis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3616-7) contains supplementary material, which is available to authorized users.

How the Knowledge of Interactions between Meningococcus and the Human Immune System Has Been Used to Prepare Effective Neisseria meningitidis Vaccines

In the last decades, tremendous advancement in dissecting the mechanisms of pathogenicity of Neisseria meningitidis at a molecular level has been achieved, exploiting converging approaches of different disciplines, ranging from pathology to microbiology, immunology, and omics sciences (such as genomics and proteomics). Here, we review the molecular biology of the infectious agent and, in particular, its interactions with the immune system, focusing on both the innate and the adaptive responses. Meningococci exploit different mechanisms and complex machineries in order to subvert the immune system and to avoid being killed. Capsular polysaccharide and lipooligosaccharide glycan composition, in particular, play a major role in circumventing immune response. The understanding of these mechanisms has opened new horizons in the field of vaccinology. Nowadays different licensed meningococcal vaccines are available and used: conjugate meningococcal C vaccines, tetravalent conjugate vaccines, an affordable conjugate vaccine against the N. menigitidis serogroup A, and universal vaccines based on multiple antigens each one with a different and peculiar function against meningococcal group B strains.

Protein Mis-Termination Initiates Genetic Diseases, Cancers, and Restricts Bacterial Genome Expansion

Protein termination is an important cellular process. Protein termination relies on the stop-codons in the mRNA interacting properly with the releasing factors on the ribosome. One third of inherited diseases, including cancers, are associated with the mutation of the stop-codons. Many pathogens and viruses are able to manipulate their stop-codons to express their virulence. The influence of stop-codons is not limited to the primary reading frame of the genes. Stop-codons in the second and third reading frames are referred as premature stop signals (PSC). Stop-codons and PSCs together are collectively referred as stop-signals. The ratios of the stop-signals (referred as translation stop-signals ratio or TSSR) of genetically related bacteria, despite their great differences in gene contents, are much alike. This nearly identical Genomic-TSSR value of genetically related bacteria may suggest that bacterial genome expansion is limited by their unique stop-signals bias. We review the protein termination process and the different types of stop-codon mutation in plants, animals, microbes, and viruses, with special emphasis on the role of PSCs in directing bacterial evolution in their natural environments. Knowing the limit of genomic boundary could facilitate the formulation of new strategies in controlling the spread of diseases and combat antibiotic-resistant bacteria.

Genomic comparative analysis and gene function prediction in infectious diseases: application to the investigation of a meningitis outbreak

Background

Next generation sequencing (NGS) is being increasingly used for the detection and characterization of pathogens during outbreaks. This technology allows rapid sequencing of pathogen full genomes, useful not only for accurate genotyping and molecular epidemiology, but also for identification of drug resistance and virulence traits.

Methods

In this study, an approach based on whole genome sequencing by NGS, comparative genomics, and gene function prediction was set up and retrospectively applied for the investigation of two N. meningitidis serogroup C isolates collected from a cluster of meningococcal disease, characterized by a high fatality rate.

Results

According to conventional molecular typing methods, all the isolates had the same typing results and were classified as outbreak isolates within the same N. meningitidis sequence type ST-11, while full genome sequencing demonstrated subtle genetic differences between the isolates. Looking for these specific regions by means of 9 PCR and cycle sequencing assays in other 7 isolates allowed distinguishing outbreak cases from unrelated cases. Comparative genomics and gene function prediction analyses between outbreak isolates and a set of reference N. meningitidis genomes led to the identification of differences in gene content that could be relevant for pathogenesis. Most genetic changes occurred in the capsule locus and were consistent with recombination and horizontal acquisition of a set of genes involved in capsule biosynthesis.

Conclusions

This study showed the added value given by whole genome sequencing by NGS over conventional sequence-based typing methods in the investigation of an outbreak. Routine application of this technology in clinical microbiology will significantly improve methods for molecular epidemiology and surveillance of infectious disease and provide a bulk of data useful to improve our understanding of pathogens biology.

Reverse vaccinology in the 21st century: improvements over the original design

Reverse vaccinology (RV), the first application of genomic technologies in vaccine research, represented a major revolution in the process of discovering novel vaccines. By determining their entire antigenic repertoire, researchers could identify protective targets and design efficacious vaccines for pathogens where conventional approaches had failed. Bexsero, the first vaccine developed using RV, has recently received positive opinion from the European Medicines Agency. The use of RV initiated a cascade of changes that affected the entire vaccine development process, shifting the focus from the identification of a list of vaccine candidates to the definition of a set of high throughput screens to reduce the need for costly and labor intensive tests in animal models. It is now clear that a deep understanding of the epidemiology of vaccine candidates, and their regulation and role in host-pathogen interactions, must become an integral component of the screening workflow. Far from being outdated by technological advancements, RV still represents a paradigm of how high-throughput technologies and scientific insight can be integrated into biotechnology research.

Attachment and invasion of Neisseria meningitidis to host cells is related to surface hydrophobicity, bacterial cell size and capsule

We compared exemplar strains from two hypervirulent clonal complexes, strain NMB-CDC from ST-8/11 cc and strain MC58 from ST-32/269 cc, in host cell attachment and invasion. Strain NMB-CDC attached to and invaded host cells at a significantly greater frequency than strain MC58. Type IV pili retained the primary role for initial attachment to host cells for both isolates regardless of pilin class and glycosylation pattern. In strain MC58, the serogroup B capsule was the major inhibitory determinant affecting both bacterial attachment to and invasion of host cells. Removal of terminal sialylation of lipooligosaccharide (LOS) in the presence of capsule did not influence rates of attachment or invasion for strain MC58. However, removal of either serogroup B capsule or LOS sialylation in strain NMB-CDC increased bacterial attachment to host cells to the same extent. Although the level of inhibition of attachment by capsule was different between these strains, the regulation of the capsule synthesis locus by the two-component response regulator MisR, and the level of surface capsule determined by flow cytometry were not significantly different. However, the diplococci of strain NMB-CDC were shown to have a 1.89-fold greater surface area than strain MC58 by flow cytometry. It was proposed that the increase in surface area without changing the amount of anchored glycolipid capsule in the outer membrane would result in a sparser capsule and increase surface hydrophobicity. Strain NMB-CDC was shown to be more hydrophobic than strain MC58 using hydrophobicity interaction chromatography and microbial adhesion-to-solvents assays. In conclusion, improved levels of adherence of strain NMB-CDC to cell lines was associated with increased bacterial cell surface and surface hydrophobicity. This study shows that there is diversity in bacterial cell surface area and surface hydrophobicity within N. meningitidis which influence steps in meningococcal pathogenesis.

Homologous recombination drives both sequence diversity and gene content variation in Neisseria meningitidis

The study of genetic and phenotypic variation is fundamental for understanding the dynamics of bacterial genome evolution and untangling the evolution and epidemiology of bacterial pathogens. Neisseria meningitidis (Nm) is among the most intriguing bacterial pathogens in genomic studies due to its dynamic population structure and complex forms of pathogenicity. Extensive genomic variation within identical clonal complexes (CCs) in Nm has been recently reported and suggested to be the result of homologous recombination, but the extent to which recombination contributes to genomic variation within identical CCs has remained unclear. In this study, we sequenced two Nm strains of identical serogroup (C) and multi-locus sequence type (ST60), and conducted a systematic analysis with an additional 34 Nm genomes. Our results revealed that all gene content variation between the two ST60 genomes was introduced by homologous recombination at the conserved flanking genes, and 94.25% or more of sequence divergence was caused by homologous recombination. Recombination was found in genes associated with virulence factors, antigenic outer membrane proteins, and vaccine targets, suggesting an important role of homologous recombination in rapidly altering the pathogenicity and antigenicity of Nm. Recombination was also evident in genes of the restriction and modification systems, which may undermine barriers to DNA exchange. In conclusion, homologous recombination can drive both gene content variation and sequence divergence in Nm. These findings shed new light on the understanding of the rapid pathoadaptive evolution of Nm and other recombinogenic bacterial pathogens.

Developing insights into the mechanisms of evolution of bacterial pathogens from whole-genome sequences

Evolution of bacterial pathogen populations has been detected in a variety of ways including phenotypic tests, such as metabolic activity, reaction to antisera and drug resistance and genotypic tests that measure variation in chromosome structure, repetitive loci and individual gene sequences. While informative, these methods only capture a small subset of the total variation and, therefore, have limited resolution. Advances in sequencing technologies have made it feasible to capture whole-genome sequence variation for each sample under study, providing the potential to detect all changes at all positions in the genome from single nucleotide changes to large-scale insertions and deletions. In this review, we focus on recent work that has applied this powerful new approach and summarize some of the advances that this has brought in our understanding of the details of how bacterial pathogens evolve.

A neutral theory of genome evolution and the frequency distribution of genes

Background

The gene composition of bacteria of the same species can differ significantly between isolates. Variability in gene composition can be summarized in terms of gene frequency distributions, in which individual genes are ranked according to the frequency of genomes in which they appear. Empirical gene frequency distributions possess a U-shape, such that there are many rare genes, some genes of intermediate occurrence, and many common genes. It would seem that U-shaped gene frequency distributions can be used to infer the essentiality and/or importance of a gene to a species. Here, we ask: can U-shaped gene frequency distributions, instead, arise generically via neutral processes of genome evolution?

Results

We introduce a neutral model of genome evolution which combines birth-death processes at the organismal level with gene uptake and loss at the genomic level. This model predicts that gene frequency distributions possess a characteristic U-shape even in the absence of selective forces driving genome and population structure. We compare the model predictions to empirical gene frequency distributions from 6 multiply sequenced species of bacterial pathogens. We fit the model with constant population size to data, matching U-shape distributions albeit without matching all quantitative features of the distribution. We find stronger model fits in the case where we consider exponentially growing populations. We also show that two alternative models which contain a "rigid" and "flexible" core component of genomes provide strong fits to gene frequency distributions.

Conclusions

The analysis of neutral models of genome evolution suggests that U-shaped gene frequency distributions provide less information than previously suggested regarding gene essentiality. We discuss the need for additional theory and genomic level information to disentangle the roles of evolutionary mechanisms operating within and amongst individuals in driving the dynamics of gene distributions.

European efforts to harmonize typing of meningococci

Neisseria meningitidis causes life-threatening disease in infants, toddlers, and adolescents. Besides representative case notification, public health management of the disease requires bacterial typing information. European reference laboratories and state epidemiologists in collaboration with European institutions have driven forward the harmonization of typing by rigorously adopting DNA sequence typing and using common reference databases. External quality assessment has been provided by supranational networks, i.e. EU-IBIS and IBD-Labnet. The recent development of novel protein-based vaccines targeting serogroup B strains highlights the necessity to complement standard typing schemes by specific vaccine antigen typing including antigen expression analysis. Although not yet feasible for routine application on hundreds of strains, novel database structures have been developed to accommodate deep sequencing data.

Transcriptome analysis of Neisseria meningitidis in human whole blood and mutagenesis studies identify virulence factors involved in blood survival

During infection Neisseria meningitidis (Nm) encounters multiple environments within the host, which makes rapid adaptation a crucial factor for meningococcal survival. Despite the importance of invasion into the bloodstream in the meningococcal disease process, little is known about how Nm adapts to permit survival and growth in blood. To address this, we performed a time-course transcriptome analysis using an ex vivo model of human whole blood infection. We observed that Nm alters the expression of ≈30% of ORFs of the genome and major dynamic changes were observed in the expression of transcriptional regulators, transport and binding proteins, energy metabolism, and surface-exposed virulence factors. In particular, we found that the gene encoding the regulator Fur, as well as all genes encoding iron uptake systems, were significantly up-regulated. Analysis of regulated genes encoding for surface-exposed proteins involved in Nm pathogenesis allowed us to better understand mechanisms used to circumvent host defenses. During blood infection, Nm activates genes encoding for the factor H binding proteins, fHbp and NspA, genes encoding for detoxifying enzymes such as SodC, Kat and AniA, as well as several less characterized surface-exposed proteins that might have a role in blood survival. Through mutagenesis studies of a subset of up-regulated genes we were able to identify new proteins important for survival in human blood and also to identify additional roles of previously known virulence factors in aiding survival in blood. Nm mutant strains lacking the genes encoding the hypothetical protein NMB1483 and the surface-exposed proteins NalP, Mip and NspA, the Fur regulator, the transferrin binding protein TbpB, and the L-lactate permease LctP were sensitive to killing by human blood. This increased knowledge of how Nm responds to adaptation in blood could also be helpful to develop diagnostic and therapeutic strategies to control the devastating disease cause by this microorganism.

Virulence evolution of the human pathogen Neisseria meningitidis by recombination in the core and accessory genome

Background

Neisseria meningitidis is a naturally transformable, facultative pathogen colonizing the human nasopharynx. Here, we analyze on a genome-wide level the impact of recombination on gene-complement diversity and virulence evolution in N. meningitidis. We combined comparative genome hybridization using microarrays (mCGH) and multilocus sequence typing (MLST) of 29 meningococcal isolates with computational comparison of a subset of seven meningococcal genome sequences.

Principal Findings

We found that lateral gene transfer of minimal mobile elements as well as prophages are major forces shaping meningococcal population structure. Extensive gene content comparison revealed novel associations of virulence with genetic elements besides the recently discovered meningococcal disease associated (MDA) island. In particular, we identified an association of virulence with a recently described canonical genomic island termed IHT-E and a differential distribution of genes encoding RTX toxin- and two-partner secretion systems among hyperinvasive and non-hyperinvasive lineages. By computationally screening also the core genome for signs of recombination, we provided evidence that about 40% of the meningococcal core genes are affected by recombination primarily within metabolic genes as well as genes involved in DNA replication and repair. By comparison with the results of previous mCGH studies, our data indicated that genetic structuring as revealed by mCGH is stable over time and highly similar for isolates from different geographic origins.

Conclusions

Recombination comprising lateral transfer of entire genes as well as homologous intragenic recombination has a profound impact on meningococcal population structure and genome composition. Our data support the hypothesis that meningococcal virulence is polygenic in nature and that differences in metabolism might contribute to virulence.

Characterization of Neisseria meningitidis isolates that do not express the virulence factor and vaccine antigen factor H binding protein

Neisseria meningitidis remains a leading cause of bacterial sepsis and meningitis. Complement is a key component of natural immunity against this important human pathogen, which has evolved multiple mechanisms to evade complement-mediated lysis. One approach adopted by the meningococcus is to recruit a human negative regulator of the complement system, factor H (fH), to its surface via a lipoprotein, factor H binding protein (fHbp). Additionally, fHbp is a key antigen in vaccines currently being evaluated in clinical trials. Here we characterize strains of N. meningitidis from several distinct clonal complexes which do not express fHbp; all strains were recovered from patients with disseminated meningococcal disease. We demonstrate that these strains have either a frameshift mutation in the fHbp open reading frame or have entirely lost fHbp and some flanking sequences. No fH binding was detected to other ligands among the fHbp-negative strains. The implications of these findings for meningococcal pathogenesis and prevention are discussed.

Neisseria meningitidis is structured in clades associated with restriction modification systems that modulate homologous recombination

Molecular data on a limited number of chromosomal loci have shown that the population of Neisseria meningitidis (Nm), a deadly human pathogen, is structured in distinct lineages. Given that the Nm population undergoes substantial recombination, the mechanisms resulting in the evolution of these lineages, their persistence in time, and the implications for the pathogenicity of the bacterium are not yet completely understood. Based on whole-genome sequencing, we show that Nm is structured in phylogenetic clades. Through acquisition of specific genes and through insertions and rearrangements, each clade has acquired and remodeled specific genomic tracts, with the potential to impact on the commensal and virulence behavior of Nm. Despite this clear evidence of a structured population, we confirm high rates of detectable recombination throughout the whole Nm chromosome. However, gene conversion events were found to be longer within clades than between clades, suggesting a DNA cleavage mechanism associated with the phylogeny of the species. We identify 22 restriction modification systems, probably acquired by horizontal gene transfer from outside of the species/genus, whose distribution in the different strains coincides with the phylogenetic clade structure. We provide evidence that these clade-associated restriction modification systems generate a differential barrier to DNA exchange consistent with the observed population structure. These findings have general implications for the emergence of lineage structure and virulence in recombining bacterial populations, and they could provide an evolutionary framework for the population biology of a number of other bacterial species that show contradictory population structure and dynamics.

Whole-genome sequence of the transformable Neisseria meningitidis serogroup A strain WUE2594

Serogroup A meningococci are a leading cause of bacterial meningitis in children and young adults worldwide. However, the genetic basis of serogroup A strains' virulence and their epidemiological properties remain poorly understood. Therefore, we sequenced the complete genome of the transformable Neisseria meningitidis serogroup A strain WUE2594.

Genomic fluidity: an integrative view of gene diversity within microbial populations

Background

The dual concepts of pan and core genomes have been widely adopted as means to assess the distribution of gene families within microbial species and genera. The core genome is the set of genes shared by a group of organisms; the pan genome is the set of all genes seen in any of these organisms. A variety of methods have provided drastically different estimates of the sizes of pan and core genomes from sequenced representatives of the same groups of bacteria.

Results

We use a combination of mathematical, statistical and computational methods to show that current predictions of pan and core genome sizes may have no correspondence to true values. Pan and core genome size estimates are problematic because they depend on the estimation of the occurrence of rare genes and genomes, respectively, which are difficult to estimate precisely because they are rare. Instead, we introduce and evaluate a robust metric - genomic fluidity - to categorize the gene-level similarity among groups of sequenced isolates. Genomic fluidity is a measure of the dissimilarity of genomes evaluated at the gene level.

Conclusions

The genomic fluidity of a population can be estimated accurately given a small number of sequenced genomes. Further, the genomic fluidity of groups of organisms can be compared robustly despite variation in algorithms used to identify genes and their homologs. As such, we recommend that genomic fluidity be used in place of pan and core genome size estimates when assessing gene diversity within genomes of a species or a group of closely related organisms.

Independent evolution of the core and accessory gene sets in the genus Neisseria: insights gained from the genome of Neisseria lactamica isolate 020-06

BACKGROUND

The genus Neisseria contains two important yet very different pathogens, N. meningitidis and N. gonorrhoeae, in addition to non-pathogenic species, of which N. lactamica is the best characterized. Genomic comparisons of these three bacteria will provide insights into the mechanisms and evolution of pathogenesis in this group of organisms, which are applicable to understanding these processes more generally.

RESULTS

Non-pathogenic N. lactamica exhibits very similar population structure and levels of diversity to the meningococcus, whilst gonococci are essentially recent descendents of a single clone. All three species share a common core gene set estimated to comprise around 1190 CDSs, corresponding to about 60% of the genome. However, some of the nucleotide sequence diversity within this core genome is particular to each group, indicating that cross-species recombination is rare in this shared core gene set. Other than the meningococcal cps region, which encodes the polysaccharide capsule, relatively few members of the large accessory gene pool are exclusive to one species group, and cross-species recombination within this accessory genome is frequent.

CONCLUSION

The three Neisseria species groups represent coherent biological and genetic groupings which appear to be maintained by low rates of inter-species horizontal genetic exchange within the core genome. There is extensive evidence for exchange among positively selected genes and the accessory genome and some evidence of hitch-hiking of housekeeping genes with other loci. It is not possible to define a 'pathogenome' for this group of organisms and the disease causing phenotypes are therefore likely to be complex, polygenic, and different among the various disease-associated phenotypes observed.

Meningococcal internalization into human endothelial and epithelial cells is triggered by the influx of extracellular L-glutamate via GltT L-glutamate ABC transporter in Neisseria meningitidis

Meningococcal internalization into human cells is likely to be a consequence of meningococcal adhesion to human epithelial and endothelial cells. Here, we identified three transposon mutants of Neisseria meningitidis that were primarily defective in the internalization of human brain microvascular endothelial cells (HBMEC), with insertions occurring in the gltT (a sodium-independent L-glutamate transporter) gene or its neighboring gene, NMB1964 (unknown function). NMB1964 was tentatively named gltM in this study because of the presence of a mammalian cell entry (MCE)-related domain in the deduced amino acid sequences. The null ΔgltT-ΔgltM N. meningitidis mutant was also defective in the internalization into human umbilical vein endothelial cells and the human lung carcinoma epithelial cell line A549, and the defect was suppressed by transcomplementation of the mutants with gltT(+)-gltM(+) genes. The intracellular survival of the ΔgltT-ΔgltM mutant in HBMEC was not largely different from that of the wild-type strain under our experimental conditions. Introduction of a1-bp deletion and amber or ochre mutations in gltT-gltM genes resulted in the loss of efficient internalization into HBMEC. The defect in meningococcal internalization into HBMEC and L-glutamate uptake in the ΔgltT-ΔgltM mutant were suppressed only in strains expressing both GltT and GltM proteins. The efficiency of meningococcal invasion to HBMEC decreased under L-glutamate-depleted conditions. Furthermore, ezrin, a key membrane-cytoskeleton linker, accumulated beneath colonies of the gltT(+)-gltM(+) N. meningitidis strain but not of the ΔgltT-ΔgltM mutant. These findings suggest that l-glutamate influx via the GltT-GltM L-glutamate ABC transporter serves as a cue for N. meningitidis internalization into host cells.

Comparative genome biology of a serogroup B carriage and disease strain supports a polygenic nature of meningococcal virulence

Neisseria meningitidis serogroup B strains are responsible for most meningococcal cases in the industrialized countries, and strains belonging to the clonal complex ST-41/44 are among the most prevalent serogroup B strains in carriage and disease. Here, we report the first genome and transcriptome comparison of a serogroup B carriage strain from the clonal complex ST-41/44 to the serogroup B disease strain MC58 from the clonal complex ST-32. Both genomes are highly colinear, with only three major genome rearrangements that are associated with the integration of mobile genetic elements. They further differ in about 10% of their gene content, with the highest variability in gene presence as well as gene sequence found for proteins involved in host cell interactions, including Opc, NadA, TonB-dependent receptors, RTX toxin, and two-partner secretion system proteins. Whereas housekeeping genes coding for metabolic functions were highly conserved, there were considerable differences in their expression pattern upon adhesion to human nasopharyngeal cells between both strains, including differences in energy metabolism and stress response. In line with these genomic and transcriptomic differences, both strains also showed marked differences in their in vitro infectivity and in serum resistance. Taken together, these data support the concept of a polygenic nature of meningococcal virulence comprising differences in the repertoire of adhesins as well as in the regulation of metabolic genes and suggest a prominent role for immune selection and genetic drift in shaping the meningococcal genome.

Biofilm formation by the human pathogen Neisseria meningitidis

The past decade has seen an increasing interest in biofilm formation by Neisseria meningitidis, a human facultative pathogen causing life-threatening childhood disease commencing from asymptomatic nasopharyngeal colonization. Studying the biology of in vitro biofilm formation improves the understanding of inter-bacterial processes in asymptomatic carriage, of bacterial aggregate formation on host cells, and of meningococcal population biology. This paper reviews publications referring to meningococcal biofilm formation with an emphasis on the role of motility and of extracellular DNA. The theory of sub-dividing the meningococcal population in settler and spreader lineages is discussed, which provides a mechanistic framework for the assumed balance of colonization efficacy and transmission frequency.

NeMeSys: a biological resource for narrowing the gap between sequence and function in the human pathogen Neisseria meningitidis

BACKGROUND

Genome sequences, now available for most pathogens, hold promise for the rational design of new therapies. However, biological resources for genome-scale identification of gene function (notably genes involved in pathogenesis) and/or genes essential for cell viability, which are necessary to achieve this goal, are often sorely lacking. This holds true for Neisseria meningitidis, one of the most feared human bacterial pathogens that causes meningitis and septicemia.

RESULTS

By determining and manually annotating the complete genome sequence of a serogroup C clinical isolate of N. meningitidis (strain 8013) and assembling a library of defined mutants in up to 60% of its non-essential genes, we have created NeMeSys, a biological resource for Neisseria meningitidis systematic functional analysis. To further enhance the versatility of this toolbox, we have manually (re)annotated eight publicly available Neisseria genome sequences and stored all these data in a publicly accessible online database. The potential of NeMeSys for narrowing the gap between sequence and function is illustrated in several ways, notably by performing a functional genomics analysis of the biogenesis of type IV pili, one of the most widespread virulence factors in bacteria, and by identifying through comparative genomics a complete biochemical pathway (for sulfur metabolism) that may potentially be important for nasopharyngeal colonization.

CONCLUSIONS

By improving our capacity to understand gene function in an important human pathogen, NeMeSys is expected to contribute to the ongoing efforts aimed at understanding a prokaryotic cell comprehensively and eventually to the design of new therapies.