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

Holly B. Bratcher, Carina Brehony, Sigrid Heuberger, Despo

Pieridou-Bagatzouni, Pavla Křížová, Steen Hoffmann, Maija Toropainen, Muhamed-Kheir Taha, Heike Claus, Georgina Tzanakaki, et al.

<figure><img src="image_1.png" /></figure>

► To cite this version:

Holly B. Bratcher, Carina Brehony, Sigrid Heuberger, Despo Pieridou-Bagatzouni, Pavla Křížová, et al. Establishment of the European meningococcal strain collection genome library (EMSC-GL) for the 2011 to 2012 epidemiological year. Eurosurveillance, 2018, 23 (20), 10.2807/1560-7917.ES.2018.23.20.17-00474. pasteur-02013211

HAL Id: pasteur-02013211

https://pasteur.hal.science/pasteur-02013211

Submitted on 10 Feb 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Establishment of the European meningococcal strain collection genome library (EMSC-GL) for the 2011 to 2012 epidemiological year

Holly B Bratcher¹, Carina Brehony¹,², Sigrid Heuberger³, Despo Pieridou-Bagatzouni⁴, Pavla Křížová⁵, Steen Hoffmann⁶, Maija Toropainen⁷, Muhamed-Kheir Taha⁸, Heike Claus⁹, Georgina Tzanakaki¹⁰, Tímea Erdősi¹¹, Jelena Galajeva¹², Arie van der Ende¹³, Anna Skoczyńska¹⁴, Marina Pana¹⁵, Alena Vaculíková¹⁶, Metka Paragi¹⁷, Martin CJ Maiden¹, Dominique A Caugant¹⁸

1. University of Oxford, Department of Zoology, Oxford, United Kingdom

2. Clinical Science Institute, National University of Ireland, Galway, Republic of Ireland

3. Austrian Agency for Health and Food Safety, Graz, Austria

4. Nicosia General Hospital, Nicosia, Cyprus

5. National Institute of Public Health, Prague, Czech Republic

6. Statens Serum Institute, Copenhagen, Denmark

7. National Institute for Health and Welfare, Helsinki, Finland

8. Pasteur Institute, Paris, France

9. University of Würzburg, Würzburg, Germany

10. National Meningitis Reference Laboratory, Athens, Greece

11. National Center for Epidemiology, Budapest, Hungary

12. Riga East University Hospital, Riga, Latvia

13. Academic Medical Center, University of Amsterdam, The Netherlands

14. National Medicines Institute, Warsaw, Poland

15. National Institute of Research and Development for Microbiology and Immunology, Bucharest, Romania

16. Public Health Authority of the Slovak Republic, Bratislava, Slovakia

17. National Institute of Public Health, Ljubljana, Slovenia

18. Norwegian Institute of Public Health, Oslo, Norway

Correspondence: Holly B. Bratcher (holly.bratcher@zoo.ox.ac.uk)

**Citation style for this article:**

Bratcher Holly B, Brehony Carina, Heuberger Sigrid, Pieridou-Bagatzouni Despo, Křížová Pavla, Hoffmann Steen, Toropainen Maija, Taha Muhamed-Kheir, Claus Heike, Tzanakaki Georgina, Erdősi Tímea, Galajeva Jelena, van der Ende Arie, Skoczyńska Anna, Pana Marina, Vaculíková Alena, Paragi Metka, Maiden Martin CJ, Caugant Dominique A. Establishment of the European meningococcal strain collection genome library (EMSC-GL) for the 2011 to 2012 epidemiological year. Euro Surveill. 2018;23(20):pii=17-00474. https://doi.org/10.2807/1560-7917.ES.2018.23.20.17-00474

Article submitted on 14 Jul 2017 / accepted on 27 Oct 2017 / published on 17 May 2018

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.

## Introduction

Invasive meningococcal disease (IMD) is a contagious, occasionally epidemic, infectious disease that causes significant morbidity and mortality with varying incidence globally [1]. Only a small proportion of individuals infected with Neisseria meningitidis, the meningococcus, develop severe disease, usually meningitis and/or septicaemia, with survivors often suffering permanent sequela [2]. While vaccination with protein-conjugated polysaccharides of serogroup A, C, W, and Y is an effective prevention strategy for the main disease-causing serogroups, currently there is no comprehensive vaccine against meningococcal disease. The structural similarity of the serogroup B polysaccharide capsule to the human neural cell adhesion molecule renders it poorly immunogenic and carries concerns over possible autoimmune reactions [3,4]. To circumvent this problem, vaccines developed to control serogroup B, the main cause of IMD in Europe, have
FIGURE 1

Breakdown by clonal complex of invasive meningococcal disease isolates submitted by 16 European countries for the epidemiological year July 2011–June 2012 (n = 799 isolates)

<figure><img src="image_3.png" /></figure>

ST: sequence type.

The most commonly observed were clonal complexes for the epidemiological year were ST-41/44, ST-32 and ST-11. The isolates unassigned to a clonal complex included four isolates without an ST and 108 isolates with a ST that was not associated with a clonal complex.

been based on outer-membrane vesicles and more
recently on recombinant proteins [5,6].

Understanding the impact of meningococcal vaccines requires surveillance that includes knowledge of disease-associated meningococci and the antigens and capsule types that they possess. However, comprehensive and timely surveillance of IMD is complicated by the high genetic and antigenic diversity of meningococcal populations [7]. While most meningococcal disease is caused by a minority of meningococcal capsule types, there are many variants of the subcapular antigens that are used in vaccines [8]. However, relatively few meningococcal genotypes, known as hyperinvasive lineages which have been identified by multilocus sequence typing (MLST) as particular clonal complexes (ccs), cause most IMD and these have been shown to be associated with particular variants of antigens

and combination of vaccine antigen variants [7,9]. Conventional isolate typing has required multiple tests involving sequencing of multiple genes [10] but whole genome sequencing (WGS) provides the potential to identify all of the salient features in a single experiment [11].

The establishment of a European meningococcal genome library has potential to enhance IMD surveillance continent-wide, providing comprehensive and comparable long-term typing data. To do this, all national reference laboratories for meningococci in Europe were asked to make available in 2014 isolates from one epidemiological year (2011–2012) to contribute to the European Meningococcal Strain Collection (EMSC) in order to provide a bank of isolate genomes representative of those causing invasive disease in Europe at that time. The EMSC Genome Library
FIGURE 2

Putative vaccine antigen coverage of invasive meningococcal disease isolates, as assessed by a direct sequence match to any one of four vaccine components, Europe, July 2011–June 2012 (n = 349 isolates<sup>a</sup>)

<figure><img src="image_4.png" /></figure>

BAST: Bexsero Antigen Sequence Type.

The vaccine components are NadA, NHBA, PorA VR1 and PorA VR2.

(EMSC-GL) was created to exploit whole genome sequence (WGS) data for surveillance, characterisation, and vaccine intervention studies to analyse longitudinal data using both historical, current and future datasets. The WGS data are hosted on the pubMLST.org/neisseria database, where the analysis tools of the BIGSdb platform can be used [12].

Methods

Specimen collection

