Interaction of myxomatosis and Rabbit Haemorrhagic Disease in wild rabbit

Increasing resistance of rabbits to myxomatosis in Australia has led to the exploration of Rabbit Haemorrhagic Disease, also called Rabbit Calicivirus Disease (RCD) as a possible control agent. While the initial spread of RCD in Australia resulted in widespread rabbit mortality in affected areas, the possible population dynamic effects of RCD and myxomatosis operating within the same system have not been properly explored. Here we present early mathematical modelling examining the interaction between the two diseases. In this study we use a deterministic compartment model, based on the classical SIR model in infectious disease modelling. We consider, here, only a single strain of myxomatosis and RCD and neglect latent periods. We also include logistic population growth, with the inclusion of seasonal birth rates. We assume there is no cross-immunity due to either disease. The mathematical model allows for the possibility of both diseases to be simultaneously present in an individual, although results are also presented for the case where co infection is not possible, since co-infection is thought to be rare and questions exist as to whether it can occur. The simulation results of this investigation show that it is a crucial issue and should be part of future field studies. A single simultaneous outbreak of RCD and myxomatosis was simulated, while ignoring natural births and deaths, appropriate for a short timescale of 20 days. Simultaneous outbreaks may be more common in Queensland. For the case where co-infection is not possible we find that the simultaneous presence of myxomatosis in the population suppresses the prevalence of RCD, compared to an outbreak of RCD with no outbreak of myxomatosis, and thus leads to a less effective control of the population. The reason for this is that infection with myxomatosis removes potentially susceptible rabbits from the possibility of infection with RCD (like a vaccination effect). We found that the reduction in the maximum prevalence of RCD was approximately 30% for an initial prevalence of 20% of myxomatosis, for the case where there was no simultaneous outbreak of myxomatosis, but the peak prevalence was only 15% when there was a simultaneous outbreak of myxomatosis. However, this maximum reduction will depend on other parameter values chosen. When co-infection is allowed then this suppression effect does occur but to a lesser degree. This is because the rabbits infected with both diseases reduces the prevalence of myxomatosis. We also simulated multiple outbreaks over a longer timescale of 10 years, including natural population growth rates, with seasonal birth rates and density dependent(logistic) death rates. This shows how both diseases interact with each other and with population growth. Here we obtain sustained outbreaks occurring approximately every two years for the case of a simultaneous outbreak of both diseases but without simultaneous co-infection, with the prevalence varying from 0.1 to 0.5. Without myxomatosis present then the simulation predicts RCD dies out quickly without further introduction from elsewhere. With the possibility of simultaneous co-infection of rabbits, sustained outbreaks are possible but then the outbreaks are less severe and more frequent (approximately yearly). While further model development is needed, our work to date suggests that: 1) the diseases are likely to interact via their impacts on rabbit abundance levels, and 2) introduction of RCD can suppress myxomatosis prevalence. We recommend that further modelling in conjunction with field studies be carried out to further investigate how these two diseases interact in the population.

Main ecological gradients and landscape matrix affect wild rabbit Oryctolagus cuniculus (Linnaeus, 1758) (Mammalia: Leporidae) abundance in a coastal region (South-East Spain)

Landscape analysis with transects, in the Marina Baja area (province of Alicante, Spain), has contributed to establish the influence of different landscape matrices and some environmental gradients on wild rabbit Oryctolagus cuniculus (Linnaeus, 1758) (Mammalia: Leporidae) abundance (kilometric abundance index, KAI). Transects (n = 396) were developed to estimate the abundance of this species in the study area from 2006 to 2008.Our analysis shows that rabbits have preferences for a specific land use matrix (irrigated: KAI = 3.47 ± 1.14 rabbits/km). They prefer the coastal area (KAI = 3.82 ± 1.71 rabbits/km), which coincides with thermo-Mediterranean (a bioclimatic belt with a tempered winter and a hot and dry summer with high human density), as opposed to areas in the interior (continental climate with lower human occupation). Their preference for the southern area of the region was also noted (KAI = 8.22 ± 3.90 rabbits/km), which coincides with the upper semi-arid area, as opposed to the northern and intermediate areas (the north of the region coinciding with the upper dry and the intermediate area with the lower dry). On the other hand, we found that the number of rabbits increased during the 3-year study period, with the highest abundance (KAI = 2.71 ± 1.30 rabbits/km) inMay. Thus, this study will enable more precise knowledge of the ecological factors (habitat variables) that intervene in the distribution of wild rabbit populations in a poorly studied area.

