Skip to main content

The impact of aggregating serogroups in dynamic models of Neisseria meningitidis transmission



Neisseria meningitidis (Nm) is a pathogen of multiple serogroups that is highly prevalent in many populations. Serogroups associated with invasive meningococcal disease (IMD) in Canada, for example, include A, B, C, W-135, X and Y. IMD is a rare but serious outcome of Nm infection, and can be prevented with vaccines that target certain serogroups. This has stimulated the development of dynamic models to evaluate vaccine impact. However, these models typically aggregate the various Nm serogroups into a small number of combined groups, instead of modelling each serogroup individually. The impact of aggregation on dynamic Nm model predictions is poorly understood. Our objective was to explore the impact of aggregation on dynamic model predictions.


We developed two age-structured agent-based models--a 2-strain model and a 4-strain model--to simulate vaccination programs in the Canadian setting. The 2-strain model was used to explore two different groupings: C, versus all other serogroups combined; and B, versus all other serogroups combined. The 4-strain model used the four groupings: C, B, Neisseria lactamica, versus all other serogroups combined. We compared the predicted impact of monovalent C vaccine, quadrivalent ACWY vaccine (MCV-4), and monovalent B vaccine (4CMenB) on the prevalence of serogroup carriage under these different models.


The 2-strain and 4-strain models predicted similar overall impacts of vaccines on carriage prevalence, especially with respect to the vaccine-targeted serogroups. However, there were some significant quantitative and qualitative differences. Declines in vaccine-targeted serogroups were more rapid in the 2-strain model than the 4-strain model, for both the C and the 4CMenB vaccines. Sustained oscillations, and evidence for multiple attractors (i.e., different types of dynamics for the same model parameters but different initial conditions), occurred in the 4-strain model but not the 2-strain model. Strain replacement was also more pronounced in the 4-strain model, on account of the 4-strain model spreading prevalence more thinly across groups and thus enhancing competitive interactions.


Simplifying assumptions like aggregation of serogroups can have significant impacts on dynamic model predictions. Modellers should carefully weigh the advantages and disadvantages of aggregation when formulating models for multi-strain pathogens.

Peer Review reports


Neisseria meningitidis (Nm) is a significant concern to public health. Prevalence of carriage varies between 5 and 10 % of the population [1]. Six main serogroups impact public health: A, B, C, W-135, X, and Y [1]. In Canada, serogroup B is most prevalent, although C, X, and Y are part of the ecological landscape as well, and are targeted by immunization programs at the time of writing [24]. Incidence of invasive meningococcal disease (IMD) in Canada ranges from 0.4 to 0.7 per 100,000 individuals [2], which is small compared to the disease burden in some countries, such as the “meningitis belt” of sub-Saharan Africa.

Nm infection is usually commensal: the pathogen typically colonizes the mucosal layer of the nasopharynx asymptomatically. Nm can be spread through the aerosol route, including coughing and sneezing, but also through close contact [2]. Therefore, closed environments, such as schools, households, and workplaces yield conditions conducive to transmission [510]. IMD incidence tends to be highest in children under 2 years of age, while Nm carriage tends to be highest in adolescent teenagers and young adults (17–21 years of age).

IMD symptoms develop when the colony penetrates the mucosal layer into the bloodstream, leading eventually to meningitis septicaemia [11]. Targeting the meninges and causing IMD, N. meningitidis can cause permanent neurological damage or fatality if not treated promptly. Many who survive IMD suffer loss of limbs or hearing. Determinants of individual susceptibility to IMD are not well understood.

There are several routes to natural immunity to Neisseria meningitis. After being infected with any serogroup of Nm, many individuals produce anti-bodies that prevent invasion by all serogroups of the pathogen, reducing the risk of becoming infected, although the duration of immunity is unknown, and probably relatively short [12]. Another route to transient Nm immunity is contracting Neisseria lactamica (Nl), a commensal pathogen. Though it is known that Nl colonizes the same loci as N. meningitidis, the mechanism and nature of the immunity is under study and currently not well known, though ideas of creating a vaccine using Nl have been proposed [13, 14]. N. lactamica colonization induces Immunoglobin A and Immunoglobin G antibody production, preventing recolonization for at least 24 weeks [15]. Though the mechanism is not clear, it is known that an immunological response to Nl causes an immunological response to Nm.

N. meningitidis and N. lactamica compete with one another for the same loci in the nasopharynx, as do the various N. meningitidis serogroups. If one pathogen is introduced to a host while the other has already colonized, the invading pathogen will not easily colonize due to competition for nutrients or adhesion [13].

Currently, vaccines reduce carriage and IMD rates in vaccinated individuals, and also unvaccinated individuals through herd immunity [16]. Previous polysaccharide vaccines were relatively ineffective, and could cause side effects [17]. However, in 2001 a conjugate vaccine that targets serogroup C was recommended for use in Canada for infants and children [18]. IMD incidence attributable to serogroup C was greatly reduced amongst the age groups that were vaccinated, as well as those that were not vaccinated [18]. However, with the reduced incidence of IMD attributable to serogroup C (and presumably, carriage as well), a small gap has been left in the ecological niche: there is less asymptomatic carriage of serogroup C in the nasopharyx, which, in principle, means more opportunities for other serogroups to colonize. This has caused some concern over potential vaccine-induced strain replacement [19]. We speculate that this ‘ecological’ strain replacement should be even stronger for vaccines against very common serogroups, such as serogroup B, if serogroups compete for space in the nasopharynx.

In 2007, Canada’s National Advisory Committee on Immunization recommended a quadrivalent conjugate Nm vaccine, meningococcal conjugate vaccine (MCV-4), that immunizes against serogroups A, C, W-135, and Y. A routine vaccine program was recommended for adolescent individuals in areas that local epidemiology warranted the prevention of the outbreaks, while high risk individuals were also recommended for vaccination [20].

Serogroup B causes significant IMD in many countries, since it dominates the ecological landscape of its environment [21]. Vaccines against serogroup B have been difficult to produce because the capsular polysaccharides for serogroup B are too similar to human neural antigens [22, 23]. Previous attempts required up to 20 different types of the antigen, PorA [24]. Though there have been difficulties, a vaccine has been developed using only four different components, 4CMenB [25]. With the use of 4CMenB in addition to MCV-4, it has been predicted that more than three quarters of the various Nm serogroups in Europe can be covered [26]. It is also predicted that 4CMenB can be used to target certain strains of serogroup X, which no vaccine is currently able to cover [27]. 4CMenB is one step closer to a universal vaccine for Nm, since the components found in the vaccine can also be found in all other capsule groups, thus 4CMenB may have an impact on other serogroups, although perhaps less efficaciously than the MCV-4 vaccine.

Previous models describing the impact of immunization on population health have aggregated all Nm serogroups into only two or three groups [2831]. Some models aggregated all serogroups of Nm while modeling the effects of Nl [28]. Other models have analyzed a single vaccine applied to different age categories [2, 30, 31], while some models compare vaccine types [2]. These models can produce valuable insights, and aggregating serogroups is a necessary model simplification. However, aggregating serogroups can also influence model projections because the ecological effects of serogroup structure and serogroup interactions cannot be fully accounted for [32, 33], and it is known that multi-strain models can produce qualitatively different dynamics from single-strain models [3438]. To our knowledge no research has directly compared the results of few-strain versus many-strain dynamic Nm models.

Here we develop, analyse and compare two age-structured, agent-based transmission models. The first model aggregates all serogroups into one of two categories while the second model aggregates all serogroups into one of four categories. Our objective is to compare the predicted impacts of monovalent C, monovalent B, and quadrivalent CAWY vaccines on Nm carriage using the two models, in order to learn more about potential biases introduced by using few-strain models to model multi-strain infectious diseases such as Nm.


4-strain model

Our 4-strain model utilizes four groupings: (1) serogroup B, (2) serogroup C, (3) serogroups A, W-135, X, Y, 29e, and ‘ungroupable’ (UG) serogroups, and (4) Neisseria lactamica. We include Neisseria lactamica since it competes for the same loci on the nasopharynx as Neisseria meningitidis serogroups [13]. Serogroups A, W-135, X, and Y are rare in Canada but 29e and the ungroupable serogroups are more prevalent. Hence, in the 4-strain model, serogroup B, serogroup C, and Nl form their own compartments, whereas other serogroups are lumped into the same compartment with one or more other serogroups.

The population is divided into annual age cohorts: <1 year old, 1 year olds, 2 year olds, 3 year olds, …, 99 year olds. Individual ages are updated monthly, moving the individual into the next age cohort after 12 simulated months. Individuals are removed each month according to all-cause, age-specific mortality rates [39] (Table S1 in Additional file 1). Individuals are born into the population at a constant rate, η, each month. A fine age structure enables comparing scenarios that distinguish between vaccinating at 12 months versus 24 months, for instance. However, for model calibration, broader age categories of 0–4, 5–9, 10–14, 15–19, and 20+ years of age were used, because 5-year age intervals are used in most available Nm carriage prevalence data as well as in the available contact data [40].

Disease transmission includes a baseline transmission rate specific to each age group and strain, βi,k, that is modified by an age specific contact matrix, Ci,j [40]. The contact matrix includes both the physical and non-physical contact of individuals as Nm is transmitted through aerosol means. Newborns are susceptible to all serogroups of Nm and to Nl. The probability per timestep that an individual of age j is infected by serogroup k is given by

$$ p\left({I}_{i,k}\right)=1-{\displaystyle {\prod}_{j=0}^{100}\left(1-\left(\frac{\beta_{i,k}{C}_{i,j}{I}_{j,k}}{N_j}\right)\right)} $$

where I i,k ,is the number of individuals of age i infected with strain k (carrying the strain asymptomatically), and N j is the number of individuals of age, j. Upon infection, individuals become a carrier of the kth serogroup.

No individual can be infected with more than one serogroup simultaneously. After the duration of carriage, τ 1k , has elapsed the individual becomes naturally immune to serogroup k that infected them, for a short duration τ 2k sampled from a gamma distribution. We assume that the average duration of carriage is the same for all ages and Nm serogroups. The average duration of carriage for Nl is also assumed to be the same for all ages, but differs from the average duration of carriage of Nm [13, 41]. While natural immunity to serogroup k, individuals also have some partial cross-protection against other serogroups, so that the probability per timestep that an individual of age j is infected by some other serogroup k is given by

$$ p\left({I}_{i,k}\right)=1-{\displaystyle {\prod}_{i=0}^{100}\left(1-\left(\frac{\beta_{i,k}{C}_{i,j}{I}_{j,k}}{N_j}\right)\left(1-{\sigma}_l\right)\right)} $$

where, 0 < σl <1 is the cross-immunity conferred by previous infection by serogroup l ≠ k. Cross-immunity lasts the same period of time as strain-specific immunity. Details on all parameter values appear in Additional file 1: Table S1.

2-strain model

The 2-strain model aggregates Nl and Nm serogroups, but it uses two groupings instead of four. The first version of the 2-strain model uses the groups: C versus Nl and all other Nm serogroups aggregated (A/B/W-135/X/Y/29e/UG/Nl). The second version uses the groups B versus Nl and all other Nm serogroups aggregated (A/C/W-135/X/Y/29e/UG/Nl). The first version is used to simulate the impact of monovalent C vaccine and compare it to the simulated impact of the monovalent C vaccine using the 4-strain model. Likewise, the second version is used to simulate the impact of monovalent B vaccine and compare it to the simulated impact using the 4-strain model. The 2-strain model is otherwise identical to the 4-strain model. In both versions, the duration of carriage and immunity was assumed to be the duration of Nm and was set by sampling from a gamma distribution.

Parameterization and uncertainty analysis

We used Canadian demographic and epidemiologic data to parameterize the model. The baseline natural history and demographic parameters were taken directly from the literature, while the baseline transmission rates βi,k were calibrated so that the age-stratified serogroup-specific carriage prevalence in the model matched the empirical data, within specified ranges of acceptability (Figure S1 in Additional file 1 and Table S2 in Additional file 1 for the 4-strain model) [2, 3]. In particular, initial values of βi,k were assumed for each simulation; every 100 simulated years, the code checked whether the modelled age-specific carriage prevalence was within the range of empirical acceptability based on seroprevalence data, for all of the ranges (20 for the 4-strain model, 10 for the 2-strain model); if so, those calibrated values of βi,k were used in the simulation, and if not, the values of βi,k were adjusted upward or downward as appropriate and the process was repeated for another 100 simulation years until the target ranges were attained.

The degree of cross-protection is not well known, so we used values from an in-vitro cross-reactivity study [42] (see Table S1 in Additional file 1 for baseline parameter values; see Tables S2, S3, S4 for the average βi,k values in Additional file 1). 50 realizations for each vaccine scenario were thereby produced, each through this calibration procedure. 50 realizations were used because this was found to be a sufficient number of realizations to produce representative averages of the (stochastic) agent-based model, at baseline parameter values. All results reported are the averages and standard deviations of the 50 realizations except where otherwise noted.

The calibration targets for the 4-strain model were obtained by adding the carriage prevalence of the constituent serogroups (Table S5 in Additional file 1). Similarly, the calibration of the 2-strain model, the acceptability ranges for prevalence of carriage were different for each age group and serogroup (Tables S6, S7 in Additional file 1).

In the case of the 4-strain model, a significant number of parameter sets yielded dynamics where serogroup prevalence oscillated over time, indicating the likely presence of multiple attractors. These parameter sets were not included in the baseline analysis, but are discussed separately at the end of the Results section.

Vaccine scenarios and assumptions

We explored nine vaccination scenarios with the 4-strain model (Table 1) and six scenarios with the 2-strain model (Table 2). Vaccination begins at t = 150 years (which provides enough “burn-in” to discard transients) and is continued for 75 years (tend = 225 years).

Table 1 Vaccine Scenarios for the 4-strain model
Table 2 Vaccine scenarios for the 2-strain model

For the 4-strain model, scenarios 13 use a monovalent C vaccine, scenarios 46 use a quadrivalent ACWY vaccine (MCV-4), and scenarios 79 use a monovalent B vaccine (4CMenB). For the 2-strain model, scenarios 13 use a monovalent C vaccine and scenarios 46 use the 4CMenB. Coverage of toddlers and infants was assumed to be 90 % and coverage for adolescents was assumed to be 80 %, based on Quebec coverage rates [38, 43]. The conjugate C and the MCV-4 vaccines were assumed to be 97 % efficacious in protecting inoculated individuals, while the B vaccine was assumed to be 64 % efficacious [44].

Vaccination operates in an “all-or-none” fashion, such that individuals who are efficaciously vaccinated receive full protection from infection and cannot transmit infection to others, while individuals who are not efficaciously vaccinated receive no protection and remain susceptible. Individuals are only protected against serogroups included in the vaccine (i.e., none of the vaccines, including 4CMenB, confer cross-immunity). Individuals are randomly selected each month to be vaccinated according to the vaccine coverage for their age group, and each vaccinated person is protected with a probability equal to the vaccine efficacy.


The overall prevalence of carriage decreases with the introduction of vaccination, regardless of the serogroup targeted by the vaccine or the age of immunization (Tables 3, 4, 5). Herd immunity was also observed for every vaccine scenario, since prevalence declined in age groups not targeted by the vaccine. We compare the 2-strain and 4-strain model predictions for each vaccine type in the following subsections.

Table 3 Reduction in prevalence of carriage with the Conjugate C vaccine
Table 4 Reduction in prevalence of carriage with the MVC-4 vaccine
Table 5 Reduction in prevalence of carriage with the 4CMenB vaccine

C Vaccine

In the 4-strain model, introducing the C vaccine causes the prevalence of serogroup C to decline quickly (Fig. 1a-c, see also Additional file 1: Figure S3a-c for time series plots with standard deviation bars on model outputs). The 2-strain model sees a more rapid reduction in the prevalence of carriage for serogroup C (Fig. 1d-f, see also Additional file 1: Figure S3d-f). In the 2-strain model, serogroup C prevalence is reduced to almost zero by approximately t = 165 years, while in the 4-strain model this occurs by approximately t = 170 years. After 15 years, the prevalence of carriage is reduced by 88.1 %, 95.0 %, and 97.3 % for the 4-strain model, compared to 96.3 %, 96.1 %, and 99.9 % for the 2-strain model, for vaccine programs 1–3 respectively (p = 0.0002, 0.0008, and 0.0001, respectively) (Table 3).

Fig. 1
figure 1

Prevalence of carriage with C vaccine. The prevalence of carriage under vaccine programs 1 (a,d,g), 2 (b,e,h) and 3 (c,f,i) for the 4-strain model (a-c, showing “C”, “B”, “Nl”, and “Other”), the 2-strain model (d-f, showing “C” and “Other”), and an overlay of model outputs for both models (g-i, showing “C” and “Other”). The initial carriage prevalence of serogroup C is different in the 2-strain and 4-strain models because of the filtering procedure used during model parameterization. The same plot with standard deviations of model outputs appears in Additional file 1: Figure S3

The 4-strain model predicts slight strain replacement caused by the C vaccine: the prevalence of “Other” increases slightly while the prevalence of Nl and B remain constant in vaccine programs 2 and 3, whereas the 2-strain model shows very little evidence of strain replacement, on account of the significant difference in prevalence between serogroup C and “Other”, which includes the very common species Neisseria lactamica. This has implications for disease control since different serogroups cause differing pathogenicity.

MCV-4 Vaccine

Vaccine programs 4 to 6 introduced a quadrivalent ACWY vaccine. In the 4-strain model, this reduced serogroup C as well as the “other” serogroup which partially includes A, W-135 and Y. Similar to the monovalent C vaccine, it takes about 20 years to reduce the prevalence of carriage serogroup C to almost zero. However, it requires an additional 10 years for the “other” serogroup to be reduced to the same level (Fig. 2a-c, see also Additional file 1: Figure S4a-c). The overall reduction in Nm prevalence was greater than for the C vaccine, on account of the MCV-4 vaccine including more serogroups. The MCV-4 vaccine causes strain replacement of serogroup B, but not of Nl, which is least affected by competition because it has a different duration of carriage than is typical for Nm serogroups.

Fig. 2
figure 2

Prevalence of carriage with MCV-4 vaccine. The prevalence of carriage with vaccine program 4 (a), 5 (b) and 6 (c) for the 4-strain model. The same plot with standard deviations of model outputs appears in Additional file 1: Figure S4

In the 4-strain model, the MCV-4 vaccine programs reduced the overall prevalence of carriage by 49.2 %, 48.7 %, and 49.9 % after 40 years, compared to the reduction of 7.0 %, 7.6 %, and 6.1 % under the C vaccine, for vaccine programs 1–3 respectively (p = <0.0001, <0.0001, and <0.0001, respectively) (Table 4). The reduction in overall prevalence of carriage is higher for the MCV-4 vaccine owing to its greater serogroup coverage than the C vaccine.

4CMenB vaccine

According to the 4-strain model, 4CMenB vaccine causes the prevalence of carriage to decline at a much slower rate than the decline in serogroups C and “Other” when C and MCV-4 vaccines were introduced. This difference is on account of the lower vaccine efficacy of 4CMenB (Fig. 3a-c, see also Additional file 1: Figure S5a-c). However, as was the case for the C and MCV-4 vaccines, the prevalence of serogroup B drops very close to zero after a sufficient amount of time.

Fig. 3
figure 3

Prevalence of carriage with 4CMenB vaccine. The prevalence of carriage with vaccine programs 7 (a,d,g), 8 (b,e,h) and 9 (c,f,i) for the 4-strain model (a-c, showing “C”, “B”, “Nl”, and “Other”), the 2-strain model (d-f, showing “B” and “Other”), and an overlay of model outputs for both models (g-i, showing “B” and “Other”). The initial carriage prevalence of serogroup “Other” is different in the 2-strain and 4-strain models because of the filtering procedure used during model parameterization. The same plot with standard deviations of model outputs appears in Additional file 1: Figure S5

In the 4-strain model, due to the high prevalence of serogroup B, the prevalence of serogroup C increases significantly when the 4CMenB vaccine is introduced on account of strain replacement, growing to a level comparable to that exhibited by serogroup B before the vaccine was introduced to the population (Fig. 3a-c, see also Additional file 1: Figure S5a-c). More interestingly, for vaccine programs 7 and 8, the prevalence of “Other” at first declines along with the prevalence of B, but then rebounds, and for vaccine program 9, the prevalence of “Other” declines steeply and remains low. This suggests that changes in the prevalence of serogroup B may be pushing the prevalence of “Other” to a basin of attraction for an alternative stable state (see next subsection). We also point out that any changes in serogroup C prevalence (caused by vaccine-induced changes to serogroup B prevalence) would also have additional impacts on “Other” serogroup prevalence, and hence the outcomes can be more complicated than would be observed with a fewer-strain model.

The 2-strain model also predicted a slower decline under the 4CMenB vaccine compared to the declines under the C and MCV-4 vaccine (Fig. 3d-f, see also Additional file 1: Figure S5d-f). However, as was observed with the C vaccine, the predicted decline in serogroup B under the 4CMenB vaccine is significantly faster under the 2-strain model than the 4-strain model (see also Fig. 3g-i, Additional file 1: Figure S5g-i). After 30 years of immunization, the 4-strain model predicts a reduction of 65.8 %, 61.7 % and 72.4 % in serogroup B for vaccine programs 7–9 respectively while the 2-strain model predicts a reduction of 93.4 %, 93.7 % and 98.5 % in serogroup B for vaccine programs 7–9 respectively, over the same 30-year time period (p = <0.0001, <0.0001, and <0.0001, respectively) (Table 5). Both 2-strain and 4-strain models predict that vaccinating at 12 months and 12 years of age causes the most rapid declines in prevalence (Table 5).

Unlike the 4-strain model, the 2-strain model predicts no strain replacement, which we again speculate is due to the much higher prevalence of the aggregated “Other” category, which includes Neisseria lactamica.

Apparent multiple attractors

During calibration of the 4-strain model it was noted that some parameter sets gave rise to sustained oscillations in serogroup prevalence. These were not included in the foregoing analysis. However, we present a few examples of these dynamics in Fig. 4. Although all serogroups oscillate in prevalence to some extent, the oscillations are particular pronounced for the “Other” serogroup. The C vaccine (programs 1–3, Fig. 4a-c) reduces the prevalence of C without having a significant impact on the oscillations in the other serogroups. The MCV-4 vaccine (programs 4–6, Fig. 4d-f). However, in some simulations (such as the one depicted in Fig. 4d), strain interactions drive serogroup B extinct even before any vaccines are introduced. For the 4CMenB vaccine, the elimination of serogroup B is accompanied by a shift in the dynamics of “Other” from sustained oscillations to a state of high, stable (non-oscillating) prevalence (Fig. 4g-i). The prevalence of C also decreases significantly, probably in response to competition from a surging prevalence of “Other”. The dynamics of a 4-strain model can be rich in ways that cannot be captured with a 2-strain model. We emphasize that these dynamics result from an interaction between forcing the system with the changes induced by the vaccine program, and the dynamical structure of systems with multiple attractors.

Fig. 4
figure 4

Evidence for multiple attractors in the 4-strain model. Examples of simulations from vaccine scenarios 1–9 (a-i) that entered the basin of attraction for oscillatory solutions. Some oscillations push serogroup prevalence to zero, causing that serogroup to go extinct in the simulation (no rescue effects were allowed)


Here we developed and compared 2-strain and 4-strain models of the impacts of several types of vaccines on Neisseria meningitidis serogroup carriage prevalence. The 2-strain and 4-strain models predicted similar overall impacts of vaccines on carriage prevalence, especially with respect to the vaccine-targeted serogroups. However, declines in vaccine-targeted serogroups were more rapid in the 2-strain model than the 4-strain model, for both the C and the 4CMenB vaccines.

Also, dynamical interactions that are present in the 4-strain model but not the 2-strain model contributed to other differences between the models. Sustained oscillations, and evidence for multiple attractors with differing basins of attraction, occurred in the 4-strain model but not the 2-strain model. In the case of the 4CMenB vaccine, the decline in serogroup B created an ecological niche, but because there are more than two other strains in a 4-strain model, it is not clear a priori whether serogroup C or “Other” would step in to fill the niche. In the model simulations, serogorup C was more successful in filling the niche, and the combined effect of changes in serogroup C and B prevalence caused the prevalence of “Other” to at first decrease, and then increase, or simply to decrease and remain low (Fig. 3d-f).

Strain replacement was more pronounced in the 4-strain model than the 2-strain model, we speculate on account of prevalence being more evenly among various serogroups in the 4-strain model, thus enhancing competitive effects. Neisseria lactamica, owing to its different natural history, did not experience the competitive effects as strongly as other Nm serogroups did from one another.

Oscillations have potentially important impacts on model calibration in multi-strain models. In our case, we excluded parameter sets that yielded oscillations in order to simplify the analysis. However, there are no empirical grounds for excluding such parameter sets, and oscillations in seroprevalence cannot be ruled out based on the existing and rather limited literature reporting cross-sectional results from various populations. Future models should permit oscillatory solutions that satisfy known empirical targets, especially since the qualitative and quantitative impacts of the vaccine could be very different for oscillatory solutions versus steady state solutions.

These types of dynamics appear exotic but they are not exclusive to mathematical models. There is evidence for such dynamics in other infectious disease systems, such as in pertussis before the introduction of vaccines [45]. Unfortunately, because high quality longitudinal data on Nm serogroup carriage prevalence do not exist, it is difficult to detect these dynamics in the case of Neisseria meningitidis. Moreover, although IMD is a notifiable disease, it is also highly stochastic, which could mask the signs of strain replacement on shorter timescales [2]. Finally, strain replacement effects can be delayed by a number of years, relative to the introduction of a vaccine program that causes them, and hence may not be immediately obvious in IMD case notification time series.

However, as more data on vaccine impacts on IMD incidence and serogroup prevalence become available over time, such new data can be used to further validate dynamic models, which can then be used to evaluate potential expansions of vaccine programs, or applications to other populations that do not currently have vaccine programs. These post-vaccine era data can—and should—be collected so that mathematical models can better inform policy recommendations. There is also value in calibrating the same model to different local circumstances. The additional validation obtained by comparing the model under different epidemiological, vaccine, and demographic circumstances can build confidence in the models so that they may better inform policy. Further work should also explore scenarios using multiple vaccines, since many jurisdictions are currently facing the choice of how best to implement vaccination programs against multiple Nm serogroups.


The phenomena explored in these models may have implications for the predicted effectiveness of Nm vaccination. Caution should be exercised when determining which type of model to use to simulate Neisseria meningitis dynamics. Models with few strains are easier to analyze, but may miss important features of disease dynamics.



Neisseria meningitidis


Invasive meningococcal disease


Quadrivalent ACWY vaccine


Monovalent B vaccine


