Impact of myxomatosis in relation to local persistence in wild rabbit populations: the role of waning immunity and the reproductive period

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.

Spatiotemporal monitoring of myxomatosis in European wild rabbit (Oryctolagus cuniculus) in Spanish Mediterranean ecosystems

A long-term active epidemiological surveillance programme was conducted to determine seroprevalence to myxoma virus (MYXV), infection prevalence and spatiotemporal patterns and factors associated with MYXV circulation in wild rabbits (Oryctolagus cuniculus) in Spanish Mediterranean ecosystems. A total of 2376 animals were sampled over four study periods: 2009-2012 (P1), 2012-2015 (P2), 2015-2018 (P3) and 2018-2021 (P4). Antibodies against MYXV were detected by a commercial indirect ELISA in 59.9% (1424/2376; 95% CI: 58.0-61.9) of wild rabbits. At least one seropositive animal was detected on 131 (96.3%) of 136 game estates sampled. MYXV infection was confirmed by PCR in 94 of 1063 (8.8%; 95% CI: 7.3-10.7) wild rabbits. Circulation of the novel recombinant MYXV (ha-MYXV) was not found in wild rabbits analysed during P4. Five statistically significant spatiotemporal clusters of high MYXV seroprevalence were identified using a Bernoulli model: one in P2 and four in P3. A generalized linear mixed model (GLMM) analysis identified sampling season (autumn), age (adult and juvenile), outbreaks of myxomatosis in the month prior to sampling, mean annual temperature, humidity and seropositivity to rabbit haemorrhagic disease virus as factors potentially linked with MYXV seropositivity. GLMM analysis identified outbreaks of myxomatosis in the month prior to sampling, MYXV seropositivity and presence of lesions compatible with myxomatosis as factors associated with MYXV infection. The results indicate high exposure, widespread but non-homogeneous distribution, and endemic circulation of MYXV in wild rabbit populations in southern Spain during the last decade. Prevalence of antibodies against MYXV showed fluctuations both within the year and over the study periods, revealing variations in the immunity of wild rabbit populations in Mediterranean ecosystems that could increase the risk of MYXV re-emergence in immunologically naïve populations. The present study highlights the importance of long-term surveillance to better understand the epidemiology of MYXV in wild lagomorphs.

Monitoring of emerging myxoma virus epidemics in Iberian hares (Lepus granatensis) in Spain, 2018-2020

Myxomatosis is an infectious disease caused by the myxoma virus (MYXV), which has very high mortality rates in European wild rabbits (Oryctolagus cuniculus). While sporadic cases of myxomatosis have also been reported in some hare species, these lagomorphs are considered to have a low susceptibility to MYXV infection. In the present study, we describe the spatiotemporal evolution and main epidemiological findings of novel hare MYXV (ha-MYXV or MYXV-Tol) epidemics in Iberian hares (Lepus granatensis) in Spain. In the period 2018-2020, a total of 487 hares from 372 affected areas were confirmed to be MYXV-infected by PCR. ha-MYXV outbreaks were detected in most of the Spanish regions where the Iberian hare is present. The spatial distribution was not homogeneous, with most outbreaks concentrated in the southern and central parts of Spain. Consecutive outbreaks reported in the last two years suggest endemic circulation in Spain of this emerging virus. A retrospective study carried out just after the first epidemic period (2018-2019) revealed that the virus could have been circulating since June 2018. The number of outbreaks started to rise in July, peaked during the first half of August and October and then decreased sharply until January 2019. The apparent mean mortality rate was 55.4% (median: 70%). The results indicated high susceptibility of the Iberian hare to ha-MYXV infection, but apparent resistance in the sympatric hare species present in Spain and less infectivity in European rabbits. The novel ha-MYXV has had significant consequences on the health status of Iberian hare populations in Spain, which is of animal health and conservation concern. The present study contributes to a better understanding of ha-MYXV emergence and will provide valuable information for the development of control strategies. Further research is warranted to assess the impact of this emerging virus on wild lagomorph populations and to elucidate its ecological implications for Iberian Mediterranean ecosystems.

Dynamic, spatial models of parasite transmission in wildlife: Their structure, applications and remaining challenges

Individual differences in contact rate can arise from host, group and landscape heterogeneity and can result in different patterns of spatial spread for diseases in wildlife populations with concomitant implications for disease control in wildlife of conservation concern, livestock and humans. While dynamic disease models can provide a better understanding of the drivers of spatial spread, the effects of landscape heterogeneity have only been modelled in a few well-studied wildlife systems such as rabies and bovine tuberculosis. Such spatial models tend to be either purely theoretical with intrinsic limiting assumptions or individual-based models that are often highly species- and system-specific, limiting the breadth of their utility. Our goal was to review studies that have utilized dynamic, spatial models to answer questions about pathogen transmission in wildlife and identify key gaps in the literature. We begin by providing an overview of the main types of dynamic, spatial models (e.g., metapopulation, network, lattice, cellular automata, individual-based and continuous-space) and their relation to each other. We investigate different types of ecological questions that these models have been used to explore: pathogen invasion dynamics and range expansion, spatial heterogeneity and pathogen persistence, the implications of management and intervention strategies and the role of evolution in host-pathogen dynamics. We reviewed 168 studies that consider pathogen transmission in free-ranging wildlife and classify them by the model type employed, the focal host-pathogen system, and their overall research themes and motivation. We observed a significant focus on mammalian hosts, a few well-studied or purely theoretical pathogen systems, and a lack of studies occurring at the wildlife-public health or wildlife-livestock interfaces. Finally, we discuss challenges and future directions in the context of unprecedented human-mediated environmental change. Spatial models may provide new insights into understanding, for example, how global warming and habitat disturbance contribute to disease maintenance and emergence. Moving forward, better integration of dynamic, spatial disease models with approaches from movement ecology, landscape genetics/genomics and ecoimmunology may provide new avenues for investigation and aid in the control of zoonotic and emerging infectious diseases.

Effects of myxoma virus and rabbit hemorrhagic disease virus on the physiological condition of wild European rabbits: Is blood biochemistry a useful monitoring tool?

Myxomatosis and rabbit hemorrhagic disease (RHD) are the major viral diseases that affect the wild European rabbit (Oryctolagus cuniculus). These diseases arrived in Europe within the last decades and have caused wild rabbit populations to decline dramatically. Both viruses are currently considered to be endemic in the Iberian Peninsula; periodic outbreaks that strongly impact wild populations regularly occur. Myxoma virus (MV) and rabbit hemorrhagic disease virus (RHDV) alter the physiology of infected rabbits, resulting in physical deterioration. Consequently, the persistence and viability of natural populations are affected. The main goal of our study was to determine if blood biochemistry is correlated with serostatus in wild European rabbits. We carried out seven live-trapping sessions in three wild rabbit populations over a two-year period. Blood samples were collected to measure anti-MV and anti-RHDV antibody concentrations and to measure biochemical parameters related to organ function, protein metabolism, and nutritional status. Overall, we found no significant relationships between rabbit serostatus and biochemistry. Our main result was that rabbits that were seropositive for both MV and RHDV had low gamma glutamyltransferase concentrations. Given the robustness of our analyses, the lack of significant relationships may indicate that the biochemical parameters measured are poor proxies for serostatus. Another explanation is that wild rabbits might be producing attenuated physiological responses to these viruses because the latter are now enzootic in the study area.

Distribution and prevalence of the Australian non-pathogenic rabbit calicivirus is correlated with rainfall and temperature

BACKGROUND

Australia relies heavily on rabbit haemorrhagic disease virus (RHDV) for the biological control of introduced European wild rabbits Oryctolagus cuniculus, which are significant economic and environmental pests. An endemic non-pathogenic rabbit calicivirus termed RCV-A1 also occurs in wild rabbits in Australian and provides partial protection against lethal RHDV infection, thus interfering with effective rabbit control. Despite its obvious importance for rabbit population management, little is known about the epidemiology of this benign rabbit calicivirus.

METHODS

We determined the continent-wide distribution and prevalence of RCV-A1 by analysing 1,805 serum samples from wild rabbit populations at 78 sites across Australia for the presence of antibodies to RCV-A1 using a serological test that specifically detects RCV-A1 antibodies and does not cross-react with co-occurring RHDV antibodies. We also investigated possible correlation between climate variables and prevalence of RCV-A1 by using generalised linear mixed effect models.

RESULTS

Antibodies to RCV-A1 were predominantly detected in rabbit populations in cool, high rainfall areas of the south-east and south-west of the continent. There was strong support for modelling RCV-A1 prevalence as a function of average annual rainfall and minimum temperature. The best ranked model explained 26% of the model structural deviance. According to this model, distribution and prevalence of RCV-A1 is positively correlated with periods of above average rainfall and negatively correlated with periods of drought.

IMPLICATIONS

Our statistical model of RCV-A1 prevalence will greatly increase our understanding of RCV-A1 epidemiology and its interaction with RHDV in Australia. By defining the environmental conditions associated with the prevalence of RCV-A1, it also contributes towards understanding the distribution of similar viruses in New Zealand and Europe.

Multi-event capture-recapture modeling of host-pathogen dynamics among European rabbit populations exposed to myxoma and Rabbit Hemorrhagic Disease Viruses: common and heterogeneous patterns

Host-pathogen epidemiological processes are often unclear due both to their complexity and over-simplistic approaches used to quantify them. We applied a multi-event capture-recapture procedure on two years of data from three rabbit populations to test hypotheses about the effects on survival of, and the dynamics of host immunity to, both myxoma virus and Rabbit Hemorrhagic Disease Virus (MV and RHDV). Although the populations shared the same climatic and management conditions, MV and RHDV dynamics varied greatly among them; MV and RHDV seroprevalences were positively related to density in one population, but RHDV seroprevalence was negatively related to density in another. In addition, (i) juvenile survival was most often negatively related to seropositivity, (ii) RHDV seropositives never had considerably higher survival, and (iii) seroconversion to seropositivity was more likely than the reverse. We suggest seropositivity affects survival depending on trade-offs among antibody protection, immunosuppression and virus lethality. Negative effects of seropositivity might be greater on juveniles due to their immature immune system. Also, while RHDV directly affects survival through the hemorrhagic syndrome, MV lack of direct lethal effects means that interactions influencing survival are likely to be more complex. Multi-event modeling allowed us to quantify patterns of host-pathogen dynamics otherwise difficult to discern. Such an approach offers a promising tool to shed light on causative mechanisms.

Early infections by myxoma virus of young rabbits (Oryctolagus cuniculus) protected by maternal antibodies activate their immune system and enhance herd immunity in wild populations

The role of maternal antibodies is to protect newborns against acute early infection by pathogens. This can be achieved either by preventing any infection or by allowing attenuated infections associated with activation of the immune system, the two strategies being based on different cost/benefit ratios. We carried out an epidemiological survey of myxomatosis, which is a highly lethal infectious disease, in two distant wild populations of rabbits to describe the epidemiological pattern of the disease. Detection of specific IgM and IgG enabled us to describe the pattern of immunity. We show that maternal immunity attenuates early infection of juveniles and enables activation of their immune system. This mechanism associated with steady circulation of the myxoma virus in both populations, which induces frequent reinfections of immune rabbits, leads to the maintenance of high immunity levels within populations. Thus, myxomatosis has a low impact, with most infections being asymptomatic. This work shows that infection of young rabbits protected by maternal antibodies induces attenuated disease and activates their immune system. This may play a major role in reducing the impact of a highly lethal disease when ecological conditions enable permanent circulation of the pathogen.

Length of intervals between epidemics: evaluating the influence of maternal transfer of immunity

The length of intervals between epidemic outbreaks of infectious diseases is critical in epidemiology. In several species of marine mammals and birds, it is pivotal to also consider the life history of the species of concern, as the contact rate between individuals can have a seasonal flux, for example, due to aggregations during the breeding season. Recently, particular interest has been given to the role of the dynamics of immunity in determining the intervals between epidemics in wild animal populations. One potentially powerful, but often neglected, process in this context is the maternal transfer of immunity. Here, we explore theoretically how the transfer of maternal antibodies can delay the recurrence of epidemics using Phocine Distemper in harbor seals as an example of a system in which epidemic outbreaks are followed by pathogen extinction. We show that the presence of temporarily protected newborns can significantly increase the predicted interval between epidemics, and this effect is strongly dependent on the degree of synchrony in the breeding season. Furthermore, we found that stochasticity in the onset of epidemics in combination with maternally acquired immunity increases the predicted intervals between epidemics even more. These effects arise because newborns with maternal antibodies temporarily boost population level immunity above the threshold of herd immunity, particularly when breeding is synchronous. Overall, our results show that maternal antibodies can have a profound influence on the dynamics of wildlife epidemics, notably in gregarious species such as many marine mammals and seabirds.

Histo-blood group antigens act as attachment factors of rabbit hemorrhagic disease virus infection in a virus strain-dependent manner

Rabbit Hemorrhagic disease virus (RHDV), a calicivirus of the Lagovirus genus, and responsible for rabbit hemorrhagic disease (RHD), kills rabbits between 48 to 72 hours post infection with mortality rates as high as 50–90%. Caliciviruses, including noroviruses and RHDV, have been shown to bind histo-blood group antigens (HBGA) and human non-secretor individuals lacking ABH antigens in epithelia have been found to be resistant to norovirus infection. RHDV virus-like particles have previously been shown to bind the H type 2 and A antigens. In this study we present a comprehensive assessment of the strain-specific binding patterns of different RHDV isolates to HBGAs. We characterized the HBGA expression in the duodenum of wild and domestic rabbits by mass spectrometry and relative quantification of A, B and H type 2 expression. A detailed binding analysis of a range of RHDV strains, to synthetic sugars and human red blood cells, as well as to rabbit duodenum, a likely gastrointestinal site for viral entrance was performed. Enzymatic cleavage of HBGA epitopes confirmed binding specificity. Binding was observed to blood group B, A and H type 2 epitopes in a strain-dependent manner with slight differences in specificity for A, B or H epitopes allowing RHDV strains to preferentially recognize different subgroups of animals. Strains related to the earliest described RHDV outbreak were not able to bind A, whereas all other genotypes have acquired A binding. In an experimental infection study, rabbits lacking the correct HBGA ligands were resistant to lethal RHDV infection at low challenge doses. Similarly, survivors of outbreaks in wild populations showed increased frequency of weak binding phenotypes, indicating selection for host resistance depending on the strain circulating in the population. HBGAs thus act as attachment factors facilitating infection, while their polymorphism of expression could contribute to generate genetic resistance to RHDV at the population level.

Rabbit hemorrhagic disease virus (RHDV), detected as late as 1984, has spread to large parts of the world, threatening rabbit populations and other species dependent on rabbits in many European countries. Mortality has been shown to be as high as 90% and rabbits are killed 48 to 72 hours after infection. Related viruses called noroviruses, infect humans in a manner dependent on the expression of histo-blood group antigens (HBGAs), which are not only expressed on red blood cells, but also on epithelial cells, in saliva and on mucins of the intestinal tract. RHDV also binds to HBGA and in this report we characterize binding of strains of all genetic groups of RHDV to different HBGAs. We also demonstrate HBGAs to function as attachment factors in a challenge experiment. As polymorphisms of genes involved in HBGA synthesis divide the rabbit population into different subgroups, we find selection of low-binding subgroups of wild rabbits in populations recovering from devastating outbreaks of RHDV. This is the first demonstration of differential HBGA specificities of RHDV strains, description of function in infection and demonstration of host selection due to RHDV infection based on HBGA phenotype.

Hantavirus infections in fluctuating host populations: the role of maternal antibodies

Infected females may transfer maternal antibodies (MatAbs) to their offspring, which may then be transiently protected against infections the mother has encountered. However, the role of maternal protection in infectious disease dynamics in wildlife has largely been neglected. Here, we investigate the effects of Puumala hantavirus (PUUV)-specific MatAbs on PUUV dynamics, using 7 years' data from a cyclic bank vole population in Finland. For the first time to our knowledge, we partition seropositivity data from a natural population into separate dynamic patterns for MatAbs and infection. The likelihood of young of the year carrying PUUV-specific MatAbs during the breeding season correlated positively with infection prevalence in the overwintered parent population in the preceding spring. The probability of PUUV infection varied between seasons (highest in spring, lowest in late summer) and depended on population structure, but was also, in late autumn, notably, negatively related to summer MatAb prevalence, as well as to infection prevalence earlier in the breeding season. Hence, our results suggest that high infection prevalence in the early breeding season leads to a high proportion of transiently immune young individuals, which causes delays in transmission. This suggests, in turn, that MatAb protection has the potential to affect infection dynamics in natural populations.

Quantitation of cytokine mRNA by real-time RT-PCR during a vaccination trial in a rabbit model of fascioliasis

Use of the rabbit as disease model has long been hampered by a lack of immunological assays specific to this species. In the present study we developed a SYBR Green-based, real-time RT-PCR protocol to quantitate cytokine mRNA in freshly harvested rabbit peripheral mononuclear cells. The method was validated in the course of a vaccination trial in which animals vaccinated with the recombinant antigen FhSAP2 were challenged with Fasciola hepatica metacercariae. Changes in the levels of rabbit interleukin (IL)-2, IL-4, IL-6, IL-10, tumor necrosis factor-alpha (TNFalpha), and interferon-gamma (IFNgamma) mRNA were determined. Messenger RNA from the universally expressed housekeeping gene GAPDH was used as an amplification control and allowed for correction of variations in the efficiencies of RNA extraction and reverse transcription. Rabbits vaccinated with FhSAP2 showed an 83.3% reduction in liver fluke burden after challenge infection when compared to non-vaccinated controls. All cytokine mRNAs were found at detectable levels; however, the levels of IFNgamma, TNFalpha, IL-2 and IL-10 were significantly higher in the vaccinated group compared to the non-vaccinated group. These results suggest that protection conferred by FhSAP2 protein could be associated with a mixed Th1/Th2 immune response in which Th1 cytokines are dominant. The real-time RT-PCR method described herein can be a useful tool for monitoring changes in basic immune functions in the rabbit model of fascioliasis and may also aid in studies of human diseases for which the rabbit is an important experimental model.