Distribution modelling of wild rabbit hunting yields in its original area (S Iberian Peninsula)

In this work we used the information of the Annual Hunting Reports (AHRs) to obtain a high-resolution model of the potential favourableness for wild rabbit harvesting in Andalusia (southern Spain), using environmental and land-use variables as predictors. We analysed 32,134 AHRs from the period 1993/2001 reported by 6049 game estates to estimate the average hunting yields of wild rabbit in each Andalusian municipality (n5771). We modelled the favourableness for obtaining good hunting yields using stepwise logistic regression on a set of climatic, orographical, land use, and vegetation variables. The favourability equation was used to create a downscaled image representing the favourableness of obtaining good hunting yields for the wild rabbit in 161 km squares in Andalusia, using the Idrisi Image Calculator. The variables that affected hunting yields of wild rabbit were altitude, dry wood crops (mainly olive groves, almond groves, and vineyards), temperature, pasture, slope, and annual number of frost days. The 161 km squares with high favourableness values are scattered throughout the territory, which seems to be caused mainly by the effect of vegetation. Finally, we obtained quality categories for the territory by combining the probability values given by logistic regression with those of the environmental favourability function.

Tracking the Origin of a Rabbit Haemorrhagic Virus 2 Outbreak in a Wild Rabbit Breeding Centre in Portugal; Epidemiological and Genetic Investigation

As key prey, the wild rabbit downsize constitutes a major drawback on the endangered Iberian lynx (Lynx pardinus) re-introduction in the Iberia. Several captive breeding units mostly located in Alentejo, endeavour the wild rabbit repopulation of depleted areas assigned for the lynx re-introduction. Here we report an RHDV2 outbreak that occurred in early 2016 in a wild rabbit captive breeding unit located in Barrancos municipality. The estimated mortality rate between March and April 2016 was approximately 8.67%. Anatomopathologic examination was carried out for 13 victimized rabbits. Molecular characterization was based on the complete vp60 capsid gene. The 13 rabbit carcasses investigated showed typical macroscopic RHD lesions testing positive to RHDV2-RNA. Comparison of the vp60 nucleotide sequences obtained from two specimens with others publically available disclosed similarities below 98.22% with RHDV2 strains originated in the Iberia and Azores and revealed that the two identical strains from Barrancos-2016 contain six unique single synonymous nucleotide polymorphisms. In the phylogenetic analysis performed, the Barrancos-2016 strains clustered apart from other known strains,meaning they may represent new evolutionary RHDV2 lineages. No clear epidemiological link could be traced for this outbreak where the mortalities were lower compared with previous years. Yet, network analysis suggested a possible connection between the missing intermediates from which the strains from Barrancos 2013, 2014 and 2016 have derived. It is therefore possible that RHDV2 has circulated endemically in the region since 2012, with periodic epizootic occurrences. Still, six years after its emergence in wild rabbits, RHDV2 continues to pose difficulties to the establishment of natural wild rabbit populations that are crucial for the self-sustainability of the local ecosystems.

Tracking the Origin of a Rabbit Haemorrhagic Virus 2 Outbreak in a Wild Rabbit Breeding Centre in Portugal; Epidemiological and Genetic Investigation

As key prey, the wild rabbit downsize constitutes a major drawback on the endangered Iberian lynx (Lynx pardinus) re-introduction in the Iberia. Several captive breeding units mostly located in Alentejo, endeavour the wild rabbit repopulation of depleted areas assigned for the lynx re-introduction. Here we report an RHDV2 outbreak that occurred in early 2016 in a wild rabbit captive breeding unit located in Barrancos municipality. The estimated mortality rate between March and April 2016 was approximately 8.67%. Anatomopathologic examination was carried out for 13 victimized rabbits. Molecular characterization was based on the complete vp60 capsid gene. The 13 rabbit carcasses investigated showed typical macroscopic RHD lesions testing positive to RHDV2- RNA. Comparison of the vp60 nucleotide sequences obtained from two specimens with others publically available disclosed similarities below 98.22% with RHDV2 strains originated in the Iberia and Azores and revealed that the two identical strains from Barrancos-2016 contain six unique single synonymous nucleotide polymorphisms. In the phylogenetic analysis performed, the Barrancos-2016 strains clustered apart from other known strains, meaning they may represent new evolutionary RHDV2 lineages. No clear epidemiological link could be traced for this outbreak where the mortalities were lower compared with previous years. Yet, network analysis suggested a possible connection between the missing intermediates from which the strains from Barrancos 2013, 2014 and 2016 have derived. It is therefore possible that RHDV2 has circulated endemically in the region since 2012, with periodic epizootic occurrences. Still, six years after its emergence in wild rabbits, RHDV2 continues to pose difficulties to the establishment of natural wild rabbit populations that are crucial for the self-sustainability of the local ecosystems.

Seropositivity to Rabbit Haemorrhagic Disease Virus in non-target mammals during periods of viral activity in a population of wild rabbits in New Zealand

As part of a longitudinal study of the epidemiology of rabbit haemorrhagic disease virus (RHDV) in New Zealand, serum samples were obtained from trapped feral animals that may have consumed European rabbit (Oryctolagus cuniculus) carcasses (non-target species). During a 21-month period when RHDV infection was monitored in a defined wild rabbit population, 16 feral house cats (Felis catus), 11 stoats (Mustela erminea), four ferrets (Mustela furo) and 126 hedgehogs (Erinaceus europaeus) were incidentally captured in the rabbit traps. The proportions of samples that were seropositive to RHDV were 38% for cats, 18% for stoats, 25% for ferrets and 4% for hedgehogs. Seropositive non-target species were trapped in April 2000, in the absence of an overt epidemic of rabbit haemorrhagic disease (RHD) in the rabbit population, but evidence of recent infection in rabbits was shown. Seropositive non-target species were found up to 2.5 months before and 1 month after this RHDV activity in wild rabbits was detected. Seropositive predators were also trapped on the site between 1 and 4.5 months after a dramatic RHD epidemic in February 2001. This study has shown that high antibody titres can be found in non-target species when there is no overt evidence of RHDV infection in the rabbit population, although a temporal relationship could not be assessed statistically owning to the small sample sizes. Predators and scavengers might be able to contribute to localised spread of RHDV through their movements.

The influence of habitat heterogeneity on patterns of connectivity among rabbit populations in southern Queensland