The EMSC was established in 2002 at the Norwegian Institute of Public Health, Oslo, Norway with the support of the European Centre for Disease Prevention and Control (ECDC), to assure preservation of meningococcal isolates causing disease in Europe. Additional funding was allocated by the ECDC in 2011 for the implementation of the EMSC-GL. Participation of each country in the creation of the EMSC-GL was voluntary. All 799 invasive meningococcal isolates received from the epidemiological year July 2011–June 2012, from 16 participating national reference laboratories, were included in the EMSC-GL.

Culture and DNA extraction

Isolates were grown on blood agar plates overnight at 37°C in an atmosphere of 5% CO₂ and DNA extracted using an epMotion 5075 pipetting robot (Eppendorf) and NucleoSpin 8 Tissue kit (Macherey-Nagel) according to the manufacturer's instructions. DNA was visualised with the FlashGel DNA system (Lonza), its quantity assessed using Qubit (Invitrogen, ThermoFisher Scientific Inc), and sent frozen in 96-well plates to the Oxford Genomics Centre (OGC), University of Oxford, United Kingdom (UK).

Whole genome sequencing, de novo assembly,
annotation and curation

The DNA samples were normalised and sequenced at the OGC. Quantification and integrity of DNA extracts were assessed using the Invitrogen PicoGreen assay on a 1% E-GelEX. DNA was fragmented using acoustic shearing and the fragment distribution determined using the Agilent Tapestation D1200 system. NEBNext DNA Sample Prep Master (New England BioLabs) was used with minor modifications to construct the libraries. Ligation of the adapters was completed using

<footnote><sup>a</sup> The graph depicts the 44 BASTs found with three or more examples, a total of 332 BAST were found in the dataset. Putative multiple component coverage BAST include: three exact sequence match vaccine components (BAST-223); two exact sequence match vaccine components (BAST-84, 2448, 1071, 1072, and 1077); and one exact sequence match vaccine component (BAST-4, 5, 10, 94, 239, 1014, 1016, 1074, 1271, and 1313).</footnote>
Neighbour-net population model analysis of invasive meningococcal disease isolates, Europe, July 2011–June 2012 (n = 799 isolates)

<figure><img src="image_5.png" /></figure>

cc: clonal complex.

Core genome multilocus sequence typing (cgMLST, vr1) and lineage assignment were performed using the pubMLST Genome Comparator Tool and visualised with SplitsTree4. Correlation of cc and lineage assignment are as follows: lineage 3 (ST-41/44 cc), lineage 8 (ST-8 cc), lineage 11 (ST-11 cc), lineage 22 (ST-22 cc), lineage 23 (ST-23 cc), lineage 2 (ST-269 cc), lineage 35 (ST-35 cc), lineage 5 (ST-32 cc), lineage 25 (ST-162 cc), lineage 13 (ST-213 cc), lineage 18 (ST-18 cc), lineage 6 (ST-60 cc), lineage 20 (ST-103 cc). Unlabelled genomes are currently unassigned to a defined lineage but may be associated with a defined cc. Scale bar represents the number of allele pairwise differences.

Illumina Adapters and PCR enriched. Libraries were pooled and quantified using the MX3005PTM instrument (Agilent) and the Agilent qPCR Library Quantification Kit. Paired-end 100 base reads were sequenced using the Illumina HiSeq 2000 sequencing platform.

Isolate records were created in the pubMLST.org/neisseria database and included provenance and phenotypic typing information [12]. Short-read sequence data were assembled de novo as described previously [13], and, uploaded to the database and then linked to their respective isolate record. Genome assembly statistics

were generated as part of this automated pipeline. The data are accessible on pubMLST.org/neisseria web page.

Each genome assembly was queried against the pubMLST.org/neisseria sequence definition database to identify the loci within individual genomes. As identified, loci were automatically annotated with a *Neisseria* (NEIS) locus identifier and an allele number assigned if an identical reference sequence was present in the database. Novel alleles with ≥ 98% sequence identity to a reference sequence in the database, were automatically assigned new reference
<caption>**TABLE 1**</caption>

Clarification of capsule type of European Meningococcal Strain Collection isolates using whole genome analysis, July 2011– June 2012 (n = 32 isolates)

<table><thead><tr><th>Capsule operon region A and C genes</th><th>Number of isolates</th><th>Serogroup result</th></tr></thead><tbody><tr><td>C</td><td>5</td><td rowspan="5">B</td></tr><tr><td>E</td><td>1</td></tr><tr><td>Cnl</td><td>1</td></tr><tr><td>Mosaic B(+E)</td><td>1</td></tr><tr><td>Incomplete assembly</td><td>1</td></tr><tr><td>Y</td><td>1</td><td rowspan="4">C</td></tr><tr><td rowspan="4">B</td><td>2</td></tr><tr><td>1</td></tr><tr><td>3</td></tr><tr><td>5</td></tr><tr><td>Cnl</td><td>2</td><td rowspan="6">ND</td></tr><tr><td>E</td><td>1</td></tr><tr><td>W</td><td>2</td></tr><tr><td>W/Y</td><td>1</td></tr><tr><td>X</td><td>2</td></tr><tr><td>Y</td><td>1</td></tr><tr><td></td><td>1</td><td>W</td></tr></tbody></table>

Cnl: capsule null; mosaic B (+ E): complete B operon plus partial E found; ND: not done; NG: non-groupable.

Assignment of capsule type using gene-by-gene analysis of capsule operon in each genome where serogroup testing did not correspond or was not done; including one mosaic capsule operons (IBD-1081, B with partial E). Two genomes with non-groupable serotype results (IBD-407, B; IBD-409, B) were found to be genetically phase variable off.

numbers. This process was repeated until all NEIS loci had been curated and as many alleles as possible had been defined and assigned in each genome assembly. The new alleles with <98% sequence identity were manually checked and curated before assigning a new reference number. The annotated loci included those defining antigens, genes known to effect antibiotic sensitivity, the capsule operon genes, and the 1,605 loci belonging to the N. meningitidis core genome MLST (cgMLST) scheme (version 1.0) [9,13,14]. Alleles with a known antibiotic sensitivity range were previously submitted to the pubMLST.org/neisseria sequence definition database and EMSC-GL alleles with the same sequence mutation pattern were inferred to convey the same antibiotic sensitivity range [15-17].

## Genomic analysis and population structure

Using the 'sequence bin breakdown' tool, the assembly statistics and genome coverage were assessed for each genome assembly. The BIGSdb Genome Comparator Tool was used to generate allele-based pairwise comparison distance matrices that were visualised as Neighbour-Net phylogenies. Cramer's V coefficient was used to measure the association of clonal complex (cc) with Bexsero Antigen Sequence Typing (BAST) and

individual vaccine antigens, using the 'cramersV' function in the 'lsr' package in R 3.3.0 [18]. The coefficient ranges from zero to one, with zero indicating no association and one indicating complete association.

## Curatorial support and database operation

The data confidentiality and security: no personal identification is recorded in the database and isolates are identified by an anonymised sample name and an automatically assigned isolate accession number. Each reference laboratory is provided immediate access to their submitted data once it has been received and assessed by a curator. The pubMLST database and its curators are currently supported by The Wellcome Trust (grant 104992) and the UK Department of Health (grant H2R00080) respectively.

## Results

### Assembly and genetic coverage

The isolates (n=799) represented IMD cases from 16 EU countries during the epidemiological year July 2011–June 2012. The genome assemblies contained an average of 466 contigs, with an average length of 5,211 nt per contig, and an average genome length of 2,143,632 nt. The larger number of contigs and unassembled genes compared to other studies was due to fragmentation of the DNA during extraction process [13,19]. The pubMLST.org/neisseria database (accessed 06/2016) contained 2,429 curated loci, of which 1,605 genes (66.1%) are considered to be a part of the core genome (cgMLST, v1.0) [13]. The curated loci are composed of coding sequences, typing fragments, antigenic peptides, intergenic promotor regions, and other sequence regions of interest. Each draft genome contained between 1,724 and 1,933 identifiable genes with an average of 1,831 per genome. The average genetic content for a *N. meningitidis* genome, based on finished genomes (NZ-05/33, FAM18, MC58, H44/76, Z2491, F2136 and Mo1-240196), was 1,967 genes.

### Capsule

Isolates had been tested for capsule phenotype by the submitting laboratory: 65.8% were serogroup B (n=526); 20.8% were serogroup C (n=166); 7.5% were serogroup Y (n=60); 4.0% were serogroup A, W, or W/Y (n=1; 30; 1); and 1.9% were non-groupable (NG) or had no serogroup data recorded (n=14; 1). Analysis of capsule operon regions A and C from the WGS data: (i) genotypically confirmed the phenotypic serogrouping of 767/799 isolates (96.0%); (ii) defined the genogroup of 14 NG isolates and of one isolate without serological data; (iii) corrected the capsule designation of 14 (1.8%) isolates; and (iv) partially confirmed one B and one C capsule and found one mosaic capsule operon (Table 1).

The isolate without serological typing information contained capsule type B genes. Two genomes had a poorly assembled region A and the gene that determines capsule type was incompletely assembled; however, the
<caption>**TABLE 2**</caption>

Putative vaccine antigen coverages of invasive meningococcal disease isolates, based on sequence matching using the Bexsero antigen sequence type (BAST) scheme and assessing presence of the fHbp peptide components in the bivalent rLP4086 (Trumemba), Europe, July 2011–June 2012 (n=799 isolates)

<table><thead><tr><th>European country</th><th>Number of isolates in dataset</th><th>NadA peptide 8</th><th>NHBA peptide 2</th><th>FHbp peptide 1</th><th>PorA VR1 peptide 7-2</th><th>PorA VR2 peptide 4</th><th>BAST total (%)</th><th>FHbp peptide 45</th><th>FHbp peptide 55</th><th>Bivalent total (%)</th></tr></thead><tbody><tr><td>Austria</td><td>30</td><td>0</td><td>1</td><td>9</td><td>7</td><td>1</td><td>18 (60)</td><td>1</td><td>0</td><td>1 (3.3)</td></tr><tr><td>Cyprus</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0 (0.0)</td><td>0</td><td>0</td><td>0 (0.0)</td></tr><tr><td>Czech Republic</td><td>36</td><td>0</td><td>1</td><td>2</td><td>3</td><td>0</td><td>6 (16.7)</td><td>0</td><td>0</td><td>0 (0.0)</td></tr><tr><td>Denmark</td><td>19</td><td>0</td><td>0</td><td>6</td><td>3</td><td>0</td><td>9 (47.4)</td><td>1</td><td>0</td><td>1 (5.3)</td></tr><tr><td>Finland</td><td>34</td><td>0</td><td>10</td><td>2</td><td>11</td><td>10</td><td>33 (97.1)</td><td>1</td><td>0</td><td>1 (2.9)</td></tr><tr><td>France</td><td>42</td><td>0</td><td>4</td><td>6</td><td>6</td><td>6</td><td>22 (52.4)</td><td>4</td><td>0</td><td>4 (9.5)</td></tr><tr><td>Germany</td><td>245</td><td>1</td><td>41</td><td>26</td><td>45</td><td>36</td><td>149 (60.8)</td><td>3</td><td>0</td><td>3 (1.2)</td></tr><tr><td>Greece</td><td>17</td><td>0</td><td>0</td><td>6</td><td>0</td><td>0</td><td>6 (35.3)</td><td>0</td><td>0</td><td>0 (0.0)</td></tr><tr><td>Hungary</td><td>15</td><td>0</td><td>1</td><td>1</td><td>1</td><td>0</td><td>3 (20.0)</td><td>0</td><td>0</td><td>0 (0.0)</td></tr><tr><td>Latvia</td><td>3</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0 (0.0)</td><td>0</td><td>0</td><td>0 (0.0)</td></tr><tr><td>Norway</td><td>27</td><td>0</td><td>3</td><td>2</td><td>5</td><td>3</td><td>13 (48.2)</td><td>0</td><td>0</td><td>0 (0.0)</td></tr><tr><td>Poland</td><td>181</td><td>0</td><td>14</td><td>26</td><td>15</td><td>15</td><td>70 (38.7)</td><td>7</td><td>0</td><td>7 (3.9)</td></tr><tr><td>Romania</td><td>14</td><td>8</td><td>0</td><td>0</td><td>0</td><td>0</td><td>8 (57.1)</td><td>0</td><td>0</td><td>0 (0.0)</td></tr><tr><td>Slovakia</td><td>22</td><td>0</td><td>0</td><td>5</td><td>2</td><td>1</td><td>8 (36.4)</td><td>2</td><td>0</td><td>2 (9.1)</td></tr><tr><td>Slovenia</td><td>12</td><td>0</td><td>4</td><td>4</td><td>1</td><td>1</td><td>10 (83.3)</td><td>0</td><td>0</td><td>0 (0.0)</td></tr><tr><td>The Netherlands</td><td>101</td><td>3</td><td>16</td><td>18</td><td>17</td><td>10</td><td>64 (63.4)</td><td>8</td><td>0</td><td>8 (7.9)</td></tr><tr><td>All countries</td><td>799 (100)</td><td>12 (1.5)</td><td>95 (11.9)</td><td>113 (14.1)</td><td>116 (14.5)</td><td>83 (10.4)</td><td>419 (52.4)</td><td>27 (3.4)</td><td>0 (0.0)</td><td>27 (3.4)</td></tr><tr><td>Total (%)</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

The Bexsero vaccine is composed of four different peptide components from four different clonal complexes and is represented by the peptide profile BAST-1 (fHbp 1, NHBA 2, NadA 8, PorA VR1 7-2, PorA VR2 4). The Trumemba vaccine is composed of two fHbp peptides, variants 45 and 55.

assembled sequence indicated that genetically one isolate was B and one isolate was C.

Of the serogroup records that did not match capsule genotype: three isolates assigned to serogroup C were capsule type B (n=2) and capsule type Y (n=1); three serogroup Y were capsule type B; seven isolates assigned to serogroup B were capsule type C (n=5), capsule type E (n=1), and one isolate did not contain the capsular regions A or C and was confirmed to be a capsule null (cnl); and one typed as serogroup W was genetically a capsule type Y.

The capsule types determined for the NG isolates included: five capsule type B; two capsule type W; two capsule type X; one capsule type E; two cnl; one capsule Y; and one with a serine residue at codon 310 of the polymerase gene, which would result in a mixed sialic acid capsule phenotype for W and Y [20]. Two of the NG isolates contained capsule B genes but were genetically phase variable off.

The isolate with a mosaic capsule contained capsule biosynthesis genes from serogroup B in addition to a partial region A for capsule type E. There was no report of secondary capsule expression nor inconclusive serogroup results reported for the isolate.

## Clonal complex

MLST was determined from the WGS data. For 402/799 isolates (50.3%) MLST data obtained by conventional means were also available. The WGS data confirmed the sequence type (ST) and cc for 377/402 (93.8%) of these isolates and determined the ST and cc for 395 of the remaining 397 (99.5%) isolates.

Concerning the 25 isolates among the 402, where the conventional MLST and the WGS data did not concur: (i) for two the *gdh* MLST locus was incomplete but the remaining six MLST alleles were identified: one isolate had a profile consistent with the ST-269cc, but was nonetheless not assigned to any specific ST, while the other was unique and did not match any currently defined ST; (ii) and for 23 isolates, the WGS-derived MLST differed from that previously reported. Most (17/23) differed at a single MLST locus (including 10 genomes differed at *pgm*, 4 at *fumC*, 2 at *abcZ*, and 1 at the *aroE* locus), three differed at two loci (*pgm* and *abcZ* or *fumC* or *gdh*, respectively), and three isolates differed at all seven loci. The isolates where the ST was one or two loci different were still assigned to the same cc, the three isolates with seven loci differences were each reassigned to a different cc. These results were consistent with the downstream effects of splitting
TABLE 3

<table><caption>Evaluation of putative antibacterial susceptibility of invasive meningococcal disease isolates by sequence comparison, Europe, July 2011–June 2012 (n=799 isolates)*</caption><thead><tr><th rowspan="2">Country</th><th colspan="4">gyrA gene: ciprofloxacin<sup>a</sup></th><th colspan="2">rpoB gene: rifampicin<sup>b</sup></th><th colspan="4">penA gene: penicillin<sup>c</sup></th></tr><tr><th>tS<sup>d</sup></th><th>tS(I)<sup>d</sup></th><th>eS<sup>e</sup></th><th>UKN</th><th>tS<sup>d</sup></th><th>UKN</th><th>tS<sup>d</sup></th><th>tS(I)<sup>d</sup></th><th>tI<sup>d</sup></th><th>tS(I)<sup>d</sup></th><th>UKN</th></tr></thead><tbody><tr><td>Austria</td><td>16</td><td>13</td><td>0</td><td>1</td><td>30</td><td>0</td><td>0</td><td>3</td><td>15</td><td>1</td><td>11</td></tr><tr><td>Cyprus</td><td>1</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>Czech Republic</td><td>7</td><td>29</td><td>0</td><td>0</td><td>34</td><td>2</td><td>1</td><td>8</td><td>15</td><td>0</td><td>12</td></tr><tr><td>Denmark</td><td>6</td><td>13</td><td>0</td><td>0</td><td>19</td><td>0</td><td>0</td><td>0</td><td>17</td><td>0</td><td>2</td></tr><tr><td>Finland</td><td>10</td><td>12</td><td>12</td><td>0</td><td>31</td><td>3</td><td>0</td><td>6</td><td>15</td><td>1</td><td>12</td></tr><tr><td>France</td><td>14</td><td>26</td><td>0</td><td>2</td><td>41</td><td>1</td><td>0</td><td>9</td><td>20</td><td>1</td><td>12</td></tr><tr><td>Germany</td><td>84</td><td>154</td><td>0</td><td>7</td><td>240</td><td>5</td><td>0</td><td>25</td><td>141</td><td>7</td><td>72</td></tr><tr><td>Greece</td><td>6</td><td>10</td><td>0</td><td>1</td><td>17</td><td>0</td><td>0</td><td>2</td><td>13</td><td>0</td><td>2</td></tr><tr><td>Hungary</td><td>5</td><td>9</td><td>0</td><td>1</td><td>15</td><td>0</td><td>0</td><td>1</td><td>10</td><td>1</td><td>3</td></tr><tr><td>Latvia</td><td>2</td><td>1</td><td>0</td><td>0</td><td>3</td><td>0</td><td>0</td><td>0</td><td>3</td><td>0</td><td>0</td></tr><tr><td>Norway</td><td>11</td><td>9</td><td>7</td><td>0</td><td>26</td><td>1</td><td>0</td><td>6</td><td>12</td><td>0</td><td>9</td></tr><tr><td>Poland</td><td>41</td><td>86</td><td>52</td><td>2</td><td>175</td><td>6</td><td>1</td><td>36</td><td>102</td><td>5</td><td>37</td></tr><tr><td>Romania</td><td>5</td><td>9</td><td>0</td><td>0</td><td>14</td><td>0</td><td>0</td><td>12</td><td>2</td><td>0</td><td>0</td></tr><tr><td>Slovakia</td><td>10</td><td>10</td><td>0</td><td>2</td><td>21</td><td>1</td><td>0</td><td>3</td><td>11</td><td>1</td><td>7</td></tr><tr><td>Slovenia</td><td>7</td><td>5</td><td>0</td><td>0</td><td>12</td><td>0</td><td>0</td><td>1</td><td>6</td><td>0</td><td>5</td></tr><tr><td>The Netherlands</td><td>26</td><td>51</td><td>22</td><td>2</td><td>98</td><td>3</td><td>1</td><td>20</td><td>52</td><td>5</td><td>23</td></tr><tr><td>Sensitive</td><td>251</td><td>0</td><td>93</td><td>0</td><td>777</td><td>0</td><td>3</td><td>0</td><td>0</td><td>22</td><td>0</td></tr><tr><td>Intermediate resistance</td><td>0</td><td>437</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>132</td><td>434</td><td>0</td><td>0</td></tr><tr><td>Unknown</td><td>0</td><td>0</td><td>0</td><td>18</td><td>0</td><td>22</td><td>0</td><td>0</td><td>0</td><td>0</td><td>208</td></tr></tbody></table>

I: intermediate; MIC: minimum inhibitory concentration; S: susceptible; S(I): allele for which additional sensitive or intermediate MIC ranges (µg/mL) have been recorded; UKN: unknown.

Isolates in this study were not directly tested for a MIC but where a genome contained a known allele with a MIC assay result the isolate was interpreted as having the same antibiotic sensitivity value.

a Ciprofloxacin MIC range in µg/mL: ≤0.03 (susceptible); >0.03 –<0.12 (intermediate); >0.12 (resistant).

b Rifampicin MIC range in µg/mL: ≤1 (susceptible); >1 (resistant).

c Penicillin MIC range in µg/mL: ≤0.06 (susceptible); >0.06 – 1 (intermediate); >1 (resistant).

d t: tested MIC (the European meningococcal strain collection genome library isolate allele is an exact sequence match to an allele previously submitted to the database with a MIC tested value).

