Hostname: page-component-8448b6f56d-tj2md Total loading time: 0 Render date: 2024-04-19T18:41:11.472Z Has data issue: false hasContentIssue false

Molecular epidemiology of Cryptosporidium spp. in calves in Estonia: high prevalence of Cryptosporidium parvum shedding and 10 subtypes identified

Published online by Cambridge University Press:  08 August 2018

Azzurra Santoro
Affiliation:
Department of Veterinary Medicine, University of Perugia, Via San Costanzo 4, 06126, Perugia, Italy
Elisabeth Dorbek-Kolin
Affiliation:
Institute of Veterinary Medicine and Animal Sciences, Estonian University of Life Sciences, Kreutzwaldi 62, 51014 Tartu, Estonia
Julia Jeremejeva
Affiliation:
Institute of Veterinary Medicine and Animal Sciences, Estonian University of Life Sciences, Kreutzwaldi 62, 51014 Tartu, Estonia
Lea Tummeleht
Affiliation:
Institute of Veterinary Medicine and Animal Sciences, Estonian University of Life Sciences, Kreutzwaldi 62, 51014 Tartu, Estonia
Toomas Orro
Affiliation:
Institute of Veterinary Medicine and Animal Sciences, Estonian University of Life Sciences, Kreutzwaldi 62, 51014 Tartu, Estonia
Pikka Jokelainen
Affiliation:
Institute of Veterinary Medicine and Animal Sciences, Estonian University of Life Sciences, Kreutzwaldi 62, 51014 Tartu, Estonia Department of Bacteria, Parasites & Fungi, Infectious Disease Preparedness, Statens Serum Institut, Artillerivej 5, 2300 Copenhagen S, Denmark Faculty of Veterinary Medicine, University of Helsinki, Helsinki, P.O. Box 66, 00014, Finland
Brian Lassen*
Affiliation:
Institute of Veterinary Medicine and Animal Sciences, Estonian University of Life Sciences, Kreutzwaldi 62, 51014 Tartu, Estonia Department of Veterinary and Animal Sciences, University of Copenhagen, Grønnegårdsvej 15, 1870 Frederiksberg C, Denmark
*
Author for correspondence: Brian Lassen, E-mail: brian.lassen@gmail.com

Abstract

We investigated the molecular epidemiology of Cryptosporidium spp. in Estonia by testing fecal samples from 486 calves aged <2 months, raised on 53 cattle farms, for the presence of Cryptosporidium DNA. The parasites were identified and characterized by sequencing of the 18S rRNA gene and of the 60 kDa glycoprotein (gp60) gene. Moreover, using a questionnaire, we surveyed factors that could be relevant for animal-to-human and human-to-animal transmission of Cryptosporidium spp. on the farms. Cryptosporidium spp. were shed by 23% of the investigated calves and at least one shedding calf was found on 66% of the farms. Cryptosporidium parvum was the most common species shed, while C. bovis and C. ryanae were also detected. More than half of the calves aged 8–14 days shed C. parvum. Nine previously described C. parvum subtypes (IIaA14G1R1, IIaA16G1R1, IIaA17G1R1, IIaA18G1R1, IIaA19G1R1, IIaA20G1R1, IIaA21G1R1, IIaA22G1R1 and IIaA16G2R1) and an apparently novel subtype IIlA21R2 were found. Calves from farms that reported spreading manure on fields during spring had 10 times higher odds to shed Cryptosporidium spp. in their feces than calves from farms that did not. Calves aged 8–14 days had higher odds to shed IIa18G1R1 as well as IIaA16G1R1 than younger calves.

Type
Research Article
Copyright
Copyright © Cambridge University Press 2018 

Introduction

Protozoan parasites of the genus Cryptosporidium can cause gastro-intestinal disease in several host species, including humans and cattle (Thompson et al., Reference Thompson, Koh and Clode2016). Cryptosporidium parvum and C. hominis are considered responsible for most cases of human cryptosporidiosis (Cacciò and Chalmers, Reference Cacciò and Chalmers2016). Cryptosporidium hominis is known as human-specific species, while C. parvum has a wider host spectrum that includes cattle. Cryptosporidium parvum has been observed as the dominant Cryptosporidium species shed by pre-weaned calves in many countries (Chako et al., Reference Chako, Tyler, Schultz, Chiguma and Beerntsen2010), but not in countries nearby Estonia (Silverlås and Blanco-Penedo, Reference Silverlås and Blanco-Penedo2013; Björkman et al., Reference Björkman, Lindström, Oweson, Ahola, Troell and Axén2015). Young calves infected with C. parvum can shed high numbers of oocysts in their feces (Xiao, Reference Xiao2010; Smith et al., Reference Smith, Clifton-Hadley, Cheney and Giles2014). In addition to animal-to-human transmission of C. parvum by direct contact, feces of infected cattle may also contaminate, e.g. water supplies (McLauchlin et al., Reference McLauchlin, Amar, Pedraza-Díaz and Nichols2000; Xiao, Reference Xiao2010; Wells et al., Reference Wells, Shaw, Hotchkiss, Gilray, Ayton, Green, Katzer, Wells and Innes2015) or ready-to-eat vegetables (Åberg et al., Reference Åberg, Sjöman, Hemminki, Pirnes, Räsänen, Kalanti, Pohjanvirta, Caccio, Pihlajasaari, Toikkanen, Huusko and Rimhanen-Finne2015).

Cryptosporidium parvum subtype families IIa and IId have been found in both humans and cattle (Xiao, Reference Xiao2010). In Sweden, one of the most common subtype of the C. parvum in cattle was IIaA16G1R1 (Silverlås et al., Reference Silverlås, Näslund, Björkman and Mattsson2010; Björkman et al., Reference Björkman, Lindström, Oweson, Ahola, Troell and Axén2015), which has been also found in humans (Silverlås et al., Reference Silverlås, Näslund, Björkman and Mattsson2010; Insulander et al., Reference Insulander, Silverlås, Lebbad, Karlsson, Mattsson and Svenungsson2013). In Estonia, the same subtype IIaA16G1R1 has been identified in both cattle and an immunocompetent human with clinical cryptosporidiosis (Lassen et al., Reference Lassen, Ståhl and Enemark2014). That case and results of a questionnaire study focusing on veterinary students (Dorbek-Kolin et al., Reference Dorbek-Kolin, Santoro, Tummeleht, Lassen and Jokelainen2018) provide evidence for zoonotic transmission of Cryptosporidium in Estonia.

In Estonia, cryptosporidiosis in humans is a notifiable but under-reported disease (Lassen et al., Reference Lassen, Ståhl and Enemark2014; Plutzer et al., Reference Plutzer, Lassen, Jokelainen, Djurkovic-Djakovic, Kucsera, Dorbek-Kolin, Šobaparl, Sréter, Imre, Omeragic, Nikolic, Bobic, Zivicnjak, Lucinger, Stefanovic, Kucinar, Sroka, Deskne, Keidane, Kvac, Huzova and Karanis2018). Surveillance data do not provide a good overview of the epidemiology of Cryptosporidium, and the need to fill the knowledge gaps with a One Health approach is evident (Plutzer et al., Reference Plutzer, Lassen, Jokelainen, Djurkovic-Djakovic, Kucsera, Dorbek-Kolin, Šobaparl, Sréter, Imre, Omeragic, Nikolic, Bobic, Zivicnjak, Lucinger, Stefanovic, Kucinar, Sroka, Deskne, Keidane, Kvac, Huzova and Karanis2018). Previous studies have shown that almost all Estonian cattle farms had cattle that were shedding Cryptosporidium spp. (Lassen et al., Reference Lassen, Viltrop, Raaperi and Järvis2009). However, the zoonotic potential of the Cryptosporidium spp. shed and the circulating C. parvum subtypes have been unknown.

The main aim of this study was to characterize Cryptosporidium spp. shed by calves in Estonia, with special emphasis on the zoonotic potential. In addition, we surveyed factors that could be relevant for transmission of Cryptosporidium spp. from animals-to-humans or from humans-to-animals on the farms.

Materials and methods

Study design

Sample size calculation was performed using OpenEpi (Dean et al., Reference Dean, Sullivan and Soe2015): 35 farms was the minimum sample size needed for this study. This calculation was based on a population size of 5572 cattle herds (Estonian Agricultural Register and Information Board, 2018a), absolute precision of 10% and an expected proportion of farms with calves shedding C. parvum of 10%. The aim was set to sample at least 50 farms. The sampling was proportionally stratified to the 15 Estonian counties according to the number of farms listed in the Estonian Animal Recording Centre (2013) in each county. Expecting that at a given moment, at least 30% of calves would be shedding Cryptosporidium spp. oocysts on a farm where Cryptosporidium spp. is present (Lassen et al., Reference Lassen, Viltrop, Raaperi and Järvis2009), it was evaluated that 10 calves per farm would be sufficient to find at least one calf shedding the parasite, if Cryptosporidium spp. was present on the farm.

Sampling

The samples were collected by veterinarians from April 2013 to May 2014 and from January to March 2015. Inclusion criteria for farms were: registration in the e-Business Register (Centre of Registers and Information Systems, 2018) and Estonian Agricultural Registers and Information Board (2018b), and herd size ⩾50 cattle to ensure a sufficient number of calves for the study. Farms were selected using a random number generator (Microsoft Excel, Microsoft Cooperation). Three farms were randomly chosen as potential replacements for each county and were included in case a farm that was originally selected opted out.

The veterinarians were instructed to collect individual fecal samples from the rectum of up to 10 calves ⩽2 months of age on each farm. Animal-level exclusion criterion was the calf being reported to be >2 months old. Samples were collected in disposable gloves and stored in a transportable cooler during transport to the laboratory. The samples were stored frozen at −18 °C until DNA extraction.

Questionnaire

A questionnaire was designed to collect information on diarrhoea in calves on the farms as well as on factors with potential relevance for animal-to-human and human-to-animal transmission of Cryptosporidium spp. on farms (Supplementary Table S1). The questionnaire was filled in by the interviewing veterinarian either on the farm at the time of sampling or by phone interview following the farm visit. The questions were asked in Estonian language.

DNA extraction, polymerase chain reaction and sequencing

We used molecular methods to detect, identify and characterize Cryptosporidium spp. from the samples. Genomic DNA was extracted from 200 µg of thawed and homogenized feces using the PSP® Spin Stool DNA Kit (STRATEC Biomedical AG, Birkenfeld, Germany) according to the manufacturer's instructions.

Three microlitres of each DNA sample were submitted to polymerase chain reaction (PCR) amplification targeting the 18S rRNA gene (Xiao et al., Reference Xiao, Morgan, Limor, Escalante, Arrowood, Shulaw, Thompson, Fayer and Lal1999). Nuclease-free water and C. parvum genomic DNA (kindly provided by the European Union Reference Laboratory for Parasites) were used as negative and positive controls. The nested reaction used 1 µL of the first-round PCR product. The thermal cycling conditions were: initial denaturation at 94 °C for 3 min, followed by 35 cycles of 94 °C for 45 s, 55 °C for 45 s, 72 °C for 1 min and a final extension of 72 °C for 7 min in both the first and second rounds. The PCR products were run on 2% ethidium bromide-stained agarose gel and visualized under an ultraviolet transilluminator. Products of the expected size (approximately 825 bp) were submitted to sequencing for species identification.

The samples that tested positive were submitted to PCR amplification targeting the 60 kDa glycoprotein (gp60) gene for subtype identification (Peng et al., Reference Peng, Matos, Gatei, Das, Stantic-Pavlinic, Bern, Sulaiman, Glaberman, Lal and Xiao2001). Three microlitres of DNA sample were used in the first PCR reaction and 1 µL of PCR product in the nested PCR reaction. Nuclease-free water and C. parvum genomic DNA were used as negative and positive controls. Thermal conditions were as follows: initial denaturation at 95 °C for 3 min, followed by 35 cycles of 94 °C for 45 s, 50 °C for 45 s, 72 °C for 1 min and a final extension of 72 °C for 7 min and 10 min in the first and second rounds, respectively. Electrophoresis was performed as described above. Products of approximately 490 bp were selected for subsequent sequencing.

The PCR products were cleaned up and sequenced with Applied Biosystems® 3130xl Genetic Analyzer by a two-directional procedure. Forward and reverse sequences were aligned with BioEdit v7.2.5 software (Hall, Reference Hall1999) to generate single consensus sequences and correct mismatches. The resulting sequences were compared with nucleotide sequences (Accession numbers KJ941147, HQ005736, AM937006, AB242226) deposited in GenBank using BLASTn (nucleotide Basic Local Alignment Search Tool, Altschul et al., Reference Altschul, Gish, Miller, Myers and Lipman1990). Gp60 subtypes were named in agreement with the system proposed by Sulaiman et al. (Reference Sulaiman, Hira, Zhou, Al-Ali, Al-Shelahi, Shweiki, Iqbal, Khalid and Xiao2005) based on the number of serine-coding trinucleotide repeats.

Statistics

Sample size and confidence intervals (CI) (Mid-P exact) were calculated using OpenEpi (Dean et al., Reference Dean, Sullivan and Soe2015). Further statistical analyses were performed using Stata IC 14.2 for Mac software (Stata Corporation, College Station, TX, USA). A calf was considered Cryptosporidium spp.-positive if its sample tested positive for Cryptosporidium spp. 18S rDNA. A farm was considered positive if at least one calf of the investigated calves tested positive. Multivariable logistic regression models were built for dichotomous animal-level outcomes: calf testing positive for Cryptosporidium spp., calf testing positive for C. parvum, calf testing positive for C. parvum subtype IIa18G1R1 and calf testing positive for C. parvum subtype IIaA16G1R1. The farm number was used as a random factor to account for clustering. Variables (Supplementary Table S1) with a P value ⩽0.20 in univariable analysis were first included in the model, followed by a stepwise backward elimination procedure. Biologically meaningful interactions and possible confounding effects were tested. P values <0.05 were considered significant.

Results

Sample

From each of the 53 farms included in the study, 3–14 (median 10) fecal samples were collected, resulting in a total of 522 individual fecal samples. A total of 36 fecal samples were excluded because the same farm had been sampled twice (the samples from second sampling were included), missing labels or insufficient amount of fecal sample available for the analysis. The final sample included in this study comprised individual fecal samples from altogether 486 calves from the 53 farms, 3–14 (median 10) per farm. Information on age was available for 454 calves, and the age of the calves ranged from 1 to 59 days (median 15, mean 18.26). The prevalence estimates were based on the results from 486 calves, i.e. including also the 32 calves with no information of age.

Questionnaire

Of the 53 farms, 49 (92.5%) answered the questionnaire. Supplementary Table S1 shows the distribution of the answers. The majority of farms had more than 150 cattle (79.2%). One (2.1%) farm had bought cattle from abroad during the previous 5 years. Almost a quarter (24.4%) of the farms were located close to natural waterbodies. Altogether, 10 (18.9%) of the farms reported that Cryptosporidium spp. had been diagnosed in calves during the previous 5 years.

Prevalence of Cryptosporidium spp.

Cryptosporidium spp. DNA was amplified and sequenced from 110 (22.63%, 95% CI 19.08–26.51) of the 486 fecal samples. Of the 110 Cryptosporidium spp.-positive fecal samples, 105 (95.45%, 95% CI 90.22–98.32) were C. parvum-positive, four (3.64%, 95% CI 1.17–8.53) were C. bovis-positive and one (0.91%, 95% CI 0.05–4.40) was C. ryanae-positive. Of the 105 C. parvum-positive fecal samples, 95 (90.48%, 95% CI 83.68–95.06) were successfully sequenced and typed by gp60 analysis.

At least one of the investigated calves was Cryptosporidium spp.-positive on 35 (66.0%, 95% CI 52.6–77.8) of the 53 farms (Table 2, Fig. 1). On 33 farms (62.3%, 95% CI 48.7–74.5), at least one of the investigated calves was C. parvum-positive. Cryptosporidium bovis was detected in fecal samples from two farms (3.8%, 95% CI 0.6–11.9), and on one of these farms, C. parvum was also detected. Cryptosporidium ryanae was detected as the only Cryptosporidium species on one farm.

Fig. 1. Map of Estonia showing 18 farms negative for Cryptosporidium spp. (grey circles), three farms positive for Cryptosporidium spp. other than C. parvum (grey diamonds), and 33 farms positive for C. parvum (black circles). We tested fecal samples from three to 14 (median 10) calves per farm, and a farm was considered positive if fecal sample of at least one of the sampled calves tested positive.

Cryptosporidium parvum DNA was almost exclusively found in fecal samples of calves ⩽28 days old, with the exception of one calf that was 36 days old, and seven calves of unknown age (Fig. 2). A total of 64 (52.03%, 95% CI 43.21–60.76) of the 123 calves aged between 8 and 14 days shed C. parvum. Cryptosporidium bovis and C. ryanae were detected in feces of calves that were >14 days old.

Fig. 2. Distribution of the 110 Cryptosporidium spp. shedding calves by age group and the Cryptosporidium species identified.

Cryptosporidium parvum subtypes

A total of 10 different subtypes were identified (Table 1 and 2). The majority (9/10) of the subtypes were in the IIa subtype family, while one subtype was identified as novel IIlA21R2 (Accession numbers MH509210–MH509219). The most common subtype was IIaA18G1R1, which was found in 35.79% (34/95) of the C. parvum-positive samples and on 15 (45.5%) of the 33 farms where C. parvum was found. The second most common subtype was IIaA16G1R1, which was identified in 16 (16.84%) of the C. parvum-positive samples and on four (12.1%) of the C. parvum-positive farms. A single C. parvum subtype per farm was found on all except one farm. Three C. parvum subtypes (IIaA18G1R1, IIaA20G1R1 and IIaA21G1R1) were identified on a farm located in southern part of the country (Fig. 1).

Table 1. Animal-level prevalence of Cryptosporidium species and subtypes in fecal samples from calves (n = 486), including 454 calves ⩽2 months of age and 32 calves of unknown age, collected from 53 cattle farms in Estonia

a Cryptosporidium parvum subtype data for 95 samples (only species level for 10 C. parvum-positive samples).

Table 2. Farm-level prevalence of Cryptosporidium species and subtypes on cattle farms (n = 53) in Estonia

A farm was considered positive if at least one of the 3–14 calves investigated from the farm tested positive.

a One farm had both C. bovis and C. parvum.

b Cryptosporidium parvum subtype data for 31 farms (only species level for two C. parvum-positive farms). A single C. parvum genotype per farm was found on all except one farm. Three C. parvum subtypes (IIaA18G1R1, IIaA20G1R1 and IIaA21G1R1) were found on that one farm.

Subtypes within the IIa family differed by the number of TCA repeats. The exception was the subtype IIaA16G2R1, which had an additional TCG triplet. The novel IIlA21R2 had no TCG triplets, an ACATCA sequence that repeated twice and several single nucleotide polymorphisms when compared with the other IIa subtypes.

Models

Based on the final multivariable model, calves from farms that spread manure on the fields during spring had 10.1 (CI 1.18–86.27) times higher odds to have Cryptosporidium spp. DNA in the feces than calves from farms that did not. Calves that were from farms that reported mortality during the first month of life in calves with severe diarrhoea, which had received veterinary treatment, had 6.2 times higher odds (CI 2.46–15.66) to shed Cryptosporidium spp. (Supplementary Table S2) and 7.4 times higher odds (2.60–21.10) to shed C. parvum (Supplementary Table S3). The odds of a calf aged 8–14 days being Cryptosporidium spp. or C. parvum-positive were 10.1 (CI 4.53–22.36) and 10.4 (CI 4.58–23.74) times higher than the odds of a calf aged up to 7 days, respectively (Supplementary Tables S2 and S3). Based on the final models, the odds of a calf being C. parvum-positive for subtypes IIa18G1R1 and IIaA16G1R1 were 4.00 (P = 0.024) and 25.47 (P = 0.018) times higher in a calf aged 8–14 days than in a calf aged up to 7 days.

Discussion

The results of this study showed that Cryptosporidium spp. were commonly shed by calves in Estonia. This result, which is based on molecular detection, is in line with the previous sample-level and farm-level microscopy-based prevalence estimates (Lassen et al., Reference Lassen, Viltrop, Raaperi and Järvis2009).

In this study, C. parvum was the predominant Cryptosporidium species detected. This is in contrast to reports from calves of comparable age groups from nearby countries Sweden (Silverlås et al., Reference Silverlås, Näslund, Björkman and Mattsson2010; Björkman et al., Reference Björkman, Lindström, Oweson, Ahola, Troell and Axén2015) and Finland (Seppä-Lassila et al., Reference Seppä-Lassila, Orro, Lassen, Lasonen, Autio, Pelkonen and Soveri2015), where the dominant species found were C. bovis, and C. bovis and C. ryanae, respectively. The results of this study resemble the results from Belgium, Slovenia and the Czech Republic, where C. parvum accounted for most of the Cryptosporidium findings from calves (Geurden et al., Reference Geurden, Berkvens, Martens, Casaert, Vercruysse and Claerebout2007; Soba and Logar, Reference Soba and Logar2008; Kváč et al., Reference Kváč, Hromadová, Květoňová, Rost and Sak2011).

Ten different C. parvum subtypes were identified in the fecal samples in this study, indicating a high genetic diversity within C. parvum in Estonia. The most frequently found subtype in this study, IIaA18G1R1, has been reported in cattle in Hungary (Plutzer and Karanis, Reference Plutzer and Karanis2007), the Czech Republic (Kváč et al., Reference Kváč, Hromadová, Květoňová, Rost and Sak2011), the Netherlands (Wielinga et al., Reference Wielinga, de Vries, van der Goot, Mank, Mars, Kortbeek and van der Giessen2008) and Serbia (Misic and Abe, Reference Misic and Abe2007). IIaA16G1R1, the second most commonly found subtype in this study has been reported in the Czech Republic (Ondrácková et al., Reference Ondrácková, Kvác, Sak, Kvetonová and Rost2009; Kváč et al., Reference Kváč, Hromadová, Květoňová, Rost and Sak2011), Hungary (Plutzer and Karanis, Reference Plutzer and Karanis2007), Romania (Imre et al., Reference Imre, Lobo, Matos, Popescu, Genchi and Dărăbuş2011; Vieira et al., Reference Vieira, Mederle, Lobo, Imre, Mederle, Xiao, Darabus and Matos2015), Slovenia (Soba and Logar, Reference Soba and Logar2008) and Sweden (Björkman et al., Reference Björkman, Lindström, Oweson, Ahola, Troell and Axén2015). IIaA17G1R1 has been reported in Poland (Kaupke and Rzeżutka, Reference Kaupke and Rzeżutka2015) and the UK (Smith et al., Reference Smith, Clifton-Hadley, Cheney and Giles2014); IIaA16G2R1 in Belgium (Geurden et al., Reference Geurden, Berkvens, Martens, Casaert, Vercruysse and Claerebout2007) and Spain (Quilez et al., Reference Quilez, Torres, Chalmers, Robinson, Del Cacho and Sanchez-Acedo2008); IIaA14G1R1 in Poland (Kaupke and Rzeżutka, Reference Kaupke and Rzeżutka2015) and Spain (Ramo et al., Reference Ramo, Quílez, Del Cacho and Sánchez-Acedo2014); IIaA20G1R1 in Serbia (Misic and Abe, Reference Misic and Abe2007); and IIaA21G1R1 and IIaA22G1R1 in Sweden (Silverlås et al., Reference Silverlås, Näslund, Björkman and Mattsson2010) and Germany (Broglia et al., Reference Broglia, Reckinger, Cacciò and Nöckler2008). At least six of the subtypes identified in this study from calves, including the two most common ones, have also been found in humans (Soba and Logar, Reference Soba and Logar2008; Chalmers et al., Reference Chalmers, Smith, Hadfield, Elwin and Giles2011; Lassen et al., Reference Lassen, Ståhl and Enemark2014; Lobo et al., Reference Lobo, Augusto, Antunes, Ceita, Xiao, Codices and Matos2014), highlighting the zoonotic potential of C. parvum shed by cattle. Subtype IIaA16G1R1, which has been connected to zoonotic transmission from a calf to a human in Estonia (Lassen et al., Reference Lassen, Ståhl and Enemark2014), was the second most common C. parvum subtype and was identified in feces of 3.3% of the calves and on 7.5% of the investigated farms.

The subtype allele family IIl (also indicated as IIj, Soba and Logar, Reference Soba and Logar2008) has been reported from calves in Serbia (Misic and Abe, Reference Misic and Abe2007), Poland (Kaupke and Rzeżutka, Reference Kaupke and Rzeżutka2015) and Lithuania (Wielinga et al., Reference Wielinga, de Vries, van der Goot, Mank, Mars, Kortbeek and van der Giessen2008), and in calves and humans from Slovenia (Soba and Logar, Reference Soba and Logar2008). To our knowledge, the IIlA21R2 identified in this study is a new subtype.

Cryptosporidium parvum is commonly regarded as a zoonotic pathogen (Cacciò and Chalmers, Reference Cacciò and Chalmers2016). However, not all C. parvum infections in humans result from zoonotic transmission. Molecular subtyping is a useful tool for determining whether human infections originate from animals. The results of studies like ours might be useful for back-tracing potential sources of Cryptosporidium infections and for evaluating the likelihood of the involvement of local cattle in outbreaks. Nevertheless, the definitive discrimination in such cases must resort to a multilocus approach (Chalmers et al., Reference Chalmers, Robinson, Hotchkiss, Alexander, May, Gilray, Connelly and Hadfield2016; Chalmers and Cacciò, Reference Chalmers and Cacciò2016).

Either IIaA16G1R1 (Iqbal et al., Reference Iqbal, Goldfarb, Slinger and Dixon2015) or IIaA16G2R1 (Ranjbar et al., Reference Ranjbar, Baghaei and Nazemalhosseini Mojarad2016), both reported to cause human cryptosporidiosis (Lassen et al., Reference Lassen, Ståhl and Enemark2014; Hijjawi et al., Reference Hijjawi, Zahedi, Kazaleh and Ryan2017), was identified in 4.3% of the fecal samples. These subtypes have also been found in river water in Romania (Imre et al., Reference Imre, Sala, Morar, Ilie, Plutzer, Imre, Hora, Badea, Herbei and Dărăbuș2017) and sewage in Portugal (Lobo et al., Reference Lobo, Xiao, Antunes and Matos2009). More studies are needed on the role of calves in shedding C. parvum in their feces and into the environment. In this study, a calf from a farm that reportedly spread cattle manure on the fields during spring had higher odds of shedding Cryptosporidium spp. as well as C. parvum. In the UK, a peak in human cryptosporidiosis cases caused by C. parvum during springtime was suspected to be linked with livestock (McLauchlin et al., Reference McLauchlin, Amar, Pedraza-Díaz and Nichols2000). In the largest human cryptosporidiosis outbreak, which occurred in Milwaukee in 1993, run-off from cattle farms was suspected as the potential cause of the water contamination that resulted in human C. parvum infections (Mac Kenzie et al., Reference Mac Kenzie, Hoxie, Proctor, Gradus, Blair, Peterson, Kazmierczak, Addiss, Fox, Rose and Davis1994). Several of the subtypes identified in the current study have been reported in wildlife (Krawczyk et al., Reference Krawczyk, van Leeuwen, Jacobs-Reitsma, Wijnands, Bouw, Jahfari, van Hoek, van der Giessen, Roelfsema, Kroes, Kleve, Dullemont, Sprong and de Bruin2015) and fish (Certad et al., Reference Certad, Dupouy-Camet, Gantois, Hammouma-Ghelboun, Pottier, Guyot, Benamrouz, Osman, Delaire, Creusy, Viscogliosi, Dei-Cas, Aliouat-Denis and Follet2015), which adds a sylvatic aspect to the epidemiology. Cryptosporidium spp. infection epidemiology is One Health epidemiology where humans, animals and the environment need to be considered.

High rates of Cryptosporidium spp. infection have been reported in calves of 1–3 weeks of age (Abeywardena et al., Reference Abeywardena, Jex and Gasser2015). In all our models, being 8–14 days old was a risk factor for the calf to shed Cryptosporidium. Young animals are usually more susceptible to Cryptosporidium spp., and may act as amplifiers and infection sources to other animals (Geurden et al., Reference Geurden, Vercruysse and Claerebout2010). In the current study, being 8–14 days old was a risk factor for shedding C. parvum subtype IIaA16G1R1 (OR 20.6), as well as for shedding IIaA18G1R1 (OR 4.0), indicating that this specific age group is a risk group for zoonotic subtypes. This information can be used to design measures that may improve animal health and reduce the occupational risks to humans: considering feces of this age group of calves as likely infective and handling them accordingly could be advisable.

Cryptosporidium parvum can cause high morbidity in calves, and the typical profuse diarrhoea can result in high mortality (Abeywardena et al., Reference Abeywardena, Jex and Gasser2015). Outbreaks with a high mortality in calves due to C. parvum have also been described in Estonia (Lassen and Talvik, Reference Lassen and Talvik2009; Niine et al., Reference Niine, Dorbek-Kolin, Lassen and Orro2017). In the current study, calves from farms reporting mortality of calves with severe diarrhoea that had received veterinary treatment had higher odds to shed Cryptosporidium spp. as well as C. parvum in feces.

The design of this study succeeded in obtaining a well-representative sample from cattle farms all over the country (Fig. 1). We chose not to concentrate the fecal samples before extracting the DNA. Consequently, calves shedding only a few oocysts may have been missed. The results are thus mainly representing calves shedding moderate-to-high numbers of oocysts, and the prevalence estimates should be considered conservative. It should be noted that PCR methods targeting the 18S rRNA and direct sequencing are likely to detect only the most abundant species and genotype in the specimen and underestimate the occurrence of mixed infections (Hadfield et al., Reference Hadfield, Robinson, Elwin and Chalmers2011; Mercado et al., Reference Mercado, Peña, Ozaki, Fredes and Godoy2015).

