Extensive genomic variation within clonal complexes of Neisseria meningitidis

Models in parasite and pathogen evolution: Genomic analysis reveals predominant clonality and progressive evolution at all evolutionary scales in parasitic protozoa, yeasts and bacteria

The predominant clonal evolution (PCE) model of pathogenic microorganisms postulates that the impact of genetic recombination in those pathogens' natural populations is not enough to erase a persistent phylogenetic signal at all evolutionary scales from microevolution till geological times in the whole ecogeographical range of the species considered. We have tested this model with a set of representative parasitic protozoa, yeasts and bacteria in the light of the most recent genomic data. All surveyed species, including those that were considered as highly recombining, exhibit similar PCE patterns above and under the species level, from macro- to micro-evolutionary scales (Russian doll pattern), suggesting gradual evolution. To our knowledge, it is the first time that such a strong common evolutionary feature among very diverse pathogens has been evidenced. The implications of this model for basic biology and applied research are exposed. These implications include our knowledge on the pathogens' reproductive mode, their population structure, the possibility to type strain and to follow up epidemics (molecular epidemiology) and to revisit pathogens' taxonomy through a flexible use of the phylogenetic species concept (Cracraft, 1983).

Rapid identification, capsular typing and molecular characterization of Streptococcus pneumoniae by using whole genome nanopore sequencing

BACKGROUND

Whole genome sequencing has emerged as a useful tool for identification and molecular characterization of pathogens. MinION (Oxford Nanopore) is a real-time third generation sequencer whose portability, affordability and speed in data production make of it an attractive device for whole genome sequencing. The objective of this study is to evaluate MinION sequencer for pathogen identification and molecular characterization of Streptococcus pneumoniae isolated at a children's Hospital. Whole genome sequencing of 32 Streptococcus pneumoniae invasive isolates, previously characterized by standard methods (Quellung reaction, Multiplex PCR and Sanger-MLST), were performed. DNA was extracted using ZymoBIOMICS DNA Microprep kit. Quantification and purity of DNA was assessed by Qubit and Nanodrop, respectively. Library preparation was performed using the Rapid Barcoding Kit. Real-time workflow EPI2ME platform "What's it in my pot" was used for species identification. Fast5 sequences were converted into FASTQ by Albacore software. Reads were assembled using CANU software. PathogenWatch, genomic epidemiology and pubmlst online tools were used for capsular typing and/or whole genome-MLST profile.

RESULTS

Rapid identification of Streptococcus pneumoniae was achieved by "What's in my pot". Capsular typing was correctly assigned with PathogenWatch in all 32 isolates at serogroup level and 24 at serotype level. Whole genome-MLST results obtained by genomic epidemiology and pubmlst were consistent with double locus variant clonal complex obtained by Sanger-MLST in 31 isolates.

CONCLUSION

MinION sequencer provides a rapid, cost-effective and promising pathway for performing WGS by a pocked-sized device for epidemiological purposes but improving its sequencing accuracy will make it more appealing to be used in clinical microbiology laboratories.

Complete genome and methylome analysis of Neisseria meningitidis associated with increased serogroup Y disease

Invasive meningococcal disease (IMD) due to serogroup Y Neisseria meningitidis emerged in Europe during the 2000s. Draft genomes of serogroup Y isolates in Sweden revealed that although the population structure of these isolates was similar to other serogroup Y isolates internationally, a distinct strain (YI) and more specifically a sublineage (1) of this strain was responsible for the increase of serogroup Y IMD in Sweden. We performed single molecule real-time (SMRT) sequencing on eight serogroup Y isolates from different sublineages to unravel the genetic and epigenetic factors delineating them, in order to understand the serogroup Y emergence. Extensive comparisons between the serogroup Y sublineages of all coding sequences, complex genomic regions, intergenic regions, and methylation motifs revealed small point mutations in genes mainly encoding hypothetical and metabolic proteins, and non-synonymous variants in genes involved in adhesion, iron acquisition, and endotoxin production. The methylation motif CACNNNNNTAC was only found in isolates of sublineage 2. Only seven genes were putatively differentially expressed, and another two genes encoding hypothetical proteins were only present in sublineage 2. These data suggest that the serogroup Y IMD increase in Sweden was most probably due to small changes in genes important for colonization and transmission.

Interconnected clusters of invasive meningococcal disease due to Neisseria meningitidis serogroup C ST-11 (cc11), involving bisexuals and men who have sex with men, with discos and gay-venues hotspots of transmission, Tuscany, Italy, 2015 to 2016

