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

Characteristics of Meningococcal Invasive Disease in Neonates and Virulence of the Corresponding Isolates

BACKGROUND

The highest incidence of invasive meningococcal disease (IMD) is observed in infants. However, its prevalence in neonates (≤28 days of age) and the characteristics of the corresponding isolates are less described. This report aimed to analyze meningococcal isolates from neonates.

METHODS

We first screened the database of the national reference center for meningococci in France for confirmed neonatal IMD cases between 1999 and 2019. We then performed whole-genome sequencing on all cultured isolates, and we evaluated their virulence in a mouse model.

RESULTS

Fifty-three neonatal cases of IMD (mainly bacteremia) were identified (50 culture-confirmed cases and 3 PCR-confirmed cases) of a total of 10,149 cases (0.5%) but represented 11% of cases among infants of under 1 year of age. Nine cases (17%) occurred among neonates of 3 days of age and younger (early onset). The neonate isolates were often of serogroup B (73.6%) and belonged to the clonal complex CC41/44 (29.4%) with at least 68.5% of coverage by vaccines against serogroup B isolates. The neonatal isolates were able to infect mice although to variable levels.

CONCLUSION

IMD in neonates is not rare and can be of early or late onsets suggesting that anti-meningococcal vaccination can target women planning to have a baby.

Unencapsulated penicillin-resistant Neisseria meningitidis ST-53 clonal complex strain causing symptomatic proctitis in an HIV-negative man who has sex with men: A case report

We report a case of Neisseria meningitidis symptomatic proctitis in an HIV-negative man who has sex with men attending the genitourinary clinic in Malta. The proctitis was caused by a rare clinical unencapsulated penicillin-resistant N. meningitidis strain of the ST-53 clonal complex.

Impact of mandatory vaccination against serogroup C meningococci in targeted and non-targeted populations in France

Since January 2018, mandatory vaccination against meningococci serogroup C has been implemented in France for children <2 years with a recommended catch-up vaccination until the age of 24 years. We aimed to analyse the impact of mandatory vaccination on populations not targeted by it (2–24 years old). We used the database of the national reference centre for meningococci to collect the number of invasive meningococcal disease (IMD) cases before (2016–2017) and after (2018–2019) the mandatory vaccination. The cultured isolates were sequenced and submitted for genomic comparison. The total number of cases was 1706, including 376 cases of IMD serogroup C. Mandatory vaccination correlated with a significant decrease among the <2 years old and a decreasing trend among the 2–14 years old group but not among 15–25 years of age. This observation may be explained by the vaccine coverage that is still low among adolescents and young adults. Moreover, the genomic analysis revealed the co-circulation of two major genotypes belonging to the clonal complex ST-11 whose distribution differed across the age groups in accord with cyclic variations of genotypes. It is important to increase specific knowledge on meningococcal epidemiology and vaccination to involve them in establishing the vaccination strategy.

Evolution of strain coverage by the multicomponent meningococcal serogroup B vaccine (4CMenB) in France

The 4CMenB, a protein-based vaccine, was licensed in Europe in 2013 against invasive meningococcal disease caused by serogroup B and is currently implemented in several countries although according to different national strategies. Isolate coverage estimation is required as vaccine-targeted antigens may vary among isolates over time. Several phenotypic and genotypic methods have been developed to predict strain coverage by scoring the expression and cross-reactivity of vaccine antigens using the Meningococcal Antigen Typing system (MATS), by the genetic correlation of alleles encoding these antigens and MATS expression data (gMATS) and by the Meningococcal Deduced Vaccine Antigen Reactivity (MenDeVAR). We applied these approaches on meningococcal B isolates in France and compared two epidemiological years, 2013-2014 and 2018-2019. A strong correlation was observed between MATS data that were generated for the year 2013-2014 and the gMATS data extracted from whole genome sequencing. gMATS and MenDeVAR were next used to compare the two years. Using gMATS, the overall coverage was 77.2% (lower limit (LL)-upper limit (UL) 66.7-87.7) and 70.7% (LL-UL 61.5-80.0) for the two years, respectively. The reduction in coverage between the two years is mainly driven by the reduction of alleles exactly matching the vaccine antigens. A high number of unpredictable isolates was observed using the MenDeVAR and was due to lack of MATS information for new or rare alleles in particular for the year 2018-2019. Our data underline the need of continuous surveillance of strain coverage and the importance of generating phenotypic MATS data to update the genetic approaches of prediction.

Bioinformatics tools used for whole-genome sequencing analysis of Neisseria gonorrhoeae: a literature review

Whole-genome sequencing (WGS) data are well established for the investigation of gonococcal transmission, antimicrobial resistance prediction, population structure determination and population dynamics. A variety of bioinformatics tools, repositories, services and platforms have been applied to manage and analyze Neisseria gonorrhoeae WGS datasets. This review provides an overview of the various bioinformatics approaches and resources used in 105 published studies (as of 30 April 2021). The challenges in the analysis of N. gonorrhoeae WGS datasets, as well as future bioinformatics requirements, are also discussed.

Predicted coverage by 4CMenB vaccine against invasive meningococcal disease cases in the Netherlands

Neisseria meningitidis serogroup B is a major cause of invasive meningococcal disease in Europe. In the absence of a conjugate serogroup B vaccine, a subcapsular 4CMenB vaccine was developed. Data on 4CMenB vaccine efficacy is still limited. Recently, genomic MATS (Meningococcal Antigen Typing System) was developed as a tool to predict strain coverage, using vaccine antigens sequence data. We characterized all invasive meningococcal isolates received by the Netherlands Reference Laboratory for Bacterial Meningitis (NRLBM) in two epidemiological years 2017-2019 using whole-genome sequencing and determined serogroup, clonal complex (cc) and estimated 4CMenB vaccine coverage by gMATS. Of 396 cases of invasive meningococcal disease, corresponding to an incidence of 1.22 cases/10⁵ inhabitants, 180 (45%) were serogroup W, 155 (39%) serogroup B, 46 (12%) serogroup Y, 10 (3%) serogroup C, 2 non-groupable (0.5%) and 3 (0.7%) unknown. The incidence was the highest among 0-4 years olds (4 cases/10⁵ inhabitants), and 57/72 (79%) of these cases were serogroup B. Serogroup W predominated among persons 45 years of age or older with 110/187 (59%) cases. Serogroup B isolates comprised 11 different clonal complexes, with 103/122 (84%) isolates belonging to 4 clonal complexes: cc32, cc41/44, cc269 and cc213. In contrast, serogroup W isolates were genetically similar with 95% belonging to cc11. Of 122 serogroup B isolates, 89 (73%; 95% CI: 64-80%) were estimated to be covered by 4CMenB and the degree of coverage varied largely by clonal complex and age. Among the 0-4 year olds, 25 of 43 (58%; 95% CI: 43-72%) MenB isolates were estimated to be covered. Since the coverage of the 4CMenB vaccine is dependent on circulating clonal complexes, our findings emphasize the need for surveillance of circulating meningococcal strains. In addition, estimation of age specific coverage is relevant to determine the right target age group for vaccination.

