Effect of Myxoma Virus Species Jump on Iberian Hare Populations.

The myxoma virus species jump from European rabbits (Oryctolagus cuniculus) to Iberian hares (Lepus granatensis) has raised concerns. We assess the decline suffered by Iberian hare populations on the Iberian Peninsula and discuss the association between the effect of myxomatosis and the average abundance index, which we estimated by using hunting bags.

Monitoring of species' genetic diversity in Europe varies greatly and overlooks potential climate change impacts.

Genetic monitoring of populations currently attracts interest in the context of the Convention on Biological Diversity but needs long-term planning and investments. However, genetic diversity has been largely neglected in biodiversity monitoring, and when addressed, it is treated separately, detached from other conservation issues, such as habitat alteration due to climate change. We report an accounting of efforts to monitor population genetic diversity in Europe (genetic monitoring effort, GME), the evaluation of which can help guide future capacity building and collaboration towards areas most in need of expanded monitoring. Overlaying GME with areas where the ranges of selected species of conservation interest approach current and future climate niche limits helps identify whether GME coincides with anticipated climate change effects on biodiversity. Our analysis suggests that country area, financial resources and conservation policy influence GME, high values of which only partially match species' joint patterns of limits to suitable climatic conditions. Populations at trailing climatic niche margins probably hold genetic diversity that is important for adaptation to changing climate. Our results illuminate the need in Europe for expanded investment in genetic monitoring across climate gradients occupied by focal species, a need arguably greatest in southeastern European countries. This need could be met in part by expanding the European Union's Birds and Habitats Directives to fully address the conservation and monitoring of genetic diversity.

The potential role of scavenging flies as mechanical vectors of Lagovirus europaeus/GI.2.

The European rabbit (Oryctolagus cuniculus) populations of the Iberian Peninsula have been severely affected by the emergence of the rabbit haemorrhagic disease virus (RHDV) Lagovirus europaeus/GI.2 (RHDV2/b). Bushflies and blowflies (Muscidae and Calliphoridae families, respectively) are important RHDV vectors in Oceania, but their epidemiological role is unknown in the native range of the European rabbit. In this study, scavenging flies were collected between June 2018 and February 2019 in baited traps at one site in southern Portugal, alongside a longitudinal capture-mark-recapture study of a wild European rabbit population, aiming to provide evidence of mechanical transmission of GI.2 by flies. Fly abundance, particularly from Calliphoridae and Muscidae families, peaked in October 2018 and in February 2019. By employing molecular tools, we were able to detect the presence of GI.2 in flies belonging to the families Calliphoridae, Muscidae, Fanniidae and Drosophilidae. The positive samples were detected during an RHD outbreak and absent in samples collected when no evidence of viral circulation in the local rabbit population was found. We were able to sequence a short viral genomic fragment, confirming its identity as RHDV GI.2. The results suggest that scavenging flies may act as mechanical vectors of GI.2 in the native range of the southwestern Iberian subspecies O. cuniculus algirus. Future studies should better assess their potential in the epidemiology of RHD and as a tool for monitoring viral circulation in the field.

The evolution of white-tailed jackrabbit camouflage in response to past and future seasonal climates.

The genetic basis of adaptive traits has rarely been used to predict future vulnerability of populations to climate change. We show that light versus dark seasonal pelage in white-tailed jackrabbits (Lepus townsendii) tracks snow cover and is primarily determined by genetic variation at endothelin receptor type B (EDNRB), corin serine peptidase (CORIN), and agouti signaling protein (ASIP). Winter color variation was associated with deeply divergent alleles at these genes, reflecting selection on both ancestral and introgressed variation. Forecasted reductions in snow cover are likely to induce widespread camouflage mismatch. However, simulated populations with variation for darker winter pelage are predicted to adapt rapidly, providing a trait-based genetic framework to facilitate evolutionary rescue. These discoveries demonstrate how the genetic basis of climate change adaptation can inform conservation.

Characterising Mitochondrial Capture in an Iberian Shrew.

Mitochondrial introgression raises questions of biogeography and of the extent of reproductive isolation and natural selection. Previous phylogenetic work on the Sorex araneus complex revealed apparent mitonuclear discordance in Iberian shrews, indicating past hybridisation of Sorex granarius and the Carlit chromosomal race of S. araneus, enabling introgression of the S. araneus mitochondrial genome into S. granarius. To further study this, we genetically typed 61 Sorex araneus/coronatus/granarius from localities in Portugal, Spain, France, and Andorra at mitochondrial, autosomal, and sex-linked loci and combined our data with the previously published sequences. Our data are consistent with earlier data indicating that S. coronatus and S. granarius are the most closely related of the three species, confirming that S. granarius from the Central System mountain range in Spain captured the mitochondrial genome from a population of S. araneus. This mitochondrial capture event can be explained by invoking a biogeographical scenario whereby S. araneus was in contact with S. granarius during the Younger Dryas in central Iberia, despite the two species currently having disjunct distributions. We discuss whether selection favoured S. granarius with an introgressed mitochondrial genome. Our data also suggest recent hybridisation and introgression between S. coronatus and S. granarius, as well as between S. araneus and S. coronatus.

Genetics of randomly bred cats support the cradle of cat domestication being in the Near East.

Cat domestication likely initiated as a symbiotic relationship between wildcats (Felis silvestris subspecies) and the peoples of developing agrarian societies in the Fertile Crescent. As humans transitioned from hunter-gatherers to farmers ~12,000 years ago, bold wildcats likely capitalized on increased prey density (i.e., rodents). Humans benefited from the cats' predation on these vermin. To refine the site(s) of cat domestication, over 1000 random-bred cats of primarily Eurasian descent were genotyped for single-nucleotide variants and short tandem repeats. The overall cat population structure suggested a single worldwide population with significant isolation by the distance of peripheral subpopulations. The cat population heterozygosity decreased as genetic distance from the proposed cat progenitor's (F.s. lybica) natural habitat increased. Domestic cat origins are focused in the eastern Mediterranean Basin, spreading to nearby islands, and southernly via the Levantine coast into the Nile Valley. Cat population diversity supports the migration patterns of humans and other symbiotic species.

A study of viral pathogens in bat species in the Iberian Peninsula: identification of new coronavirus genetic variants.

Bats have long been associated with multiple pathogens, including viruses affecting humans such as henipaviruses, filoviruses, bunyaviruses and coronaviruses. The alpha and beta coronaviruses genera can infect most mammalian species. Among them, betacoronavirus SARS-CoV, MERS-CoV and SARS-CoV-2, which have caused the three major pandemics in the last two decades, have been proposed to originate in bats. In this study, 194 oral swabs from 22 bats species sampled in 19 locations of the Iberian Peninsula were analysed and characterized by three different PCR tests (coronavirus generic real-time RT-PCR, multiplex conventional PCR, and SARS-CoV-2 specific real-time RT-PCR) to detect bat coronaviruses. Screening with coronavirus generic PCR showed 102 positives out of 194 oral swabs analysed. Then, metabarcoding with multiplex PCR amplified 15 positive samples. Most of the coronaviruses detected in this study belong to alphacoronavirus (α-CoV) genus, with multiple alphacoronaviruses identified by up to five different genetic variants coexisting in the same bat. One of the positive samples identified in a Miniopterus schreibersii bat positive for the generic coronavirus PCR and the specific SARS-CoV-2 PCR was classified as betacoronavirus (-CoV) through phylogenetic analysis. These results support the rapid evolution of coronaviruses to generate new genomic potentially pathogenic variants likely through co-infection and recombination.

Absence of anti-Besnoitia spp. specific antibodies in European wild lagomorphs from the Iberian Peninsula.

Besnoitia besnoiti is an apicomplexan parasite whose life cycle is not completely understood. It is assumed that this parasite might have an indirect life cycle with a carnivore as a definitive host able to shed oocysts after the ingestion of mature cysts in tissues of an infected intermediate host. Cattle and wild cervids on the Iberian Peninsula can act as intermediate hosts of B. besnoiti, and exposure to the parasite has been demonstrated in equids. In this study, we aimed to assess the presence of members of the genera Besnoitia in wild lagomorphs from the Iberian Peninsula and the potential role of these host species in the life cycle of B. besnoiti, as all the animals were sampled from 19 regions of the Iberian Peninsula where cases of bovine besnoitiosis have been previously detected. Serum samples (Oryctolagus cuniculus: n = 552; Lepus europaeus: n = 122) were first analysed by ELISA and subsequently confirmed by Western blot (WB). Specific antibodies against B. besnoiti were not found in any sampled animal by WB. In addition, lung samples from a subset of wild rabbits (n = 16) were tested by PCR and Besnoitia spp. DNA was not detected. These results suggest that Besnoitia spp. are unlikely to circulate in wild lagomorphs in the Iberian Peninsula. Thus, lagomorphs are not expected to play a key role in the biological cycle of B. besnoiti. Further studies are necessary to assess whether different micromammal species, such as rodents, can serve as natural reservoirs of Besnoitia spp. in other European regions.

The evolutionary pathways for local adaptation in mountain hares.

Understanding the evolution of local adaptations is a central aim of evolutionary biology and key for the identification of unique populations and lineages of conservation relevance. By combining RAD sequencing and whole-genome sequencing, we identify genetic signatures of local adaptation in mountain hares (Lepus timidus) from isolated and distinctive habitats of its wide distribution: Ireland, the Alps and Fennoscandia. Demographic modelling suggested that the split of these mountain hares occurred around 20 thousand years ago, providing the opportunity to study adaptive evolution over a short timescale. Using genome-wide scans, we identified signatures of extreme differentiation among hares from distinct geographic areas that overlap with area-specific selective sweeps, suggesting targets for local adaptation. Several identified candidate genes are associated with traits related to the uniqueness of the different environments inhabited by the three groups of mountain hares, including coat colour, ability to live at high altitudes and variation in body size. In Irish mountain hares, a variant of ASIP, a gene previously implicated in introgression-driven winter coat colour variation in mountain and snowshoe hares (L. americanus), may underlie brown winter coats, reinforcing the repeated nature of evolution at ASIP moulding adaptive seasonal colouration. Comparative genomic analyses across several hare species suggested that mountain hares' adaptive variants appear predominantly species-specific. However, using coalescent simulations, we also show instances where the candidate adaptive variants have been introduced via introgressive hybridization. Our study shows that standing adaptive variation, including that introgressed from other species, was a crucial component of the post-glacial dynamics of species.

Stepping up from wildlife disease surveillance to integrated wildlife monitoring in Europe.

In a context of disease emergence and faced with the ever-growing evidence of the role of wildlife in the epidemiology of transmissible diseases, efforts have been made to develop wildlife disease surveillance (WDS) programs throughout Europe. Disease monitoring is ideally composed of "numerator data" (number of infected individuals) and "denominator data" (size of the target population). Too often however, information is available for only one. Hence, there is a need for developing integrated and harmonized disease and population monitoring tools for wildlife: integrated wildlife monitoring (IWM). IWM should have three components. Passive disease surveillance improves the likelihood of early detection of emerging diseases, while active surveillance and population monitoring are required to assess epidemiological dynamics, freedom of disease, and the outcome of interventions. Here, we review the characteristics of ongoing WDS in Europe, observe how pathogens have been ranked, and note a need for ranking host species, too. Then, we list the challenges for WDS and draw a roadmap for stepping up from WDS to IWM. There is a need to integrate and maintain an equilibrium between the three components of IWM, improve data collection and accessibility, and guarantee the adaptability of these schemes to each epidemiological context and temporal period. Methodological harmonization and centralization of information at a European level would increase efficiency of national programs and improve the follow-up of eventual interventions. The ideal IWM would integrate capacities from different stakeholder; allow to rapidly incorporate relevant new knowledge; and rely on stable capacities and funding.

Assessing changes in stream macroinvertebrate communities across ecological gradients using morphological versus DNA metabarcoding approaches.

Freshwater macroinvertebrates provide valuable indicators for biomonitoring ecosystem change in relation to natural and anthropogenic drivers. DNA metabarcoding is an efficient approach for estimating such indicators, but its results may differ from morphotaxonomic approaches traditionally used in biomonitoring. Here we test the hypothesis that despite differences in the number and identity of taxa recorded, both approaches may retrieve comparable patterns of community change, and detect similar ecological gradients influencing such changes. We compared results obtained with morphological identification at family level of macroinvertebrates collected at 80 streams under a Water Framework Directive biomonitoring program in Portugal, with results obtained with metabarcoding from the ethanol preserving the bulk samples, using either single (COI-M19BR2, 16S-Inse01, 18S-Euka02) or multiple markers. Metabarcoding recorded less families and different communities compared to morphotaxonomy, but community sensitivities to disturbance estimated with the IASPT index were more similar across approaches. Spatial variation in local community metrics and the factors influencing such variation were significantly correlated between morphotaxonomy and metabarcoding. After reducing random noise in the dissimilarity matrices, the spatial variation in community composition was also significantly correlated across methods. A dominant gradient of community change was consistently retrieved, and all methods identified a largely similar set of anthropogenic stressors strongly influencing such gradient. Overall, results confirm our initial hypothesis, suggesting that morphotaxonomy and metabarcoding can estimate consistent spatial patterns of community variation and their main drivers. These results are encouraging for macroinvertebrate biomonitoring using metabarcoding approaches, suggesting that they can be intercalibrated with morphotaxonomic approaches to recover equivalent spatial and temporal gradients of ecological change.

Sex and Age-Specific Hematology and Biochemistry Reference Intervals of Live Iberian Hares (Lepus granatensis) and Comparison with Postmortem Sampling.

Hematology and serum biochemistry reference intervals were estimated for the Iberian hare (Lepus granatensis). Most parameters differed significantly between hunted and livetrapped Iberian hares. Significant differences were found for sex (red blood cell count, mean corpuscular volume, total protein, albumin, uric acid, triglycerides, cholesterol, chloride) and age classes (red blood cell count, hematocrit, mean corpuscular hemoglobin, glucose, calcium, and sodium). Sex- and age-specific reference intervals were estimated for these parameters. Red blood cell count, hemoglobin concentration, mean corpuscular hemoglobin, urea, and potassium show seasonal variations, with the lowest values in summer and the highest in winter. Creatinine, calcium, sodium, and phosphorus achieve their highest values in summer and stable baseline values throughout the rest of the year. These reference intervals can be used as baseline to monitor health, physiology, ecology, and nutrition of Iberian hare populations.

The Legacy of Recurrent Introgression during the Radiation of Hares.

Hybridization may often be an important source of adaptive variation, but the extent and long-term impacts of introgression have seldom been evaluated in the phylogenetic context of a radiation. Hares (Lepus) represent a widespread mammalian radiation of 32 extant species characterized by striking ecological adaptations and recurrent admixture. To understand the relevance of introgressive hybridization during the diversification of Lepus, we analyzed whole exome sequences (61.7 Mb) from 15 species of hares (1-4 individuals per species), spanning the global distribution of the genus, and two outgroups. We used a coalescent framework to infer species relationships and divergence times, despite extensive genealogical discordance. We found high levels of allele sharing among species and show that this reflects extensive incomplete lineage sorting and temporally layered hybridization. Our results revealed recurrent introgression at all stages along the Lepus radiation, including recent gene flow between extant species since the last glacial maximum but also pervasive ancient introgression occurring since near the origin of the hare lineages. We show that ancient hybridization between northern hemisphere species has resulted in shared variation of potential adaptive relevance to highly seasonal environments, including genes involved in circadian rhythm regulation, pigmentation, and thermoregulation. Our results illustrate how the genetic legacy of ancestral hybridization may persist across a radiation, leaving a long-lasting signature of shared genetic variation that may contribute to adaptation. [Adaptation; ancient introgression; hybridization; Lepus; phylogenomics.].

Quantification of the Animal Tuberculosis Multi-Host Community Offers Insights for Control.

Animal tuberculosis (TB) is a multi-host zoonotic disease whose prevalence in cattle herds in Europe has been increasing, despite a huge investment in eradication. The composition of the host community is a fundamental driver of pathogen transmission, and yet this has not been formally quantified for animal TB in Europe. We quantified multi-host communities of animal TB, using stochastic models to estimate the number of infected domestic and wild hosts in three regions: officially TB-free Central-Western Europe, and two largely TB-endemic regions, the Iberian Peninsula and Britain and Ireland. We show that the estimated number of infected animals in the three regions was 290,059-1,605,612 and the numbers of infected non-bovine domestic and wild hosts always exceeded those of infected cattle, with ratios ranging from 3.3 (1.3-19.6):1 in Britain and Ireland to 84.3 (20.5-864):1 in the Iberian Peninsula. Our results illustrate for the first time the extent to which animal TB systems in some regions of Europe are dominated by non-bovine domestic and wild species. These findings highlight the need to adapt current strategies for effective future control of the disease.

Ecological traits and the spatial structure of competitive coexistence among carnivores.

Competition is a widespread interaction among carnivores, ultimately manifested through one or more dimensions of the species' ecological niche. One of the most explicit manifestations of competitive interactions regards spatial displacement. Its interpretation under a theoretical context provides an important tool to deepen our understanding of biological systems and communities, but also for wildlife management and conservation. We used Bayesian multispecies occupancy models on camera-trapping data from multiple sites in Southwestern Europe (SWE) to investigate competitive interactions within a carnivore guild, and to evaluate how species' ecological traits are shaping coexistence patterns. Seventeen out of 26 pairwise interactions departed from a hypothesis of independent occurrence, with spatial association being twice as frequent as avoidance. Association behaviors were only detected among mesocarnivores, while avoidance mainly involved mesocarnivores avoiding the apex predator (n = 4) and mesocarnivore-only interactions (n = 2). Body mass ratios, defined as the dominant over the subordinate species body mass, revealed an important negative effect ( ß^=-0.38;CI95=-0.81to-0.06 ) on co-occurrence probability, and support that spatially explicit competitive interactions are mostly expressed by larger species able to dominate over smaller ones, with a threshold in body mass ratios of ~4, above which local-scale intraguild coexistence is unlikely. We found a weak relationship between pairwise trophic niche overlap and the probability of coexistence ( ß^=-0.19;CI95=-0.58to0.21 ), suggesting that competition for feeding resources may not be a key driver of competition, at least at the scale of our analysis. Despite local-scale avoidance, regional-scale coexistence appears to be maintained by the spatial structuring of the competitive environment. We provide evidence that SWE ecosystems consist of spatially structured competitive environments, and propose that coexistence among near-sized species is likely achieved through the interplay of "facultative" and "behavioral" character displacements. Factors influencing carnivore coexistence likely include context-dependent density and trait-mediated effects, which should be carefully considered for a sound understanding of the mechanisms regulating these communities.

Transcriptomic regulation of seasonal coat color change in hares.

Color molts from summer brown to winter white coats have evolved in several species to maintain camouflage year-round in environments with seasonal snow. Despite the eco-evolutionary relevance of this key phenological adaptation, its molecular regulation has only recently begun to be addressed. Here, we analyze skin transcription changes during the autumn molt of the mountain hare (Lepus timidus) and integrate the results with an established model of gene regulation across the spring molt of the closely related snowshoe hare (L. americanus). We quantified differences in gene expression among three stages of molt progression-"brown" (early molt), "intermediate," and "white" (late molt). We found 632 differentially expressed genes, with a major pulse of expression early in the molt, followed by a milder one in late molt. The functional makeup of differentially expressed genes anchored the sampled molt stages to the developmental timeline of the hair growth cycle, associating anagen to early molt and the transition to catagen to late molt. The progression of color change was characterized by differential expression of genes involved in pigmentation, circadian, and behavioral regulation. We found significant overlap between differentially expressed genes across the seasonal molts of mountain and snowshoe hares, particularly at molt onset, suggesting conservatism of gene regulation across species and seasons. However, some discrepancies suggest seasonal differences in melanocyte differentiation and the integration of nutritional cues. Our established regulatory model of seasonal coat color molt provides an important mechanistic context to study the functional architecture and evolution of this crucial seasonal adaptation.

Applying genomic data in wildlife monitoring: Development guidelines for genotyping degraded samples with reduced single nucleotide polymorphism panels.

The genomic era has led to an unprecedented increase in the availability of genome-wide data for a broad range of taxa. Wildlife management strives to make use of these vast resources to enable refined genetic assessments that enhance biodiversity conservation. However, as new genomic platforms emerge, problems remain in adapting the usually complex approaches for genotyping of noninvasively collected wildlife samples. Here, we provide practical guidelines for the standardized development of reduced single nucleotide polymorphism (SNP) panels applicable for microfluidic genotyping of degraded DNA samples, such as faeces or hairs. We demonstrate how microfluidic SNP panels can be optimized to efficiently monitor European wildcat (Felis silvestris S.) populations. We show how panels can be set up in a modular fashion to accommodate informative markers for relevant population genetics questions, such as individual identification, hybridization assessment and the detection of population structure. We discuss various aspects regarding the implementation of reduced SNP panels and provide a framework that will allow both molecular ecologists and practitioners to help bridge the gap between genomics and applied wildlife conservation.

An Annotated Draft Genome of the Mountain Hare (Lepus timidus).

Hares (genus Lepus) provide clear examples of repeated and often massive introgressive hybridization and striking local adaptations. Genomic studies on this group have so far relied on comparisons to the European rabbit (Oryctolagus cuniculus) reference genome. Here, we report the first de novo draft reference genome for a hare species, the mountain hare (Lepus timidus), and evaluate the efficacy of whole-genome re-sequencing analyses using the new reference versus using the rabbit reference genome. The genome was assembled using the ALLPATHS-LG protocol with a combination of overlapping pair and mate-pair Illumina sequencing (77x coverage). The assembly contained 32,294 scaffolds with a total length of 2.7 Gb and a scaffold N50 of 3.4 Mb. Re-scaffolding based on the rabbit reference reduced the total number of scaffolds to 4,205 with a scaffold N50 of 194 Mb. A correspondence was found between 22 of these hare scaffolds and the rabbit chromosomes, based on gene content and direct alignment. We annotated 24,578 protein coding genes by combining ab-initio predictions, homology search, and transcriptome data, of which 683 were solely derived from hare-specific transcriptome data. The hare reference genome is therefore a new resource to discover and investigate hare-specific variation. Similar estimates of heterozygosity and inferred demographic history profiles were obtained when mapping hare whole-genome re-sequencing data to the new hare draft genome or to alternative references based on the rabbit genome. Our results validate previous reference-based strategies and suggest that the chromosome-scale hare draft genome should enable chromosome-wide analyses and genome scans on hares.

Genomic approaches to identify hybrids and estimate admixture times in European wildcat populations.

The survival of indigenous European wildcat (Felis silvestris silvestris) populations can be locally threatened by introgressive hybridization with free-ranging domestic cats. Identifying pure wildcats and investigating the ancestry of admixed individuals becomes thus a conservation priority. We analyzed 63k cat Single Nucleotide Polymorphisms (SNPs) with multivariate, Bayesian and gene-search tools to better evaluate admixture levels between domestic and wild cats collected in Europe, timing and ancestry proportions of their hybrids and backcrosses, and track the origin (wild or domestic) of the genomic blocks carried by admixed cats, also looking for possible deviations from neutrality in their inheritance patterns. Small domestic ancestry blocks were detected in the genomes of most admixed cats, which likely originated from hybridization events occurring from 6 to 22 generations in the past. We identified about 1,900 outlier coding genes with excess of wild or domestic ancestry compared to random expectations in the admixed individuals. More than 600 outlier genes were significantly enriched for Gene Ontology (GO) categories mainly related to social behavior, functional and metabolic adaptive processes (wild-like genes), involved in cognition and neural crest development (domestic-like genes), or associated with immune system functions and lipid metabolism (parental-like genes). These kinds of genomic ancestry analyses could be reliably applied to unravel the admixture dynamics in European wildcats, as well as in other hybridizing populations, in order to design more efficient conservation plans.

Gastrointestinal parasite infestation in the alpine mountain hare (Lepus timidus varronis): Are abiotic environmental factors such as elevation, temperature and precipitation affecting prevalence of parasite species?

Information concerning factors regulating Alpine mountain hare (Lepus timidus varronis) populations such as host-parasite interactions is missing as only a few parasitological surveys exist of this subspecies. Parasites are not only dependent on their host but also on suitable environmental conditions for infestation. Abiotic environmental factors have an important regulating role on parasites in mammals. It is estimated that the elevation range of parasites is likely to shift in response to alternate host movement and changes in climate. Here we assess the parasitic infestation in the Alpine mountain hare by analysing the parasites in faeces and comparing the parasite infestation at different elevation ranges and at varied weather conditions for two years in the Austrian Alps. Almost half of the faecal samples were free of parasites (46.2%, n = 52). Most frequent was the infection by Coccidia (46.2%), whereas stomach intestine strongylids, Trichuris spp, and Cestoda were only found in 9.6% of all faeces. Hence, only Coccidia may be prevalent enough to regulate Alpine mountain hare populations in the Austrian Alps. Elevation had a significant positive effect on the infection of animals by Trichuris spp, whereas temperature had a significant negative effect on the infection by any parasite traceable in faeces and, when looking at the parasite groups individually, on Coccidia.