HAT: de novo variant calling for highly accurate short-read and long-read sequencing data.

MOTIVATION: de novo variants (DNVs) are variants that are present in offspring but not in their parents. DNVs are both important for examining mutation rates as well as in the identification of disease-related variation. While efforts have been made to call DNVs, calling of DNVs is still challenging from parent-child sequenced trio data. We developed Hare And Tortoise (HAT) as an automated DNV detection workflow for highly accurate short-read and long-read sequencing data. Reliable detection of DNVs is important for human genomics and HAT addresses this need. RESULTS: HAT is a computational workflow that begins with aligned read data (i.e. CRAM or BAM) from a parent-child sequenced trio and outputs DNVs. HAT detects high-quality DNVs from Illumina short-read whole-exome sequencing, Illumina short-read whole-genome sequencing, and highly accurate PacBio HiFi long-read whole-genome sequencing data. The quality of these DNVs is high based on a series of quality metrics including number of DNVs per individual, percent of DNVs at CpG sites, and percent of DNVs phased to the paternal chromosome of origin. AVAILABILITY AND IMPLEMENTATION: https://github.com/TNTurnerLab/HAT.

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.

Patterns of evolution in MHC class II DQA and DQB exon 2 genes of Alpine mountain hares, Lepus timidus varronis, and sympatric and parapatric brown hares, L. europaeus, from Switzerland.

In natural populations, hybridization is known to occur between a wide range of species. However, its evolutionary significance is less clear. Genes involved in fighting pathogens are considered excellent candidates for studying adaptive introgression, although both introgression and balancing selection can generate similar patterns of diversity and differentiation. Here, we compared DQA and DQB MHC class II and microsatellite allelic diversity of sympatric and parapatric mountain (Lepus timidus) and brown hare (L. europaeus) populations from Switzerland. We detected higher genetic diversity in brown hares compared to mountain hares at both MHC and microsatellite loci. We consider the observed patterns of microsatellite diversity both for L. europaeus and L. timidus as result of stochastic demographic processes while the pattern of MHC polymorphism of the studied hare populations can be explained by pathogen-driven selection. Rare bidirectional gene flow between both hare species seems to occur specifically for MHC alleles. However, the high number of shared alleles showing similar high frequency in both species suggests that reciprocally exchanged MHC alleles are being maintained via balancing selection. Adaptation to similar pathogen communities can also lead to parallel selection of MHC alleles. Positive selection, recombination and mutations have played different roles in shaping the patterns of MHC allelic diversity in and differentiation between both species. Results for the latter evolutionary forces do not show a better matching between the sympatric populations compared to the parapatric ones, suggesting a minor role of introgression for the observed evolutionary patterns of the studied hare species.

Vegetation biomass and topography are associated with seasonal habitat selection and fall translocation behavior in Arctic hares.

Habitat selection theory suggests that environmental features selected at coarse scales reveal fundamental factors affecting animal fitness. When these factors vary across seasons, they may lead to large-scale movements, including long-distance seasonal migrations. We analyzed the seasonal habitat selection of 25 satellite-tracked Arctic hares from a population on Ellesmere Island (Nunavut, Canada) that relocated over 100 km in the fall. Since no other lagomorph is known to perform such extensive movements, this population offered an ideal setting to test animal movement and habitat selection theory. On summer grounds hares selected low elevation areas, while on winter grounds they selected high vegetation biomass, high elevation, and steep slopes. During fall relocation, they alternated between stopover and traveling behavioral states (ratio 2:1). Stopover locations were characterized by higher vegetation heterogeneity and lower rugosity than traveling locations, while vegetation biomass and elevation interacted to explain stopover locations in a more complex way. The selected combination of environmental features thus varied across seasons and behavioral states, in a way broadly consistent with predictions based on the changing food and safety needs of hares. Although causality was not demonstrated, our results improve our understanding of long-distance movements and habitat selection in Arctic hares, as well as herbivore ecology in the polar desert. Results also provide strong support to animal movement and habitat selection theory, by showing how some important hypotheses hold when tested in a species phylogenetically distinct from most animal models used in this research field.

Typing of Leishmania isolates from vectors and leporids of the Madrid (Spain) outbreak.