Patterns of connectivity among local populations influence the dynamics of regional systems, but most ecological models have concentrated on explaining the effect of connectivity on local population structure using dynamic processes covering short spatial and temporal scales. In this study, a model was developed in an extended spatial system to examine the hypothesis that long term connectivity levels among local populations are influenced by the spatial distribution of resources and other habitat factors. The habitat heterogeneity model was applied to local wild rabbit populations in the semi-arid Mitchell region of southern central Queensland (the Eastern system). Species' specific population parameters which were appropriate for the rabbit in this region were used. The model predicted a wide range of long term connectivity levels among sites, ranging from the extreme isolation of some sites to relatively high interaction probabilities for others. The validity of model assumptions was assessed by regressing model output against independent population genetic data, and explained over 80% of the variation in the highly structured genetic data set. Furthermore, the model was robust, explaining a significant proportion of the variation in the genetic data over a wide range of parameters. The performance of the habitat heterogeneity model was further assessed by simulating the widely reported recent range expansion of the wild rabbit into the Mitchell region from the adjacent, panmictic Western rabbit population system. The model explained well the independently determined genetic characteristics of the Eastern system at different hierarchic levels, from site specific differences (for example, fixation of a single allele in the population at one site), to differences between population systems (absence of an allele in the Eastern system which is present in all Western system sites). The model therefore explained the past and long term processes which have led to the formation and maintenance of the highly structured Eastern rabbit population system. Most animals exhibit sex biased dispersal which may influence long term connectivity levels among local populations, and thus the dynamics of regional systems. When appropriate sex specific dispersal characteristics were used, the habitat heterogeneity model predicted substantially different interaction patterns between female-only and combined male and female dispersal scenarios. In the latter case, model output was validated using data from a bi-parentally inherited genetic marker. Again, the model explained over 80% of the variation in the genetic data. The fact that such a large proportion of variability is explained in two genetic data sets provides very good evidence that habitat heterogeneity influences long term connectivity levels among local rabbit populations in the Mitchell region for both males and females. The habitat heterogeneity model thus provides a powerful approach for understanding the large scale processes that shape regional population systems in general. Therefore the model has the potential to be useful as a tool to aid in the management of those systems, whether it be for pest management or conservation purposes.

A real time Taqman RT-PCR for the detection of rabbit hemorrhagic disease virus 2 (RHDV2)

As key prey, the wild rabbit downsize constitutes a major drawback on the endangered Iberian lynx (Lynx pardinus) re-introduction in the Iberia. Several captive breeding units mostly located in Alentejo, endeavour the wild rabbit repopulation of depleted areas assigned for the lynx re-introduction. Here we report an RHDV2 outbreak that occurred in early 2016 in a wild rabbit captive breeding unit located in Barrancos municipality. The estimated mortality rate between March and April 2016 was approximately 8.67%. Anatomopathologic examination was carried out for 13 victimized rabbits. Molecular characterization was based on the complete vp60 capsid gene. The 13 rabbit carcasses investigated showed typical macroscopic RHD lesions testing positive to RHDV2- RNA. Comparison of the vp60 nucleotide sequences obtained from two specimens with others publically available disclosed similarities below 98.22% with RHDV2 strains originated in the Iberia and Azores and revealed that the two identical strains from Barrancos-2016 contain six unique single synonymous nucleotide polymorphisms. In the phylogenetic analysis performed, the Barrancos-2016 strains clustered apart from other known strains, meaning they may represent new evolutionary RHDV2 lineages. No clear epidemiological link could be traced for this outbreak where the mortalities were lower compared with previous years. Yet, network analysis suggested a possible connection between the missing intermediates from which the strains from Barrancos 2013, 2014 and 2016 have derived. It is therefore possible that RHDV2 has circulated endemically in the region since 2012, with periodic epizootic occurrences. Still, six years after its emergence in wild rabbits, RHDV2 continues to pose difficulties to the establishment of natural wild rabbit populations that are crucial for the self-sustainability of the local ecosystems.

A method for mapping the distribution and density of rabbits and other vertebrate pests in Australia

The European wild rabbit has been considered Australia’s worst vertebrate pest and yet little effort appears to have gone into producing maps of rabbit distribution and density. Mapping the distribution and density of pests is an important step in effective management. A map is essential for estimating the extent of damage caused and for efficiently planning and monitoring the success of pest control operations. This paper describes the use of soil type and point data to prepare a map showing the distribution and density of rabbits in Australia. The potential for the method to be used for mapping other vertebrate pests is explored. The approach used to prepare the map is based on that used for rabbits in Queensland (Berman et al. 1998). An index of rabbit density was determined using the number of Spanish rabbit fleas released per square kilometre for each Soil Map Unit (Atlas of Australian Soils). Spanish rabbit fleas were released into active rabbit warrens at 1606 sites in the early 1990s as an additional vector for myxoma virus and the locations of the releases were recorded using a Global Positioning System (GPS). Releases were predominantly in arid areas but some fleas were released in south east Queensland and the New England Tablelands of New South Wales. The map produced appears to reflect well the distribution and density of rabbits, at least in the areas where Spanish fleas were released. Rabbit pellet counts conducted in 2007 at 54 sites across an area of south east South Australia, south eastern Queensland, and parts of New South Wales (New England Tablelands and south west) in soil Map Units where Spanish fleas were released, provided a preliminary means to ground truth the map. There was a good relationship between mean pellet count score and the index of abundance for soil Map Units. Rabbit pellet counts may allow extension of the map into other parts of Australia where there were no Spanish rabbit fleas released and where there may be no other consistent information on rabbit location and density. The recent Equine Influenza outbreak provided a further test of the value of this mapping method. The distribution and density of domestic horses were mapped to provide estimates of the number of horses in various regions. These estimates were close to the actual numbers of horses subsequently determined from vaccination records and registrations. The soil Map Units are not simply soil types they contain information on landuse and vegetation and the soil classification is relatively localised. These properties make this mapping method useful, not only for rabbits, but also for other species that are not so dependent on soil type for survival.

Preferencia del conejo silvestre Oryctolagus cuniculus por diferentes tipos de sustratos

Asignatura. Etología de los recursos pesqueros (Licenciatura Ciencias del mar)

Ingeniería de la conservación: Investigación aplicada para la conservación de la fauna mediterránea amenazada