e e: expected MIC (the European meningococcal strain collection genome library isolate allele has key known mutations allowing to infer antibacterial susceptibility).

samples in seven-locus PCR-based MLST and mislabeling of samples [13].

In total, 795 isolates had an assigned ST, among which 108 with no cc association. Two-hundred-ninety STs were defined (795/799; 99.5%) and twenty-nine cc were identified (687/799; 86.0%). The three most common cc, cc41/44 (n=177 isolates), cc32 (n=132), and cc11 (n=92), represented 50.2% of the 799 isolates. Eight additional cc had at least 10 representatives each and comprised 29.0% (232) of the isolates (Figure 1). The remaining 18 cc were represented by one to nine isolates each. The 108 isolates (13.5%) without a cc association were less than 1% of the 290 STs defined. Overall the most prominent STs were ST-5238 from Poland (n=6) and Germany (n=1) and ST-10958 from Romania (n=8).

## Antigenic typing and vaccine coverage

The antigenic fine-type data (PorA VR1:PorA VR2:FetA VR) were assessed for each of the isolates' genome assemblies. The finetype antigen profile was complete

in 778 of the isolates (97.4%). The porA gene (NEIS1364) was absent in 10 genomes (1.3%), while four genomes had a deleted PorA VR1 loop, two genomes had a deleted PorA VR2 loop, and one genome had an incompletely assembled NEIS1364 gene and its PorA VR1 loop could not be assigned. Additionally, six genomes (0.8%) were lacking FetA VR assignment: five were missing the NEIS1963 (fetA) gene and one had a non-coding protein sequence. The most prevalent antigen fine-types in the dataset (n≥10) included: 7-2:4:F1-5 (8.0%); 7:16:F3-3 (5.2%); 5:2:F3-3 (4.7%); 22:14:F5-5 (3.7%); 18-1:3:F3-9 (3.1%); 5-1:10-8:F3-6 (3.0%); and 5-2:10-1:F4-1 (3.0%).

The meningococcal genomes were characterised with the BAST scheme. The 4CMenB (Bexsero) vaccine contains four different peptide components and is represented by the peptide profile BAST-1 (fHbp 1, NHBA 2, NadA 8, PorA VR1 7-2, PorA VR2 4) [9]. There were no occurrences of BAST-1 in the EMSC-GL. BAST antigens were diverse and 687 genomes (86.0%) were assigned to one of 332 BASTs. The genomes without a BAST
assigned (n=112; 14.0%) had an incomplete sequence assembly in the region containing the NEIS2109 (nhba) or the NEIS1969 (nadA) gene. The most common 12 BASTs were present at a frequency of at least 1.0% (n=8) and an additional 32 BASTs were represented by at least three to seven genomes each (Figure 2).

Among the total 332 BASTs, the remaining BASTs included 50 profiles (15.1%) with two representative genomes each and 238 profiles (71.7%) that occurred once. Overall the two most common profiles were BAST-223 (FHbp 14, NHBA 2, NadA 0, PorA VR1 7-2, PorA VR2 4) and BAST-4 (FHbp 1, NHBA 3, NadA 1, PorA VR1 7, PorA VR2 16). The gene encoding nadA was absent in 64.2% (n=513) of 799 genomes and 229 of the BASTs. The BAST of each genome was assessed for direct matches to the BAST-1 components for prediction of potential vaccine coverage. Predicted coverage, based on sequence identity, was most frequently associated with the PorA VR1 7-2 (n=116; 14.5%) and fHbp 1 (n=113; 14.1%) vaccine components (Table 2).

The genomes were also assessed for exact peptide matches for the bivalent rLP2086 (Trumemba) vaccine composed of fHbp variants 45 and 55. The fHbp variant 55 was not observed, while variant 45 was present in 27 (3.4%) genomes from eight countries: the Netherlands (n=8), Poland (n=7), France (n=4), Germany (n=3), Slovakia (n=2), and one each in Austria, Denmark, and Finland. There was no overlap of coverage between the vaccines, 52.4% (n=419) of the genomes had exact matches to BAST component and an additional 3.4% (n=27) of the genomes contained the fHbp variant 45.

## Antimicrobial resistance

The pubMLST.org/neisseria database contained some phenotypic antimicrobial data in the form of minimum inhibitory concentration (MIC) ranges. These data were limited to a subset of isolates used for targeted gene sequencing to define susceptibility for ciprofloxacin (NEIS1320, *gyrA* gene), penicillin (NEIS1753, *penA* gene), and rifampicin (NEIS0123, *rpoB* gene) [15,16]. Alleles without MIC data were compared with those with a MIC value to determine if the allele sequence contained known mutations associated with a phenotypic MIC range (Table 3).

Forty-one *rpoB* alleles were present in the dataset, 23 of which were associated with a susceptible MIC (≤1µg/mL). Most genomes (n=501; 62.7%) contained one of four *rpoB* alleles: 18, 2, 4, or 34, and an additional eight alleles were present in at least 1% of the genomes (n=256; 32.0%). The next most prevalent, allele 184 (n=5; 0.63%), had no MIC associated testing data and was found in Germany (n=3) and Poland (n=2). There were 18 additional alleles, each present in the dataset once, with no MIC associated testing data.

There were 141 genomes (17.7%) with an incompletely assembled *penA* gene and no allele assignment and 44 genomes (5.5%) did not have the *penA* gene. There

were 67 *penA* alleles in the dataset, 13 of which had MIC values that were generally sensitive, ≤0.06µg/mL (two alleles were always sensitive and 11 alleles were sensitive but had at least one isolate record where the MIC value was intermediate, >0.06 to 1µg/mL). Thirty-eight *penA* alleles were associated with MIC values of intermediate sensitivity (≥0.06 to 1µg/mL). There were 16 additional alleles with no MIC associated testing data, 15 of which were found once in the dataset and one allele was found eight times, from Germany (n=4), Poland (n=3) and Hungary (n=1).

There were a total of 19 *gyrA* alleles and 642 genomes (80.4%) contained allele 4, 2, or 3. Eight genomes (1.0%) had an incompletely assembled *gyrA* gene. Five additional alleles were found in at least 1% of the isolate genomes and represented an additional 16.5% of the genomes (n=132). All eight alleles had an associated susceptible MIC test result (≤0.03µg/mL) and half of the alleles had at least one isolate record where the MIC result was intermediate (≥0.03 but <0.12µg/mL). Only three of the eleven remaining alleles had MIC associated data and sequence comparison of the eight alleles with known MIC data was not discriminatory.

## Population structure

Meningococcal lineage nomenclature was used as the fundamental unit of analysis for the genetic relationships at the genome level [13]. This is based on ribosomal sequence types (rST) derived from the 53 ribosomal protein genes (rMLST) and supported with cgMLST analysis [13,21]. There were 788 complete rSTs, and 11 rSTs unassigned due to the incomplete assembly of one (n=6), two (n=3), or three (n=2) of the 53 loci. A total of 13 lineages, corresponding to a known cc were identified: lineage 3 (cc41/44); lineage 8 (cc8); lineage 11 (cc11); lineage 22 (cc22); lineage 23 (cc23); lineage 2 (cc269); lineage 35 (cc35); lineage 5 (cc32); lineage 25 (cc162); lineage 13 (cc213); lineage 18 (cc18); lineage 6 (cc60); and lineage 20 (cc103) (Figure 3).