In 2009, a large outbreak of leishmaniasis, associated with environmental changes, was declared near Madrid (Spain), in which Phlebotomus perniciosus was the vector, whereas the main reservoirs were hares and rabbits. Analysis of isolates from humans, vectors and leporids from the focus identified the Leishmania infantum ITS-Lombardi genotype. However, multilocus enzyme electrophoresis (MLEE), the reference technique for Leishmania typing, and sequencing of the hsp70 gene, a commonly used marker, were not performed. In the present study, 19 isolates from P. perniciosus (n = 11), hares (n = 5) and rabbits (n = 3) from the outbreak area, all characterized as ITS-Lombardi in previous studies, were analysed by MLEE and hsp70 sequencing. The hsp70 results confirmed that all the analysed strains are L. infantum. However, by MLEE, 4 different zymodemes of L. infantum were identified based on variable mobilities of the NP1 enzyme: MON-34 (NP1100, n = 11), MON-80 (NP1130, n = 6), MON-24 (NP1140, n = 1) and MON-331 (NP1150, n = 1). The relative frequency of these zymodemes does not correspond to their usual occurrence in Spain. Moreover, MON-34 and MON-80 were found in P. perniciosus, hares and rabbits for the first time. These findings continue to provide insights into the outbreak and call for further studies with a higher number of strains.

Estrogen regulated transcription of the non-estrogen-regulated hamster uteroglobin gene in MCF-7 cells.

To study the estrogen regulated transcription of the uteroglobin (UG) gene, the founding member of the secretoglobin family widely expressed in many different mammalian species, we re-created functional estrogen response elements (EREs) in the UG gene promoter from a species where UG expression is not regulated by estrogens: the hamster Mesocricetus auratus (Ma), to ascertain if the lack of functional EREs is the real cause of its estrogen insensitivity. Functional EREs in the hamster promoter, including the consensus ERE (cERE), failed to respond to an appropriate estrogen stimulus compared with its estrogen regulated ortholog from the brown hare Lepus capensis (Lc). As the nucleotide sequence is the only difference between genetic constructs from these two species, we suspected that the UG promoter from the hamster probably contains cis-acting genetic elements that negatively impairs the estrogen-regulated transcription mediated by the functional ERE. Accordingly, we prepared chimeric DNA constructs which eventually allowed to identify a region located 29 base pairs (bp) downstream of the ERE as responsible for the lack of estrogen-responsiveness of the Ma-UG gene in the breast cancer cell line MCF-7. This region contains the sequence ACACCCC which has been identified as the core sequence of the Sp/ Krüppel-like factor (KLF) family of transcription factors. This finding is relevant, not only due to the observation on a novel mechanism that control estrogen-induced transcription, but also because it may encourage further investigation for better defining specific genes with an ERE that do not respond to estrogen signaling in MCF-7 cells, a cell line widely employed as an in vitro model in breast cancer research.

Selective sweep analysis of the adaptability of the Yarkand hare (Lepus yarkandensis) to hot arid environments using SLAF-seq.

The Yarkand hare (Lepus yarkandensis) inhabits arid desert areas and is endemic to China. It has evolved various adaptations to survive in hot arid environments, including stress responses, the ability to maintain water homeostasis and heat tolerance. Here, we performed a selective sweep analysis to identify the candidate genes for adaptation to hot arid environments in the Yarkand hare. A total of 397 237 single-nucleotide polymorphisms were obtained from 80 Yarkand hares, which inhabit hot arid environments, and 36 Tolai hares (Lepus tolai), which inhabit environments with a mild climate, via specific-locus amplified fragment sequencing. We identified several candidate genes that were associated with the heat stress response (HSPE1), oxidative stress response (SLC23A and GLRX2), immune response (IL1R1 and IRG1), central nervous system development (FGF13, THOC2, FMR1 and MECP2) and regulation of water homeostasis (CDK1) according to fixation index values and θπ ratios in the selective sweep analysis, and six of these genes (GLRX2, IRG1, FGF13, FMR1, MECP2 and CDK1) are newly discovered genes. To the best of our knowledge, this is the first study to identify candidate genes for adaptation to hot arid environments in the Yarkand hare. The results of this study enhance our understanding of the adaptation of the Yarkand hare to hot arid environments and will aid future studies aiming to functionally verify these candidate genes.

Toxoplasma gondii infection in European brown hares (Lepus europaeus) in Brazil.

Toxoplasmosis is a worldwide zoonosis that affects warm-blooded animals, including humans. Wild animals can act as intermediate hosts of this pathogen; thus, this study aims to detect Toxoplasma gondii infection in invasive European brown hares in Brazil. For this, 72 wild European brown hares were captured from July 2020 to June 2022 in three Brazilian states: São Paulo, Paraná, and Rio Grande do Sul. The diagnostic of Toxoplasma gondii infection was performed by bioassay in mouse, histopathology in Hematoxylin-Eosin-stained tissue sections (brain, liver, lungs, kidneys, and small intestine), serology by IFAT, and molecular techniques by conventional PCR and qPCR. The combined prevalence of the different diagnostic methods was 51.4% (37/72, CI= 40.1 - 62.6 %), and there was no statistical difference between sexes, age range, or geographical region of the hosts. Mouse bioassay was the technique that detected more positive hares. To our knowledge, this is the first confirmation of Toxoplasma gondii infection in invasive European brown hares in Brazil. These animals act as reservoirs and potential infection source for carnivores and other wild and domestic animals, including humans, thus contributing to perpetuate the disease cycle in São Paulo, Paraná, and Rio Grande do Sul States. Research such as the present study is necessary to raise awareness about the role of animals in the disease cycle.

Potentially toxic metals in seawater, sediment and seaweeds: bioaccumulation, ecological and human health risk assessment.

This study assesses the bioaccumulation, ecological, and health risks associated with potentially toxic metals (PTMs), including Pb, Hg, Cd, As, and Cr in Hare Island, Thoothukudi. The results revealed that the concentration of PTMs in sediment, seawater, and S. wightii ranged from 0.095 to 2.81 mg kg-1, 0.017 to 1.515 mg L-1, and 0.076 to 5.713 mg kg-1, respectively. The highest concentrations of PTMs were found in the S. wightii compared to seawater and sediment. The high bioaccumulation of Hg and As in S. wightii suggests that it can be used as a bioindicator for these elements in this region. The ecological risk indices, which include individual, complex, biological, and ecological pollution indices, suggest that Hare Island had moderate contamination with Hg and Cd. However, there are no human health risks associated with PTMs. This study examines the current ecological and health risks associated with PTMs and emphasizes the importance of regular monitoring.

New Data on the Ecology of the Fossil Don Hare Lepus tanaiticus Gureev, 1964.

A spore and pollen research was performed for the first time to study the contents of the gastrointestinal tracts (GITs) and sediments containing frozen mummies of the fossil Don hare from the Upper Pleistocene ice complex of the Verkhoyansk district of Yakutia. Radiocarbon dating (C14) revealed that the hares lived during the Karginian Interstadial of the Late Pleistocene, 32.5 thousand years ago (calibrated date). The results expanded the understanding of the ecology of extinct Lepus tanaiticus. The species was assumed to live in cold steppes dominated by xerophytic communities, as well as in grass-forb and sedge-forb meadows. Herbaceous plants mostly constituted the winter diet of Don hare in contrast to the modern mountain hare L. timidus, which feeds mainly on branches and bark of trees and shrubs in winter.

Chromosome-scale genome assembly of Lepus oiostolus (Lepus, Leporidae).

Lepus oiostolus (L. oiostolus) is a species endemic to the Qinghai-Tibet Plateau. However, the absence of a reference genome limits genetic studies. Here, we reported a high-quality L. oiostolus genome assembly, with scaffolds anchored to 24 chromosomes and a total assembled length of 2.80 Gb (contig N50 = 64.25 Mb). Genomic annotation uncovered 22,295 protein-coding genes and identified 49.84% of the sequences as transposable elements. Long interspersed nuclear elements (LINEs) constitute a high proportion of the genome. Our study is at the first time to report the chromosome-scale genome for the species of the L. oiostolus. It provides a valuable genomic resource for future research on the evolution of the Leporidae.

Serologic Survey of Selected Arthropod-Borne Pathogens in Free-Ranging Snowshoe Hares (Lepus americanus) Captured in Northern Michigan, USA.

Snowshoe hares (Lepus americanus) in the Upper Peninsula (UP) of Michigan, USA, occupy the southern periphery of the species' range and are vulnerable to climate change. In the eastern UP, hares are isolated by the Great Lakes, potentially exacerbating exposure to climate-change-induced habitat alterations. Climate change is also measurably affecting distribution and prevalence of vector-borne pathogens in North America, and increases in disease occurrence and prevalence can be one signal of climate-stressed wildlife populations. We conducted a serosurvey for vector-borne pathogens in snowshoe hares that were captured in the Hiawatha National Forest in the eastern UP of Michigan, USA, 2016-2017. The most commonly detected antibody response was to the mosquito-borne California serogroup snowshoe hare virus (SSHV). Overall, 24 (51%) hares screened positive for SSHV antibodies and of these, 23 (96%) were confirmed positive by plaque reduction neutralization test. We found a positive association between seroprevalence of SSHV and live weight of snowshoe hares. Additionally, we detected a significant effect of ecological land type group on seroprevalence of SSHV, with strong positive support for a group representing areas that tend to support high numbers of hares (i.e., acidic mineral containing soils with cedar, mixed swamp conifers, tamarack and balsam fir as common overstory vegetation). We also detected and confirmed antibodies for Jamestown Canyon virus and Silverwater virus in a single hare each. We did not detect antibodies to other zoonotic vector-borne pathogens, including Lacrosse encephalitis virus, West Nile virus, Borrelia burgdorferi, Powassan virus, and Francisella tularensis. These results provide a baseline for future serological studies of vector-transmitted diseases that may increase climate vulnerability of snowshoe hares in the UP of Michigan, as well as pose a climate-related zoonotic risk.

Use of hare bone for the manufacture of a Clovis bead.

A tubular bone bead dating to ~ 12,940 BP was recovered from a hearth-centered activity area at the La Prele Mammoth site in Converse County, Wyoming, USA. This is the oldest known bead from the Western Hemisphere. To determine the taxonomic origin of the bead, we extracted collagen for zooarchaeology by mass spectrometry (ZooMS). We also used micro-CT scanning for morphological analysis to determine likely skeletal elements used for its production. We conclude that the bead was made from a metapodial or proximal phalanx of a hare (Lepus sp.). This find represents the first secure evidence for the use of hares during the Clovis period. While the use of hare bone for the manufacture of beads was a common practice in western North America during the Holocene, its origins can now be traced back to at least the terminal Pleistocene.

Species-specific variation in mitochondrial genome tandem repeat polymorphisms in hares (Lepus spp., Lagomorpha, Leporidae) provides insight into their evolution.

The non-coding regions of the mitochondrial DNAs (mtDNAs) of hares, rabbits, and pikas (Lagomorpha) contain short (∼20 bp) and long (130-160 bp) tandem repeats, absent in related mammalian orders. In the presented study, we provide in-depth analysis for mountain hare (Lepus timidus) and brown hare (L. europaeus) mtDNA non-coding regions, together with a species- and population-level analysis of tandem repeat variation. Mountain hare short tandem repeats (SRs) as well as other analyzed hare species consist of two conserved 10 bp motifs, with only brown hares exhibiting a single, more variable motif. Long tandem repeats (LRs) also differ in sequence and copy number between species. Mountain hares have four to seven LRs, median value five, while brown hares exhibit five to nine LRs, median value six. Interestingly, introgressed mountain hare mtDNA in brown hares obtained an intermediate LR length distribution, with median copy number being the same as with conspecific brown hare mtDNA. In contrast, transfer of brown hare mtDNA into cultured mtDNA-less mountain hare cells maintained the original LR number, whereas the reciprocal transfer caused copy number instability, suggesting that cellular environment rather than the nuclear genomic background plays a role in the LR maintenance. Due to their dynamic nature and separation from other known conserved sequence elements on the non-coding region of hare mitochondrial genomes, the tandem repeat elements likely to represent signatures of ancient genetic rearrangements. clarifying the nature and dynamics of these rearrangements may shed light on the possible role of NCR repeated elements in mitochondria and in species evolution.

The hare syphilis agent is related to, but distinct from, the treponeme causing rabbit syphilis.

The treponemes infecting lagomorphs include Treponema paraluisleporidarum ecovar Cuniculus (TPeC) and ecovar Lepus (TPeL), infecting rabbits and hares, respectively. In this study, we described the first complete genome sequence of TPeL, isolate V3603-13, from an infected mountain hare (Lepus timidus) in Sweden. In addition, we determined 99.0% of the genome sequence of isolate V246-08 (also from an infected mountain hare, Sweden) and 31.7% of the genome sequence of isolate Z27 A77/78 (from a European hare, Lepus europeaus, The Netherlands). The TPeL V3603-13 genome had considerable gene synteny with the TPeC Cuniculi A genome and with the human pathogen T. pallidum, which causes syphilis (ssp. pallidum, TPA), yaws (ssp. pertenue, TPE) and endemic syphilis (ssp. endemicum, TEN). Compared to the TPeC Cuniculi A genome, TPeL V3603-13 contained four insertions and 11 deletions longer than three nucleotides (ranging between 6 and2,932 nts). In addition, there were 25 additional indels, from one to three nucleotides long, altogether spanning 36 nts. The number of single nucleotide variants (SNVs) between TPeC Cuniculi A and TPeL V3603-13 were represented by 309 nucleotide differences. Major proteome coding differences between TPeL and TPeC were found in the tpr gene family, and (predicted) genes coding for outer membrane proteins, suggesting that these components are essential for host adaptation in lagomorph syphilis. The phylogeny revealed that the TPeL sample from the European brown hare was more distantly related to TPeC Cuniculi A than V3603-13 and V246-08.

The low survival rate of European hare leverets in arable farmland: evidence from the predation experiment.

The low survival rate of leverets may significantly contribute to steep population declines and slow recovery of European hares (Lepus europaeus). However, the leveret survival rate in farmlands with different landscape structures is poorly understood, and the existing evidence comes mainly from Western Europe. In this study, we explored the survival of leveret hare dummies along linear semi-natural habitats in homogeneous Central European arable farmland during the main part of the European hare reproduction period (March-April) in 2019 and 2020. The survival rate of hare leverets during the 14-day period was only 22.2%, and all predation events were recorded during the first six days of the experiment. Mammalian predators were responsible for 53.1% of predation events, avian predators for 40.8%, and agricultural operations for 6.1%. The red fox (Vulpes vulpes) was the dominant predator in our study area and was the primary cause of leveret dummy mortality (32.7%), but it also had the highest use-intensity and visit frequency of all of the study plots. Predation by avian predators was associated with patches of lower vegetation height and cover (such as plowed fields) and during daylight hours, whereas the opposite was true for mammalian predators. We propose that improving the habitat quality of arable landscapes by increasing the proportion and quality of extensively used non-farmed habitats (e.g., set-asides, wildflower areas, extensive meadows, fallow land, and semi-natural habitats on arable land) providing cover and shelter for leverets could be an effective management measure for reducing predation risk on leverets.

Identification and Characterisation of a Myxoma Virus Detected in the Italian Hare (Lepus corsicanus).

Myxoma virus (MYXV) is a Leporipoxvirus (genus) belonging to the family Poxviridae; it is characterised by a genome of approximately 161 kb dsDNA encoding for several proteins that play an essential role in both host spectrum determination and immunomodulation. The healthy reservoir of the virus is Sylvilagus spp. At the same time, in wild and domestic European rabbits (Oryctolagus cuniculus), MYXV is the etiologic agent of myxomatosis, a disease with an extremely high mortality rate. In 2014, an interspecies jump of MYXV was reported in Lepus europaeus in the UK. In 2018, myxomatosis induced by a new recombinant strain called MYXV-To was identified during a large outbreak in Iberian hares (Lepus granatensis) in Spain. Here, we describe the case of myxomatosis in another hare species: an adult male Italian hare (Lepus corsicanus) found dead in 2018 in Sicily with lesions suggestive of myxomatosis and treponema infection. Laboratory tests, e.g., end-point PCR and negative staining electron microscopy, confirmed the presence of both pathogens. MYXV was then isolated from tissue samples in permissive cells and sequenced using NGS technology. Main genomic differences concerning known MYXV strains are discussed.

A new HaCV-EBHSV recombinant lagovirus circulating in European brown hares (Lepus europaeus) from Catalonia, Spain.

In 2020/2021, several European brown hare syndrome virus (EBHSV) outbreaks were recorded in European hares (Lepus europaeus) from Catalonia, Spain. Recombination analysis combined with phylogenetic reconstruction and estimation of genetic distances of the complete coding sequences revealed that 5 strains were recombinants. The recombination breakpoint is located within the non-structural protein 2C-like RNA helicase (nucleotide position ~ 1889). For the genomic fragment upstream of the breakpoint, a non-pathogenic EBHSV-related strain (hare calicivirus, HaCV; GII.2) was the most closely related sequence; for the rest of the genome, the most similar strains were the European brown hare syndrome virus (EBHSV) strains recovered from the same 2020/2021 outbreaks, suggesting a recent origin. While the functional impact of the atypical recombination breakpoint remains undetermined, the novel recombinant strain was detected in different European brown hare populations from Catalonia, located 20-100 km apart, and seems to have caused a fatal disease both in juvenile and adult animals, confirming its viability and ability to spread and establish infection. This is the first report of a recombination event involving HaCV and EBHSV and, despite the recombination with a non-pathogenic strain, it appears to be associated with mortality in European brown hares, which warrants close monitoring.

Gut diazotrophs in lagomorphs are associated with season but not altitude and host phylogeny.

Invertebrates such as termites feeding on nutrient-poor substrate receive essential nitrogen by biological nitrogen fixation of gut diazotrophs. However, the diversity and composition of gut diazotrophs of vertebrates such as Plateau pikas living in nutrient-poor Qinghai-Tibet Plateau remain unknown. To fill this knowledge gap, we studied gut diazotrophs of Plateau pikas (Ochotona curzoniae) and its related species, Daurian pikas (Ochotona daurica), Hares (Lepus europaeus) and Rabbits (Oryctolagus cuniculus) by high-throughput amplicon sequencing methods. We analyzed whether the gut diazotrophs of Plateau pikas are affected by season, altitude, and species, and explored the relationship between gut diazotrophs and whole gut microbiomes. Our study showed that Firmicutes, Spirochaetes, and Euryarchaeota were the dominant gut diazotrophs of Plateau pikas. The beta diversity of gut diazotrophs of Plateau pikas was significantly different from the other three lagomorphs, but the alpha diversity did not show a significant difference among the four lagomorphs. The gut diazotrophs of Plateau pikas were the most similarly to that of Rabbits, followed by Daurian pikas and Hares, which was inconsistent with gut microbiomes or animal phylogeny. The dominant gut diazotrophs of the four lagomorphs may reflect their living environment and dietary habits. Season significantly affected the alpha diversity and abundance of dominant gut diazotrophs. Altitude had no significant effect on the gut diazotrophs of Plateau pikas. In addition, the congruence between gut microbiomes and gut diazotrophs was low. Our results proved that the gut of Plateau pikas was rich in gut diazotrophs, which is of great significance for the study of ecology and evolution of lagomorphs.

Occurrence of Coxiella burnetii in wild lagomorphs and their ticks in Spanish Mediterranean ecosystems.

BACKGROUND: Coxiella burnetii, the causative agent of Q fever, is a zoonotic multi-host vector-borne pathogen of major public health importance. Although the European Food Safety Authority has recently made the monitoring of this bacterium in wildlife a priority, the role of wild lagomorphs in the transmission and maintenance of C. burnetii is poorly understood. AIMS: The aims of this study were to determine the prevalence and risk factors associated with C. burnetii circulation in European wild rabbits (Oryctolagus cuniculus) and Iberian hares (Lepus granatensis) and to assess the presence of this pathogen in ticks that feed on them in Mediterranean ecosystems in Spain, the country with the highest number of reported cases of Q fever in Europe. METHODS: A total of 574 spleen samples were collected from 453 wild rabbits and 121 Iberian hares, and 513 ticks (processed in 120 pools) between the 2017/2018 and 2021/2022 hunting seasons. RESULTS: C. burnetii DNA was detected in 103 (17.9%; 95% CI: 14.8-21.1) of the 574 wild lagomorphs tested. By species, prevalence was 16.3% (74/453; 95% CI: 12.9-19.7) in the European wild rabbit and 24.0% (29/121; 95% CI: 16.4-31.6) in the Iberian hare. At least one positive lagomorph was found on 47.9% of the 96 hunting estates sampled and in every hunting season since 2018/2019. Two risk factors associated with C. burnetii infection were as follows: outbreak of myxomatosis on the hunting estate in the month prior to sampling and high tick abundance observed by gamekeepers on the hunting estate. C. burnetii DNA was also found in 33 of the 120 (27.5%; 95% CI: 19.5-35.5) tick pools tested. The pathogen was detected in 66.7% (4/6), 29.2% (26/89) and 21.4% (3/14) of Haemaphysalis hispanica, Rhipicephalus pusillus and Hyalomma lusitanicum pools respectively. CONCLUSIONS: This study provides new epidemiological data on C. burnetii in European wild rabbits and is the first survey on this zoonotic pathogen performed in Iberian hares. Our results indicate widespread endemic circulation of C. burnetii and highlight the importance of both wild lagomorph species as natural reservoirs of this zoonotic bacterium in Mediterranean ecosystems in southern Spain, which may be of public and animal health concern. The high prevalence and wide diversity of positive tick species suggest the possible role of ticks in the epidemiological cycle of C. burnetii, with the potential risk of transmission to sympatric species, including humans.