In 2015 an increased incidence of invasive meningococcal disease due to serogroup-C (MenC) occurred in Tuscany, Italy. This led the Regional Health Authority of Tuscany to implement a reactive immunisation campaign and to launch an epidemiological field investigation aiming to address targeted immunisation interventions. In 2011–14, 10 MenC cases had been reported compared with 62 cases in 2015–16. The case fatality rate was 21% (n = 13) and 51 cases (82.3%) were confirmed as C:P1.5–1,10–8:F3–6:ST-11(cc11). Overall, 17 clusters were recognised. Six discos and four gay-venues were found to have a role as transmission-hotspots, having been attended by 20 and 14 cases in the 10 days before symptoms onset. Ten and three cases occurred, respectively, among men who have sex with men (MSM) and bisexual individuals, who were involved in 11 clusters. In addition, heterosexual cases (n = 5) attending gay-venues were also found. Secondary cases were not identified. Molecular typing indicated close relationship with MenC clusters recently described among gay, bisexual and other MSM in Europe and the United States, suggesting a possible international spread of the serogroup-C-variant P1.5–1,10–8:F3–6:ST-11(cc11) in this population-group; however, epidemiological links were not identified. In December 2016, a targeted vaccination campaign involving discos and lesbian, gay, bisexual, and transgender (LGBT) associations was implemented. During 2017, 10 cases of MenC occurred, compared with 32 and 30 cases reported in 2015 and 2016 respectively, suggesting the effectiveness of the reactive and targeted immunisation programmes.

The comparative population genetics of Neisseria meningitidis and Neisseria gonorrhoeae

Neisseria meningitidis and N. gonorrhoeae are closely related pathogenic bacteria. To compare their population genetics, we compiled a dataset of 1,145 genes found across 20 N. meningitidis and 15 N. gonorrhoeae genomes. We find that N. meningitidis is seven-times more diverse than N. gonorrhoeae in their combined core genome. Both species have acquired the majority of their diversity by recombination with divergent strains, however, we find that N. meningitidis has acquired more of its diversity by recombination than N. gonorrhoeae. We find that linkage disequilibrium (LD) declines rapidly across the genomes of both species. Several observations suggest that N. meningitidis has a higher effective population size than N. gonorrhoeae; it is more diverse, the ratio of non-synonymous to synonymous polymorphism is lower, and LD declines more rapidly to a lower asymptote in N. meningitidis. The two species share a modest amount of variation, half of which seems to have been acquired by lateral gene transfer and half from their common ancestor. We investigate whether diversity varies across the genome of each species and find that it does. Much of this variation is due to different levels of lateral gene transfer. However, we also find some evidence that the effective population size varies across the genome. We test for adaptive evolution in the core genome using a McDonald–Kreitman test and by considering the diversity around non-synonymous sites that are fixed for different alleles in the two species. We find some evidence for adaptive evolution using both approaches.

Neisseria genomics: current status and future perspectives

High-throughput whole genome sequencing has unlocked a multitude of possibilities enabling members of the Neisseria genus to be examined with unprecedented detail, including the human pathogens Neisseria meningitidis and Neisseria gonorrhoeae. To maximise the potential benefit of this for public health, it is becoming increasingly important to ensure that this plethora of data are adequately stored, disseminated and made readily accessible. Investigations facilitating cross-species comparisons as well as the analysis of global datasets will allow differences among and within species and across geographic locations and different times to be identified, improving our understanding of the distinct phenotypes observed. Recent advances in high-throughput platforms that measure the transcriptome, proteome and/or epigenome are also becoming increasingly employed to explore the complexities of Neisseria biology. An integrated approach to the analysis of these is essential to fully understand the impact these may have in the Neisseria genus. This article reviews the current status of some of the tools available for next generation sequence analysis at the dawn of the ‘post-genomic’ era.

Unraveling Neisseria meningitidis pathogenesis: from functional genomics to experimental models

Neisseria meningitidis is a harmless commensal bacterium finely adapted to humans. Unfortunately, under “privileged” conditions, it adopts a “devious” lifestyle leading to uncontrolled behavior characterized by the unleashing of molecular weapons causing potentially lethal disease such as sepsis and acute meningitis. Indeed, despite the lack of a classic repertoire of virulence genes in N. meningitidis separating commensal from invasive strains, molecular epidemiology and functional genomics studies suggest that carriage and invasive strains belong to genetically distinct populations characterized by an exclusive pathogenic potential. In the last few years, “omics” technologies have helped scientists to unwrap the framework drawn by N. meningitidis during different stages of colonization and disease. However, this scenario is still incomplete and would benefit from the implementation of physiological tissue models for the reproduction of mucosal and systemic interactions in vitro. These emerging technologies supported by recent advances in the world of stem cell biology hold the promise for a further understanding of N. meningitidis pathogenesis.

Chromosomal rearrangements and protein globularity changes in Mycobacterium tuberculosis isolates from cerebrospinal fluid

Background

