Skip to main content

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

Abstract

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.

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, cross-protection 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 serosubtype-specific 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.

Methods

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.

Identification and typing of N meningitidis

Isolates of N. meningitidis were identified at the local hospitals and their identity was confirmed and they were serogrouped at the Laboratoire de santé publique du Québec using rabbit anti-meningococcal antisera produced by the NML. Identification of serotype and serosubtype was done by a monoclonal antibody kit containing antibodies to serotypes 1, 2a, 2b, 4, 14, and 15; and serosubtypes P1.1, P1.2, P1.4, P1.5, P1.6, P1.7, P1.9, P1.10, P1.12, P1.13, P1.14, P1.15, and P1.16 (Rijksinstitut voor Volksgezondheid en Milieu, National Institute of Public Health, Bethoven, The Netherlands) using indirect whole cell ELISA [22]. Mouse monoclonal antibodies to serotype 2c, 17, 19, and serosubtype P1.19 were kind gifts from Dr. Wendell Zollinger of the Walter Reed Army Institute of Research (Maryland, U.S.A.). Clonal analysis was done by MLST [23] and ST/CC were assigned using tools available in the Neisseria MLST website (http://pubmlst.org/neisseria/). PorA genotype was determined by PCR amplification of the variable region gene fragments followed by standard DNA sequencing technique using published methods [24, 25]. The nomenclature given in the N. meningitidis PorA variable region database (http://pubmlst.org/neisseria/PorA/) and by de Filippis et al. [26] was followed. The fHbp, nhba, and nadA gene sequences were determined by PCR amplification and standard sequencing reactions following the protocols described by Lucidarme et al. [27]; and their peptide types were determined by the online tools available from the Neisseria.org website (http://pubmlst.org/neisseria). The nomenclature scheme for fHbp peptide according to Novartis was used: e.g. variant 1, peptide 1 was described as 1.1, etc.

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.).

Results

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.

Table 1 Distribution of Neisseria meningitidis serogroups among culture-positive invasive meningococcal disease cases in Québec, Canada, 2009-2013

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 most commonly occurring CCs were ST-269 CC (n = 159 or 69.4 %), followed by ST-41/44 CC (n = 43 or 18.8 %), ST-32 CC (n = 5 or 2.2 %), ST-35 CC (n = 4 or 1.8 %), and ST-461 CC and ST-1157 CC (each with 2 isolates or 0.9 %). There were also 5 other CCs (ST-37, ST-60, ST-103, ST-213, and ST-254 CCs) each with only a single isolate (see Additional file 1: Table S1).

Despite being the most frequently occurring CC with 159 isolates, 91 % (n = 145) of the isolates in the ST-269 CC belonged to a single ST (ST-269). There were 9 other STs with 7 of them (ST-565, ST-1986, ST-5494, ST-8767, ST-8772, ST-8924, and ST-8880) being single locus variants (SLVs) of ST-269 while ST-13 and ST-2738 were double, and triple locus variants, respectively. Most (n = 150 or 94 %) expressed the serotype antigen 17, with 8 isolates being non-serotypeable (NT), and 1 isolate of serotype 14,19. One hundred and twenty two isolates (77 %) expressed the serosubtype antigen P1.19 and 25 isolates (16 %) expressed the serosubtype antigen P1.9. There were 6 isolates that were non-serosubtypeable (P1.-), and 1 isolate each of serosubtypes P1.6; P1.7; P1.12; P1.16; P1.5,2; and P1.15,29.

In contrast to the ST-269 CC, the 43 MenB case isolates that belonged to the ST-41/44 CC were grouped into 23 different STs (Additional file 1: Table S1). Although 1 ST (ST-571) was common (responsible for 16 isolates or 37 % of the ST-41/44 CC), 19 other STs were represented by a single isolate only, and another 3 STs with multiple isolates (ST-944, n = 3; ST-282, n = 3; and ST-41, n = 2). The diversity of the ST-41/44 CC isolates was also reflected in their serotype and serosubtype antigens. Altogether, there were 7 different combinations of serotype antigens, including seventeen isolates of serotype 19; fourteen of serotype 4; four of serotype 15,19; three of serotype 17; two of serotype 1,19; and one each of serotypes 15, and 14,19. There was only one isolate that was non-serotypeable (NT). There were 8 different combinations of serosubtype antigens, including twenty of P1.9; five of P1.6; three of P1.13; two each of P1.4; P1.5; P1.12; and one each of P1.14; and P1.5, 2. Seven isolates were non-serosubtypeable (P1.-). The antigenic formula of the ST-571 isolates appeared to be generally uniform with thirteen isolates typed as B:19:P1.9; and one isolate each of B:1,19:P1.9; B:15,19:P1.9; and B:15,19:P1.-.

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.

