Dynamics of Humoral Immunity to Myxoma and Rabbit Hemorrhagic Disease Viruses in Wild European Rabbits Assessed by Longitudinal Semiquantitative Serology

Myxoma virus (MYXV) and rabbit hemorrhagic disease virus (RHDV) are important drivers of the population decline of the European rabbit, an endangered keystone species. Both viruses elicit strong immune responses, but the long-term dynamics of humoral immunity are imperfectly known. This study aimed to assess the determinants of the long-term dynamics of antibodies to each virus based on a longitudinal capture-mark-recapture of wild European rabbits and semiquantitative serological data of MYXV and RHDV GI.2-specific IgG. The study included 611 indirect enzyme-linked immunosorbent assay (iELISA) normalized absorbance ratios for each MYXV and RHDV GI.2 from 505 rabbits from 2018 to 2022. Normalized absorbance ratios were analyzed using log-linear mixed models, showing a significant positive relationship with the time since the first capture of individual rabbits, with monthly increases of 4.1% for antibodies against MYXV and 2.0% against RHDV GI.2. Individual serological histories showed fluctuations over time, suggesting that reinfections boosted the immune response and likely resulted in lifelong immunity. Normalized absorbance ratios significantly increased with the seroprevalence in the population, probably because of recent outbreaks, and with body weight, highlighting the role of MYXV and RHDV GI.2 in determining survival to adulthood. Juvenile rabbits seropositive for both viruses were found, and the dynamics of RHDV GI.2 normalized absorbance ratios suggest the presence of maternal immunity up to 2 months of age. Semiquantitative longitudinal serological data provide epidemiological information, otherwise lost when considering only qualitative data, and support a lifelong acquired humoral immunity to RHDV GI.2 and MYXV upon natural infection. IMPORTANCE This study addresses the long-term dynamics of humoral immunity to two major viral pathogens of the European rabbit, an endangered keystone species of major ecological relevance. Such studies are particularly challenging in free-ranging species, and a combination of longitudinal capture-mark-recapture and semiquantitative serology was used to address this question. Over 600 normalized absorbance ratios of iELISA, obtained from 505 individual rabbits in 7 populations over 5 years, were analyzed using linear mixed models. The results support a lifelong acquired humoral immunity to myxoma virus and rabbit hemorrhagic disease virus upon natural infection and suggest the presence of maternal immunity to the latter in wild juvenile rabbits. These results contribute to understanding the epidemiology of two viral diseases threatening this keystone species and assist in developing conservation programs.

Keeping an eye on the use of eye-lens weight as a universal indicator of age for European wild rabbits

Accurate methods for age determination are critical to the knowledge of wildlife populations' age structure and, therefore, to their successful management. The reliability of age estimation may have profound economic and ecological consequences on the management of the European wild rabbits, Oryctolagus cuniculus, in its native and introduced range, where it is a keystone species and a major pest, respectively. As in other mammal species, European rabbits' age is often estimated using the Gompertz relationship between age and lens' weight. The growth rate formula has been developed based on data collected from European rabbits introduced in Australia, where a single subspecies (O. cuniculus cuniculus, Occ) is present. However, this curve has never been validated in the species native range, the Iberian Peninsula, where two subspecies (Occ, and O. c. algirus, Oca) coexist naturally. In this study, we tested the relationship between age and lens' weight using 173 Occ and 112 Oca wild rabbits that were surveyed in two experimental facilities in Spain. Our findings show that, in the native range, the published growth curve formula fits well Occ but not Oca data. Therefore, we recommend using the formula reported in this study to estimate the age of Oca (Lens dry weight = 240 × 10^{(-64.9/(Age+32))}). This study supports Oca rabbits' distinctiveness revealed by previous studies, which suggests that management interventions should be applied to protect this subspecies whose distribution range is very narrow and whose populations seem to be declining. More broadly, our findings point to the importance of testing the suitability of growth curves defined for other species with different genetic forms as occurs in the European wild rabbit case.