An Improved Approach to Identify Bacterial Pathogens to Human in Environmental Metagenome

The identification of bacterial pathogens to humans is critical for environmental microbial risk assessment. However, current methods for identifying pathogens in environmental samples are limited in their ability to detect highly diverse bacterial communities and accurately differentiate pathogens from commensal bacteria. In the present study, we suggest an improved approach using a combination of identification results obtained from multiple databases, including the multilocus sequence typing (MLST) database, virulence factor database (VFDB), and pathosystems resource integration center (PATRIC) databases to resolve current challenges. By integrating the identification results from multiple databases, potential bacterial pathogens in metagenomes were identified and classified into eight different groups. Based on the distribution of genes in each group, we proposed an equation to calculate the metagenomic pathogen identification index (MPII) of each metagenome based on the weighted abundance of identified sequences in each database. We found that the accuracy of pathogen identification was improved by using combinations of multiple databases compared to that of individual databases. When the approach was applied to environmental metagenomes, metagenomes associated with activated sludge were estimated with higher MPII than other environments (i.e., drinking water, ocean water, ocean sediment, and freshwater sediment). The calculated MPII values were statistically distinguishable among different environments (p<0.05). These results demonstrate that the suggested approach allows more for more accurate identification of the pathogens associated with metagenomes.

Toward a Global Genomic Epidemiology of Meningococcal Disease

Whole-genome sequencing (WGS) is invaluable for studying the epidemiology of meningococcal disease. Here we provide a perspective on the use of WGS for meningococcal molecular surveillance and outbreak investigation, where it helps to characterize pathogens, predict pathogen traits, identify emerging pathogens, and investigate pathogen transmission during outbreaks. Standardization of WGS workflows has facilitated their implementation by clinical and public health laboratories (PHLs), but further development is required for metagenomic shotgun sequencing and targeted sequencing to be widely available for culture-free characterization of bacterial meningitis pathogens. Internet-accessible servers are being established to support bioinformatics analysis, data management, and data sharing among PHLs. However, establishing WGS capacity requires investments in laboratory infrastructure and technical knowledge, which is particularly challenging in resource-limited regions, including the African meningitis belt. Strategic WGS implementation is necessary to monitor the molecular epidemiology of meningococcal disease in these regions and construct a global view of meningococcal disease epidemiology.

Difference in virulence between Neisseria meningitidis serogroups W and Y in transgenic mice

BACKGROUND

Neisseria meningitidis serogroups W and Y are the most common serogroups causing invasive meningococcal disease in Sweden. The majority of cases are caused by the serogroup W UK 2013 strain of clonal complex (cc) 11, and subtype 1 of the serogroup Y, YI strain of cc23. In this study, virulence factors of several lineages within cc11 and cc23 were investigated in transgenic BALB/c mice expressing human transferrin. Transgenic mice were infected intraperitoneally with serogroup W and Y isolates. Levels of bacteria and the proinflammatory cytokine CXCL1 were determined in blood collected 3 h and 24 h post-infection. Apoptosis was investigated in immune cells from peritoneal washes of infected mice. Adhesion and induction of apoptosis in human epithelial cells were also scored.

RESULTS

The levels of bacteraemia, CXCL1, and apoptosis were higher in serogroup W infected mice than in serogroup Y infected mice. Serogroup W isolates also induced higher levels of apoptosis and adhesion in human epithelial cells. No significant differences were observed between different lineages within cc11 and cc23.

CONCLUSIONS

N. meningitidis Serogroup W displayed a higher virulence in vivo in transgenic mice, compared to serogroup Y. This was reflected by higher bacteremia, proinflammatory activity, and ability to induce apoptosis in mouse immune cells and human epithelial cells.

Establishment of the European meningococcal strain collection genome library (EMSC-GL) for the 2011 to 2012 epidemiological year

Invasive meningococcal disease surveillance in Europe combines isolate characterisation and epidemiological data to support public health intervention. A representative European Meningococcal Strain Collection (EMSC) of IMD isolates was obtained, and whole genome sequenced to characterise 799 EMSC isolates from the epidemiological year July 2011-June 2012. To establish a genome library (GL), the isolate information was deposited in the pubMLST.org/neisseria database. Genomes were curated and annotated at 2,429 meningococcal loci, including those defining clonal complex, capsule, antigens, and antimicrobial resistance. Most genomes contained genes encoding B (n = 525; 65.7%) or C (n = 163; 20.4%) capsules; isolates were genetically highly diverse, with >20 genomic lineages, five of which comprising 60.7% (n = 485) of isolates. There were >350 antigenic fine-types: 307 were present once, the most frequent (P1.7-2,4:F5-1) comprised 8% (n = 64) of isolates. Each genome was characterised for Bexsero Antigen Sequence Typing (BAST): 25.5% (n = 204) of isolates contained alleles encoding the fHbp and/or the PorA VR1 vaccine component, but most genomes (n = 513; 64.2%) did not contain the NadA component. EMSC-GL will support an integrated surveillance of disease-associated genotypes in Europe, enabling the monitoring of hyperinvasive lineages, outbreak identification, and supporting vaccine programme implementation.

Typing complex meningococcal vaccines to understand diversity and population structure of key vaccine antigens

Background: Protein-conjugate capsular polysaccharide vaccines can potentially control invasive meningococcal disease (IMD) caused by five (A, C, W, X, Y) of the six IMD-associated serogroups.  Concerns raised by immunological similarity of the serogroup B capsule to human neural cell carbohydrates, meant that 'serogroup B substitute' vaccines target more variable subcapsular protein antigens.  A successful approach using outer membrane vesicles (OMVs) as major vaccine components had limited strain coverage. In 4CMenB (Bexsero ^®), recombinant proteins have been added to ameliorate this problem.  Methods: Scalable, portable, genomic techniques were used to investigate the Bexsero ^® OMV protein diversity in meningococcal populations. Shotgun proteomics identified 461 proteins in the OMV, defining a complex proteome. Amino acid sequences for the 24 proteins most likely to be involved in cross-protective immune responses were catalogued within the PubMLST.org/neisseria database using a novel OMV peptide Typing (OMVT) scheme. Results: Among these proteins there was variation in the extent of diversity and association with meningococcal lineages, identified as clonal complexes (ccs), ranging from the most conserved peptides (FbpA, NEISp0578, and putative periplasmic protein, NEISp1063) to the most diverse (TbpA, NEISp1690).  There were 1752 unique OMVTs identified amongst 2492/3506 isolates examined by whole-genome sequencing (WGS). These OMVTs were grouped into clusters (sharing ≥18 identical OMVT peptides), with 45.3% of isolates assigned to one of 27 OMVT clusters. OMVTs and OMVT clusters were strongly associated with cc, genogroup, and Bexsero ^® antigen variants, demonstrating that combinations of OMV proteins exist in discrete, non-overlapping combinations associated with genogroup and Bexsero ^® Antigen Sequence Type. This highly structured population of IMD-associated meningococci is consistent with strain structure models invoking host immune and/or metabolic selection. Conclusions: The OMVT scheme facilitates region-specific WGS investigation of meningococcal diversity and is an open-access, portable tool with applications for vaccine development, especially in the choice of antigen combinations, assessment and implementation.