Clonal analysis of non-group B meningococci

All 7 MenC isolates belonged to ST-11 (ST-11 CC) and expressed the serotype antigen 2a. Three different serosubtype antigen combinations were detected including four isolates of P1.5; two isolates of P1.7,1; and one isolate of P1.2.

The four different STs found among the 10 MenW isolates all belonged to the ST-22 CC. Besides the single isolate of ST-22, there were 7 isolates of ST-184, and 1 isolate each of ST-1617 and ST-2625; with the latter 3 STs being SLVs of ST-22. All ten MenW were NT; five were positive for P1.6; one was positive for P1.14; and four were P1.-.

Twelve of the 16 MenY were grouped into 4 CCs, with ST-23 CC being the most common, made up of 6 isolates, followed by ST-167 CC with 4 isolates, and 1 isolate each of ST-174 CC and ST-22 CC. Four MenY isolates were not classified into any known CC. Half of the MenY isolates were NT; and the remaining 8 isolates presented 4 different combinations of serotype antigens 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).

In contrast to isolates in the ST-269 CC; isolates of ST-41/44 CC were also more diverse in their fHbp and nhba gene sequences (Additional file 2: Table S2). Ten fHbp gene alleles together with five nhba gene alleles were found among isolates in the ST-41/44 CC. They predicted fHbp peptides 1.4 (n = 4), 1.14 (n = 1), 2.19 (n = 26), 2.23 (n = 2), 2.24 (n = 2), 3.30 (n = 1), 1.100 (n = 1), 1.410 (n = 4), 2.632 (n = 1), and 1.687 (n = 1); and NHBA peptides 2 (n = 6),10 (n = 4), 29 (n = 4), 47 (n = 3), 112 (n = 25), and 286 (n = 1).

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).

Table 2 Factor H binding protein (fHbp) peptide types in invasive Neisseria meningitidis from Québec, Canada, 2009-2013
Table 3 Neisserial Heparin Binding Antigen (NHBA) peptide types in invasive Neisseria meningitidis from Québec, Canada, 2009-2013

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 isolates belonging to the ST-32 CC (n = 5); ST-1157 CC (n = 2); ST-213 CC (n = 1), and ST-336 (n = 1). However, the nadA genes in isolates of ST-1157 CC (allele 20) and ST-213 CC (allele 34) did not predict to produce any peptide due to frame-shift mutations. In the single isolate of ST-336, the nadA gene was found to predict NadA-2/3 peptide 3. Four of the 5 ST-32 CC isolates had the nadA gene which predicted NadA-1 peptide 1; while the nadA gene in the fifth isolate was predicted to encode NadA-1 peptide 85.

Among the 229 MenB case isolates, only two expressed the serosubtype antigen P1.4; with the PorA genotype of P1.7-2,4,37. Both isolates belonged to the ST-41/44 CC (one being of ST-41 and the other was ST-154).

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).

Table 4 Invasive Neisseria meningitidis from Québec, Canada in 2009-2013: predicted vaccine coverage based on at least one antigen showing an exact match to the 4CMenB vaccine components or with Novartis variant 1 factor H binding protein (fHbp)/Neisseria adhesin A (NadA) peptides belonging to variants 1 to 3

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).

Fig. 1
figure 1

