Surveillance design for hare-associated myxoma virus in Europe: arthropod evidence, seasonal concordance, and detection power

# MYX3 figures and figure source data Figures and the numerical source data behind them, for the manuscript: **Surveillance design for hare-associated myxoma virus in Europe: arthropod evidence, seasonal concordance, and detection power** Štěpán Ryba (ORCID 0009-0005-7126-6409)Department of Animal Husbandry Sciences, Faculty of Agriculture and Technology, University of South Bohemia in České Budějovice, Studentská 1668, 370 05 České Budějovice, Czech Republic. Contact: [e-mail] ## Scope of this deposit This deposit contains the five published figures as vector PDF and one tab-separated file per figure panel set. Each source-data file holds the values that the corresponding figure draws. The files are sufficient to redraw every figure and to check every plotted value against the manuscript text. The analysis code, input data, model configuration files, software manifest and claim ledger are deposited separately and are cited in the manuscript. Raster and editable exports (SVG, TIFF, JPG) are not included. The PDF files are vector and are the authoritative version of each figure. ## Contents ### figures/ | File | Figure | Size (mm) | Content ||---|---|---|---|| F1v1_vector_evidence_architecture.pdf | Fig. 1 | 193.9 x 145.2 | Architecture of the curated MYXV arthropod evidence set, 49 evidence units, four panels || F2v1_host_vector_transferability.pdf | Fig. 2 | 180.0 x 127.6 | Schematic of the rabbit-flea mechanism, transferability to hares, and the dipteran contact route || F3v7_phenology_concordance.pdf | Fig. 3 | 170.2 x 219.6 | ha-MYXV onset timing and phenological concordance of candidate dipterans || F4v1_negative_surveillance_power.pdf | Fig. 4 | 172.9 x 204.2 | Detection power and limits of negative arthropod surveillance || F5v1_suitability_frontier.pdf | Fig. 5 | 381.7 x 223.9 | Environmental suitability for the three focal arthropods and the documented 2025 front | Figure 2 is a schematic and has no numerical source data. Its panel content follows the host-vector mechanism matrix described in the manuscript Materials and Methods. The animal illustrations in Figure 2 are original assets produced with an image-generation model and are visual depiction only. They carry no data. ### source_data/ **F1v1_vector_evidence_architecture_source_data.tsv** (38 records) Columns: `metric_key`, `statement`, `value`, `unit`, `used_in`. One record per plotted or stated quantity, with the sentence it supports. **F3v7_phenology_concordance_source_data.tsv** (1136 records) Columns: `panel`, `series`, `x_type`, `x`, `y_type`, `y`, `lower`, `upper`, `onset_id`, `country_iso2`, `note`. Panel A holds the standardized daily activity percentile for each taxon. Panel B holds the observed onset interval of each region-year onset unit. Panel C holds the mean dipteran-minus-flea activity-percentile difference at lags of 0, 14 and 28 days, with the event-level bootstrap interval and the country-year cluster bootstrap interval at lag 0. **F4v1_negative_surveillance_power_source_data.tsv** (1905 records) Columns: `panel`, `scenario`, `design_prevalence`, `target_probability`, `n`, `detection_probability`, `assay_sensitivity`, `clustering`, `note`. Design prevalence is a planning quantity. It is not an estimate of field positivity in arthropods. **F5v1_suitability_frontier_event_source_data.tsv** (77 records) Columns: `event_id`, `country`, `event_date_start`, `year`, `latitude`, `longitude`, `x3035`, `y3035`, `documented_2025_front_point`, `front_geometry_interpolated`. Coordinates are given in WGS84 and in EPSG:3035. The column `front_geometry_interpolated` is `no` for every record. The 2025 front is represented by 13 observed points. No line, hull or interpolated corridor was constructed. **F5v1_suitability_frontier_summary_source_data.tsv** (3 records, one per taxon) Columns include `primary_q10_threshold`, `M2_novel_environment_fraction`, `documented_2025_front_equal_country_fraction_ge_q10`, `front_adjacent_100km_equal_country_area_fraction_ge_q10`, the analysis output path behind each value, and an interpretation field. Country values are averaged with equal weight, so that Hungary, which contributes 10 of the 13 front points, does not dominate the summary. ## How to read the suitability values Suitability is the cloglog output of a MaxEnt-type model fitted to occurrence records and bioclimatic predictors. It describes environmental conditions for the arthropod. It is not abundance, vectorial capacity, transmission probability or disease risk. Cloglog values are model-specific, so absolute suitability is not comparable among taxa. Each taxon is assessed against its own 10th-percentile training-presence threshold, given as `primary_q10_threshold`. ## Third-party data Occurrence records come from three DOI-bearing GBIF downloads and are not redistributed here: https://doi.org/10.15468/dl.zcr7fb, https://doi.org/10.15468/dl.ja6ke4 and https://doi.org/10.15468/dl.8ef2rx. Environmental predictors come from CHELSA-bioclim v2.1, https://doi.org/10.16904/envidat.228. Administrative and coastline reference layers come from Natural Earth 1:10m, version 5.1.1. The digitized source figures behind the three phenology reference series are third-party material and are not redistributed. The numerical series derived from them are included in the Figure 3 source data. ## File formats PDF 1.4 vector, generated with Matplotlib 3.11.1. Tab-separated UTF-8 text with one header row. The decimal separator is a point. Empty cells are empty and are not coded. ## Licence CC BY 4.0. ## Checksums SHA-256 checksums for every file are in SHA256SUMS.txt.