Meningitis is a major cause of mortality in tuberculosis (TB). It is not clear what factors promote central nervous system invasion and pathology but it has been reported that certain strains of Mycobacterium tuberculosis (Mtb) might have genetic traits associated with neurotropism.

Methods

In this study, we generated whole genome sequences of eight clinical strains of Mtb that were isolated from the cerebrospinal fluid (CSF) of patients presenting with tuberculous meningitis (TBM) in Malaysia, and compared them to the genomes of H37Rv and other respiratory Mtb genomes either downloaded from public databases or extracted from local sputum isolates. We aimed to find genomic features that might be distinctly different between CSF-derived and respiratory Mtb.

Results

Genome-wide comparisons revealed rearrangements (translocations, inversions, insertions and deletions) and non-synonymous SNPs in our CSF-derived strains that were not observed in the respiratory Mtb genomes used for comparison. These rearranged segments were rich in genes for PE (proline-glutamate)/PPE (proline-proline-glutamate), transcriptional and membrane proteins. Similarly, most of the ns SNPs common in CSF strains were noted in genes encoding PE/PPE proteins. Protein globularity differences were observed among mycobacteria from CSF and respiratory sources and in proteins previously reported to be associated with TB meningitis. Transcription factors and other transcription regulators featured prominently in these proteins. Homologs of proteins associated with Streptococcus pneumoniae meningitis and Neisseria meningitidis virulence were identified in neuropathogenic as well as respiratory mycobacterial spp. examined in this study.

Discussion

The occurrence of in silico genetic differences in CSF-derived but not respiratory Mtb suggests their possible involvement in the pathogenesis of TBM. However, overall findings in this comparative analysis support the postulation that TB meningeal infection is more likely to be related to the expression of multiple virulence factors on interaction with host defences than to CNS tropism associated with specific genetic traits.

Genomic Investigation Reveals Highly Conserved, Mosaic, Recombination Events Associated with Capsular Switching among Invasive Neisseria meningitidis Serogroup W Sequence Type (ST)-11 Strains

Neisseria meningitidis is an important cause of meningococcal disease globally. Sequence type (ST)-11 clonal complex (cc11) is a hypervirulent meningococcal lineage historically associated with serogroup C capsule and is believed to have acquired the W capsule through a C to W capsular switching event. We studied the sequence of capsule gene cluster (cps) and adjoining genomic regions of 524 invasive W cc11 strains isolated globally. We identified recombination breakpoints corresponding to two distinct recombination events within W cc11: A 8.4-kb recombinant region likely acquired from W cc22 including the sialic acid/glycosyl-transferase gene, csw resulted in a C→W change in capsular phenotype and a 13.7-kb recombinant segment likely acquired from Y cc23 lineage includes 4.5 kb of cps genes and 8.2 kb downstream of the cps cluster resulting in allelic changes in capsule translocation genes. A vast majority of W cc11 strains (497/524, 94.8%) retain both recombination events as evidenced by sharing identical or very closely related capsular allelic profiles. These data suggest that the W cc11 capsular switch involved two separate recombination events and that current global W cc11 meningococcal disease is caused by strains bearing this mosaic capsular switch.

Distribution of the type III DNA methyltransferases modA, modB and modD among Neisseria meningitidis genotypes: implications for gene regulation and virulence

Neisseria meningitidis is a human-specific bacterium that varies in invasive potential. All meningococci are carried in the nasopharynx, and most genotypes are very infrequently associated with invasive meningococcal disease; however, those belonging to the ‘hyperinvasive lineages’ are more frequently associated with sepsis or meningitis. Genome content is highly conserved between carriage and disease isolates, and differential gene expression has been proposed as a major determinant of the hyperinvasive phenotype. Three phase variable DNA methyltransferases (ModA, ModB and ModD), which mediate epigenetic regulation of distinct phase variable regulons (phasevarions), have been identified in N. meningitidis. Each mod gene has distinct alleles, defined by their Mod DNA recognition domain, and these target and methylate different DNA sequences, thereby regulating distinct gene sets. Here 211 meningococcal carriage and >1,400 disease isolates were surveyed for the distribution of meningococcal mod alleles. While modA11-12 and modB1-2 were found in most isolates, rarer alleles (e.g., modA15, modB4, modD1-6) were specific to particular genotypes as defined by clonal complex. This suggests that phase variable Mod proteins may be associated with distinct phenotypes and hence invasive potential of N. meningitidis strains.

Global epidemiology of capsular group W meningococcal disease (1970-2015): Multifocal emergence and persistence of hypervirulent sequence type (ST)-11 clonal complex

