Detection and Molecular Diversity of Cryptosporidium spp. and Giardia duodenalis in the Endangered Iberian Lynx (Lynx pardinus), Spain

Simple Summary The Iberian lynx is an iconic feline species endemic to the Iberian Peninsula. Since the second half of the past century, its global population has decreased dramatically to the brink of extinction as a consequence of human-driven activities (habitat reduction and transformation, illegal hunting, road kills, density decrease in natural preys) and infectious diseases. Fortunately, the successful implementation of conservation programs has reversed this gloomy trend, allowing for an increase in the Iberian lynx population to over 1600 free-ranging animals in 2022. Regarding infectious diseases, very little is known on the epidemiology and health impact of the diarrhoea-causing intestinal protozoan parasites Cryptosporidium and Giardia in the Iberian lynx. To tackle these questions, we investigated the presence and molecular diversity of both pathogens in 256 collected faecal samples from 251 free-ranging and captive Iberian lynxes in Spain during the period 2017–2023. Our results demonstrate that Cryptosporidium (2.4%) and Giardia (27.9%) are present at different frequencies in the surveyed individuals. Our molecular analyses also indicate that a significant proportion of the Cryptosporidium infections detected are caused by strains that are typically found in the preys the Iberian lynx feed on. Interestingly, we also found that the Iberian lynx can harbour genetic variants of Cryptosporidium and Giardia with the potential to infect humans, although the likelihood of such events is judged low due to the light infections detected in the investigated animals. Abstract Cryptosporidium spp. and Giardia duodenalis are the main non-viral causes of diarrhoea in humans and domestic animals globally. Comparatively, much less information is currently available in free-ranging carnivore species in general and in the endangered Iberian lynx (Lynx pardinus) in particular. Cryptosporidium spp. and G. duodenalis were investigated with molecular (PCR and Sanger sequencing) methods in individual faecal DNA samples of free-ranging and captive Iberian lynxes from the main population nuclei in Spain. Overall, Cryptosporidium spp. and G. duodenalis were detected in 2.4% (6/251) and 27.9% (70/251) of the animals examined, respectively. Positive animals to at least one of them were detected in each of the analysed population nuclei. The analysis of partial ssu rRNA gene sequences revealed the presence of rodent-adapted C. alticolis (n = 1) and C. occultus (n = 1), leporid-adapted C. cuniculus (n = 2), and zoonotic C. parvum (n = 2) within Cryptosporidium, and zoonotic assemblages A (n = 5) and B (n = 3) within G. duodenalis. Subgenotyping analyses allowed for the identification of genotype VaA19 in C. cuniculus (gp60 locus) and sub-assemblages AI and BIII/BIV in G. duodenalis (gdh, bg, and tpi loci). This study represents the first molecular description of Cryptosporidium spp. and G. duodenalis in the Iberian lynx in Spain. The presence of rodent/leporid-adapted Cryptosporidium species in the surveyed animals suggests spurious infections associated to the Iberian lynx’s diet. The Iberian lynx seems a suitable host for zoonotic genetic variants of Cryptosporidium (C. parvum) and G. duodenalis (assemblages A and B), although the potential risk of human transmission is regarded as limited due to light parasite burdens and suspected low excretion of infective (oo)cysts to the environment by infected animals. More research should be conducted to ascertain the true impact of these protozoan parasites in the health status of the endangered Iberian lynx.


