Characterization of invasive Neisseria meningitidis strains from Québec, Canada, during a period of increased serogroup B disease, 2009-2013: phenotyping and genotyping with special emphasis on the non-carbohydrate protein vaccine targets

Background The epidemiology of invasive meningococcal disease (IMD) in Québec, Canada, has been dominated in the past decade by a clone of serogroup B (MenB) Neisseria meningitidis defined by multi-locus sequence typing (MLST) as sequence type (ST)-269. With the licensure of a new MenB vaccine Bexsero (4CMenB) in Canada, this study characterized invasive N. meningitidis recovered in Québec from 2009 to 2013, with an objective to examine the diversity of the 4CMenB vaccine antigens. Isolates were serogrouped by antisera and genogrouped by PCR, and further typed by whole cell ELISA for serotype and serosubtype antigens. Clonal analysis was done by MLST. Isolates were genotyped by analysis of their 4CMenB vaccine antigen genes of PorA, factor H binding protein (fHbp), Neisserial Heparin Binding Antigen (NHBA), and Neisseria Adhesin A (NadA). Results Of the 263 IMD isolates analysed, 229, 16, 10, 7, and 1 belonged to MenB, MenY, MenW, MenC, and MenX, respectively. Of the 229 MenB, 159 (69.4 %) were typed as ST-269 clonal complex (CC); and they possessed a restricted number of three fHbp and five nhba gene alleles. Nine N. meningitidis isolates (eight MenB and one MenY) were found to possess at least one gene that encoded for an antigen that matched exactly with protein variants in the 4CMenB vaccine. Two MenB expressed PorA antigen P1.4 and possessed the nhba gene for peptide 2; four other MenB were predicted to have NHBA peptide 2; another two MenB were predicted to encode fHbp peptide 1.1; and a single MenY was found to have nadA gene for NadA peptide 8. In addition, another 172 isolates were found to possess genes for variant 1 fHbp peptides other than peptide 1.1 or NadA variant 1-2/3 peptides other than peptide 8; and therefore, may potentially be covered by 4CMenB. Conclusion The most prevalent clone of N. meningitidis in Quebec was ST-269 CC; and 96 % of the isolates in this CC were predicted to be covered by 4CMenB vaccine. Extensive genetic diversity was found in the other IMD isolates in Québec which might suggest a lower coverage by the vaccine when compared to the ST-269 MenB. Electronic supplementary material The online version of this article (doi:10.1186/s12866-015-0469-6) contains supplementary material, which is available to authorized users.


