Prevalence of APV lesions was analyzed using General Linear Mixed Model with binomial distributed error and logit link function. Presence/absence of lesions was included in the model as the response variable, trios and site within each trio as random factors and habitat, sex, age, month and year as fixed independent factors. Only individuals with known age and sex were included in the analyses. Statistical analyses focus only on southern Spain data due to the larger dataset available and standardized sampling. Data from central Spain were used for comparison of prevalences and phylogenetic analyses. Models were fit using the GLIMMIX procedure in SAS 9.4.