Open Access Research article

Cohort effects in dynamic models and their impact on vaccination programmes: an example from Hepatitis A

Arni SR Srinivasa Rao1, Maggie H Chen2, Ba' Z Pham5, Andrea C Tricco5, Vladimir Gilca4, Bernard Duval4, Murray D Krahn3 and Chris T Bauch1*

Author Affiliations

1 Department of Mathematics and Statistics, University of Guelph, Guelph, Canada

2 Toronto General Hospital and University Health Network, Toronto, Canada

3 Department of Health Policy, Management and Evaluation, University of Toronto and University Health Network, Toronto, Canada

4 Institut national de santé publique Québec, Québec City, Canada, and Research Centre, Centre hospitalier universitaire de Québec, Québec City, Canada

5 BioMedical Data Sciences, GlaxoSmithKline, Toronto, Canada

For all author emails, please log on.

BMC Infectious Diseases 2006, 6:174  doi:10.1186/1471-2334-6-174

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2334/6/174


Received:14 July 2006
Accepted:5 December 2006
Published:5 December 2006

© 2006 Srinivasa Rao et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Infection rates for many infectious diseases have declined over the past century. This has created a cohort effect, whereby older individuals experienced a higher infection rate in their past than younger individuals do now. As a result, age-stratified seroprevalence profiles often differ from what would be expected from constant infection rates.

Methods

Here, we account for the cohort effect by fitting an age-structured compartmental model with declining transmission rates to Hepatitis A seroprevalence data for Canadian-born individuals. We compare the predicted impact of universal vaccination with and without including the cohort effect in the dynamic model.

Results

We find that Hepatitis A transmissibility has declined by a factor of 2.8 since the early twentieth century. When the cohort effect is not included in the model, incidence and mortality both with and without vaccination are significantly over-predicted. Incidence (respectively mortality) over a 20 year period of universal vaccination is 34% (respectively 90%) higher than if the cohort effect is included. The percentage reduction in incidence and mortality due to vaccination are also over-predicted when the cohort effect is not included. Similar effects are likely for many other infectious diseases where infection rates have declined significantly over past decades and where immunity is lifelong.

Conclusion

Failure to account for cohort effects has implications for interpreting seroprevalence data and predicting the impact of vaccination programmes with dynamic models. Cohort effects should be included in dynamic modelling studies whenever applicable.

Background

Hepatitis A

Hepatitis A virus (HAV) is endemic in most countries [1-4], although infection rates continue to decline as hygiene and sanitation improve [5]. Infection is often subclinical in children but the majority of infected adults develop clinical symptoms, and often do not seek medical attention [5,6]. The severity of HAV infection increases with age and is particularly severe for those with underlying chronic liver disease [7-10]. The rate of mortality attributable to HAV increases by age from a 0.2% in symptomatic young adults to 3.9% in symptomatic adults over the age of 80 [11]. The clinical illness usually lasts 4 weeks and symptoms include nausea, loss of appetite, fatigue, fever, abdominal pain and jaundice [6].

HAV is transmitted by the fecal-oral route. Chronic infection does not occur [12], and immunity appears to be lifelong [7]. Traditionally, children have played a prominent role in transmission due to poor hygiene and the high probability of subclinical infection [13]. The reported incidence in Canada from 1980–1994 was twice as high in the 5–9 age groups as in other age groups [14], and once the reduced probablility of clinical infection in children is taken into account [6,15-19], true infection rates in children are seen to be significantly higher than in adolescents and adults. Approximately 25–50% of clinical cases have unidentified risk factors [20]. Many of these may have been contacts of children with subclinical infection.

In Canada, HAV is not endemic and transmission is primarily person-to-person [5] with rare foodborne outbreaks [21]. Risk factors for contracting HAV are often unidentified [20,21], but known risk factors include membership in high risk groups (including men who have sex with men, injection drug users and the homeless), low socioeconomic conditions [22], membership in socially-contained religious communities [23-25], and travel to endemic areas [26,27]. The latter is a particularly significant risk factor in Canada. Secondary infections from index cases infected abroad [25,28,29] and intra-household transmission [30] are also significant. From 1980 to 1994 (the year of vaccine licensure in Canada), the average incidence of reported cases was 6.3 per 100,000 per year [14] (Figure 1).

thumbnailFigure 1. Hepatitis A incidence in Canada (number of cases per 100,000 per year), 1978–2003.

Case reports reveal when clinical infections occur, but they do not account for subclinical infection and under-reporting [1,5]. By comparison, seroprevalence data can accurately determine what proportion of individuals have been infected by a given time, but cannot say when the infection occurred. Estimates of true infection rates usually rely upon seroprevalence surveys and/or case reporting data. Catalytic modelling uses integral equations to reconcile case reporting and seroprevalence data via statistical regression, yielding estimates of the true infection rates and how they have evolved over time. According to estimates from catalytic modelling, infections are under-reported by a factor of 10 in the United States [31,32].

Cohort effects