Introduction
Cryptosporidium spp.and Giardia duodenalis are major causative agents of diarrheal diseases in humans and a wide diversity of animals with a worldwide distribution [1,2].Human cryptosporidiosis is the leading protozoan cause of diarrheal mortality worldwide [3].In contrast, human giardiasis is rarely mortal but is associated with malabsorptive diarrhoea and impaired childhood growth [4,5].Both cryptosporidiosis and giardiasis also cause diarrhoea in neonatal ruminants, leading to high morbidity and mortality rates in the first three weeks [6][7][8][9] and significant economic losses for farmers [10,11].Cryptosporidium and Giardia infections are typically asymptomatic in free-living animals, raising concerns about their true health impact in wildlife and the role of wildlife in the epidemiology of these parasites [6,12].
The Iberian lynx (Lynx pardinus) is an emblematical felid species endemic to the Iberian Peninsula.It is listed as "endangered" by the International Union for Conservation of Nature's Red List of Threatened Species [42].Since the second half of the twentieth century, a sharp decrease in the number of Iberian lynxes brought the species to the brink of extinction due to habitat loss/transformation, illegal hunting, road kills, reduction in the density of its primary prey, the European rabbit (Oryctolagus cuniculus), and infectious diseases [43,44].Among the latter, clinical cases and mortality reported during the last two decades have been associated to bacterial (e.g., Mycobacterium bovis, Streptococcus canis) [45,46], viral (e.g., feline leukaemia virus, feline herpes virus, feline calicivirus, pseudorabies virus) [47,48], and parasitic (e.g., Neospora caninum, Toxoplasma gondii, Cystoisospora spp.) [49][50][51][52] pathogens.Although the development of conservation programs has reversed the trend, allowing for an increase in the Iberian lynx population to over 1600 free-ranging animals in 2022 [53], the monitoring of pathogens that could affect captive and free-ranging animals is still a key component of ongoing conservation programs [54,55].Following this line of action, this study aims to investigate the occurrence, genetic diversity, and zoonotic potential of the diarrhoea-causing enteric protozoan Cryptosporidium spp.and G. duodenalis in the Iberian lynx, a host species for which this information is currently lacking.

Study Area and Sampling
Faecal samples (n = 251) from Iberian lynxes were collected between 2017 and 2023.These included a total of 223 free-ranging animals from the three major population nuclei of this species in Spain (central, n = 63; south, n = 125; southwest, n = 33; unknown, 2), whereas 20 were lynxes maintained in captivity, including 14 animals from three captive breeding centres (BC1-BC3) belonging to the Iberian lynx ex situ conservation program and six from four zoo/conservation centres (ZC1-ZC4).The breeding and zoo/conservation centres were located in southern (n = 9) and southwestern (n = 10) Spain, respectively (Figure 1).Status information was not available for eight animals.In addition, five (three free-living, two captive) animals were longitudinally sampled during the study period.All faecal samples were taken from biological banks or animals subjected to medical check-ups, health programs, or surgical interventions during the study period.Faecal samples were obtained from the ground or the intestinal content of examined animals.Epidemiological information, including habitat status (free-living vs. captivity), sampling date, age (yearlings: <1 year old; subadults: 1 to 3 years old; adults: 3 to 10 years old; senile: >10 years old), sex, and sampling georeferenced location, was collected from each animal, whenever possible.All faecal samples studied were formed.This survey expands and complements those previously conducted on the very same Iberian lynx population that investigated the presence of other intestinal protists, including Microsporidia [56]

DNA Extraction and Purification of Faecal and Tissue Samples
Genomic DNA was isolated from approximately 100 mg of each faecal sample by using the IndiSpin Pathogen Kit (Indical Bioscience, Leipzig, Germany) according to the manufacturer's instructions.Extracted and purified DNA samples were eluted in 90 µL of PCR-grade water and kept at 4 °C until further molecular analysis.

Molecular Detection and Characterisation of Cryptosporidium spp.
Cryptosporidium spp.presence was investigated using a nested PCR protocol, amplifying a 587 bp fragment of the small subunit of the rRNA (ssu RNA) gene of the parasite [57].A subtyping tool based on the amplification of partial sequences of the 60 kDa gly-

DNA Extraction and Purification of Faecal and Tissue Samples
Genomic DNA was isolated from approximately 100 mg of each faecal sample by using the IndiSpin Pathogen Kit (Indical Bioscience, Leipzig, Germany) according to the manufacturer's instructions.Extracted and purified DNA samples were eluted in 90 µL of PCR-grade water and kept at 4 • C until further molecular analysis.

Molecular Detection and Characterisation of Cryptosporidium spp.
Cryptosporidium spp.presence was investigated using a nested PCR protocol, amplifying a 587 bp fragment of the small subunit of the rRNA (ssu RNA) gene of the parasite [57].A subtyping tool based on the amplification of partial sequences of the 60 kDa glycoprotein (gp60) [58] gene was used to ascertain intra-species genetic diversity in the samples that tested positive for C. parvum and C. cuniculus with ssu-PCR.

