- Open Access
A multi-event capture-recapture analysis of Toxoplasma gondii seroconversion dynamics in farm cats
Parasites & Vectors volume 11, Article number: 339 (2018)
Domestic cats play a key role in the epidemiology of the parasite Toxoplasma gondii by excreting environmentally-resistant oocysts that may infect humans and other warm-blooded animals. The dynamics of Toxoplasma gondii seroconversion, used as a proxy for primo-infection dynamics, was investigated in five cat populations living on farms.
Serological tests on blood samples from cats were performed every three months over a period of two years, for a total of 400 serological tests performed on 130 cats. Variations in seroconversion rates and associated factors were investigated using a multi-event capture-recapture modelling approach that explicitly accounted for uncertainties in cat age and serological status.
Seroprevalence varied between farms, from 15 to 73%, suggesting differential exposure of cats to T. gondii. In farms with high exposure, cats could become infected before reaching the age of six months. Seroconversion rates varied from 0.42 to 0.96 seroconversions per cat per year and were higher in autumn and winter than in spring and summer.
Our results suggest inter-farm and seasonal variations in the risks of exposure to T. gondii oocysts for humans and livestock living on farms. The paper also discusses the role of young cats in the maintenance of environmental contamination by T. gondii oocysts on farms.
Estimating parasitic infection rates in natural host populations often requires longitudinal surveys with repeated sampling of individuals in order to gather information on their infection history . These surveys present methodological challenges and can result in a range of uncertainties that researchers must deal with. The first source of uncertainty is inherent to longitudinal surveys, as a previously sampled individual may not be detected, captured or tested at each sampling occasion, resulting in partial observations of its infectious states over time (i.e. incomplete infection history; ). A second source of uncertainty may result from the lack of reliability of certain diagnostic tests or observations in assigning the infectious state of an individual: a given state may be assigned erroneously (false positive or negative), resulting in misclassification and mistaken interpretations [2,3,4]. These are particularly encountered when pathogen exposure is inferred from qualitative measures, such as the presence of specific antibodies . The classification of an individual as ‘infected’ or ‘susceptible’ then results in an approximation of the continuous distribution of the immune response [6, 7].
Recently, modelling methodologies used for population dynamics, such as site-occupancy modelling [4, 8] and multi-event capture-mark-recapture (CMR) , have been proposed to deal with this state uncertainty issue in epidemiologic studies. Multi-event CMR approach especially allows the estimation of infection rates corrected by individuals detection probabilities by relating the way an individual is observed and recorded in the field (observational process) to the unobserved biological states of individuals that is inferred by the model (biological process; [9,10,11]). Parameters of these two processes are estimated using probabilistic statements describing the probability of each possible detection and serological histories combined, based on the likelihood estimation of the whole histories of all the surveyed individuals. Multi-event capture-mark-recapture has especially been used to study the exposure dynamics of black-legged kittiwakes (Rissa tridactyla) to Lyme disease , the immune dynamic of European rabbits (Oryctolagus cuniculus) to myxoma and rabbit hemorrhagic disease virus  and more recently the probability of true Mycobacterium bovis infection in badger (Meles meles) . In this study, we used this approach to investigate the dynamics of Toxoplasma gondii infection in farm cat populations.
Toxoplasma gondii is a ubiquitous parasite that causes toxoplasmosis in humans and other warm-blooded animals . This protozoan is considered of high medical importance as it can cause severe illness that can be life-threatening in immunocompromised individuals or in foetuses when acquired congenitally . Toxoplasma gondii also represents a veterinary issue as a major source of reproductive failure in small ruminants [16,17,18], and a cause of fatal infections in some wildlife species [19,20,21]. Wild and domestic felids are the only known definitive hosts of T. gondii and they play a major role in spreading the parasite by shedding oocysts in faeces. These oocysts become infectious after sporulation in the environment and can survive and remain infectious for months in soil and water [22,23,24]. Accidental ingestion of oocysts contained in water, soil and vegetables is a major source of infection for animal intermediate hosts  as well as for humans [26, 27].
Among definitive hosts, the domestic cat (Felis silvestris catus) plays a key role in the epidemiology of the parasite [14, 28], especially the free ranging cats settled in livestock farms (and so called “farm cats”), since these farms are considered as hot-spots of environmental contamination by T. gondii in rural areas [29,30,31]. It is believed that cats become infected most often through the consumption of a contaminated prey . After a primary-infection, cats can excrete millions of T. gondii oocysts in their faeces over a period of 7–20 days and then develop a long-lasting humoral immune response against the parasite. Bradyzoites developed following infection may persist within cysts in tissues for the life of the host and IgG antibodies probably do as well . In cats with reactivation of chronic toxoplasmosis from immune suppression, IgG titres only rarely increase .
Estimating primary infection rates in cats is directly relevant to predicting the oocyst burden in the environment, and thus in evaluating the infection risk to humans and animals. Most of the studies that have attempted to assess the determinants of infection by T. gondii in natural hosts populations were based on serological data by estimating seroprevalence, i.e. proportion of individuals within a population that demonstrate T. gondii-specific antibody in the serum (see reviews in [14, 28, 32]). However, if seroprevalence provides evidence of past exposure, it does not inform about the timing or frequency of infection. These latter are usually explored by estimating incidence rates . To our knowledge, only [35, 36] have focused on estimating incidence rates in domestic cat populations, but none of these studies focused on farm cats.
In this study, the determinants of T. gondii seroconversion dynamics (used as a proxy of primo-infection dynamics) were investigated in five farm cat populations based on data from a two-year survey. Blood samples were collected and tested every three months to detect specific antibodies against T. gondii. Seroconversion rates were estimated in relation to the cats’ age and gender and to the season using multi-event capture-mark-recapture (CMR) to account for the uncertain age and serological status of cats that were not captured and not tested for T. gondii antibodies. The probabilities of misclassifications in the assignation of serological status (i.e. apparent false-negative and false-positive results) was also estimated since multi-event CMR models do not require a one-to-one correspondence between a test result and the serological status, but rather estimate the probabilities of each test results conditional on each serological status (e.g. a seropositive cat may be tested as having no T. gondii antibodies but with a different - much lower - probability than a seronegative cat). Although these values do not indicate the intrinsic bias due to the performance of diagnostic test, they allow to examine the apparent misclassification bias and its related factors.
Study site and cat populations
The study was conducted on five populations of domestic cats living on five dairy farms (farms A, B, C, F and T) located in five villages in the Ardennes region of north eastern France (49°27'3.49"N, 4°47'0.7"E to 49°28'09"N, 5°00'52"E; Fig. 1). A recent study conducted on these farms found the soil to be highly contaminated with T. gondii , let expecting a high incidence of T. gondii in host populations. The five studied cat populations were considered as independent units since the sampled farms were between 3.4 and 14 km distant from each other, while the range of domestic cat movement around farms in the study area has been shown not to exceed 2.5 km . The studied cats were free-ranging and mostly dependent on predation to survive as they were not fed or barely fed by farm owners. Their reproduction was not controlled, except on farm A where the first litters of the reproductive period (in spring) were sometimes removed by the farm owner.
Cat capture and anaesthesia
The five cat populations were surveyed from April 2014 to January 2016. Captures were performed four times per year: in January (winter), April (spring), July (summer) and October (autumn) for a total of eight sampling sessions. On each sampling session, 8–10 cage traps (100 × 30 × 30 cm) were baited with dry cat food and distributed over the farmyard for a period of 2–4 days. The number of traps and the duration of the trapping session per farm were determined on the first session, depending on the specific requirements, in order to capture all the individuals in the population (visually observed and/or reported by farm owners). After the first session, these parameters remained unchanged for a given population to ensure similar capture effort over the study period. Cage traps were mainly used for unsociable cats that could not be approached, while sociable cats that could be approached by humans were caught by hand. On each session, captured cats were anaesthetised directly on the farm with isoflurane gas in order to take a blood sample. All cats were individually identified by the colour pattern of their coat or by a coloured plastic collar when having a uniform coat colour. They were also permanently marked with passive integrated transponder (PIT) tags implanted subcutaneously when first anaesthetised. Kittens less than two months old and cats in bad health (e.g. presenting symptoms of upper respiratory illness or extremely thin) were released without being tested. For each captured cat, the age class and gender were recorded. Three age classes were considered: ‘kittens’, from two to three months-old; ‘juveniles’, from four to six months-old; and ‘adults’, over 6 months-old. The age class was estimated based on birth observation or other information collected from farm owners and/or on teeth development . Cat sociability towards humans was also recorded: ‘easy to catch by hand’ or ‘need trap to be caught’. Lastly, information on cat mortality events was also collected from farm owners to complement the dataset.
Assessing cat immunity against Toxoplasma gondii
On each capture, a blood sample was taken from the marginal ear vein of the cat using a sterile needle. The drops of blood were applied on blotting paper (Whatman 3MM CHR) until a 2–3 cm2 surface area was saturated with blood. The blotting paper was then placed in a closed drying box with individual compartments and stored overnight at ambient temperature. When the blood had completely dried, blotting papers were placed in an individual envelope, labelled with the animal’s number and the collection date, and stored at -20 °C for a maximum of one month before performing serological tests. For these tests, a 1 × 1 cm strip of blotting paper was cut out using clean scissors, eluted in 300 μl of phosphate-buffered saline solution (PBS, pH 7.2) and incubated while being agitated overnight. After centrifugation, the supernatant (i.e. the eluate) was collected and tested for the presence of IgG antibodies to T. gondii using the modified agglutination test (MAT). Results were recorded as titres, corresponding to the inverse of the highest dilution for which an agglutination reaction was observed. All serological analyses were performed by the same person to avoid bias in agglutination detection.
Based on experimental results on cats [33, 39, 40], most epidemiological studies have used a cut-off value of 20 or 25 to discriminate between negative and positive results obtained from MAT tests (see reviews in [28, 32]). However, false negatives based on the cut-off value have been observed in experimentally infected cats even though the parasite was still present in their tissues . Furthermore, the elution step performed for serological tests on blotting paper can weakly dilute antibody quantities, which may result in antibody titres inferior to those normally obtained with a standard serum [41, 42] and a misclassification in the assignment of serological status, in particular, producing false-negative results (titre < 25 from a seropositive individual). As multi-event models allow a serological status to give rise to any titre but with different probabilities , three titre classes were considered in the models: ‘titre = 0’ (most likely with seronegative individuals), ‘titre ≥ 25’ (most likely with seropositive individuals) and ‘0 < titre < 25’ (intermediate titres).
As per [35, 36], T. gondii seroprevalence in the five cat populations was estimated over the study period by randomly selecting one antibody measurement for each cat. The proportion of “seropositive” cats and the associated Wilson’s confidence intervals at 95% were calculated considering two cut-off values: (i) titre = 25 (the value generally used with serum samples); and (ii) titre = 10 (one titre lower than the value generally used with serum samples). In order to simulate a large number of the possible combinations using one measurement per cat, this calculation was repeated 1000 times. Overall seroprevalence in a population was then estimated as the mean proportion of seropositive cats in the 1000 re-samplings.
Multi-event capture-mark-recapture modelling
Because unsociable cats were not captured during some sampling sessions, we performed analyses under the general framework of multi-event capture-mark-recapture models (CMR) accounting for uncertainties about state of some individuals (sensu ‘partial observations’) and for heterogeneity in capture rates between sociable and unsociable cats. Two sources of uncertainty in state assessment were observed in the dataset and considered in the models: (i) uncertainty about the serological status of some re-sighted individuals whose blood was not collected and thus not tested for the presence of T. gondii antibodies; and (ii) uncertainty about the age class of some individuals due to incomplete information about birth date and teeth. These uncertainties were considered explicitly in the capture histories by taking into account 15 events defined as follows: 0, not detected; 1, kitten tested with titre = 0; 2, kitten tested with titre ≥ 25; 3, kitten detected but not tested; 4, juvenile tested with titre = 0; 5, juvenile tested with titre ≥ 25; 6, juvenile detected but not tested; 7, adult tested with titre = 0; 8, adult tested with 0 < titre < 25; 9, adult tested with titre ≥ 25; 10, adult detected but not tested; 11, individual of undetermined age tested with titre = 0; 12, individual of undetermined age tested with titre ≥ 25; 13, individual of undetermined age detected but not tested;14, cat recovered dead.
The events ‘kitten/juveniles/undetermined age tested with 0 < titre < 25’ were not considered since uncertainty about age class only concerned the assignment of kittens or juveniles and both of these classes presented titre > 25 when ≠ 0. The event ‘0, not detected’ could correspond to a living cat present on the farm but not observed, a living cat absent from the farm or a dead undetected cat.
The 15 events might correspond to six biological states: ‘seronegative kitten’, ‘seropositive kitten’, ‘seronegative juvenile’, ‘seropositive juvenile’, ‘seronegative adult’ or ‘seropositive adult’. In addition, as information about dead cat recoveries was available from farm owners, the states ‘newly dead’ and ‘dead’ were also considered.
The unobserved biological process in CMR models is described by initial state probability (I), which is the estimated probability of an individual being in a specific state when first captured, and conditional transition probability (Φ), which is the estimated probability of an individual moving from one state to another between two successive sessions following its first capture (Fig. 2a). We decomposed initial states in two parameters: (i) initial age class (Ia), corresponding to the probability of an individual belonging to one of the three age classes when first captured; and (ii) initial serological status (Iss), corresponding to the probability of an individual being seronegative or seropositive when first captured (Fig. 2a and see Additional file 1 for details regarding the matrices constructed for parameters estimation). Conditional transitions were broken down into three steps corresponding to three parameters: (i) survival (S), the estimated probability of an individual surviving between two successive sessions; (ii) growth (c), the estimated probability of an individual moving to the next age class between two successive sessions (conditional on its being alive); and (iii) seroconversion (Ψ), the estimated probability of an individual moving from the state ‘seronegative’ to the state ‘seropositive’ between two successive sessions, i.e. during a three months period (Fig. 2a). The three cat age classes (2–3 months, 4–6 months, and > 6 months) were chosen in such a way that an individual recorded in a given age class at session i systematically moved to the next age class at session i+1 (which occurred three months after session i), except for adult cats which remain in the same age class. Consequently, growth rates were considered constant and fixed to 1 (c = 1). Furthermore, as T. gondii antibodies have lifetime persistence in domestic cats , the transition probability from ‘seropositive’ to ‘seronegative’ was fixed at zero (see Additional file 1 for details regarding the matrices constructed for parameters estimation).
Event probabilities conditional on underlying states (B) were broken down into four steps (Fig. 2b and see Additional file 1 for details regarding the matrices constructed for parameters estimation): detection (β), the estimated probability of an individual being detected (i.e. observed or captured) at a given session; age assignment (α), the estimated probability of an individual being correctly assigned to a given age class at a given session when detected; blood collection (η), the estimated probability of collecting a blood sample from an individual and testing it for the presence of T. gondii antibodies (conditional on being detected); and serological result (τ), the estimated probability of an individual being assigned to a given class of titre values (conditional on blood collection and testing). This last parameter especially allows to estimate the probability a cat assigned as seropositive by the model shows null or low antibody titres (i.e. apparent false negative results) and the probability a cat assigned seronegative by the model shows positive antibody titres (i.e. apparent false positive results).
Cat gender (gender) and sociability towards humans (soc) were integrated in the models as individual covariates. The gender of 23 kittens that died prematurely before being captured was randomly attributed based on a sex ratio at birth of 1:1 . It has been previously reported that survival in free-ranging cats varies according to gender . The variable gender was then tested on survival, as well as on detection and blood collection rates since sexual behaviour may affect the probability of cat presence on farms depending on the season . Gender was also tested on initial serological status and seroconversion rates, although only a few studies have found a consistent relationship between gender and T. gondii seroprevalence in cat populations [46,47,48]. As the prevalence of T. gondii infection in cats and oocyst excretion vary according to season [49,50,51], this effect (season) was tested on seroconversion rates, as well as on survival, detection and blood collection rates. In addition, a group effect of the farm (farm) was also tested on almost all the parameters. This factor should reflect the effects of the farm’s environment, management and cat population characteristics known to influence the infection risk for animals on farms . Finally, as the quantity of blood collected on the blotting papers was increased after the third sampling session, thanks to acquired technical skills, and potentially resulting in an increase in MAT accuracy to detect T. gondii antibodies on blotting papers, an effect of the period (before the third session versus after the third session) was tested on the serological results obtained. All the variables tested on the estimated parameters are summarised in Fig. 2.
Goodness-of-fit tests and model selection
No specific test exists to detect a structural lack of fit of multi-event models to data . In this study, the goodness-of-fit of the multi-event model was assessed after pooling all ‘detected alive’ events, resulting in a multi-state dataset containing only alive and dead encounters and ignoring uncertain states. The goodness-of-fit for the fully time-dependent Jolly-Move (JMV) multi-state model  was then applied using the U-CARE program version 2.3.2 .
Multi-event models were fitted using the E-SURGE program version 1.9.0 . The model selection procedure was performed following a bottom-up approach (see Additional file 1: Table S1 for details on the starting model and the selection procedure). First, the best structure for initial states, i.e. detection, age assignment, blood sample and survival was selected. The last selection was conducted for serological results and seroconversion rates as these were the parameters of primary interest. In a final check, all models neighbouring the model with the best fit were examined. The relative support of competing models was assessed using the Akaike information criterion corrected for small sample size and possible overdispersion (QAICc) . Models with the lowest QAICc were retained as good candidates, based on the criteria that a model fits the data significantly better when having two QAICc-points less than competing models . If the difference in QAICc between models was less than two, the one with the smallest number of parameters was retained.
Model assumptions and starting model structure
Rather than considering all possible models, only those related to a priori biological hypotheses were specified and compared. The following assumptions were used to build the starting model: (i) initial age (Ia) structure varies according to time t (i.e. sampling session); a higher proportion of young individuals in a population was expected from sessions performed during the cat breeding season (March-October) [45, 57]; (ii) initial serological status (Iss) and seroconversion rates (Ψ) vary between farms due to differences in soil contamination with T. gondii ; (iii) survival rates (S) vary according to cat age class; generally being lower for kittens than for juveniles and adults [44, 58]; (iv) survival rates (S) are considered equal for seronegative and seropositive individuals since T. gondii infection rarely induces acute toxoplasmosis and death in domestic cats ; (v) detection (β) and blood collection (η) rates are not affected by the serological status (ss) of individuals but vary with age, since kittens and juveniles are more systematically detected and caught for blood analysis and sociability towards humans since cats caught by hand can be detected and submitted to blood analysis at each session, contrary to unsociable cats; and (vi) serological results (τ) vary according to the inferred serological status (ss); we considered a seropositive cat has a greater probability to show positive antibody titres than a seronegative one and that a seronegative cat has a greater probability to have no antibody than a seropositive one. Based on these assumptions, the starting model had the following structure: Ia(t), Iss(farm), S(age), c, Ψ(farm+soc ad ), β(age), α, η(age+soc ad ), τ(ss), where parentheses indicate an effect of the variables on the parameter and soc ad is the effect of sociability on adults only.
The size of the cat populations varied between 7–25 individuals per farm depending on the period (i.e. kitten season or not). The average number of cats (± SD) observed per farm and per session varied from 10.37 ± 3.25 to 21.62 ± 4.41 (Table 1). A total of 154 cats were detected (observed and/or captured) during the study period. Of these cats, 130 (84.41%) were tested at least once for the presence of T. gondii-specific antibodies. The remaining 24 individuals were never tested because they were kittens that died before weaning. In total, 400 serological tests were carried out (one to eight per individual; see Table 1 for details). During the study period, 91 cats were born: 32 on farm A, 22 on farm B, 16 on farm C, 17 on farm F and 4 on farm T. The age was known to the nearest week for 79 of these individuals (86.81%). The 12 cats for which it was impossible to assign an age class with certainty represented 13.19% of the cats born during the study and 7.79% of the total cats.
Observed antibody titres from the MAT
Transitions from ‘titre = 0’ to ‘titre ≠ 0’ (i.e. intermediate titre or titre ≥ 25) were observed for 28 individuals (see Additional file 2: Table S3 for details of the data): 15 young cats born during the study (for which post-transition titres varied from 25 to 12,800) and 13 adult cats (for which post-transition titres varied from 10 to 2500). Of the 91 cats born during the study, 15 (16.48%) showed titres ≥ 25 at the first capture. Of these 15 cats, two were less than three months and showed antibody titres of 400 and 1600, respectively, nine were 4–6 months old, three could not be assigned in an age class and one was adult. Eleven of these cats originated from farm A (34.37% of kittens born during the study on this farm) and four from farm F (23.53% of the kittens born during the study on this farm).
Based on a cut-off value of titre = 25, the estimated seroprevalence over the study period varied from 15.38% in farm T to 73.08% in farm F (see Table 2). When a cut-off value of titre = 10 was considered, seroprevalence varied from 29.63% in farm B to 73.08% in farm F (Table 2).
The fully time-dependent Jolly-Move model, which ignores uncertainty and serological status, did not fit the data (χ2 = 36.04, df = 21, P = 0.02). Based on the details of the tests (test 3G.Sm: χ2 =19.67, df = 6, P = 0.003), we suspected that a major cause of lack of fit could be a rather strong heterogeneity of detection among individuals. When sociability towards humans (‘caught only with traps’ versus ‘caught by hand’) was taken into account, there remained no evidence of structural failure in the application of the fully time-dependent Jolly-Move model (χ2 = 21.19, df = 19, P = 0.33 and χ2 = 13.77, df = 9, P = 0.13 for group 1 and 2, respectively). An effect of sociability on detection rate was thus considered in the starting model. Cat age was also a variable which may favour the detection of heterogeneity in our dataset, since cats ≤ 6 months were almost systematically detected at each occasion contrary to adults. This effect was taken into account in the models and tested on each parameter.
Seroconversion rates estimated by multi-event CMR models
To investigate seroconversion rates, 95 competing models were performed (Additional file 1: Table S1). In the most supported model, estimated quarterly seroconversion rates (Ψ) varied according to age class, farm and season (Table 3 and Additional file 1: Table S1). A general seasonal pattern was observed on all the farms: quarterly seroconversion rates were higher from October to April than from April to October (Fig. 3). Seroconversion rates and their distribution among age classes also varied according to farm (Fig. 3 and Additional file 1: Table S2 for details of the estimated parameters). These rates were higher on farms A and F than on the other farms and did not differ between age classes (Ψ = 0.27 for the periods April-July and July-October, and Ψ = 0.73 for the periods October-January and January-April). Lower seroconversion rates were observed on farm C and occurred only in cats ≤ 6 months (Ψ = 0.17 and Ψ = 0.59). Seroconversion rates were lowest on farms B and T and occurred only in adult cats (Ψ = 0.04 and Ψ = 0.21). Based on these parameters, the annual seroconversion rates were estimated based on the following formula: Ψannual = 1 - [(1 - ΨJanuary-April) (1 - ΨApril-July) (1 - ΨJuly-October) (1 - ΨOctober-January)] and were equal to 0.42 seroconversions/cat/year on farms B and T, 0.88 seroconversions/cat/year on farm C and 0.96 seroconversions/cat/year on farms A and F (Table 2).
Distribution and estimated occurrence probabilities of antibody titre classes
The distribution of antibody titres ≠ 0 differed between cat populations. On farms C and T, most cats showing titres ≠ 0 exhibited titres between 10 and 100, whereas on farms A, B and F, titres were most often > 100 (see Table 2 and Additional file 2: Table S3 for details of the data). Similar results were also observed in the retained model (Table 3), in which cats assigned as seropositive by the model and living on farms C and T, had a probability of showing intermediate titres (i.e. 0 < titres < 25) of τ = 0.23 and τ = 0.32, and a probability of showing titres ≥ 25 of τ = 0.47 and τ = 0.68, during sessions 1–3 and sessions 4–8, respectively. In contrast, cats assigned as seropositive by the model and living on farms A, B and F, had a probability of showing intermediate titres almost null (τ = 0.05 and τ = 0.06), whereas they had a probability of showing titres ≥ 25 of τ = 0.71 and τ = 0.94, in sessions 1–3 and sessions 4–8, respectively (Fig. 4).
In addition, according to the selected model, antibody titres ≥ 25 or intermediate titres rarely occurred in cats assigned as seronegative by the model. The estimated frequencies of these titres did not reach 0.04 (Additional file 1: Table S2). In contrast, titres = 0 were frequently observed in cats assigned as seropositive by the model at the beginning of the study (τ = 0.24 to 0.30 in sessions 1–3) whereas no titres = 0 were detected after the third session when the quantity of blood collected was increased (i.e. sessions 4–8).
To our knowledge, this study provides for the very first time a detailed description of intra-annual T. gondii seroconversion dynamics in cats living on farms.
Seroprevalences and seroconversion rates among farm cat populations
The average annual seroconversion rates estimated on the five studied farms (0.42 for the farms B and T, 0.88 for the farm C and 0.96 for the farms A and F) are high compared to those previously reported in an urban population of stray cats (0.17 seroconversions/year)  and in three rural populations of owned domestic cats (0.26, 0.36 and 0.39 seroconversions/year) . These results support previous findings showing that livestock farms are hotspots of T. gondii infection in rural areas [29, 31, 59,60,61]. However, the seroconversion rates, as well as the estimated seroprevalences were highly variable between farms (from 15.38 to 73.08% or from 29.63 to 73.08% depending on the selected cut-off value) and the range of seroprevalences we obtained corresponds to the range of seroprevalences in farm cats previously reported in the literature: 33.3% in Spain ; 41.9% ; 68%  and 53.19%  in the USA; 62% in the UK ; and 75% in Poland . The spatial heterogeneity we observed in our study could result from variations in many factors such as: (i) precipitation [35, 36]; (ii) the number of surrounding farms ; (iii) prey availability and access to food ; and/or (iv) the number of kittens in the cat population .
Seasonal variations of seroconversion rates
Seroconversion rates in farm cats vary temporally according to the season, with significantly higher rates during autumn and winter (from 0.21 to 0.73 between October and April) than during spring and summer (from 0.04 to 0.27 between April and October). These variations are in line with the predictions of  suggesting that cat seroprevalence and oocyst excretion are the result of two non-exclusive factors: (i) host population dynamics following intra-annual cycles; and (ii) oocyst survival depending on seasonal climatic variations. Seroconversion probability in cats would be low in spring and summer because newborn kittens, which represent the most susceptible part of the population, are protected from infection by maternally derived antibodies [35, 70]. Furthermore, oocyst survival is expected to be lower during dry and hot periods [71, 72] corresponding to summer in the study area. In contrast, moisture and mild temperatures in autumn, combined with an increasing proportion of infected intermediate host prey (birds and small mammals), would favour the infection risk for cats .
In previous studies carried out in Germany [50, 51], the proportion of owned cat faeces with T. gondii oocysts was lower in winter (December to March) than during the rest of the year while in our study, seroconversion rates in farm cats were high between January and April. The latency between oocyst shedding and seroconversion in cats  may partly explain this difference but it could also originate from differences in prey availability between cat populations. A low availability of infected prey in winter could explained the low risk of infection for free-ranging owned cats in Germany [50, 51], while the concentration of commensal rodent species in farm buildings  could maintain a high risk of farm cats infection during winter in our study area. In addition, oocysts located inside barns may remain viable and infectious during winter since farms are considered to be moist, shady environments favourable to their survival .
Age of infection
Considering the short time span of the survey and the high seroprevalences on farms, surprisingly high seroconversion rates were observed in adult cats on some farms (up to 0.59 on farm C and up to 0.73 in farms A and F). It cannot be excluded that these high seroconversion rates partly result from a previous decrease in immune response in some individuals as has been reported for experimentally infected cats with low T. gondii antibody titres . However, apparent seroconversions may also arise from false-negative results obtained with MAT, especially for the three first sampling sessions when the probability to observe titres = 0 in cats assigned as seropositive by the model were 0.24 to 0.30. As a consequence, some of the transitions from negative to positive titres occurring in adult cats may have been erroneously considered as seroconversions in the models, resulting in an overestimation of the seroconversion rate in cat populations.
In contrast, specific T. gondii antibodies were detected in very young individuals: two kittens aged 2–3 months and one young cat of around three months (assigned in the uncertain age state). Seropositivity in these kittens may result from maternal antibodies which are passively transferred from the colostrum and that can be detected in kitten serum until the age of 10–12 weeks [46, 70]. However, it may also result from an actively acquired infection since titres in the three cats were high (1600, 400 and 3200, respectively) and increased or remained high (3200 for the three individuals) three months later. As a mother starts bringing her kittens dead prey when they are about four weeks of age , kittens can become infected through the consumption of tissue cysts from the end of the weaning period (after eight weeks). For example, in a natural population of stray cats,  reported that five kittens aged 5–9 weeks had acquired T. gondii infection.
Effect of young cats on infection risk in farms
Despite the potential bias reported above (“false” seroconversions in adults and presence of maternal derived antibodies in kittens), different age patterns of seroconversion were observed between farms. These patterns seem to depend on the infection risk on the farm since cats can get infected before reaching six months of age on farms A, F and C, where seroconversion rates and seroprevalence are high, whereas infections only occurred in adults on farms B and T where seroconversion rates and seroprevalences are lower. Infections in young farm cats have been previously reported based on the detection of both T. gondii-specific antibodies  and T. gondii-like oocysts in faeces . The renewal of individuals in farm cat populations is generally high since their reproduction is not, or poorly, controlled by farm owners. This results in a high proportion of susceptible new individuals each year. Only some of these individuals survive to the point of getting infected, as kittens of 2–3 months in our study have an average survival rate of 0.57 (see the results for survival rates in Additional file 3). However, on farms where the risk of exposure to T. gondii is high, young cats can get infected early, so even if they die young they can contribute to environmental contamination by shedding oocysts. In this way, kittens may contribute to the maintenance of oocysts on farms with high T. gondii contamination. In this case, a control strategy of sterilising cats to limit the annual number of new susceptible cats should help reduce environmental contamination.
Variability of serological titres
High heterogeneity in T. gondii antibody titres among seropositive cats has been observed both experimentally  and in natural cat populations . Cats can exhibit four to eight-fold variations in maximal titres after ingesting the same number of bradyzoites . Our study also found that the pattern of antibody titres varied according to the different cat populations, most of the cats on farms B and T showed low to medium titres (mostly from 10 to 100), whereas the cats on farms A, C and F showed high titres (> 100, up to 12,800). As noted in other studies, these variations in antibody titres may be explained both by host factors, such as genotypes (as observed in mice ), health status and/or coinfection with other parasites and virus [76,77,78] and/or by the strain of T. gondii that may impact the level of produced antibodies in cats [33, 39]. However, further studies are needed to understand how immune response varies from a cat to another and at a larger range from a cat population to another.
Based on longitudinal serological data collected in five cat populations every three months over a two-year period, variations in seroconversion rates were observed between farms but also between seasons, suggesting that risks of exposure to T. gondii oocysts for humans and livestock differ during the year. Thanks to the high frequency of blood sampling, our study allowed the detection of seroconversion occurring early in young cats on farms where parasite exposure seems to be greatest. It suggests that kittens, often present in a large number on a farm, may be responsible for contributing to a significant part of oocyst shedding and thus for maintaining high levels of environmental contamination by T. gondii. These results provide important insight for understanding T. gondii infection dynamic in natural cat populations. In addition to these biological findings, this study highlights the utility of multi-event CMR modelling approaches in order to estimate seroconversion rates from partial serological data collected on untamed outdoor cat populations, and the apparent risk of misclassification in serological status based on serological tests.
Modified agglutination test
Cooch EG, Conn PB, Ellner SP, Dobson AP, Pollock KH. Disease dynamics in wild populations: modeling and estimation: a review. J Ornithol. 2012;152(2):485–509.
Conn PB, Cooch EG. Multistate capture-recapture analysis under imperfect state observation: an application to disease models. J Appl Ecol. 2009;46:486–92.
McClintock BT, Nichols JD, Bailey LL, MacKenzie DI, Kendall WL, Franklin AB. Seeking a second opinion: uncertainty in disease ecology. Ecol Letters. 2010;13:659–74.
Lachish S, Gopalaswamy AM, Knowles SCL, Sheldon BC. Site-occupancy modelling as a novel framework for assessing test sensitivity and estimating wildlife disease prevalence from imperfect diagnostic tests. Meth Ecol Evol. 2012;3:339–48.
Gilbert AT, Fooks AR, Hayman DTS, Horton DL, Müller T, Plowright R, et al. Deciphering serology to understand the ecology of infectious diseases in wildlife. EcoHealth. 2013;10:298–313.
Chambert T, Staszewski V, Lobato E, Choquet R, Carrie C, McCoy KD, et al. Exposure of black-legged kittiwakes to Lyme disease spirochetes: dynamics of the immune status of adult hosts and effects on their survival. J Anim Ecol. 2012;81:986–95.
Choquet R, Carrié C, Chambert T, Boulinier T. Estimating transitions between states using measurements with imperfect detection: application to serological data. Ecology. 2013;94:2160–5.
Colvin ME, Peterson JT, Kent ML, Schreck CB. Occupancy modeling for improved accuracy and understanding of pathogen prevalence and dynamics. PLoS One. 2015;10:e0116605.
Pradel R. Multievent, an extension of multistate capture-recapture models to uncertain states. Biometrics. 2005;61:442–7.
Lebreton J-D, Burnham KP, Clobert J, Anderson DR. Modeling survival and testing biological hypotheses using marked animals: a unified approach with case studies. Ecol Monogr. 1992;62:67–118.
Lebreton J-D, Nichols JD, Barker RJ, Pradel R, Spendelow JA. Modeling individual animal histories with multistate capture-recapture models. Adv Ecol Res. 2009;41:87–173.
Santoro S, Pacios I, Moreno S, Bertó-Moran A, Rouco C. Multi-event capture-recapture modeling of host-pathogen dynamics among European rabbit populations exposed to myxoma and rabbit hemorrhagic disease viruses: common and heterogeneous patterns. Vet Res. 2014;45:39.
Buzdugan SN, Vergne T, Grobois V, Delahay RJ, Drewe JA. Inference of the infection status of individuals using longitudinal testing data from cryptic populations: towards a probabilistic approach to diagnosis. Sci Rep. 2017;7:1111.
Tenter AM, Heckeroth AR, Weiss LM. Toxoplasma gondii: from animals to humans. Int J Parasitol. 2000;30:1217–58.
Robert-Gangneux F, Sterkers Y, Yera H, Accoceberry I, Menotti J, Cassaing S, et al. Molecular diagnosis of toxoplasmosis in immunocompromised patients: a 3-year multicenter retrospective study. J Clin Microbiol. 2015;53(5):1677–84.
Dubey JP. Toxoplasmosis in goats. Agri Pract. 1987;8:43–52.
Innes EA, Bartley PM, Buxton D, Katzer F. Ovine toxoplasmosis. Parasitology. 2009;136:1887–94.
Abu-Dalbou MA, Ababneh MM, Giadinis ND, Lafi SQ. Ovine and caprine toxoplasmosis (Toxoplasma gondii). Iran J Vet Sci Tech. 2010;2(2):61–76.
Parameswaran N, Thompson RCA, Sundar N, Pan S, Johnson M, Smith NC, et al. Nonarchetypal Type II-like and atypical strains of Toxoplasma gondii infecting marsupials of Australia. Int J Parasitol. 2010;40(6):635–40.
Shapiro K, Miller M, Mazet J. Temporal association between land-based runoff events and California sea otter (Enhydra lutris nereis) protozoal mortalities. J Wild Dis. 2012;48(2):394–404.
Roe WD, Howe L, Baker EJ, Burrows L, Hunter SA. An atypical genotype of Toxoplasma gondii as a cause of mortality in Hector’s dolphins (Cephalorhynchus hectori). Vet Parasitol. 2013;192:67–74.
Yilmaz SM, Hopkins SH. Effects on different conditions on duration on infectivity of Toxoplasma gondii oocysts. J Parasitol. 1972;58:938–9.
Frenkel JK, Ruiz A, Chinchilla M. Soil survival of Toxoplasma oocysts in Kansas and Costa Rica. Am J Trop Med Hyg. 1975;24:439–43.
Lindsay DS, Dubey JP. Long-term survival of Toxoplasma gondii sporulated oocysts in seawater. J Parasitol. 2009;95:1019–20.
VanWormer E, Fritz H, Shapiro K, Mazet JAK, Conrad PA. Molecules to modeling: Toxoplasma gondii oocysts at the human-animal-environment interface. Comp Immunol Microbiol Infect Dis. 2013;36:217–31.
Muñoz-Zanzi CA, Fry P, Lesina B, Hill D. Toxoplasma gondii oocyst-specific antibodies and source of infection. Emerg Infect Dis. 2010;16:1591–3.
Boyer K, Hill D, Mui E, Wroblewski K, Karrison T, Dubey JP, et al. Unrecognized ingestion of Toxoplasma gondii oocysts leads to congenital toxoplasmosis and causes epidemics in North America. Clin Infect Dis. 2011;53:1081–9.
Dabritz HA, Conrad PA. Cats and Toxoplasma: implications for public health. Zoonoses Public Health. 2010;57:34–52.
Afonso E, Germain E, Poulle M-L, Ruette S, Devillard S, Say L, et al. Environmental determinants of spatial and temporal variations in the transmission of Toxoplasma gondii in its definitive hosts. Int J Parasitol. 2013;2:278–85.
Gotteland C, Chaval Y, Villena I, Galan M, Geers R, Aubert D, et al. Species or local environment, what determines the infection of rodents by Toxoplasma gondii? Parasitology. 2014;141:259–68.
Simon JA, Kurdzielewicz S, Jeanniot E, Dupuis E, Marnef F, Aubert D, et al. Spatial distribution of soil contaminated with Toxoplasma gondii in relation to the distribution and use of domestic cat defecation sites on dairy farms. Int J Parasitol. 2017;47:357–67.
Dubey JP. Toxoplasmosis of Animals and Man. 2nd ed. Boca, Raton: CRC Press; 2010.
Dubey JP, Lappin MR, Thulliez P. Long-term antibody responses of cats fed Toxoplasma gondii tissue cysts. J Parasitol. 1995;81:887–93.
Lappin MR. Update on the diagnosis and management of Toxoplasma gondii infection in cats. Top Companion Anim Med. 2010;25:136–41.
Afonso E, Thulliez P, Gilot-Fromont E. Transmission of Toxoplasma gondii in an urban population of domestic cats (Felis catus). Int J Parasitol. 2006;36:1373–82.
Afonso E, Thulliez P, Gilot-Fromont E. Local meteorological conditions, dynamics of seroconversion to Toxoplasma gondii in cats (Felis catus) and oocyst burden in a rural environment. Epidemiol Infect. 2010;138:1105–3.
Germain E, Benhamou S, Poulle M-L. Spatio-temporal sharing between the European wildcat, the domestic cat and their hybrids. J Zool Lond. 2008;276:195–203.
Courchamp F, Say L, Pontier D. Transmission of feline immunodeficiency virus in a population of cats (Felis catus). Wildl Res. 2000;27:603–11.
Dubey JP, Lappin MR, Thulliez P. Diagnosis of induced toxoplasmosis in neonatal cats. J Am Vet Med Assoc. 1995;207:179–85.
Dubey JP, Thulliez P. Serologic diagnosis of toxoplasmosis in cats fed Toxoplasma gondii tissue cysts. J Am Vet Med Assoc. 1989;194:1297–9.
Simon J. Déterminants spatio-temporels de la transmission de Toxoplasma gondii dans les fermes d'élevage: dynamique d'infection et patron de dépôt des fèces chez un hôte définitif le chat domestique. PhD Thesis. Reims: Université de Reims Champagne-Ardenne; 2017.
Portejoie Y, Faure E, Georges F, Artois M, Peroux R, Guitton JS. Investigation of specific European brown hare syndrome antibodies in wild hares using blood samples dried on blotting paper. Eur J Wildl Res. 2009;55:53.
van Aarde RJ. Demographic parameters of the feral cat Felis catus population at Marion Island. S Afr J Wildl Res. 1983;13:12–6.
Devillard S, Say L, Pontier D. Dispersal pattern of domestic cats (Felis catus) in a promiscuous urban population: do females disperse or die? J Anim Ecol. 2003;72:203–11.
Turner DC, Bateson P. The domestic cat. The biology of its behaviour. 3rd ed. Cambridge: Cambridge University Press; 2013.
Dubey JP, Bhatia CR, Lappin MR, Ferreira LR, Thorn A, Kwok OCH. Seroprevalence of Toxoplasma gondii and Bartonella spp. antibodies in cats from Pennsylvania. J Parasitol. 2009;95:578–80.
Sharif M, Daryani A, Nasrolahei M, Ziapour SP. Prevalence of Toxoplasma gondii antibodies in stray cats in Sari, northern Iran. Trop Anim Health Prod. 2009;41:183–7.
Rahimi MT, Daryani A, Sarvi S, Shokri A, Ahmadpour E, Teshnizi SH, et al. Cats and Toxoplasma gondii: a systematic review and meta-analysis in Iran. Onderstepoort J Vet Res. 2015;82:823.
Salant H, Spira DT. A cross-sectional survey of anti-Toxoplasma gondii antibodies in Jerusalem cats. Vet Parasitol. 2004;124:167–77.
Herrmann DC, Pantchev N, Globokar Vrhovec M, Barutzki D, Wilking H, Fröhlich A, et al. Atypical Toxoplasma gondii genotypes identified in oocysts shed by cats in Germany. Int J Parasitol. 2010;40:285–92.
Schares G, Ziller M, Herrmann DC, Globokar MV, Pantchev N, Conraths FJ. Seasonality in the proportions of domestic cats shedding Toxoplasma gondii or Hammondia hammondi oocysts is associated with climatic factors. Int J Parasitol. 2016;46:263–73.
Guo M, Dubey JP, Hill D, Buchanan RL, Gamble HR, Jones JL, et al. Prevalence and risk factors for Toxoplasma gondii infection in meat animals and meat products destined for human consumption. J Food Prot. 2015;78:457–76.
Pradel R, Wintrebert CMA, Gimenez O. A proposal for a goodness-of-fit test to the Arnason-Schwarz multisite capture-recapture model. Biometrics. 2003;59:43–53.
Choquet R, Lebreton J-D, Gimenez O, Reboulet AM, Pradel R. U-CARE: utilities for performing goodness-of-fit tests and manipulating CApture-REcapture data. Ecography. 2009;32:1071–4.
Choquet R, Rouan L, Pradel R. Program E-Surge: a software application for fitting multievent models. In: Thomson DL, Cooch EG, Conroy MJ, editors. Modeling demographic processes in marked populations, vol. 3. New York: Springer; 2009. p. 845–65.
Burnham KP, Anderson DR. Model selection and multimodel inference: a practical information-theoretic approach. 2nd ed. New York: Springer; 2002.
Bradshaw JWS, Casey RA, Brown SL. The behaviour of the domestic cat. 2nd ed. Wallingford: CAB International; 2012.
Nutter FB, Levine JF, Stoskopf MK. Reproductive capacity of free-roaming domestic cats and kitten survival rate. J Am Vet Med Assoc. 2004;225:1399–402.
Lehmann T, Graham DH, Dahl E, Sreekumar C, Launer F, Corn JL, et al. Transmission dynamics of Toxoplasma gondii on a pig farm. Infect Genet Evol. 2003;3:135–41.
Langlais M, Lélu M, Avenet C, Gilot-Fromont E. A simplified model system for Toxoplasma gondii spread within a heterogeneous environment. Nonlinear Dyn. 2012;68:381–99.
Gotteland C, McFerrin BM, Zhao X, Gilot-Fromont E, Lélu M. Agricultural landscape and spatial distribution of Toxoplasma gondii in rural environment: an agent-based model. Int J Health Geogr. 2014;13:45.
Miró G, Montoya A, Jiménez S, Frisuelos C, Mateo M, Fuentes I. Prevalence of antibodies to Toxoplasma gondii and intestinal parasites in stray, farm and household cats in Spain. Vet Parasitol. 2004;126:249–55.
Smith KE, Zimmerman JJ, Patton S, Beran GW, Hill HT. The epidemiology of toxoplasmosis on Iowa swine farms with an emphasis on the roles of free-living mammals. Vet Parasitol. 1992;42:199–211.
Dubey JP, Weigel RM, Siegel AM, Thulliez P, Kitron UD, Mitchell MA, et al. Sources and reservoirs of Toxoplasma gondii infection on 47 swine farms in Illinois. J Parasitol. 1995;81:723–9.
Dubey JP, Saville WJA, Stanek JF, Reed SM. Prevalence of Toxoplasma gondii antibodies in domestic cats from rural Ohio. J Parasitol. 2002;88:802–3.
Yamaguchi N, Macdonald DW, Passanisi WC, Harbour DA, Hopper CD. Parasite prevalence in free-ranging farm cats, Felis sivestris catus. Epidemiol Infect. 1996;116:217–23.
Sroka J, Wójcik-Fatla A, Szymańska J, Dutkiewicz J, Zajac V, Zwoliński J. The occurrence of Toxoplasma gondii infection in people and animals from rural environment of Lublin region - estimate of potential role of water as a source of infection. Ann Agric Environ Med. 2010;17:125–32.
Lélu M, Langlais M, Poulle M-L, Gilot-Fromont E. Transmission dynamics of Toxoplasma gondii along an urban-rural gradient. Theor Pop Biol. 2010;78:139–47.
Gilot-Fromont E, Lélu M, Dardé M-L, Richomme C, Aubert D, Afonso E, et al. The life cycle of Toxoplasma gondii in the natural environment. In: Djurkovic-Djakovic O, editor. Toxoplasmosis - Recent Advances. Rijeka: InTech; 2012. p. 1–36.
Omata Y, Oikawa H, Kanda M, Mikazuki K, Dilorenzo C, Claveria FG, et al. Transfer of antibodies to kittens from mother cats chronically infected with Toxoplasma gondii. Vet Parasitol. 1994;52:211–8.
Dubey JP. Toxoplasma gondii oocyst survival under defined temperatures. J Parasitol. 1998;84:862–5.
Lélu M, Villena I, Dardé M-L, Aubert D, Geers R, Dupuis E, et al. Quantitative estimation of the viability of Toxoplasma gondii oocysts in soil. Appl Environ Microbiol. 2012;78:5127–32.
Crowell-Davis SL, Curtis TM, Knowles RJ. Social organization in the cat: a modern understanding. J Feline Med Surg. 2004;6:19–28.
Beelitz P, Göbel E, Fauna GR. incidence of endoparasites in kittens and their mothers from different husbandry situations in south Germany. Tierarztl Prax. 1992;20:297–300.
Suzuki Y, Sher A, Yap G, Park D, Neyer LE, Liesenfeld O, et al. IL-10 is required for prevention of necrosis in the small intestine and mortality in both genetically resistant BALB/c and susceptible C57BL/6 mice following peroral infection with Toxoplasma gondii. J Immunol. 2010;164:5375–82.
Dorny P, Speybroeck N, Verstraete S, Baeke M, De Becker A, Berkvens D, et al. Serological survey of Toxoplasma gondii, feline immunodeficiency virus and feline leukaemia virus in urban stray cats in Belgium. Vet Rec. 2002;151:626–9.
Khan IA, Hakak R, Eberle K, Sayles P, Weiss LM, Urban JF. Coinfection with Heligmosomoides polygyrus fails to establish CD8+ T-cell immunity against Toxoplasma gondii. Infect Immun. 2008;76:1305–13.
Stoicov C, Whary M, Rogers AB, Lee FS, Klucevsek K, Li H, et al. Coinfection modulates inflammatory responses and clinical outcome of Helicobacter felis and Toxoplasma gondii infections. J Immunol. 2004;173:3329–36.
The authors would like to thank all the farm owners who permitted access to their properties and their cat populations and who provided very useful information on cat population dynamics. The authors also warmly thank the two anonymous referees for helpful comments on the previous version of the article. All authors are grateful to Alain Simon, Pierre Belle-Isle, Marion Dellinger, Sabine Kurdzielewicz, Céline Sakonda, Candice Gagnaison and Marie Bonneric for their assistance in the field work.
This research was supported by Julie Simon's grant from the French Ministry of National Education, Higher Education and Research.
Availability of data and materials
All data generated or analysed during this study are included in this published article and its additional files.
Ethics approval and consent to participate
All the cat experiments in this study were approved by the French Ministry of National Education, Higher Education and Research and by the Ethics Committee of Champagne-Ardenne (Approval no. 5387) and were conducted according to French regulations on ethics in animal experimentation. Each owner of the respective farms gave consent to perform the study on the property.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Text. Matrices constructed for parameters estimation in the E-SURGE program. Table S1. Detailed structure of the set of models and the model selection procedure from the multi-event analysis. Table S2. Estimated parameters from the selected model. (DOCX 53 kb)
Table S3. Dataset showing Toxoplasma gondii antibody titres obtained from the modified agglutinated test (MAT) and the characteristics of the sampled cat. (DOCX 51 kb)
Description of survival rates, detection rates, age assignment and blood collection rates estimated from the best retained model. (DOCX 12 kb)
About this article
Cite this article
Simon, J.A., Pradel, R., Aubert, D. et al. A multi-event capture-recapture analysis of Toxoplasma gondii seroconversion dynamics in farm cats. Parasites Vectors 11, 339 (2018). https://doi.org/10.1186/s13071-018-2834-4
- Domestic cats
- Multi-event modelling
- Infection dynamics
- Serological titres
- Blotting paper