Typing complex meningococcal vaccines to understand diversity and population structure of key vaccine antigens

Background: Protein-conjugate capsular polysaccharide vaccines can potentially control invasive meningococcal disease (IMD) caused by five (A, C, W, X, Y) of the six IMD-associated serogroups.  Concerns raised by immunological similarity of the serogroup B capsule to human neural cell carbohydrates, meant that ‘serogroup B substitute’ vaccines target more variable subcapsular protein antigens.  A successful approach using outer membrane vesicles (OMVs) as major vaccine components had limited strain coverage. In 4CMenB (Bexsero ^®), recombinant proteins have been added to ameliorate this problem. 

Methods: Scalable, portable, genomic techniques were used to investigate the Bexsero ^® OMV protein diversity in meningococcal populations. Shotgun proteomics identified 461 proteins in the OMV, defining a complex proteome. Amino acid sequences for the 24 proteins most likely to be involved in cross-protective immune responses were catalogued within the PubMLST.org/neisseria database using a novel OMV peptide Typing (OMVT) scheme.

Results: Among these proteins there was variation in the extent of diversity and association with meningococcal lineages, identified as clonal complexes (ccs), ranging from the most conserved peptides (FbpA, NEISp0578, and putative periplasmic protein, NEISp1063) to the most diverse (TbpA, NEISp1690).  There were 1752 unique OMVTs identified amongst 2492/3506 isolates examined by whole-genome sequencing (WGS). These OMVTs were grouped into clusters (sharing ≥18 identical OMVT peptides), with 45.3% of isolates assigned to one of 27 OMVT clusters. OMVTs and OMVT clusters were strongly associated with cc, genogroup, and Bexsero ^® antigen variants, demonstrating that combinations of OMV proteins exist in discrete, non-overlapping combinations associated with genogroup and Bexsero ^® Antigen Sequence Type. This highly structured population of IMD-associated meningococci is consistent with strain structure models invoking host immune and/or metabolic selection.

Conclusions:The OMVT scheme facilitates region-specific WGS investigation of meningococcal diversity and is an open-access, portable tool with applications for vaccine development, especially in the choice of antigen combinations, assessment and implementation.

Gen2Epi: an automated whole-genome sequencing pipeline for linking full genomes to antimicrobial susceptibility and molecular epidemiological data in Neisseria gonorrhoeae

BACKGROUND

Recent adva1nces in whole genome sequencing (WGS) based technologies have facilitated multi-step applications for predicting antimicrobial resistance (AMR) and investigating the molecular epidemiology of Neisseria gonorrhoeae. However, generating full scaffolds of N. gonorrhoeae genomes from short reads, and the assignment of molecular epidemiological information (NG-MLST, NG-MAST, and NG-STAR) to multiple assembled samples, is challenging due to required manual tasks such as annotating antimicrobial resistance determinants with standard nomenclature for a large number of genomes.

RESULTS

We present Gen2Epi, a pipeline that assembles short reads into full scaffolds and automatically assigns molecular epidemiological and AMR information to the assembled genomes. Gen2Epi is a command-line tool integrating third-party software and tailored specifically for N. gonorrhoeae. For its evaluation, the Gen2Epi pipeline successfully assembled the WGS short reads from 1484 N. gonorrhoeae samples into full-length genomes for both chromosomes and plasmids and was able to assign in silico molecular determinant information to each dataset automatically. The assemblies were generated using raw as well as trimmed short reads. The median genome coverage of full-length scaffolds and "N" statistics (N50, NG50, and NGA50) were higher than, or comparable to, previously published results and the scaffolding process improved the quality of the draft genome assemblies. Molecular antimicrobial resistant (AMR) determinants identified by Gen2Epi reproduced information for the 1484 samples as previously reported, including NG-MLST, NG-MAST, and NG-STAR molecular sequence types.

CONCLUSIONS

Gen2Epi can be used to assemble short reads into full-length genomes and assign accurate molecular marker and AMR information automatically from NG-STAR, NG-MAST, and NG-MLST. Gen2Epi is publicly available under "CC BY-NC 2.0 CA" Creative Commons licensing as a VirtualBox image containing the constituent software components running on the LINUX operating system (CentOS 7). The image and associated documentation are available via anonymous FTP at ftp://www.cs.usask.ca/pub/combi or ftp://ftp.cs.usask.ca/pub/combi.

The Impact of Nucleotide Sequence Analysis on Meningococcal Vaccine Development and Assessment

Since it became available as a routine tool in biology, the determination and analysis of nucleotide sequences has been applied to the design of vaccines and the investigation of their effectiveness. As vaccination is primarily concerned with the interaction of biological molecules with the immune system, the utility of sequence data is not immediately obvious and, indeed, nucleotide sequence data are most effective when used to complement more conventional immunological approaches. Here, the impact of sequencing on the field of vaccinology will be illustrated with reference to the development and implementation of vaccines against Neisseria meningitidis (the meningococcus) over the 30-year period from the late-1980s to the late-2010s. Nucleotide sequence-based studies have been important in the fight against this aggressive pathogen largely because of its high genetic and antigenic diversity, properties that were only fully appreciated because of sequence-based studies. Five aspects will be considered, the use of sequence data to: (i) discover vaccine antigens; (ii) assess the diversity and distribution of vaccine antigens; (iii) determine the evolutionary and population biology of the organism and their implications for immunization; and (iv) develop molecular approaches to investigate pre- and post-vaccine pathogen populations to assess vaccine impact. One of the great advantages of nucleotide sequence data has been its scalability, which has meant that increasingly large data sets have been available, which has proved invaluable in the investigation of an organism as diverse and enigmatic as the meningococcus.

Multilocus Sequence Typing (MLST) of Chlamydiales

Developed two decades ago as a molecular method to provide definite characterization of a bacterial isolate, Multilocus Sequence Typing (MLST) is today globally adopted as a universal fine-detailed molecular typing tool and has been applied to numerous pathogenic and nonpathogenic bacterial as well eukaryotic organisms. MLST utilizes DNA sequence of several conserved housekeeping (HK) genes which are assigned an allelic number, which then collectively constitute an allelic profile or sequence type (ST), a "molecular barcode" of the interrogated bacterial strain or a eukaryotic organism. Here, we describe the principles and molecular approaches for generating MLST data for an analysis of a bacteria in the order Chlamydiales, using a Chlamydia pecorum-specific MLST scheme as an example.

Genomic epidemiology of multidrug-resistant Gram-negative organisms

The emergence and spread of antibiotic-resistant Gram-negative bacteria (rGNB) across global healthcare networks presents a significant threat to public health. As the number of effective antibiotics available to treat these resistant organisms dwindles, it is essential that we devise more effective strategies for controlling their proliferation. Recently, whole-genome sequencing has emerged as a disruptive technology that has transformed our understanding of the evolution and epidemiology of diverse rGNB species, and it has the potential to guide strategies for controlling the evolution and spread of resistance. Here, we review specific areas in which genomics has already made a significant impact, including outbreak investigations, regional epidemiology, clinical diagnostics, resistance evolution, and the study of epidemic lineages. While highlighting early successes, we also point to the next steps needed to translate this technology into strategies to improve public health and clinical medicine.

Open-access bacterial population genomics: BIGSdb software, the PubMLST.org website and their applications

The PubMLST.org website hosts a collection of open-access, curated databases that integrate population sequence data with provenance and phenotype information for over 100 different microbial species and genera.  Although the PubMLST website was conceived as part of the development of the first multi-locus sequence typing (MLST) scheme in 1998 the software it uses, the Bacterial Isolate Genome Sequence database (BIGSdb, published in 2010), enables PubMLST to include all levels of sequence data, from single gene sequences up to and including complete, finished genomes.  Here we describe developments in the BIGSdb software made from publication to June 2018 and show how the platform realises microbial population genomics for a wide range of applications.  The system is based on the gene-by-gene analysis of microbial genomes, with each deposited sequence annotated and curated to identify the genes present and systematically catalogue their variation.  Originally intended as a means of characterising isolates with typing schemes, the synthesis of sequences and records of genetic variation with provenance and phenotype data permits highly scalable (whole genome sequence data for tens of thousands of isolates) means of addressing a wide range of functional questions, including: the prediction of antimicrobial resistance; likely cross-reactivity with vaccine antigens; and the functional activities of different variants that lead to key phenotypes.  There are no limitations to the number of sequences, genetic loci, allelic variants or schemes (combinations of loci) that can be included, enabling each database to represent an expanding catalogue of the genetic variation of the population in question.  In addition to providing web-accessible analyses and links to third-party analysis and visualisation tools, the BIGSdb software includes a RESTful application programming interface (API) that enables access to all the underlying data for third-party applications and data analysis pipelines.

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.

Distinct evolutionary patterns of Neisseria meningitidis serogroup B disease outbreaks at two universities in the USA

Neisseria meningitidis serogroup B (MnB) was responsible for two independent meningococcal disease outbreaks at universities in the USA during 2013. The first at University A in New Jersey included nine confirmed cases reported between March 2013 and March 2014. The second outbreak occurred at University B in California, with four confirmed cases during November 2013. The public health response to these outbreaks included the approval and deployment of a serogroup B meningococcal vaccine that was not yet licensed in the USA. This study investigated the use of whole-genome sequencing(WGS) to examine the genetic profile of the disease-causing outbreak isolates at each university. Comparative WGS revealed differences in evolutionary patterns between the two disease outbreaks. The University A outbreak isolates were very closely related, with differences primarily attributed to single nucleotide polymorphisms/insertion-deletion (SNP/indel) events. In contrast, the University B outbreak isolates segregated into two phylogenetic clades, differing in large part due to recombination events covering extensive regions (>30 kb) of the genome including virulence factors. This high-resolution comparison of two meningococcal disease outbreaks further demonstrates the genetic complexity of meningococcal bacteria as related to evolution and disease virulence.

Persistent Pandemic Lineages of Uropathogenic Escherichia coli in a College Community from 1999 to 2017

The incidence of drug-resistant community-acquired urinary tract infections (CA-UTI) continues to increase worldwide. In 1999 to 2000, a single lineage of uropathogenic Escherichia coli (UPEC) sequence type 69 (ST69) caused 51% of trimethoprim-sulfamethoxazole-resistant UTI in a Northern California university community. We compared the clonal distributions of UPEC and its impact on antimicrobial resistance prevalence in the same community during two periods separated by 17 years. We analyzed E. coli isolates from urine samples from patients with symptoms of UTI who visited a health service between September 2016 and May 2017 and compared them to UPEC isolates collected similarly between October 1999 and March 2000. Isolates were tested for antimicrobial drug susceptibility and genotyped by multilocus sequence typing. In 1999 to 2000, strains belonging to ST95, ST127, ST73, ST69, ST131, and ST10 caused 125 (56%) of 225 UTI cases, while the same STs caused 148 (64%) of 233 UTI cases in 2016 to 2017. The frequencies of ampicillin resistance and ciprofloxacin resistance rose from 24.4% to 41.6% (P < 0.001) and from 0.9% to 5.1% (P < 0.003), respectively. The six STs accounted for 78.6% and 72.7% of these increases, respectively. Prevalence of drug-resistant UTI in this community appears to be largely influenced by a small set of dominant UPEC STs circulating in the same community 17 years apart. Further research to determine the origin and reasons for persistence of these dominant genotypes is necessary to combat antimicrobial-resistant CA-UTI.

Survey on the Use of Whole-Genome Sequencing for Infectious Diseases Surveillance: Rapid Expansion of European National Capacities, 2015-2016

Whole-genome sequencing (WGS) has become an essential tool for public health surveillance and molecular epidemiology of infectious diseases and antimicrobial drug resistance. It provides precise geographical delineation of spread and enables incidence monitoring of pathogens at genotype level. Coupled with epidemiological and environmental investigations, it delivers ultimate resolution for tracing sources of epidemic infections. To ascertain the level of implementation of WGS-based typing for national public health surveillance and investigation of prioritized diseases in the European Union (EU)/European Economic Area (EEA), two surveys were conducted in 2015 and 2016. The surveys were designed to determine the national public health reference laboratories’ access to WGS and operational WGS-based typing capacity for national surveillance of selected foodborne pathogens, antimicrobial-resistant pathogens, and vaccine-preventable diseases identified as priorities for European genomic surveillance. Twenty-eight and twenty-nine out of the 30 EU/EEA countries participated in the survey in 2015 and 2016, respectively. National public health reference laboratories in 22 and 25 countries had access to WGS-based typing for public health applications in 2015 and 2016, respectively. Reported reasons for limited or no access were lack of funding, staff, and expertise. Illumina technology was the most frequently used followed by Ion Torrent technology. The access to bioinformatics expertise and competence for routine WGS data analysis was limited. By mid-2016, half of the EU/EEA countries were using WGS analysis either as first- or second-line typing method for surveillance of the pathogens and antibiotic resistance issues identified as EU priorities. The sampling frame as well as bioinformatics analysis varied by pathogen/resistance issue and country. Core genome multilocus allelic profiling, also called cgMLST, was the most frequently used annotation approach for typing bacterial genomes suggesting potential bioinformatics pipeline compatibility. Further capacity development for WGS-based typing is ongoing in many countries and upon consolidation and harmonization of methods should enable pan-EU data exchange for genomic surveillance in the medium-term subject to the development of suitable data management systems and appropriate agreements for data sharing.