El ser humano nunca ha tenido una capacidad semejante a la actual para alterar el medio en el que vive. A través fundamentalmente de cambios en el uso del hábitat, alteración de los ciclos biogeoquímicos, cambio climático y presencia de especies invasoras, la tasa de extinción de especies se ha acelerado enormemente. Al crecer más rápido la población humana dentro de las áreas en las que se concentra una mayor biodiversidad (biodiversity hotspots), las tasas de transformación son mayores, por lo que dicha biodiversidad se encuentra más amenazada. Entre los principales 25 biodiversity hotspots globales se encuentra la cuenca mediterránea. Las afecciones a la biodiversidad hacen por tanto que en dicha cuenca mediterránea sea necesario poner en práctica iniciativas innovadoras para su conservación. Pero además de una cuidadosa selección, es necesario el seguimiento y la evaluación de la eficacia y la eficienca de dichas medidas, para comprobar su idoneidad y mejorarlas cuando sea posible. Esta es la tarea de la ingeniería de la conservación. Esta Memoria busca analizar y proponer mejoras en cuatro acciones experimentales de conservación de la biodiversidad que tienen un marco común, la conservación de especies subrogadas del monte mediterráneo. En concreto de cinco especies amenazadas: conejo de monte (Oryctolagus cuniculus, especie llave), águila imperial ibérica (Aquila adalberti), águila perdicera (Aquila fasciata), cigüeña negra (Ciconia nigra) y lince ibérico (Lynx pardinus), todas ellas especies paraguas y bandera. Casi todas estas especies se enfrentan a amenazas comunes, como son la pérdida de la calidad del hábitat, la mortalidad de origen humano, la escasez de alimento o la falta de eficiencia en las técnicas de seguimiento. Las acciones destinadas a revertir o minimizar estas amenazas han tenido un marco geográfico y ecológico común, 49 fincas privadas de monte mediterráneo situadas en 9 espacios de la Red Natura 2000 del centro-occidental peninsular, en las que se ha podido trabajar merced a acuerdos de custodia del territorio. Dichas actividades han sido financiadas casi por completo con cargo a fondos LIFE Naturaleza de la Comisión Europea. El presente documento contribuye a aportar el conocimiento científico imprescindible para mejorar la eficacia y eficiencia de las medidas técnicas de gestión y conservación: la ingeniería de la conservación. Como resultados más relevantes se pueden citar que la selección del conejo de las zonas sembradas aumenta con la presencia tanto de isletas de matorral como de refugios y resulta inferior en las zonas de matorral más cerrado; que la electrocución es una de las principales causas de mortalidad de las rapaces amenazadas, así como que se ha detectado que la mortalidad se concentra en determinadas zonas que deben ser prioritarias para su corrección y que éstas deben pasar por las modificaciones estructurales; que la selección de zonas de alimentación por parte de la cigüeña negra se ve afectada por numerosos factores y que una adecuada ubicación y gestión de las charcas para ganadería y caza mayor puede favorecer no sólo la alimentación de esta especie, sino suponer un lugar de cría para diversas especies de anfibios y peces autóctonos, en franco declive; finalmente, hay que destacar que la mejora en el seguimiento de las poblaciones de lince ibérico pasa por la incorporación de métodos de seguimiento novedosos, como las cámaras infrarrojas en letrinas, dada la capacidad de aprendizaje de la especie. Abstract Nowadays the human being has reached its top ability to exchange its environment. Mainly through habitat change, altering biogeochemical cycles, climatic change and alien species, global extinction rate has increased dramatically. Population increases faster in biodiversity hotspots, thus change rates are higher. Mediterranean basin is included within the 25 most valuable biodiversity hotspots. Affections on biodiversity turns necessary to develop innovative actions to maintain the most threatened species in this hotspot. But further than a careful selection and implementarion, survey and evaluation of those practices is necessary to check its suitability. Those are actions for the conservation engineering. This Memory looks for to analyze and to put forward improvements in four experimental biodiversity conservation actions with a common framework: conservation of surrogate Mediterranean species. Those actions were specifically directed to five threatened species: the European wild rabbit (Oryctolagus cuniculus, key species), Spanish imperial eagle (Aquila adalberti) and Bonelli’s eagle (Aquila fasciata, both of the umbrella and flagship species), black stork (Ciconia nigra, umbrella and flagship species) and the Iberian lynx (Lynx pardinus, umbrella and flagship species). Those species have common threats: habitat loss, maninduced mortality, food shortening and lack of surveying efficiency. Actions aimed to minimize those threats also had a common geographical frame: 49 privately owned states within 9 Natura 2000 areas in Southwestern Spain. Actions were developed through lands tewardship agreements, and were even wholly undertaken through European Commision LIFE’s fund. This document contributes to develop cientific knowledge necessary for increasing efficacy and efficiency for implementing technical measures devoted to conservation, the conservation engineering. One of the most outstanding result might be the necessity of including bushy islets or rabbit shelter within sowings when enhacing rabbit populations is the goal. Electrocution is one of the commonest man-induced death causes for raptors, especially for the threatened ones. Areas of concentration of this mortality have been detected, and when mitigation measures have been evaluated, only structural modifications had results. The black stork’s selection of feeding areas is influenced by several factors. An adequate location and management of livestock’s ponds can favour not only this species, but also creating new breeding places for amphibians and indigenous fish species, most of them in a sharp decline. Finally, improving the survey of the Iberian lynx might be achieved through implementing newly developed surveying systems, such as the infrared cameras in latrines.

Influência da gestão do habitat na distribuição e abundância do coelho-bravo no Sítio Natura 2000 Monchique

O coelho-bravo, devido à sua importância ecológica e económica, tem sido alvo de diversos planos de gestão e vários esforços têm sido empreendidos no sentido de contrariar o decréscimo das suas populações. Este estudo foi realizado em três zonas de caça do Sítio Monchique e o principal objectivo é determinar se as medidas de gestão implementadas influenciam a distribuição e abundância da espécie na área de estudo. A abundância relativa foi interpolada com o método "Inverso do Peso da Distância" {IDW), e as relações entre presença de coelho e os descritores ambientais foram analisadas através de Modelos Lineares Generalizados (GLM). Os resultados da modelação estatística mostraram que as medidas de melhoria de habitat parecem ter sido determinantes para um aumento da área de distribuição do coelho-bravo nos locais intervencionados. São propostas novas medidas de gestão, cujo objectivo será promover a continuação do aumento da ocorrência e abundância da espécie neste local. /ABSTRACT: The wild rabbit, due to its ecological and economical role, has been the target of several management plans and considerable efforts have been made to enhance its populations. This study was held in three game estates located inside Monchique Natura 2000. Site and aims to determine if the habitat management actions implemented in the study area influence rabbit distribution and abundance. The relative abundance was interpolated to all study area with lnverse Distance Weight method {IDW} and the relationships between rabbit presence and the environmental descriptors were evaluated with Generalized Linear Models (GLM). The results of the statistical modelling showed that the management actions seem to have contributed significantly to an enhancement on the rabbit occurrence in the studied game estates. Several new management actions are proposed with the aim to continue to increase rabbit occurrence and abundance in this site.