Additional groups were evident, including a cluster of genomes from Poland (n=7), Germany (n=7), and Hungary (n=2), containing eight STs not associated with a cc. Members of this new lineage, while sharing up to five of seven MLST alleles with members of cc174 (lineage 14), were distinct using cgMLST (data not shown). A number of lineages formed larger groups or clades including: lineages 11 and 8; and lineages 2 and 35 with lineage 5. Lineage 11 and lineage 23 contained distinct sub-lineages, (11.1 and 11.2; 23.1 and 23.2) as previously described [19]. The vaccine antigen variants and combinations were associated with particular ccs and lineages, with Cramer's V values for cc association ranging from 0.458 (NadA) to 0.703 (NHBA) and 0.94 (BAST).

## Lessons learnt

*Neisseria* genomes contain repeat regions, homopolymeric tracts, and paralogous loci [22,23]; all of which may affect sequencing fidelity, assembly synteny, and
genome annotation. In addition, whole chromosome
extraction and library preparation methods will have
downstream effects on the quality of the assembled
genome data. A pilot run comprising a minimum of 96
samples is therefore recommended so that the impact
of such issues can be ameliorated before protocols are
implemented on larger studies, which can be costly.
The curation of genomic data within the pubMLST.org/
neisseria database was developed using a population
annotation approach which places allelic diversity at
the centre of analysis [12,13], and therefore required
a sequence definition database of proportional size
and variability to the organism's genetic variation. The
allelic annotation against a locus specific definition,
where accurate start and stop codons were described
while capturing the population variability in a meaning-
ful and usable way, requires dedicated manual over-
sight of curation. Sustaining the computational and
comparative value of a genome library will require reg-
ular manual checks, data cleaning, and periodic locus
definition updates over the life of the database.

Discussion

The sporadic nature of IMD, combined with the unpre-
dictable occurrence of epidemics and the diversity of
meningococci, and especially of meningococcal anti-
gens, makes accurate surveillance an essential ele-
ment of effective public health measures [24]. The
well-documented propensity of particular hyperinva-
sive lineages to cause IMD outbreaks across continents
makes international surveillance of particular impor-
tance [7]. Within Europe, there has been coordinated
surveillance of IMD since the early 1990s [25]. During
this time, meningococcal isolate characterisation was
greatly improved by the introduction of molecular
methods [26]. The latest enhancement, the use of WGS
data was first deployed for routine and comprehensive
IMD surveillance on a national scale in the UK, with the
establishment of the Meningitis Research Foundation
Meningococcus Genome Library (MRF-MGL) from 2010
[19]. The WGS approach provides high-resolution index-
ing of bacterial diversity potentially replacing phe-
notypic methods and enabling multiple features of a
clinical isolate to be characterised in a single test [13].
With WGS technology becoming increasingly accessible
and affordable the EMSC-GL becomes a practical tool
for IMD surveillance, complementing and enhancing
the information collected by the European Surveillance
System (TESSy) database at ECDC and contributing to
studies of population biology, evolution, and the con-
sequences of vaccine implementation [27,28].

To maximise international public health impact, the
EMSC-GL is hosted within the pubMLST.org/neisseria
web-accessible database, which also hosts the MRF-
MGL, and allows members to submit genome data and
the isolates associated provenance information. This
provides publicly-accessible sequence data, described
with a uniform nomenclature and linked to relevant
provenance and phenotype information. As of early
2017, the continually-expanding pubMLST.org/neisseria

database comprised over 2,400 Neisseria loci, including
38 sequence-based typing loci such as antigen genes,
seven locus MLST, and cgMLST loci. This permits the
rapid genetic characterisation and comparison of >93%
of the content of future meningococci genomes added
to the genome library with minimal additional curation.
Within pubMLST.org/neisseria these loci are grouped
into different schemes: currently 40 such curated
schemes are publicly available, including those for
metabolic pathways, plasmid and phage components,
typing profiles, and genomic regions [13]. The WGS
approach provides a number of practical advantages
in addition to its inherent high-resolution and com-
parability of nomenclature [13,29]. It is a single test
that eliminates the multiple testing of a single speci-
men, which can introduce errors [30]. For example, the
incorrect grouping of alleles into a conventional MLST
profile can affect the investigation of an outbreak of
invasive meningococcal disease [31]. While errors can
occur during de novo WGS assembly, cgMLST analysis
ameliorates these problems as, in a high-quality draft
genome, such problems affect a minority of loci which
are unlikely to materially alter an analysis [13].

Given its importance in invasion and as a vaccine anti-
gen the meningococcal capsule, which defines the
serogroup, is one of the most important isolate char-
acteristics required for disease surveillance. The cap-
sule encoding genome regions are analysable by WGS,
resolving anomalies that can arise from phenotypic
data and enabling the characterisation of meningo-
cocci that do not express capsule. The EMSC-GL con-
tained a small number of isolates with indeterminate
W/Y expression and isolates containing genes of more
than one capsule type [19]. The difficulty in determin-
ing the phenotype of phenotypically W/Y isolates was
the result of a rare serine residue at codon position 310
of the NEIS2162 (csw) and NEIS2163 (csy) genes and,
while there is no evidence of multiple capsule expres-
sion, horizontal genetic transfer of region A capsule
genes has been described for other isolates [14,32].
The presence of different capsule genes highlights
recombination potential among closely related menin-
gococci, which can lead to novel phenotypes of clinical
and epidemiological significance [14]. Finally, although
the production of a capsule is considered to be a viru-
lence factor, there have been reports of un-encapsu-
lated meningococci, one of which was present in the
EMSC genomes, causing IMD [33,34].

Analysis of cgMLST data defined 13 lineages, had good
ccs concordance as defined by conventional MLST, and
identified a new lineage. Vaccine antigens were also
rapidly extracted, including those used to assign BAST
types, demonstrating a strong association between
vaccine antigen variants and cc among the EMSC-GL.
This association was consistent with previous studies
in other countries such as Ireland, the Netherlands and
the UK, and is important when considering the likely
impact of protein-based vaccines because prevalent
lineages change over time likely causing changes in disease incidence [9,19,35].

The deposition of future meningococcal genome data in the EMSC-GL will permit the resolution of the epidemiology of IMD at the whole genome across the entire continent, and as long as the European Union (EU) maintains funding the infrastructure costs are minimal. The inclusion of additional isolates in close to real time will enable the rapid identification of outbreaks. The past 40 years have seen the spread of successive hyperinvasive meningococci throughout Europe and, while disease levels remain low in most European countries most of the time, epidemics, hyperendemics, and localised disease outbreaks occur regularly across the continent, including transnational outbreaks [19,26,36-38]. The advantage of WGS data is that it comprehensively characterises isolates for their genetic lineage, vaccine antigens, antimicrobial sensitivity, and other properties rapidly and effectively and is able to resolve outbreaks to a high level of discrimination [13,31]. These advantages greatly outweigh perceived high costs of the WGS approach. Future addition of genomes to the library will enable monitoring of secular changes in circulating meningococci over time and an assessment of the possible effects of the introduction of vaccines, such as the protein-based 4MenB (Bexsero) and LP2086 (Trumenba) vaccines on IMD epidemiology [9]. Effective meningococcal vaccine development and implementation depends on an understanding of the population biology and the effects of herd immunity in an endemic setting, which in turn relies on high-quality surveillance data and isolate characterisation [39].

## Conclusion

In conclusion, the detailed meningococcal isolate characterisation available thought the EMSC-GL will be essential for European IMD surveillance, leading to the improvement of control strategies and disease prevention across the continent. This model is generalisable to other bacterial pathogens, many examples can be found on the pubMLST.org website (https://pubmlst.org/databases/).

## *Erratum

In Table 3, the column headings were misaligned. This was corrected on 31 May 2018.

## Acknowledgements

The project was funded by the ECDC within the framework of the IBD-labnet program through the Specific Contract No3 ECD.4051 implementing Framework Contract No ECDC/2011/014.

Database curators are currently supported by the UK Department of Health (grant H2R00080).

This publication made use of the Neisseria Multi Locus Sequence Typing website (https://pubmlst.org/neisseria/) developed by Keith Jolley and sited at the University of Oxford (Jolley & Maiden 2010, BMC Bioinformatics,

11:595). The development of this site has been funded by the Wellcome Trust and European Union, Wellcome Trust Biomedical Resource Grant (104992).

The authors wish to thank Charlene Rodrigues for her work in assigning the BAST designations.

## Conflict of interest

None declared.

## Authors' contributions

SH, DPB, PK, StH, MT, MKT, HC, GT, TE, JG, AVDE, AS, MP, AV, MPa, DAC: contributed isolates to the project; DAC: coordinated the collection of the strain collection and provided the DNA for the isolates; HBB: annotated and curated the genomes, analysed the data, and wrote the manuscript; CB: coordinated the genome sequencing and calculated Cramer's V coefficient; MCJM and DAC designed the study and contributed to the manuscript content and editing.

## References

1. Jafri RZ, Ali A, Messonnier NE, Tevi-Benissan C, Durrheim D, Eskola J, et al. Global epidemiology of invasive meningococcal disease. Popul Health Metr. 2013;11(1):17. https://doi.org/10.1186/1478-7954-11-17 PMID: 24016339

2. Pace D, Pollard AJ. Meningococcal disease: clinical presentation and sequelae. Vaccine. 2012;30(Suppl 2):B3-9. https://doi.org/10.1016/j.vaccine.2011.12.062 PMID: 22607896

3. Finne J, Leinonen M, Mäkelä PH. Antigenic similarities between brain components and bacteria causing meningitis. Implications for vaccine development and pathogenesis. Lancet. 1983;2(8346):355-7. https://doi.org/10.1016/S0140-6736(83)90340-9 PMID: 6135869

4. Jódar L, Feavers IM, Salisbury D, Granoff DM. Development of vaccines against meningococcal disease. Lancet. 2002;359(9316):1499-508. https://doi.org/10.1016/S0140-6736(02)08416-7 PMID: 11988262

5. Holst J, Martin D, Arnold R, Huergo CC, Oster P, O'Hallahan J, et al. Properties and clinical performance of vaccines containing outer membrane vesicles from Neisseria meningitidis. Vaccine. 2009;27(Suppl 2):B3-12. https://doi.org/10.1016/j.vaccine.2009.04.071 PMID: 19481313

6. Pizza M, Rappuoli R. Neisseria meningitidis: pathogenesis and immunity. Curr Opin Microbiol. 2015;23:68-72. https://doi.org/10.1016/j.mib.2014.11.006 PMID: 25461575

7. Caugant DA, Maiden MC. Meningococcal carriage and disease-population biology and evolution. Vaccine. 2009;27(Suppl 2):B64-70. https://doi.org/10.1016/j.vaccine.2009.04.061 PMID: 19464092

8. Feavers IM, Maiden MCJ. Recent progress in the prevention of serogroup B meningococcal disease. Clin Vaccine Immunol. 2017;24(5):e00566-16. https://doi.org/10.1128/CVI.00566-16 PMID: 28356256

9. Brehony C, Rodrigues CMC, Borrow R, Smith A, Cunney R, Moxon ER, et al. Distribution of Bexsero® Antigen Sequence Types (BASTs) in invasive meningococcal disease isolates: Implications for immunisation. Vaccine. 2016;34(39):4690-7. https://doi.org/10.1016/j.vaccine.2016.08.015 PMID: 27521232

10. Jolley KA, Brehony C, Maiden MC. Molecular typing of meningococci: recommendations for target choice and nomenclature. FEMS Microbiol Rev. 2007;31(1):89-96. https://doi.org/10.1111/j.1574-6976.2006.00057.x PMID: 17168996

11. Jolley KA, Maiden MC. Automated extraction of typing information for bacterial pathogens from whole genome sequence data: Neisseria meningitidis as an exemplar. Euro Surveill. 2013;18(4):20379. https://doi.org/10.2807/ese.18.04.20379-en PMID: 23369391

12. Jolley KA, Maiden MC. BIGSdb: Scalable analysis of bacterial genome variation at the population level. BMC Bioinformatics. 2010;11(1):595. https://doi.org/10.1186/1471-2105-11-595 PMID: 21143983

13. Bratcher HB, Corton C, Jolley KA, Parkhill J, Maiden MC. A gene-by-gene population genomics platform: de novo assembly, annotation and genealogical analysis of 108
representative Neisseria meningitidis genomes. BMC
Genomics. 2014;15(1):1138. https://doi.org/10.1186/1471-2164-
15-1138 PMID: 25523208

14. Harrison OB, Claus H, Jiang Y, Bennett JS, Bratcher HB, Jolley KA, et al. Description and nomenclature of Neisseria meningitidis capsule locus. Emerg Infect Dis. 2013;19(4):566-73. https://doi.org/10.3201/eid1904.111799 PMID: 23628376

15. Taha MK, Hedberg ST, Szatanik M, Hong E, Ruckly C, Abad R, et al. Multicenter study for defining the breakpoint for rifampin resistance in Neisseria meningitidis by rpoB sequencing. Antimicrob Agents Chemother. 2010;54(9):3651-8. https://doi.org/10.1128/AAC.00315-10 PMID: 20606072

16. Taha MK, Vázquez JA, Hong E, Bennett DE, Bertrand S, Bukovski S, et al. Target gene sequencing to characterize the penicillin G susceptibility of Neisseria meningitidis. Antimicrob Agents Chemother. 2007;51(8):2784-92. https://doi.org/10.1128/AAC.00412-07 PMID: 17517841

17. Hong E, Thulin Hedberg S, Abad R, Fazio C, Enríquez R, Deghmane AE, et al. Target gene sequencing to define the susceptibility of Neisseria meningitidis to ciprofloxacin. Antimicrob Agents Chemother. 2013;57(4):1961-4. https://doi.org/10.1128/AAC.02184-12 PMID: 23357770

18. R Core Team RD. R: A Language and Environment for Statistical Computing. In: Vienna, Austria: R Foundation for Statistical Computing; 2006.

19. Hill DMC, Lucidarme J, Gray SJ, Newbold LS, Ure R, Brehony C, et al. Genomic epidemiology of age-associated meningococcal lineages in national surveillance: an observational cohort study. Lancet Infect Dis. 2015;15(12):1420-8. https://doi.org/10.1016/S1473-3099(15)00267-4 PMID: 26515523

20. Claus H, Borrow R, Achtman M, Morelli G, Kantelberg C, Longworth E, et al. Genetics of capsule O-acetylation in serogroup C, W-135 and Y meningococci. Mol Microbiol. 2004;51(1):227-39. https://doi.org/10.1046/j.1365-2958.2003.03819.x PMID: 14651624

21. Jolley KA, Bliss CM, Bennett JS, Bratcher HB, Brehony C, Colles FM, et al. Ribosomal multilocus sequence typing: universal characterization of bacteria from domain to strain. Microbiology. 2012;158(Pt 4):1005-15. https://doi.org/10.1099/mic.0.055459-0 PMID: 22282518

22. Bentley SD, Vernikos GS, Snyder LA, Churcher C, Arrowsmith C, Chillingworth T, et al. Meningococcal genetic variation mechanisms viewed through comparative analysis of serogroup C strain FAM18. PLoS Genet. 2007;3(2):e23. https://doi.org/10.1371/journal.pgen.0030023 PMID: 17305430

23. Parkhill J, Achtman M, James KD, Bentley SD, Churcher C, Klee SR, et al. Complete DNA sequence of a serogroup A strain of Neisseria meningitidis Z2491. Nature. 2000;404(6777):502-6. https://doi.org/10.1038/35006655 PMID: 10761919

24. German RR, Lee LM, Horan JM, Milstein RL, Pertowski CA, Waller MNGuidelines Working Group Centers for Disease Control and Prevention (CDC). Updated guidelines for evaluating public health surveillance systems: recommendations from the Guidelines Working Group. MMWR Recomm Rep. 2001;50(RR-13):1-35, quiz CE1-7. PMID: 18634202

25. Trotter CL, Chandra M, Cano R, Larrauri A, Ramsay ME, Brehony C, et al. A surveillance network for meningococcal disease in Europe. FEMS Microbiol Rev. 2007;31(1):27-36. https://doi.org/10.1111/j.1574-6976.2006.00060.x PMID: 17168995

26. Brehony C, Jolley KA, Maiden MC. Multilocus sequence typing for global surveillance of meningococcal disease. FEMS Microbiol Rev. 2007;31(1):15-26. https://doi.org/10.1111/j.1574-6976.2006.00056.x PMID: 17168997

27. Reuter S, Ellington MJ, Cartwright EJ, Köser CU, Török ME, Gouliouris T, et al. Rapid bacterial whole-genome sequencing to enhance diagnostic and public health microbiology. JAMA Intern Med. 2013;173(15):1397-404. https://doi.org/10.1001/jamainternmed.2013.7734 PMID: 23857503

28. Croucher NJ, Finkelstein JA, Pelton SI, Mitchell PK, Lee GM, Parkhill J, et al. Population genomics of post-vaccine changes in pneumococcal epidemiology. Nat Genet. 2013;45(6):656-63. https://doi.org/10.1038/ng.2625 PMID: 23644493

29. Didelot X, Urwin R, Maiden MC, Falush D. Genealogical typing of Neisseria meningitidis. Microbiology. 2009;155(Pt 10):3176-86. https://doi.org/10.1099/mic.0.031534-0 PMID: 19643763

30. Bratcher HB, Bennett JS, Maiden MCJ. Evolutionary and genomic insights into meningococcal biology. Future Microbiol. 2012;7(7):873-85. https://doi.org/10.2217/fmb.12.62 PMID: 22827308

31. Jolley KA, Hill DM, Bratcher HB, Harrison OB, Feavers IM, Parkhill J, et al. Resolution of a meningococcal disease outbreak from whole-genome sequence data with rapid Web-based analysis methods. J Clin Microbiol. 2012;50(9):3046-53. https://doi.org/10.1128/JCM.01312-12 PMID: 22785191

32. Tsang RS, Tsai CM, Henderson AM, Tyler S, Law DK, Zollinger W, et al. Immunochemical studies and genetic background of two Neisseria meningitidis isolates expressing unusual capsule polysaccharide antigens with specificities of both serogroup Y and W135. Can J Microbiol. 2008;54(3):229-34. https://doi.org/10.1139/W07-132 PMID: 18388994

33. Vogel U, Claus H, von Müller L, Bunjes D, Elias J, Frosch M. Bacteremia in an immunocompromised patient caused by a commensal Neisseria meningitidis strain harboring the capsule null locus (cnl). J Clin Microbiol. 2004;42(7):2898-901. https://doi.org/10.1128/JCM.42.7.2898-2901.2004 PMID: 15243035

34. Findlow H, Vogel U, Mueller JE, Curry A, Njanpop-Lafourcade BM, Claus H, et al. Three cases of invasive meningococcal disease caused by a capsule null locus strain circulating among healthy carriers in Burkina Faso. J Infect Dis. 2007;195(7):1071-7. https://doi.org/10.1086/512084 PMID: 17330799

35. Bambini S, Piet J, Muzzi A, Keijzers W, Comandi S, De Tora L, et al. An analysis of the sequence variability of meningococcal fHbp, NadA and NHBA over a 50-year period in the Netherlands. PLoS One. 2013;8(5):e65043. https://doi.org/10.1371/journal.pone.0065043 PMID: 23717687

36. Harrison OB, Bray JE, Maiden MC, Caugant DA. Genomic analysis of the evolution and global spread of hyper-invasive meningococcal lineage 5. EBioMedicine. 2015;2(3):234-43. https://doi.org/10.1016/j.ebiom.2015.01.004 PMID: 25984558

37. Taha MK, Achtman M, Alonso JM, Greenwood B, Ramsay M, Fox A, et al. Serogroup W135 meningococcal disease in Hajj pilgrims. Lancet. 2000;356(9248):2159. https://doi.org/10.1016/S0140-6736(00)03502-9 PMID: 11191548

38. Lucidarme J, Hill DM, Bratcher HB, Gray SJ, du Plessis M, Tsang RS, et al. Genomic resolution of an aggressive, widespread, diverse and expanding meningococcal serogroup B, C and W lineage. J Infect. 2015;71(5):544-52.; Epub ahead of print. https://doi.org/10.1016/j.jinf.2015.07.007 PMID: 26226598

39. Trotter CL, Maiden MC. Meningococcal vaccines and herd immunity: lessons learned from serogroup C conjugate vaccination programs. Expert Rev Vaccines. 2009;8(7):851-61. https://doi.org/10.1586/erv.09.48 PMID: 19538112

License and copyright

This is an open-access article distributed under the terms of the Creative Commons Attribution (CC BY 4.0) Licence. You may share and adapt the material, but must give appropriate credit to the source, provide a link to the licence, and indicate if changes were made.

This article is copyright of the authors, 2018.