Microbiome at the Frontier of Personalized Medicine

The genomic revolution promises to transform our approach to treat patients by individualizing treatments, reducing adverse events, and decreasing health care costs. The early advances using this have been realized primarily by optimizing preventive and therapeutic approaches in cancer using human genome sequencing. The ability to characterize the microbiome, which includes all the microbes that reside within and upon us and all their genetic elements, using next-generation sequencing allows us to now incorporate this important contributor to human disease into developing new preventive and therapeutic strategies. In this review we highlight the importance of the microbiome in all aspects of human disease, including pathogenesis, phenotype, prognosis, and response to treatment, as well as their role as diagnostic and therapeutic biomarkers. We provide a role for next-generation sequencing in both precise microbial identification of infectious diseases and characterization of microbial communities and their function. Taken together, the microbiome is emerging as an integral part of precision medicine approach as it not only contributes to interindividual variability in all aspects of a disease but also represents a potentially modifiable factor that is amenable to targeting by therapeutics.

Meningococcal serogroup B strain coverage of the multicomponent 4CMenB vaccine with corresponding regional distribution and clinical characteristics in England, Wales, and Northern Ireland, 2007-08 and 2014-15: a qualitative and quantitative assessment

BACKGROUND

The UK introduced 4CMenB-a multicomponent vaccine against serogroup B meningococcal disease-into the national infant immunisation programme in September, 2015. The Meningococcal Antigen Typing System (MATS) was used to estimate coverage by 4CMenB of invasive meningococcal group B isolates obtained during 2007-08 in England and Wales (MATS coverage). We aimed to repeat the MATS survey for invasive meningococcal group B isolates obtained during 2014-15, before 4CMenB introduction; compare strain coverage between 2007-08 and 2014-15; and investigate associations between MATS coverage, age, region, and disease outcomes.

METHODS

Invasive serogroup B meningococcal isolates from cases in England, Wales, and Northern Ireland during 2014-15 were assayed using MATS and compared with 2007-08 data. MATS coverage was assessed by geographical region and age group. Clinical characteristics, risk factors, and outcomes were assessed according to MATS coverage for 2014-15 English cases.

FINDINGS

In 2014-15, 165 of 251 (66%; 95% CI 52-80) meningococcal group B isolates were estimated by MATS to be covered by 4CMenB, compared with 391 of 535 (73%; 95% CI 57-87) in 2007-08. The proportion of MATS-positive isolates with one vaccine antigen increased from 23% (122 of 535) in 2007-08 to 31% (78 of 251) in 2014-15, whereas the proportion with more than one antigen fell from 50% (269 of 535) to 35% (87 of 251). This effect reflected changes in circulating strains, particularly ST-269 clonal complex strains. MATS coverage increased with age, varied by geographical region, and was associated with more severe disease.

INTERPRETATION

In 2014-15, two-thirds of meningococcal group B isolates were predicted to be covered by 4CMenB. Temporal changes in MATS coverage underscore the need for continued monitoring of antigen expression and diversity, particularly in countries with 4CMenB programmes.

FUNDING

Public Health England, GlaxoSmithKline.

MOST: a modified MLST typing tool based on short read sequencing

Multilocus sequence typing (MLST) is an effective method to describe bacterial populations. Conventionally, MLST involves Polymerase Chain Reaction (PCR) amplification of housekeeping genes followed by Sanger DNA sequencing. Public Health England (PHE) is in the process of replacing the conventional MLST methodology with a method based on short read sequence data derived from Whole Genome Sequencing (WGS). This paper reports the comparison of the reliability of MLST results derived from WGS data, comparing mapping and assembly-based approaches to conventional methods using 323 bacterial genomes of diverse species. The sensitivity of the two WGS based methods were further investigated with 26 mixed and 29 low coverage genomic data sets from Salmonella enteridis and Streptococcus pneumoniae. Of the 323 samples, 92.9% (n = 300), 97.5% (n = 315) and 99.7% (n = 322) full MLST profiles were derived by the conventional method, assembly- and mapping-based approaches, respectively. The concordance between samples that were typed by conventional (92.9%) and both WGS methods was 100%. From the 55 mixed and low coverage genomes, 89.1% (n = 49) and 67.3% (n = 37) full MLST profiles were derived from the mapping and assembly based approaches, respectively. In conclusion, deriving MLST from WGS data is more sensitive than the conventional method. When comparing WGS based methods, the mapping based approach was the most sensitive. In addition, the mapping based approach described here derives quality metrics, which are difficult to determine quantitatively using conventional and WGS-assembly based approaches.

Global Meningococcal Initiative: guidelines for diagnosis and confirmation of invasive meningococcal disease

The Global Meningococcal Initiative (GMI) is an international group of scientists and clinicians with recognized expertise in meningococcal disease including microbiology, immunology, epidemiology, public health and vaccinology. The GMI was established to promote the global prevention of meningococcal disease through education, research and international cooperation. The GMI held its second summit meeting in 2013 to discuss the different aspects of existing meningococcal immunization programmes and surveillance systems. Laboratory confirmation and characterization were identified as essential for informing evidence-based vaccine implementation decisions. The relative merits of different confirmatory methodologies and their applications in different resource settings were a key component of the discussions. This paper summarizes the salient issues discussed, with special emphasis on the recommendations made and any deficiencies that were identified.

Meningococcal vaccine antigen diversity in global databases

The lack of an anti-capsular vaccine against serogroup B meningococcal disease has necessitated the exploration of alternative vaccine candidates, mostly proteins exhibiting varying degrees of antigenic variation. Analysis of variants of antigen-encoding genes is facilitated by publicly accessible online sequence repositories, such as the Neisseria PubMLST database and the associated Meningitis Research Foundation Meningococcus Genome Library (MRF-MGL). We investigated six proposed meningococcal vaccine formulations by deducing the prevalence of their components in the isolates represented in these repositories. Despite high diversity, a limited number of antigenic variants of each of the vaccine antigens were prevalent, with strong associations of particular variant combinations with given serogroups and genotypes. In the MRF-MGL and globally, the highest levels of identical sequences were observed with multicomponent/multivariant vaccines. Our analyses further demonstrated that certain combinations of antigen variants were prevalent over periods of decades in widely differing locations, indicating that vaccine formulations containing a judicious choice of antigen variants have potential for long-term protection across geographic regions. The data further indicated that formulations with multiple variants would be especially relevant at times of low disease incidence, as relative diversity was higher. Continued surveillance is required to monitor the changing prevalence of these vaccine antigens.

Population and Functional Genomics of Neisseria Revealed with Gene-by-Gene Approaches

Rapid low-cost whole-genome sequencing (WGS) is revolutionizing microbiology; however, complementary advances in accessible, reproducible, and rapid analysis techniques are required to realize the potential of these data. Here, investigations of the genus Neisseria illustrated the gene-by-gene conceptual approach to the organization and analysis of WGS data. Using the gene and its link to phenotype as a starting point, the BIGSdb database, which powers the PubMLST databases, enables the assembly of large open-access collections of annotated genomes that provide insight into the evolution of the Neisseria, the epidemiology of meningococcal and gonococcal disease, and mechanisms of Neisseria pathogenicity.

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.