Molecular Detection and Characterisation of Giardia duodenalis
For the identification of G. duodenalis, a real-time PCR (qPCR) method was set-up to amplify a 62 bp fragment of the ssu RNA gene of the parasite [59].Samples that yielded cycle threshold (C T ) values < 35 in qPCR were then analysed through a nested PCR, used to amplify a 300 bp fragment of the ssu RNA gene [60,61] to assess G. duodenalis molecular diversity at the assemblage level.Samples that yielded qPCR C T values < 32 were additionally assessed using a sequence-based multilocus genotyping (MLST) scheme targeting the genes encoding for the glutamate dehydrogenase (gdh), β-giardin (bg), and triose phosphate isomerase (tpi) proteins to assess G. duodenalis molecular diversity at the sub-assemblage level.A 432 bp fragment of the gdh gene was amplified using a seminested PCR [62], while 511 and 530 bp fragments of the bg and tpi genes, respectively, were amplified through nested PCRs [63,64].

General Procedures
Detailed information on the PCR cycling conditions and oligonucleotides used for molecular identification and/or characterisation of the abovementioned parasites can be found in Tables S1 and S2, respectively.The previously described PCR protocols were conducted on a 2720 Thermal Cycler (Applied Biosystems, Foster City, CA, USA).The reaction mixes included 2.5 units of MyTAQ TM DNA polymerase (Bioline GmbH, Luckenwalde, Germany) and 5-10 µL 5× MyTAQ TM Reaction Buffer containing five mM deoxynucleotide triphosphates and 15 mM MgCl 2 .Negative and positive controls were included in all PCR runs.The PCR amplicons obtained were examined on a 1.5% D5 agarose gel stained with Pronasafe (Conda, Madrid, Spain) and sized using a 100 bp DNA ladder (Boehringer Mannheim GmbH, Mannheim, Germany).

