Modelling sexually transmitted infections: The effect of partnership activity and number of partners on
Introduction
In the present paper we investigate, by means of mathematical modelling, the potential spread of a sexually transmitted infection (STI) in a community of interest, with particular focus on heterogeneities in terms of number of sex partners and sexual activity within partnerships.
Similar questions have been addressed elsewhere in the literature using various approaches. A complete survey is beyond the scope of the present paper, but here are some examples. Inspired by Hethcote and Yorke (1984), May and Anderson (1987) acknowledge variations in sexual activity by dividing the community into subgroups, where the subgroup an individual belongs to is defined by the average number of new sex partners per unit of time s/he has. By approximating the initial stages of the outbreak by a set of differential equations they derive an expression for the basic reproduction number, an expression containing the variance of the number of new sex partners. In Diekmann et al. (1991) the forming and separation of partnerships is modelled dynamically in time—partnerships break up and new partnerships are formed with specified intensities. This leads to a Markov-type model and the number of partners during the infectious period will typically follow the Poisson distribution (or a sum of Poisson variables when there is more than one disease state). The initial phase of the disease outbreak is approximated by a suitable branching process assuming a large population and an expression for the basic reproduction number is obtained. A similar model, but where variation in sexual activity within partnerships is acknowledged by distinguishing between “steady” and “casual” partnerships, is analysed by Kretzschmar et al. (1994) and is derived. A somewhat different approach is taken by Diekmann et al. (1998) where each individual has k fixed “acquaintances” chosen uniformly within the community, and contacts between each pair of acquaintances occur independently and randomly in time. For such a model they derive expressions for but also the final size of the epidemic in case a major outbreak occurs. Eames and Keeling (2004) study a model in which each individual has at most one active partnership among a set of fixed potential partnerships but these active partnerships may change over time. The number of sexual partners over a given period can thus be fitted to empirical data. The spread of the disease in the dynamic network of active partnerships is approximated by differential equations relying on a large population and the moment closure method. Newman (2002) studies a general model for infectious diseases where the focus is in modelling social relationships by means of static network models (an approach used also in the present paper). A model for a heterosexually transmitted disease is also presented and briefly analysed. Nordvik and Liljeros (2006) look at an STI-model where the transmission probability depends on the number of sex acts, fits the model to the data reanalysed in the present paper, and derives an expression for the expected number of new infections caused by a randomly selected individual in the community.
In the present paper we study a model which takes into account both individual variation of number of sex partners and variation in sexual activity between partnerships. Empirical evidence shows that the distribution of the number of sexual partners, the degree distribution, is usually heavy-tailed (Colgate et al., 1989, Liljeros et al., 2001). For this reason our model allows an arbitrary degree distribution. Further, in the model we incorporate two types of partnerships, “steady” and “casual” (cf. Kretzschmar et al., 1994), having distinct transmission probabilities. This distinction refers to the number of sex acts over a period corresponding to the length of the infectious period. We distinguish between steady and casual (or high and low transmission risk) partnerships for two reasons. First, the number of sex acts affects the probability of disease transmission in such a way that the more sex acts the higher over-all probability of disease transmission (e.g. Rottingen and Garnett, 2002). Secondly, individuals having many sex partners tend to have fewer sex acts per partner compared to individuals with one or few partners (e.g. Giesecke et al., 1992, Blower and Boe, 1993, Nordvik and Liljeros, 2006).
Assuming a large population we approximate the initial phase of the epidemic by a suitable multitype branching process and derive an expression for the basic reproduction number which determines whether a major outbreak is possible or not. We then neglect the fact that there are steady and casual partnerships, simply treating all partnerships identically, and derive under this assumption. By calibrating parameters in the two models we can compare under the different assumptions. These questions are first addressed for a homosexual community (Section 2) and thereafter for a heterosexual community in which the degree distributions as well as transmission probabilities may differ between sexes (Section 3).
We model the sexual partnerships in the community by a static random network. It would of course be more realistic to allow new partnerships to be formed and old to break up, as is done in some of the papers cited above. However, if in reality new partnerships are formed and old break up in a time-stationary way, then the fixed set of partners of a given individual in the present model can be interpreted as the partners of that individual during a time period corresponding to the typical length of the infectious period in a dynamic network model. With this interpretation the present model can approximate a time-dynamic network model.
In Section 4 the heterosexual model is applied to data on sexual activity collected in the island of Gotland, Sweden (Giesecke et al., 1992, recently reanalysed in Nordvik and Liljeros, 2006). The data are from a representative community sample and contain the number of sex partners of the individuals and also information about the sexual activity in each partnership. is then computed numerically in terms of transmission parameters for two distinctions between steady and casual partnerships, as well as for the case when not distinguishing between partnerships. A comparison is made between the three partnership definitions. The comparison is “fair” because the transmission probabilities of the different partnership definitions are calibrated by computing them as means, within the specified partnership, of a simple per-sex-act transmission model. The main empirical finding from the analysis is that, for relevant transmission probabilities, is higher (and hence over-estimated) when neglecting differences in partnerships as compared to the case when partnerships are separated into steady and casual. The paper is concluded with a discussion on limitations of the present work and important future problems.
Section snippets
An STI model for a homosexual population
Consider a homosexual community of size n where we assume n to be large.
An STI model for a heterosexual population
Consider now a population consisting of males and females and assume that all sexual contacts are heterosexual.
Application to a study on sexual partnerships on the island of Gotland, Sweden
The calculations in this study are based on data from a study on sexual behaviour which was conducted in 1988 on Gotland, a Swedish island in the Baltic Sea (Giesecke et al., 1992).
Discussion
In the paper a model for the spread of an STI was defined allowing arbitrary partnership distribution and distinguishing between steady and casual partnerships. An approximation using branching processes was derived thus giving an expression for determining whether a major outbreak is possible or not. When fitted to data the main conclusion was that neglecting differences between partnerships has the effect that is systematically over-estimated.
The model can be made more realistic in
Acknowledgments
Tom Britton is supported by the Swedish Research Council and Fredrik Liljeros and Monica K. Nordvik are supported by the Swedish National Institute of Public Health and the European research NEST project DYSONET 012911. All three authors acknowledge Stockholm Group for Epidemic Modelling (www.s-gem.se).
References (23)
- et al.
The basic reproduction ratio for sexually transmitted diseases: I. Theoretical considerations
Math. Biosci.
(1991) - et al.
Monogamous networks and the spread of sexually transmitted diseases
Math. Biosci.
(2004) - et al.
The basic reproduction ratio for a sexually transmitted disease in a pair formation model with two types of pairs
Math. Biosci.
(1994) Epidemic models and social networks
Math. Scientist
(1999)- et al.(2000)
- et al.
Sex acts, sex partners, and sex budgets: implications for risk factor analysis and estimation of HIV transmission probabilities
J. AIDS
(1993) - Bollobás, B., Janson, S., Riordan, O., 2007. The phase transition in inhomogeneous random graphs. Rand. struct. alg.,...
- et al.
Generating simple random graphs with prescribed degree distribution
J. Stat. Phys.
(2006) - et al.
Risk behavior-based model of the cubic growth of acquired immunodeficiency syndrome in the United States
PNAS
(1989) - et al.
A deterministic epidemic model taking into account of repeated contacts between the same individuals
J. Appl. Probab.
(1998)
Sexual behavior related to the spread of sexually transmitted disease—a population-based survey
Int. J. STD AIDS
Cited by (21)
Case and partnership reproduction numbers for a curable sexually transmitted infection
2013, Journal of Theoretical BiologyCitation Excerpt :These studies did not, however, distinguish between the number of secondary infections and the number of secondary cases. Other studies have derived R0 using models that take into account partnership duration explicitly (Diekmann et al., 1991; Kretzschmar and Dietz, 1998; Diekmann et al., 1998; Britton et al., 2007) but these are for Susceptible-Infected (SI) infections (for example HIV) where an index case can infect the susceptible partner only once. In order to account for an infection in which an index case can infect a partner more than once we define the average number of secondary cases a primary case will infect during his or her infectious period as the case reproduction number Rc (Table 1).
Weighting for sex acts to understand the spread of STI on networks
2012, Journal of Theoretical BiologyCitation Excerpt :The number of RSA between two individuals i and j depends on the number of PSA that i attributed to j and vice versa (see the Methods section). Our framework stands out compared to earlier studies because we use a finite size network inferred from contact tracing data, contrary to earlier models that are based on infinite ‘scale-free’ networks (Joo and Lebowitz, 2004; Wang et al., 2007) or that use degree distributions without explicit network topology (Röttingen and Garnett, 2002; Britton et al., 2007). This allows us to model explicitly how individual sexual behaviours affect the network weighted topology and, hence disease spread.
Global asymptotic stability and a property of the SIS model on bipartite networks
2012, Nonlinear Analysis: Real World ApplicationsReverse-feeding effect of epidemic by propagators in two-layered networks
2015, Chinese Physics BImproving concordance in environmental epidemiology: A three-part proposal
2015, Journal of Toxicology and Environmental Health - Part B: Critical ReviewsMathematical Models for Estimating the Risks of Bovine Spongiform Encephalopathy (BSE)
2015, Journal of Toxicology and Environmental Health - Part B: Critical Reviews