|
|
||||||||
Department of Population Medicine and Diagnostic Sciences, College of Veterinary Medicine, Cornell University, Ithaca, NY 14853
1 Corresponding author: hom1{at}cornell.edu
| ABSTRACT |
|---|
|
|
|---|
Key Words: prevalence cross-sectional Cryptosporidium risk factor
| INTRODUCTION |
|---|
|
|
|---|
Transmission of Cryptosporidium species takes place via the fecaloral route, and susceptible hosts contract the infection by ingesting the infective form of the parasite, the sporulated oocyst. Human cryptosporidiosis is most often associated with the consumption of contaminated drinking water or exposure to contaminated recreational water, such as public swimming pools. A small number of cases have also been attributed to contaminated foods such as salads (Fayer et al., 2000).
In cattle, C. parvum typically causes disease in calves and has been identified as one of the primary etiologic agents of neonatal calf diarrhea (Naciri et al., 1999). Herd-level prevalence of C. parvum has been reported to range from 13 to 100% (Wade et al., 2000; Santin et al., 2004), although a portion of this variability is likely attributable to study design and the geographical region under investigation. A new species, Cryptosporidium bovis, has recently been described based on both genomic and host range differences (Fayer et al., 2005). To date, this organism has been identified only in cattle and is thus considered unlikely to be zoonotic. Cryptosporidium bovis is indistinguishable from C. parvum using traditional diagnostic methods such as flotation and ELISA techniques; thus, the term "C. parvum-like" has been used to describe infections identified using such methods (Starkey et al., 2005).
Research efforts were increased after Cryptosporidium species were found to be responsible for several large waterborne disease outbreaks in North America, the United Kingdom, and Europe (Patel et al., 1998; Ong et al., 1999; Glaberman et al., 2002). The largest such outbreak occurred in Milwaukee in 1993 (Mac Kenzie et al., 1994). In their 1994 paper, Mac Kenzie et al. speculated that the source of Cryptosporidium oocysts may have been contamination of the body of water from which a failed water treatment plant drew its intakeMilwaukee harborby cattle herds or slaughterhouses located along 2 rivers leading to the harbor. In addition to animal sources, the authors also commented on the possible role of human sewage and mixed-source runoff attributable to spring rains and thawing snow. Subsequent molecular work indicated that C. hominis was responsible for the Milwaukee outbreak (Sulaiman et al., 2001). This finding essentially removed cattle from consideration as the source of the outbreak because C. hominis has never been isolated from cattle, and it is considered an anthroponotic pathogen (Morgan-Ryan et al., 2002). However, C. parvum has been identified at the molecular level in a number of smaller waterborne outbreaks (Ong et al., 1999; Glaberman et al., 2002). This finding, along with the fact that cattle are quite common within watersheds in the United States, has led to further speculation about their role in the contamination of drinking water supplies (Xiao and Ryan, 2004). Furthermore, human doseresponse studies have shown that C. parvum of bovine origin can have an infectious dose 50 (ID50) as low as 87 oocysts (Okhuysen et al., 1999). Therefore, there is little doubt that cattle pose a potential, albeit imperfectly quantified, risk to human health with regard to C. parvum. An additional source of this imperfect risk quantification stems from the relative lack of information about the prevalence of the nonzoonotic C. bovis, especially among cattle populations within watersheds.
Dairy farms represent a significant industry within New York State, with cash receipts from milk totaling over $1.56 billion in 2002, more than half of all agricultural receipts collected that year (New York State Department of Agriculture and Markets, 2006). Such operations are often located within the prime agricultural lands that form watersheds. Academic institutions as well as state and local government agencies have conducted extensive research on the dynamics of infection of Cryptosporidium in dairy cattle as a means to better understand the level of risk these animals pose to the environment. By way of continuation of such research activities, the present study aimed to determine the prevalence of C. parvum in the Upper Susquehanna Watershed and investigate the potential risk factors associated with the prevalence of C. parvum (as identified by traditional methods) in the target population. Additionally, the potential role of cattle in the risk to water quality was examined by differentiating between zoonotic and nonzoonotic isolates using DNA sequencing.
| MATERIALS AND METHODS |
|---|
|
|
|---|
The study population consisted of 19 dairy farms selected from the target population of approximately 100 farms. A block sampling design was adopted to ensure representation of all geographic areas covered by the target population. In this design, farms were grouped into clusters according to their geographic location within the county (i.e., towns). Within this framework, a proportional sampling scheme was also used. Thus, although all areas were represented among the selected farms, more study farms were included from those areas with a relatively higher proportion of dairy operations. While following the above study design, we used a participatory approach for the selection of study farms. This approach involved collaboration with the Delaware County Soil and Water Conservation District and representatives of the local farmers association. The collaborators provided planning and logistical support and facilitated study farm selection by granting access to their ongoing, countywide farm census data. In addition to providing census data, the collaborators also provided detailed local knowledge during the study farm selection phase. We believe that the participatory approach, in conjunction with the proportional block design, led to the selection of a representative study population.
Sample and Data Collection
An age-stratified sampling design was used to collect fecal samples on the study farms (Wade et al., 2000; Santin et al., 2004). Animals less than 6 mo of age were differentially targeted to improve the chances of detecting those shedding C. parvum-like oocysts. In the present study, the fecal sampling protocol called for the sampling on study farms of all calves less than 6 mo of age, up to a maximum of 15 animals. If more than 15 such animals were present, 80% of all calves less than 6 mo of age were to be sampled. Fecal samples were also collected from 5 heifers from 6 mo of age until first freshening, 5 lactating cows, and 5 dry cows at each study farm.
At the time of sampling, study personnel collected detailed farm management and demographic data through the administration of a comprehensive management questionnaire. The questionnaire was conducted in the form of an in-person interview with the farm owner or manager. A single member of the study team administered the questionnaire on all farms to eliminate interobserver error. The questionnaire covered a large range of physical and management factors pertaining to each farm (Table 1
). Emphasis was placed on factors hypothesized to relate to the prevalence of this parasite, such as colostrum management, nature and frequency of bedding changes, nature of manure management, open or closed herd status, and existence of rodent or wild bird problems.
|
Nested PCR
A nested PCR targeting the 18S rRNA gene was used to allow for the differentiation between C. parvum and C. bovis, allowing elucidation of the zoonotic risk posed by cattle within watersheds. All samples testing positive by flotation, ELISA, or both were subjected to the PCR protocol. Additionally, 100 flotation-negative samples were analyzed with this method to determine the degree of agreement between the techniques, as measured by the kappa statistic.
Immediately on return to the laboratory, a 1-mL sub-sample of feces of approximately 50% solids was obtained from each individual fecal sample and stored at 20°C in a sealed tube prior to DNA extraction. Upon thawing the feces, we used a modified mechanical disruption method for DNA extraction (Zhu et al., 1998). Samples were examined at the 18S rRNA gene loci using a nested PCR protocol. Two sets of primers amplifying a final fragment of approximately 830 bp were used. The primary reaction consisted of 1 µL of 1:10 diluted DNA solution obtained from the extraction procedure, added to a mixture consisting of 10.8 µL of reverse osmosis water, 2 µL of 1x PCR buffer, 4.8 µL of MgCl2 (50 mM), 0.4 µL of deoxyribonucleoside triphosphates, 0.4 µL of the forward primer (SSU 1) 5'-GAT AAC CGT GGT AAT TCT AGA GCTA-3' (10 µM), 0.4 µL of the reverse primer (SSU 2) 5'- TAA GGT GCT GAA GGA GTA AGG -3' (10 µM), and 0.2 µL of Taq DNA polymerase. The secondary reaction consisted of 1 µL of the product from the primary reaction added to a mixture consisting of 13.2 µL of water, 2 µL of 1x PCR buffer, 2.4 µL of MgCl2 (50 mM), 0.4 µL of deoxyribonucleoside triphosphates, 0.4 µL of the forward nested primer (SSU 3) 5'-GAA RGG TYG TAT TTA TTA GAT AAA GGAAC -3' (10 µM), 0.4 µL of the reverse nested primer (SSU 4) 5'-AAG GAG TAA GGA ACA ACC TCC A-3' (10 µM), and 0.2 µL of Taq DNA polymerase. Both the primary and secondary reactions were run under the same conditions: 35 cycles of 96°C for 45 s, 55°C for 45 s, and 72°C for 1 min.
DNA Sequencing
Polymerase chain reaction products were purified using Exonuclease I/Shrimp Alkaline Phosphatase (Exo-SAP-IT; USB Corporation, Cleveland, OH). Purified products were sequenced using the internal primers described above in 9-µL reactions using an automated sequencer (3730 DNA Analyzer; Applied Biosystems, Foster City, CA). Samples were sequenced in both directions and sequence chromatograms from each strand were aligned and inspected using MEGA 3.1 (Kumar et al., 2004).
Data Management and Statistical Analysis
Data management, coalition, and validation were performed using Microsoft Access 2000 and Microsoft Excel 2000 (Microsoft, Redbank, CA). The SAS statistical program (Version 9.1 for Windows; SAS Institute, Cary, NC) was used to generate descriptive statistics and manage data. Prevalence was determined as the proportion of samples testing positive for C. parvum-like oocysts relative to those testing negative. Because all samples were analyzed with the flotation method, prevalence is reported based on the results yielded by this method. Prevalence was examined in the population as a whole, at the herd level, and by age.
Fecal flotation results were used in the statistical analyses to identify risk factors associated with the prevalence of C. parvum-like oocyst shedding. These analyses were restricted to animals less than 60 d of age. Such a restriction was used to focus on the population perceived to be at greatest risk, and to increase the power of the statistical analyses. This cut-off was based on the results of a study conducted in an adjacent watershed, which identified an age range of 3 to 60 d among C. parvum-like-positive animals (mean 15 d, standard deviation 6.6 d; Starkey et al., 2005). To this end, a logistic regression model was used to identify and examine the significance of various risk factors when examined simultaneously. During this process, appropriate transformations, including polynomial transformations, were assessed for continuous variables in an attempt to improve model fit. Individual animal and herd management factors obtained from the management questionnaire were first examined for association with the prevalence using appropriate bivariate techniques. Discrete factors were examined by
2 test of independence or Fishers exact test when expected cell counts were less than 5. Continuous variables were examined by unconditional logistic regression. The model-building step continued with the investigation of correlation among significant variables (P
0.2) grouped under each of the management categories (general farm management, maternity management, preweaning calf management, postweaning calf management). Pearson correlation coefficients were examined, and when a correlation was found between variables, the decision regarding which to retain was based on biological plausibility. When 3 or more variables were retained in a management category after the above steps, a best subset selection procedure was run, using Mallows Cp statistic, to determine which to advance to the final model-building step. Age was included as an independent variable in these steps to control for this important risk factor throughout the model-building process. Ultimately, all the significant variables chosen at each management level were submitted, along with age, to a final best subset selection procedure for the selection of the best model of association between parameters of interest and the prevalence of C. parvum-like oocyst shedding.
Because the sampling units (the animals) in this study are clustered into herds, we assumed that this clustering would lead to a correlation in the likelihood of infection within the study population. This correlation between responses occurs because they are dependent on exogenous factors that are associated with these responses (i.e., infection with the organism). Conditioning on an observed set of these factors by controlling for their effect in the analysis and including them as covariates in the logistic regression analysis will sometimes achieve approximate conditional independence. However, more often this correlation in the response arises from both observed and unobserved risk factors. We assumed that the unobserved risk factors were randomly distributed among farms, and the overall significance of this assumption was evaluated by using a mixed-effect logistic regression model (Rosner, 1989). The mixed-effect logistic regression model was specified as follows:
![]() |
where P(CP/
, ß,
) is the probability that an animal within a herd level µi would have been shedding C. parvum-like oocysts given a set of fixed factors Zi with an effect of ßi. The likelihood ratio test was used to evaluate the significance of the farm random effect parameter in the mixed model. The mixed-effect logistic regression analysis was performed using the EGRET statistical software (Cytel Statistical Software, Cambridge, MA). The effect of each factor on the likelihood of infection with the organism was quantified by the odds ratio (OR), which was computed as the exponent of the respective regression coefficient.
| RESULTS |
|---|
|
|
|---|
The average age of animals shedding oocysts was 16 d, with a standard deviation of 9.2 d and a range of 3 to 31 d. When examining animals less than or equal to 60 d of age, prevalence was found to be 20% across the study subjects. The prevalence was higher still when animals of 31 d of age or younger were examined, with 32% of such animals testing positive at the time of sampling.
Risk Factors
After variables within in each management category were examined using bivariate techniques and significant variables were investigated with Pearson correlation coefficients, several variables remained and were subsequently analyzed using a best subset selection procedure (Table 2
). The variables under consideration during the development of the final best subset model were age, type of water supply for the barn, type of housing for milking cows, nature of calf housing (tied vs. not tied), total number of unweaned calves on a farm, and contact of postweaned calves with adult cows (binary).
|
|
| DISCUSSION |
|---|
|
|
|---|
The overall prevalence of 3.9% observed in this study via the flotation method is higher than the 0.9% reported in a cross-sectional study conducted in an adjacent watershed (Wade et al., 2000). Although similar age stratification was used in both study designs, additional attention was paid to the youngest animals within the less than 6-mo-old category in the present study, potentially explaining the increase in observed prevalence. The prevalence of C. parvum-like oocyst shedding among unweaned animals in this study was 20%, a level significantly higher than the overall prevalence and in keeping with other reports in the literature that this age group has the highest prevalence of the organism (Santin et al., 2004). The present study identified a relatively high herd-level prevalence of the organism of 42%. This finding is also in keeping with previous reports in the literature (Garber et al., 1994; Wade et al., 2000).
As indicated by the kappa statistic of 0.71, there was a high degree of agreement between the flotation and the PCR methods used in this study. For the purposes of reporting prevalence and investigating risk factors in the present study, we chose to use the results of flotation. This decision was based on the fact that all samples were subjected to this test and thus all had an equal chance of becoming positive. Additionally, this test remains in common use among practitioners and in diagnostic laboratories throughout the United States. This decision subsequently precluded the investigation of factors that would have led to a differential prevalence of C. parvum and C. bovis in the present study. Such work would require a larger number of positive animals to draw conclusions, and could initially be performed as a case-control study or a larger retrospective cross-sectional study.
For the purpose of model building, we focused our analysis on animals that were less than 60 d of age. The rationale for this strategy is 2-fold. First, from our previous studies we learned that among the populations of cattle in New York State, the risk of infection is limited to younger calves that are typically less than 60 d of age (Wade et al., 2000; Starkey et al., 2005). Second, other researchers found similar results in other cattle populations outside New York State (Santin et al., 2004). Given these findings and the potential for increased statistical power, we opted to restrict the analysis to cattle below 60 d of age.
Among the findings of interest in the final logistic regression model was the estimated 37% increase in the odds of infection among those animals whose barn water source was not a well. In the study population, non-well sources of water were either springs or streams. It is possible that these water sources are associated with an increased prevalence because they may have a greater risk of fecal contamination relative to well water. In the final model, there was also a significant increase in the risk of infection in calves that were housed in a tied manner, with a point estimate indicating a 434% increase in the odds of infection. Most other unweaned calves in the study population were housed in greenhouses. One possible explanation for this increased risk among tied animals is the relative ease of contact with immediate neighbors in most tied environments, thereby increasing the risk of fecaloral contact. However, it is important to note that even moderate physical separation may not eliminate the risk of transmission of this organism owing to the frequently propulsive nature of the diarrhea associated with the infection. The number of unweaned animals on a property was significantly associated with an increase in the risk of C. parvum-like oocyst shedding. This finding is logical, because an increased number of at-risk (i.e., unweaned) animals would be expected to increase both the chances of detecting infection and the likelihood of an uninfected animal encountering an infected animal or infectious oocyst in its environment.
The study was designed to increase both internal and external validity; however, a large-scale longitudinal study needs to be conducted to confirm the identity of the identified risk factors and to allow causal inferences to be made. Considering increasing evidence regarding the relatively high prevalence of C. bovis, further work is also needed to better quantify the risk cattle pose to water supplies and downstream water users. An integrative approach combining stochastic and deterministic reasoning with data in the literature may aid in identifying a useful determination of risk.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication January 27, 2006. Accepted for publication May 8, 2006.
| REFERENCES |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |