# First identification of Cryptosporidium parvum virus 1 (CSpV1) in various subtypes of Cryptosporidium parvum from diarrheic calves, lambs and goat kids from France

Karim Tarik Adjou, Aurélie Chevillot, Pierrick Lucas, Yannick Blanchard, Houria Louifi, Razika Arab, Mohamed Mammeri, Myriam Thomas, Bruno Polack, Grégory Karadjian, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Karim Tarik Adjou, Aurélie Chevillot, Pierrick Lucas, Yannick Blanchard, Houria Louifi, et al. First identification of Cryptosporidium parvum virus 1 (CSpV1) in various subtypes of Cryptosporidium parvum from diarrheic calves, lambs and goat kids from France. Veterinary Research, 2023, 54 (1), pp.66. 10.1186/s13567-023-01196-4. anses-04185936

HAL Id: anses-04185936

https://anses.hal.science/anses-04185936

Submitted on 23 Aug 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

Open Access

<figure><img src="image_3.png" /></figure>

Check for
updates

# First identification of Cryptosporidium parvum virus 1 (CSpV1) in various subtypes of *Cryptosporidium parvum* from diarrheic calves, lambs and goat kids from France

Karim Tarik Adjou¹*, Aurélie Chevillot¹, Pierrick Lucas³, Yannick Blanchard³, Houria Louifi¹, Razika Arab¹, Mohamed Mammeri¹, Myriam Thomas¹, Bruno Polack¹, Grégory Karadjian¹ and Nolwenn M. Dheilly²*

## Abstract

*Cryptosporidium* spp. remain a major cause of waterborne diarrhea and illness in developing countries and represent a significant burden to farmers worldwide. Cryptosporidium parvum virus 1 (CSpV1), of the genus *Cryspovirus*, was first reported to be present in the cytoplasm of *C. parvum* in 1997. Full-length genome sequences have been obtained from *C. parvum* from Iowa (Iowa), Kansas (KSU) and China. We aimed at characterizing the genome of CSpV1 from France and used sequence analysis from *Cryptosporidium* isolates to explore whether CSpV1 genome diversity varies over time, with geographical sampling location, *C. parvum* genetic diversity, or ruminant host species. A total of 123 fecal samples of cattle, sheep and goats were collected from 17 different French departments (57 diseased animal fecal samples and 66 healthy animal fecal samples). Subtyping analysis of the *C. parvum* isolates revealed the presence of two zoonotic subtype families Ila and Ild. Sequence analysis of CSpV1 revealed that all CSpV1 from France, regardless of the subtype of *C. parvum* (IIaA15G2R1, IIaA17G2R1 and IldA18G1R1) are more closely related to CSpV1 from Turkey, and cluster on a distinct branch from CSpV1 collected from *C. parvum* subtype IIaA15G2R1 from Asia and North America. We also found that samples collected on a given year or successive years in a given location are more likely to host the same subtype of *C. parvum* and the same CSpV1 strain. Yet, there is no distinct clustering of CSpV1 per French department or ruminants, probably due to trade, and transmission of *C. parvum* among host species. Our results point towards (i) a close association between CSpV1 movement and *C. parvum* movement, (ii) recent migrations of *C. parvum* among distantly located departments and (iii) incidental transmission of *C. parvum* between ruminants. All together, these results provide insightful information regarding CSpV1 evolution and suggest the virus might be used as an epidemiological tracer for *C. parvum*. Future studies need to investigate CSpV1's role in *C. parvum* virulence and on subtype ability to infect different species.

**Keywords** *Cryptosporidium parvum*, *cryspovirus*, *calves*, *lambs*, *kids*, *France*

Handling editor: Frank Katzer

*Correspondence:
Karim Tarik Adjou
karim.adjou@vet-alfort.fr

Nolwenn M. Dheilly
nolwenn.dheilly@anses.fr

Full list of author information is available at the end of the article
**Introduction**

Species of the genus *Cryptosporidium* are protozoan parasites specific of vertebrates that are responsible for mild to severe diarrhea in both human and animals [1, 2]. In ruminants, diarrhea is particularly severe in newborn and young. It causes dehydration, weakness, anorexia, weight-loss and a delayed growth. Adults can remain asymptomatic and serve as reservoir host for the parasite [3]. Thus, Cryptosporidiosis represents a significant economic loss for farmers. *Cryptosporidium* is also a public health issue because contamination of surface water can lead to human infections. *Cryptosporidium* is among the most important cause of diarrhea in developing countries, and remain a major cause of waterborne outbreaks and illness worldwide in both immunocompetent and immunocompromised people. In healthy immunocompetent people, cryptosporidiosis usually resolves itself within a couple of weeks [4–7]. However, epidemic outbreaks are still associated with fatalities, mostly in young children, and a chronic and life-threatening cryptosporidiosis can develop in immunocompromised individuals, such as those affected by HIV, transplant recipients or infants [7].

The genus *Cryptosporidium* is composed of 44 recognized species, and 120 genotypes remain of unknown species status [8]. They represent different public health significance depending on their host specificity and zoonotic potential. Humans are mainly infected by *C. parvum* and *C. hominis*. *C. parvum* is by far the *Cryptosporidium* species with the broadest host range and most important zoonotic potential. Nearly 20 *C. parvum* subtypes have been described using the 60 kDa glycoprotein gene (*gp60*) as genetic locus. The subtype family IIc appears to have adapted to humans whereas the subtype family IIa is mostly found in cattle and subtype family IID is mostly found in sheep and goats [9, 10]. The highly transmissible subtype IIaA15G2R1 is the most prevalent in most industrialized nations [10–12]. In France, this highly transmissible subtype has been reported in feces from cattle, goat and sheep [13–15]. However, *C. parvum* has a panmictic/epidemic population structure resulting in a high genetic diversity, with frequent disagreement between *gp60* results and other genetic markers due to frequent genetic recombination within this subtype [11]. Multilocus typing tools have allowed a better characterization of *C. parvum* subtype IIaA15G2R1, revealing a geographic segregation with country-specific clusters and a genetic distance that correlates with geographic distance [16, 17]. Different subpopulations have also been found in different host species due to reduced gene flow [18–20]. The global expansion of the IIaA15G2R1 zoonotic subtype presents a major challenge. In particular, there is a

need to develop novel, easy to use methods to track *C. parvum*'s genetic structure, host-adaptation, virulence and threat to public health.

The *C. parvum* virus 1, of the genus *Cryspovirus*, family *Partitiviridae*, was first reported in *C. parvum* cytoplasm in 1997 [21, 22]. The viral genome is bisegmented and dsRNA1 and dsRNA2 each encode a single protein. Closely related strains of the cryspovirus were found in distinct *C. parvum* strains infecting humans and calves, as well as in *C. hominis*, *C. felis* and *C. meleagridis*, with >92% amino-acid similarities [23, 24]. A single study has assessed the cryspovirus prevalence in parasite populations, in Japan [25]. It showed that CSpV1 has 100% apparent prevalence and appears to co-diversify with *C. parvum* so that viral sequences might be used for sensitive detection of *C. parvum* subtype IIaA15G2R1 and to determine the regional origin of the parasite [25–27].

The aims of the present study were the following: (1) to determine by PCR (*18S rRNA*, PCR-RFLP, *gp60*) the presence of *Cryptosporidium* in calves, lambs and kid goats in different regions of France, (2) to characterize the presence and prevalence of cryspovirus in these French *Cryptosporidium* isolates, (3) to genotype the virus in order to test whether CSpV1 genome diversity varies over time, with geographical sampling location, with *C. parvum* genetic diversity, or with ruminant host species.

**Materials and methods**

**Sample collection**

A total of 123 fecal samples were collected voluntarily by veterinarians or departmental veterinarian laboratories, between 2018 and 2022, and originated from 17 different French departments (Table 1). Among those, 57 samples were collected from diseased animals that had been naturally infected by *Cryptosporidium* and presented a diarrhea at the time of sampling. An additional 66 fecal samples were collected from apparently healthy animals. Samples were collected with sterile gloves, mailed at room temperature, and conserved at 4 °C until use.

**Cryptosporidium detection and enrichment**

*Cryptosporidium* spp. were initially detected microscopically with direct immunofluorescence assay (DFA) (MeriFluor® *Cryptosporidium/Giardia*, Meridian Bioscience Europe, Milano, Italy) as previously described [14]. *Cryptosporidium* oocysts were purified from samples with positive DFA using Dynabeads™ anti-*Cryptosporidium* kit per manufacturer protocol. The success of *Cryptosporidium* purification was controlled using DFA.
<table><caption>Table 1 Number and origin of samples used in this study.</caption><thead><tr><th rowspan="2">Departement of origin</th><th colspan="2">Sheep</th><th colspan="2">Cattle</th><th colspan="2">Goats</th></tr><tr><th>Diseased lambs</th><th>Healthy</th><th>Diseased calves</th><th>Healthy</th><th>Diseased Goat kids</th><th>Healthy</th></tr></thead><tbody><tr><td>Tarn (81)</td><td>15</td><td>–</td><td>–</td><td>–</td><td>4</td><td>–</td></tr><tr><td>Haute-Vienne(87)</td><td>7</td><td>–</td><td>–</td><td>–</td><td>–</td><td>–</td></tr><tr><td>Loiret (45)</td><td>1</td><td>9 adults</td><td>–</td><td>–</td><td>–</td><td>–</td></tr><tr><td>Loire (42)</td><td>–</td><td>8 lambs</td><td>5</td><td>–</td><td>–</td><td>–</td></tr><tr><td>Allier (3)</td><td>–</td><td>–</td><td>8</td><td>–</td><td>–</td><td>–</td></tr><tr><td>Puy-de-Dôme (63)</td><td>–</td><td>–</td><td>4</td><td>–</td><td>–</td><td>–</td></tr><tr><td>Moselle (57)</td><td>–</td><td>–</td><td>2</td><td>–</td><td>–</td><td>–</td></tr><tr><td>Nièvre (58)</td><td>–</td><td>2 adults</td><td>1</td><td>–</td><td>–</td><td>8 adults + 6 kids</td></tr><tr><td>Yonne (89)</td><td>–</td><td>6 lambs</td><td>–</td><td>–</td><td>–</td><td>–</td></tr><tr><td>Torcy (77)</td><td>–</td><td>16 adults</td><td>–</td><td>–</td><td>–</td><td>2 adults</td></tr><tr><td>Côte d'or (21)</td><td>–</td><td>4 adults</td><td>–</td><td>–</td><td>–</td><td>–</td></tr><tr><td>Cher (18)</td><td>–</td><td>1 adult</td><td>–</td><td>–</td><td>–</td><td>–</td></tr><tr><td>Aveyron (12)</td><td>–</td><td>1 adult</td><td>–</td><td>–</td><td>–</td><td>–</td></tr><tr><td>Aisne (2)</td><td>–</td><td>2 adults</td><td>–</td><td>–</td><td>–</td><td>–</td></tr><tr><td>Lot (64)</td><td>1</td><td>–</td><td>–</td><td>–</td><td>–</td><td>–</td></tr><tr><td>Pyrénées-Atlantiques (46)</td><td>6</td><td>–</td><td>–</td><td>–</td><td>–</td><td>–</td></tr><tr><td>Ardennes (08)</td><td>–</td><td>–</td><td>3</td><td>1</td><td>–</td><td>–</td></tr><tr><td>Total</td><td>30</td><td>49</td><td>23</td><td>1</td><td>4</td><td>16</td></tr></tbody></table>

### DNA extraction, *Cryptosporidium* species determination and *C. parvum* subtyping

Disruption of oocyst walls was achieved using ten freeze-thaw cycles as previously described [28]. Then, DNA extraction was conducted using the QIAamp DNA Stool Mini Kit (Qiagen, France), according to the manufacturer's instructions. The *Cryptosporidium* genus was determined using a nested PCR amplifying a 840 bp fragment of 18S rRNA gene [29]. PCR products were digested with SspI and MboII endonucleases and restriction fragment length polymorphism (RFLP) analysis was used to determine the *Cryptosporidium* species [30]. *C. parvum* samples were subtyped by nested PCR-sequence analysis of the partial 60 kDa glycoprotein locus (gp60) [31] and Sanger sequencing of both strands (Genoscreen) (Genbank numbers in Table 2). *C. parvum* subtypes were named using the recommended nomenclature system [10, 32].

### RNA extraction, and Cryspovirus sequencing

Disruption of oocyst walls was achieved using ten freeze-thaw cycles and a proteinase K treatment for 1 h at 55 °C. Total RNA were extracted with RNeasy mini kit (Qiagen) following the manufacturer's recommendations. Total RNA was used for cDNA synthesis with Maximus H minus reverse transcriptase (Thermofisher). PCR was carried out with primer sets amplifying the