Neisseria lactamica


  1. WHO | Meningococcal Meningitis []

  2. Vickers, DM, Anonychuk, AM, De Wals, P, Demarteau, N, Bauch, CT. Evaluation of serogroup C and ACWY meningococcal vaccine programs: projected impact on disease burden according to a stochastic 2-strain dynamic model. Vaccine 2014, in press.

  3. Christensen H, May M, Bowen L, Hickman M, Trotter CL. Meningococcal carriage by age: a systematic review and meta-analysis. Lancet Infect Dis. 2010;10(12):853–61.

    Article  PubMed  Google Scholar 

  4. PHAC Invasive Meningococcal Disease for Health Professionals []

  5. Caugant DA, Høiby EA, Magnus P, Scheel O, Hoel T, Bjune G, et al. Asymptomatic carriage of Neisseria meningitidis in a randomly sampled population. J Clin Microbiol. 1994;32(2):323–30.

    CAS  PubMed  PubMed Central  Google Scholar 

  6. Gauld JR, Nitz RE, Hunter DH, Rust JH, Gauld RL. Epidemiology of meningococcal meningitis at Fort Ord. Am J Epidemiol. 1965;82(1):56–72.

    Google Scholar 

  7. Hudson PJ, Vogt RL, Heun EM, Brondum J, Coffin RR, Plikaytis BD, et al. Evidence for school transmission of Neisseria meningitidis during a Vermont outbreak. Pediatr Infect Dis J. 1986;5(2):213–7.

    Article  CAS  Google Scholar 

  8. Munford RS, de Vasconcelos ZJS, Phillips CJ, Gelli DS, Gorman GW, Risi JB, et al. Eradication of carriage of Neisseria meningitidis in families: A study in Brazil. J Infect Dis. 1974;129(6):644–9.

    Article  CAS  PubMed  Google Scholar 

  9. Olcen P, Kjellander J, Danielsson D, Lindquist BL. Epidemiology of Neisseria meningitidis; prevalence and symptoms from the upper respiratory tract in family members to patients with meningococcal disease. Scand J Infect Dis. 1981;13(2):105.

    Article  CAS  PubMed  Google Scholar 

  10. Zangwill KM, Schuchat A, Riedo FX, Pinner RW, Koo DT, Reeves MW, et al. School-based clusters of meningococcal disease in the united states. JAMA. 1997;277(5):389–95.

    Article  CAS  PubMed  Google Scholar 

  11. Kremastinou J, Tzanakaki G, Pagalis A, Theodondou M, Weir DM, Blackwell C. Detection of igg and igm to meningococcal outer membrane proteins in relation to carriage of Neisseria meningitidis or Neisseria lactamica. FEMS Immunol Med Microbiol. 1999;24(1):73–8.

    Article  CAS  PubMed  Google Scholar 

  12. Troncoso G, Sanchez S, Moreda M, Criado MT, Ferreiros CM. Antigenic cross-reactivity between outer membrane proteins of Neisseria meningitidis and commensal neisseria species. FEMS Immunol Med Microbiol. 2000;27(2):103–9.

    Article  CAS  PubMed  Google Scholar 

  13. Gold R, Goldschneider I, Lepow ML, Draper TF, Randolph M. Carriage of Neisseria meningitidis and Neisseria lactamica in infants and children. J Infect Dis. 1978;137(2):112–21.

    Article  CAS  PubMed  Google Scholar 

  14. Gorringe A, Halliwell D, Matheson M, Reddin K, Finney M, Hudson M. The development of a meningococcal disease vaccine based on Neisseria lactamica outer membrane vesicles. Vaccine. 2005;23(17):2210–3.

    Article  CAS  PubMed  Google Scholar 

  15. Evans CM, Pratt CB, Matheson M, Vaughan TE, Findlow J, Borrow R, et al. Nasopharyngeal colonization by Neisseria lactamica and induction of protective immunity against Neisseria meningitidis. Clin Infect Dis. 2011;52(1):70–7.

    Article  CAS  PubMed  Google Scholar 

  16. Tzeng YL, Stephens DS. Epidemiology and pathogenesis of Neisseria meningitidis. Microbes Infect. 2000;2(6):687–700.

    Article  CAS  PubMed  Google Scholar 

  17. MacDonald NE, Halperin SA, Law BJ, Forrest B, Danzig LE, Granoff DM. Induction of immunologic memory by conjugated vs plain meningococcal c polysaccharide vaccine in toddlers. JAMA. 1998;280(19):1685–9.

    Article  CAS  PubMed  Google Scholar 

  18. MacDonald N. Immunization update 2005: Stepping forward. Can J Infect Dis Med Microbiol. 2005;16(4):219.

    PubMed  PubMed Central  Google Scholar 

  19. Kinlin M, Jamieson F, Brown EM, Brown S, Rawte P, Dolman S, et al. Rapid identification of herd effects with the introduction of serogroup c meningococcal conjugate vaccine in ontario, canada, 2000–2006. Vaccine. 2009;27(11):1735–40.

    Article  PubMed  Google Scholar 

  20. PHAC Statement on conjugate meningococcal vaccine for serogroups A, C, Y, and W-135 []

  21. PHAC The recommended use of the multicomponent meningococcal B (4CMenB) vaccine in Canada []

  22. Tan LKK, Carlone GM, Borrow R. Advances in the development of vaccines against Neisseria meningitidis. New Engl J Med. 2010;362(16):1511–20.

    Article  CAS  PubMed  Google Scholar 

  23. Sadarangani M, Pollard AJ. Serogroup B meningococcal vaccines—an unfinished story. Lancet Infect Dis. 2010;10(2):112–24.

    Article  CAS  PubMed  Google Scholar 

  24. Pillai S, Howell A, Alexander K, Bentley BE, Jiang HQ, Ambrose K, et al. Outer membrane protein (omp) based vaccine for Neisseria meningitidis serogroup b. Vaccine. 2005;23(17):2206–9.

    Article  CAS  PubMed  Google Scholar 

  25. Giuliani MM, Adu-Bobie J, Comanducci M, Aricò B, Savino S, Santini L, et al. A universal vaccine for serogroup B meningococcus. Proc Natl Acad Sci. 2006;103(29):10834–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Vogel U, Taha MK, Vazquez JA, Findlow J, Claus H, Stefanelli P, et al. Predicted strain coverage of a meningococcal multicomponent vaccine (4CMenB) in Europe: a qualitative and quantitative assessment. Lancet Infect Dis. 2013;13(5):416–25.

    Article  PubMed  Google Scholar 

  27. Hong E, Giuliani MM, Deghmane AE, Comanducci M, Brunelli B, Dull P, et al. Could the multicomponent meningococcal serogroup B vaccine (4CMenB) control Neisseria meningitidis capsular group X outbreaks in Africa? Vaccine. 2013;31(7):1113–6.

    Article  CAS  PubMed  Google Scholar 

  28. Coen PG, Cartwright K, Stuart J. Mathematical modelling of infection and disease due to Neisseria meningitidis and Neisseria lactamica. Int J Epidemiol. 2000;29(1):180–8.

    Article  CAS  PubMed  Google Scholar 

  29. Caro JJ, Möller J, Getsios D, Coudeville L, El-Hadi W, Chevat C, et al. Invasive meningococcal disease epidemiology and control measures: a framework for evaluation. BMC Public Health. 2007;7(1):130.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Trotter CCL, Andrews NJ, Kaczmarski EB, Miller E, Ramsay ME. Effectiveness of meningococcal serogroup C conjugate vaccine 4 years after introduction. Lancet. 2004;364(9431):365–7.

    Article  CAS  PubMed  Google Scholar 

  31. Trotter CL, Gay NJ, Edmunds WJ. Dynamic models of meningococcal carriage, disease, and the impact of serogroup C conjugate vaccination. Am J Epidemiol. 2005;162(1):89–100.

    Article  PubMed  Google Scholar 

  32. Stollenwerk N, Maiden MCJ, Jansen VAA. Diversity in pathogenicity can cause outbreaks of meningococcal disease. Proc Natl Acad Sci U S A. 2004;101(27):10229–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Van de Velde N, Brisson M, Boily MC. Understanding differences in predictions of HPV vaccine effectiveness: A comparative model-based analysis. Vaccine. 2010;28(33):5473–84.

    Article  PubMed  Google Scholar 

  34. Wearing HJ, Rohani P. Ecological and immunological determinants of dengue epidemics. Proc Nat Acad Sci USA. 2006;103(31):11802–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Poletto C, Meloni S, Colizza V, Moreno Y, Vespignani A. Host mobility drives pathogen competition in spatially structured population. PLoS Comput Biol. 2013;9(8):e1003169.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Ballesteros S, Vergu E, Cazelles B. Influenza A gradual and epochal evolution: insight from simple models. PLoS One. 2009;4(10):e7426.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Nuño M, Feng Z, Martcheva M, Castillo-Chavez C. Dynamics of two-strain influenza with isolation and partial cross-immunity. SIAM J Appl Math. 2005;65(3):964–82.

    Article  Google Scholar 

  38. Gumel AB. Global dynamics of a two-strain avian influenza model. Int J Comput Math. 2009;86(1):85–108.

    Article  Google Scholar 

  39. Statistics Canada Table 3–1 Deaths by single year of age and geography – Both sexes []

  40. Mossong J, Hens N, Jit M, Beutels P, Auranen K, Mikolajczyk R, et al. Social contacts and mixing patterns relevant to the spread of infectious diseases. PLoS Med. 2008;5(3):e74.

    Article  PubMed  PubMed Central  Google Scholar 

  41. De Wals P, Bouckaert A. Methods of estimating the duration of bacterial carriage. Int J Epidemiol. 1985;14(4):628–43.

    Article  PubMed  Google Scholar 

  42. Goldschneider I, Gotschlich EC, Artenstein MS. Human immunity to meningococcus: development of natural immunity. J Exp Med. 1969;129(6):1327–48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Bradet R, Audet D, Deceuninck G. Enquête sur la couverture vaccinale des enfants de 1 an et 2 ans au Québec en 2008. Institut national de santé publique du Québec, Direction des risques biologiques et de la santé au travail, 2009.

  44. An Advidory Committee Statement (ACS) National Advisory Committee on Immunization (NACI). Advice for the use of the Multicomponent Meningococcal Serogroup B (4CMenB) Vaccine. 2014.

    Google Scholar 

  45. Cherry J. The epidemiology of pertussis and pertussis immunization in the United Kingdom and the United State: A Comparative study. Curr Probl Pediatr. 1984;14(2):7–77.

    Google Scholar 

  46. White CP, Scott J. Meningococcal serogroup C conjugate vaccination in Canada: how far have we progressed? How far do we have to go? Can J Public Health. 2010;101(1):12–4.

    PubMed  Google Scholar 

Download references


The authors are grateful to support from the Natural Sciences and Engineering Research Council of Canada (NSERC) and the Canada Foundation for Innovation (CFI), and to Maria Major, Chiara Poletto, and Alberto Eugenio Tozzi for helpful comments.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Chris T Bauch.

Additional information

Competing interests

CTB has received research contracts from GlaxoSmithKline Vaccines for developing models of meningococcal vaccination.

Authors’ contributions

CTB conceived of the study. CTB and KP developed the model and wrote the manuscript. KP coded and simulated the model and analyzed the model outputs. Both authors read and approved the final manuscript.

Additional files

Additional file 1:

Supplementary Material. Additional file descriptions text (including details of how to view the file, if it is in a non-standard format).

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Poore, K.D., Bauch, C.T. The impact of aggregating serogroups in dynamic models of Neisseria meningitidis transmission. BMC Infect Dis 15, 300 (2015).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: