key: cord-0935152-zexzeflz
authors: Delafosse, Arnaud; Castro-Hermida, José Antonio; Baudry, Christian; Ares-Mazás, Elvira; Chartier, Christophe
title: Herd-level risk factors for Cryptosporidium infection in dairy-goat kids in western France
date: 2006-11-17
journal: Prev Vet Med
DOI: 10.1016/j.prevetmed.2006.07.001
sha: e76bbce895ecc53e1ddb7fe2840101dca1128e13
doc_id: 935152
cord_uid: zexzeflz

We conducted a cross-sectional study of risk factors for herd-level kid positivity for Cryptosporidium parvum oocysts in dairy-goat farms (Deux-Sèvres, western France). From January to March 2003, faeces from a convenient sample of 879 5- to 30-day-old goat kids from 60 herds were examined microscopically after staining with carbol fuschin. Oocyst shedding was scored semi-quantitatively (0 to 4+) allowing us to obtain a cumulative score per herd. Standardized questionnaires with information about management practices were collected in each farm. We found positive kids in 32 of 60 herds (53.3%) and in 142 animals out of 879 (16.2%). We used logistic regression for two risk-factor model: (1) simple positive (case: herd score ≥1+, at least one positive kid in the herd, versus control: herd score = 0), (2) strongly positive (case: overall herd score ≥3+ versus control: herd score <3+). Risk factors associated with simple positive herds were period of sampling compared to the peak of births (After versus Before, OR = 4.2, 95% CI 1.2, 15.3) and practice of kid grouping by age or weight (Yes versus No, OR = 4.4, 95% CI 1.0, 19.1). Risk factors associated with strongly positive herds were period of investigation (February/March versus January, OR = 12.7, 95% CI 2.1, 76.6), exposure to graminaceous plants in forage (OR = 11.6, 95% CI 1.7, 81.0) and type of ventilation in the goat premises (Vertical versus Wind effect, OR = 14.7, 95% CI 2.1, 106.1). No important association was found between kid-management practices and herd positivity. These results suggest a major role of the environment of kids during their first hours of life in the adult-goat premises regarding the transmission of C. parvum infection.

Cryptosporidium parvum infection is one of the principal enteropathogens in neonatal goats (De Graaf et al., 1999) and has been reported in French dairy goats since 1983 (Polack et al., 1983) . In both natural and experimental conditions, cryptosporidial infection is associated with diarrhoea and mortality in kids aged 1-2 weeks (Nagy et al., 1984; Thamsborg et al., 1990; Koudela and Jiri, 1997; Vieira et al., 1997) . As demonstrated in natural conditions in calves, oocyst excretion mainly occurs beginning on day 3-6 postinfection, continuing for 6-9 days and usually persisting at a detectable but asymptomatic level until 1 month of age (Olson et al., 1997; Uga et al., 2000; Castro-Hermida et al., 2002) . However, clinical cryptosporidiosis occasionally occurs in goats >4 week old (Johnson et al., 1999) . Because no fully satisfactory chemoprophylaxis is available to control neonatal cryptosporidiosis (De Graaf et al., 1999) , a better knowledge of the main risk factors of C. parvum infection in kids is needed to allow the implementation of preventive hygienic measures both reducing the environmental oocyst burden and preventing the transmission to neonates. Such epidemiological information is relatively scarce in dairy-goat farms compared to farms with other ruminants except for some data from Spain (Matos-Fernandez et al., 1994) and from Poland (Majewska et al., 2000) . Our aim was to assess herd-level risk factors for two different levels of oocyst shedding from kids.

We did this study between January and March, 2003.

The target population consisted of all dairy-goat farms of Deux-Sèvres, western France. The sampling frame contained herds from 850 owners (90% of all dairy-goat farms in Deux-Sèvres in 2003) who had a kidding period in winter.

A computed-generated list of 60 farms was selected (Epi Info Version 6.04) to be 95% confident that estimated herd-level prevalence was 40% AE 12.

Fifteen 5-to 30-day-old goat kids were selected from each herd as a convenience sample, based on a 95% confidence in selecting at least one positive kid, should the withinherd prevalence be 20% and test sensitivity (Se) be 50% (in Toma et al., 1996) .

Faecal samples were collected directly from the rectum, identified and examined by microscopy after staining the smear with carbol fuschin. Slides were examined by a single experienced microscopist who did not know the kid's values on any of the other variables. Oocysts were counted in 20 fields (1000Â) and samples were recorded as negative when no oocysts was found.

Cryptosporidium oocyst shedding was scored semi quantitatively (0 to 4+) according to the number of oocysts per microscopic field: 0: absence of oocyst in 20 fields, 1+: <1 oocyst/ field, 2+: 1-10 oocysts/field, 3+: 11-20 oocysts/field, 4+: >20 oocysts/field (Heine, 1982) .

Although some faecal samples contain small numbers of non C. parvum particulates (Maldonado-Camargo et al., 1998) , specificity (Sp) of Heine technique was assumed to have perfect animal-level specificity (leading also to perfect herd-level specificity). In contrast, the animal-level sensitivity is low because the method cannot detect light infections (Farrington et al., 1995) .

Standardized questionnaires with information about management practices were systematically collected. A personal interview of the farm owner/manager was performed during the visit. Factors hypothesized to be associated with the risk of C. parvum infection in kids were selected after a review of the scientific literature (Garber et al., 1994; Atwill et al., 1998; Maldonado-Camargo et al., 1998; De Graaf et al., 1999; Mohammed et al., 1999; Sischo et al., 2000; Huetink et al., 2001; Castro-Hermida et al., 2002) . Factors were divided into three management categories: general, adult-goat and kid-goat management (Tables 1 and 2 ). The questionnaire was tested on five farms (not among the 60) to identify potential sources for misinterpretation of the questions and to further refine the questions. Individual data (age, sex, breed, presence of diarrhoea) were collected for each goat sampled.

Individual results of the Heine technique were added to obtain a score at the herd level (cumulative score). The outcome variables in the study were the case-control status of a herd in two situations: (1) simple positive (case: cumulative herd score !1+, at least one positive kid in the herd, versus control: cumulative herd score = 0); (2) strongly positive (case: cumulative score !3+, i.e. three kids with 1+ or one with 1+ and one with 2+ or one with 3+, versus control: cumulative herd score <3+).

The approach used to deal with large numbers of independent variables was to investigate potential associations fully between the independent and dependent variables. Factors that were associated ( p < 0.10) with the likelihood of infection in the bivariable analysis were further considered in the multivariable analysis. The choice of a screening criterion of 0.10 was designed to limit multicollinearity but to ensure that potentially important variables were included in the next analytical step, the multivariable analysis (Dohoo et al., 1996) .

The bivariable association between each hypothesized factor and each C. parvum infection situation (1, simple positivity; 2, strong positivity) was obtained from Student t-test, Wilcoxon test or Chi-square test depending on whether the independent variable had a continuous distribution (Gaussian or non-Gaussian) or categorical distributions. Normality of the continuous variables was tested using the Shapiro-Wilk test (Royston, 1995) . Multivariable analysis consisted in a backward stepwise logistic regression with a procedure minimising the Akaike Criterion (AIC) (Akaike, 1974) . The stepwise method based on AIC that we used was based on the stepAIC function (Venables and Ripley, 1999) (STEP procedure for glm models, R statistical software). The R function calculates the AIC according to the formula (R Documentation): AIC = À2  log-likelihood + k  npar, where npar represents the number of parameters in the fitted model, and k = 2 for the usual AIC. The starting candidate model was based on all of the predictors with p < 0.10 in the screening. Subsequent models were based on omitting a variable from the current candidate model or adding a variable that was not in the model, with the choice based on minimizing AIC. The final model was found when adding or omitting a variable did not reduce AIC further.

The linearity assumption for continuous variables selected for testing in the models was assessed by categorizing each continuous variable into multiple dichotomous variables of equal units and plotting each variable's coefficient against the midpoint of the variable. We also performed the Chi-squared test for trend in proportions (R statistical software).

For dealing with potential multicollinearity, the bivariable associations between factors retained and not retained in the logistic models were studied from Student t-test, Wilcoxon test, Chi-square test or Rank Correlation test depending on the distributions of the variables (Table 4 ).

In total, 879 goat kids were tested by Heine's technique. Of these, 142 (16.2%) were positive from 32 herds, out of a total of 60 herds sampled (53.3%). The number of goat kids sampled was <15 in 7 farms (10 in 1 farm; 11 in 2 farms; 12 in 1 farm; 13 in 2 farms and 14 in 1 farm). One hundred and three (11.7%), 24 (2.7%), 13 (1.5%) and 2 kids (0.2%) had respectively a score of 1+, 2+, 3+ and 4+.

Twenty three herds (38.3%) were classified as strongly positive (cumulative herd score !3+).

The highest prevalence for diarrhoea (51.9%) was detected in strongly positive herds compared with 17.0 and 8.1% in the low-positive (cumulative herd score <3+) and negative herds, respectively (x 2 = 442, p < 0.001).

Categorical and continuous variables are presented in Tables 1 and 2. Four variables were retained in model 1 (simple positivity) and 5 in model 2 (strong positivity). Only one variable was selected in both models (Table 3) .

Mainly retained factors were general management or adult goat management factors. Herds sampled in February/March (versus January) or After (versus Before) peak of births were at higher risk of infection. Risk also increased with characteristics of goat housing (wall in cement-stone-brick, vertical ventilation, large surface available by goat) and with breeding practices (birth grouping, forage with graminaceous plants) ( Table 3 ). The only kid goat management factor identified was the practice of kid grouping by age or weight (Table 3) .

Although water and kid change were not retained in model 1 (Table 3) , each was associated with another risk factor that was retained (Table 4) ; so, these have not truly been ruled out as having association with simple herd positivity. Similarly, birth peak, autumn kidding, and kid grouping (Table 4) were associated with factors retained in the final version of model 2. 

Extrapolation of our results to kidding in the autumn should be done only with care. Many other agents (Escherichia coli, rotavirus, coronavirus) are involved in diarrhoea of neonatal ruminants. It would have been relevant to include these factors in our statistical analysis. However, this information in goat is extremely scarce (Millemann et al., 2003) . In a few herds, the sample of 15 5-to 30-day-old kids was not achieved with an increased risk of misclassification in the negative group (sensitivity of the diagnosis at the herd level lower than the expected value of 0.95).

One of the objectives of our study was the identification and quantification of risk factors for C. parvum infection at the herd level. Many factors were investigated. In such a situation, the possibility of finding associations 'due to chance alone' goes up substantially (Dohoo et al., 1996) . The other problem is multicollinearity which occurs when predictor variables are not statistically independent and which results in unstable estimates of regression coefficients in logistic-regression models (Dohoo et al., 1996) .

Several variables eliminated by the backward-elimination process were highly collinear with variables retained by the modelling (Table 4 ). The information from the variables ''Birth peak'' and ''Period'' was probably redundant because the birth peak generally occurs before the end of the winter period (February/March). In the same way, ''Autumn kidding'' and ''Birth grouping'' variables included similar information because birth in autumn is allowed by non natural methods (photoperiod control or hormonal treatments) on groups of goats. To simplify our statistical approach, we could have preserved only one of the redundant variables, and discarded the other one, before the backward-elimination process. On the other hand, no obvious relationship appeared between the variables ''Water'' and ''Birth peak'' on the one hand and ''Kid grouping'' and ''Graminaceous'' on the other hand.

Our strategy for reducing the number of independent variables was to screen potential predictor variables using simple (unconditional) statistics and then select a subset of independent variables for inclusion in the final analysis. Another approach would be to create indices or scores which combine data from multiple factors into a single variable or to create a smaller set of independent variables through the use of multivariable techniques (Dohoo et al., 1996) . However, those techniques do not assess the statistical significance of the direct associations between specific independent variables and the dependent variable. Variable-selection methods attempt to balance the goodness-of-fit of a model with considerations of parsimony. However this is done, in theory one could look at all possible logistic-regression models to find the ''best'' one, but this becomes computationally prohibitive when p is large. For this reason, it is more classical to use a stepwise procedure, where candidate models are based on adding or removing a term from the current ''best'' model.

A flaw in this approach is that it does not directly address the crucial balance of goodness-of-fit versus parsimony. Akaike (1974) proposed a measure (AIC) based on information considerations that explicitly quantifies this balance. Model selection using AIC does better, more often resulting in improved performance relative to using the full logistic-regression model, particularly for smaller samples sizes (Perlich et al., 2003) .

We acknowledge that we did not have random sampling at each step of our selection process. However, the specificity was good enough that we are confidant that, where positive kids were found, the herds truly were infected. The problem remaining is that because of uncertain sensitivity and the convenient sampling of kids, there might also be infected herds among those used as controls in our models. The effect, we feel, would have been to diminish the contrast between the control and the positive herds. With such large ORs as we found especially in the model for strong positivity, we believe that we have identified a few risk factors worthy of additional examination.

The period of sampling was linked to the risk of C. parvum infection with February/ March period being more at risk than January (OR = 12.7). Monthly variations in prevalence were also recorded by Bourgoin (1996) and Lefay et al. (2000) in cattle in France and by Causape et al. (2002) in lambs in Spain. One of the likely explanations is the hygienic degradation of the neonates' environment and the simultaneous increase in disease incidence during the whole kidding or calving period especially when parturition is highly seasonal. This is particularly the case in dairy-goat breeding where kidding period mainly extends from November to March with the neonates' premises not being cleaned during this period. Otherwise, changes in crowding of premises or multiple peaks of newborn animals may account for seasonal variations in Cryptosporidium oocyst shedding (Huetink et al., 2001) .

When examining the risk factors of Cryptosporidium infection in farms, none of the parameters related to the management of goat neonates, except the practice of kid grouping by age or weight, was significant in our study. This surprised us, because an extensive previous study of Atwill et al. (1998) performed in a similar system where dairy calves are immediately separated from their dams brought some evidence of the key role of environment in the contamination with Cryptosporidium. Those authors demonstrated the role of the floors and the walls of calf hutches in the transmission of cryptosporidial infection to neonates whereas the dam excretion was undetectable. Similarly, significant relationships between some calf-environment factors such as ventilation, frequency of disposal or addition of bedding, feeding, calf-to-calf contact and the odds of shedding C. parvum oocysts have been described (Maldonado-Camargo et al., 1998; Mohammed et al., 1999; Sischo et al., 2000) .

In contrast, other results suggest an important role of the dam in the direct contamination of their offsprings. According to both the prepatent period of C. parvum infection and the breeding conditions of dairy calves in the Netherlands (immediate removal of calves from their dams, rearing in individual hutches), Huetink et al. (2001) considered that direct or indirect cow-to-calf transmission was probably the most important route of infection followed by indirect calf-to-calf transmission through vectors (caretakers, insects). Asymptomatic adult cattle or sheep can be inapparent carriers especially around parturition when an increase in oocyst output generally occurs (Xiao et al., 1994; Ortega-Mora et al., 1999; Ralston et al., 2003; Faubert and Litvinsky, 2000) . Similarly, >25% of the goats >1 year old shed oocysts (Noordeen et al., 2001) . Thus the risk factors we found (which mainly involved adult-goat management) have to be interpreted in relation with oocyst contamination, i.e output and persistency in adult premises. The vertical ventilation (ridge tile) appeared as a risk factor and this might be related to a more efficient system (compared to horizontal ventilation) leading to a paradoxically more favourable environment for the oocyst survival due to a lower level of ammonia in the litter (Jenkins et al., 1999; Walker et al., 1998) . We also found that the type of the diet given to the adult goat was associated with an increased risk of C. parvum infection prevalence in kids. Forages might be involved through a direct (contaminated feed) or an indirect effect. Very little information is available on survival of oocysts in forage except the data of Merry et al. (1997) on ryegrass ensilage showing a 30% survival after 3 months. Except for contamination with rodent or pet faeces after collection, the survival of occysts on hay is expected to be very low if any because of the dryness of such a material. The impact of graminaceous versus leguminous in the diet on oocyst excretion remains to be investigated.

Knowledge of transmission routes of C. parvum infection remains essential for the development of a control program on farms. Our results suggest a major role of the environment of the first hours of kids life regarding the transmission of C. parvum infection in the adult-goat premises.

A new look at statistical model identification

Evaluation of periparturient dairy cows and contact surfaces as a reservoir of Cryptosporidium parvum for calfhood infection

La place de la cryptosporidiose dans les maladies néo-natales du veau en Corrèze

Prevalence and analysis of potential risk factors for Cryptosporidium parvum infection in lambs in Zaragoza (northeastern Spain)

A study of cryptosporidiosis in a cohort of neonatal calves

A review of the importance of cryptosporidiosis in farm animals

An overview of techniques for dealing with large numbers of independent variables in epidemiologic studies

Patterns of Cryptosporidium antigen and oocyst excretion in calves studied by reverse passive haemagglutination and light microscopy

Natural transmission of Cryptosporidium parvum between dams and calves on a dairy farm

Potential risk factors for Cryptosporidium infection in dairy calves

An easy technique for the demonstration of Cryptosporidia in faeces

Epidemiology of Cryptosporidium spp. and Giardia duodenalis on a dairy farm

Inactivation of Cryptosporidium parvum oocysts by ammonia

Atypical outbreak of caprine cryptosporidiosis in the Sultanate of Oman

Experimental cryptosporidiosis in kids

Prevalence of Cryptosporidium infection in calves in France

Prevalence of and risk factors for shedding of Cryptosporidium parvum in Holstein Freisian dairy calves in central México

Prevalence of Cryptosporidium in sheep and goats bred on five farms in west-central region of Poland

Epidemiologica de la criptosporidiosis en el ganado, ovino y caprino de la montana de Leon

Viability of Cryptosporidium parvum during ensilage of perennial ryegrass

Les diarrhées néonatales des agneaux et des chevreaux

Risk factors associated with Cryptosporidium parvum infection in dairy cattle in southeastern New York State

Studies on cryptosporidial infection of goat kids. Les maladies de la chèvre

Excretion of Cryptosporidium oocysts by goats in relation to age and season in the dry zone of Sri Lanka

Giardia and Cryptosporidium in dairy calves in British Columbia

Role of adult sheep in transmission of infection by Cryptosporidium parvum to lambs: confirmation of periparturient rise

Tree induction versus logistic regression: a learning-curve analysis

Les cryptosporidies en France. Techniques usuelles d'identification et résultats préliminaires d'enquêtes épidémiologiques. Point Vet

Prevalence and infection pattern of naturally acquired giardiasis and cryptosporidiosis in range beef calves and their dams

A remark on algorithm AS 181: the W test for normality

Cryptosporidia on dairy farms and the role these farms may have in contaminating surface water supplies in the northeastern United States

Cryptosporidiosis in kids of dairy goats

Epidémiologie appliquée à la lutte collective contre les maladies transmissibles majeures

Prevalence of Cryptosporidium parvum infection and pattern of oocyst shedding in calves in Japan

Modern Applied Statistics with S-Plus

Outbreak of cryptoporidiosis in dairy goats in Brazil

Source water assessment and nonpoint sources of acutely toxic contaminants: a review of research related to survival and transport of Cryptosporidium parvum

Periparturient rise in the excretion of Giardia sp. cysts and Cryptosporidium parvum oocysts as a source of infection for lambs

We thank Frantz Jénot (Chambre d'Agriculture of Deux-Sèvres, France) for his technical assistance. J.A. Castro-Hermida is supported by Xunta de Galicia (Spain) by means of the Conselleria de Politica Agroalimentaria e Desenvolvemento Rural (PGIDIT02RAG20301 PR) and Post-Doctoral fellowship of the Conselleria de Presidencia e Administracion Publica (PR409A2002/7-0). We are indebted to Conseil Régional Poitou-Charentes (France) for additional financial support.