Analyzing infectious disease data presents unique challenges for epidemiologists and biostatisticians. Unlike chronic diseases, for which a person's risk depends only on his or her personal exposures and risk factors, the risk that someone will acquire an infection is inherently dependent on whether others in the population are infected. As such, traditional statistical methods that assume outcomes are independent cannot always be applied to infectious disease data, and novel statistical methods are often needed. Furthermore, to understand the potential impact of interventions, one must account for the nonlinear feedbacks that give rise to population patterns of infection and disease. This requires blending epidemiologic methods with ecologic principles to develop models for transmission dynamics, which play an integral role in our understanding of infectious disease epidemiology (Figure).
An excellent illustration of the synthesis of data collection, analysis, and transmission dynamic models is provided by Lipsitch and colleagues.1 The authors analyze longitudinal data of Streptococcus pneumoniae nasopharyngeal carriage among children from Kilifi, Kenya. The study design and data collection themselves represent a remarkable feat, with baseline assessment of the carrier status of 2840 children and multiple follow-up assessments of 1868 children who were positive at baseline.1 – 3 Analyzing these data to obtain unbiased estimates of the serotype-specific rates of acquisition, clearance, and resistance to competition for 27 unique pneumococcal serotypes presents a considerable added challenge. In addition to individual risk factors such as age, the children's ability to clear colonization could depend on the resident serotype, how often they are exposed to someone shedding a different serotype, how resistant the resident serotype is to competition from other serotypes, and how good those serotypes are at competing. Thus, all the rates of interest are dependent on one another.
To address this challenge, Lipsitch et al apply a Markov transition model.1 This approach allows them to estimate the competing risks of clearance of the resident serotype and switching to each of the other serotypes. It also allows for estimation of the rate of acquisition for each of the 27 serotypes from the baseline prevalence distribution. This approach has considerable advantages over previous attempts to estimate serotype-specific parameters for S. pneumoniae, which were limited because they grouped vaccine serotypes and nonvaccine serotypes, estimated parameters for only a limited number of serotypes, assumed certain parameters were the same for all serotypes, or defined separate models for each serotype.3 – 9 The large longitudinal cohort study conducted in Kilifi presents a unique opportunity to estimate these parameters, providing data that are essential for understanding competition and coexistence among pneumococcal serotypes.
Complexities Underlying Pneumococcal Epidemiology
A key question for pneumococcal epidemiology is how vaccinations that protect against only certain serotypes will affect the overall diversity of serotypes in the population. Serotype replacement following the introduction of heptavalent pneumococcal conjugate vaccine (PCV7) into national immunization schedules has become a concern for prospects of controlling pneumococcal disease.10 Underlying this is the question of why so many serotypes of S. pneumoniae are able to coexist in the first place.
Transmission dynamic models can aid our understanding of the complex biologic processes that drive epidemiologic patterns such as the coexistence of pneumococcal serotypes, and they can provide an environment for testing hypotheses by simulating counterfactual scenarios. Models for the transmission dynamics are also increasingly used to inform policy decisions, such as whether to introduce a new vaccine into the national immunization schedule. However, to generate meaningful insight with real-world applications, such models must accurately portray the natural history of infection. Too often, transmission dynamic models are structured and parameterized based on few or no data, or on assumptions made by others, and then used to generate “predictions” about future trends in incidence or the impact of interventions. When this occurs, there is often a blurring of the distinction between what is an assumption and what is a prediction of the model. Similarly, models that rely too heavily on parameterizations obtained through fitting to observed incidence data as a substitute for a priori knowledge of the natural history are at best difficult to extrapolate to other settings and at worst may be misleading. The study by Lipsitch et al1 plays a crucial role in elucidating the natural history of various serotypes of S. pneumoniae and provides a basis for modeling the transmission dynamics. Consequently, their findings helped inform the parameter assumptions for an individual-based model for the transmission dynamics of the various pneumococcal serotypes.11
The coexistence of so many different pneumococcal serotypes is surprising, particularly given the evidence of Lipsitch et al1 that the fittest serotypes tend to be the best across all dimensions of fitness, ie, tend to have longer durations of infection and higher rates of acquisition, as well as to be more resistant to competition. By modeling the transmission dynamics and interactions of 25 pneumococcal serotypes, Cobey et al11 provided an explanation for this apparent paradox. They found that both weak levels of serotype-specific immunity generated by previous infection and a nonspecific decrease in the duration of colonization with each subsequent infection regardless of serotype (as observed by Lipsitch et al1) were necessary to generate patterns of diversity similar to those observed.11 The serotype-specific immunity helps the less-fit serotypes gain a competitive advantage among individuals previously infected by the more fit types, whereas the nonspecific immunity helps to equalize the fitness differences among serotypes by disproportionately penalizing the more fit types. Without the empirically derived characteristics of pneumococcal natural history obtained by Lipsitch et al,1 such a well-grounded understanding of the dynamics of coexistence would not have been possible.
Next Steps: Modeling Serotype Replacement
Transmission dynamic models are most useful when they are able to shed light on how incidence patterns depend on factors such as the transmission rate, which may be changing through time or as a result of interventions. The next step for pneumococcal epidemiology is to understand the patterns of serotype replacement that have emerged since the introduction of PCV7 vaccines in the United States and other countries, and to predict the potential impact of the new 13-valent PCV (PCV13). Although a few modeling studies have attempted to predict patterns of serotype replacement following vaccination,11 – 19 all of these have had shortcomings. Static models have attempted to predict which serotypes are expected to become more prevalent following vaccination based on their structural properties19 or observed changes in carriage and disease following PCV7 introduction (and assuming similar effects for PCV13),16,17 but these models do not explain why this replacement occurs or how the level of replacement is expected to scale with vaccination coverage. Dynamic models have been largely theoretical13,14 or have considered only vaccine versus nonvaccine serotypes, grouping together bacterial strains that may in fact be described by different natural histories.12,15 Although Cobey et al11 make some predictions about the potential impact of vaccination on the distribution and prevalence of serotypes associated with carriage, they do not attempt to extrapolate their findings to the incidence of invasive and noninvasive disease. Furthermore, they note that the precise impact of vaccination will depend on the strength of natural and vaccine-induced cross-immunity,11 which is currently not well understood. Additional epidemiologic studies are needed to assess how the correlates of fitness that determine the prevalence of carriage relate to those that determine a serotype's invasiveness.
Collecting, analyzing, and using data to inform models for the transmission dynamics and ultimately to provide guidance for policymakers is an iterative process (Figure). Transmission dynamic models are essential for understanding the expected indirect (ie, “herd immunity”) effects of vaccination, which can make a difference in whether a vaccine is deemed cost-effective. Public health officials also rely on dynamic models to interpret patterns such as changes in the age distribution of cases following vaccine introduction. Acknowledging the limitations of existing transmission dynamic models helps to identify key gaps in our understanding and highlight which data should be collected in future epidemiologic and laboratory studies. Once the data are collected, analysis using statistical methods to appropriately deal with the complexities of data about infectious diseases can lead to better parameter estimates. These parameters can then be used to refine dynamic models and further our understanding of disease transmission.
One of the overarching challenges in infectious disease epidemiology is to understand the comparative multistrain dynamics of various pathogens and how each is affected by vaccination. Serotype replacement is a concern not only for pneumococcal vaccination, but is also an issue, or potential issue, for other newly introduced vaccines (eg, rotavirus, human papillomavirus, and meningococcal vaccines), vaccines in development (eg, malaria, HIV, and dengue), and existing vaccines (eg, influenza). Integrating empirical data analysis with transmission dynamic models will enable us to understand the differences and similarities that drive the patterns of strain replacement occurring naturally and under the selective pressures of vaccination. Ultimately, this will help us to design better vaccines and more effective vaccination strategies for staying one step ahead of the pathogens, rather than always struggling to catch up.
We thank Bryan Grenfell for helpful comments.
1. Lipsitch M, Abdullahi O, D'Amour A, et al.. Estimating rates of carriage acquisition and clearance and competitive ability for pneumococcal serotypes in Kenya with a Markov transition model. Epidemiology. 2012;23:510–519.
2. Abdhullahi O, Karani A, Tigoi CC, et al.. The rates of acquisition and clearance of pneumococcal serotypes in the nasopharynges of children in Kilifi district, Kenya. J Infect Dis. In press.
3. Abdullahi O, Karani A, Tigoi CC, et al.. The prevalence and risk factors for pneumococcal colonization of the nasopharynx among children in Kilifi district, Kenya. PLoS ONE. 2012;7:e30787.
4. Auranen K, Mehtala J, Tanskanen A, S Kaltoft M. Between-strain competition in acquisition and clearance of pneumococcal carriage–epidemiologic evidence from a longitudinal study of day-care children. Am J Epidemiol. 2010;171:169–176.
5. Cauchemez S, Temime L, Valleron AJ, et al.. S. pneumoniae transmission according to inclusion in conjugate vaccines: Bayesian analysis of a longitudinal follow-up in schools. BMC Infect Dis. 2006;6:14.
6. Domenech de Celles M, Opatowski L, Salomon J, et al.. Intrinsic epidemicity of Streptococcus pneumoniae depends on strain serotype and antibiotic susceptibility pattern. Antimicrob Agents Chemother. 2011;55:5255–5261.
7. Erasto P, Hoti F, Granat SM, Mia Z, Makela PH, Auranen K. Modelling multi-type transmission of pneumococcal carriage in Bangladeshi families. Epidemiol Infect. 2010;138:861–872.
8. Melegaro A, Choi Y, Pebody R, Gay N. Pneumococcal carriage in United Kingdom families: estimating serotype-specific transmission parameters from longitudinal data. Am J Epidemiol. 2007;166:228–235.
9. Melegaro A, Gay NJ, Medley GF. Estimating the transmission parameters of pneumococcal carriage in households. Epidemiol Infect. 2004;132:433–441.
10. Weinberger DM, Malley R, Lipsitch M. Serotype replacement in disease after pneumococcal vaccination. Lancet. 2011;378:1962–1973.
11. Cobey S, Lipsitch M. Niche and neutral effects of acquired immunity permit coexistence of pneumococcal serotypes. Science. 2012;335:1376–1380.
12. Choi YH, Jit M, Gay N, et al.. 7-Valent pneumococcal conjugate vaccination in England and Wales: is it still beneficial despite high levels of serotype replacement? PLoS One. 2011;6:e26190.
13. Lipsitch M. Vaccination against colonizing bacteria with multiple serotypes. Proc Natl Acad Sci U S A. 1997;94:6571–6576.
14. Lipsitch M. Bacterial vaccines and serotype replacement: lessons from Haemophilus influenzae and prospects for Streptococcus pneumoniae. Emerg Infect Dis. 1999;5:336–345.
15. Melegaro A, Choi YH, George R, Edmunds WJ, Miller E, Gay NJ. Dynamic models of pneumococcal carriage and the impact of the heptavalent pneumococcal conjugate vaccine on invasive pneumococcal disease. BMC Infect Dis. 2010;10:90.
16. Rubin JL, McGarry LJ, Strutton DR, et al.. Public health and economic impact of the 13-valent pneumococcal conjugate vaccine (PCV13) in the United States. Vaccine. 2010;28:7634–7643.
17. Shea KM, Weycker D, Stevenson AE, Strutton DR, Pelton SI. Modeling the decline in pneumococcal acute otitis media following the introduction of pneumococcal conjugate vaccines in the US. Vaccine. 2011;29:8042–8048.
18. Van Effelterre T, Moore MR, Fierens F, et al.. A dynamic model of pneumococcal infection in the United States: implications for prevention through vaccination. Vaccine. 2010;28:3650–3660.
19. Weinberger DM, Trzcinski K, Lu YJ, et al.. Pneumococcal capsular polysaccharide structure predicts serotype prevalence. PLoS Pathog. 2009;5:e1000476.
ABOUT THE AUTHORS
VIRGINIA E. PITZER is a Postdoctoral Research Associate in the Department of Ecology and Evolutionary Biology at Princeton University, and will start as an Assistant Professor at Yale School of Public Health in July. Her research focuses on the transmission dynamics of imperfectly immunizing infections, particularly rotavirus. NICOLE E. BASTA is an Associate Research Scholar in Ecology and Evolutionary Biology at Princeton University and an Affiliate Investigator at the Hutchinson Cancer Research Center. Her research focuses on understanding vaccine effects, particularly for influenza and meningitis.