Sequence and Phylogenetic Analysis
All amplicons of the expected size were directly sequenced in both directions with the internal primer pair in 10 µL reactions using Big Dye TM chemistries and an ABI 3730xl sequencer analyser (Applied Biosystems).The raw sequencing data were examined with Chromas Lite version 2.1 software ( http://chromaslite.software.informer.com/2.1,accessed on 18 January 2023) to generate consensus sequences.These sequences were compared with reference sequences deposited at the National Center for Biotechnology Information (NCBI) using the BLAST tool ( http://blast.ncbi.nlm.nih.gov/Blast.cgi,accessed on 18 January 2023).
To analyse the phylogenetic relationship among Cryptosporidium species and genotypes at the ssu rRNA locus, a maximum-likelihood tree was constructed using MEGA version 10 [65], based on substitution rates calculated with the general time reversible model and gamma distribution with invariant sites (G+I).Bootstrapping with 1000 replicates was used to determine support for the clades.The representative nucleotide sequences obtained in the present study were deposited in the GenBank public repository database under accession numbers OR916202-OR916206 and OR921171 (Cryptosporidium spp.) and OR916207-OR916209 and OR921172-OR921177 (G.duodenalis).

Statistics Analysis
Prevalence rates were estimated by dividing the number of positive animals by the total number of animals tested using two-sided exact binomial 95% confidence intervals (95% CI).Pearson's chi-squared test or Fisher's exact test was used to assess differences in the Cryptosporidium spp.and G. duodenalis infection rates, according to habitat, sex, age, sampling areas, and sampling period (categorised by terciles), using the R Statistical Package version 2.15.3 [66].A p-value < 0.05 was considered as statistically significant.

Results
The full dataset of this study, showing sampling, epidemiological, diagnostic, and molecular data, can be found in Table S3.

Occurrence of Cryptosporidium spp. and Giardia duodenalis
Table 3 summarises the occurrence of Cryptosporidium spp.and G. duodenalis in the Iberian lynx population (n = 251) under investigation according to the main epidemiological variables considered in this study.All faecal samples analysed (n = 256) had a formed consistency, suggestive of an apparent absence of gastrointestinal manifestations.a Missing values (number of samples with unknown data).b yearlings: <1 year old; sub-adults: 1 to 3 years old; adults: 3 to 10 years old; senile: >10 years old.
Cryptosporidium spp.DNA was detected in 2.4% (6/251; 95% CI: 0.9-5.1) of the individuals tested.On the other hand, G. duodenalis DNA was detected in 27.9% (70/251; 95% CI: 22.4-33.9) of the individuals tested.Giardia infections were observed in animals of all age groups, whereas no Cryptosporidium infections were detected in senile individuals.Three Iberian lynxes (two free-living, one captive) were co-infected with Cryptosporidium spp.and G. duodenalis.
None of the epidemiological variables considered were significantly associated with a higher likelihood of Giardia or Cryptosporidium infection except the sampling period for the latter (p = 0.042).The highest prevalence was detected in individuals sampled during the 2017-2020 period (6.8%), followed by 2022-2023 (1.8%) and 2021 (0.0%).Both Cryptosporidium spp.and G. duodenalis were detected in the three free-ranging areas sampled with frequencies varying from 1.5 to 7.3% and 23.9 to 33.3%, respectively.

Molecular Characterisation of Cryptosporidium spp.
The sequence analyses of the ssu rRNA region revealed the presence of four distinct Cryptosporidium species (C.alticolis, C. cuniculus, C. occultus, and C. parvum) in the Iberian lynx populations under study (Table 4).Cryptosporidium alticolis was identified in a freeliving animal from south Spain.The sequences generated at the ssu rRNA locus differed by five single nucleotide polymorphisms (SNPs, including three indels) from reference sequence MH145330 originally isolated from a common vole in the Czech Republic.Cryptosporidium cuniculus was identified in a free-living and a captive Iberian lynx, both in south Spain (Table 4 and Figure 1).Ssu rRNA sequences had 100% identity with reference sequence AY120901.One of the two isolates was successfully genotyped at the gp60 locus, revealing the presence of genotype VaA19.Cryptosporidium occultus was identified in a free-living animal in southwest Spain.Two additional isolates were assigned to C. parvum at the ssu rRNA marker: one belonged to a captive Iberian lynx in central Spain and the other to a free-living animal in southwest Spain.Both ssu rRNA sequences differed by 6-7 SNPs from the reference sequence AF112571 (Table 4).These include a hallmark deletion of 3-4 nucleotides at positions 686 to 689 of AF112571.Attempts to amplify these sequences at the gp60 marker failed.Phylogenetic analysis of ssu rRNA sequences revealed that all sequences generated in the present study belonging to C. alticolis, C. cuniculus, and C. parvum grouped together with appropriate reference sequences in well-defined clusters (Figure 2).

Molecular Characterisation of Giardia duodenalis
Giardia-positive samples with qPCR yielded C T values ranging from 20.0 to 39.7 (median: 34.5; standard deviation: 3.5).Approximately half of them (53.0%,35/66) had C T values > 34 and were not further investigated for genotyping purposes.All 31 Giardiapositive samples with qPCR C T values ≤ 34 were subjected to nested ssu-PCR to ascertain the assemblage of the parasite involved.Of them, 25.8% (8/31) were successfully genotyped at this locus (Table 5).Sequence analyses revealed that assemblage A (62.5%, 5/8) was more prevalent than assemblage B (37.5%, 3/8).Overall, MLST data at the four assessed loci were available for 3.0% (2/66) of samples, whereas subtyping data at a single locus (ssu rRNA) were available for 9.1% (6/66) of samples.No mixed infections nor host-adapted assemblages of canine (C, D), feline (F), or livestock (E) origin were detected.Out of the three assemblage A sequences at the ssu rRNA locus, two showed 100% identity with reference sequence M54878 with the remaining one differing from it by three SNPs in the form of ambiguous (double peak) positions.A single assemblage A sequence was confirmed as sub-assemblage AI at the gdh, bg, and tpi loci.The sequences generated at the three markers were identical to their respective reference sequences (Table 6).All three assemblage B sequences at the ssu rRNA locus showed 100% identity with reference AF113898.One of them was successfully genotyped at the three markers used, being identified as sub-assemblage BIV at the gdh marker and as sub-assemblage BIII at the tpi marker.This sample was, therefore, considered as an ambiguous BIII/BIV isolate (Table 6).

Discussion
This study shows that Cryptosporidium spp.and G. duodenalis are present at very different rates (2.4% vs. 27.9%) in faecal samples from Iberian lynxes without apparent gastrointestinal manifestations.The strengths of this study include (i) the use of molecular (PCR and Sanger sequencing) methods for accurate detection and genotyping of the two pathogens under investigation, (ii) a large sample size that includes a significant proportion (15-20%) of the estimated total population of free-living Iberian lynxes, (iii) representativeness of all three major distribution areas where the Iberian lynx is naturally present in Spain, (iv) the first report describing the molecular diversity of Cryptosporidium spp.and G. duodenalis in this carnivore host species, and (v) molecular evidence suggesting that a significant proportion of the positive samples might correspond to spurious infections as a direct consequence of predation on infected preys.
Cryptosporidiosis is regarded as a high-risk and often fatal opportunistic infection for undernourished young children and immunocompromised individuals as well as a major cause of neonatal diarrhoea in livestock [1][2][3].Comparatively, much less information is available on the epidemiology of Cryptosporidium spp. in wildlife with most studies conducted globally indicating low-to-medium infection rates and an apparent absence of gastrointestinal manifestations [2].This trend is particularly manifest in wild carnivore species.In the European scenario, Cryptosporidium infections have been reported in badgers (2.8-20.0%),foxes (6.1-13.3%),genets (16.6%),Eurasian lynxes (4.2%), martens (29.2-29.4%),minks (6.2%), otters (4.0%), raccoons (3.9-43.7%),raccoon dogs (24.1%), and wolves (35.7%), mostly with PCR (Table 1).Only two previous studies conducted in the Iberian Peninsula attempted to identify the presence of Cryptosporidium spp. in Iberian lynxes, but the limited number of samples analysed did not allow for the detection of the protozoa [20,32].In the present survey, Cryptosporidium spp. was detected in 2.4% (6/251) of the faecal samples from the Iberian lynxes examined, a figure in the lower range of those reported for other free-living carnivore species in Spain, Portugal, and other European countries.Despite the limited prevalence, positive animals were detected in the three sampling areas.These findings, together with the statistically significant differences among sampling periods, denote a wide but temporally heterogeneous circulation of Cryptosporidium in the Iberian lynx populations.
Molecular analyses of the six Cryptosporidium-positive isolates successfully genotyped revealed interesting data.First, four of the six infections detected were caused by Cryptosporidium species (C.alticolis, C. cuniculus, and C. occultus) with a strong preference for hosts that are common preys of the Iberian lynx.In this regard, although the Iberian lynx diet is mainly based on European rabbit, they can sporadically consume birds, wild ungulates, and also small mammals [67].Rodent-adapted Cryptosporidium alticolis and C. occultus were initially described in common voles and rats [68,69], whereas leporids, including rabbits and hares, are the preferred host species for C. cuniculus [70].Interestingly, C. alticolis has been previously reported in two red foxes in Poland [17].To our knowledge, this is the first report of C. cuniculus and C. occultus in free-living carnivores (including the Iberian lynx) globally.Taken together, these data seem to indicate that the presence of C. alticolis, C. cuniculus, and C. occultus in faecal samples from Iberian lynxes might be the consequence of spurious (mechanical carriage) rather than true infections.Second, the identification of generalist C. parvum allows for a wider interpretation.This Cryptosporidium species is characterised by a loose host specificity and great cross-species potential [71], making difficult the distinction between spurious and true infections.Regardless the case, the failure to amplify the two C. parvum isolates at the gp60 marker might be indicative of a low number of oocysts in faeces, compatible with a subclinical infection.Cryptosporidium parvum infections have been described in other European free-living carnivores, including wolves in Poland [15] and red foxes in Spain [20,21] and the UK [22].Third, we managed to characterise one of our two C. cuniculus isolates as genotype VaA19.Of note, previous studies conducted in Spain reported the presence of VaA16 (n = 1), VaA18 (n = 2), VbA24 (n = 1), VbA26 (n = 1), and VbA31 (n = 1) in wild populations of European rabbits and Iberian hares [72,73].These data expand our knowledge on the epidemiology of C. cuniculus in the country and support the spurious nature of our findings in Iberian lynxes.And fourth, the assignment of one of our Cryptosporidium-positive isolates as C. occultus should be interpreted with caution, as the generated ssu sequence was relatively short (214 bp) and this species is closely related to C. suis [69].We based our decision on two facts: (i) Our C. occultus sequence differed by two SNPs (688DelA, and T692A) with C. suis reference sequence AF115377, and (ii) the predator-prey relationship makes more likely that Iberian lynxes fed on small rodents than on suids, including domestic pigs and wild boars (the natural host species for C. suis).
Our molecular findings on the frequency and diversity of Cryptosporidium species in the Iberian lynx could also have public health implications.Whereas C. alticolis is not considered a zoonotic pathogen and only sporadic cases of human cryptosporidiosis by C. occultus have been reported in China [74], both C. cuniculus and C. parvum are able to cause significant morbidity in humans.Cryptosporidium cuniculus is typically identified at low infection (<1.5%) rates in European countries, including Spain [75], Sweden [76], and the UK [77,78].However, because C. cuniculus is closely related to C. hominis, its potential to cause human infections if the opportunity arises should not be underestimated [79].The finding of C. parvum has more relevance as this Cryptosporidium species causes one in four human cryptosporidiosis cases in Spain [80][81][82][83][84][85].
In contrast with cryptosporidiosis, giardiasis is widely regarded as a debilitating rather than a fatal condition in both human [5,86] and animal [9,10] hosts.Giardiasis in freeranging animals has only been investigated opportunistically, and relatively little is known about the epidemiology and health impact of the infection in wildlife populations [87].At the European level, Giardia infections have been reported in several wild carnivores, including badgers (25.6%), jackals (12.5%), lynxes (16.7%), martens (12.5-15.8%),otters (3.1-6.8%),raccoons (29.2-33.3%),red foxes (2.2-44.2%),wildcats (10.0%), and wolves (5.0-28.6%),mostly with PCR (Table 2).An infection rate of 26.7% was reported in 30 Iberian lynxes sampled from Portugal in a previous study [32], a figure very similar to that (27.9%) found in the present study also with PCR.These data suggest a high circulation of this parasite among the Iberian lynx populations and denote that the Iberian lynx could be a suitable host for G. duodenalis.The fact that neither geographical origin, sex, age, status, nor sampling year have an effect on the likelihood of having the parasite seems to support this hypothesis.
In the present study, the effort to assess the genetic diversity of G. duodenalis was hampered by the limited amount of parasitic DNA present in most positive samples, as indicated by the median qPCR C T value (34.5).This fact compromised the performance of our genotyping PCRs and explains why only a low proportion (25.8%, 8/31) of the tested G. duodenalis-positive samples were successfully characterised at one or more of the four (ssu, gdh, bg, and tpi) genetic markers used for this purpose.Our sequence analyses revealed the presence of two assemblages with assemblage A being more prevalent than assemblage B (62.5% vs. 37.5%, respectively).Remarkably, no feline-specific assemblage F was identified in the surveyed Iberian lynx populations.Considering that both assemblages A and B have zoonotic potential, these findings deserve attention.Out of the six assemblage A sequences, only one could be resolved at the sub-assemblage level as AI.This sub-assemblage is the most frequently found in animals [13], although it has also been reported at non-negligible rates in some human communities, primarily in low-income countries [88].The finding of assemblage B is somehow more worrying as this genetic variant is the most predominantly found circulating in the Spanish human population regardless of clinical status [89-92].Of note, in the only survey reporting molecular data on G. duodenalis infections in European free-living felids, assemblage B was identified in a single wildcat in Luxembourg [36].Taken together, these findings indicate that felids including the Iberian lynx can act as suitable hosts and spreaders of zoonotic variants of G. duodenalis.However, the finding that G. duodenalis infections are most likely associated with light parasite burdens (and, therefore, low cyst count in faeces) might limit the environmental contamination with infective cysts and reduce human exposure to them.
This study has some limitations that should be considered when interpreting the results obtained.First, it is possible that long-term storage of faecal samples has affected the quality/quantity of parasitic DNA, reducing the sensitivity and compromising the performance of the PCR protocols used for detection and genotyping purposes.Second, light parasitic infections leading to low (oo)cyst counts in faecal samples together with the limited sensibility of our genotyping PCRs have negatively impacted our ability to determine intra-species molecular variability in some Cryptosporidiumand G. duodenalispositive samples.And third, low Cryptosporidium infection rates might have compromised the accuracy of the statistical analyses conducted.

Conclusions
This study describes for the first time the occurrence and genetic diversity on Cryptosporidium spp.and G. duodenalis in the endangered Iberian lynx.The large sample size available, including animals from the main distribution areas, guarantee that the results obtained are representative of the whole free-living Iberian lynx population in Spain.Our results denote a limited but wide circulation of Cryptosporidium and a high wide and endemic distribution of Giardia among these individuals, which could be of animal health concern.The generated molecular data suggest that most Cryptosporidium species found correspond to rodent-or leporid-adapted strains that very likely cause spurious rather than true infection in the surveyed Iberian lynxes.However, the finding of zoonotic C. parvum and G. duodenalis assemblages A and B indicates that the Iberian lynx can act as a suitable host and spreader of these pathogens.Although the role of the Iberian lynx as a source of human cryptosporidiosis and giardiasis is regarded as low, this possibility should not be underestimated.Individuals (researchers, veterinarians, hunters) in close contact with infected animals or their faeces should be aware of the potential risk of zoonotic transmission of these protozoan parasites.The information provided in this study expands our knowledge on the epidemiology and public health relevance of Cryptosporidium spp.and G. duodenalis in Spain.
Institutional Review Board Statement: This study was carried out in accordance with Spanish legislation guidelines (RD 8/2003 of Animal Health and RD 53/2013 of Guiding Principles for Biomedical Research Involving Animals).The animal study protocol was approved in 2012 as part of the LIFE 10NAT/ES/570 IBERLINCE project ( https://webgate.ec.europa.eu/life/publicWebsite/project/LIFE10-NAT-ES-000570/recovering-the-historic-distribution-range-of-the-iberian-lynx-lynx-pardinu s-in-spain-and-portugal; accessed on 21 December 2023) and was adopted in the LIFE 19NAT/ ES001005LINXCONNECT project ( https://lifelynxconnect.eu/en/project/; accessed on 21 December 2023) in which sampling was carried out on live animals.Animals killed in road accidents and collected by the Wildlife Recovery Centres of the Regional Environmental Departments of the Autonomous Communities where this study was conducted were analysed in accordance with the corresponding authorisations (reference DGPFEN/SEN/avp_21_103_bis in Castilla-La Mancha, AB/is Exp.AUES/CYL/001/2021 in Castilla y León).
Informed Consent Statement: Not applicable.

Figure 1 .
Figure 1.Spatial distribution and molecular results of Iberian lynx samples.Total number of faecal samples analysed (n = 256) and frequency of positivity of Giardia duodenalis (*) and Cryptosporidium spp.(**) in each sampling area and captivity centre are shown in brackets.

Figure 1 .
Figure 1.Spatial distribution and molecular results of Iberian lynx samples.Total number of faecal samples analysed (n = 256) and frequency of positivity of Giardia duodenalis (*) and Cryptosporidium spp.(**) in each sampling area and captivity centre are shown in brackets.

Figure 2 .
Figure 2. Phylogenetic relationship among Cryptosporidium species and genotypes revealed with a maximum likelihood analysis of the partial ssu rDNA gene.Substitution rates were calculated by using the general time reversible model.Numbers on branches are percent bootstrapping values over 50% using 1000 replicates.The filled red circle indicates the nucleotide sequence generated in the present study.The filled green triangle indicates selected nucleotide sequences previously reported in wild carnivore species globally used for comparative purposes.

Figure 2 .
Figure 2. Phylogenetic relationship among Cryptosporidium species and genotypes revealed with a maximum likelihood analysis of the partial ssu rDNA gene.Substitution rates were calculated by using the general time reversible model.Numbers on branches are percent bootstrapping values over 50% using 1000 replicates.The filled red circle indicates the nucleotide sequence generated in the present study.The filled green triangle indicates selected nucleotide sequences previously reported in wild carnivore species globally used for comparative purposes.

Table 1 .
Infection rates and molecular diversity of Cryptosporidium spp. in European wild carnivore species, 2007-2023.

Table 2 .
Infection rates and molecular diversity of Giardia duodenalis in European wild carnivore species, 2007-2023.

Table 3 .
Infection rates by Cryptosporidium spp.and Giardia duodenalis in Iberian lynxes (n = 251) according to distribution area, sex, age, status, and sampling year of the animals.95% confidence intervals (95% CI) are indicated.

Table 4 .
Diversity, frequency, and molecular features of Cryptosporidium spp.isolates identified in the Iberian lynx population investigated in the present study.

Table 5 .
Multilocus sequence typing results of the eight G. duodenalis-positive samples successfully genotyped at any of the four loci investigated in the present survey.The age and gender of the infected Iberian lynxes are also shown.

Table 6 .
Diversity, frequency, and molecular features of G. duodenalis isolates identified in the Iberian lynx population investigated in the present study.