Phylogeographic Triangulation: Using Predator-Prey-Parasite Interactions to Infer Population History from Partial Genetic Information

Phylogeographic studies, which infer population history and dispersal movements from intra-specific spatial genetic variation, require expensive and time-consuming analyses that are not always feasible, especially in the case of rare or endangered species. On the other hand, comparative phylogeography of species involved in close biotic interactions may show congruent patterns depending on the specificity of the relationship. Consequently, the phylogeography of a parasite that needs two hosts to complete its life cycle should reflect population history traits of both hosts. Population movements evidenced by the parasite’s phylogeography that are not reflected in the phylogeography of one of these hosts may thus be attributed to the other host. Using the wild rabbit (Oryctolagus cuniculus) and a parasitic tapeworm (Taenia pisiformis) as an example, we propose comparing the phylogeography of easily available organisms such as game species and their specific heteroxenous parasites to infer population movements of definitive host/predator species, independently of performing genetic analyses on the latter. This may be an interesting approach for indirectly studying the history of species whose phylogeography is difficult to analyse directly.

Increased virulence of rabbit haemorrhagic disease virus associated with genetic resistance in wild Australian rabbits (Oryctolagus cuniculus)

The release of myxoma virus (MYXV) and Rabbit Haemorrhagic Disease Virus (RHDV) in Australia with the aim of controlling overabundant rabbits has provided a unique opportunity to study the initial spread and establishment of emerging pathogens, as well as their co-evolution with their mammalian hosts. In contrast to MYXV, which attenuated shortly after its introduction, rapid attenuation of RHDV has not been observed. By studying the change in virulence of recent field isolates at a single field site we show, for the first time, that RHDV virulence has increased through time, likely because of selection to overcome developing genetic resistance in Australian wild rabbits. High virulence also appears to be favoured as rabbit carcasses, rather than diseased animals, are the likely source of mechanical insect transmission. These findings not only help elucidate the co-evolutionary interaction between rabbits and RHDV, but reveal some of the key factors shaping virulence evolution.

Temporal dynamics of rabbit haemorrhagic disease virus infection in a low-density population of wild rabbits (Oryctolagus cuniculus) in New Zealand

A longitudinal capture-mark-recapture study was conducted to determine the temporal dynamics of rabbit haemorrhagic disease (RHD) in a European rabbit (Oryctolagus cuniculus) population of low to moderate density on sand-hill country in the lower North Island of New Zealand. A combination of sampling ( trapping and radio-tracking) and diagnostic (cELISA, PCR and isotype ELISA) methods was employed to obtain data weekly from May 1998 until June 2001. Although rabbit haemorrhagic disease virus ( RHDV) infection was detected in the study population in all 3 years, disease epidemics were evident only in the late summer or autumn months in 1999 and 2001. Overall, 20% of 385 samples obtained from adult animals older than 11 weeks were seropositive. An RHD outbreak in 1999 contributed to an estimated population decline of 26%. A second RHD epidemic in February 2001 was associated with a population decline of 52% over the subsequent month. Following the outbreaks, the seroprevalence in adult survivors was between 40% and 50%. During 2000, no deaths from RHDV were confirmed and mortalities were predominantly attributed to predation. Influx of seronegative immigrants was greatest in the 1999 and 2001 breeding seasons, and preceded the RHD epidemics in those years. Our data suggest that RHD epidemics require the population immunity level to fall below a threshold where propagation of infection can be maintained through the population.