Figure captions
Figure 1 . Van diagram with the iELISA and qPCR results for the 451 serum samples analysed.
Figure 2. Seroprevalence in Iberian hares by year (A) and season (B). Seroprevalence includes only results from hares randomly sampled (hunted or live captured), showing the proportion of the sample with antibodies recognizing myxoma virus antigens with 95% confidence interval.
Figure 3. Dot graph with the RI10 values from serum samples collected in different years. The left graph (A) refers to data from serum samples collected from hunted or live trapped hares. The right graph (B) refers to data from found dead hares’ serum samples. The two horizontal lines point out the IR value, the lower line representing the original cut-off of iELISA (IR=2.0) and upper line (IR 6.1) indicates the adjusted ELISA cutoff value. M- macroscopic lesions compatible with MYXV infection (no qPCR results available), N-negative in qPCR, P-positive in qPCR, ND-no data available on the qPCR test.
Figure 4. Sample size and seropositivity by administrative units in 1994-1999 and 2017-2019. A) Number of Iberian hares sampled; B) Administrative units with seropositive Iberian hares, considering the whole sample. Administrative units are districts in Portugal and provinces in Spain.
Figure 5. Hypotheses for Iberian hare susceptibility to MYXV and MYXV-related viruses. A-Hares are not susceptible  to infection by MYXV strains B-Hares are susceptible to MYXV if exposed to high infectious doses; C- ha-MYXV emerged in hares after a species jump event, involving or not, an intermediate host; D- Hares are infected by MYXV, but develop subclinical disease; E-ha-MYXV emerged in hares infected with rabbit strains, by a recombination event that increased virulence of the virus towards hares;  F-ha-MYXV emerged in hares from a strain that has circulated unnoticed in the population.