Phylogenetic tree to show the relationship of factor H binding protein peptide types identified in invasive Neisseria meningitidis from Québec, Canada, 2009-2013. *Peptide Id from the Neisseria.org website (http://pubmlst.org/neisseria/fHbp/). Factor H binding protein peptide Id 1.1, 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, and 1.687 belong to Novartis variant 1 or subfamily B. Factor H binding protein peptide Id 2.16, 2.19, 2.21, 2.22, 2.23, 2.24, 2.25, 2.102, 2.106, 2.527, 2.629, and 2.632 belong to Novartis variant 2 or subfamily A. Factor H binding protein peptide Id 3.30, 3.47, and 3.555 belong to Novartis variant 3 or subfamily A. Dotted line indicates a negative branch length

Fig. 2
figure 2

Phylogenetic tree to show the relationship of Neisseria Heparin Binding Antigen peptide types in invasive Neisseria meningitidis from Québec, Canada, 2009-2013. *Peptide Id from the Neisseria.org website (http://pubmlst.org/neisseria/NHBA). Variants were labelled with the peptide Id. Dotted line indicates a negative branch length

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 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 vaccine component (variant 1 fHbp peptide 1.1; two isolates) or with genes that encoded for NadA variants 2/3 peptides (3 isolates). This finding is in line with other observations that many MenB of ST-32 CC are likely to be covered by the 4CMenB vaccine [19]. In contrast a small number of MenB isolates of the ST-41/44 CC (14 % or 6 out of 43 isolates) were predicted to be covered, possibly because of the heterogeneity of strains found in this CC. Of the 27 MenB that did not belong to ST-269 and ST-41/44 CCs, only 12 or 44 % were predicted to be covered by the 4CMenB vaccine. Therefore, aside from ST-269, ST-32 and ST-41/44 CCs, MenB of other clonal types in Quebec might be only poorly covered by the vaccine according to the strains analysed in this study. Besides MenB, 6 of the 7 ST-11 MenC were predicted to be covered by the 4CMenB vaccine by virtue of either their nadA genes encoding for variant 2/3 NadA peptides (3 isolates) or fHbp genes encoding for variant 1 fHbp (3 isolates). In contrast, none of the 10 MenW and only 2 of the 16 MenY were predicted to be covered by the vaccine (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.

References

  1. Rosenstein NE, Perkins BA, Stephens DS, Popovic T, Huges JM. Meningococcal disease. N Engl J Med. 2001;344:1378–88.

    Article  CAS  PubMed  Google Scholar 

  2. Al-Tawfiq JA, Clark TA, Memish ZA. Meningococcal disease: the organism, clinical presentation, and worldwide epidemiology. J Trop Med. 2010;17(Suppl):3–8.

    Google Scholar 

  3. Xie Q, Pollard AJ, Mueller JE, Norheim G. Emergence of serogroup X meningococcal disease in Africa: need for a vaccine. Vaccine. 2013;31:2852–61.

    Article  PubMed  Google Scholar 

  4. Caugant DA. Population genetics and molecular epidemiology of Neisseria meningitidis. APMIS. 1998;106:505–25.

    Article  CAS  PubMed  Google Scholar 

  5. Yazdankhah SP, Kriz P, Tzanakaki G, Kremastinou J, Kalmusova J, Musilek M, et al. Distribution of serogroups and genotypes among disease-associated and carried isolates of Neisseria meningitidis from the Czech Republic, Greece, and Norway. J Clin Microbiol. 2014;42:5146–53.

    Article  Google Scholar 

  6. Harrison LH, Trotter CL, Ramsay ME. Global epidemiology of meningococcal disease. Vaccine. 2009;27:B51–63.

    Article  PubMed  Google Scholar 

  7. Harrison OB, Brueggemann AB, Caugant DA, van der Ende A, Frosch M, Tray S, et al. Molecular typing methods for outbreak detection and surveillance of invasive disease caused by Neisseria meningitidis, Haemophilus influenzae, and Streptococcus pneumoniae, a review. Microbiolology. 2011;157:2181–95.

    CAS  Google Scholar 

  8. Thompson EA, Feavers IM, Maiden MC. Antigenic diversity of meningococcal enterobactin receptor FetA, a vaccine component. Microbiology. 2003;149:1849–58.

    Article  CAS  PubMed  Google Scholar 

  9. Warshawsky B. Update on the use of quadrivalent meningococcal vaccines. An Advisory Committee Statement (ACS) National Advisory Committee on Immunization. Can Comm Dis Rep. 2013;39(ACS-1):1–40.

    Google Scholar 

  10. Finne J, Leinonen M, Makela PH. Antigenic similarities between brain components and bacteria causing meningitis. Implications for vaccine development and pathogenesis. Lancet. 1983;322:355–7.

    Article  Google Scholar 

  11. Tettelin H, Saunders NJ, Heidelberg J, Jeffries AC, Nelson KE, Eisen JA, et al. Complete genome sequence of Neisseria meningitidis serogroup B strain MC58. Science. 2000;287:1809–15.

    Article  CAS  PubMed  Google Scholar 

  12. Sette A, Rappuoli R. Reverse vaccinology: developing vaccines in the era of genomics. Immunity. 2010;33:530–41.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  13. Perez-Trallero E, Esnal O, Marimon JM. Progressive decrease in the potential useflness of meningococcal serogroup B vaccine (4CMenB, Bexsero®) in Gipuzkoa, Northers Spain. PLoS One. 2014;9, e116024.

    Article  PubMed Central  PubMed  Google Scholar 

  14. Serruto D, Bottomley MJ, Ram S, Giuliani MM, Rappuloi R. The new multicomponent vaccine against meningococcal serogroup B, 4CMenB: immunological, functional and structural characterization of the antigens. Vaccine. 2012;30 Suppl 2:B87–97.

    Article  CAS  PubMed  Google Scholar 

  15. Richmond PC, Marshall HS, Nissan MD, Jiang Q, Jansen KJ, Garces-Sanchez M, et al. Safety, immunogenicity, and tolerability of meningococcal serogroup B bivalent recombinant lipoprotein 2086 vaccine in healthy adolescents: a randomized, single-blind, placebo-controlled, phase 2 trial. Lancet Infect Dis. 2012;12:597–607.

    Article  CAS  PubMed  Google Scholar 

  16. Oster P, Lennon D, O’Hallahan MK, Reid S, Martin D. MeNZB™: a safe and highly immunogenic tailor-made vaccine against the New Zealand Neisseria meningitidis serogroup B disease epidemic strain. Vaccine. 2005;23:2191–6.

    Article  CAS  PubMed  Google Scholar 

  17. Fletcher LD, Bernfield L, Barniak V, Farley JE, Howell A, Knauf M, et al. Vaccine potential of the Neisseria meningitidis 2086 lipoprotein. Infect Immun. 2004;72:2088–100.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  18. Masignani V, Comanducci M, Giuliani MM, Bambini S, Adu-Bobie J, Arico B, et al. Vaccination against Neisseria meningitidis using three variants of the lipoprotein GNA1870. J Exp Med. 2003;197:789–99.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  19. Giuliana MM, Adu-Bobie J, Commanducci M, Arico B, Savino S, Santini L, et al. A universal vaccine for serogroup B meningococcus. Proc Natl Acad Sci U S A. 2006;103:10834.

    Article  Google Scholar 

  20. Bambini S, De Chiara M, Muzzi A, Mora M, Lucidarme J, Brehony C, et al. Neisseria Adhesin A variation and revised nomenclature scheme. Clin Vaccine Immunol. 2014;21:966–71.

    Article  PubMed Central  PubMed  Google Scholar 

  21. Gilca R, Deceuninck G, Lefebvre B, Tsang R, Amin R, Gilca V, et al. The changing epidemiology of meningococcal disease in Quebec, Canada, 1991-2011: potential implications of emergence of new strains. PLoS One. 2012;7, e50659.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  22. Abdillahi H, Poolman JT. Whole-cell ELISA for typing Neisseria meningitidis with monoclonal antibodies. FEMS Microbiol Lett. 1987;48:367–71.

    Article  Google Scholar 

  23. Maiden MCJ, Bygraves JA, Feil E, Morelli E, Russell JE, Urwin R, et al. Multilocus sequence typing: a portable approach to the identification of clones within populations of pathogenic microorganisms. Proc Natl Acad Sci U S A. 1998;95:3140–5.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  24. Sacchi CT, Lemos AP, Brandt ME, Whitney AM, Melles CEA, Solari CA, et al. Proposed standardization of Neisseria meningitidis PorA variable-region typing nomemclature. Clin Diagn Lab Immunol. 1998;5:845–55.

    CAS  PubMed Central  PubMed  Google Scholar 

  25. Clark SC, Diggle MA, Molling P, Unemo M, Olcen P. Analysis of PorA variable region 3 in menincococci: implications for vaccine policy. Vaccine. 2003;21:2468–73.

    Article  Google Scholar 

  26. De Filippis I, Gopalan V, Huyen Y. PorA VR3 typing database: a web-based resource for the determination of PorA VR3 alleles of Neisseria meningitidis. Infect Gene Evol. 2011;11:248–9.

    Article  CAS  Google Scholar 

  27. Lucidarme J, Comanducci M, Findlow J, Gray SJ, Kaczmarski EB, Guiver M, et al. Characterization of fHbp, nhba (gna2132), nadA, porA, sequence type (ST), and genomic presence of IS1301 in group B meningococcal ST-269 clonal complex isolates from England and Wales. J Clin Microbiol. 2009;47:3577–85.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  28. Lucidarme J, Tan L, Exley RM, Findlow J, Borrow R, Tang CM. Characterization of Neisseria meningitidis isolates that do not express the virulence factor and vaccine antigen factor H binding protein. Clin Vaccine Immunol. 2011;18:1002–14.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  29. Law DKS, Lorange M, Ringuette L, Dion R, Giguere M, Henderson AM, et al. Invasive meningococcal disease in Quebec, Canada, due to an emerging clone of ST-269 serogroup B meningococci with serotype antigen 17 and serosubtype antigen P1.19 (B:17:P1.19). J Clin Microbiol. 2006;44:2743–9.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  30. Zhou J, Lefebvre B, Deng S, Gilca R, Deceuninck G, Law DKS, et al. Invasive serogroup B Neisseria meningitidis in Quebec, Canada, 2003 to 2010: persistence of the ST-269 clone since it first emerged in 2003. J Clin Microbiol. 2012;50:1545–51.

    Article  PubMed Central  PubMed  Google Scholar 

  31. Tsang RSW, Jamieson FB, Lefebvre B, Gilca R. Emerging sequence type-269 serogroup B Neisseria meningitidis has the potential to cause outbreaks. Can J Infect Dis Med Microbiol. 2014;25:8.

    PubMed Central  PubMed  Google Scholar 

  32. Bettinger JA, Scheifele DW, Halperin SA, Vaudry W, Findlow J, Borrow R, et al. Diversity of Canadian meningococcal serogroup B isolates and estimated coverage by an investigational meningococcal serogroup B vaccine (4CMenB). Vaccine. 2014;32:124–30.

    Article  Google Scholar 

  33. Law DKS, Zhou J, Deng S, Hoang L, Tyrrell G, Horsman G, et al. Determination of serotyping antigens, clonal analysis and genetic characterization of the 4CMenB vaccine antigen genes in invasive Neisseria meningitidis from Western Canada, 2009-2013. J Med Microbiol. 2014;63:1490–9.

    Article  PubMed  Google Scholar 

  34. Tsang RSW, Hoang L, Tyrrell G, Horsman G, Wylie J, Jamieson F, et al. Genetic and antigenic characterization of Canadian Neisseria meningitidis serogroup C (MenC) case isolates in the post-MenC conjugate vaccine era, 2009-2013. J Med Microbiol. 2015;64:174–9.

    Article  PubMed  Google Scholar 

  35. Elias J, Vogel U. IS1301 fingerprint analysis of Neisseria meningitidis strains belonging to the ET-15 clone. J Clin Microbiol. 2007;45:159–67.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  36. Erickson L, De Wals P. Complications and sequelae of meningococcal disease in Québec, Canada, 1990-1994. Clin Infect Dis. 1998;26:1159–64.

    Article  CAS  PubMed  Google Scholar 

  37. COMITÉ SUR L’IMMUNISATION DU QUÉBEC. Avis du Comité sur l’immunisation du Québec concernant l’utilité d’une dose de rappel de vaccin contre le méningocoque chez les adolescents au Québec. INSTITUT NATIONAL DE SANTÉ PUBLIQUE DU QUÉBEC, Québec: 2010. Available at: http://www.inspq.qc.ca/pdf/publications/1178_AvisCIQDoseRappelMeningocoqueAdos.pdf.

Download references

Acknowledgements

We thank the staff of the National Microbiology Laboratory’s DNA Core Facility for providing primers and help with DNA sequencing. We also would like to thank staff at the local and regional hospitals in Québec for providing the isolates for this study. RSWT thanks Dr. Wendell Zollinger, Walter Reed Army Institute of Research for his gifts of the serotyping and serosubtyping monoclonal antibodies.

This study made use of the Neisseria MLST website (http://pubmlst.org/neisseria/) developed by Keith Jolley and Man-Suen Chan and sited at the University of Oxford, UK. The development of this site has been funded by the Wellcome Trust and European Union.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Brigitte Lefebvre or Raymond S.W. Tsang.

Additional information

Competing interests

Dennis Law, Brigitte Lefebvre, Saul Deng, Jianwei Zhou, and Raymond Tsang do not have any competing interest to declare. Rodica Gilca received research grants from Pfizer and Sanofi for unrelated studies. Philippe De Wals has received research grants and reimbursement of travel expenses from vaccine manufactures including Pfizer, Novartis and GlaxoSmithKline.

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.

Additional files

Additional file 1: Table S1.

A list of invasive meningococcal disease isolates in Quebec, Canada from 2009 to 2013.

Additional file 2: Table S2.

4CMenB protein antigensΔ predicted in serogroup B Neisseria meningitidis belonging to ST-269 and ST-41/44 claonal complexes (CCs) and isolated in Québec, Canada, 2009-2013. These additional files as well as Figs. 1 and 2 showing phylogenetic trees of factor H binding protein and Neisserial Heparin Binding Antigen peptide types were also deposited in Dryad (http://datadryad.org/).

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Law, D.K., Lefebvre, B., Gilca, R. et al. 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. BMC Microbiol 15, 143 (2015). https://doi.org/10.1186/s12866-015-0469-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12866-015-0469-6

Keywords