To date, there is no representative national seroprevalence survey for HAV in Canada. However, an approximate age-stratified seroprevalence profile has recently been synthesized from various independent surveys by the method of systematic review [33]. The resulting profile is both reasonable and consistent with the profile observed in the NHANES (National Health and Nutrition Examination Survey) II and III surveys in the United States [31]. The observed seroprevalence in Canada and the United States is too low in the youngest age classes and too high in older age classes to be consistent with the assumption of constant infection rates over time. Instead, it is necessary to posit that infection rates were higher in past decades due to a lower level of sanitation and hygiene [34]. This is usually described as a cohort effect [5]. Cohort effects have also been observed in other populations [34,35].

The presence of cohort effects in seroprevalence data is widely recognized. However it is difficult to characterize cohort effects without a mathematical model since the two dimensions involved – age class and calendar year – are easily confounded. Figure 2 illustrates how two cohorts born at different times experience different infection rate histories, when infection rates decline over time. The expected seroprevalence for a given cohort can be thought of as an integration of the infection rate history over the cohort's life history path in Figure 2. When infection rates have been constant over time, age-stratified seroprevalence profiles have a convex shape, increasingly monotonically with age but with progressively smaller increases for increasing age. However, when a strong cohort effect is present due to rapidly declining infection rates, an S-shaped age-stratified seroprevalence profile should be observed.

thumbnailFigure 2. Illustration of the cohort effect for an idealized situation of no age dependence in infection rates. The top arrow represents time evolution in the infection rate on a cohort born in 1920 and the bottom arrow represents the same for a cohort born in 1970. The infection rate is the rate at which a susceptible person becomes infected, and here it has been expressed in an idealized unit of measurement.

Goals of study

Few dynamic models have been developed for HAV, compared to other common infectious diseases [34,36-38]. Dynamic modelling studies and other studies which estimate infection rates for many diseases often do not take the cohort effect into account, resulting in inaccurate predictions. Here we develop and analyze a dynamic model for Hepatitis A in Canada. Our goal is to quantify how infection rates have declined over the past century, to illustrate how dynamic models which include the cohort effect can reconcile case reporting data and seroprevalence data, and to illustrate how the predicted impact of vaccination on incidence and mortality varies depending on whether the cohort effect is accounted for in the dynamic model.

Methods

Compartmental model

We develop an age-structured compartmental model to study HAV transmission and vaccination in Canada. The literature on age-structured compartmental models is well-established and has developed extensively since the mid-1980s [39]. Such models have long been applied to infectious diseases where natural immunity is lifelong and transmission is primarily person-to-person, as is the case for HAV in Canada. They have been validated against epidemiologic data [39-42] and have been used by health authorities in decision making processes for vacination policy [43]. Our SEIRV (Susceptible-Exposed-Infectious-Recovered-Vaccinated) compartmental model allocates all members of the population into mutually exclusive compartments according to age class and epidemiologic status: the number susceptible (Si), exposed (infected but not yet infectious, Ei), infectious (Ii), recovered (Ri), and vaccinated (Vi) in each age class i. The age classes are 0–4, 5–9, 10–19, 20–29, 30–39, 40–59 and 60+, chosen to reflect age categories in available sources of demographic and epidemiologic data.

Flow rates between compartments are defined by model parameters. Exposed individuals enter the infectious compartment at rate δ, and infectious individuals enter the recovered compartment at rate γi, thereafter retaining lifelong immunity. Individuals in age class i are vaccinated at per capita rate gi thereby entering the vaccinated compartment, and vaccinated individuals lose their immunity at per capita rate f, re-entering the susceptible compartment. Individuals are born susceptible at rate b, and when a given birth cohort enters the next highest age class, a proportion di of them die. Maternal immunity is short-lived and affects relatively few individuals in a non-endemic country such as Canada, so we do not include it [5]. HAV in Canada is spread primarily person-to-person [5]. Hence we do not model foodborne or waterborne outbreaks. The transmission rates and the model for the cohort effect are described in the following subsections. The model is parameterized using data from the clinical literature [44-50], demographic [51] and travel data [26,52-55], case reports [14], and seroprevalence suveys [33]. Model equations appear in 1 and the parameterization is described in 2. A diagram of the model appears in Figure 3.

Additional File 1. Model Equations. This file defines the model equations, parameters and variables.

Format: PDF Size: 114KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional File 2. Parameter Estimation. This file describes the model parameterization.

Format: PDF Size: 350KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

thumbnailFigure 3. Diagram of the SEIRV model.

Travel transmission rates

Estimates of the proportion of HAV cases attributable to travel in endemic countries range from 26% to 40%, depending upon the year and location [26,53-55]. A considerably higher proportion of cases are attributable to travel in younger age classes than in older age classes [55]. The predicted impact of vaccination, and the fit of the model to the data, can vary significantly depending on whether travel-related incidence is accounted for. Here we include travel-related incidence in our model.

If λi denotes the total infection rate in age class i (the rate at which a susceptible individual in age class i becomes infected), and if κi denotes the proportion of infections in age class i attributable to travel in endemic countries, then clearly τi = κi λi where τi is the travel transmission rate, ie, the rate at which a susceptible person in age class i becomes infected due to travel in an endemic country. Let <a onClick="popup('http://www.biomedcentral.com/1471-2334/6/174/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2334/6/174/mathml/M1">View MathML</a>i be the average travel transmission rate for 1980–1994. We estimate <a onClick="popup('http://www.biomedcentral.com/1471-2334/6/174/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2334/6/174/mathml/M1">View MathML</a>i from the average force of infection λi from 1980–1994 as estimated by catalytic modelling, and from average κi values for 1980–1994, estimated from surveys on HAV risk factors (2) [26,53-55].

Domestic transmission rates

Since the contribution to the infection rate from travel-related transmission is κi λi, the remainder (1 - κi)λi is the contribution from domestic transmission. The domestic transmission rate βij is the rate at which a susceptible individual of age class i in Canada is infected by infectious individuals of age class j in Canada. The resulting matrix of 7 × 7 domestic transmission rates form a WAIFW (Who Acquires Infection From Whom) matrix. By specifying a structure for the WAIFW matrix, the average values of λi and κi for 1980–1994 can be used along with estimated true incidence Ii and data from seroprevalence surveys [33] to determine the domestic transmission rates <a onClick="popup('http://www.biomedcentral.com/1471-2334/6/174/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2334/6/174/mathml/M2">View MathML</a>ij for 1980–1994 by solving the equation <a onClick="popup('http://www.biomedcentral.com/1471-2334/6/174/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2334/6/174/mathml/M3">View MathML</a> (2).

Modelling the cohort effect

We have estimated the average travel transmission rate <a onClick="popup('http://www.biomedcentral.com/1471-2334/6/174/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2334/6/174/mathml/M1">View MathML</a>i and the average domestic transmission rate <a onClick="popup('http://www.biomedcentral.com/1471-2334/6/174/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2334/6/174/mathml/M2">View MathML</a>ij for the years 1980–1994. The model could be simulated without including the cohort effect by assuming constant values βij = <a onClick="popup('http://www.biomedcentral.com/1471-2334/6/174/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2334/6/174/mathml/M2">View MathML</a>ij and τi = <a onClick="popup('http://www.biomedcentral.com/1471-2334/6/174/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2334/6/174/mathml/M1">View MathML</a>i for all time. To include the cohort effect, these constant values must be replaced by time-varying values βi = βi (t) and τi = τi (t).

To a good approximation, the per capita annual volume of travel by Canadian residents to HAV-endemic countries has increased quadratically since 1952, following the regression y = 2.38 × 10-5 t2 + 1.52 × 10-3 t + 2.98 × 10-3 (R2 = 0.98) [52]. Hence, it is possible to estimate a time-varying travel transmission rate, τi (t), by assuming that τi (t) increased over the past half century according to the same quadratic regression, rescaled to match the known average value <a onClick="popup('http://www.biomedcentral.com/1471-2334/6/174/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2334/6/174/mathml/M1">View MathML</a>i for 1980–1994 (2).

Additionally, the domestic transmission rate βi (t) can be assumed to obey βij (t) = <a onClick="popup('http://www.biomedcentral.com/1471-2334/6/174/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2334/6/174/mathml/M2">View MathML</a>ij F (t) = <a onClick="popup('http://www.biomedcentral.com/1471-2334/6/174/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2334/6/174/mathml/M2">View MathML</a>ij [L + H (1 - tanh(A(t - T)))], where tanh(t) = (exp(2t) - 1)/(exp(2t) + 1) is the hyperbolic tangent function. The function F(t) describes a continuous decline from an upper bound to a lower bound. The parameter L determines the lower bound attained for sufficiently large time t, L+H determines the upper bound attained for sufficiently small t, T determines the timing of the transition between these upper and lower bounds, and A controls how abrupt the transition is.

We estimated L, H, T and A by exploring a plausible region of LHTA parameter space with fine resolution and minimizing the least-squares error between predicted seroprevalence values and the observed values from the 11 seroprevalence surveys for Canadian-born individuals of distinct age classes and survey years, between 1980 and 1994 [27,33,56-58]. It was assumed that τi = τi (t) as estimated above. With this method we obtained L = 0.73, H = 0.84, A = 0.038, T = 1960.

A previous study has also captured cohort effects by fitting a function for declining infection rates to age-stratified seroprevalence data [34]. The main differences with the previous study are that the current study includes infection due to travel in endemic countries and also explores how including the cohort effect changes the predicted impact of vaccination.

Uncertainty analysis

To derive uncertainty intervals we applied the Latin hypercube method [59] to those parameters with the least certain values: the infection rates as estimated from catalytic modelling (λi), the duration of infectiousness (γi) and the proportion of travel-related incidence (κi). For each sampled parameter set, the model was fitted as before to the seroprevalence data to obtain values for L, H, T, and A and then seroprevalence and incidence were computed. Details appear in 3.

Additional File 3. Uncertainty Analysis. This file describes the methods of uncertainty analysis.

Format: PDF Size: 293KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Results

Declining transmissibility

If infection rates are assumed to be unchanged over the past century (ie, constant transmission rates <a onClick="popup('http://www.biomedcentral.com/1471-2334/6/174/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2334/6/174/mathml/M2">View MathML</a>ij and <a onClick="popup('http://www.biomedcentral.com/1471-2334/6/174/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2334/6/174/mathml/M1">View MathML</a>i are assumed) then the modelled seroprevalence profile exhibits a large error relative to the observed data from the 11 seroprevalence surveys for Canadian-born individuals (Table 1) [27,33,56-58]. The predicted seroprevalence is too large in the youngest age classes and too small in the oldest age classes. The predicted incidence is significantly larger than observed incidence for all age classes (Table 2).

Table 1. Observed seroprevalence versus fitted seroprevalence for the the case of constant domestic (βij =

<a onClick="popup('http://www.biomedcentral.com/1471-2334/6/174/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2334/6/174/mathml/M2">View MathML</a>

ij) and travel (τi =

<a onClick="popup('http://www.biomedcentral.com/1471-2334/6/174/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2334/6/174/mathml/M1">View MathML</a>

i) transmission rates, and for the case of time-varying domestic (βij = βij (t)) and travel (τi = τi (t)) transmission rates.

Table 2. Average observed incidence of reported cases (per 100,000 per year) versus predicted incidence of reported cases for the case of constant domestic (βij =

<a onClick="popup('http://www.biomedcentral.com/1471-2334/6/174/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2334/6/174/mathml/M2">View MathML</a>

ij) and travel (τi =

<a onClick="popup('http://www.biomedcentral.com/1471-2334/6/174/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2334/6/174/mathml/M1">View MathML</a>

i) transmission rates, and for the case of time-varying domestic (βij = βij (t)) and travel (τi = τi (t)) transmission rates, from 1980–1994.

By comparison, there is good agreement between fitted and observed seroprevalence when infection rates βij (t) and τi (t) evolve over time as described in Methods (Table 1). Accounting for the cohort effect through a declining transmissibility function explains most of the difference between fitted and observed seroprevalence (R2 = 0.97). The remaining discrepancies may partly be due to other heterogeneities not accounted for, such as outbreaks in high risk groups during specific periods, or constraints inherent to the hyperbolic tangent model. The fitted seroprevalence remains somewhat high in the younger age classes and somewhat low in the older age classes. The agreement between predicted and observed incidence also improves dramatically once the cohort effect is included (Table 2). Observed incidence values fall within 96% uncertainty intervals of the predicted values. A piecewise exponential function for F(t) yields similar results.

Figure 4 illustrates the time evolution of transmissibility F(t) over the past century from the fitted model. Transmissibility has dropped by a factor of 2.8 since the 1920s. The period of most rapid change was from 1950 to 1970, which is consistent with an era of rapid improvements in hygiene and sanitation in Canada, due ultimately to post-war economic growth.

thumbnailFigure 4. HAV transmissibility over time, as determined from fitting the dynamic model to seroprevalence data. The vertical axis shows the function F(t) which describes time evolution of transmissibility.

Universal vaccination

Here we compare the predicted effects of universal vaccination in Canada over 20-year and 50-year periods starting in 2006, with and without inclusion of the cohort effect in the dynamic model. For the case of constant transmissibility (cohort effect is not included), transmission rates in the vaccine era are assumed to be constant (βij = <a onClick="popup('http://www.biomedcentral.com/1471-2334/6/174/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2334/6/174/mathml/M2">View MathML</a>ij τi = <a onClick="popup('http://www.biomedcentral.com/1471-2334/6/174/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2334/6/174/mathml/M1">View MathML</a>i, as in Methods). For the case of declining transmissibility (cohort effect is included), transmission rates in the vaccine era are given by βij = βij (t), τi = τi (t), the same functions used to fit the model to the seroprevalence surveys but extrapolated to the vaccine era. The effects of the targeted policy in place from 1995 to the present are accounted for by reducing the travel transmission rate until the reported age-specific incidence from 1995–2005 is recovered, and thereafter the targeted policy is discontinued. Other parameters are unchanged from 1980–1994 values. Vaccination is applied to 4-year-olds and calls for two doses of vaccine with 70% compliance and 97% vaccine efficacy [49,50]. Vaccination occurs at age 4 because this would allow Hepatitis A vaccination to be efficiently scheduled along with other currently-administered vaccines in Canada. The assumed rate of waning of vaccine-derived immunity was 1.65% per year [47], which is a conservative assumption since cell-mediated immunity may actually result in lifetime protection.

Over both 20- and 50-year periods, failure to include the cohort effect over-estimates incidence and mortality, both with and without vaccination (Tables 3, 4). If the cohort effect is not included, the average reported incidence over 20 years is predicted to be 7.8 per 100,000 per year. If the cohort effect is included, it is 5.8 per 100,000 per year. Likewise, the cumulative number of deaths attributable to HAV over 20 years is predicted to be 262 instead of 138. Hence, from the perspective of absolute incidence in the vaccine era, not including the cohort effect will under-estimate the effectiveness of vaccination.

Table 3. Predicted annual incidence of reported cases, with and without including the cohort effect in the dynamic model.

Table 4. Predicted cumulative number of deaths attributable to HAV, with and without including the cohort effect in the dynamic model.

On the other hand, from the perspective of percentage reductions in incidence and mortality brought on by vaccination, the effectiveness of vaccination can be significantly over-estimated if the cohort effect is not included. In this case the percentage reduction in incidence (respectively mortality) due to vaccination is 43% (respectively 83%) over a 20-year period if the cohort effect is not included. These figures are 23% (respectively 73%) when the cohort effect is included. Patterns are similar over 50 years. The drop in mortality brought on by vaccination is largest in the oldest age classes since older individuals tend to be at a higher risk for mortality attributable to HAV.

The absolute and relative differences between incidence and mortality with and without vaccination, over 20-year and 50-year periods, are also larger when the cohort effect is not included.

We note that the average annual reported cases and deaths are somewhat higher over the 50-year period than the 20-year period (Tables 3 and 4). This occurs because the travel transmission rate after 2006 continues to climb according to our fitted function τi = τi (t) that reflects long-term historical trends in travel. Hence, the number of individuals infected through travel to HAV-endemic countries continues to rise, while the fall in the domestic transmission rate βij = βij (t) after 2006 is much more moderate. Another possible reason for this effect is that a decrease in incidence in the younger, recently-vaccinated age classes is offset by an increase in incidence in the older age classes, who are increasingly susceptible as more time passes since the cessation of the (all-ages) targeted policy in 2006. We have assumed that vaccine-derived immunity wanes, and so not all individuals vaccinated as children retain their immunity when older. Older age classes also experience higher rates of mortality attributable to HAV, so the net number of deaths can increase over time.

Discussion

These results indicate that the presence of a cohort effect can significantly alter the impact of vaccination. The implicit assumption that infection rates have always been constant is frequently made in modelling studies, but it is incorrect for many common infectious diseases. Failure to include the cohort effect overestimates incidence and mortality both before and after introduction of universal vaccination, as well as the percentage reduction achieved by universal vaccination.

Cohort effects are not immediately obvious in case reporting data or seroprevalence survey data. In seroprevalence data, it is often difficult to disentangle the two dimensions – calendar year and age – that are present in any set of seroprevalence surveys (Figure 2). In case reporting data, there is not a simple one-to-one relationship between reported incidence and true infection rates, particularly in cases where the probability of clinically apparent infection varies with age. In some situations, a rise in reported incidence may actually signal reduced infection rates in the general population or in certain age cohorts. This may have happened in Canada in the last decades of the twentieth century. Reported incidence in Hepatitis A in Canada exhibits an increasing trend between the years of 1980 and 1996 (linear regression: y = 0.21x + 4.61, R2 = 0.24; the R2 value is low because of the recurrent outbreaks). The dynamic model during this time predicts an increase in reported incidence (linear regression: y = 0.18x + 5.92, R2 = 0.99) while simultaneously predicting a decrease in true (actual) incidence due to a declining transmission rate (Figure 4). This occurs because declining transmissibility yields a decrease in the infection rates and hence an increase in the mean age at infection. Because of age dependence in the probability of clinical infection, this translates into more reported HAV infections. However, in other situations, variation in reported incidence over time may also be due to changing case reporting methods.

Because of such subtleties, dynamic modelling provides a rigorous and transparent framework for the interpretation and synthesis of seroprevalence and case reporting, and especially the clarification of cohort effects (eg, Figure 4). Modelling is valuable not only for predicting the impact of proposed policies but also for understanding and interpreting epidemiologic data. When developing dynamic models of disease transmission and vaccination, seroprevalence data should be interpreted carefully and incorporated into model parameterization whenever possible.

In our parameterization, some of the seroprevalence survey data were used twice. Firstly, they were used in the catalytic modelling to determine the proportion of seropositive individuals and hence the infection rates. Secondly, they were used to fit the function F(t) for declining transmissibility. This two-step approach may introduce bias, since it assumes transmissibility has declined in each age class in the same way over the past century. An alternative approach might have been to fit functions for declining transmissibility for each age class separately. However, this would mean estimating 10 parameters using 11 data points, which constitutes a poorly determined problem. In situations where more data points are available, this latter approach could be desirable.

In Canada, Hepatitis A transmission is primarily person-to-person. For other countries where environmental transmission is a significant source of Hepatitis A infection, the model would have to be modified to take environmental transmission into account [60].

Although this study addresses HAV specifically, the results should apply to any infectious disease in which there is lifelong immunity, person-to-person transmission, and where infection rates have changed significantly in past decades. Many pediatric infectious diseases fall into this category, with one difference being that disease severity declines with age instead of increasing, as with HAV. Hence, for many dynamic modelling studies, it may be necessary to account for cohort effects.

Conclusion

Failure to account for cohort effects has implications for interpreting seroprevalence survey data and parameterizing dynamic models. Moreover, for Hepatitis A, the predicted impact of universal vaccination can vary widely depending on whether the cohort effect is included in the dynamic model. This is likely to be true for any infectious disease with lifelong immunity in populations with pronounced cohort effects. Hence, cohort effects should be accounted for in dynamic modelling studies whenever applicable.

Competing interests

This study was funded by a Canadian Institutes of Health Research Industry-Partnered Operating Grant in connection with GlaxoSmithKline (GSK). BZP and ACT were employed by GSK. CTB was a consultant for GSK. CTB, BD and VG have received research funding from GSK.

Authors' contributions

All authors provided feedback on the modelling and the writing of manuscript as well as background to mathematical, epidemiological, clinical and/or public health literature. ASRSR programmed and helped construct the dynamic model, and wrote early drafts of the manuscript; MHC conducted the catalytic modelling under supervision of BZP; ACT reviewed the literature on Hepatitis A epidemiology; VG, BD and MDK provided guidance and feedback on epidemiological and public health aspects of HAV transmission and vaccination; CTB conceived of the study, constructed and analyzed the dynamic model, and completed the writing of the manuscript. All authors read and approved the final manuscript.

Acknowledgements

This study was funded by a Canadian Institutes of Health Research Industry-Partnered Operating Grant in connection with GlaxoSmithKline. The funding bodies had no role in the collection, analysis or interpretation of data, the dynamic modelling, the writing of the manuscript, or the decision to submit the manuscript for publication. The authors are grateful to Heather Whitaker and Paolo Bonanni for acting as reviewers for the manuscript.

References

  1. Tapia-Conyer R, Santos JI, Cavalcanti AM, Urdaneta E, Rivera L, Manterola A, Potin M, Ruttiman R, Tanaka Kido J: Hepatitis A in Latin America: a changing epidemiologic pattern.

    Am J Trop Med Hyg 1999, 61:825-829. PubMed Abstract | Publisher Full Text OpenURL

  2. Gust ID: Epidemiological patterns of hepatitis A in different parts of the world.

    Vaccine 1992, 10(Suppl 1):S56-S58. PubMed Abstract | Publisher Full Text OpenURL

  3. Cuthbert JA: Hepatitis A: old and new.

    Clin Microbiol Rev 2001, 14:38-58. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Wasley A, Samandari T, Bell BP: Incidence of hepatitis A in the United States in the era of vaccination.

    J Am Med Soc 2005, 294(2):194-201. OpenURL

  5. Scheifele DW: Hepatitis A vaccines: the growing case for universal immunisation of children.

    Expert Opin Pharmacother 2005, 6:157-164. PubMed Abstract | Publisher Full Text OpenURL

  6. Lednar WM, Lemon SM, Kirkpatrick JW, Redfield RR, Fields ML, Kelley PW: Frequency of illness associated with epidemic hepatitis A virus infection in adults.

    Am J Epidemiol 1985, 122:226-233. PubMed Abstract OpenURL

  7. Bell BP, Feinstone SM: Hepatitis A vaccine. In Vaccines. 4th edition. Edited by Plotkin SA, OWA. Philadelphia, PA: Saunders; 2004:269-298. OpenURL

  8. Koff RS: Seroepidemiology of hepatitis A in the United States.

    J Infect Dis 1995, 171(Suppl 1):S19-23. PubMed Abstract OpenURL

  9. Almasio PL, Amoroso P: HAV infection in chronic liver disease: a rationale for vaccination.

    Vaccine 2003, 21:2238-2241. PubMed Abstract | Publisher Full Text OpenURL

  10. Yao G: Clinical spectrum and natural history of viral hepatitis A in a 1988 Shanghai epidemic. In Viral Hepatitis and Liver Disease. Edited by Hollinger FB, Lemon SM, Margolis HS. Baltimore, MD, USA: Williams & Wilkins; 1991:76-78. OpenURL

  11. Jacobs RJ, Margolis HS, Coleman PJ: The cost-effectiveness of adolescent hepatitis A vaccination in states with the highest disease rates.

    Arch Pediatr Adolesc Med 2000, 154:763-770. PubMed Abstract | Publisher Full Text OpenURL

  12. Glickson M, Galun E, Orun R, Tur-Kaspa R, Shouval D: Relapsing hepatitis A. Review of 14 cases and literature survey.

    Medicine 1992, 71:14-23. PubMed Abstract OpenURL

  13. Advisory Committee on Immunization Practices (ACIP), Centers for Disease Control and Prevention: Prevention of hepatitis A through active or passive immunization.

    Morb Mort Wkly Rep: Recommendations and Reports 1999, 48:1-37. OpenURL

  14. Public Health Agency of Canada (PHAC): Notifiable Diseases Online. [http://dsol-smed.hc-sc.gc.ca/dsol-smed/ndis/c_time_e.html] webcite

  15. Benenson MW, Takafuji ET, Bancroft WH, Lemon SM, Callahan MC, Leach DA: A military community outbreak of hepatitis type A related to transmission in a child care facility.

    Am J Epidemiol 1980, 112:471-481. PubMed Abstract OpenURL

  16. Hadler SC, Webster HM, Erben JJ, Swanson JE, Maynard JE: Hepatitis A in day-care centers. A community-wide assessment.

    N Engl J Med 1980, 302:1222-1227. PubMed Abstract OpenURL

  17. Gingrich GA, Hadler SC, Elder HA, Ash KO: Serologic investigation of an outbreak of hepatitis A in a rural day-care center.

    Am J Public Health 1983, 73:1190-1193. PubMed Abstract | PubMed Central Full Text OpenURL

  18. Pavia AT, Nielsen L, Armington L, Thurman DJ, Tierney E, Nichol CR: A community-wide outbreak of hepatitis A in a religious community: impact of mass administration of immune globulin.

    Am J Epidemiol 1990, 131:1085-1093. PubMed Abstract OpenURL

  19. Staes CJ, Schlenker TL, Risk I, Cannon KG, Harris H, Pavia AT, Shapiro CN, Bell BP: Sources of infection among persons with acute hepatitis A and no identified risk factors during a sustained community-wide outbreak.

    Pediatrics 2000, 106(4):E54. PubMed Abstract | Publisher Full Text OpenURL

  20. Centers for Disease Control and Prevention (CDC):

    Hepatitis Surveillance Report. 2000., 57 OpenURL

  21. Public Health Agency of Canada (PHAC): Canadian Immunization Guide. [http://www.phac-aspc.gc.ca/publicat/cig-gci/] webcite

    6th edition. 2000.

  22. Harb J, Lem M, Fyfe M, Patrick D, Ochnio J, Dobson S, Hockin J: Hepatitis A in the Northern interior of British Columbia: An outbreak among members of a first nations community.

    Health Canada Canada Comm Dis Rep 26(19)

    1 October 2000

    OpenURL

  23. Hockin J: Hepatitis A outbreak in a socially-contained religious community in rural southern Ontario.

    Health Canada Canada Comm Dis Rep 23(21)

    1 November 1997

    OpenURL

  24. Naus M, Everett W, Davies S, Coutts J: A school outbreak of hepatitis A in southwestern Ontario.

    Can Dis Wkly Rep 15(45):225-228.

    1 November 1989

    PubMed Abstract OpenURL

  25. Deshaies D: Immunization against Hepatitis A during an Outbreak in a Jewish Orthodox Community-Quebec, 1997–1998.

    Health Canada Canada Comm Dis Rep 24(18)

    15 September 1998

    OpenURL

  26. De Serres G, Duval B, Shadmani R, Boulianne N, Pohani G, Naus M, Fradet MD, Rochette L, Ward BJ, Kain KC: Ineffectiveness of the current strategy to prevent hepatitis A in travelers.

    J Travel Med 2002, 9(1):10-16. PubMed Abstract OpenURL

  27. Ochnio JJ, Scheifele DW, Ho M: Hepatitis A virus in urban children – are preventative opportunities being missed?

    J Inf Dis 1997, 176:1610-1613. OpenURL

  28. Honish L: Hepatitis A Infected Food Handler at an Edmonton, Alberta Retail Food facility: Public Health Protection Strategies.

    Heath Canada CCDR 27(21)

    1 November 2001

    OpenURL

  29. Thompson L, McArthur B, Dalich-Adair T, Bell A: Hepatitis A outbreak in Grand Forks, British Columbia, 1990.

    Can Dis Wkly Rep 17(37):201-204.

    14 September 1991

    PubMed Abstract OpenURL

  30. Minuk GY, Ding LX, Hannon C, Sekla K: The risks of transmission of acute hepatitis A and B virus infection in an urban centre.

    J Hepatol 1994, 21:118-121. PubMed Abstract | Publisher Full Text OpenURL

  31. Armstrong GL, Bell BP: Hepatitis A virus infections in the United States: model-based estimates and implications for childhood immunization.

    Pediatrics 2002, 109:839-845. PubMed Abstract | Publisher Full Text OpenURL

  32. Muensch H: Catalytic models in epidemiology. Cambridge, Massachussetts: Harvard University Press; 1959. OpenURL

  33. Pham B, Duval B, DeSerres G, Gilca V, Tricco AC, Ochnio J, Scheifele DW: Seroprevalence of hepatitis A infection in a lower endemicity country: a systematic review.

    BMC Inf Dis 2005, 5:56-67. BioMed Central Full Text OpenURL

  34. Jacobsen KJ, Koopman JS: The effects of socioeconomic development on worldwide hepatitis A virus seroprevalence patterns.

    Int J Epidem 2005, 34:600-609. Publisher Full Text OpenURL

  35. Wong KH, Liu YM, Ng PSP, Young BWY, Lee SS: Epidemiology of Hepatitis A and Hepatitis E infection and their determinants in adult Chinese community in Hong Kong.

    J Med Virology 2004, 72:538-544. Publisher Full Text OpenURL

  36. Sattenspiel L, Simon CP: The spread and persistence of disease in structured populations.

    Math Biosci 1998, 90:341-366. Publisher Full Text OpenURL

  37. Diel R, Rappenhoener B, Schneider S: Cost-effectiveness of hepatitis A immunization of children and adolescents in Germany.

    Health Econ in Prev and Care 2001, 2:96-103. Publisher Full Text OpenURL

  38. Samandari T, Bell BP, Armstrong GL: Quantifying the impact of hepatitis A immunization in the United States, 1995–2001.

    Vaccine 2004, 22:4342-4350. PubMed Abstract | Publisher Full Text OpenURL

  39. Anderson RM, May RM: Infectious diseases of humans: dynamics and control. Oxford: Oxford University press; 1991. OpenURL

  40. Bauch CT, Earn DJD: Transients and attractors in epidemics.

    Proc R Soc Lond B 2003, 270:1573-1578. Publisher Full Text OpenURL

  41. Grassly NC, Fraser C, Garnett G: Host immunity and synchronized epidemics of syphilis across the United States.

    Nature 2005, 433:417-421. PubMed Abstract | Publisher Full Text OpenURL

  42. Babad HR, Nokes DJ, Gay NJ, Miller E, Morgan-Capner P, Anderson RM: Predicting the impact of measles vaccination in England and Wales: model validation and analysis of policy options.

    Epidemiol Infect 1995, 114(2):319-344. PubMed Abstract OpenURL

  43. Ministry of Health of New Zealand: Modelling Measles: Predicting and preventing measles epidemics in New Zealand. Application of a mathematical model. Wellington: Ministry of Health; 1998. OpenURL

  44. Goldrick BA: Another vaccine preventable disease continues to emerge.

    Am J Nursing 2004, 104:27-28. OpenURL

  45. Leach C: Hepatitis A in the United States.

    Pediatric Infec Dis J 2004, 23:551-2. Publisher Full Text OpenURL

  46. Craig AS, Schaffner W: Prevention of hepatitis A with the hepatitis A vaccine.

    New Engl J Med 2004, 350:476-481. PubMed Abstract | Publisher Full Text OpenURL

  47. Jacobs RJ, Gibson GA, Meyerhoff AS: Cost-effectiveness of Hepatitis A-B Vaccine Versus Hepatitis B Vaccine for Healthcare and Public Safety Workers in the Western United States.

    Inf Control and Hosp Epidem 2004, 25(7):563-569. Publisher Full Text OpenURL

  48. Centers for Disease Control and Prevention:

    Viral Hepatitis Surveillance Programme, 1983–1989. OpenURL

  49. Innis BL, Snitbhan R, Kunasol P, Laorakpongse T, Poopatanakool W, Kozik CA, Suntayakorn S, Suknuntapong T, Safary A, Tang DB: Protection against hepatitis A by an inactivated vaccine.

    J Am Med Assoc 1994, 271:1328-1334. Publisher Full Text OpenURL

  50. Werzberger A, Mensch B, Kuter B, Brown L, Lewis J, Sitrin R, Miller W, Shouval D, Wiens B, Calandra G: A controlled trial of a formalin-inactivated hepatitis A vaccine in healthy children.

    N Engl J Med 1992, 327:453-457. PubMed Abstract OpenURL

  51. Statistics Canada: Canadian Statistics Section, Vital Statistics. [http://www.statscan.ca] webcite

  52. Statistics Canada: CANSIM Table 427–0003 and historical supplement 66–201. [http://cansim2.statcan.ca] webcite

  53. Duval B, Boucher F, Dion J, Dupont M, Robinson E, Tremblay M, Venne S, Willems B: Controle de l'hépatite A par l'immunisation au Quebec (Final Report, in French), Quebec Immunization Committee, December 1997.

  54. el Saadany S, Gully P, Giulivi A: Hepatitis A, B, and C in Canada. Results from the National Sentinel Health Unit Surveillance System, 1993–1995.

    Can J Public Health 2002, 93(6):435-438. PubMed Abstract OpenURL

  55. Gosselin C, De Serres G, Rouleau I, Duval B, Shadmani R, Naus M, Ward BJ: Comparison of trip characteristics of children and adults with travel-acquired Hepatitis A infection.

    Ped Inf Dis J 2006, 25(12):1184-1186. Publisher Full Text OpenURL

  56. McFarlane ES, Embil JA, Manuel FR, Gorelick M: Prevalence of antibodies to hepatitis A antigen in patients attending a clinic for treatment of sexually transmitted disease.

    Sex Trans Dis 1980, 7:87-89. Publisher Full Text OpenURL

  57. Crewe MD, Embil JA, Garner JB: Prevalence of antibodies to hepatitis A virus in Nova Scotia children.

    Can Med Assoc J 1983, 25:1195-1197. OpenURL

  58. Payment P: Antibody levels to selected enteric virus in a French-Canadian population in the province of Quebec.

    Immunology and Inf Dis 1991, 1:317-322. OpenURL

  59. Blower SM, Dowlatabadi H: Sensitivity and uncertainty analysis of complex models of disease transmission: an HIV model, as an example.

    International Statistical Review 1994, 62:229-243. OpenURL

  60. Golinski MR, Barany E, Ballyk M: Ecological conditions that favor the evolution of intermediate-virulence in an environmentally transmitted parasite.

    Journal of Mathematical Biology 2005, 51:389-402. PubMed Abstract | Publisher Full Text OpenURL

Pre-publication history

The pre-publication history for this paper can be accessed here:

http://www.biomedcentral.com/1471-2334/6/174/prepub