The gp60 sequence analysis we used is a common approach employed to characterize C. parvum (Xiao, Reference Xiao2010). The findings of this study indicate that subtypes of C. parvum that have also been found in humans were the rule, not the exception, in calves raised in Estonia. It would be important to characterize Cryptosporidium spp. from humans in the country as well, to evaluate the proportion attributable to zoonotic transmission.

Supplementary material

The supplementary material for this article can be found at https://doi.org/10.1017/S0031182018001348.

Acknowledgements

We thank the farms that participated in the study and the European Union Reference Laboratory for Parasites, Istituto Superiore di Sanità, Rome, Italy, for providing us with the positive control. We also thank Piret Kalmus and Kaisa Velström for their contributions in collecting samples, Rachel Chalmers and Annetta Zintl for technical guidance, Pille Paats and Milvi Jallajas for their help in the laboratory, and Heli Kirik for sharing expertise.

Financial support

This work was supported by the Estonian Research Council (T.O., grant number IUT8-1), Estonian Science Foundation (B.L., grant number ETF9433), the Strategic Development Fund of the Estonian University of Life Sciences (P.J., grant number M14143VLVP) and Base Financing of the Estonian University of Life Sciences (P.J., grant number 8P160014VLVP).

Conflict of interest

None.

Ethical standards

The study was approved by the Ethical Committee of Ministry of Agriculture (currently, Ministry of Rural Affairs, permit nr. 7.2-11/46). The participating farms were recruited during farm visits for unrelated reasons, or contacted by phone by a veterinarian. Participation was voluntary and oral informed consent was given.

Footnotes

*

Contributed equally.

References

Åberg, R, Sjöman, M, Hemminki, K, Pirnes, A, Räsänen, S, Kalanti, A, Pohjanvirta, T, Caccio, SM, Pihlajasaari, A, Toikkanen, S, Huusko, S and Rimhanen-Finne, R (2015) Cryptosporidium parvum caused a large outbreak linked to Frisée salad in Finland, 2012. Zoonoses and Public Health 62, 618624.Google Scholar
Abeywardena, H, Jex, AR and Gasser, RB (2015) A perspective on Cryptosporidium and Giardia, with an emphasis on bovines and recent epidemiological findings. Advances in Parasitology 88, 243301.Google Scholar
Altschul, SF, Gish, W, Miller, W, Myers, EW and Lipman, DJ (1990) Basic local alignment search tool. Journal of Molecular Biology 215, 403410.Google Scholar
Björkman, C, Lindström, L, Oweson, C, Ahola, H, Troell, K and Axén, C (2015) Cryptosporidium infections in suckler herd beef calves. Parasitology 142, 11081114.Google Scholar
Broglia, A, Reckinger, S, Cacciò, SM and Nöckler, K (2008) Distribution of Cryptosporidium parvum subtypes in calves in Germany. Veterinary Parasitology 154, 813.Google Scholar
Cacciò, SM and Chalmers, RM (2016) Human cryptosporidiosis in Europe. Clinical Microbiology and Infection 22, 471480.Google Scholar
Centre of Registers and Information Systems (2018) e-Business Register. Retrieved from Centre of Registers and Information Systems website. Available at http://www.rik.ee/en/e-business-register (Accessed 8 March 2018).Google Scholar
Certad, G, Dupouy-Camet, J, Gantois, N, Hammouma-Ghelboun, O, Pottier, M, Guyot, K, Benamrouz, S, Osman, M, Delaire, B, Creusy, C, Viscogliosi, E, Dei-Cas, E, Aliouat-Denis, CM and Follet, J (2015) Identification of Cryptosporidium species in fish from Lake Geneva (Lac Léman) in France. PLoS ONE 10, e0133047.Google Scholar
Chako, CZ, Tyler, JW, Schultz, LG, Chiguma, L and Beerntsen, BT (2010) Cryptosporidiosis in people: it's not just about the cows. Journal of Veterinary Internal Medicine 24, 3743.Google Scholar
Chalmers, R and Cacciò, S (2016) Towards a consensus on genotyping schemes for surveillance and outbreak investigations of Cryptosporidium, Berlin, June 2016. Eurosurveillance 21, 30338.Google Scholar
Chalmers, RM, Smith, RP, Hadfield, SJ, Elwin, K and Giles, M (2011) Zoonotic linkage and variation in Cryptosporidium parvum from patients in the United Kingdom. Parasitology Research 108, 13211325.Google Scholar
Chalmers, RM, Robinson, G, Hotchkiss, E, Alexander, C, May, S, Gilray, J, Connelly, L and Hadfield, SJ (2016) Suitability of loci for multiple-locus variable-number of tandem-repeats analysis of Cryptosporidium parvum for inter-laboratory surveillance and outbreak investigations. Parasitology 144, 3747.Google Scholar
Dean, AG, Sullivan, KM and Soe, MM (2015) OpenEpi: Open Source Epidemiologic Statistics for Public Health. Version 3.03a. Available at http://www.openepi.com (Accessed 8 March 2018).Google Scholar
Dorbek-Kolin, E, Santoro, A, Tummeleht, L, Lassen, B and Jokelainen, P (2018) Occurrence of diarrhoea in undergraduate veterinary students in Estonia, 2010–2016: questionnaire study focusing on cryptosporidiosis and giardiasis, and a case report. Submitted manuscript.Google Scholar
Estonian Agricultural Registers and Information Board (2018 a) Data enquiry retrieved from Estonian Agricultural Registers and Information Board website. Available at https://www.pria.ee (Accessed 8 March 2018a).Google Scholar
Estonian Agricultural Registers and Information Board (2018 b) Retrieved from Estonian Agricultural Registers and Information Board website. Available at https://kls.pria.ee/kaart/ (Accessed 8 March 2018b).Google Scholar
Estonian Animal Recording Centre (2013) Results of Animal Recording in Estonia 2012. Tartu, Estonia: Kirjastus ELMATAR, p. 52.Google Scholar
Geurden, T, Berkvens, D, Martens, C, Casaert, S, Vercruysse, J and Claerebout, E (2007) Molecular epidemiology with subtype analysis of Cryptosporidium in calves in Belgium. Parasitology 134, 19811987.Google Scholar
Geurden, T, Vercruysse, J and Claerebout, E (2010) Is Giardia a significant pathogen in production animals? Experimental Parasitology 124, 98106.Google Scholar
Hadfield, SJ, Robinson, G, Elwin, K and Chalmers, RM (2011) Detection and differentiation of Cryptosporidium spp. in human clinical samples by use of real-time PCR. Journal of Clinical Microbiology 49, 918924.Google Scholar
Hall, TA (1999) Bioedit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symposium Series 41, 9598.Google Scholar
Hijjawi, N, Zahedi, A, Kazaleh, M and Ryan, U (2017) Prevalence of Cryptosporidium species and subtypes in paediatric oncology and non-oncology patients with diarrhoea in Jordan. Infection, Genetics and Evolution 55, 127130.Google Scholar
Imre, K, Lobo, LM, Matos, O, Popescu, C, Genchi, C and Dărăbuş, G (2011) Molecular characterisation of Cryptosporidium isolates from pre-weaned calves in Romania: is there an actual risk of zoonotic infections? Veterinary Parasitology 181, 321324.Google Scholar
Imre, K, Sala, C, Morar, A, Ilie, MS, Plutzer, J, Imre, M, Hora, , Badea, C, Herbei, MV and Dărăbuș, G (2017) Giardia duodenalis and Cryptosporidium spp. as contaminant protozoa of the main rivers of western Romania: genetic characterization and public health potential of the isolates. Environmental Science and Pollution Research 24, 1867218679.Google Scholar
Insulander, M, Silverlås, C, Lebbad, M, Karlsson, L, Mattsson, JG and Svenungsson, B (2013) Molecular epidemiology and clinical manifestations of human cryptosporidiosis in Sweden. Epidemiology and Infection 141, 10091020.Google Scholar
Iqbal, A, Goldfarb, DM, Slinger, R and Dixon, BR (2015) Prevalence and molecular characterization of Cryptosporidium spp. and Giardia duodenalis in diarrhoeic patients in the Qikiqtani Region, Nunavut, Canada. International Journal of Circumpolar Health 74, 27713.Google Scholar
Kaupke, A and Rzeżutka, A (2015) Emergence of novel subtypes of Cryptosporidium parvum in calves in Poland. Parasitology Research 114, 47094716.Google Scholar
Krawczyk, AI, van Leeuwen, AD, Jacobs-Reitsma, W, Wijnands, LM, Bouw, E, Jahfari, S, van Hoek, AH, van der Giessen, JW, Roelfsema, JH, Kroes, M, Kleve, J, Dullemont, Y, Sprong, H and de Bruin, A (2015) Presence of zoonotic agents in engorged ticks and hedgehog faeces from Erinaceus europaeus in (sub) urban areas. Parasites & Vectors 8, 210.Google Scholar
Kváč, M, Hromadová, N, Květoňová, D, Rost, M and Sak, B (2011) Molecular characterization of Cryptosporidium spp. in pre-weaned dairy calves in the Czech Republic: absence of C. ryanae and management-associated distribution of C. andersoni, C. bovis and C. parvum subtypes. Veterinary Parasitology 177, 378382.Google Scholar
Lassen, B and Talvik, H (2009) Parasitic protozoans in livestock and pets in Estonia. Review. Veterinarija ir Zootechnika 46, 3036.Google Scholar
Lassen, B, Viltrop, A, Raaperi, K and Järvis, T (2009) Eimeria and Cryptosporidium in Estonian dairy farms in regard to age, species, and diarrhoea. Veterinary Parasitology 166, 212219.Google Scholar
Lassen, B, Ståhl, M and Enemark, HL (2014) Cryptosporidiosis – an occupational risk and a disregarded disease in Estonia. Acta Veterinaria Scandinavica 56, 13.Google Scholar
Lobo, ML, Xiao, L, Antunes, F and Matos, O (2009) Occurrence of Cryptosporidium and Giardia genotypes and subtypes in raw and treated water in Portugal. Letters in Applied Microbiology 48, 732737.Google Scholar
Lobo, ML, Augusto, J, Antunes, F, Ceita, J, Xiao, L, Codices, V and Matos, O (2014) Cryptosporidium spp., Giardia duodenalis, Enterocytozoon bieneusi and other intestinal parasites in young children in Lobata province, Democratic Republic of São Tomé and Principe. PLoS ONE 9, e97708.Google Scholar
Mac Kenzie, WR, Hoxie, NJ, Proctor, ME, Gradus, MS, Blair, KA, Peterson, DE, Kazmierczak, JJ, Addiss, DG, Fox, KR, Rose, JB and Davis, JP (1994) A massive outbreak in Milwaukee of Cryptosporidium infection transmitted through the public water supply. The New England Journal of Medicine 331, 161167.Google Scholar
McLauchlin, J, Amar, C, Pedraza-Díaz, S and Nichols, GL (2000) Molecular epidemiological analysis of Cryptosporidium spp. in the United Kingdom: results of genotyping Cryptosporidium spp. in 1,705 fecal samples from humans and 105 fecal samples from livestock animals. Journal of Clinical Microbiology 38, 39843990.Google Scholar
Mercado, R, Peña, S, Ozaki, LS, Fredes, F and Godoy, J (2015) Multiple Cryptosporidium parvum subtypes detected in a unique isolate of a Chilean neonatal calf with diarrhea. Parasitology Research 114, 19851988.Google Scholar
Misic, Z and Abe, N (2007) Subtype analysis of Cryptosporidium parvum isolates from calves on farms around Belgrade, Serbia and Montenegro, using the 60 kDa glycoprotein gene sequences. Parasitology 134, 351358.Google Scholar
Niine, T, Dorbek-Kolin, E, Lassen, B and Orro, T (2017) Cryptosporidium outbreak in calves on a large dairy farm: effect of treatment and the association with the inflammatory response and short-term weight gain. Research in Veterinary Science 117, 200208.Google Scholar
Ondrácková, Z, Kvác, M, Sak, B, Kvetonová, D and Rost, M (2009) Prevalence and molecular characterization of Cryptosporidium spp. in dairy cattle in South Bohemia, the Czech Republic. Veterinary Parasitology 165, 141144.Google Scholar
Peng, MM, Matos, O, Gatei, W, Das, P, Stantic-Pavlinic, M, Bern, C, Sulaiman, IM, Glaberman, S, Lal, AA and Xiao, L (2001) A comparison of Cryptosporidium subgenotypes from several geographic regions. The Journal of Eukaryotic Microbiology 48, 2831.Google Scholar
Plutzer, J and Karanis, P (2007) Genotype and subtype analyses of Cryptosporidium isolates from cattle in Hungary. Veterinary Parasitology 146, 357362.Google Scholar
Plutzer, J, Lassen, B, Jokelainen, P, Djurkovic-Djakovic, O, Kucsera, I, Dorbek-Kolin, E, Šobaparl, B, Sréter, T, Imre, K, Omeragic, J, Nikolic, A, Bobic, B, Zivicnjak, T, Lucinger, S, Stefanovic, LL, Kucinar, J, Sroka, J, Deskne, G, Keidane, D, Kvac, M, Huzova, Z and Karanis, P (2018) Cryptosporidium spp. and Giardia intestinalis in Eastern Europe. Eurosurveillance 23, pii=16-00825. doi: https://doi.org/10.2807/1560-7917.ES.2018.23.4.16-00825.Google Scholar
Quilez, J, Torres, E, Chalmers, RM, Robinson, G, Del Cacho, E and Sanchez-Acedo, C (2008) Cryptosporidium species and subtype analysis from dairy calves in Spain. Parasitology 135, 16131620.Google Scholar
Ramo, A, Quílez, J, Del Cacho, E and Sánchez-Acedo, C (2014) Intra-species genetic diversity and clonal structure of Cryptosporidium parvum in sheep farms in a confined geographical area in Northeastern Spain. PLoS ONE 11, e0155336.Google Scholar
Ranjbar, R, Baghaei, K and Nazemalhosseini Mojarad, E (2016) Genetic characterization of Cryptosporidium spp. among patients with gastrointestinal complaints. Gastroenterology and Hepatology from Bed to Bench 9, 301307.Google Scholar
Seppä-Lassila, L, Orro, T, Lassen, B, Lasonen, R, Autio, T, Pelkonen, S and Soveri, T (2015) Intestinal pathogens, diarrhoea and acute phase proteins in naturally infected dairy calves. Comparative Immunology, Microbiology and Infectious Diseases 41, 1016.Google Scholar
Silverlås, C and Blanco-Penedo, I (2013) Cryptosporidium spp. in calves and cows from organic and conventional dairy herds. Epidemiology & Infection 141, 529539.Google Scholar
Silverlås, C, Näslund, K, Björkman, C and Mattsson, JG (2010) Molecular characterisation of Cryptosporidium isolates from Swedish dairy cattle in relation to age, diarrhoea and region. Veterinary Parasitology 169, 289295.Google Scholar
Smith, RP, Clifton-Hadley, FA, Cheney, T and Giles, M (2014) Prevalence and molecular typing of Cryptosporidium in dairy cattle in England and Wales and examination of potential on-farm transmission routes. Veterinary Parasitology 29, 111119.Google Scholar
Soba, B and Logar, J (2008) Genetic classification of Cryptosporidium isolates from humans and calves in Slovenia. Parasitology 135, 12631270.Google Scholar
Sulaiman, IM, Hira, PR, Zhou, L, Al-Ali, FM, Al-Shelahi, FA, Shweiki, HM, Iqbal, J, Khalid, N and Xiao, L (2005) Unique endemicity of cryptosporidiosis in children in Kuwait. Journal of Clinical Microbiology 43, 28052809.Google Scholar
Thompson, RCA, Koh, WH and Clode, PL (2016) Cryptosporidium – what is it? Food and Waterborne Parasitology 4, 5461.Google Scholar
Vieira, PM, Mederle, N, Lobo, ML, Imre, K, Mederle, O, Xiao, L, Darabus, G and Matos, O (2015) Molecular characterisation of Cryptosporidium (apicomplexa) in children and cattle in Romania. Folia Parasitologica (Praha) 62, 002.Google Scholar
Wells, B, Shaw, H, Hotchkiss, E, Gilray, J, Ayton, R, Green, J, Katzer, F, Wells, A and Innes, E (2015) Prevalence, species identification and genotyping Cryptosporidium from livestock and deer in a catchment in the Cairngorms with a history of a contaminated public water supply. Parasites & Vectors 8, 66.Google Scholar
Wielinga, PR, de Vries, A, van der Goot, TH, Mank, T, Mars, MH, Kortbeek, LM and van der Giessen, JW (2008) Molecular epidemiology of Cryptosporidium in humans and cattle in The Netherlands. International Journal for Parasitology 38, 809817.Google Scholar
Xiao, L (2010) Molecular epidemiology of cryptosporidiosis: an update. Experimental Parasitology 124, 8089.Google Scholar
Xiao, L, Morgan, UM, Limor, J, Escalante, A, Arrowood, M, Shulaw, W, Thompson, RC, Fayer, R and Lal, AA (1999) Genetic diversity within Cryptosporidium parvum and related Cryptosporidium species. Applied and Environmental Microbiology 65, 33863391.Google Scholar
Figure 0

Fig. 1. Map of Estonia showing 18 farms negative for Cryptosporidium spp. (grey circles), three farms positive for Cryptosporidium spp. other than C. parvum (grey diamonds), and 33 farms positive for C. parvum (black circles). We tested fecal samples from three to 14 (median 10) calves per farm, and a farm was considered positive if fecal sample of at least one of the sampled calves tested positive.

Figure 1

Fig. 2. Distribution of the 110 Cryptosporidium spp. shedding calves by age group and the Cryptosporidium species identified.

Figure 2

Table 1. Animal-level prevalence of Cryptosporidium species and subtypes in fecal samples from calves (n = 486), including 454 calves ⩽2 months of age and 32 calves of unknown age, collected from 53 cattle farms in Estonia

Figure 3

Table 2. Farm-level prevalence of Cryptosporidium species and subtypes on cattle farms (n = 53) in Estonia

Supplementary material: File

Santoro et al. supplementary material

Santoro et al. supplementary material 1

Download Santoro et al. supplementary material(File)
File 50 KB