Genomic Analysis of Serogroup Y Neisseria meningitidis Isolates Reveals Extensive Similarities Between Carriage-Associated and Disease-Associated Organisms

BACKGROUND

Neisseria meningitidis is a frequent colonizer of the human nasopharynx, with asymptomatic carriage providing the reservoir for invasive, disease-causing strains. Serogroup Y (MenY) strains are a major cause of meningococcal disease. High-resolution genetic analyses of carriage and disease isolates can establish epidemiological relationships and identify potential virulence factors.

METHODS

Whole-genome sequence data were obtained for 99 MenY carriage isolates recovered in the United Kingdom during 1997-2010. Sequences were compared to those of 73 MenY invasive isolates recovered during 2010-2011, using a gene-by-gene approach.

RESULTS

Comparisons across 1605 core genes resolved 91% of isolates into one of 8 clusters containing closely related disease and carriage isolates. Six clusters contained carried meningococci isolated during 1997-2001, suggesting temporal stability. One cluster of isolates, predominately sharing the designation Y: P1.5-1,10-1: F4-1: ST-1655 (cc23), was resolved into one subcluster with 86% carriage isolates and a second with 90% invasive isolates. These subclusters were defined by specific allelic differences in 5 core genes encoding glycerate kinase (glxK), valine-pyruvate transaminase (avtA), superoxide dismutase (sodB), and 2 hypothetical proteins.

CONCLUSIONS

High-resolution genetic analyses detected long-term temporal stability and temporally overlapping carriage and disease populations for MenY clones but also evidence of a disease-associated clone.

Whole genome multilocus sequence typing as an epidemiologic tool for Yersinia pestis

Human plague is a severe and often fatal zoonotic disease caused by Yersinia pestis. For public health investigations of human cases, nonintensive whole genome molecular typing tools, capable of defining epidemiologic relationships, are advantageous. Whole genome multilocus sequence typing (wgMLST) is a recently developed methodology that simplifies genomic analyses by transforming millions of base pairs of sequence into character data for each gene. We sequenced 13 US Y. pestis isolates with known epidemiologic relationships. Sequences were assembled de novo, and multilocus sequence typing alleles were assigned by comparison against 3979 open reading frames from the reference strain CO92. Allele-based cluster analysis accurately grouped the 13 isolates, as well as 9 publicly available Y. pestis isolates, by their epidemiologic relationships. Our findings indicate wgMLST is a simplified, sensitive, and scalable tool for epidemiologic analysis of Y. pestis strains.

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.

BacWGSTdb, a database for genotyping and source tracking bacterial pathogens

Whole genome sequencing has become one of the routine methods in molecular epidemiological practice. In this study, we present BacWGSTdb (http://bacdb.org/BacWGSTdb), a bacterial whole genome sequence typing database which is designed for clinicians, clinical microbiologists and hospital epidemiologists. This database borrows the population structure from the current multi-locus sequence typing (MLST) scheme and adopts a hierarchical data structure: species, clonal complex and isolates. When users upload the pre-assembled genome sequences to BacWGSTdb, it offers the functionality of bacterial genotyping at both traditional MLST and whole-genome levels. More importantly, users are told which isolates in the public database are phylogenetically close to the query isolate, along with their clinical information such as host, isolation source, disease, collection time and geographical location. In this way, BacWGSTdb offers a rapid and convenient platform for worldwide users to address a variety of clinical microbiological issues such as source tracking bacterial pathogens.

Draft Genome Sequence of Dichelobacter nodosus ATCC 25549, Strain VPI 2340 [11342], a Bacterium Causing Footrot in Sheep

We report a draft genome sequence for Dichelobacter nodosus ATCC 25549, strain VPI 2340 [11342], a causative agent of ovine footrot. The draft genome shares ~98% gene similarity with the available genome of D. nodosus strain VCS1703A but is differentiated by extensive gene duplication and the absence of 13 particular genes.

Comparison of Phenotypic and Genotypic Approaches to Capsule Typing of Neisseria meningitidis by Use of Invasive and Carriage Isolate Collections

Neisseria meningitidis serogroup B (MnB) is a leading cause of bacterial meningitis; however, MnB is most commonly associated with asymptomatic carriage in the nasopharyngeal cavity, as opposed to the disease state. Two vaccines are now licensed for the prevention of MnB disease; a possible additional benefit of these vaccines could be to protect against disease indirectly by disrupting nasopharyngeal carriage (e.g., herd protection). To investigate this possibility, accurate diagnostic approaches to characterize MnB carriage isolates are required. In contrast to invasive meningococcal disease (IMD) isolates, which can be readily serogrouped, carriage isolates often lack capsule expression, making standard phenotypic assays unsuitable for strain characterization. Several antibody-based methods were evaluated for their abilities to serogroup isolates and were compared with two genotyping methods (real-time PCR [rt-PCR] and whole-genome sequencing [WGS]) to identify which approach would most accurately ascertain the polysaccharide groups associated with carriage isolates. WGS and rt-PCR were in agreement for 99% of IMD isolates, including those with coding sequences for MnB, MnC, MnW, and MnY, and the phenotypic methods correctly identified serogroups for 69 to 98% of IMD isolates. In contrast, only 47% of carriage isolates were groupable by genotypic methods, due to mutations within the capsule operon; of the isolates identified by genotypic methods, ≤43% were serogroupable with any of the phenotypic methods tested. These observations highlight the difficulties in the serogrouping and capsular genogrouping of meningococcal carriage isolates. Based on our findings, WGS is the most suitable approach for the characterization of meningococcal carriage isolates.

Genome-Based Characterization of Emergent Invasive Neisseria meningitidis Serogroup Y Isolates in Sweden from 1995 to 2012

Invasive meningococcal disease (IMD) caused by Neisseria meningitidis serogroup Y has increased in Europe, especially in Scandinavia. In Sweden, serogroup Y is now the dominating serogroup, and in 2012, the serogroup Y disease incidence was 0.46/100,000 population. We previously showed that a strain type belonging to sequence type 23 was responsible for the increased prevalence of this serogroup in Sweden. The objective of this study was to investigate the serogroup Y emergence by whole-genome sequencing and compare the meningococcal population structure of Swedish invasive serogroup Y strains to those of other countries with different IMD incidence. Whole-genome sequencing was performed on invasive serogroup Y isolates from 1995 to 2012 in Sweden (n = 186). These isolates were compared to a collection of serogroup Y isolates from England, Wales, and Northern Ireland from 2010 to 2012 (n = 143), which had relatively low serogroup Y incidence, and two isolates obtained in 1999 in the United States, where serogroup Y remains one of the major causes of IMD. The meningococcal population structures were similar in the investigated regions; however, different strain types were prevalent in each geographic region. A number of genes known or hypothesized to have an impact on meningococcal virulence were shown to be associated with different strain types and subtypes. The reasons for the IMD increase are multifactorial and are influenced by increased virulence, host adaptive immunity, and transmission. Future genome-wide association studies are needed to reveal additional genes associated with serogroup Y meningococcal disease, and this work would benefit from a complete serogroup Y meningococcal reference genome.

Using multilocus sequence typing to study bacterial variation: prospects in the genomic era

Multilocus sequence typing (MLST) indexes the sequence variation present in a small number (usually seven) of housekeeping gene fragments located around the bacterial genome. Unique alleles at these loci are assigned arbitrary integer identifiers, which effectively summarizes the variation present in several thousand base pairs of genome sequence information as a series of numbers. Comparing bacterial isolates using allele-based methods efficiently corrects for the effects of lateral gene transfer present in many bacterial populations and is computationally efficient. This 'gene-by-gene' approach can be applied to larger collections of loci, such as the ribosomal protein genes used in ribosomal MLST (rMLST), up to and including the complete set of coding sequences present in a genome, whole-genome MLST (wgMLST), providing scalable, efficient and readily interpreted genome analysis.

The Diversity of Meningococcal Carriage Across the African Meningitis Belt and the Impact of Vaccination With a Group A Meningococcal Conjugate Vaccine

Background. Study of meningococcal carriage is essential to understanding the epidemiology of Neisseria meningitidis infection.

Methods. Twenty cross-sectional carriage surveys were conducted in 7 countries in the African meningitis belt; 5 surveys were conducted after introduction of a new serogroup A meningococcal conjugate vaccine (MenAfriVac). Pharyngeal swab specimens were collected, and Neisseria species were identified by microbiological and molecular techniques.

Results. A total of 1687 of 48 490 participants (3.4%; 95% confidence interval [CI], 3.2%–3.6%) carried meningococci. Carriage was more frequent in individuals aged 5–14 years, relative to those aged 15–29 years (adjusted odds ratio [OR], 1.41; 95% CI, 1.25–1.60); in males, relative to females (adjusted OR, 1.17; 95% CI, 1.10–1.24); in individuals in rural areas, relative to those in urban areas (adjusted OR, 1.44; 95% CI, 1.28–1.63); and in the dry season, relative to the rainy season (adjusted OR, 1.54; 95% CI, 1.37–1.75). Forty-eight percent of isolates had genes encoding disease-associated polysaccharide capsules; genogroup W predominated, and genogroup A was rare. Strain diversity was lower in countries in the center of the meningitis belt than in Senegal or Ethiopia. The prevalence of genogroup A fell from 0.7% to 0.02% in Chad following mass vaccination with MenAfriVac.

Conclusions. The prevalence of meningococcal carriage in the African meningitis belt is lower than in industrialized countries and is very diverse and dynamic, even in the absence of vaccination.

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.

Transoceanic spreading of pathogenic strains of Vibrio parahaemolyticus with distinctive genetic signatures in the recA gene

Vibrioparahaemolyticus is an important human pathogen whose transmission is associated with the consumption of contaminated seafood. Consistent multilocus sequence typing for V. parahaemolyticus has shown difficulties in the amplification of the recA gene by PCR associated with a lack of amplification or a larger PCR product than expected. In one strain (090–96, Peru, 1996), the produced PCR product was determined to be composed of two recA fragments derived from different Vibrio species. To better understand this phenomenon, we sequenced the whole genome of this strain. The hybrid recA gene was found to be the result of a fragmentation of the original lineage-specific recA gene resulting from a DNA insertion of approximately 30 kb in length. This insert had a G+C content of 38.8%, lower than that of the average G+C content of V. parahaemolyticus (45.2%), and contained 19 ORFs, including a complete recA gene. This new acquired recA gene deviated 24% in sequence from the original recA and was distantly related to recA genes from bacteria of the Vibrionaceae family. The reconstruction of the original recA gene (recA3) identified the precursor as belonging to ST189, a sequence type reported previously only in Asian countries. The identification of this singular genetic feature in strains from Asia reveals new evidence for genetic connectivity between V. parahaemolyticus populations at both sides of the Pacific Ocean that, in addition to the previously described pandemic clone, supports the existence of a recurrent transoceanic spreading of pathogenic V. parahaemolyticus with the corresponding potential risk of pandemic expansion.

Using surface-enhanced Raman spectroscopy and electrochemically driven melting to discriminate Yersinia pestis from Y. pseudotuberculosis based on single nucleotide polymorphisms within unpurified polymerase chain reaction amplicons

The development of sensors for the detection of pathogen-specific DNA, including relevant species/strain level discrimination, is critical in molecular diagnostics with major impacts in areas such as bioterrorism and food safety. Herein, we use electrochemically driven denaturation assays monitored by surface-enhanced Raman spectroscopy (SERS) to target single nucleotide polymorphisms (SNPs) that distinguish DNA amplicons generated from Yersinia pestis, the causative agent of plague, from the closely related species Y. pseudotuberculosis. Two assays targeting SNPs within the groEL and metH genes of these two species have been successfully designed. Polymerase chain reaction (PCR) was used to produce Texas Red labeled single-stranded DNA (ssDNA) amplicons of 262 and 251 bases for the groEL and metH targets, respectively. These amplicons were used in an unpurified form to hybridize to immobilized probes then subjected to electrochemically driven melting. In all cases electrochemically driven melting was able to discriminate between fully homologous DNA and that containing SNPs. The metH assay was particularly challenging due to the presence of only a single base mismatch in the middle of the 251 base long PCR amplicon. However, manipulation of assay conditions (conducting the electrochemical experiments at 10 °C) resulted in greater discrimination between the complementary and mismatched DNA. Replicate data were collected and analyzed for each duplex on different days, using different batches of PCR product and different sphere segment void (SSV) substrates. Despite the variability introduced by these differences, the assays are shown to be reliable and robust providing a new platform for strain discrimination using unpurified PCR samples.

Cronobacter, the emergent bacterial pathogen Enterobacter sakazakii comes of age; MLST and whole genome sequence analysis

BACKGROUND

Following the association of Cronobacter spp. to several publicized fatal outbreaks in neonatal intensive care units of meningitis and necrotising enterocolitis, the World Health Organization (WHO) in 2004 requested the establishment of a molecular typing scheme to enable the international control of the organism. This paper presents the application of Next Generation Sequencing (NGS) to Cronobacter which has led to the establishment of the Cronobacter PubMLST genome and sequence definition database (http://pubmlst.org/cronobacter/) containing over 1000 isolates with metadata along with the recognition of specific clonal lineages linked to neonatal meningitis and adult infections

RESULTS

Whole genome sequencing and multilocus sequence typing (MLST) has supports the formal recognition of the genus Cronobacter composed of seven species to replace the former single species Enterobacter sakazakii. Applying the 7-loci MLST scheme to 1007 strains revealed 298 definable sequence types, yet only C. sakazakii clonal complex 4 (CC4) was principally associated with neonatal meningitis. This clonal lineage has been confirmed using ribosomal-MLST (51-loci) and whole genome-MLST (1865 loci) to analyse 107 whole genomes via the Cronobacter PubMLST database. This database has enabled the retrospective analysis of historic cases and outbreaks following re-identification of those strains.

CONCLUSIONS

The Cronobacter PubMLST database offers a central, open access, reliable sequence-based repository for researchers. It has the capacity to create new analysis schemes 'on the fly', and to integrate metadata (source, geographic distribution, clinical presentation). It is also expandable and adaptable to changes in taxonomy, and able to support the development of reliable detection methods of use to industry and regulatory authorities. Therefore it meets the WHO (2004) request for the establishment of a typing scheme for this emergent bacterial pathogen. Whole genome sequencing has additionally shown a range of potential virulence and environmental fitness traits which may account for the association of C. sakazakii CC4 pathogenicity, and propensity for neonatal CNS.

SRST2: Rapid genomic surveillance for public health and hospital microbiology labs

Rapid molecular typing of bacterial pathogens is critical for public health epidemiology, surveillance and infection control, yet routine use of whole genome sequencing (WGS) for these purposes poses significant challenges. Here we present SRST2, a read mapping-based tool for fast and accurate detection of genes, alleles and multi-locus sequence types (MLST) from WGS data. Using >900 genomes from common pathogens, we show SRST2 is highly accurate and outperforms assembly-based methods in terms of both gene detection and allele assignment. We include validation of SRST2 within a public health laboratory, and demonstrate its use for microbial genome surveillance in the hospital setting. In the face of rising threats of antimicrobial resistance and emerging virulence among bacterial pathogens, SRST2 represents a powerful tool for rapidly extracting clinically useful information from raw WGS data. Source code is available from http://katholt.github.io/srst2/.

SRST2: Rapid genomic surveillance for public health and hospital microbiology labs

Rapid molecular typing of bacterial pathogens is critical for public health epidemiology, surveillance and infection control, yet routine use of whole genome sequencing (WGS) for these purposes poses significant challenges. Here we present SRST2, a read mapping-based tool for fast and accurate detection of genes, alleles and multi-locus sequence types (MLST) from WGS data. Using >900 genomes from common pathogens, we show SRST2 is highly accurate and outperforms assembly-based methods in terms of both gene detection and allele assignment. We include validation of SRST2 within a public health laboratory, and demonstrate its use for microbial genome surveillance in the hospital setting. In the face of rising threats of antimicrobial resistance and emerging virulence among bacterial pathogens, SRST2 represents a powerful tool for rapidly extracting clinically useful information from raw WGS data. Source code is available from http://katholt.github.io/srst2/.

Clostridium difficile surveillance: harnessing new technologies to control transmission

Clostridium difficile surveillance allows outbreaks of cases clustered in time and space to be identified and further transmission prevented. Traditionally, manual detection of groups of cases diagnosed in the same ward or hospital, often followed by retrospective reference laboratory genotyping, has been used to identify outbreaks. However, integrated healthcare databases offer the prospect of automated real-time outbreak detection based on statistically robust methods, and accounting for contacts between cases, including those distant to the ward of diagnosis. Complementary to this, rapid benchtop whole genome sequencing, and other highly discriminatory genotyping, has the potential to distinguish which cases are part of an outbreak with high precision and in clinically relevant timescales. These new technologies are likely to shape future surveillance.

Real-time whole-genome sequencing for routine typing, surveillance, and outbreak detection of verotoxigenic Escherichia coli

Fast and accurate identification and typing of pathogens are essential for effective surveillance and outbreak detection. The current routine procedure is based on a variety of techniques, making the procedure laborious, time-consuming, and expensive. With whole-genome sequencing (WGS) becoming cheaper, it has huge potential in both diagnostics and routine surveillance. The aim of this study was to perform a real-time evaluation of WGS for routine typing and surveillance of verocytotoxin-producing Escherichia coli (VTEC). In Denmark, the Statens Serum Institut (SSI) routinely receives all suspected VTEC isolates. During a 7-week period in the fall of 2012, all incoming isolates were concurrently subjected to WGS using IonTorrent PGM. Real-time bioinformatics analysis was performed using web-tools (www.genomicepidemiology.org) for species determination, multilocus sequence type (MLST) typing, and determination of phylogenetic relationship, and a specific VirulenceFinder for detection of E. coli virulence genes was developed as part of this study. In total, 46 suspected VTEC isolates were characterized in parallel during the study. VirulenceFinder proved successful in detecting virulence genes included in routine typing, explicitly verocytotoxin 1 (vtx1), verocytotoxin 2 (vtx2), and intimin (eae), and also detected additional virulence genes. VirulenceFinder is also a robust method for assigning verocytotoxin (vtx) subtypes. A real-time clustering of isolates in agreement with the epidemiology was established from WGS, enabling discrimination between sporadic and outbreak isolates. Overall, WGS typing produced results faster and at a lower cost than the current routine. Therefore, WGS typing is a superior alternative to conventional typing strategies. This approach may also be applied to typing and surveillance of other pathogens.

Genome sequence analyses show that Neisseria oralis is the same species as 'Neisseria mucosa var. heidelbergensis'

Phylogenies generated from whole genome sequence (WGS) data provide definitive means of bacterial isolate characterization for typing and taxonomy. The species status of strains recently defined with conventional taxonomic approaches as representing Neisseria oralis was examined by the analysis of sequences derived from WGS data, specifically: (i) 53 Neisseria ribosomal protein subunit (rps) genes (ribosomal multi-locus sequence typing, rMLST); and (ii) 246 Neisseria core genes (core genome MLST, cgMLST). These data were compared with phylogenies derived from 16S and 23S rRNA gene sequences, demonstrating that the N. oralis strains were monophyletic with strains described previously as representing 'Neisseria mucosa var. heidelbergensis' and that this group was of equivalent taxonomic status to other well-described species of the genus Neisseria. Phylogenetic analyses also indicated that Neisseria sicca and Neisseria macacae should be considered the same species as Neisseria mucosa and that Neisseria flavescens should be considered the same species as Neisseria subflava. Analyses using rMLST showed that some strains currently defined as belonging to the genus Neisseria were more closely related to species belonging to other genera within the family; however, whole genome analysis of a more comprehensive selection of strains from within the family Neisseriaceae would be necessary to confirm this. We suggest that strains previously identified as representing 'N. mucosa var. heidelbergensis' and deposited in culture collections should be renamed N. oralis. Finally, one of the strains of N. oralis was able to ferment lactose, due to the presence of β-galactosidase and lactose permease genes, a characteristic previously thought to be unique to Neisseria lactamica, which therefore cannot be thought of as diagnostic for this species; however, the rMLST and cgMLST analyses confirm that N. oralis is most closely related to N. mucosa.

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.