Identifying the Age Cohort Responsible for Transmission in a Natural Outbreak of
Identifying the major routes of disease transmission and reservoirs of infection are needed to increase our understanding of disease dynamics and improve disease control. Despite this, transmission events are rarely observed directly. Here we had the unique opportunity to study natural transmission of Bordetella bronchiseptica – a directly transmitted respiratory pathogen with a wide mammalian host range, including sporadic infection of humans – within a commercial rabbitry to evaluate the relative effects of sex and age on the transmission dynamics therein. We did this by developing an a priori set of hypotheses outlining how natural B. bronchiseptica infections may be transmitted between rabbits. We discriminated between these hypotheses by using forceofinfection estimates coupled with random effects binomial regression analysis of B. bronchiseptica ageprevalence data from within our rabbit population. Forceofinfection analysis allowed us to quantify the apparent prevalence of B. bronchiseptica while correcting for age structure. To determine whether transmission is largely within social groups (in this case litter), or from an external group, we used randomeffect binomial regression to evaluate the importance of social mixing in disease spread. Between these two approaches our results support young weanlings – as opposed to, for example, breeder or maternal cohorts – as the age cohort primarily responsible for B. bronchiseptica transmission. Thus ageprevalence data, which is relatively easy to gather in clinical or agricultural settings, can be used to evaluate contact patterns and infer the likely agecohort responsible for transmission of directly transmitted infections. These insights shed light on the dynamics of disease spread and allow an assessment to be made of the best methods for effective longterm disease control.
Published in the journal:
. PLoS Pathog 6(12): e32767. doi:10.1371/journal.ppat.1001224
Category:
Research Article
doi: 10.1371/journal.ppat.1001224
Summary
Identifying the major routes of disease transmission and reservoirs of infection are needed to increase our understanding of disease dynamics and improve disease control. Despite this, transmission events are rarely observed directly. Here we had the unique opportunity to study natural transmission of Bordetella bronchiseptica – a directly transmitted respiratory pathogen with a wide mammalian host range, including sporadic infection of humans – within a commercial rabbitry to evaluate the relative effects of sex and age on the transmission dynamics therein. We did this by developing an a priori set of hypotheses outlining how natural B. bronchiseptica infections may be transmitted between rabbits. We discriminated between these hypotheses by using forceofinfection estimates coupled with random effects binomial regression analysis of B. bronchiseptica ageprevalence data from within our rabbit population. Forceofinfection analysis allowed us to quantify the apparent prevalence of B. bronchiseptica while correcting for age structure. To determine whether transmission is largely within social groups (in this case litter), or from an external group, we used randomeffect binomial regression to evaluate the importance of social mixing in disease spread. Between these two approaches our results support young weanlings – as opposed to, for example, breeder or maternal cohorts – as the age cohort primarily responsible for B. bronchiseptica transmission. Thus ageprevalence data, which is relatively easy to gather in clinical or agricultural settings, can be used to evaluate contact patterns and infer the likely agecohort responsible for transmission of directly transmitted infections. These insights shed light on the dynamics of disease spread and allow an assessment to be made of the best methods for effective longterm disease control.
Introduction
Containing and ultimately eliminating infectious disease remains a central goal for many animal and public health officials. Dissecting disease transmission – in terms of identifying the routes and potentially heterogeneous rates of disease spread [1] – is an essential step in devising or optimizing intervention strategies aimed at pathogen eradication [1], [2]. This is because heterogeneities in transmission that arise due to for example age or sexspecific differences among individuals [2], [3] can greatly affect invasion and eradication criteria [1]. Unfortunately, precise measurements of transmission remain elusive due to the immense difficulties associated with identifying the nature of a potential contact, the probability of infection given a contact [4] and important drivers of heterogeneities in transmission [2], [3], [5]. A key reason for these difficulties is that transmission events are rarely observed directly, with some notable exceptions [6].
One useful approach that can shed partial light on the transmission process is to measure the forceofinfection (FOI: λ), or the per capita conversion rate of susceptible hosts [7]. The simplest way to think about the FOI, is that over a short interval of time – say from time t to t+Δ – the probability that a disease negative individual becomes disease positive is λΔ. The most popular way to estimate λ is through use of the observed agespecific prevalence (or the proportion of individuals that are disease positive in a crosssectional sample), due to the ease with which it is measured in most populations [8]. Indeed, FOI estimates have been calculated from ageprevalence data for human, and to a lesser extent, wildlife infections [3], [9], [10], [11], [12], [13], [14], [15], [16], [17], [18], [19]. Estimating the FOI can help identify those ageclasses responsible for transmission and evaluate the relative effects of each group on overall transmission. Here we evaluate the relative effects of sex, age and social structure on the transmission dynamics of the respiratory pathogen Bordetella bronchiseptica within a commercial rabbitry of New Zealand White (NZW) rabbits. In doing so we illustrate how analysis of ageprevalence data can be used to estimate the agespecific FOI. The importance of social organization in B. bronchiseptica transmission is also considered. To test for litterbased transmission events – for example, from mother to offspring or between siblings – we checked for significant correlation among the fate of siblings by using a litterbased randomeffects binomial regression to estimate the importance of horizontal versus pseudovertical transmission [20], [21]. The statistical tools we employ here are general and can be applied to a range of directly transmitted medical and veterinary diseases to help shed light on the dynamics of disease spread and allow an assessment to be made of the best methods for effective longterm disease control.
The Bordetella genus contains three closely related gramnegative bacteria that cause respiratory infections in humans and other mammals [22]. Whereas B. pertussis and B. parapertussis largely infect humans and cause the acute respiratory disease known as whooping cough [23], B. bronchiseptica typically causes chronic infections in a wide range of mammals [24]. Indeed B. bronchiseptica infection is often endemic in agricultural settings – including commercial rabbitries [25], [26] – where rapid spread and persistent infection make it difficult to control [23]. Despite its widespread nature, there is a paucity of data describing the epidemiology of B. bronchiseptica in terms of both the main route(s) of, and likely cohort(s) responsible for disease transmission. As a respiratory infection, the major physical route of transmission is oralnasal via direct aerosol droplets [27], [28]. Based on the published literature [26], [27], [29], we propose a set of plausible routes of transmission within a commercial rabbitry would include:

the pseudovertical oralnasal route from mother to new born offspring before the age of weaning (4 to 5 weeks of age), with transmission occurring during suckling or grooming (figure 1a);

horizontal between sibling transmission when rehoused at age of independence (4 to 5 weeks until 5 to 6 months of age) via routine social behavior such as densharing or grooming (figure 1b);

horizontally at time of breeding (5 to 6 months onward: figure 1c)

constant from birth via environmental contamination (figure 1d).
These (not mutually exclusive) possible routes of transmission may result in the prevalence of infection changing with age in different ways which can be related to different underlying hazard models/FOI patterns (Figure 1). As we will explore below, each of these transmission possibilities translates into a specific prediction which can be tested using our statistical framework. To identify the parsimonious hypothesis, we applied a piecewise constant model for the agespecific FOI [9], [11], [15]. Since B. bronchiseptica is an endemic persistent infection [23], we used a catalytic framework (which assumes a oneway flow from susceptible to infected). The importance of sex and location (facility building) and time of sampling on FOI estimates was also determined. We considered the importance of social mixing and organization in B. bronchiseptica transmission using random effect logistic regression estimates to control for litter as confounding variable in transmission models. In parallel, we took a molecular epidemiological approach to investigate whether strainspecific differences existed in the epidemiological pattern of infection [30], [31].
Methods
Ethics statement
All protocols involving rabbits were approved by the Institutional Animal Care and Use Committee (IACUC) at the Pennsylvania State University according to the guidelines of the American Association for Laboratory Animal Science.
Rabbit hosts and facility
This study was conducted at a commercial rabbitry which breeds NZW rabbits. The rabbitry comprised of three separate animal breeding buildings (Table 1; buildings A  C). Background health checks – in the form of comprehensive monthly pathology reports testing for >17 pathogens – have been carried out since January 2003 (n = 2 to 4 rabbits/month/building). These reports show that B. bronchiseptica has been endemic in the rabbitry since testing began and that of the other pathogens screened, only nonpathogenic Eimeria species (intestinal coccidia) are occasionally isolated. Importantly, our rabbitry is Pasteurella multocida free – infection with this respiratory pathogen has long been associated with upper respiratory disease (URD) in rabbits [32] – with no URD reported in the last 30 years.
Kits were weaned at 4–5 weeks of age and ‘weanlings’ segregated by sex and cohoused in sibling pairs. Rabbits of good breeding stock were selected as ‘future breeders’ and housed in pairs. The remaining ‘stock’ rabbits were housed singly and sold at 8–10 weeks old. ‘Breeder’ rabbits were initially bred at 5 or 6 months of age for females and males respectively. Breeders were, housed individually and rebred when litters were weaned. For all rabbits included in this study, the date of sampling, building and rabbit identification number were ascertained along with nasal swab (BD sterile swab, product # 220518).
Sampling strategies
A total of eight sampling efforts were carried out from November 2006 to September 2008, culminating in the collective nasal swabbing of 602 rabbits total. All nasal swabs were streaked onto BordetGengou (BG) agar (Difco) containing 10% sheep's blood (Hema Resources) with 20 µg/mL streptomycin (Sigma) as soon as possible after collection and incubated at 35°C for ∼3–5 days. The sampling strategies were as follows:
Sampling strategy one: forceofinfection
To estimate the FOI, dateofbirth and sex information was obtained for each rabbit along with a nasal swab for culture from 4 crosssectional sampling efforts across two buildings (Table 1, buildings A & B) that comprised 214 rabbits; 73 does, 44 bucks, 32 male kits, 30 female kits, 16 female future breeders and 19 male future breeders. Litter information was only available for kits. Exact dateofbirth and kit sex information was not collected during sampling strategies two and three (see below). This data could not be obtained retrospectively from the rabbitry and hence prevented the inclusion of these rabbits into FOI analyses.
Sampling strategy two: siblingtosibling transmission
Nasal swabs of kits (and their does for inclusion in the maternal transmission analyses below) were taken before weaning (∼2 and 4weeks of age) to ascertain infection status and comprised 160 kits total. At weaning, kits were cohoused in sibling pairs based on their B. bronchiseptica infection status as follows; 23 pairs with one sibling positive and one negative for B. bronchiseptica and 54 pairs where both siblings were negative for B. bronchiseptica. As a control, six B. bronchiseptica negative siblings were housed solitarily. 4weeks later, nasal swabs were collected to calculate the number of disease conversion events: specifically, cages containing rabbits converting from B. bronchiseptica negative to positive were scored ‘1’, with no change in infection status scored ‘0’.
Sampling strategy three: maternal transmission
Nasal swabs of does and their kits taken at weaning (∼4 weeks old) across 3 rabbitry buildings (A, B & C) were used in these analyses and included data from 3 new sampling efforts (n = 208), as well as relevant doekit data from the FOI analysis (n = 106 rabbits) and siblingsibling analysis (n = 180) described above. Thus a total of 86 does and 408 kits were included in the maternal transmission analyses.
Bordetella diversity
We used MultiLocus Sequence Typing (MLST) analysis [33], [34] to determine the phylogenetic relationships among 90 B. bronchiseptica isolates from 4 sampling efforts across three rabbitry buildings, as previously described for Bordetella. Briefly, genomic DNA from each isolate was obtained using a DNAeasy Tissue Kit (Qiagen) and nucleotide sequences were determined for internal regions of seven housekeeping genes for all 90 isolates (see the Bordetella MLST database at http://pubmlst.org/bordetella). All alleles were double stranded sequenced at The Pennsylvania State University Genomic Sequencing Center and an allele number was assigned to each unique allele sequence. The combination of the allele numbers at the seven loci defines the sequence types (ST) or allelic profile of each strain [33], [35].
Statistical analysis
(i) ForceofInfection models
All analyses were performed in R version 2.7.0 (http://www.Rproject.org) and the fully annotated R code is available as Supporting Information (Text S1). There are several methods for establishing disease foothold in a population; one measure is the FOI (λ), or the per capita infection rate of susceptible hosts. For persistent infections the most popular way to estimate λ is through use of the observed agespecific prevalence, largely due to the ease with which it is measured in most populations. It should be noted that for nonpersistent fullyimmunizing infections, the FOI can be inferred from ageseroprevalence profiles [8]. Calculating the FOI from ageprofiled data is generally referred to as the ‘catalytic framework’ in mathematical epidemiology [7]. In our study, agespecific prevalence data was collected from rabbit nasal swabs and is intervalcensored infectiontime data such that each individual is either infected (Y = 1) or not (Y = 0) within a set interval of time.
To estimate the FOI, we make a number of necessary and common assumptions about our data: that no portion of the rabbit population is free from disease exposure; that a perfect test is used (i.e. in our case, that nasal swabbing always detects infection when present); that diseaseinduced mortality is negligible; and that infection is lifelong (or in the case of ageseroprevalence data, that immunity is lifelong). As B. bronchiseptica infections in our rabbit population are chronic and nonvirulent, we assume that the infection is irreversible – meaning that infection is lifelong – and that the infectioninduced mortality is negligible and can be ignored. For data sets which may violate any of the above assumptions, alternative methods [8], as well as extensions to existing models to incorporate infection processes such as death and periods of passive maternal protection are now available for estimating the agespecific FOI, λ(a) [15]. For a nonimmunizing persistent infection such as B. bronchiseptica, the agespecific prevalence, P(a), is defined by the differential equation (1)where (1P(a)) is the agespecific proportion of susceptible hosts. In epidemiology, age is considered one of the main risk factors for infection and disease prevalence often shows strong and distinct relationships with age [15]. To illustrate this, suppose individuals are born susceptible. Total time of disease exposure therefore increases with age and thus, disease prevalence appears to be directly linked to population agestructure. Our population is an animal breeding facility and so is skewed towards younger ageclasses, which suffer less diseaseexposure time and therefore possibly lower overall disease prevalence. Equation (1) leads to the following prediction for how disease prevalence should depend on age: (2)
Equation (2) presents a catalytic model [7] that incorporates how the cumulative FOI up to age a will have acted on susceptible hosts. The probability of past infection is thus the cumulative distribution function of the time to infection (or 1 – the survival function). The underlying idea here is that any given individual of age a will have experienced a cumulative FOI throughout its lifetime, represented by the integral in Equation (2). The exponential term thus represents the probability of being uninfected and its complement (Equation 2), is thus the probability of being infected by age a. The parameters of Equation (2) can be estimated from the observed ageprevalence data recording the change in function P(a) with age (see Text S1 for detailed R code).
We use the so called ‘piecewise constant’ model to estimate agedependent λ(a), which assumes a fixed FOI within predetermined age intervals which should be inspired by some prior knowledge of ageheterogeneities in the population of interest. For example, schoolage versus preschool children and measles dynamics [11]. We chose a 3 ageinterval model which corresponds to the hypothesized routes of transmission (see Figure 1) and has a comparable expected number of events in each interval (i.e. ∼ equal numbers of rabbits in each of the 3 ageclasses) which is important for model integrity [36]. It should be noted that in our initial explorations we used a flexible piecewise constant model for λ(a) with 8 ageintervals (in months) [0, 1), [1, 4), [4, 8), [8, 12), [12, 18), [18, 22), [22, 26) and [26+), but such a highly segmented model was illdetermined and is not reported here. The probability that a subject has converted by age a is P(a), and the probability that the subject has not converted is 1  P(a). We use maximum likelihood techniques to find values of FOI elements that best fit the agestructured prevalence data. The individual infection status is considered a binary outcome variable Y, with Y = 1 for infected and 0 for uninfectedindividuals and we assume a conditional Bernoulli likelihood for the binary outcome variable (i.e. disease incidence), where the probability, P, of an individual being infected at age a is:(3)where y_{j} denotes the infectious status of the jth individual. We can estimate constant λ(a) the agevarying FOI by integrating equation (2) and numerically maximize the likelihood of the coefficients according to equation (3). This is straight forward within the R statistical programming language (http://www.Rproject.org). The detailed numerical recipe is given in the R code Text S1. We use the quasiNewton method to estimate agevarying λ(a) by minimizing the negative loglikelihood. We computed standard errors using partial profile likelihood (see Text S1 and [37]).
For the hypothesized transmission routes (a)–(c) (see Figure 1) we predict the FOI curve to peak at young, intermediate and high ageclasses, respectively, whereas (d) should result in an ageinvariant curve. We tested these predictions using the FOI estimates and their associated standard errors and used Akaike's Information Criterion (AIC) and the ΔAIC statistic for model comparison [38].
(ii) Generalized linear models (GLMs) and generalized linear mixed models (GLMMs)
To test for evidence of significant siblingtosibling transmission we used a binomial regression to investigate whether the infection status of cohoused siblings could help explain disease transmission events (see Sampling Strategy Two in M&M for data collection details and Text S1 for R code). Transition of a previously uninfected animal to infected within a sibling pair during their time of cohousing (4 weeks) was merited a ‘disease conversion’ event and scored ‘1’, with no change in infection status scored a ‘0’.
The fate of siblings in the motheroffspring cages is unlikely to be independent (see Sampling Strategy Three). For example, if motheroffspring infection is important, all offspring will have a simultaneously increased risk during periods of enhanced maternal shedding. Indeed most litterbased transmission events should result in significant correlation among the fate of siblings. Testing for such correlation is therefore important in its own right. However, additionally nonindependence of siblings violates a basic assumption of generalized linear regression. We can deal with the statistical challenge of the nonindependence by modeling littermembership as a random effect in a generalized linear mixedmodel (see [21] for a general introduction to GLMM's for biologists). Specifically we used random effect binomial regression analysis (with a complementary loglog link) with litter as a random variable. The algorithm we used for fitting this model was Penalized quasilikelihood (PQL) as implemented in the glmmPQLfunction of MASS Rpackage [39] (see Text S1 for more details and R code). The betweenlitter random effect is assumed to follow a multivariate normal distribution. Qualitative differences due to age, sex and date of sampling were consistent across rabbitry buildings A–C and quantitative differences were controlled for by including facility building as a factor in our analyses. For all analyses, maximal models were first fit to the data and minimal models then obtained by removing nonsignificant terms, beginning with the interaction terms (with p = 0.05 as a threshold), to arrive at the final models for which parameter estimates are reported.
Results
Forceofinfection
Two different catalytic models were fitted to the data: a piecewise constant FOI with 3 ageintervals (Figure 2b), corresponding to the hypothesized routes of transmission outlined in the Introduction and Figure 1 and a constant FOI corresponding to the null hypothesis. The 3 ageinterval model fits the data better (having the lowest AIC value of 219.39) than the null hypothesis of a constant FOI (ΔAIC of model (c) versus (a) = 39.6). The model results show close correspondence between the observed and expected prevalence data (Figure 2a). Both data and model fits exhibit a rapid increase in prevalence during the first and second ageclasses (i.e. in rabbits up to 5 months; Figure 2a). During the first month of life, the estimated FOI is substantial (Figure 2b; FOI = 0.16 month^{−1}) and peaks in the second ageclass (Figure 2b; FOI = 0.20 month^{−1}), with the older age classes (from 5–30 months of age) having the lowest FOI estimates (Figure 2b; from 3.3×10^{−4} to virtually zero).
Next we examined whether gender differences existed for FOI estimates. No differences between sexes were found in the FOI estimates. In the younger ageclasses, the prevalence data and model estimates peaked in the second ageclass resulting in positive FOI estimates in young weaned kits (1 to 4 month olds: Figure 3a–d). In the older ageclasses, B. bronchiseptica prevalence asymptoted and subsequently fell for both sexes, with a concomitant decline in the FOI estimates (Figure 3a–d).
Sibling versus maternal routes of infection
To examine the likelihood of becoming infected from an infected sibling, we ran a binomial regression on the experimental siblingtosibling transmission experiment (see Sampling Strategy Two in M&M for details). Being cohoused with an infected sibling increased the probability of becoming B. bronchiseptica positive (Figure 4a: cohoused with infected sibling: Z = 2.42, p = 0.016), such that uninfected kits were 3.85 times more likely to become infected when they were cohoused with an infected compared to an uninfectedkit (Figure 4a: 95% C.I. for odds ratio 3.85: 1.29–11.46). None of the solitary Bordetellafree rabbits (housed alone in isolation) converted to diseasepositive during this time.
Using the maternal transmission data (Sampling Strategy Three in M&M for details), the importance of siblingtosibling versus motheroffspring routes of transmission was investigated. First, the data revealed substantial correlation (0.53) among the infection fate of siblings and a highly significant litterrandom effect (litter variance = 4.2±0.7), demonstrating the importance of withinlitter transmission. Although the prevalence of B. bronchiseptica was significantly higher in does compared to kits (Figure 4b; Z = 5.03 p<0.0001), having an infected mother did not significantly increase the probability of kits being infected (infected mother: Z = 1.74, p = 0.09). Nor was there any significant relationship between the litter random effect and the mothers prevalence status (Z = −1.05×10^{−15}, p = 1.0).
Bordetella bronchiseptica diversity
We used MLST analysis to characterize the relationship between 90 isolates collected from four sampling efforts across rabbitry buildings. All isolates were of sequence type (ST)14, which is a member of the B. bronchiseptica complex I [35]. Thus, one circulating strain appears to dominate in our rabbit population.
Discussion
This study demonstrates how FOI estimates coupled with random effects binomial regression analyses represent powerful tools for discerning between alternative modes of transmission for a directly transmitted pathogen. Specifically, our results support a role for siblingtosibling transmission among young weaned kits as a major route of B. bronchiseptica spread in the rabbit population studied (Figure 2). That the FOI reached a maximum value between 1 to 4 months of age – a time period when kits are rehoused in sibling pairs – followed by a sharp decline in the older ageclasses, is consistent with high betweensibling transmission in young weanlings (Figure 2), regardless of host sex (Figure 3). Results from the binomial regression analyses further support a major role for siblingtosibling transmission in driving B. bronchiseptica dynamics in the rabbitry; being cohoused with an infected sibling increased the risk of infection almost 4fold (Figure 4a). In comparison, the data did not support all other potential transmission routes; namely maternal, breeder or environmental routes. These insights shed light on the dynamics of disease spread and allow an assessment to be made of the best method(s) for effective longterm disease control, discussed more fully below.
A basic motivation for this study was to demonstrate how robust statistical tools can be used to disentangle routes and modes of transmission in humans and social animals from infectionatage data (within family groups), which is of broad medical, ecological and veterinary interest. The FOI analyses we present may have greatest application for analyzing disease dynamics in medical and agricultural settings because here one often has direct access to dateofbirth information, knowledge of the distinct mixing patterns over the lifetime of the host, as well as host infection status (for example, by detecting a serological response in the live animal, by the polymerase chain reaction (PCR) or by pathogen isolation). One complexity which often arises in analyses of medical and agricultural diseases is clustering in the data; hosts live in families, litters or herds and once an infection is introduced, hosts within that cluster have a higher instantaneous rate of becoming infected than those outside the cluster. Our use of random effect binomial regression analysis allows us to estimate the subjectspecific measure of the effect [20] and evaluate the importance of social mixing in disease spread. Thus, using the following protocol, the transmission dynamics of a range of directly transmitted infections can be analyzed by: (1) using the catalytic model and associated FOI analysis to determine the core susceptible ageclass(es); (2) using randomeffect binomial regression to inform on whether transmission is largely within the social group (family/litter/herd etc) or from an external social group; (3) carefully constructing transmission experiments, whenever possible, to test whether withingroup versus betweengroup individuals are the dominant source of infection.
What might explain heterogeneities in rabbit susceptibility to B. bronchiseptica infection; for example, the decline in B. bronchiseptica prevalence in olderage classes (in rabbits ∼20 months of age)? Betweenrabbit variation in protective antiB. bronchiseptica immunity – and hence resistance to infection – is likely to at least partly explain differences in host susceptibility to infection. Indeed, recent work has shown that the protective immune response against B. bronchiseptica varies between individual rabbits, with robust serum IgG detected in some hosts for up to 5 months post infection, which correlated with clearance from the respiratory tract [28]. Given the persistent nature of B. bronchiseptica infections in rabbits – infections of 5 months were routinely recorded [28] – and other mammals [23], the decline in prevalence we observe is unlikely to be driven by bacteria clearance and recovery. Rather, some level of enhanced immune protection in older ageclasses may be responsible for conferring some level of antibordetella resistance. Thus, the low attack rates (or number of reported cases per unit time in a given ageclass, divided by the number in that age class) in olderage classes likely reflect low proportions of rabbits susceptible to infection – i.e. immune, diseasenegative hosts – rather than a real decline in the rate at which susceptible rabbits acquire infection. In addition, betweenrabbit heterogeneities in protective anti B. bronchiseptica immunity might also help explain differences in rabbit susceptibility to infection in the maternal and cohoused sibling transmission studies reported here.
Is there any epidemiological support for the major route of B. bronchiseptica spread (siblingtosibling) identified using our statistical framework? B. bronchiseptica is known to pass efficiently and spread rapidly between populations of young weaned pigs [40], consistent with a siblingtosibling route for B. bronchiseptica transmission amongst young farmed animals. This would be particularly true in agricultural systems where an allin/allout (the facility is completely emptied and cleaned between groups of agematched animals which move together between phases of production) policy of animal breeding is not practised, as is the case in the rabbitry under study. However, that our FOI estimates were above 0.1 before 1 month of age suggests some maternal or environmental transmission is occurring in young weanlings and may be key to initiating the siblingtosibling transmission which follows. Indeed, a maternal route of transmission is thought initiate B. bronchiseptica infections in swine and rabbits [26], [27], but that infection only becomes endemic when passed horizontally between different batches of susceptible young [27]. Interestingly, the time when FOI values peaked in young weanlings, coincided with a period where maternal protection wanes in kits – antibodies against B. bronchiseptica decreased between 2 6 weeks of age in rabbits [25] – and could also contribute to increased susceptibility to infection observed in this age class. Thus, based on our findings and the published literature, we propose that the cycle of B. bronchiseptica infection in our rabbitry is maintained by a proportion of chronically infected breeder females and males (the infectious reservoir) with the majority of transmission occurring between young weaned siblings.
One important application for the analytical tools presented here is in the implementation of targeted disease control programs. Given that targeting those highrisk subgroups identified as playing key roles in transmission – rather than applying disease control measures randomly – is one efficient strategy to control disease [2], [6], a precautionary management approach might rely on the selective removal of infected weanlings to reduce siblingtosibling transmission. Selective removal of breeder animals – which may represent potential maintenance hosts for B. bronchiseptica – may also improve disease control by eliminating the infectious reservoir. Indeed, preemptive culling based on predetermined patterns of disease spread has been successfully used to combat the spread of footandmouth disease in cattle [41], [42]. The relationship between culling intensity and the resulting disease prevalence can be estimated when knowledge on population density and disease prevalence is available [43]. This allows estimates to be made regarding the level of culling needed to produce significant reductions in disease prevalence.
The analyses presented here can be applied to a range of medical and veterinary diseases to better understand the dynamics and mechanisms of disease spread, provided they are directly transmitted and induce lifelong immunity to reinfection. For example, the disease caused by mycobacterium – the etiological agent of tuberculosis in animals including bovine and humans – is largely directly transmitted, causes a subacute or chronic disease state which is irreversible [14], [15], [18] and can be routinely confirmed via culture, making it a tractable disease for application of FOI analyses. Indeed, the tools of infectious disease quantitative epidemiology have successfully been applied to further understand Mycobacterium bovis infection dynamics in wildlife population of badgers [44], ferrets [14] and bison [15], [45] and has shed light on likely patterns of mycobacteria transmission in the wild. However, these tools have not been used to the same effect in agricultural settings despite the debilitating effects of this disease and the potential to improve disease control therein. Other veterinary diseases which are tractable for this type of analyses include brucellosis, bovine herpes infection, classical swine fever, bovine mastitis and atrophic rhinitis in swine, to name but a few. Finally, the FOI model presented here can be extended to include diseases with reversion to nondiseased state or nonbenign diseases (i.e. associated with increasing death rate), or indeed to include a period where hosts are not exposed to infection (for example, when maternal antibodies are known to provide protection against specific diseases early in life) similar to a guarantee time in survival analysis (see Caley & Hone 2002 for examples of such extensions).
Our study has some limitations. Although the method we outline can clearly reveal the ageclass for which most of the new infection occurs, it cannot easily discern whether that infection is mainly within an ageclass versus from a different ageclass. However, once the high FOI ageclass is identified, careful design of transmission experiments could confirm the likely source of infection, and such studies are underway in our University. To control and possibly eradicate infectious diseases we need a better understanding of pathogen population dynamics and structure. Indeed, only when HIV population structure was understood did the requirement for a threecocktail HIV drug therapy become clear [46]. Knowledge of pathogen population structure is also needed to determine which diseaseassociated genes are under directional selection change. To this end we used MLST analysis to investigate whether strainspecific differences existed in the epidemiological pattern of infection [30], [31]. However, only one major circulating sequence type – ST14 – was identified in our rabbits regardless of rabbit age, sex or facility building. The dominance of ST14 across our facility may be due to the successful expansion of this single serotype over time. Alternatively, a limitation in sampling could have potentially biased our results; the sequence type of only 1 colony per swabbed plate (i.e. per rabbit) was determined at each sampling round. Therefore if the rabbit was colonized with multiple strains we most likely detected the dominant type (ST14). More intensive sequencing typing is required to test whether the lack of genetic variation we report is real and such studies are ongoing.
This study demonstrates the ease with which potential routes and reservoirs of infection can be discriminated amongst from ageprevalence data in medical, agricultural, and wildlife setting when we have access to fundamental ageprevalence data. Much remains to be done to achieve a better understanding of the complex dynamics of chronic infections and to extend this model to incorporate factors such as host immunity and parasite genetic variation.
Supporting Information
Zdroje
1. WoolhouseME
DyeC
EtardJF
SmithT
CharlwoodJD
1997 Heterogeneities in the transmission of infectious agents: implications for the design of control programs. Proc Natl Acad Sci U S A 94 338 342
2. AndersonRM
MayRM
1991 Infectious Diseases of Humans: Dynamics and Control: Oxford Oxford University Press
3. KlepacP
PomeroyLW
BjornstadON
KuikenT
OsterhausAD
2009 Stagestructured transmission of phocine distemper virus in the Dutch 2002 outbreak. Proc Biol Sci 276 2469 2476
4. McCallumH
BarlowN
HoneJ
2001 How should pathogen transmission be modelled? Trends Ecol Evol 16 295 300
5. PerkinsSE
CattadoriIM
TagliapietraV
RizzoliAP
HudsonPJ
2003 Empirical evidence for key hosts in persistence of a tickborne disease. Int J Parasitol 33 909 917
6. EamesKT
KeelingMJ
2003 Contact tracing and disease control. Proc Biol Sci 270 2565 2571
7. MuenchH
1959 Catalytic Models in Epidemiology. Harvard University Press
8. HensN
AertsM
FaesC
ShkedyZ
LejeuneO
Seventyfive years of estimating the force of infection from current status data. Epidemiol Infect 138 802 812
9. GriffithsDA
1974 A catalytic model of infection for measles. Applied Statistics 23 330 339
10. AndersonRM
MayRM
1985 Agerelated changes in the rate of disease transmission: implications for the design of vaccination programmes. Journal of Hygiene 94 365 436
11. GrenfellBT
AndersonRM
1985 The estimation of agerelated rates of infection from case notifications and serological data. J Hygiene 95 419 436
12. FarringtonCP
1990 Modelling forces of infection for measles, mumps and rubella. Statistics in Medicine 9 953 967
13. KeidingN
1991 Agespecific incidence and prevalence: a statistical perspective. Journal of the Royal Statistical Society, Series B 154 371 412
14. CaleyP
HoneJ
2002 Estimating the forceofinfection; Mycobacterium bovis infection in feral ferrets Mustela furo in New Zealand. J Animal Ecol 71 44 54
15. HeiseyDM
JolyDO
MessierFO
2006 The fitting of general ForceOfInfection models to wildlife disease prevalence data Ecology 87 2356 2365
16. CornellSJ
BjornstadON
CattadoriIM
BoagB
HudsonPJ
2008 Seasonality, cohortdependence and the development of immunity in a natural hostnematode system. Proc Biol Sci 275 511 518
17. CattadoriIM
BoagB
BjornstadON
CornellSJ
HudsonPJ
2005 Peak shift and epidemiology in a seasonal hostnematode system. Proc Biol Sci 272 1163 1169
18. GauthierDT
LatourRJ
HeiseyDM
BonzekCF
GartlandJ
2008 Mycobacteriosisassociated mortality in wild striped bass (Morone saxatilis) from Chesapeake Bay, U.S.A. Ecol Appl 18 1718 1727
19. OzgulA
OliMK
BolkerBM
PerezHeydrichC
2009 Upper respiratory tract disease, force of infection, and effects on survival of gopher tortoises. Ecol Appl 19 786 798
20. LairdNM
WareJH
1982 Randomeffects models for longitudinal data. Biometrics 38 963 974
21. BolkerBM
BrooksME
ClarkCJ
GeangeSW
PoulsenJR
2009 Generalized Linear Mixed Models: a Practical Guide for Ecology and Evolution. Trends Ecol Evol 24 127 135
22. BjornstadON
HarvillET
2005 Evolution and emergence of Bordetella in humans. Trends Microbiol 13 355 359
23. MattooS
CherryJD
2005 Molecular pathogenesis, epidemiology, and clinical manifestations of respiratory infections due to Bordetella pertussis and other Bordetella subspecies. Clin Microbiol Rev 18 326 382
24. GoodnowRA
1980 Biology of Bordetella bronchiseptica. Microb Rev 44 722 738
25. GlassLS
BeasleyJN
1989 Infection with and antibody response to Pasteurella multocida and Bordetella bronchiseptica in immature rabbits. Lab Anim Sci 39 406 410
26. DeebBJ
DiGiacomoRF
BernardBL
SilbernagelSM
1990 Pasteurella multocida and Bordetella bronchiseptica infections in rabbits. J Clin Microbiol 28 70 75
27. de JongMF
1992 Progressive atrophic rhinitis. In: Diseases of Swine 1: pp. 414435, Iowa State University Press, Ames, IA
28. PathakAK
CreppageKE
WernerJR
CattadoriIM
2010 Immune regulation of a chronic bacteria infection and consequences for pathogen transmission. BMC Microbiol 10 226
29. PathakAK
BoagB
PossM
HarvillET
CattadoriIM
in press Seasonal Breeding drives the incidence of a chronic bacterial infection in a freeliving herbivore population. Infect Epid
30. GuptaS
FergusonN
AndersonR
1998 Chaos, persistence, and evolution of strain structure in antigenically diverse infectious agents. Science 280 912 915
31. ThompsonRCA
2000 Molecular epidemiology of infectious diseases Arnold: London
32. WebsterLT
1924 The epidemiology of a rabbit respiratory infection Il. Clinical, pathological, and bacteriological study of snuffles. J Exp Med 39 843 856
33. JolleyKA
ChanMS
MaidenMC
2004 mlstdbNet  Distributed multilocus sequence typing (MLST) databases. BMC Bioinformatics 5 86
34. MaidenMC
BygravesJA
FeilE
MorelliG
RussellJE
1998 Multilocus sequence typing: A portable approach to the identification of clones within populations of pathogenic microorganisms. Proc Natl Acad Sci U S A 95 3140 3145
35. DiavatopoulosDA
CummingsCA
SchoulsLM
BrinigMM
RelmanDA
2005 Bordetella pertussis, the causative agent of whooping cough, evolved from a distinct, humanassociated lineage of B. bronchiseptica. PLoS Path 1 e45
36. FriedmanM
1982 Piecewise Exponential Models for Survival Data with Covariates. Annals of Statistics 10 101 113
37. BolkerBM
2008 Ecological Models and Data in R: Princeton University Press
38. BurnhamKP
AndersonDR
2002 Model Selection and MultiModel Inference; A Practical Information  Theoretic Approach. 2: SpringerVerlag, New York
39. VenablesWN
RipleyBD
1994 Modern applied statistics with Splus. Statistics and Computing
40. SmithIM
GilesCJ
BaskervilleAJ
1982 Immunisation of pigs against experimental infection with Bordetella bronchiseptica. Vet Rec 110 488 494
41. KeelingMJ
WoolhouseME
ShawDJ
MatthewsL
ChaseToppingM
2001 Dynamics of the 2001 UK foot and mouth epidemic: stochastic dispersal in a heterogeneous landscape. Science 294 813 817
42. WoolhouseM
DonaldsonA
2001 Managing footandmouth. Nature 410 515 516
43. DobsonA
MeagherM
1996 The Population Dynamics of Brucellosis in Yellowstone National Park. Ecology 77 1026 1036
44. AndersonRM
TrewhellaW
1985 Population dynamics of the badger (Meles meles) and the epidemiology of bovine tuberculosis (Mycobacterium bovis). Philos Trans R Soc Lond B Biol Sci 310 327 381
45. JolyDO
MessierF
2004 Factors affecting apparent prevalence of tuberculosis and brucellosis in wood bison. Journal of Animal Ecology 73 623 631
46. SimonV
HoDD
Abdool KarimQ
2006 HIV/AIDS epidemiology, pathogenesis, prevention, and treatment. Lancet 368 489 504
Štítky
Hygiena a epidemiologie Infekční lékařství LaboratořČlánek vyšel v časopise
PLOS Pathogens
2010 Číslo 12
Nejčtenější v tomto čísle
 Glacial Refugia in Pathogens: European Genetic Structure of Anther Smut Pathogens on and
 Structural and Functional Studies of Nonstructural Protein 2 of the Hepatitis C Virus Reveal Its Key Role as Organizer of Virion Assembly
 Molecular Architectures of Trimeric SIV and HIV1 Envelope Glycoproteins on Intact Viruses: StrainDependent Variation in Quaternary Structure
 HIV Capsid is a Tractable Target for Small Molecule Therapeutic Intervention