RNA-dependent RNA polymerase (RdRP) (dsRNA1: CPVL_ORF_F 5'-AAGTTTGTCAATATCTATGAGATA C-3', CPVL_ORF_R 5'-TCCATAAATTTTGTGACT CCTG-3') and capsid (dsRNA2: CPVS_ORF_F 5'-ATT ACAAGTTTTGAATCAATAGAG-3', CPVS_ORF_R 5'-ATGGGAGCGATCTGCGCTACAC-3') genes as initially described by Murakoshi et al. [25]. The resulting 1468 bp fragments and 867 bp fragments represent 80% of dsRNA1 and 58% of dsRNA2 genome fragments. They were visualized using gel electrophoresis before Sanger sequencing (Genoscreen). Genbank accession numbers of successful sequencing are provided in Table 2. Additional PCR were carried out using a second set of primers amplifying smaller fragments of the RdRP (RdrpV-F: 5'-TGGGCATATGGTGCTCACTA-3'; RdrpV-R: 5'-GCTAAGAGAT CGT AGATGTCCA-3) and capsid (CPV-F: 5'-TGGTTCCGATTTTACCGGAA-3'; CPV-R: 5'-ACGACAATTAGGACTCAAATGACC-3') as described by Jenkins et al. [33]. CSpV1 was considered detected if at least one of the two PCR approaches returned positive results.

### Full-length genome sequencing

CSpV1 genome was obtained from a pool of samples B44 and B45, using RNA sequencing. cDNA libraries were prepared using an Ion total RNAseq Kit (Life technologies, Carlsbad, CA, USA) according to the supplier's instructions. The cDNA libraries were sequenced using
<table><caption>Table 2 CSpV1 detection and partial genome sequencing in C. parvum isolates.</caption><thead><tr><th>Spl</th><th>Year</th><th>Dpt</th><th>GP60 subtyping</th><th>Genbank acc #</th><th>RdRpV</th><th>VL</th><th>Genbank acc #</th><th>CPV</th><th>VS</th><th>Genbank acc #</th></tr></thead><tbody><tr><td>B36</td><td>2020</td><td>42</td><td>IIaA15G2R1</td><td>OQ722152</td><td>pos</td><td>pos</td><td>OQ749437</td><td>pos</td><td>pos</td><td>OQ749467</td></tr><tr><td>B38</td><td>2021</td><td>42</td><td>IIaA17G1R1</td><td>OQ722153</td><td>pos</td><td>pos</td><td>OQ749438</td><td>pos</td><td>pos</td><td>OQ749468</td></tr><tr><td>B40</td><td>2021</td><td>42</td><td>IIaA17G1R1</td><td>OQ722172</td><td>pos</td><td>pos</td><td>OQ749439</td><td>pos</td><td>pos</td><td>OQ749469</td></tr><tr><td>B41</td><td>2021</td><td>3</td><td>IIaA15G2R1</td><td>OQ722154</td><td>pos</td><td>pos</td><td>OQ749440</td><td>pos</td><td>pos</td><td>OQ749470</td></tr><tr><td>B42</td><td>2021</td><td>3</td><td>IIaA15G2R1</td><td>OQ722155</td><td>pos</td><td>pos</td><td>OQ749441</td><td>pos</td><td>pos</td><td>OQ749471</td></tr><tr><td>B43</td><td>2021</td><td>3</td><td>IIaA15G2R1</td><td>OQ722156</td><td>pos</td><td>pos</td><td>OQ749442</td><td>pos</td><td>pos</td><td>OQ749472</td></tr><tr><td>B44</td><td>2021</td><td>58</td><td>IIaA15G2R1</td><td>OQ722157</td><td>pos</td><td>pos</td><td>OQ749443</td><td>pos</td><td>pos</td><td>OQ749473</td></tr><tr><td>B45</td><td>2021</td><td>3</td><td>IIaA15G2R1</td><td>OQ722158</td><td>pos</td><td>pos</td><td>OQ749444</td><td>pos</td><td>pos</td><td>OQ749474</td></tr><tr><td>B46</td><td>2021</td><td>63</td><td>IIaA15G2R1</td><td>OQ722159</td><td>pos</td><td>pos</td><td>OQ749445</td><td>pos</td><td>pos</td><td>OQ749475</td></tr><tr><td>B47</td><td>2021</td><td>63</td><td>IIaA15G2R1</td><td>OQ722160</td><td>pos</td><td>pos</td><td>OQ749446</td><td>pos</td><td>pos</td><td>OQ749476</td></tr><tr><td>B48</td><td>2021</td><td>3</td><td>IIaA15G2R1</td><td>OQ722161</td><td>pos</td><td>pos</td><td>OQ749447</td><td>pos</td><td>pos</td><td>OQ749477</td></tr><tr><td>B49</td><td>2021</td><td>3</td><td>IIaA15G2R1</td><td>OQ722162</td><td>pos</td><td>pos</td><td>OQ749448</td><td>pos</td><td>pos</td><td>OQ749478</td></tr><tr><td>B50</td><td>2021</td><td>3</td><td>IIaA15G2R1</td><td>OQ722163</td><td>neg</td><td>neg</td><td></td><td>neg</td><td>neg</td><td></td></tr><tr><td>B51</td><td>2021</td><td>63</td><td>IIaA15G2R1</td><td>OQ722164</td><td>pos</td><td>pos</td><td>OQ749449</td><td>pos</td><td>pos</td><td>OQ749479</td></tr><tr><td>B52</td><td>2021</td><td>63</td><td>IIaA15G2R1</td><td>OQ722165</td><td>pos</td><td>pos</td><td>OQ749450</td><td>pos</td><td>pos</td><td>OQ749480</td></tr><tr><td>B53</td><td>2021</td><td>3</td><td>IIaA15G2R1</td><td>OQ722166</td><td>pos</td><td>pos</td><td></td><td>neg</td><td>pos</td><td>OQ749481</td></tr><tr><td>B54</td><td>2021</td><td>57</td><td>unknown</td><td></td><td>pos</td><td>pos</td><td>OQ749451</td><td>pos</td><td>pos</td><td>OQ749482</td></tr><tr><td>B55</td><td>2021</td><td>57</td><td>IIaA15G2R1</td><td>OQ722167</td><td>pos</td><td>pos</td><td>OQ749452</td><td>pos</td><td>pos</td><td>OQ749483</td></tr><tr><td>B65</td><td>2022</td><td>8</td><td>IIaA15G2R1</td><td>OQ722168</td><td>pos</td><td>pos</td><td>OQ749434</td><td>pos</td><td>pos</td><td></td></tr><tr><td>B66</td><td>2022</td><td>8</td><td>IIaA15G2R1</td><td>OQ722169</td><td>pos</td><td>pos</td><td>OQ749435</td><td>pos</td><td>pos</td><td></td></tr><tr><td>B67</td><td>2022</td><td>8</td><td>IIaA15G2R1</td><td>OQ722170</td><td>pos</td><td>pos</td><td>OQ749436</td><td>pos</td><td>pos</td><td></td></tr><tr><td>O1</td><td>2018</td><td>81</td><td>IIaA15G2R1</td><td>MN037849</td><td>pos</td><td>pos</td><td>OQ749453</td><td>pos</td><td>neg</td><td></td></tr><tr><td>O10</td><td>2019</td><td>81</td><td>IIaA16G3R1</td><td>MN037858</td><td>pos</td><td>pos</td><td></td><td>pos</td><td>pos</td><td>OQ749488</td></tr><tr><td>O11</td><td>2019</td><td>81</td><td>IIaA15G2R1</td><td>MN037859</td><td>pos</td><td>pos</td><td>OQ749456</td><td>pos</td><td>pos</td><td>OQ749489</td></tr><tr><td>O12</td><td>2019</td><td>81</td><td>IIaA15G2R1</td><td>MN037860</td><td>pos</td><td>neg</td><td></td><td>pos</td><td>pos</td><td>OQ749490</td></tr><tr><td>O13</td><td>2019</td><td>81</td><td>IIaA15G2R1</td><td>MN037861</td><td>pos</td><td>pos</td><td></td><td>neg</td><td>pos</td><td>OQ749491</td></tr><tr><td>O14</td><td>2019</td><td>81</td><td>IIaA15G2R1</td><td>MN037862</td><td>pos</td><td>pos</td><td></td><td>neg</td><td>pos</td><td>OQ749492</td></tr><tr><td>O15</td><td>2019</td><td>81</td><td>IIaA15G2R1</td><td>MN037863</td><td>neg</td><td>neg</td><td></td><td>neg</td><td>pos</td><td>OQ749493</td></tr><tr><td>O16</td><td>2019</td><td>81</td><td>IIaA16G3R1</td><td>MN037864</td><td>pos</td><td>pos</td><td>OQ749457</td><td>pos</td><td>pos</td><td>OQ749494</td></tr><tr><td>O17</td><td>2019</td><td>81</td><td>IIaA15G2R1</td><td>MN037865</td><td>neg</td><td>neg</td><td></td><td>neg</td><td>pos</td><td></td></tr><tr><td>O18</td><td>2019</td><td>81</td><td>IIaA15G2R1</td><td>MN037866</td><td>pos</td><td>pos</td><td>OQ749458</td><td>pos</td><td>pos</td><td>OQ749495</td></tr><tr><td>O19</td><td>2019</td><td>81</td><td>IIaA15G2R1</td><td>MN037867</td><td>pos</td><td>pos</td><td>OQ749459</td><td>pos</td><td>pos</td><td>OQ749496</td></tr><tr><td>O2</td><td>2018</td><td>81</td><td>IIaA15G2R1</td><td>MN037850</td><td>pos</td><td>neg</td><td></td><td>pos</td><td>pos</td><td></td></tr><tr><td>O20</td><td>2019</td><td>81</td><td>IIaA15G2R1</td><td>MN037868</td><td>pos</td><td>neg</td><td></td><td>pos</td><td>pos</td><td>OQ749497</td></tr><tr><td>O21</td><td>2019</td><td>81</td><td>IIaA15G2R1</td><td>MN037869</td><td>pos</td><td>pos</td><td>OQ749460</td><td>pos</td><td>pos</td><td>OQ749498</td></tr><tr><td>O22</td><td>2019</td><td>81</td><td>IIaA15G2R1</td><td>MN037870</td><td>neg</td><td>neg</td><td></td><td>neg</td><td>neg</td><td></td></tr><tr><td>O23</td><td>2020</td><td>45</td><td>IIaA17G2R1</td><td>MN037871</td><td>pos</td><td>pos</td><td>OQ749461</td><td>pos</td><td>neg</td><td></td></tr><tr><td>O3</td><td>2019</td><td>87</td><td>IIaA15G2R1</td><td>MN037851</td><td>pos</td><td>pos</td><td>OQ749454</td><td>pos</td><td>pos</td><td>OQ749484</td></tr><tr><td>O4</td><td>2019</td><td>87</td><td>IIaA15G2R1</td><td>MN037852</td><td>neg</td><td>neg</td><td></td><td>pos</td><td>neg</td><td></td></tr><tr><td>O5</td><td>2019</td><td>87</td><td>IIaA15G2R1</td><td>MN037853</td><td>pos</td><td>pos</td><td>OQ749455</td><td>pos</td><td>pos</td><td>OQ749485</td></tr><tr><td>O6</td><td>2019</td><td>87</td><td>IIaA15G2R1</td><td>MN037854</td><td>pos</td><td>neg</td><td></td><td>pos</td><td>pos</td><td>OQ749486</td></tr><tr><td>O7</td><td>2019</td><td>87</td><td>IIaA15G2R1</td><td>MN037855</td><td>neg</td><td>neg</td><td></td><td>pos</td><td>pos</td><td></td></tr><tr><td>O73</td><td>2021</td><td>64</td><td>IIaA15G2R1</td><td>OQ722148</td><td>pos</td><td>pos</td><td>OQ749462</td><td>pos</td><td>pos</td><td>OQ749499</td></tr><tr><td>O74</td><td>2021</td><td>46</td><td>IIaA18G1R1</td><td>OQ722149</td><td>pos</td><td>pos</td><td>OQ749463</td><td>pos</td><td>pos</td><td>OQ749500</td></tr><tr><td>O75</td><td>2021</td><td>46</td><td>IIaA18G1R1</td><td>OQ722150</td><td>pos</td><td>pos</td><td>OQ749464</td><td>pos</td><td>pos</td><td>OQ749501</td></tr><tr><td>O76</td><td>2021</td><td>46</td><td>IIaA18G1R1</td><td>OQ722151</td><td>pos</td><td>pos</td><td>OQ749465</td><td>pos</td><td>pos</td><td>OQ749502</td></tr><tr><td>O77</td><td>2021</td><td>46</td><td>IIaA18G1R1</td><td>OQ747365</td><td>pos</td><td>pos</td><td>OQ749466</td><td>pos</td><td>pos</td><td></td></tr><tr><td>O78</td><td>2021</td><td>46</td><td>IIaA18G1R1</td><td>OQ747366</td><td>neg</td><td>pos</td><td></td><td>pos</td><td>pos</td><td></td></tr><tr><td>O79</td><td>2021</td><td>46</td><td>IIaA15G2R1</td><td>OQ747367</td><td>pos</td><td>pos</td><td></td><td>pos</td><td>pos</td><td></td></tr></tbody></table>
<table><caption>**Table 2** (continued)</caption><thead><tr><th>Spl</th><th>Year</th><th>Dpt</th><th>GP60 subtyping</th><th>Genbank acc #</th><th>RdRpV</th><th>VL</th><th>Genbank acc #</th><th>CPV</th><th>VS</th><th>Genbank acc #</th></tr></thead><tbody><tr><td>O8</td><td>2019</td><td>87</td><td>IIaA15G2R1</td><td>MN037856</td><td>pos</td><td>neg</td><td></td><td>neg</td><td>pos</td><td></td></tr><tr><td>O9</td><td>2019</td><td>87</td><td>IIaA15G2R1</td><td>MN037857</td><td>pos</td><td>neg</td><td></td><td>pos</td><td>pos</td><td>OQ749487</td></tr><tr><td>C1</td><td>2019</td><td>81</td><td>IIaA15G2R1</td><td>MN037844</td><td>pos</td><td>pos</td><td></td><td>neg</td><td>pos</td><td>OQ749503</td></tr><tr><td>C2</td><td>2019</td><td>81</td><td>IIaA15G2R1</td><td>MN037845</td><td>pos</td><td>neg</td><td></td><td>neg</td><td>pos</td><td></td></tr><tr><td>C3</td><td>2019</td><td>81</td><td>IIaA15G2R1</td><td>MN037847</td><td>pos</td><td>neg</td><td></td><td>neg</td><td>pos</td><td>OQ749504</td></tr><tr><td>C4</td><td>2019</td><td>81</td><td>IIaA15G2R1</td><td>MN037848</td><td>neg</td><td>pos</td><td></td><td>neg</td><td>neg</td><td></td></tr></tbody></table>

an Ion Proton Sequencer and an Ion PI Chip v2 (Life technologies).

Sequence reads were cleaned and trimmed for adapter removal using fastP version 0.20.1 [34] and sequence quality was verified using FastQC version 0.11.8 [35]. Reads were assembled with rnaSPADES de novo assembler as implemented in SPAdes assembler version 3.10.0 [36]. Resulting contigs were aligned on local nt database with Megablast version 2.10.1 to identify viral references. Then, both sequence reads and assembled contigs were aligned using Burroughs-Wheeler Aligner (BWA, version 0.7.8) [37] against the CSpV1 strain Iowa genome fragments (NC_038843 dsRNA1 and NC_038844 dsRNA2) and visualized in Integrative Genome Viewer (IGV) [38, 39] to control the quality of the consensus sequence extracted using Samtools pileup version 1.8 [40].

## Phylogenetic analyses

The CSpV1 dsRNA1 and dsRNA2 sequences were aligned against all other CSpV1 sequences known to date. Nucleotide and protein percentage of identity were calculated with Clustal Omega [41]. Phylogenetic analysis were conducted on the full-length open reading frame. Phylogenetic tree were then inferred using the maximum likelihood method implemented in PhyML (version 3.0) [42] using the best-fit model and best of NNI and Subtree Pruning and Regrafting (SPR) branch swapping. Support for nodes on the trees were assessed using an approximate likelihood ratio test (aLRT) with the Shimodaira-Hasegawa-like procedure. Trees generated using the Neighbor-Joining and Maximum Parsimony methods gave identical results.

## Data availability

Sequencing data were submitted to Genbank under Bio-project ID PRJNA947786 using the MIUVIG symbiont-associated package to provide metadata on the host of the virus (*C. parvum*) and on the host of the host (cattle, sheep and goat) [43, 44].

## Results

None of the healthy goats, sheep and cattle were infected by *Cryptosporidium*. Among the 57 diseased animals presenting diarrhea, 55 were infected with *Cryptosporidium* and one individual was infected with *Giardia*. 18S PCR and sequencing revealed the presence of *C. parvum* only. This was in agreement with previous studies in France showing that *C. parvum* is responsible for most infections in pre-weaned calves, lambs and goats [45]. In France, the two most prevalent *C. parvum* belong to subtype families IIa and IId [14]. For one cattle sample with low parasite load, the 18S gene could not be sequenced. For all other samples, *C. parvum* subtyping revealed that cattle were infected by subtypes IIaA15G2R1 (18 ind; 90%) and IIaA17G1R1 (2 ind, 10%), sheep were infected by subtypes IIaA15G2R1 (22 ind; 73%), IIaA16G3R1 (2 ind; 7%), IIaA17G2R1 (1 ind, 3%) and IIdA18G1R1 (5 ind; 17%), and all 4 goat kids were infected by IIaA15G2R1 (Figure 1, Table 2). Thus, IIaA15G2R1 *C. parvum* subtype was the most prevalent, and was found in 11 departments on all four successive years (2018 to 2022). IIaA17G1R1 and IIaA17G2R1 were found in one department in 2020 and 2021, while IIaA16G3R1 and IIdA18G1R1 were found only once, in one department.

The full-length genome of CSpV1 from France was obtained from a pool of *C. parvum* of subtype IIaA15G2R1 collected from cattle in 2021 (Samples B44 and B45; Genbank acc numbers OQ686777 and OQ686778 for dsRNA1 and dsRNA2 respectively). The representative dsRNA1 and dsRNA2 FR strain genome was obtained with a coverage depth of 50 and 70 reads per position and with 0.6 and 0.8% iupac codes for dsRNA1 and dsRNA2, respectively. The genome sequence was compared to the three available full-length Cryspovirus genomes from Iowa (Iowa), Kansas (KSU) and China. The FR strain dsRNA1 shows a maximum of 97.3% nucleotide (nc) identity to the Iowa strain and a minimum of 96.4% nc identity to the KSU strain from North America. The RdRP protein shows 99% and 97% amino-acid (aa)

<footnote>* Spl: Sample. Letter (B, bovine / cattle; C, Caprine / goats and O, ovine / sheep); Dpt: Department; Cryspovirus PCR was carried out using primers RdRPV and CPV from Jenkins et al. [33] primers VL and VS from Murakoshi et al. [25]. GP60 *C. parvum* subtyping and VL/VS Cryspovirus amplicon were followed with Sanger sequencing. Genbank accession numbers are provided.</footnote>
<figure><img src="image_4.png" /><figcaption>**Figure 1** *C. parvum* subtyping in France between 2019 and 2021. Gp60 based subtyping of *C. parvum* was conducted on isolates from 11 different French departments collected from Cattle, Sheep and Goats. A pie chart is provided for each department (dpt) within which *C. parvum* was detected, and provide the *C. parvum* subtype and the host. The size of the pie charts is proportional to the number of samples collected.</figcaption></figure>

identity to the Iowa and KSU strains, respectively. The FR strain dsRNA2 shows a maximum of 98.4% nc identity to the Iowa strain and a minimum of 98% nc identity to KSU strain and was more distantly related to the related viruses discovered in *Cryptosporidium meleagridis* (85.6% nc identity), *Cryptosporidium felis* (86.7% nc identity), and *Cryptosporidium hominis* (90.8% nc identity). The capsid protein shows 98.6% and 99.1% aa identity to the Iowa and KSU strains, respectively.

Using two different sets of primers per virus genome fragment, we confirmed the presence of CSpV1 within all but two isolates of *C. parvum* (Table 2), including in the sample for which the *C. parvum* subtype could not be determined. This high prevalence was in accordance with previous reports [21, 25]. Good quality partial sequences of CSpV1 virus dsRNA1 encoding the RdRP, and dsRNA2 encoding the Capsid proteins were obtained for 30 and

42 samples, respectively (Genbank accession numbers are provided in Table 2). CSpV1 partial genome sequencing provided evidence of co-infection by two different genotypes of CSpV1 in *C. parvum* isolate O77, depicted by the presence of two overlapping peaks on chromatograms. Interestingly, point variations observed co-localize with nucleotide positions that show variations when comparing all sequences with each other, further supporting the co-infection hypothesis (Additional file 1). A greater genetic diversity was observed within dsRNA1 that shows an average of 95.6% nc identity to CSpV1 strain Iowa whereas dsRNA2 presents an average of 98% nc identity to CSpV1 strain Iowa. We conducted phylogenetic analyses on the nucleotide sequences of dsRNA1 and dsRNA2 and yielded similar results. As previously described, all sequences from Japan clustered closely together on a single branch. Similarly, sequences from
France clustered closely together and are most closely related to the only other CSpV1 partial virus sequence from the European continent (origin: Turkey) (Figures 2 and 3).

Discussion

The parasite population structure and rate of CSpV1 ver-
tical and horizontal transmission is expected to influence
virus evolution. Here, phylogenetic analyses conducted
on both genome fragments (dsRNA1 and dsRNA2) show

that CSpV1 from France clustered together with the
sequence of CSpV1 from Turkey, and on a branch dis-
tinct from CSpV1 collected from North America or Asia
(Figures 2 and 3). This result indicates that the CSpV1
dsRNA1 and dsRNA2 can distinguish sampling area at
the continental level. Furthermore, we provide herein the
first characterization of CSpV1 from different subtypes
of *C. parvum*. Indeed, studies conducted on CSpV1 from
North America and Japan focused on *C. parvum* subtype
IIaA15G2R1 [25, 46], while the subtype of *C. parvum*

<figure><img src="image_5.png" /><figcaption>**Figure 2** Tracking CSpV1 movements and evolution through phylogenetic analysis of dsRNA1. Phylogenetic tree of a 1468 nc long fragment of dsRNA1 of CSpV1 using all available CSpV1 sequences. The tree was inferred in PhyML using the LG substitution model. Branch points indicate that the results of Shimodaira-Hasgawa branch test > 0.8. Scale bar shows the number of nucleotide changes. Genbank virus isolate names are given as follows: accession number_**C. parvum** subtype_strain_name. French virus isolate names are provided as follows: Letter (B: bovine/cattle and O; ovine / sheep) followed with individual number, department number, C. parvum subtype, and sampling year. A map of France providing the localisation of departments (dpt) is provided next to the tree.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Figure 3** Tracking CSpV1 movements and evolution through phylogenetic analysis of dsRNA2. Phylogenetic tree of a 867 nt long fragment of dsRNA2 of CSpV1 using all available CSpV1 sequences. The tree was inferred in PhyML using the LG substitution model. Branch points indicate that the results of Shimodaira-Hasgawa branch test > 0.8. Scale bar shows the number of nucleotide changes. Genbank virus isolate names are given as follows: accession number_C. parvum subtype_strain name. French virus isolates names are provided as follows: Letter (B: bovine / cattle, C, Caprine / goats and O; ovine / sheep) followed with individual number, department number, C. parvum subtype, and sampling year. A map of France providing the localisation of departments (dpt) is provided next to the tree.</figcaption></figure>

from which the CSpV1 from Turkey and China were sequenced had not been characterized [47, 48]. The fact that regardless of the subtype of C. parvum (IIaA15G2R1, IIaA17G2R1 and IIdA18G1R1), all European sequences cluster together indicates that CSpV1 host shift between different parasite subtypes does occur. The rate of such host shift might be high considering that, focusing on sequencing from France only, we did not observe any clustering of CSpV1 depending on parasite subtypes. Experimental studies would be needed to estimate the rate of virus host shift and characterize the underlying

mechanism. Host shift can only occur within a host co-infected by two subtypes of C. parvum. Coinfections at the species level have frequently been reported and recent studies have demonstrated that different C. parvum subtypes recombine and can give rise to more highly virulent progeny [49, 50]. Indeed, recent reports indicate that classical methods based on Sanger sequencing of PCR amplified *gp60* alleles strongly underestimate the intra-host genotype diversity [51–53]. Moreover, single oocysts can harbor a mixed population of sporozoites, which provide opportunities for CSpV1 host shift
between *gp60* subtypes of *C. parvum* in the same animal. Host co-infection by different strains of *C. parvum* could also explain why we found some evidence of co-infection by different CSpV1 in isolate O77. We did not conduct clonal isolation of *C. parvum*, which means that if the host was co-infected by different *C. parvum* strains, they were present in the mixture that was analyzed. Therefore, the working hypothesis that each *C. parvum* strain carries only a single CSpV1 remains to be tested through more advanced analyses of CSpV1 genotype and *C. parvum* genotyping.

It has been suggested that CSpV1 evolution can be used to track movement of *C. parvum* [25]. Thus, we compared the phylogenetic position of CSpV1 dsRNA1 and dsRNA2 with the ruminant host, sampling year, and geographic locations. Our results indicate that samples collected on a given year and given location are more likely to host the same subtype of *C. parvum* and the same CSpV1 strain (Figures 1 and 2). In addition, when sampling occurred on successive years in the same location, we found the same subtype of *C. parvum* and the same CSpV1 strain (Figure 2). Yet, there is no distinct clustering of viruses per department or ruminant host: closely related CSpV1 were found in distant departments (i.e. viruses of IIaA15G2R1 from cattle of departments 63, 03, 57 and 58, Figures 1 and 2); also, closely related CSpV1 are often found in different host species (Figures 2 and 3). Our results point towards (i) a close association between CSpV1 movement and *C. parvum* movement, (ii) recent migrations of *C. parvum* among distantly located departments and (iii) incidental transmission of *C. parvum* between ruminants. Inter-species local transmission of *Cryptosporidium* spp. is characteristic of the parasite epidemiology [54–56]. Thus, the presence of *C. parvum* in many animal species, inter-species transmission, movements of animals due to trade among different regions of the country, and the potential role of wild animals in parasite dispersion [56], implies that the sources of *C. parvum* and CSpV1 movements are multiple.

In conclusion, we provide the first genomic data of CSpV1 in France. We also provide the first comparative analysis of the genome of CSpV1 from different *C. parvum* subtypes, from different ruminant hosts and over successive years. This analysis provides insightful information regarding both *C. parvum* and CSpV1 transmission and evolution. In agreement with Murakoshi et al. [25], our results suggest that CSpV1 is vertically transmitted in *C. parvum*, and that transmission between subtypes likely occurs as a result of recombination during the parasite's sexual life-cycle. While *gp60* gene coding is the most common marker used for epidemiological surveillance of *C. parvum*, the limitations associated with this approach are recognized

and novel techniques such as whole genome sequencing and multilocus genotyping are being investigated as alternative approaches [49, 57]. Herein, we identified CSpV1 in all but two *C. parvum* isolates, which can either suggest the absence of the virus within these two *C. parvum*-positive samples, or illustrate that our sensitivity was suboptimal. Given that CSpV1 was often identified using a set of primers and not others suggest that neither of our primers provided sufficient sensitivity. Moreover, we failed at obtaining good quality CSpV1 sequences from many samples, which limited the scope of our analyses. In the future, amplicon-based sequencing approaches are likely to provide better sensitivity, while allowing for the sequencing of multiple variants of CSpV1 within the same sample. CSpV1 genome sequencing may represent a low-cost alternative approach to track *C. parvum* epidemiology and will help identify co-infections, and characterize the rates of outcrossing and inter-species transmission. Coupling amplicon-based sequencing of CSpV1 with multilocus genotyping of *C. parvum* would most likely provide a more accurate description of the parasite and virus co-diversification. Whether the virus' ability to jump host between *C. parvum* belonging to different subclades can modulate the protozoan parasite's ability to infect different ruminants, its pathogenicity, and zoonotic potential remains to be investigated through experimental approaches.

## Supplementary Information

The online version contains supplementary material available at https://doi.org/10.1186/s13567-023-01196-4.

**Additional file 1: Coinfection of O77 by different variants of CSpV1.**

A Example of point mutation in capillary sequencing chromatogram confirming the presence of different variants of CSpV1 in sample O77. B Alignment of all sequences to the same region showing that the point mutation co-localize with a region where variation occur often between different CSpV1.

## Acknowledgements

This project is part of the Parasite Microbiome Project and was funded by the French Agency for Food, Environmental and occupational Health Safety project PARAVIR to NMD and GK. We would like to thank the veterinarians who participated in fecal sample collection.

## Authors' contributions

KTA, NMD and GK participated in the conception and design of the study. KTA coordinated sampling of the fecal samples. KTA, NMD, BP, MM and GK designed protocol for molecular assay and participated in the analysis of results. AC, HL, RA, MT carried out molecular assays. PL, YB and NMD carried out the viral sequencing and genome assembly. NMD conducted phylogenetic analyses. NMD drafted the manuscript, with help from KTA, and GK. All authors read and approved the final manuscript.

## Declarations

### Competing interests

The authors declare that they have no competing interests.
Author details

¹Laboratoire de Santé Animale, Ecole Nationale Vétérinaire d'Alfort, Anses, INRAE, UMR BIPAR, 94700 Maisons-Alfort, France. ²Anses Animal Health Laboratory, UMR1161 Virology, INRAE, Anses, ENVA, Maisons-Alfort, France.

³Laboratoire de Ploufragan-Plouzané-Niort, Unité Génétique virale et biosécurité, ANSES, Agence Nationale de Sécurité Sanitaire de l'Alimentation, de l'Environnement et du Travail, Ploufragan, France.

Received: 12 May 2023 Accepted: 23 June 2023

Published online: 22 August 2023

References

1. Meganck V, Hoflack G, Opsomer G (2014) Advances in prevention and therapy of neonatal dairy calf diarrhoea: a systematical review with emphasis on colostrum management and fluid therapy. Acta Vet Scand 56:75

2. Kotloff KL, Nataro JP, Blackwelder WC, Nasrin D, Farag TH, Panchalingam S, Wu Y, Sow SO, Sur D, Breiman RF, Faruque AS, Zaidi AK, Saha D, Alonso PL, Tamboura B, Sanogo D, Onwuchekwa U, Manna B, Ramamurthy T, Kanungo S, Ochieng JB, Omore R, Oundo JO, Hossain A, Das SK, Ahmed S, Qureshi S, Quadri F, Adegbola RA, Antonio M et al (2013) Burden and aetiology of diarrhoeal disease in infants and young children in developing countries (the Global Enteric Multicenter Study, GEMS): a prospective, case-control study. Lancet 382:209–222

3. Thompson RC, Palmer CS, O'Handley R (2008) The public health and clinical significance of Giardia and *Cryptosporidium* in domestic animals. Vet J 177:18–25

4. Mac Kenzie WR, Hoxie NJ, Proctor ME, Gradus MS, Blair KA, Peterson DE, Kazmierczak JJ, Addiss DG, Fox KR, Rose JB, Davis JP (1994) A massive outbreak in milwaukee of *Cryptosporidium* infection transmitted through the public water supply. N Engl J Med 331:161–167

5. Smith A, Reacher M, Smerdon W, Adak GK, Nichols G, Chalmers RM (2006) Outbreaks of waterborne infectious intestinal disease in England and Wales, 1992–2003. Epidemiol Infect 134:1141–1149

6. Beaudeau P, de Valk H, Vaillant V, Mannschott C, Tillier C, Mouly D, Ledrans M (2008) Lessons learned from ten investigations of waterborne gastroenteritis outbreaks, France, 1998–2006. J Water Health 6:491–503

7. Tzipori S, Widmer G (2008) A hundred-year retrospective on cryptosporidiosis. Trends Parasitol 24:184–189

8. Feng Y, Ryan UM, Xiao L (2018) Genetic diversity and population structure of *Cryptosporidium*. *Trends Parasitol* 34:997–1011

9. Xiao L, Feng Y (2017) Molecular epidemiologic tools for waterborne pathogens *Cryptosporidium* spp. and *Giardia duodenalis*. Food Water-borne Parasitol 8:9–14–32

10. Xiao L (2010) Molecular epidemiology of cryptosporidiosis: an update. Exp Parasitol 124:80–89

11. Feng Y, Torres E, Li N, Wang L, Bowman D, Xiao L (2013) Population genetic characterization of dominant *Cryptosporidium parvum* subtype IIaA15G2R1. Int J Parasitol 43:1141–1147

12. McKerr C, Adak GK, Nichols G, Gorton R, Chalmers RM, Kafatos G, Cosford P, Charlett A, Reacher M, Pollock KG, Alexander CL, Morton S (2015) An outbreak of *Cryptosporidium parvum* across England & Scotland associated with consumption of fresh pre-cut salad leaves, May 2012. PLoS One 10:e0125955

13. Mammeri M, Cartou L, Chevillot A, Thomas M, Julien C, Vallée I, Polack B, Follet J, Adjou KT (2019) First identification of *Cryptosporidium parvum* zoonotic subtype IIaA15G2R1 in diarrheal lambs in France. Vet Parasitol Reg Stud Reports 18:100355

14. Mammeri M, Chevillot A, Chenafi I, Thomas M, Julien C, Vallée I, Polack B, Follet J, Adjou KT (2019) Molecular characterization of Cryptosporidium isolates from diarrheal dairy calves in France. Vet Parasitol Reg Stud Reports 18:100323

15. Castro-Hermida JA, Pors I, Poupin B, Ares-Mazás E, Chartier C (2005) Prevalence of *Giardia duodenalis* and *Cryptosporidium parvum* in goat kids in western France. Small Rumin Res 56:259–264

16. Tanriverdi S, Grinberg A, Chalmers RM, Hunter PR, Petrovic Z, Akiyoshi DE, London E, Zhang L, Tzipori S, Tumwine JK, Widmer G (2008) Inferences about the global population structures of *Cryptosporidium parvum* and *Cryptosporidium hominis*. Appl Environ Microbiol 74:7227–7234

17. Cacciò SM, de Waele V, Widmer G (2015) Geographical segregation of *Cryptosporidium parvum* multilocus genotypes in Europe. Infect Genet Evol 31:245–249

18. Quilez J, Vergara-Castiblanco C, Monteagudo L, Del Cacho E, Sánchez-Acedo C (2011) Multilocus fragment typing and genetic structure of *Cryptosporidium parvum* Isolates from diarrheic preweaned calves in Spain. Appl Environ Microbiol 77:7779–7786

19. Quilez J, Vergara-Castiblanco C, Monteagudo L, Cacho ED, Sánchez-Acedo C (2013) Host association of *Cryptosporidium parvum* populations infecting domestic ruminants in Spain. Appl Environ Microbiol 79:5363–5371

20. Drumo R, Widmer G, Morrison LJ, Tait A, Grelloni V, D'Avino N, Pozio E, Cacciò SM (2012) Evidence of host-associated populations of *Cryptosporidium parvum* in Italy. Appl Environ Microbiol 78:3523–3529

21. Khramtsov NV, Woods KM, Nesterenko MV, Dykstra CC, Upton SJ (1997) Virus-like, double-stranded RNAs in the parasitic protozoan *Cryptosporidium parvum*. Mol Microbiol 26:289–300

22. Nibert ML, Woods KM, Upton SJ, Ghabrial SA (2009) Cryspovirus: a new genus of protozoan viruses in the family Partitiviridae. Arch Virol 154:1959–1965

23. Leoni F, Gallimore CI, Green J, McLauchlin J (2006) Characterisation of small double stranded RNA molecule in *Cryptosporidium hominis*, *Cryptosporidium felis* and *Cryptosporidium meleagridis*. Parasitol Int 55:299–306

24. Khramtsov NV, Chung PA, Dykstra CC, Griffiths JK, Morgan UM, Arrowood MJ, Upton SJ (2000) Presence of double-stranded RNAs in human and calf isolates of *Cryptosporidium parvum*. J Parasitol 86:275–282

25. Murakoshi F, Ichikawa-Seki M, Aita J, Yaita S, Kinami A, Fujimoto K, Nishikawa Y, Murakami S, Horimoto T, Kato K (2016) Molecular epidemiological analyses of *Cryptosporidium parvum* virus 1 (CSpV1), a symbiotic virus of *Cryptosporidium parvum*, in Japan. Virus Res 211:69–72

26. de Souza MS, O'Brien C, Santin M, Jenkins M (2019) A highly sensitive method for detecting *Cryptosporidium parvum* oocysts recovered from source and finished water using RT-PCR directed to Cryspovirus RNA. J Microb Meth 156:77–80

27. Tai L, Li J, Yin J, Zhang N, Yang J, Li H, Yang Z, Gong P, Zhang X (2018) A novel detection method of *Cryptosporidium parvum* infection in cattle based on Cryptosporidium parvum virus 1. Acta Biochim Biophys Sin 51:104–111

28. Sahraoui L, Thomas M, Chevillot A, Mammeri M, Polack B, Vallée I, Follet J, Ain-Baaziz H, Adjou KT (2019) Molecular characterization of zoonotic *Cryptosporidium* spp. and *Giardia duodenalis* pathogens in Algerian sheep. Vet Parasitol Reg Stud Rep 16:100280

29. Xiao L, Escalante L, Yang C, Sulaiman I, Escalante AA, Montali RJ, Fayer R, Lal AA (1999) Phylogenetic analysis of Cryptosporidium parasites based on the small-subunit rRNA gene locus. Appl Environ Microbiol 65:1578–1583

30. Feng Y, Ortega Y, He G, Das P, Xu M, Zhang X, Fayer R, Gatei W, Cama V, Xiao L (2007) Wide geographic distribution of Cryptosporidium bovis and the deer-like genotype in bovines. Vet Parasitol 144:1–9

31. Alves M, Xiao L, Sulaiman I, Lal AA, Matos O, Antunes F (2003) Subgenotype analysis of *Cryptosporidium* isolates from humans, cattle, and zoo ruminants in Portugal. J Clin Microbiol 41:2744–2747

32. Sulaiman IM, Hira PR, Zhou L, Al-Ali FM, Al-Shelahi FA, Shweiki HM, Iqbal J, Khalid N, Xiao L (2005) Unique endemicity of cryptosporidiosis in children in Kuwait. J Clin Microbiol 43:2805–2809

33. Jenkins M, O'Brien C, Fetterer R, Santin M (2016) RT-PCR specific for Cryspovirus is a highly sensitive method for detecting *Cryptosporidium parvum* oocysts. Food Waterborne Parasitol 5:14–20

34. Chen S, Zhou Y, Chen Y, Gu J (2018) fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34:i884–i890

35. Andrews S (2010) FASTQC. A quality control tool for high throughput sequence data. Babraham Bioinformatics, Babraham Institute, Cambridge, United Kingdom.

36. Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, Lesin VM, Nikolenko SI, Pham S, Prjibelski AD, Pyshkin AV, Sirokin AV, Vyahhi N, Tesler G, Alekseyev MA, Pevzner PA (2012) SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol 19:455–477

37. Li H, Durbin R (2009) Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25:1754–1760
38. Robinson JT, Thorvaldsdóttir H, Winckler W, Guttman M, Lander ES, Getz G, Mesirov JP (2011) Integrative genomics viewer. Nat Biotechnol 29:24–26

39. Thorvaldsdóttir H, Robinson JT, Mesirov JP (2012) Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform 14:178–192

40. Danecek P, Bonfield JK, Liddle J, Marshall J, Ohan V, Pollard MO, Whittwham A, Keane T, McCarthy SA, Davies RM, Li H (2021) Twelve years of SAMtools and BCFtools. GigaScience. https://doi.org/10.1093/gigascience/giab008

41. Madeira F, Park YM, Lee J, Buso N, Gur T, Madhusoodanan N, Basutkar P, Tivey ARN, Potter SC, Finn RD, Lopez R (2019) The EMBL-EBI search and sequence analysis tools APIs in 2019. Nucleic Acids Res 47:W636–W641

42. Lefort V, Longueville J-E, Gascuel O (2017) SMS: smart model selection in PhyML. Mol Biol Evol 34:2422–2424

43. Jorge F, Brealey JC, Brindley PJ, Buysse M, Cantacessi C, Duron O, Fichorova R, Fitzpatrick CR, Hahn M, Hunter C, Hervé V, Knoll LJ, Kohl KD, Lalle M, Lukeš J, Martínez JM, Perkins SL, Poulin R, Rosario K, Schneider AC, Schriml LM, Thompson LR, Walls RL, Dheilly NM (2022) MlXS-SA: a MlXS extension defining the minimum information standard for sequence data from symbiont-associated micro-organisms. ISME Commun 2:9

44. Roux S, Adriaenssens EM, Dutilh BE, Koonin EV, Kropinski AM, Krupovic M, Kuhn JH, Lavigne R, Brister JR, Varsani A, Amid C, Aziz RK, Bordenstein SR, Bork P, Breitbart M, Cochrane GR, Daly RA, Desnués C, Duhaime MB, Emerson JB, Enault F, Fuhrman JA, Hingamp P, Hugenholtz P, Hurwitz BL, Ivanova NN, Labonté JM, Lee K-B, Malmstrom RR, Martinez-Garcia M et al (2019) Minimum Information about an Uncultivated Virus Genome (MIUVIG). Nat Biotechnol 37:29–37

45. Santín M, Trout JM, Fayer R (2008) A longitudinal study of cryptosporidiosis in dairy cattle from birth to 2 years of age. Vet Parasitol 155:15–23

46. Vong M, Ludington JG, Ward HD, Nibert ML (2017) Complete cryspovirus genome sequences from *Cryptosporidium parvum* isolate Iowa. Arch Virol 162:2875–2879

47. Berber E, Şimşek E, Çanakoğlu N, Sürsal N, Gençay Göksu A (2021) Newly identified *Cryptosporidium parvum* virus-1 from newborn calf diarrhoea in Turkey. Transbound Emerg Dis 68:2571–2580

48. Li W, Zhang N, Liang X, Li J, Gong P, Yu X, Ma G, Ryan UM, Zhang X (2009) Transient transfection of *Cryptosporidium parvum* using green fluorescent protein (GFP) as a marker. Mol Biochem Parasitol 168:143–148

49. Tanriverdi S, Blain JC, Deng B, Ferdig MT, Widmer G (2007) Genetic crosses in the apicomplexan parasite *Cryptosporidium parvum* define recombination parameters. Mol Microbiol 63:1432–1439

50. Feng X, Rich SM, Tzipori S, Widmer G (2002) Experimental evidence for genetic recombination in the opportunistic pathogen *Cryptosporidium parvum*. Mol Biochem Parasitol 119:55–62

51. Grinberg A, Biggs PJ, Dukkipati VSR, George TT (2013) Extensive intra-host genetic diversity uncovered in *Cryptosporidium parvum* using Next Generation Sequencing. Infect Genet Evol 15:18–24

52. Dettwiler I, Troell K, Robinson G, Chalmers RM, Basso W, Rentería-Solís ZM, Daugschies A, Mühlethaler K, Dale MI, Basapathi Raghavendra J, Ruf M-T, Poppert S, Meylan M, Olias P (2022) TIDE Analysis of Cryptosporidium infections by gp60 typing reveals obscured mixed infections. J Infect Dis 225:686–695

53. Zahedi A, Gofton AW, Jian F, Paparini A, Oskam C, Ball A, Robertson I, Ryan U (2017) Next Generation Sequencing uncovers within-host differences in the genetic diversity of Cryptosporidium gp60 subtypes. Int J Parasitol 47:601–607

54. Thomson S, Hamilton CA, Hope JC, Katzer F, Mabbott NA, Morrison LJ, Innes EA (2017) Bovine cryptosporidiosis: impact, host-parasite interaction and control strategies. Vet Res 48:42

55. Krumkamp R, Conraths FJ, Caccio S, Scharés G, Hogan B, Winter D, Jaeger A, Melhem S, Rakotozandrindrainy N, May J, Rakotozandrindrainy R, Eibach D (2022) Clustering of *Cryptosporidium* species infections among sheep and cattle but not children in remote highland communities of Madagascar. Parasit Vectors 15:304

56. Hatam-Nahavandi K, Ahmadpour E, Carmena D, Spotin A, Bangoura B, Xiao L (2019) Cryptosporidium infections in terrestrial ungulates with focus on livestock: a systematic review and meta-analysis. Parasit Vectors 12:453

57. Robinson G, Pérez-Cordón G, Hamilton C, Katzer F, Connelly L, Alexander CL, Chalmers RM (2022) Validation of a multilocus genotyping scheme for

subtyping *Cryptosporidium parvum* for epidemiological purposes. Food Waterborne Parasitol 27:e00151

**Publisher's Note**

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Ready to submit your research? Choose BMC and benefit from:

* fast, convenient online submission

* thorough peer review by experienced researchers in your field

* rapid publication on acceptance

* support for research data, including large and complex data types

* gold Open Access which fosters wider collaboration and increased citations

* maximum visibility for your research: over 100M website views per year

At BMC, research is always in progress.

Learn more biomedcentral.com/submissions

<figure><img src="image_7.png" /></figure>