(Continued from previous page)
Conclusion: The most prevalent clone of N. meningitidis in Quebec was ST-269 CC; and 96 % of the isolates in this CC were predicted to be covered by 4CMenB vaccine. Extensive genetic diversity was found in the other IMD isolates in Québec which might suggest a lower coverage by the vaccine when compared to the ST-269 MenB.
Keywords: Invasive Neisseria meningitidis, 4CMenB protein vaccine targets Background Neisseria meningitidis is a Gram-negative diplococcus which infects only humans. On average 10 % of healthy individuals may carry this bacterium in their nasopharynx, and due to reasons not completely understood, the bacterium sometimes gains access into the blood stream and causes systemic diseases such as meningitis, septicemia, and occasionally localised infections in the joint, and heart [1]. Six of the 12 known serogroups (A, B, C, W, X, and Y) of meningococci are responsible for most of the invasive diseases which occur worldwide [2,3]. Although all 6 serogroups may cause outbreaks, serogroups A, B, C, and W have caused epidemics that may spread globally [4]. Case fatality rates associated with the different serogroups and in particular clonal groups may also vary [5,6].
Besides the capsule, typing of N. meningitidis is usually based on surface antigens or their genes, such as the PorB and PorA (corresponding to serotype and serosubtype antigens, respectively), as well as non-cell surface components such as housekeeping enzymes or their genes using the method of multilocus enzyme electrophoresis or multilocus sequence typing (MLST), respectively [7]. The surface antigens and their genes are under constant selective pressure and therefore, resulting in greater diversity which made them suitable only for short term or localised epidemiological tracking. In contrast, housekeeping enzymes and their genes are neutral (not under selection) and are therefore, ideal targets for long term and global epidemiology [4]. In MLST, isolates with unique allelic profiles are assigned to sequence types (STs) and related STs are grouped together into clonal complexes (CCs). Most IMD, especially those occurring in clusters or epidemics are caused by strains recognised as hypervirulent clones [4]. Other additional targets may include the iron-regulated outer membrane protein FetA [8], as well as other potential non-capsule protein-based vaccine components (see below).
The capsules of meningococci are well known virulence factor by offering the pathogen resistance to innate immunity including activation of the alternative pathway of complement. Therefore, vaccines based on capsules of serogroups A, C, W, and Y (MenA, C, W, and Y) have been successfully used to control these serogroups of meningococci [9]. However, the capsule of serogroup B meningococci (MenB) is non-immunogenic due to its cross-reactivity to host tissue and therefore, development of a MenB vaccine based on its capsule is not feasible [10]. Genome sequencing of N. meningitidis [11] has opened a new approach, called reverse vaccinology, to vaccine development [12]. Two new MenB vaccines have been licensed in some countries. Protein-based meningococcal vaccines developed against MenB may have the potential to protect also against many clones belonging to any serogroups, including MenC, MenY, and MenW [13] or even non-capsular clones. The 4CMenB vaccine or Bexsero® (Novartis) [14] is licensed for use in Europe, Canada, the United States, and Australia, while the recombinant bi-valent vaccine Trumenba® based on factor H binding protein (fHbp), also known as lipoprotein LP2086 [15], is licensed for use in the United States. The Bexsero® vaccine is made up of three major components, fHbp, subfamily B or variant 1, peptide 1; Neisserial Heparin Binding Antigen (NHBA) peptide 2; Neisseria Adhesin A (NadA) peptide 8, combined with the outer membrane vesicle (OMV) vaccine prepared from the New Zealand MenB strain NZ98/254 which has been used for control of the MenB epidemic in New Zealand caused by the MenB strain of B:4:P1.4 (Lineage 3 or ST-41/44 CC) [16]. Meningococcal fHbp is divided into 2 subfamilies, A and B which correspond to variant 2/3 and 1, respectively [17,18]. Proteins within a subfamily share a relatively high degree of similarity and therefore, crossprotection between peptides within a subfamily is feasible while peptides between subfamilies are more divergent and do not provide cross-protection [18]. Variants of NHBA have also been shown in animal studies to offer cross-protection [19]. NadA is made up of at least 4 variants, 1, 2/3, 4/5, and 6; with proteins within variants 1, and 2/3 shown to provide cross-protection [20]. The protection afforded by OMV vaccine is PorA or serosubtypespecific and the major immunological epitope resides on the PorA VR2 region. Thus OMV vaccine made from MenB strain NZ98/254 offers only protection against meningococcal strains expressing the PorA antigen P1.4 [16].
The epidemiology of invasive meningococcal disease in Québec, Canada, has been dominated by the ST-269 clone of MenB in the last decade associated with outbreaks or clusters of cases in certain regions of the province [21]. The aim of the current study is to document the trend since our last report as well as to expand on our previous studies to include an analysis of the genetic diversity of the newer MenB vaccine antigens, i.e. fHbp, NHBA, NadA, and PorA. This study will provide good baseline data before the newer MenB vaccines are introduced in the population.

Isolates of N. meningitidis
Individual IMD case isolates (n = 263), defined by their isolation from normally sterile body sites of individual patients were included in this study. These isolates constituted all the IMD isolates received at the National Microbiology Laboratory (NML) of the Public Health Agency of Canada between the period of 1 January 2009 and 31 December 2013.
Phylogenetic trees were constructed based on the predicted amino acid sequences of fHbp and NHBA peptides using software provided in the DNASTAR Lasergene 12.1 (Madison, Wisconsin, U.S.A.).

Serogroup distribution per year
During the period 2009-2013, 352 IMD cases were reported to Public Health authorities in Québec. Our study included 263 individual case isolates obtained from culture-confirmed IMD patients which constituted 75 % of all IMD reported in Québec. Overall, there were 229 (87 %) MenB, 16 (6 %) MenY, 10 (4 %) MenW, 7 (3 %) MenC, and 1 isolate of MenX (Table 1). This study did not include the culture-negative and PCR diagnosed IMD cases, and this might be the major reason for the smaller number of case isolates included in this study compared to the total number of cases reported in the province.

Clonal analysis of MenB
Among the 229 individual MenB case isolates, 220 or 96 % were grouped by MLST into 11 CCs, while 9 isolates did not belong to any known CC according to the current MLST typing scheme (http://pubmlst.org/ neisseria/; accessed April 30, 2015).    The other MenB case isolates (n = 27), whether belonging to a known CC or that did not belong to any known CC, were in general presenting diverse STs, serotype or serosubtype antigens. No one ST, nor serotype/ serosubtype was found in large numbers. including four with serotype 14,19; two with serotype 2c; and one each of serotype 19, and 15,19. Four different combinations of serosubtype antigens were found among the MenY isolates, including six isolates with P1.5; five with P1.5,2; and one each of P1.6 and P1.9. Three MenY were P1.-. Nine different combinations of serotype:serosubtype combinations were found among the 16 MenY isolates (data provided in Supplementary Table S1, see Additional file 1).

Characteristics of non-carbohydrate protein-based vaccine targets in MenB
Similar to the clonal data of MenB isolates, a limited number of fHbp and NHBA peptide types were predicted among isolates in the ST-269 CC (Additional file 2: Table S2). There were only three fHbp peptide types identified with 152 isolates predicted to have fHbp peptide 1.15; five for peptide 2.19 and one for peptide 1.249. In one isolate, no peptide was predicted based on its fHbp gene sequence identified as allele 755 which contained a frame-shift mutation. The five NHBA peptide types predicted to be produced by isolates of the ST-269 CC included peptide 21 (149 isolates), peptide 770 which differed from peptide 21 by only 1 amino acid (4 isolates), peptide 6 (3 isolates), peptide 768 (2 isolates), and peptide 122 (1 isolate).
The remaining MenB isolates were made up of many different STs belonging to 9 CCs as well as some that were not assigned to any known CC. Because the number of isolates in each CC or ST were very small, they would not be described individually but were grouped together for an overview of MenB strains. Of the 229 MenB IMD isolates in Québec, 25 fHbp peptide types were identified based on their fHbp gene sequences, including 12 peptide types with 175 isolates (76.4 %) belonging to variant 1 or subfamily B; and 13 peptide types with 52 isolates (22.7 %) belonging to subfamily A (10 peptides with 47 isolates belonged to variant 2; and 3 peptides with 5 isolates belonged to variant 3) ( Table 2). In addition, there was 1 isolate with a fHbp gene allele 755 that did not suggest a peptide would be produced and another isolate (ST-6977) whose fHbp gene sequence suggested a genetic recombination, similar to that described for strain M99 241177 (GenBank JF916571) [28]. Based on their nhba gene sequences, 21 NHBA peptide types were expected to be present in the 229 MenB isolates in Québec. The most common NHBA peptide was peptide 21 with 154 isolates and another 4 isolates predicted to synthesise peptide 770, which was related to NHBA peptide 21 but with only 1 amino acid difference. Other common NHBA peptide types included peptide 112 (25 isolates), peptide 2 (6 isolates), and peptides 6, 10, 47, and 122, each with four isolates. The remaining NHBA peptide types were made up of either 1 or 2 isolates each ( Table 3).
The nadA gene could not be detected in MenB strains belonging to ST-269 and ST-41/44 CCs. In the other MenB strains, the nadA gene could only be found in

Characteristics of non-carbohydrate protein-based vaccine targets in non-serogroup B meningococci
In two of the seven MenC isolates, their fHbp gene allele 669 contained a frame-shift mutation that led to a premature stop codon and therefore, no fHbp peptide would be expected in them. The other five isolates each have a different fHbp gene allele predicting the following fHbp peptides (1.12, 2.16, 2.102, 1.129, and 1.218); three of which belonged to variant 1 and two to variant 2. The nadA genes (allele 29) in four MenC isolates were disrupted by presence of an IS1301 element, and therefore, no peptide were expected in them. In one MenC isolate, its nadA gene predicted a NadA-2/3 peptide 2; and in two MenC isolates, their nadA gene predicted a NadA-2/3 peptide 3. All seven MenC isolates were found to have nhba genes that predicted a NHBA peptide 20.
Of the 16 MenY isolates, seven fHbp peptide types were predicted based on their fHbp gene sequences. Other than the one isolate that was predicted to produce a variant 1 fHbp peptide 1.278, the remaining 15 MenY isolates were predicted to encode variant 2 fHbp peptides 2.16, 2.21, 2.23, 2.24, 2.25, and 2.629; with peptide 2.25 being the most common peptide found in seven isolates. Seven NHBA peptide types (6,7,8,9,20, 680, and 769) were identified based on their nhba gene sequences. Peptide 680 was related to peptide 6 but with 1 amino acid difference. Only the single ST-1466 (ST-174 CC) MenY was found to have a nadA gene allele 80, which predicted a NadA-2/3 peptide 8. The nadA gene was not detected in all the other MenY isolates.
All ten MenW isolates were found to have fHbp genes that encoded for variant 2 peptides, with nine of them predicted to encode for peptide 2.16, and only one isolate was predicted to produce a fHbp peptide 2.25. All ten MenW were found to have nhba gene allele encoding for NHBA peptide 20. None of the MenW was found to harbor the nadA gene and therefore, no NadA protein was expected to be present in these MenW isolates.
The single MenX was found to have fHbp and nhba genes to encode for fHbp peptide 2.19 (variant 2), and NHBA peptide 767. No nadA gene was detected in the MenX isolate.
Predicted coverage by the 4CMenB vaccine based on vaccine target gene sequence analysis There were altogether 30 fHbp peptide types identified among the invasive N. meningitidis strains found in Québec. Their frequencies of occurrence in different serogroups are described in Table 2. In addition, there were 3 isolates with fHbp gene mutations that would prevent synthesis of a mature fHbp peptide. In another isolate, there was a recombination resulting in an interruption of its fHbp gene by a foreign unrelated gene (GenBank JF916571) [28]. There were 27 different NHBA peptide types and their distribution according to serogroup is described in Table 3.
For MenB isolates, none was found to have more than two vaccine antigens that matched exactly with the 4CMenB components and only two isolates were found to have PorA P1.4 and NHBA peptide 2 that matched exactly with the respective components found in the vaccine. Another six MenB isolates were found to have one antigen, either fHbp peptide 1.1 (two isolates) or NHBA peptide 2 (four isolates), that matched exactly with the corresponding component found in the 4CMenB vaccine (Table 4). Assuming all isolates with fHbp variant 1 and NadA variants 1 and 2/3 peptides would be covered by the vaccine, two MenB isolates might be regarded to be covered by three components of the vaccine by virtue of their PorA P1.4, fHbp peptide 1.4, and NHBA peptide 2 antigens; and five MenB isolates might be regarded to be covered by two components of the vaccine by virtue of their variant 1 fHbp peptides and NHBA peptide 2 (four isolates), and variant 1 fHbp peptide and NadA peptide 85 (one isolate). Another 166 MenB isolates were considered to be covered by one component of the vaccine: variant 1 fHbp (163 isolates) or NadA variants 1, 2/3 peptides (3 isolates).
For the non-serogroup B IMD strains in Québec, none was found to have fHbp and nhba genes encoding for fHbp peptide 1.1 and NHBA peptide 2. However 1 MenY isolate (ST-1466) was found to possess a nadA gene for NadA peptide 8 found in the 4CMenB vaccine. Another 4 isolates (3 MenC and 1 MenY) were found to have fHbp genes that predicted variant 1 fHbp peptides. Also 3 MenC were found to have nadA genes encoding for NadA-2/3 peptides, likely to have enough similarity with the NadA peptide component in the 4CMenB to be regarded as covered by the vaccine.
Regardless of serogroup, a total of 9 isolates could be regarded to be covered by the 4CMenB vaccine by possession of at least 1 antigen showing an exact match to the vaccine components; and an additional 172 isolates might be covered by virtue of possession of at least 1 antigen showing a partial match with the vaccine components, e.g. fHbp variant 1 peptides.
The diversity of the fHbp and NHBA peptides found in the invasive N. meningitidis strains in Québec are presented in Figs. 1 and 2, respectively. With respect to fHbp peptides, even among the variant 1 peptides, considerable amino acid diversity was found among IMD isolates in Québec ( Fig. 1; bottom half of the figure). In addition, it was evident that peptides of the variants 2 and 3 were more distantly related to peptides that belong to variant 1. Similarly, NHBA peptides found in the Québec IMD isolates were also grouped into clusters, with some clusters of peptides showing various degree of diversity with the vaccine peptide 2 (Fig. 2).

Discussion
Factor H binding protein (fHbp) is one of Neisseria meningitidis virulence factors and is found in two meningococcal vaccines, Bexsero® and Trumenbal®. Characterization of fHbp from various N. meningitidis strains have identified 2 subfamilies (A and B) or 3 variants (variants 1, 2, 3) of fHbp. Subfamily A fHbp is equivalent to variants 2 and 3 while subfamily B fHbp is equivalent to variant 1. Immunological cross-reactivity and cross-protection between different fHbp within a subfamily have been observed but this is not extended to fHbp from different subfamilies [18]. Bexsero® also known as 4CMenB, is made up of multi-components with fHbp being only one of them. The fHbp present in this vaccine is restricted to peptide 1.1 from subfamily B, which presumably will cover strains expressing subfamily B fHbp. Broader protective nature of Bexsero® is extended by its other components, including NHBA, NadA and PorA P1.4. On the other hand, Trumenbal® is composed of only fHbp but it is a bivalent fHbp vaccine with 2 peptides present representing both subfamily A and subfamily B that allows it to cover strains expressing either subfamily of fHbp. This study provided a comprehensive genetic analysis of this virulence factor found in invasive N. meningitidis strains in the province of Quebec, Canada over the five years period of 2009-2013, just before the licensure of either of these vaccines.
This study was based on our previous reports of a changing epidemiology of IMD in the province of Quebec [21,29,30] but the focus was on characteristics of invasive N. meningitidis strains and the genetic diversity of vaccine antigens found in the newer MenB vaccines. Results of this study documented the continued dominance of MenB as a cause of IMD in Québec, with 92 % of the IMD case isolates in 2013 belonging to this serogroup (Table 1), and in particular the ST-269 CC. Compared to the study that examined IMD case isolates in Québec from 2003 to 2010 which consisted of 180 MenB isolates belonging to the ST-269 CC with 91.7 % due to a strain typed as ST-269 [30], the present study showed that the % of ST-269 isolates in the ST-269 CC has not changed and remained as 91.2 % (145 of 159 isolates of MenB ST-269 CC). In the original report of the emergence of the ST-269 clone, there were only 3 STs identified among the MenB isolates in this CC in Québec [29]. Then the number of STs within this CC expanded to involve 11 STs [30], with 8 STs newly identified when compared to the previous finding [29]. The  1.4, 1.12, 1.13, 1.14, 1.15, 1.100, 1.110, 1.129, 1.144, 1.218, 1.249, 1.278, 1.410 present study also identified ten STs within the Quebec ST-269 CC but four of the STs found previously [30] were not found in the present study, and three new STs not previously found were included in here. Also of interest is ST-275 and its SLV ST-1163 found in the period of 2003-2010 did not appear to have expanded to cause IMD in Québec. This was in contrast to the endemic MenB disease reported in England/Wales, where the ST-275 clone has expanded while ST-269 clone has decreased in frequency of detection among IMD case isolates [27]. The difference between the molecular epidemiology of MenB disease in Québec and England/ Wales might be in Québec MenB ST-269 has been associated with clusters of cases whereas in England/ Wales, the MenB disease has been characterized as endemic disease without any single clone dominating. Our data may suggest that ST-269 strain in Québec might have a unique advantage in establishing itself in the population to cause significant disease over an extended period of time [30]. Indeed we speculated that MenB strain of ST-269 has the potential to cause outbreaks [31].
Among the MenB ST-41/44 CC isolates, ST-571 continued to be a common ST found, responsible for 37 % of the MenB in this CC, similar to the finding reported previously [30]. Isolates of ST-571 also appeared to have a uniform antigenic formula expressing serotype antigen 19 (in all 16 isolates), and serosubtype antigen P1.9 (in 15 isolates) together with the PorA genotype of P1. 18-7, 9, 35-1 (in 15 isolates). Other MenB isolates within the ST-41/44 CC or other CCs or those not belonging to any known CC can essentially be described as diverse, not restricted to any particular ST or antigenic form.
Analysis of the 4CMenB vaccine antigen genes among the various meningococcal serogroups and CCs have provided some interesting observations (Table 4). Although MenB that belonged to the ST-269 CC might not have any antigens that matched exactly with the 4CMenB vaccine components, a significant percentage (96 % or 153 out of 159 isolates) were predicted to be covered based on the majority of them possessed fHbp genes that encoded for variant 1 fHbp peptide 1.15. The coverage of this clone by the vaccine was further supported by an in vitro assay Meningococcal Antigen Typing System (MATS) that measures expression, and cross-reactivity of the vaccine antigens on clinical isolates. Using MATS, the majority of strains within the ST-269 CC were covered by fHbp and NHBA; and the coverage for strains typed as ST-269 was 95 % [32]. Also all five MenB isolates of the ST-32 CC were also predicted to be covered by the vaccine, either because they possessed genes that encoded for a matching 4CMenB  Table 4).
The number of MenB isolates in this study predicted to synthesize variant 1 fHbp was greatly influenced by the presence of a large number of the ST-269 CC isolates which predominantly predicted to have fHbp peptide 15 (variant 1). If this clone was removed from the calculation, the number and percentage of MenB isolates predicted to have variant 1, 2, 3 fHbp would be 23, 47, and 5 isolates or 30.7, 62.7, and 6.7 %, respectively. The percentage of MenB in Québec without the ST-269 CC would be more comparable to the percentage of MenB in Western Canada predicted to have variant 1 fHbp (37.7 %) [33]. In Western Canada, the percentage of MenB predicted to encode for variants 2 and 3 fHbp were 34.8 % and 27.5 %, respectively The higher percentage of variant 3 fHbp found among MenB in Western Canada (27.5 %) was due to the presence of the ST-213 cc, which were predominantly (92 %) predicted to have fHbp peptide 45 (variant 3) [33].
The number of MenC, Y, W, and X were relatively small and therefore, prediction of how well the vaccine may or may not work against these serogroups of meningococi in Québec may change when more isolates are tested. The 34 MenC, Y, W, and X isolates were found to possess diverse alleles of fHbp and nhba genes predicted to encode for 13 different fHbp and eight NHBA peptide types. However, only 4 (3 MenC and 1 MenY) isolates were predicted to have variant 1 fHbp while 28 isolates were predicted to have variant 2 fHbp; and 2 MenC were predicted not to produce fHbp due to a frame-shift mutation in their fHbp genes [28]. Among MenC, no two isolates were found to have identical fHbp gene alleles; a phenomenon linked to the ET-15 variant of MenC in the ST-11 CC [34]. Also 4 of the 7 MenC isolates were found to have their nadA genes disrupted by presence of IS1301 in the coding region, thereby leading to no NadA protein predicted in these isolates [35]. NadA-2/3 peptides were predicted in 3 of the 34 non-B meningococci, including 2 MenC (NadA peptide 3) and one MenY (NadA peptide 8).
Despite conjugated MenC vaccination programs implemented since 2001-2002 to cover all infants and adolescents in the province, it appears invasive MenC disease still occurred at a low level in the adult population. Six of the seven MenC cases were in the age range of 30 to 83 years (with a median age of 59.5 years); and in one case the age information is not available. Also the responsible MenC strain appears to be identical to the hyper-virulent ET-15 clone prominent in prevalence before the MenC conjugate vaccine was implemented. The continuing circulation of this MenC ET-15 clone with the potential to cause outbreaks with severe outcome [36], justifies the existence of an immunisation program consisting of two MenC conjugate vaccine doses offered at 12 months and 12 years of age, respectively [37], in order to maintain the immunity in the population and to keep the level of invasive MenC disease to a minimum in the province. Genetic analysis of new MenB vaccine targets suggested that the MenC ET-15 clone may be covered by the fact that they all have nhba genes for NHBA peptide 20 (which may cross-react with NHBA peptide 2), or either presence of fHbp genes for variant 1 fHbp peptides in three isolates or nadA genes that encoded for NadA peptide 2 or 3 in another three isolates.
With the implementation of the Bexsero® vaccine in regions of Quebec with the highest incidence of IMD due to the ST-269 clone (authors' unpublished data), attention are now being focused on gathering data that would provide information on the efficacy of the vaccine, its safety profile in the field, its acceptability by the population, its effect on potential carriage and herd immunity, as well as its potential effect on the circulating population of N. menngitidis strains. This post-implementation surveillance data would help to improve our knowledge of this new vaccine and may assist in refining the dosage and/or schedule of the vaccine, as well as its true strain coverage.

Conclusion
This study provided an overall molecular analysis of the invasive N. meningitidis in Québec in recent years with an intent to examine the genetic diversity and distribution of the newer meningococcal vaccine targets. For MenB, eight strains were predicted to have one (six isolates) or two (two isolates) antigens that matched exactly with the respective component in the 4CMenB vaccine. Six isolates (including two that also expressed the PorA P1.4 antigen) were predicted to encode NHBA peptide 2, and two isolates predicted to encode variant 1 fHbp peptide 1.1. Another 165 MenB isolates might be regarded to be covered by the vaccine by virtue of either possessing genes for variant 1 fHbp or NadA variants 1, 2/3 peptides. For the 34 non-MenB isolates, only one MenY was found to possess a nadA gene for the NadA peptide 8 found in the 4CMenB vaccine. In addition seven other isolates predicted to be potentially covered by the vaccine by virtue of either possession of a fHbp gene for variant 1 fHbp peptide (four isolates) or a nadA gene for variant 2/3 NadA peptides (three isolates).
The major limitation is we did not perform assay to examine expression of the newer vaccine targets in our clinical isolates and their cross-reactivities with those found in the vaccine. Also the number of isolates as a whole were small which might affect the confidence level of the finding. Furthermore, we did not have clinical data on hand to compare the virulence of the different genotypes identified, in particular those in the ST-269 CC. Therefore, further prospective studies employing these tools and approaches may allow for a more in-depth analysis of IMD in Québec, Canada.

Additional files
Additional file 1: Table S1. A list of invasive meningococcal disease isolates in Quebec, Canada from 2009 to 2013. Authors' contributions RSWT and DKSL designed the study; BL, SD, JZ, and DKSL carried out the laboratory study and the interpretation of the data; BL, RG, PDW and RSWT participated in the analysis of the data and preparation of the manuscript, and all authors read, and approved the content of the manuscript.