Following an outbreak in Mecca Saudi Arabia in 2000, meningococcal strains expressing capsular group W (W) emerged as a major cause of invasive meningococcal disease (IMD) worldwide. The Saudi Arabian outbreak strain (Hajj clone) belonging to the ST-11 clonal complex (cc11) is similar to W cc11 causing occasional sporadic disease before 2000. Since 2000, W cc11 has caused large meningococcal disease epidemics in the African meningitis belt and endemic disease in South America, Europe and China. Traditional molecular epidemiologic typing suggested that a majority of current W cc11 burden represented global spread of the Hajj clone. However, recent whole genome sequencing (WGS) analyses revealed significant genetic heterogeneity among global W cc11 strains. While continued spread of the Hajj clone occurs in the Middle East, the meningitis belt and South Africa have co-circulation of the Hajj clone and other unrelated W cc11 strains. Notably, South America, the UK, and France share a genetically distinct W cc11 strain. Other W lineages persist in low numbers in Europe, North America and the meningitis belt. In summary, WGS is helping to unravel the complex genomic epidemiology of group W meningococcal strains. Wider application of WGS and strengthening of global IMD surveillance is necessary to monitor the continued evolution of group W lineages.

Genomic epidemiology of age-associated meningococcal lineages in national surveillance: an observational cohort study

Background

Invasive meningococcal disease (IMD) is a worldwide health issue that is potentially preventable with vaccination. In view of its sporadic nature and the high diversity of Neisseria meningitidis, epidemiological surveillance incorporating detailed isolate characterisation is crucial for effective control and understanding the evolving epidemiology of IMD. The Meningitis Research Foundation Meningococcus Genome Library (MRF-MGL) exploits whole-genome sequencing (WGS) for this purpose and presents data on a comprehensive and coherent IMD isolate collection from England and Wales via the internet. We assessed the contribution of these data to investigating IMD epidemiology.

Methods

WGS data were obtained for all 899 IMD isolates available for England and Wales in epidemiological years 2010–11 and 2011–12. The data had been annotated at 1720 loci, analysed, and disseminated online. Information was also available on meningococcal population structure and vaccine (Bexsero, GlaxoSmithKline, Brentford, Middlesex, UK) antigen variants, which enabled the investigation of IMD-associated genotypes over time and by patients' age groups. Population genomic analyses were done with a hierarchical gene-by-gene approach.

Findings

The methods used by MRF-MGL efficiently characterised IMD isolates and information was provided in plain language. At least 20 meningococcal lineages were identified, three of which (hyperinvasive clonal complexes 41/44 [lineage 3], 269 [lineage 2], and 23 [lineage 23]) were responsible for 528 (59%) of IMD isolates. Lineages were highly diverse and showed evidence of extensive recombination. Specific lineages were associated with IMD in particular age groups, with notable diversity in the youngest and oldest individuals. The increased incidence of IMD from 1984 to 2010 in England and Wales was due to successive and concurrent epidemics of different lineages. Genetically, 74% of isolates were characterised as encoding group B capsules: 16% group Y, 6% group W, and 3% group C. Exact peptide matches for individual Bexsero vaccine antigens were present in up to 26% of isolates.

Interpretation

The MRF-MGL represents an effective, broadly applicable model for the storage, analysis, and dissemination of WGS data that can facilitate real-time genomic pathogen surveillance. The data revealed information crucial to effective deployment and assessment of vaccines against N meningitidis.

Funding

Meningitis Research Foundation, Wellcome Trust, Public Health England, European Union.

Phylogenetic study of clonal complex (CC)198 capsule null locus (cnl) genomes: A distinctive group within the species Neisseria meningitidis

Capsule null locus (cnl) strains, one type of specific unencapsulated Neisseria spp., only have regions D and E of the capsule gene cluster which encodes the genes for capsule biosynthesis, modification, and transportation. Compared with encapsulated strains, regions A and C of cnl strains have been replaced by 113 or 114 bp conserved non-coding sequences. Cnl strains include multiple clonal complexes (CC). According to previous studies, CC198 is the major clonal lineage in both cnl patients and healthy cnl carriers. We hypothesized that CC198 possesses different genome characteristics compared with other cnl strains. In this study, we obtained the draft genomes of two CC198 strains from healthy carriers. Using 75071 single nucleotide polymorphisms located in 1163 core genes, we constructed the phylogenetic relationships between a batch of representative Neisseria meningitidis genomes. CC198 and CC1136 clustered together, but apart from other N. meningitidis strains including CC53. We also aligned the sequences of genes located in regions D and E of the capsule gene locus from encapsulated and unencapsulated strains. A number of possible recombination events were identified in the galE and tex genes between different serogroups of encapsulated N. meningitidis and CC53 strains, especially in tex. In contrast, there is almost no recombination in N. meningitidis CC198 strains. These results showed that CC198 belongs to a phylogenetically distinct group within the species N. meningitidis, which may be directly derived from the cnl-type ancestor of N. meningitidis. The encapsulated strains may acquire other necessary genes for capsule formation by horizontal transfer.

Genomic Analysis of the Evolution and Global Spread of Hyper-invasive Meningococcal Lineage 5

Background

The predominant model for bacterial pandemics is the emergence of a virulent variant that diversifies as it spreads in human populations. We investigated a 40-year meningococcal disease pandemic caused by the hyper-invasive ET-5/ST-32 complex.

Methods

A global collection of Neisseria meningitidis isolates dating from 1969 to 2008 was whole genome sequenced (WGS) and analysed using a gene-by-gene approach at http://pubmlst.org/neisseria.

Findings

Analysis of WGS data identified a ‘Lineage 5 pan genome’ of 1940 genes, 1752 (92%) of which were present in all isolates (Lineage 5 ‘core genome’). Genetic diversity, which was mostly generated by horizontal gene transfer, was unevenly distributed in the genome; however, genealogical analysis of diverse and conserved core genes, accessory genes, and antigen encoding genes, robustly identified a star phylogeny with a number of sub-lineages. Most European and American isolates belonged to one of two closely related sub-lineages, which had diversified before the identification of the pandemic in the 1970s. A third, genetically more diverse sub-lineage, was associated with Asian isolates. Several isolates had acquired DNA from the related gonococcus.

Interpretation

These data were inconsistent with a single point of origin followed by pandemic spread, rather suggesting that the sub-lineages had diversified and spread by asymptomatic transmission, with multiple distinct strains causing localised hyperendemic outbreaks.

•

Genomic analysis of hyper-invasive N. meningitidis lineage 5 distinguishing sub-lineages and localised outbreaks.

•

Identification of a gonococcal conjugative plasmid consistent with horizontal genetic transfer.

•

Discovery of a distinct type 4 secretion system previously unidentified in Neisseria.

•

Gene-by-gene analysis of conserved and diverse loci providing essential tools in serogroup B vaccine development.

A multi-component meningococcal serogroup B vaccine (4CMenB): the clinical development program

Recently approved in Europe and Australia, the multi-component meningococcal B vaccine, 4CMenB (Bexsero^®, Novartis Vaccines and Diagnostics), contains three surface-exposed recombinant proteins (fHbp, NadA, and NHBA) and New Zealand strain outer membrane vesicles (NZ OMV) with PorA 1.4 antigenicity. This comprehensive review of the 4CMenB clinical development program covers pivotal phase I/IIb/III studies in over 7,000 adults, adolescents, and infants. The immunological correlate for clinical protection used was human complement-mediated serum bactericidal activity titers ≥4 or 5 against indicator strains for individual antigens. Based on achievement of protective titers, a four-dose schedule (three primary doses and one booster dose) for infants and a two-dose schedule for adolescents provided the best results. Observed increases in injection site pain/tenderness and fever in infants, and injection site pain, malaise, and headache in adolescents compared with routine vaccines, were mostly mild to moderate; frequencies of rare events (Kawasaki disease, juvenile arthritis) were not significantly different from non-vaccinated individuals. 4CMenB is conservatively estimated to provide 66–91 % coverage against meningococcal serogroup B strains worldwide.

HGT turbulence: Confounding phylogenetic influence of duplicative horizontal transfer and differential gene conversion

Horizontal gene transfer (HGT) often leads to phylogenetic incongruence. When "duplicative HGT" introduces a second copy of a pre-existing gene, the two copies may then engage in gene conversion, leading to phylogenetically mosiac genes. When duplicative HGT is followed by differential gene conversion among descendant lineages, as under the DH-DC model, phylogenetic analysis is further complicated. To explore the effects of DH-DC on phylogeny reconstruction, we analyzed two sets of sequences: (1) an augmented set of plant mitochondrial atp1 sequences for which we recently published evidence of DH-DC; and (2) a set of simulated sequences for which we varied the extent of chimerism, the number of chimeric genes and nucleotide substitution rates. We show that the phylogenetic behavior of evolutionarily chimeric genes is highly volatile and depends on both the degree of chimerism and the number of differentially chimeric genes present in the analysis. Furthermore, we show that the presence of chimeric genes in gene trees can spuriously affect the phylogenetic position of purely native sequences, especially by attracting these sequences toward basal positions in trees. We propose the term "HGT turbulence" to describe these complex effects of evolutionarily chimeric genes on phylogenetic results.

Extensive genomic variation within clonal bacterial groups resulted from homologous recombination

Due to divergence, genetic variation is generally believed to be high among distantly related strains, low among closely related ones and little or none within the same classified clonal groups. Several recent genome-wide studies, however, revealed that significant genetic variation resides in a considerable number of genes among strains with identical MLST (Multilocus sequence typing) types and much of the variation was introduced by homologous recombination. Recognizing and understanding genomic variation within clonal bacterial groups could shed new light on the evolutionary path of infectious agents and the emergence of particularly pathogenic or virulent variants. This commentary presents our recent contributions to this line of work.

Metabolism and virulence in Neisseria meningitidis

A longstanding question in infection biology addresses the genetic basis for invasive behavior in commensal pathogens. A prime example for such a pathogen is Neisseria meningitidis. On the one hand it is a harmless commensal bacterium exquisitely adapted to humans, and on the other hand it sometimes behaves like a ferocious pathogen causing potentially lethal disease such as sepsis and acute bacterial meningitis. Despite the lack of a classical repertoire of virulence genes in N. meningitidis separating commensal from invasive strains, molecular epidemiology suggests that carriage and invasive strains belong to genetically distinct populations. In recent years, it has become increasingly clear that metabolic adaptation enables meningococci to exploit host resources, supporting the concept of nutritional virulence as a crucial determinant of invasive capability. Here, we discuss the contribution of core metabolic pathways in the context of colonization and invasion with special emphasis on results from genome-wide surveys. The metabolism of lactate, the oxidative stress response, and, in particular, glutathione metabolism as well as the denitrification pathway provide examples of how meningococcal metabolism is intimately linked to pathogenesis. We further discuss evidence from genome-wide approaches regarding potential metabolic differences between strains from hyperinvasive and carriage lineages and present new data assessing in vitro growth differences of strains from these two populations. We hypothesize that strains from carriage and hyperinvasive lineages differ in the expression of regulatory genes involved particularly in stress responses and amino acid metabolism under infection conditions.

Role of factor H binding protein in Neisseria meningitidis virulence and its potential as a vaccine candidate to broadly protect against meningococcal disease

Neisseria meningitidis is a Gram-negative microorganism that exists exclusively in humans and can cause devastating invasive disease. Although capsular polysaccharide-based vaccines against serogroups A, C, Y, and W135 are widely available, the pathway to a broadly protective vaccine against serogroup B has been more complex. The last 11 years has seen the discovery and development of the N. meningitidis serogroup B (MnB) outer membrane protein factor H binding protein (fHBP) as a vaccine component. Since the initial discovery of fHBP, a tremendous amount of work has accumulated on the diversity, structure, and regulation of this important protein. fHBP has proved to be a virulence factor for N. meningitidis and a target for functional bactericidal antibodies. fHBP is critical for survival of meningococci in the human host, as it is responsible for the primary interaction with human factor H (fH). Binding of hfH by the meningococcus serves to downregulate the host alternative complement pathway and helps the organism evade host innate immunity. Preclinical studies have shown that an fHBP-based vaccine can elicit serum bactericidal antibodies capable of killing MnB, and the vaccine has shown very encouraging results in human clinical trials. This report reviews our current knowledge of fHBP. In particular, we discuss the recent advances in our understanding of fHBP, its importance to N. meningitidis, and its potential role as a vaccine for preventing MnB disease.

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.

Controlling serogroup B invasive meningococcal disease: the Canadian perspective

With publically funded meningococcal immunization programs established in infants, children and adolescents, Canada is at the forefront of invasive meningococcal disease prevention. The advent of two new serogroup B vaccines that may protect against multiple disease-causing strains offers the potential to reduce endemic disease to very low levels in Canada. Canada likely will be one of the first countries with approval to use recombinant serogroup B vaccine. However, inclusion of these new vaccines into public immunization programs will be decided at the provincial/territorial level, rather than nationally, and may result initially in different immunization schedules throughout the country as we have seen with conjugate meningococcal vaccines. Such heterogeneous use and adoption of new vaccines complicates disease control, but may assist in evaluation of effectiveness. Minimally, it requires regionally specific information. In this article, the authors provide an overview of the Canadian epidemiology, serogroup B vaccine characteristics, potential strain coverage, immunization strategies and remaining postmarketing research questions.

Routine vaccination against MenB: considerations for implementation

Effective polysaccharide(conjugate) vaccines against Neisseria meningitidis serogroups A, C, W, and Y have been widely used, but serogroup B meningococci remain a major cause of severe invasive meningococcal disease (IMD) worldwide, especially in infants. Recently, a vaccine, 4CMenB (Bexsero^®), containing three recombinant proteins, and outer membrane vesicles (OMV) derived from a serogroup B meningococcal strain (MenB) has been licensed in Europe and Australia and is indicated for persons aged 2 mo or older. This article discusses what should be considered to enable a successful implementation of a broad coverage MenB vaccine in national immunization programs. Epidemiology data, vaccine characteristics including vaccine coverage, immunogenicity, post-implementation surveillance and costs are relevant aspects that should be taken into account when selecting an appropriate immunization strategy. The potential impact on strain variation and carriage, as well as monitoring vaccine effectiveness, and rare but potentially serious adverse events are points that need to be included in a post-implementation surveillance plan.

Constraints on genome dynamics revealed from gene distribution among the Ralstonia solanacearum species

Because it is suspected that gene content may partly explain host adaptation and ecology of pathogenic bacteria, it is important to study factors affecting genome composition and its evolution. While recent genomic advances have revealed extremely large pan-genomes for some bacterial species, it remains difficult to predict to what extent gene pool is accessible within or transferable between populations. As genomes bear imprints of the history of the organisms, gene distribution pattern analyses should provide insights into the forces and factors at play in the shaping and maintaining of bacterial genomes. In this study, we revisited the data obtained from a previous CGH microarrays analysis in order to assess the genomic plasticity of the R. solanacearum species complex. Gene distribution analyses demonstrated the remarkably dispersed genome of R. solanacearum with more than half of the genes being accessory. From the reconstruction of the ancestral genomes compositions, we were able to infer the number of gene gain and loss events along the phylogeny. Analyses of gene movement patterns reveal that factors associated with gene function, genomic localization and ecology delineate gene flow patterns. While the chromosome displayed lower rates of movement, the megaplasmid was clearly associated with hot-spots of gene gain and loss. Gene function was also confirmed to be an essential factor in gene gain and loss dynamics with significant differences in movement patterns between different COG categories. Finally, analyses of gene distribution highlighted possible highways of horizontal gene transfer. Due to sampling and design bias, we can only speculate on factors at play in this gene movement dynamic. Further studies examining precise conditions that favor gene transfer would provide invaluable insights in the fate of bacteria, species delineation and the emergence of successful pathogens.

Dialects of the DNA uptake sequence in Neisseriaceae

In all sexual organisms, adaptations exist that secure the safe reassortment of homologous alleles and prevent the intrusion of potentially hazardous alien DNA. Some bacteria engage in a simple form of sex known as transformation. In the human pathogen Neisseria meningitidis and in related bacterial species, transformation by exogenous DNA is regulated by the presence of a specific DNA Uptake Sequence (DUS), which is present in thousands of copies in the respective genomes. DUS affects transformation by limiting DNA uptake and recombination in favour of homologous DNA. The specific mechanisms of DUS-dependent genetic transformation have remained elusive. Bioinformatic analyses of family Neisseriaceae genomes reveal eight distinct variants of DUS. These variants are here termed DUS dialects, and their effect on interspecies commutation is demonstrated. Each of the DUS dialects is remarkably conserved within each species and is distributed consistent with a robust Neisseriaceae phylogeny based on core genome sequences. The impact of individual single nucleotide transversions in DUS on meningococcal transformation and on DNA binding and uptake is analysed. The results show that a DUS core 5'-CTG-3' is required for transformation and that transversions in this core reduce DNA uptake more than two orders of magnitude although the level of DNA binding remains less affected. Distinct DUS dialects are efficient barriers to interspecies recombination in N. meningitidis, N. elongata, Kingella denitrificans, and Eikenella corrodens, despite the presence of the core sequence. The degree of similarity between the DUS dialect of the recipient species and the donor DNA directly correlates with the level of transformation and DNA binding and uptake. Finally, DUS-dependent transformation is documented in the genera Eikenella and Kingella for the first time. The results presented here advance our understanding of the function and evolution of DUS and genetic transformation in bacteria, and define the phylogenetic relationships within the Neisseriaceae family.

Unrecognized fine-scale recombination can mimic the effects of adaptive radiation

Gene sequences can undergo accelerated nucleotide changes and rapid diversification. The rapid sequence changes can then potentially lead to phylogenetic incongruence. Recently, Bodilis et al. (2011) observed artificial phylogenetic incongruence using the Pseudomonas surface protein gene oprF, and hypothesized that it was the result of a long-branch attraction artifact ultimately caused by adaptive radiation. In this study, an alternative hypothesis, namely fine-scale recombination, was tested on the same dataset. The results reveal that regions in oprF are of different evolutionary origins, and the mosaic gene structure resulted in confounding phylogenetic signals. These findings demonstrate that unrecognized fine-scale recombination can confound the phylogenetic interpretation and emphasize the limitation of using whole genes as the unit of phylogenetic analysis.

Automated extraction of typing information for bacterial pathogens from whole genome sequence data: Neisseria meningitidis as an exemplar

Whole genome sequence (WGS) data are increasingly used to characterise bacterial pathogens. These data provide detailed information on the genotypes and likely phenotypes of aetiological agents, enabling the relationships of samples from potential disease outbreaks to be established precisely. However, the generation of increasing quantities of sequence data does not, in itself, resolve the problems that many microbiological typing methods have addressed over the last 100 years or so; indeed, providing large volumes of unstructured data can confuse rather than resolve these issues. Here we review the nascent field of storage of WGS data for clinical application and show how curated sequence-based typing schemes on websites have generated an infrastructure that can exploit WGS for bacterial typing efficiently. We review the tools that have been implemented within the PubMLST website to extract clinically useful, strain-characterisation information that can be provided to physicians and public health professionals in a timely, concise and understandable way. These data can be used to inform medical decisions such as how to treat a patient, whether to instigate public health action, and what action might be appropriate. The information is compatible both with previous sequence-based typing data and also with data obtained in the absence of WGS, providing a flexible infrastructure for WGS-based clinical microbiology.

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.

Unique secreted-surface protein complex of Lactobacillus rhamnosus, identified by phage display

Proteins are the most diverse structures on bacterial surfaces; hence, they are candidates for species- and strain-specific interactions of bacteria with the host, environment, and other microorganisms. Genomics has decoded thousands of bacterial surface and secreted proteins, yet the function of most cannot be predicted because of the enormous variability and a lack of experimental data that would allow deduction of function through homology. Here, we used phage display to identify a pair of interacting extracellular proteins in the probiotic bacterium Lactobacillus rhamnosus HN001. A secreted protein, SpcA, containing two bacterial immunoglobulin-like domains type 3 (Big-3) and a domain distantly related to plant pathogen response domain 1 (PR-1-like) was identified by screening of an L. rhamnosus HN001 library using HN001 cells as bait. The SpcA-"docking" protein, SpcB, was in turn detected by another phage display library screening, using purified SpcA as bait. SpcB is a 3275-residue cell-surface protein that contains general features of large glycosylated Serine-rich adhesins/fibrils from gram-positive bacteria, including the hallmark signal sequence motif KxYKxGKxW. Both proteins are encoded by genes within a L. rhamnosus-unique gene cluster that distinguishes this species from other lactobacilli. To our knowledge, this is the first example of a secreted-docking protein pair identified in lactobacilli.

Genetic and antigenic characterization of invasive endemic serogroup B Neisseria meningitidis from Ontario, Canada, in 2001-2010

This study examined the antigenic and genetic diversity of serogroup B Neisseria meningitidis (MenB) recovered from invasive meningococcal disease (IMD) cases in Ontario, Canada, over the period 2001-2010 during which no MenB outbreaks had occurred. MenB was found to be responsible for 39 % of all IMD cases, with the remaining cases caused mainly by serogroups Y (28 %), C (23.5 %) and W135 (8 %). One hundred and ninety-three individual MenB case isolates were collected and characterized. Of the 88 sequence types (STs) identified, 75 were grouped into 14 known clonal complexes (CCs), whilst 13 STs were not assigned to any known CC. Fifty-seven different PorA genotypes and 88 STs defined the diversity of invasive MenB in Ontario, which supported the endemic nature of MenB disease in Ontario. Despite the presence of the hypervirulent ST-41/44 and ST-32 CCs, no single ST was predominant and responsible for a large number of IMD cases. Although the Québec outbreak clone of ST-269 was also found in Ontario, the 20 case isolates were genetically diverse: they grouped into seven STs and did not have a predominant PorA genotype. eburst analysis identified a new CC responsible for 14.5 % of the MenB case isolates. The six most common PorA variable region 2 (VR2) genotypes (VR2-9, -4, -14, -16, -13-1 and -16-3) were found in 67 % of invasive MenB isolates.

Phylogenetic incongruence in E. coli O104: understanding the evolutionary relationships of emerging pathogens in the face of homologous recombination

Escherichia coli O104:H4 was identified as an emerging pathogen during the spring and summer of 2011 and was responsible for a widespread outbreak that resulted in the deaths of 50 people and sickened over 4075. Traditional phenotypic and genotypic assays, such as serotyping, pulsed field gel electrophoresis (PFGE), and multilocus sequence typing (MLST), permit identification and classification of bacterial pathogens, but cannot accurately resolve relationships among genotypically similar but pathotypically different isolates. To understand the evolutionary origins of E. coli O104:H4, we sequenced two strains isolated in Ontario, Canada. One was epidemiologically linked to the 2011 outbreak, and the second, unrelated isolate, was obtained in 2010. MLST analysis indicated that both isolates are of the same sequence type (ST678), but whole-genome sequencing revealed differences in chromosomal and plasmid content. Through comprehensive phylogenetic analysis of five O104:H4 ST678 genomes, we identified 167 genes in three gene clusters that have undergone homologous recombination with distantly related E. coli strains. These recombination events have resulted in unexpectedly high sequence diversity within the same sequence type. Failure to recognize or adjust for homologous recombination can result in phylogenetic incongruence. Understanding the extent of homologous recombination among different strains of the same sequence type may explain the pathotypic differences between the ON2010 and ON2011 strains and help shed new light on the emergence of this new pathogen.