The value of having no wild rabbits in south-east Queensland

Rabbits released in Australia in 1859 spread to most areas of suitable habitat by 1910 causing great damage to the environment and primary industries. Measurement of damage is essential to justify spending money and utilising resources to remove rabbits. Damage to pasture and biodiversity may be irreversible and therefore difficult to measure without comparison with an area that has never suffered such damage. A rabbit proof fence completed in 1906 protected a large part of south east Queensland from rabbits. The Darling Downs Moreton Rabbit Board (DDMRB) continues to maintain the fence and keep the area relatively free of rabbits. This area is unique because it is highly suitable for rabbits and yet it has never ‘experienced’ the damage caused by plagues of uncontrolled rabbits. A study site was established where the DDMRB fence separates an area heavily used by rabbits (‘dirty side’) from an area that has never been infested by rabbits (‘clean side’). The number and location of all rabbit warrens and log piles were recorded. The absence of warrens from the ‘clean side’ shows clearly that the rabbit proof fence has prevented rabbits from establishing warren systems. The ‘dirty side’ is characterised by a high number of warrens, a high density of rabbits, fewer pasture species and low macropod activity. Future work will determine whether the rabbit populations are viable in the absence of rabbit warrens. We plan to radio collar rabbits on both sides of the fence to measure their survival rate. In selected warrens and log piles of varying degrees of complexity and size, rabbits will be trapped and information on reproduction and age structure will be collected. This will allow better targeting of the source of rabbits during control operations. Once the initial comparative analysis of the site has been completed, all rabbit warrens will be destroyed on the dirty side of the fence. After rabbits are removed from this area, monitoring will continue to determine if pasture and biodiversity on opposite sides of the fence begin to mirror each other.

Co-evolution of wild rabbits (Oryctolagus cuniculus) and RHDV: resistance and virulence

Rabbit Haemorrhagic Disease Virus (RHDV) was introduced to Australia in 1995 for the control of wild rabbits. Initial outbreaks greatly reduced rabbit numbers and the virus has continued to control rabbits to varying degrees in different parts of Australia. However, recent field evidence suggests that the virus may be becoming less effective in those areas that have previously experienced repeated epizootics causing high mortality. There are also reports of rabbits returning to pre-1995 density levels, Virus and host can be expected to co-evolve. The host will develop resistance to the virus with the virus subsequently changing to overcome that resistance. It has been 12 years since the release of RHDV and it is an opportune time to examine where the dynamic currently stands between RHDV and rabbits. Laboratory challenge tests have indicated that resistance to RHDV has developed to different degrees in populations throughout Australia. In one population a low dose (1:25 dilution) of Czech strain RHDV failed to infect a single susceptible rabbit, yet infected a low to high (up to 73%) percentage across other populations tested. Different selection pressures are present in these populations and will be driving the level of resistance being seen. The mechanisms and genetics behind the development of resistance are also important as the on-going use of RHDV as a control tool in the management of rabbits relies on our understanding of factors influencing the efficacy of the virus. Understanding how resistance has developed may provide clues on how best to use the virus to circumvent these mechanisms. Similarly, it will help in managing populations that have yet to develop high levels of resistance.

How can warren destruction by ripping control European wild rabbits (Oryctolagus cuniculus) on large properties in the Australian arid zone?

Context: For over 100 years, control efforts have been unable to stop rabbits causing damage to cattle production and native plants and animals on large properties in arid parts of Australia. Warren destruction by ripping has shown promise, but doubts about long-term success and the perceived expense of treating vast areas have led to this technique not being commonly used. Aims: This study measured the long-term reduction in rabbit activity and calculated the potential cost saving associated with treating just the areas where rabbits are believed to survive drought. Wealso considered whether ripping should be used in a full-scale rabbit control program on a property where rabbits have been exceptionally resilient to the influence of biological and other control measures. Methods: Rabbits were counted along spotlight transects before warrens were ripped and during the two years after ripping, in treated and untreated plots. Rabbit activity was recorded to determine the immediate and long-term impact of ripping, up to seven years after treatment. The costs of ripping warrens within different distances from drought refuge areas were calculated. Key results: Destroying rabbit warrens by ripping caused an immediate reduction in rabbit activity and there were still 98% fewer rabbits counted by spotlight in ripped plots five months after ripping. Seven years after ripping no active warrens were found in ripped plots, whereas 57% of warrens in unripped plots showed signs of rabbit activity. The cost of ripping only the areas where rabbits were likely to seek refuge from drought was calculated to be less than 4%of the cost of ripping all warrens on the property. Conclusions: Destroying rabbit warrens by ripping is a very effective way of reducing rabbit numbers on large properties in arid Australia. Ripping should commence in areas used by rabbits to survive drought. It is possible that no further ripping will be required. Implications: Strategic destruction of warrens in drought refuge areas could provide an alternative to biological control for managing rabbits on large properties in the Australian arid zone.

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.

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.

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.

Rabbit haemorrhagic disease: are Australian rabbits (Oryctolagus cuniculus) evolving resistance to infection with Czech CAPM 351 RHDV?

Rabbit haemorrhagic disease is a major tool for the management of introduced, wild rabbits in Australia. However, new evidence suggests that rabbits may be developing resistance to the disease. Rabbits sourced from wild populations in central and southeastern Australia, and domestic rabbits for comparison, were experimentally challenged with a low 60 ID50 oral dose of commercially available Czech CAPM 351 virus - the original strain released in Australia. Levels of resistance to infection were generally higher than for unselected domestic rabbits and also differed (0-73% infection rates) between wild populations. Resistance was lower in populations from cooler, wetter regions and also low in arid regions with the highest resistance seen within zones of moderate rainfall. These findings suggest the external influences of non-pathogenic calicivirus in cooler, wetter areas and poor recruitment in arid populations may influence the development rate of resistance in Australia.

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.

Survival of rabbit haemorrhagic disease virus (RHDV) in the environment

A study was conducted to investigate the persistence of rabbit haemorrhagic disease virus (RHDV) in the environment. Virus was impregnated onto two carrier materials (cotton tape and bovine liver) and exposed to environmental conditions on pasture during autumn in New Zealand. Samples were collected after 1, 10, 44 and 91 days and the viability of the virus was determined by oral inoculation of susceptible 11- to 14-week-old New Zealand White rabbits. Evidence of RHDV infection was based on clinical and pathological signs and/or seroconversion to RHDV. Virus impregnated on cotton tape was viable at 10 days of exposure but not at 44 days, while in bovine liver it was still viable at 91 days. The results of this study suggest that RHDV in animal tissues such as rabbit carcasses can survive for at least 3 months in the field, while virus exposed directly to environmental conditions, such as dried excreted virus, is viable for a period of less than I month. Survival of RHDV in the tissues of dead animals could, therefore, provide a persistent reservoir of virus, which could initiate new outbreaks of disease after extended delays.

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.

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.

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.