Next Article in Journal
Contrasting Patterns of Pomacea maculata Establishment and Dispersal in an Everglades Wetland Unit and a Central Florida Lake
Next Article in Special Issue
Age-Independent Adult Mortality in a Long-Lived Herb
Previous Article in Journal
Alterations of Arbuscular Mycorrhizal Fungal Diversity in Soil with Elevation in Tropical Forests of China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Lifetime Mortality Trajectories in Wildlife Disease Research: BaSTA and Beyond

by
Dave W. Hudson
1,
Richard Delahay
2,
Robbie A. McDonald
3,
Trevelyan J. McKinley
4 and
Dave J. Hodgson
1,*
1
Centre for Ecology and Conservation, University of Exeter, Penryn TR10 9FE, UK
2
National Wildlife Management Centre, Animal and Plant Health Agency, Sand Hutton YO41 1LZ, UK
3
Environment and Sustainability Institute, University of Exeter, Penryn TR10 9FE, UK
4
College of Engineering, Mathematics and Physical Sciences, University of Exeter, Penryn TR10 9FE, UK
*
Author to whom correspondence should be addressed.
Submission received: 30 July 2019 / Revised: 9 September 2019 / Accepted: 24 September 2019 / Published: 1 October 2019
(This article belongs to the Special Issue Bayesian Survival Trajectory Analysis in Wildlife)

Abstract

:
Wildlife hosts are important reservoirs of a wide range of human and livestock infections worldwide, and in some instances, wildlife populations are threatened by disease. Yet wildlife diseases are difficult to monitor, and we often lack an understanding of basic epidemiological parameters that might inform disease management and the design of targeted interventions. The impacts of disease on host survival are generally associated with age, yet traditional epidemiological models tend to use simplistic categories of host age. Mortality trajectory analysis provides the opportunity to understand age-specific impacts of disease and uncover epidemiological patterns across complete life histories. Here, we use Bayesian survival trajectory analysis (BaSTA) software to analyse capture-mark-recapture data from a population of wild badgers Meles meles naturally infected with Mycobacterium bovis, the causative agent of tuberculosis in badgers and cattle. We reveal non-constant mortality trajectories, and show that infection exaggerates an age-dependent increase in late-life mortality. This study provides evidence for actuarial senescence in badgers, a species previously believed to display constant mortality throughout life. Our case study demonstrates the application of mortality trajectory analysis in wildlife disease research, but also highlights important limitations. We recommend BaSTA for mortality trajectory analysis in epidemiological research, but also suggest combining approaches that can include diagnostic uncertainty and the movement of hosts between disease states as they age. We recommend future combinations of multi-state and multi-event modelling frameworks for complex systems incorporating age-varying disease states.

Graphical Abstract

1. Introduction

Investigating the epidemiology of any disease in wild populations is challenging because of the practical difficulties in monitoring both infections and their wild hosts [1], yet wild animals are important contributors to many emergent and widespread infections of humans and livestock worldwide [2]. Consequently, any improvements in our understanding of wildlife disease epidemiology can be beneficial to human health, animal welfare and productivity, as well as to biodiversity conservation [1,3]. Disease-related mortality is a critical parameter in any epidemiological model [4] and methods for its estimation, developed for the study of human populations, are now commonly used in studies of wildlife [5]. Disease can increase rates of host mortality [6,7] and, in extreme cases, can contribute to decline and risks of species extinction [8,9]. Disease-induced mortality is commonly modelled as a simple increase in otherwise fixed mortality rates at particular life stages [10,11] but the reality can be much more complex. Mortality trajectories reveal patterns in age-specific mortality that are often missed when using fixed rates in discrete age classes, yet the inclusion of these trajectories in epidemiological models is rare. The ability to accurately model factors that influence mortality is fundamental to any demographic investigation and the inclusion of mortality trajectories could offer greater precision whilst uncovering subtle variations in the patterns of age-specific mortality. Here we show that the inclusion of disease status, alongside full mortality trajectory analysis, can reveal disease-induced changes in the shape of lifetime schedules of mortality.
Predicted increases in the proportion of the global human population living over the age of 80 [12] have attracted substantial economic investment and scientific interest [13] in age-associated diseases in humans. Disease effects have been shown to vary with age in both human [14] and now non-human animals [15]. While such variation can have a direct impact on population size [16], empirical studies of age-specific causes of mortality in wild populations are scarce [17]. The decline of physiological function with age, i.e. senescence, can affect fecundity [18], morphological traits [19], behaviour [20], and physiology [21] but most commonly it refers to an increase in the rate of mortality with age, known as actuarial senescence [22]. A long-standing belief that environmental factors would result in an animal’s death long before the impact of senescence became manifest has now been shown to be false [23]. The emergence of several high-quality, long-term demographic studies, combined with theoretical and statistical advances, have since revealed senescence in many wild populations of various species, and have promoted the study of factors affecting senescence in the wild. Traditional approaches in wildlife epidemiology have focused on estimating fixed mortality parameters for a finite and predetermined number of age categories (e.g., cubs versus adults [24]), which results in the loss of information and risks unnecessarily coarse conclusions. Mortality trajectory analysis, which considers rates of mortality through the entire lifespan, is a means of mitigating some of these problems, but until recently has been difficult to employ with wildlife populations. This is because of several statistical challenges [25], as well as the requirement to monitor marked individuals through entire life histories [23]. Furthermore, the accurate interpretation of mortality trajectories and their associated parameters requires ecological understanding of a species’ behaviour and life strategies [26], as well as suitable model selection through comparisons of model fit [27]. If an inappropriate model is fitted to the data, then subsequent conclusions are often specious.
Several mathematical functions have been used to describe lifetime trajectories of age-specific mortality. Historically, the preferred choice was the Gompertz curve [28], which describes mortality as increasing exponentially with age from an intercept representative of an initial baseline risk of death. The more flexible Weibull distribution [29] has recently been used in survivorship analysis, and can accommodate accelerating increase, decelerating increase, decreasing, and constant mortality [30], but can fail to capture early decelerations in mortality rates. Makeham [31] proposed that death could be separated into age-dependent and age-independent sources, resulting in the addition of a “Makeham” term (a constant) to established mortality functions (e.g., Gompertz–Makeham, Weibull–Makeham, etc.). Criticism of these models has focused on the omission of individual heterogeneity, leading to potentially flawed statistical inference [32] and development of the logistic model has been proposed as a solution [33]. This model accommodates individual variation as a decrease in mortality at advancing ages (i.e. individual frailty; frail individuals tend to die younger, leaving more robust individuals as survivors) [34,35,36]. Although these functions have been applied across a variety of species [37,38], support for them is mixed and it is often difficult to deduce whether uncertainty in parameter estimates is evidence of absence (the parameter is not important) or absence of evidence (the data do not provide a sufficient signal to help infer the parameter). Life histories can be split into distinct phases, each with an individual mortality signature. Caughley [39] identified three stages: juvenile (characterised by relatively high mortality), adult (with relatively low mortality), and senescent (increasing rates of mortality in later life). When combined, these stages form a “U” shape or “bathtub” mortality curve, which can capture complex patterns of mortality over complete life histories (e.g. the Siler function [40]). The use of mortality trajectories and the separation of life stages allows us to unravel subtle changes in mortality and identify actuarial senescence, hereafter referred to as senescence. This combining of mortality models to accommodate different life stages has found support in a variety of species [41,42,43], and has helped to reveal some of the drivers [44] of age-specific patterns of mortality and senescent variation.
To evaluate the use of mortality trajectories and investigate patterns of age-specific mortality, we use data from a long-term study of a population of European badgers (Meles meles) naturally infected with Mycobacterium bovis, the causative agent of bovine tuberculosis (bTB). In the United Kingdom, badgers are the primary wild maintenance host for bTB [45], and in some locations contribute to the persistence of infection in cattle populations [46]. Being able to discern patterns of mortality and distinguish between groups in terms of disease susceptibility may help inform approaches to managing infection in badgers, thereby reducing risks of onward transmission to cattle. The epidemiology of bTB in badgers has been the focus of many studies (for a review, see [47]), yet the use of mortality trajectories, using age as a continuous variable, has to our knowledge been implemented only once. McDonald et al. [44] used Bayesian survival trajectory analysis (BaSTA) [48] to identify sex differences in mortality from the point of infection and infer the mechanisms underpinning them. Increased mortality as a result of bTB infection in male badgers is already known [11,24], but this was further dissected to identify where in the infection process the sex differences arose [44]. These analyses, however, do not capture the full lifetime trajectory of mortality in badgers, and cannot reveal actuarial senescence. Here, we study mortality trajectories by using cubs of known age and known infection status to analyse disease-related mortality trajectories across entire life histories. Previous analyses of badger mortality have suggested that badgers show no evidence of actuarial senescence [11], but these conclusions were based on static life table data, and so may not be an authentic representation of badger life histories.
In the following analyses, we use Bayesian approaches to demonstrate the benefits of using age-specific mortality trajectories to reveal when and how gender differences in disease-induced mortality occur. Research employing mortality trajectory analysis commonly recognises sex and individual heterogeneity [49,50] as critical drivers of mortality, although environmental variation has also been considered [23,51]. Males are often more susceptible to infection [52], have a weaker immune response [53], and greater disease-induced mortality [11,24], stemming from physiological variation [54] or behavioural and ecological differences that result in males being more likely to become infected [44]. Gender differences are not consistent across all species, and our understanding of the complex nature of sex influences on ageing patterns is far from complete [55]. Long-lived species are more likely to display negligible or negative patterns of senescence [56], but by using BaSTA to assess mortality trajectories, we uncover complex age- and sex-specific patterns of mortality that vary with infection status. We also reveal important limitations of standard mortality trajectory analysis for epidemiological research, including age-dependent predictors (disease state changes with age) and diagnostic uncertainty (tests for infection vary in their sensitivity and specificity towards the target pathogen [57]). We discuss alternative modelling strategies that will help to deal with these limitations.

2. Materials and Methods

2.1. Ecological Data

We used capture-mark-recapture (CMR) data from a population of wild badgers naturally infected with Mycobacterium bovis in Woodchester Park, Gloucestershire. Data for the present study consisted of badgers trapped from 1982 to 2015 inclusive. The badger population was sampled using live traps on (usually) four occasions per year. All trapped badgers were anaesthetized and subjected to several diagnostic tests for bTB before being released (for a more detailed account of the trapping and testing procedures, see [58,59]). On first capture, each badger was given a unique identifying tattoo, so that it could be identified at subsequent captures. Data available in Supplementary material—Hudson_2019_MTA_Data.

2.2. Diagnostic Tests

Samples of sputum, faeces, urine, and swabs of any abscesses or wounds were taken for M. bovis culture [60], as well as blood samples obtained for antibody tests as follows:
  • Brock ELISA (enzyme-linked immunosorbent assay) [61]—used from 1982 to 2006
  • BrockTB Stat-Pak lateral flow immunoassay (Chembio Diagnostics Systems, USA) [62]—used from 2006 to 2015
We used badgers of known age (i.e. badgers caught as cubs) and categorised them as “cub-positive” if they tested positive to either the Brock ELISA, Stat-Pak, or culture during the first year of their lives, and “never-positive” if they never tested positive to any test throughout their lives (sex was also included as a categorical covariate for each group). The distinction between these two groups allows for the comparison of disease effects on mortality trajectories across entire life histories, and identification of where variations occurred. The analysis of age-dependent acquisition of disease requires a more complex modelling framework, but for this analysis we filtered out any badgers that acquired test-positive status beyond the first year of life (see discussion for commentary on the risk of bias). We also recognise that the diagnostic tests employed in this study have limitations in terms of their sensitivity (Brock ELISA: 40.7% [63]; Stat-Pak: 49%; culture: 10% [64], yet all the tests used are highly specific (Brock ELISA: 94%–98% [61]; Stat-Pak: 97%; culture: 100% [64]). As a result, we are confident that individuals diagnosed as infected are very rarely truly uninfected, but can be rather less confident that individuals diagnosed as uninfected are truly uninfected. Our classification of “never-positive” removes some of this uncertainty, as an individual is required to have tested negative throughout their capture history and is not reliant on a single diagnostic test result (see discussion for implications).
Ongoing research seeks to clarify links between diagnostic outcomes and infection status [65,66], but for the purposes of our analyses here, we worked with diagnostic results rather than true infection status.
Badgers caught and identified as cubs were assigned a birth occasion as the first trapping season for the year of first capture, as the majority of cubs are born between February and March each year [67]. Badgers recovered dead were assigned a known death occasion as the time of the post-mortem examination.

2.3. Analysis

To understand age-specific mortality patterns in badgers, and how their mortality trajectories are affected by diagnosis of bTB infection, we used the package BaSTA [48] in R [68]. This package implements a hierarchical Bayesian model (for full details and likelihood functions, see [30]) and draws inference on age-specific survival from CMR data when large portions of the data consist of unknown birth and death years. Age-specific survival analysis requires the definition of mortality or hazard rate: we define a random variable X for ages at death, where X > 0 with any given age, represented by (lower case) x . Standard theory of survival models defines the hazard function ( μ ) (i.e., the instantaneous rate of death [69], given survival to age x ) as
μ ( x | b ) = lim Δ x 0 Pr ( x < X < x + Δ x | X x , b ) Δ x ,
where b is a vector of mortality parameters to be estimated, and Δ x is some very small time period. The cumulative hazard function ( H ) at age x , t years later is then defined as
H ( x | b ) = 0 x μ ( t | b )   d t ,
From Equations (1) and (2), it is possible to derive the survival function ( S ) as
S ( x | b ) = Pr { X > x | b } = exp   [ H ( x | b ) ] ,
and the probability density function (pdf) of ages at death as
f ( x | b ) = Pr   ( x X < ( x + Δ x ) ) = S ( x | b )   μ ( x | b )
BaSTA allows for the comparison of four different functional forms of the mortality function in Equation (1) (Exponential, Gompertz, Weibull, and Logistic), as well as the extension of the latter three to incorporate more complex shapes. The inclusion of a “Makeham” [31] term models the effect of age-independent mortality ( c ) . The further addition of a declining Gompertz function then allows the exploration of “bathtub” mortality shapes. In total, we compared ten different forms of mortality function (Table 1).
The Makeham structure consists of adding a constant ( c ) to the mortality ( μ b ) and survival rates ( S b ) when b is a vector of parameters to be estimated and x is age:
μ b ( x | b , c ) = c + μ b ( x | b ) ,
S b ( x | b , c ) = e c x S b ( x | b ) ,
with c > 0 if μ ( x ) is declining or c > μ ( 0 ) otherwise. The bathtub structures are constructed by adding a declining Gompertz function and a constant to the basic mortality:
μ b ( x | b , a , c ) = e a 0 a 1 x + c + μ b ( x | b ) ,
S 0 ( x | b , a ,   c ) = exp [ e a 0 a 1 ( e a 1 x 1 ) c x ] S b ( x | b ) ,
with < a 0   < , a 1 > 0 , and c > [ e a 0 a 1 x m i n + μ b ( x m i n | b ) ] , where x m i n is the age at which Equation (5) is at its lowest value.
Examples of the shapes of mortality trajectory that can be defined by these models are shown in Figure 1.
The models in Table 1 define parameters and describe mortality in four different ways: (1) exponential, in which the model assumes mortality to be constant and independent of age; (2) Gompertz [28], in which the model describes baseline mortality ( b 0 ) and an exponential increase with age ( b 1 ) ; (3) Weibull [29], in which b 0 is the shape parameter and b 1 is the scale parameter, the model assumes that mortality increases (or decreases) as a power function of age; and (4) logistic [33], in which the model consists of an initial exponential increase in mortality that decelerates to a plateau after a particular age. The b 2 parameter describes the degree of deceleration in mortality with age. Bathtub variants of the models include a declining Gompertz function e a 0 a 1 x , where a 0 represents initial mortality rate at birth and a 1 is the exponential decrease in age-dependent mortality from birth.
Badger sex and infection status (“cub-positive” or “never-positive”) were included as categorical covariates and incorporated into the models as linear functions of the survival parameters—for example, in a Gompertz mortality function:
μ ( x | θ ) = exp [ ( α T z ) b 0 + ( β T z ) b 1 x ] ,
where α T and β T are transposed vectors of linear coefficients that link the covariate z with the survival parameters b 0 and b 1 , such that θ = ( α ,   β ) . Covariates are included for all parameters in each model structure.
To ensure model convergence, initial trials of four Markov chain Monte Carlo (MCMC) chains were run for each model, followed by 1,000,000 iterations, with a burn-in of 10,001 iterations and thinning every 100 iterations [70]. Convergence was assessed visually, ensuring mixing of the chains, and formally within each model using the potential scale reduction ( R ^ ) [48]. Convergence is reasonable when R ^ 1 . We assessed the sensitivity of the mortality parameters to the choice of prior distributions [71] by running the analysis under four different prior structures, and found there to be no differences in the selection or identification of parameters as a result.
The fits of the ten models were compared using their deviance information criterion (DIC) [72], which is a measure of predictive power and criterion for model fit, akin to the Akaike information criterion and the Bayesian information criterion (for a review of its use, see [73]). To evaluate the impact of categorical covariates, BaSTA uses an adapted version of the Kullback–Liebler discrepancy [74,75], which estimates the degree of overlap in the posterior distributions of the parameter estimates. This indicator provides a value for k β [ 0.5 , 1 ] for a given parameter β (0.5 indicating full overlap, 1 meaning no overlap); k β 0.65 is generally interpreted as indicative of a difference between means that is unlikely to occur if the distributions of the two variables are the same [76].
Recapture probabilities were modelled as fully time-dependent, allowing the parameter estimate to vary for each occasion. Sample sizes for each type of datum are provided in Table 2. R code available in Supplementary material—RMarkdown_BaSTA_code.

3. Results

The Gompertz bathtub or Siler [40] function was the most supported mortality model across both the “cub-positive” (Figure 2a) and “never-positive” (Figure 2b) badger data sets, with substantial support as the “best” model and with a clear difference in DIC to the nearest rival model (Table 3). The Siler model is the sum of three different mortality models: the first describing a decrease in mortality over the initial phase of life, with e a 0 being the initial level and a 1 modelling the rate of decrease. The central “Makeham” term is a constant hazard, which is independent of age, and the final term is a Gompertz function, which describes mortality as increasing exponentially with a rate of b 1 from an initial level e b 0 .
Having revealed sex-, disease-, and age-specific patterns of mortality, we then analysed the entire data set using the Siler function, with disease status and sex as individual-level categorical covariates, creating four different groups (cub-positive male, cub-positive female, never-positive male, never-positive female). We assessed model fit by producing Kaplan–Meier plots of observed survival against predicted survival trajectories (Figure S1), and were satisfied that the Siler function and posterior parameter estimates were appropriate. We compared the posterior distributions between groups for each mortality parameter, both visually and using an adapted version of the Kullback–Liebler discrepancy [74] proposed by McCulloch [75] (Figure 3). The discrepancy measure works on the log of estimated probabilities, and due to the posterior densities of the a 1 parameters between “cub-positive” and “never-positive” female badgers being so different, the estimated probability produced zeroes, which in turn produced an NA in the Kullback–Liebler calculations.
Badgers diagnosed as infected shared very similar patterns of age-dependent mortality, with the only sex difference in mortality trajectories lying within the Makeham/age-independent parameter ( c ). “Never-positive” badgers displayed a greater degree of dissimilarity between the sexes. Male badgers had greater age-independent mortality ( c ) and initial “at birth” mortality rates ( a 0 ), which follows previous research [11,76]; however, we also uncovered sex variation in the senescent phase ( b 0 and b 1 ) of the mortality trajectory. “Never-positive” male badgers fail to display any evidence of a senescent decline in survival, whereas female badgers spend the majority of their lives with a comparatively lower mortality rate until the later phase of their lives, when they display a much clearer pattern of senescence. Infection raised the initial mortality rate ( a 0 ) and intensified the senescent decline ( b 1 ) in survival for both sexes, inducing a pattern of senescence in male badgers.
The parameter estimates also uncovered a lack of any real decline in mortality during the juvenile phase of the life history ( a 1 ) for “never-positive” individuals (male a 1 = 0.05 , female a 1 = 0.06 ), as well as extremely flat posteriors for “cub-positive” individuals, suggesting a lack of certainty. This pattern of mortality is further supported by the second “best” model, being the simple exponential (“never-positive”) and simple Gompertz (“cub-positive”), neither of which define an initial phase of reducing mortality for uninfected badgers.

4. Discussion

Mortality trajectory analysis has revealed complex, age-specific patterns of mortality, and provided evidence of actuarial senescence in badgers, a species previously believed to display near-constant rates of mortality throughout life. Our results discern the impact of bTB infection on wild badgers, and has allowed us to estimate where in a species’ life history disease- and sex-differences occur.
The inclusion of age as a continuous variable enabled us to generate mortality trajectories and clarify age-specific sex differences. In our case study, we found “never positive” male badgers to have a higher initial mortality rate (compared to females), which appears to decline with age, yet “never-positive” females show a more intense increase in later life (senescence). Higher rates of mortality in males resemble the sex bias evidenced across other mammalian species [77]. Explanations for higher rates of mortality among males include differences in intra-sexual reproductive competition [78], physiology [54], behaviour [79], and ecology [80]. In addition to an age-independent sex difference, the male bias in juvenile mortality hints at differences in physiology being the key driver: although behavioural differences could appear, ecological explanations are unlikely to have come into effect at this early stage of a badger’s life. An unanticipated result was the more intense senescent increase in mortality with age found in female badgers. Males do not live as long as females (11% reduction in lifespan irrespective of infection status), and it is possible that this reduced lifespan removes or reduces the potential to exhibit detectable signs of senescence. There is a growing body of evidence reporting sex differences in senescence (e.g., [81,82]), with the majority indicating a greater effect in males. A number of reasons have been proposed, mostly associated with a greater impact of sexual selection in males, although this has been challenged in a study of herbivores [83]. The pattern of “never-positive” male mortality found here is more in line with previous findings that show negligible actuarial senescence in badgers [44].
It was already known that bTB infection, as indicated from several diagnostic tests, results in increased male mortality in badgers [11,24], but our analysis provides evidence for a similar effect in females. Mortality trajectories of the “cub-positive” individuals revealed novel patterns of senescence and indicated that infection promotes a senescent increase in late-life mortality; this effect is more pronounced in male badgers, thus homogenizing the shape of the mortality trajectories and removing any signature of age-dependent sex differences. Recent work has revealed a diversity of age-specific patterns of mortality, both within and amongst species [35], with sex and individual differences in frailty being highlighted as the most important sources of variation in survival parameters [84]. Our analyses found only limited support for individual heterogeneity amongst badgers, indicated by the low ranking of the logistic model (third for infected badgers and fifth for uninfected). Males and females, on the other hand, do display differences in individual heterogeneity, and may in fact be better described by sex-specific models [27]. Our combined sex analysis may have masked such variation, suggesting that any reduction in mortality with advancing age could be representative of real change [34], or be an artefact of heterogeneity among individuals in the population [33]. Previous studies have suggested that ignoring frailty may result in biased parameter estimates, but with only moderate support for the logistic mortality model it is difficult to draw any firm conclusions. Should individual heterogeneity in frailty exist, then the more frail individuals with higher mortality should be removed from the population as a result of within-cohort selection [85]. The lack of evidence for individual heterogeneity in the present study may in part be due to aspects of badger ecology. Most badgers give birth between mid-January and mid-March [86], then cubs spend their first few weeks underground and are unavailable for capture. Previous research has estimated pre-capture mortality at 24% [67], and perhaps it is at this stage when the majority of weaker individuals are lost from the population. Pre-capture mortality has not been included in our analysis, and with trapping not taking place during March and April, cubs may be upwards of 6 months old when they are first captured and tested for bTB. Consequently, our estimates of initial mortality are likely to be underestimates, and in reality, the downward curve assumed in the juvenile phase of the Siler function would be more pronounced than presented here.
We used the package BaSTA, with age-specific survival data and a hierarchical structure within a Bayesian framework, to draw inference from censored CMR data, in order to investigate survival and mortality in a badger population naturally infected with bTB. Our study has demonstrated the benefits of this approach to mortality trajectory analysis for the study of wildlife disease epidemiology, and clarified age-specific sex and disease variation in patterns of mortality. There are, however, inherent weaknesses in our modelling framework, and to overcome them will require a more flexible modelling approach. BaSTA dictated our categorisation of badgers into two fixed groups of “cub-positive” and “never positive” individuals, as it is unable to assess the influence of time-varying disease states on mortality parameters. This restriction meant filtering out badgers that first tested positive after their first year of life, and leaves unanswered questions regarding the mortality trajectories of such individuals. The inability to model time-varying disease states also forced us to use a binary disease classification. Recent research into the epidemiology of bTB in badgers has identified at least four different states of infection (susceptible, test positive, single site excretor, multi-site/occasion excretor), with differing mortality rates at each stage [11], and high levels of individual heterogeneity in disease progression over time [24]. This adds further complexity, and suggests the potential for differing mortality trajectories at each disease state, as well as for each age-at-infection/state—if, for example, older badgers are more likely to become infected and experience more rapid disease progression. To model this system more completely will require a more flexible multi-state framework that allows for time-varying transition through multiple disease states. Further complications relevant to studies of bTB are the limitations of diagnostic tests and changes in the prevailing tests employed in studies over time [63,87]. False-positive diagnoses would potentially weaken the signal of mortality in “cub-positive” badgers, whereas the more likely false negatives could strengthen the signal in the “never-positive” individuals. Our results may therefore be more conservative for “cub-positive”, yet potentially overstated for “never-positive” individuals. We are confident that any impact is minimal, due to tests being highly specific (low chance of false-positive) and recapture being common, with an average of 5.9 capture/testing occasions per badger (therefore reducing the chance of any individual being wrongly included in either category). There was no bias in the recapture rate between categories, which reduces any impact of selective disappearance. The issue of test performance could be explicitly addressed by the further development of mixture models that can incorporate diagnostic uncertainty generated from false-positive and false-negative test results [57]. Multi-event methods have been developed to deal with such problems [17,88,89], often implementing the software package E-Surge [90]. To fully address all the difficulties outlined here requires the combination of these methods and a need to develop more complex, hierarchical state-space models developed from first principles, using MCMC software such as BUGS, JAGS, or Stan, which also allow for the inclusion of age as a continuous covariate.
We have demonstrated the use of relatively simple mortality trajectory analysis in the study of wildlife disease epidemiology. We recommend widespread use of this approach in systems that have sufficient information to infer patterns of mortality across the whole life cycle. Despite senescence being notoriously difficult to detect in wild populations [25], the present study has provided credible evidence of actuarial senescence in badgers. Furthermore, it has revealed that senescence can be intensified in diseased individuals. Understanding the physiological, evolutionary, and ecological drivers of these mortality trajectories remains an ongoing challenge [23]. BaSTA is a powerful tool to begin exploration whilst accommodating many of the problems associated with CMR data. However, there are limitations associated with age-dependence and diagnostic uncertainty of disease states, and we are now developing similar Bayesian methods with more complex hierarchical and multi-state frameworks.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1424-2818/11/10/182/s1, Figure S1: Kaplan–Meier plots (in red) of observed survival on top of predicted survival trajectories and 95% confidence intervals for a population of European badgers naturally infected with bovine tuberculosis, split by sex and infection status: (a) “cub-positive” females; (b) “cub-positive” males; (c) “never-positive” females; (d) “never-positive” males. Table S1: Posterior parameter estimates, standard deviations, and lower/upper 95% confidence intervals of mortality trajectories for a population of wild European badgers naturally infected with bovine tuberculosis. R Code: RMarkdown_BaSTA_code. Data: Hudson_2019_MTA_Data.

Author Contributions

Conceptualization, D.W.H., D.J.H., and R.A.M.; methodology, D.W.H., R.D., D.J.H., and T.J.M.; writing—original draft preparation, D.W.H.; writing—review and editing, D.W.H., D.J.H., R.D., R.A.M. and T.J.M.; supervision, D.J.H. and T.J.M.; funding acquisition, D.J.H., R.A.M., T.J.M., and R.D.

Funding

This research was funded by the Natural Environment Research Council (Research council grant ref: NE/M010260/1), the Animal and Plant Health Agency, and the College of Life and Environmental Sciences, University of Exeter. The Woodchester Park badger study is funded by the Department for Environment, Food and Rural Affairs.

Acknowledgments

We are grateful to the Animal and Plant Health Agency’s Woodchester Park team for data collection and managing the long-term dataset and DEFRA for funding fieldwork activities at Woodchester Park. We also wish to thank the many land owners and occupiers for consent to work on their land. Finally we would like to thank the anonymous reviewers for their comments on the manuscript.

Ethics Statement

The work was carried out under Home Office license and approved by the APHA Ethics Panel.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; and in the decision to publish the results.

References

  1. Delahay, R.J.; Smith, G.C.; Hutchings, M.R. The science of wildlife disease management. In Management of Disease in Wild Mammals; Springer: Tokyo, Japan, 2009; pp. 1–8. [Google Scholar]
  2. Gortázar, C.; Ferroglio, E.; Höfle, U.; Frölich, K.; Vicente, J. Diseases shared between wildlife and livestock: A European perspective. Eur. J. Wildl. Res. 2007, 53, 241–256. [Google Scholar] [CrossRef]
  3. Wiethoelter, A.K.; Beltrán-Alcrudo, D.; Kock, R.; Mor, S.M. Global trends in infectious diseases at the wildlife-livestock interface. Proc. Natl. Acad. Sci. USA 2015, 112, 9662–9667. [Google Scholar] [CrossRef] [PubMed]
  4. Heisey, D.M.; Joly, D.O.; Messier, F. The fitting of general force-of-infection models to wildlife disease prevalence data. Ecology 2006, 87, 2356–2365. [Google Scholar] [CrossRef]
  5. Samuel, M.D.; Woodworth, B.L.; Atkinson, C.T.; Hart, P.J.; Lapointe, D.A. Avian malaria in Hawaiian forest birds: Infection and population impacts across species and elevations. Ecosphere 2015, 6, 1–21. [Google Scholar] [CrossRef]
  6. Kilpatrick, A.M.; Briggs, C.J.; Daszak, P. The ecology and impact of chytridiomycosis: An emerging disease of amphibians. Trends Ecol. Evol. 2010, 25, 109–118. [Google Scholar] [CrossRef] [PubMed]
  7. Frick, W.F.; Pollock, J.F.; Hicks, A.C.; Langwig, K.E.; Reynolds, D.S.; Turner, G.G.; Butchkoski, C.M.; Kunz, T.H. An emerging disease causes regional population collapse of a common North American bat species. Science 2010, 329, 679–682. [Google Scholar] [CrossRef]
  8. van Riper, C.; van Riper, S.G.; Goff, M.L.; Laird, M. The epizootiology and ecological significance of malaria in Hawaiian land birds. Ecol. Monogr. 1986, 56, 327–344. [Google Scholar] [CrossRef]
  9. Berger, L.; Roberts, A.A.; Voyles, J.; Longcore, J.E.; Murray, K.A.; Skerratt, L.F. History and recent progress on chytridiomycosis in amphibians. Fungal Ecol. 2016, 19, 89–99. [Google Scholar] [CrossRef]
  10. Samuel, M.D.; Storm, D.J. Chronic wasting disease in white-tailed deer: Infection, mortality, and implications for heterogeneous transmission. Ecology 2016, 97, 3195–3205. [Google Scholar] [CrossRef]
  11. Wilkinson, D.; Smith, G.C.; Delahay, R.J.; Rogers, L.M.; Cheeseman, C.L.; Clifton-Hadley, R.S. The effects of bovine tuberculosis (Mycobacterium bovis) on mortality in a badger (Meles meles) population in England. J. Zool. 2000, 250, 389–395. [Google Scholar] [CrossRef]
  12. Fontana, L.; Kennedy, B.K.; Longo, V.D.; Seals, D.; Melov, S. Medical research: Treat ageing. Nature 2014, 511, 405–407. [Google Scholar] [CrossRef] [PubMed]
  13. Hayflick, L. The future of ageing. Nature 2000, 408, 267–269. [Google Scholar] [CrossRef] [PubMed]
  14. Haas, G.P.; Delongchamps, N.; Brawley, O.W.; Wang, C.Y.; de la Roza, G. The worldwide epidemiology of prostate cancer: Perspectives from autopsy studies. Can. J. Urol. 2008, 15, 3866–3871. [Google Scholar] [PubMed]
  15. Jorgenson, J.T.; Festa-Bianchet, M.; Gaillard, J.-M.; Wishart, W.D. Effects of age, sex, disease, and density on survival of bighorn sheep. Ecology 1997, 78, 1019–1032. [Google Scholar] [CrossRef]
  16. Larsen, D.G.; Gauthier, D.A.; Markel, R.L. Causes and rate of moose mortality in the Southwest Yukon. J. Wildl. Manag. 1989, 53, 548–557. [Google Scholar] [CrossRef]
  17. Koons, D.N.; Gamelon, M.; Gaillard, J.-M.; Aubry, L.M.; Rockwell, R.F.; Klein, F.; Choquet, R.; Gimenez, O. Methods for studying cause-specific senescence in the wild. Methods Ecol. Evol. 2014, 5, 924–933. [Google Scholar] [CrossRef] [Green Version]
  18. Reid, J.M.; Bignal, E.M.; Bignal, S.; McCracken, D.I.; Monaghan, P. Age-specific reproductive performance in red-billed choughs Pyrrhocorax pyrrhocorax: Patterns and processes in a natural population. J. Anim. Ecol. 2003, 72, 765–776. [Google Scholar] [CrossRef]
  19. Reimers, E.; Holmengen, N.; Mysterud, A. Life-history variation of wild reindeer (Rangifer tarandus) in the highly productive North Ottadalen region, Norway. J. Zool. 2005, 265, 53–62. [Google Scholar] [CrossRef]
  20. Chen, G.H.; Wang, Y.J.; Wang, X.M.; Zhou, J.N.; Liu, R.Y. Effect of aging on species-typical behaviors in senescence-accelerated mouse. Physiol. Behav. 2005, 85, 536–545. [Google Scholar] [CrossRef]
  21. Angelier, F.; Weimerskirch, H.; Dano, S.; Chastel, O. Age, experience and reproductive performance in a long-lived bird: A hormonal perspective. Behav. Ecol. Sociobiol. 2007, 61, 611–621. [Google Scholar] [CrossRef]
  22. Ricklefs, R.E. Evolutionary theories of aging: Confirmation of a fundamental prediction, with implications for the genetic basis and evolution of life span. Am. Nat. 1998, 152, 24–44. [Google Scholar] [CrossRef] [PubMed]
  23. Nussey, D.H.; Froy, H.; Lemaitre, J.-F.; Gaillard, J.-M.; Austad, S.N. Senescence in natural populations of animals: Widespread evidence and its implications for bio-gerontology. Ageing Res. Rev. 2013, 12, 214–225. [Google Scholar] [CrossRef] [PubMed]
  24. Graham, J.; Smith, G.C.; Delahay, R.J.; Bailey, T.; McDonald, R.A.; Hodgson, D. Multi-state modelling reveals sex-dependent transmission, progression and severity of tuberculosis in wild badgers. Epidemiol. Infect. 2013, 141, 1429–1436. [Google Scholar] [CrossRef] [PubMed]
  25. Nussey, D.H.; Coulson, T.; Festa-Bianchet, M.; Gaillard, J.-M. Measuring senescence in wild animal populations: Towards a longitudinal approach. Funct. Ecol. 2008, 22, 393–406. [Google Scholar] [CrossRef]
  26. Ricklefs, R.E.; Scheuerlein, A. Biological implications of the Weibull and Gompertz models of aging. J. Gerontol. Ser. A Biol. Sci. Med. Sci. 2002, 57, B69–B76. [Google Scholar] [CrossRef] [PubMed]
  27. Wilson, D.L. The analysis of survival (mortality) data: Fitting Gompertz, Weibull, and logistic functions. Mech. Ageing Dev. 1994, 74, 15–33. [Google Scholar] [CrossRef]
  28. Gompertz, B. On the nature of the function expressive of the law of human mortality, and on a new. Philos. Trans. R. Soc. Lond. 1825, 115, 513–583. [Google Scholar]
  29. Pinder, J.E.; Wiener, J.G.; Smith, M.H. The Weibull distribution: A new method of summarizing survivorship data. Ecology 1978, 59, 175–179. [Google Scholar] [CrossRef]
  30. Colchero, F.; Clark, J.S. Bayesian inference on age-specific survival for censored and truncated data. J. Anim. Ecol. 2012, 81, 139–149. [Google Scholar] [CrossRef] [PubMed]
  31. Makeham, W.M. On the law of mortality. J. Inst. Actuar. 1867, 13, 325–358. [Google Scholar] [CrossRef]
  32. Cam, E.; Aubry, L.M.; Authier, M. The conundrum of heterogeneities in life history studies. Trends Ecol. Evol. 2016, 31, 872–886. [Google Scholar] [CrossRef] [PubMed]
  33. Vaupel, J.W.; Manton, K.G.; Stallard, E. The impact of heterogeneity in individual frailty on the dynamics of mortality. Demography 1979, 16, 439–454. [Google Scholar] [CrossRef] [PubMed]
  34. Pletcher, S.D. Pletcher model fitting and hypothesis testing for age-specific mortality data. J. Evol. Biol. 1999, 12, 430–439. [Google Scholar] [CrossRef]
  35. Jones, O.R.; Scheuerlein, A.; Salguero-Gómez, R.; Camarda, C.G.; Schaible, R.; Casper, B.B.; Dahlgren, J.P.; Ehrlén, J.; García, M.B.; Menges, E.S.; et al. Diversity of ageing across the tree of life. Nature 2014, 505, 169–173. [Google Scholar] [CrossRef] [PubMed]
  36. Vaupel, J.W.; Carey, J.R.; Christensen, K.; Johnson, T.E.; Yashin, A.I.; Holm, N.V.; Iachine, I.A.; Kannisto, V.; Khazaeli, A.A.; Liedo, P.; et al. Biodemographic trajectories of longevity. Science 1998, 280, 855–860. [Google Scholar] [CrossRef] [PubMed]
  37. Ricklefs, R.E.; Scheuerlein, A. Comparison of aging-related mortality among birds and mammals. Exp. Gerontol. 2001, 36, 845–857. [Google Scholar] [CrossRef]
  38. Barthold, J.A.; Loveridge, A.J.; Macdonald, D.W.; Packer, C.; Colchero, F. Bayesian estimates of male and female African lion mortality for future use in population management. J. Appl. Ecol. 2016, 53, 295–304. [Google Scholar] [CrossRef] [Green Version]
  39. Caughley, G. Mortality patterns in mammals. Ecology 1966, 47, 906–918. [Google Scholar] [CrossRef]
  40. Siler, W. A competing-risk model for animal mortality. Ecology 1979, 60, 750–757. [Google Scholar] [CrossRef]
  41. Klutke, G.A.; Kiessler, P.C.; Wortman, M.A. A critical look at the bathtub curve. IEEE Trans. Reliab. 2003, 52, 125–129. [Google Scholar] [CrossRef]
  42. Bebbington, M.; Lai, C.-D.; Zitikis, R. Modeling human mortality using mixtures of bathtub shaped failure distributions. J. Theor. Biol. 2007, 245, 528–538. [Google Scholar] [CrossRef] [PubMed]
  43. Snoke, M.S.; Promislow, D.E.L. Quantitative genetic tests of recent senescence theory: Age-specific mortality and male fertility in Drosophila melanogaster. Heredity 2003, 91, 546–556. [Google Scholar] [CrossRef]
  44. McDonald, J.L.; Smith, G.C.; McDonald, R.A.; Delahay, R.J.; Hodgson, D. Mortality trajectory analysis reveals the drivers of sex-specific epidemiology in natural wildlife-disease interactions. Proc. R. Soc. B Biol. Sci. 2014, 281, 20140526. [Google Scholar] [CrossRef]
  45. Delahay, R.J.; Walker, N.; Smith, G.S.; Wilkinson, D.; Clifton-Hadley, R.S.; Cheeseman, C.L.; Tomlinson, A.J.; Chambers, M.A. Long-term temporal trends and estimated transmission rates for Mycobacterium bovis infection in an undisturbed high-density badger (Meles meles) population. Epidemiol. Infect. 2013, 141, 1445–1456. [Google Scholar] [CrossRef] [PubMed]
  46. Donnelly, C.A.; Woodroffe, R.; Cox, D.R.; Bourne, F.J.; Cheeseman, C.L.; Clifton-Hadley, R.S.; Wei, G.; Gettinby, G.; Gilks, P.; Jenkins, H.; et al. Positive and negative effects of widespread badger culling on tuberculosis in cattle. Nature 2006, 439, 843–846. [Google Scholar] [CrossRef]
  47. Gallagher, J.; Clifton-Hadley, R.S. Tuberculosis in badgers; a review of the disease and its significance for other animals. Res. Vet. Sci. 2000, 69, 203–217. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Colchero, F.; Jones, O.R.; Rebke, M. BaSTA: An R package for Bayesian estimation of age-specific survival from incomplete mark-recapture/recovery data with covariates. Methods Ecol. Evol. 2012, 3, 466–470. [Google Scholar] [CrossRef]
  49. Colchero, F.; Aliaga, A.E.; Jones, O.R.; Conde, D.A. Individual heterogeneity determines sex differences in mortality in a monogamous bird with reversed sexual dimorphism. J. Anim. Ecol. 2017, 86, 899–970. [Google Scholar] [CrossRef] [PubMed]
  50. Tidière, M.; Gaillard, J.-M.; Müller, D.W.H.; Lackey, L.B.; Gimenez, O.; Clauss, M.; Lemaître, J.-F. Does sexual selection shape sex differences in longevity and senescence patterns across vertebrates? A review and new insights from captive ruminants. Evolution 2015, 69, 3123–3140. [Google Scholar] [CrossRef]
  51. Lemaître, J.-F.; Gaillard, J.-M.; Lackey, L.B.; Clauss, M.; Müller, D.W.H. Comparing free-ranging and captive populations reveals intra-specific variation in aging rates in large herbivores. Exp. Gerontol. 2013, 48, 162–167. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Klein, S.L. The effects of hormones on sex differences in infection: From genes to behavior. Neurosci. Biobehav. Rev. 2000, 24, 627–638. [Google Scholar] [CrossRef]
  53. Møller, A.P. Advantages and disadvantages of coloniality in the swallow, Hirundo rustica. Anim. Behav. 1987, 35, 819–832. [Google Scholar] [CrossRef]
  54. Zuk, M.; McKean, K.A. Sex differences in parasite infections: Patterns and processes. Int. J. Parasitol. 1996, 26, 1009–1024. [Google Scholar] [CrossRef]
  55. Berger, V.; Lemaître, J.-F.; Dupont, P.; Allainé, D.; Gaillard, J.-M.; Cohas, A. Age-specific survival in the socially monogamous alpine marmot (Marmota marmota): Evidence of senescence. J. Mammal. 2016, 97, 992–1000. [Google Scholar] [CrossRef]
  56. Baudisch, A. The pace and shape of ageing. Methods Ecol. Evol. 2011, 2, 375–382. [Google Scholar] [CrossRef]
  57. McDonald, J.L.; Hodgson, D.J. Prior precision, prior accuracy, and the estimation of disease prevalence using imperfect diagnostic tests. Front. Vet. Sci. 2018, 5, 83. [Google Scholar] [CrossRef]
  58. McDonald, J.L.; Robertson, A.; Silk, M.J. Wildlife disease ecology from the individual to the population: Insights from a long-term study of a naturally infected European badger population. J. Anim. Ecol. 2018, 87, 101–112. [Google Scholar] [CrossRef] [PubMed]
  59. Delahay, R.J.; Langton, S.; Smith, G.C.; Clifton-Hadley, R.S.; Cheeseman, C.L. The spatio-temporal distribution of Mycobacterium bovis (bovine tuberculosis) infection in a high-density badger population. J. Anim. Ecol. 2000, 69, 428–441. [Google Scholar] [CrossRef]
  60. Gallagher, J.; Horwill, D.M. A selective oleic acid albumin agar medium for the cultivation of Mycobacterium bovis. Epidemiol. Infect. 1977, 79, 155–160. [Google Scholar] [CrossRef]
  61. Goodger, J.; Nolan, A.; Russell, W.P.; Dalley, D.J.; Thorns, C.J.; Stuart, F.A.; Croston, P.; Newell, D.G. Serodiagnosis of Mycobacterium bovis infection in badgers: Development of an indirect ELISA using a 25 kDa antigen. Vet. Rec. 1994, 135, 82–85. [Google Scholar] [CrossRef]
  62. Chambers, M.A.; Crawshaw, T.; Waterhouse, S.; Delahay, R.; Hewinson, R.G.; Lyashchenko, K.P. Validation of the BrockTB stat-pak assay for detection of tuberculosis in Eurasian badgers (Meles meles) and influence of disease severity on diagnostic accuracy. J. Clin. Microbiol. 2008, 46, 1498–1500. [Google Scholar] [CrossRef] [PubMed]
  63. Clifton-Hadley, R.S.; Sayers, A.R.; Stock, M.P. Evaluation of an ELISA for Mycobacterium bovis infection in badgers (Meles meles). Vet. Rec. 1995, 137, 555–558. [Google Scholar] [CrossRef] [PubMed]
  64. Drewe, J.A.; Tomlinson, A.J.; Walker, N.J.; Delahay, R.J. Diagnostic accuracy and optimal use of three tests for tuberculosis in live badgers. PLoS ONE 2010, 5, e11196. [Google Scholar] [CrossRef] [PubMed]
  65. Buzdugan, S.N.; Chambers, M.A.; Delahay, R.J.; Drewe, J.A. Diagnosis of tuberculosis in groups of badgers: An exploration of the impact of trapping efficiency, infection prevalence and the use of multiple tests. Epidemiol. Infect. 2016, 144, 1717–1727. [Google Scholar] [CrossRef]
  66. Wawegama, N.K.; Markham, P.F.; Kanci, A.; Schibrowski, M.; Oswin, S.; Barnes, T.S.; Firestone, S.M.; Mahony, T.J.; Browning, G.F. Evaluation of an IgG enzyme-linked immunosorbent assay as a serological assay for detection of mycoplasma bovis infection in feedlot cattle. J. Clin. Microbiol. 2016, 54, 1269–1275. [Google Scholar] [CrossRef]
  67. Rogers, L.M.; Cheeseman, C.L.; Mallinson, P.J.; Clifton-Hadley, R. The demography of a high-density badger (Meles meles) population in the west of England. J. Zool. 1997, 242, 705–728. [Google Scholar] [CrossRef]
  68. R Development Core Team. R: A language and environment for statistical computing. R Found. R Found. Stat. Comput. 2011. Available online: https://www.r-project.org/ (accessed on 12 February 2019).
  69. Kalbfleisch, J.D.; Prentice, R.L. The Statistical Analysis of Failure Time Data; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2002; ISBN 9781118032985. [Google Scholar]
  70. Gelman, A.; Carlin, J.B.B.; Stern, H.S.S.; Dunson, D.B.; Vehtari, A.; Rubin, D.B.B. Bayesian Data Analysis, 3rd ed.; Texts in Statistical Science; CRC Press: Boca Raton, FL, USA, 2014; ISBN 9781439840955. [Google Scholar]
  71. Kass, R.E.; Raftery, A.E. Bayes factors. J. Am. Stat. Assoc. 1995, 90, 773–795. [Google Scholar] [CrossRef]
  72. Spiegelhalter, D.J.; Best, N.G.; Carlin, B.P.; van der Linde, A. Bayesian measures of model complexity and fit. J. R. Stat. Soc. Ser. B 2002, 64, 583–639. [Google Scholar] [CrossRef] [Green Version]
  73. Spiegelhalter, D.J.; Best, N.G.; Carlin, B.P.; van der Linde, A. The deviance information criterion: 12 years on. J. R. Stat. Soc. Ser. B 2014, 76, 485–493. [Google Scholar] [CrossRef]
  74. Kullback, S.; Leibler, R.A. On information and sufficiency. Ann. Math. 1951, 22, 79–86. [Google Scholar] [CrossRef]
  75. McCulloch, R.E. Local model influence. J. Am. Stat. Assoc. 1989, 84, 473–478. [Google Scholar] [CrossRef]
  76. Larson, S.M.; Colchero, F.; Jones, O.R.; Williams, L.; Fernandez-Duque, E. Age and sex-specific mortality of wild and captive populations of a monogamous pair-bonded primate (Aotus azarae). Am. J. Primatol. 2016, 78, 315–325. [Google Scholar] [CrossRef] [PubMed]
  77. Bronikowski, A.M.; Altmann, J.; Brockman, D.K.; Cords, M.; Fedigan, L.M.; Pusey, A.; Stoinski, T.; Morris, W.F.; Strier, K.B.; Alberts, S.C. Aging in the natural world: Comparative data reveal similar mortality patterns across primates. Science 2011, 311, 1325–1328. [Google Scholar] [CrossRef] [PubMed]
  78. Beirne, C.; Delahay, R.; Young, A. Sex differences in senescence: The role of intra-sexual competition in early adulthood. Proc. R. Soc. B Biol. Sci. 2015, 282, 20151086. [Google Scholar] [CrossRef] [PubMed]
  79. Byrnes, J.P.; Miller, D.C.; Schafer, W.D. Gender differences in risk taking: A meta-analysis. Psychol. Bull. 1999, 125, 367. [Google Scholar] [CrossRef]
  80. Delahay, R.J.; Walker, N.J.; Forrester, G.J.; Harmsen, B.; Riordan, P.; Macdonald, D.W.; Newman, C.; Cheeseman, C.L. Demographic correlates of bite wounding in Eurasian badgers, Meles meles L., in stable and perturbed populations. Anim. Behav. 2006, 71, 1047–1055. [Google Scholar] [CrossRef]
  81. Greiner, S.; Nagy, M.; Mayer, F.; Knörnschild, M.; Hofer, H.; Voigt, C.C. Sex-biased senescence in a polygynous bat species. Ethology 2014, 120, 197–205. [Google Scholar] [CrossRef]
  82. Descamps, S.; Boutin, S.; Berteaux, D.; Gaillard, J.M. Age-specific variation in survival, reproductive success and offspring quality in red squirrels: Evidence of senescence. Oikos 2008, 117, 1406–1416. [Google Scholar] [CrossRef]
  83. Lemaitre, J.F.; Gaillard, J.M. Male survival patterns do not depend on male allocation to sexual competition in large herbivores. Behav. Ecol. 2012, 24, 421–428. [Google Scholar] [CrossRef] [Green Version]
  84. Noonburg, E.G.; Chen, A.; Shima, J.S.; Swearer, S.E. Demographic heterogeneity and the dynamics of open populations. Ecology 2015, 96, 1159–1165. [Google Scholar] [CrossRef] [PubMed]
  85. Dahlgren, J.P.; Colchero, F.; Jones, O.R.; Øien, D.-I.; Moen, A.; Sletvold, N. Actuarial senescence in a long-lived orchid challenges our current understanding of ageing. Proc. R. Soc. B Biol. Sci. 2016, 283, 20161217. [Google Scholar] [CrossRef] [PubMed]
  86. Kynaston, S.; Neal, E.; Cheeseman, C. Badgers. J. Anim. Ecol. 2006, 65, 533. [Google Scholar] [CrossRef]
  87. Dalley, D.; Davé, D.; Lesellier, S.; Palmer, S.; Crawshaw, T.; Hewinson, R.G.; Chambers, M. Development and evaluation of a gamma-interferon assay for tuberculosis in badgers (Meles meles). Tuberculosis 2008, 88, 235–243. [Google Scholar] [CrossRef]
  88. Conn, P.B.; Cooch, E.G. Multistate capture-recapture analysis under imperfect state observation: An application to disease models. J. Appl. Ecol. 2009, 46, 486–492. [Google Scholar] [CrossRef]
  89. Pradel, R.; Hines, J.E.; Lebreton, J.-D.; Nichols, J.D. Capture-recapture survival models taking account of transients. Biometrics 1997, 53, 60–72. [Google Scholar] [CrossRef]
  90. Choquet, R.; Rouan, L.; Pradel, R. Program E-surge: A software application for fitting multievent models. In Modeling Demographic Processes in Marked Populations; Springer: Boston, MA, USA, 2009; pp. 845–865. [Google Scholar]
Figure 1. Possible mortality trajectories, μ(x|θ) when θ is a vector of mortality parameters to be estimated, resulting from the four models included in BaSTA. Mortality or hazard rate is shown on the y-axis, with age as a continuous variable on the x-axis. The line styles represent examples of each model shape that are possible to test: simple, Makeham, and bathtub. The exponential model is only applicable in a simple format.
Figure 1. Possible mortality trajectories, μ(x|θ) when θ is a vector of mortality parameters to be estimated, resulting from the four models included in BaSTA. Mortality or hazard rate is shown on the y-axis, with age as a continuous variable on the x-axis. The line styles represent examples of each model shape that are possible to test: simple, Makeham, and bathtub. The exponential model is only applicable in a simple format.
Diversity 11 00182 g001
Figure 2. Age-specific survival and mortality trajectories of European badgers from a wild population naturally infected with bovine tuberculosis: (a) “cub-positive” individuals (badgers that tested positive in the first year of their lives) and (b) “never-infected” individuals (badgers that never tested positive throughout their lives). Female (red) and male (blue) estimated survival and mortality curves of the Siler function. Coloured areas surrounding the curve represent 95% credible intervals. Age in years.
Figure 2. Age-specific survival and mortality trajectories of European badgers from a wild population naturally infected with bovine tuberculosis: (a) “cub-positive” individuals (badgers that tested positive in the first year of their lives) and (b) “never-infected” individuals (badgers that never tested positive throughout their lives). Female (red) and male (blue) estimated survival and mortality curves of the Siler function. Coloured areas surrounding the curve represent 95% credible intervals. Age in years.
Diversity 11 00182 g002
Figure 3. Posterior distributions of mortality trajectory parameter estimates for a population of wild European badgers naturally infected with bovine tuberculosis. Female (red) and male (blue) posterior distributions of mortality parameters of the Siler function, split by infection status. “Cub-positive” refers to individuals who tested positive in the first year of their life. “Never-positive” refers to individuals who never tested positive throughout their lives. Adapted Kullback–Liebler discrepancy measures are shown in black, indicating the degree of distribution overlap between classifications of sex and disease status. Measures over 0.65 are considered important (0.5 = distributions are identical, 1 = no overlap). Iterations = 1,000,000; burn in = 10,001; thinning = 200; number of chains = 4. (For table of estimates and 95% confidence intervals, see Table S1).
Figure 3. Posterior distributions of mortality trajectory parameter estimates for a population of wild European badgers naturally infected with bovine tuberculosis. Female (red) and male (blue) posterior distributions of mortality parameters of the Siler function, split by infection status. “Cub-positive” refers to individuals who tested positive in the first year of their life. “Never-positive” refers to individuals who never tested positive throughout their lives. Adapted Kullback–Liebler discrepancy measures are shown in black, indicating the degree of distribution overlap between classifications of sex and disease status. Measures over 0.65 are considered important (0.5 = distributions are identical, 1 = no overlap). Iterations = 1,000,000; burn in = 10,001; thinning = 200; number of chains = 4. (For table of estimates and 95% confidence intervals, see Table S1).
Diversity 11 00182 g003
Table 1. Basic mortality and survival probability functions available to test in Bayesian survival trajectory analysis (BaSTA)—namely, exponential, Gompertz, Weibull, and logistic. Additional terms can be added to each model (except exponential) to test different shapes, namely “Makeham” and “bathtub”. Parameter constraints are also indicated.
Table 1. Basic mortality and survival probability functions available to test in Bayesian survival trajectory analysis (BaSTA)—namely, exponential, Gompertz, Weibull, and logistic. Additional terms can be added to each model (except exponential) to test different shapes, namely “Makeham” and “bathtub”. Parameter constraints are also indicated.
Model Mortality   Rate   μ b ( x | b ) Survival   Probability   S b ( x | b ) Parameters
Exponential b 0 e b x b 0 > 0
Gompertz [28] exp ( b 0 + b 1 x ) exp [ e b 0 b 1   ( 1 e b 1 x ) ] < b 0 ,   b 1 <
Weibull [29] b 0 b 1 ( b 1 x ) b 0 1 exp [ ( b 1 x ) b 0 ] b 0 , b 1 > 0
Logistic [33] exp ( b 0 + b 1 x ) 1 +   ( e b 0 b 1 )   b 2 ( e b 1 x 1 ) ( 1 + b 2 e b 0 b 1 ( e b 1 x 1 ) ) 1 b 2 b 0 , b 1 ,   b 2 > 0
Table 2. Summary of capture-mark-recapture data sets describing badgers in one of two health states: “cub-positive” (badgers that tested positive in the first year of their lives) and “never-positive” (badgers that have never tested positive throughout their lives). Known birth years refers to cubs captured and identified as being within the first year of their lives. Known death years refers to badgers who were recovered dead and subjected to post-mortem. Detections are the total number of capture events that took place over the duration of the study.
Table 2. Summary of capture-mark-recapture data sets describing badgers in one of two health states: “cub-positive” (badgers that tested positive in the first year of their lives) and “never-positive” (badgers that have never tested positive throughout their lives). Known birth years refers to cubs captured and identified as being within the first year of their lives. Known death years refers to badgers who were recovered dead and subjected to post-mortem. Detections are the total number of capture events that took place over the duration of the study.
Summary StatisticCub-PositiveNever-Positive
Total number badgers428 (M 191; F 237)1768 (M 833; F 935)
Number of known birth years4281768
Number of known death years13323
Total number of detections25157588
Table 3. Ranked list of tested mortality functions, fitted to data from a wild population of European badgers naturally infected with bovine tuberculosis. Badgers were separated into two health states (“cub-positive”: tested positive in the first year of their lives; “never-positive”: never tested positive throughout their lives). Deviance information criterion (DIC) values are given for each model, as well as corresponding differences (ΔDIC) with respect to the “best” model. Substantial support for the “best” model is indicated when alternative candidate models have ΔDIC > 3 [73].
Table 3. Ranked list of tested mortality functions, fitted to data from a wild population of European badgers naturally infected with bovine tuberculosis. Badgers were separated into two health states (“cub-positive”: tested positive in the first year of their lives; “never-positive”: never tested positive throughout their lives). Deviance information criterion (DIC) values are given for each model, as well as corresponding differences (ΔDIC) with respect to the “best” model. Substantial support for the “best” model is indicated when alternative candidate models have ΔDIC > 3 [73].
Cub-Positive (in Rank Order by DIC)Never-Positive
ModelShapeDIC Δ DICModelShapeDIC Δ DIC
GompertzBathtub46220GompertzBathtub25,6780
GompertzSimple464220ExponentialSimple25,69315
LogisticBathtub466139WeibullBathtub25,69517
WeibullBathtub466947WeibullMakeham25,954276
WeibullMakeham467553LogisticMakeham25,975297
LogisticMakeham468260LogisticSimple25,982304
WeibullSimple468967GompertzMakeham26,004326
LogisticSimple469775LogisticBathtub26,048370
GompertzMakeham471088GompertzSimple26,136458
ExponentialSimple4741119WeibullSimple26,235557

Share and Cite

MDPI and ACS Style

Hudson, D.W.; Delahay, R.; McDonald, R.A.; McKinley, T.J.; Hodgson, D.J. Analysis of Lifetime Mortality Trajectories in Wildlife Disease Research: BaSTA and Beyond. Diversity 2019, 11, 182. https://0-doi-org.brum.beds.ac.uk/10.3390/d11100182

AMA Style

Hudson DW, Delahay R, McDonald RA, McKinley TJ, Hodgson DJ. Analysis of Lifetime Mortality Trajectories in Wildlife Disease Research: BaSTA and Beyond. Diversity. 2019; 11(10):182. https://0-doi-org.brum.beds.ac.uk/10.3390/d11100182

Chicago/Turabian Style

Hudson, Dave W., Richard Delahay, Robbie A. McDonald, Trevelyan J. McKinley, and Dave J. Hodgson. 2019. "Analysis of Lifetime Mortality Trajectories in Wildlife Disease Research: BaSTA and Beyond" Diversity 11, no. 10: 182. https://0-doi-org.brum.beds.ac.uk/10.3390/d11100182

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop