Favourable areas for expansion and reintroduction of Iberian lynx accounting for distribution trends and genetic diversity of the European rabbit

Although on a local scale Iberian lynx distribution is determined by the availability of prey rabbits, recent modelling analyses have uncovered broad-scale disagreements between these two species’ distribution trends. These analyses showed also that the lynx had become restricted to only a fraction of the rabbit’s genetic diversity, and that this could be jeopardising its survival in the face of environmental hazards and uncertainty. In the present paper, a follow-up was carried out through the building of lynx and rabbit distribution models based on the most recent Spanish mammal atlas. Environmental favourability values for lynx and rabbit were positively correlated within the lynx’s current distribution area, but they were negatively correlated within the total Spanish area where lynx used to occur in the 1980’s. Environmental favourability for rabbits was significantly higher where lynx maintains reproductive populations than where it recently disappeared, indicating that rabbit favourability plays an important role and can be a good predictor of lynx persistence. The lynx and rabbit models were extrapolated to predict favourable areas for both species in Spain as well as in Portugal, on the original scale of the distribution data (10x10 km) and on a 100 times finer spatial resolution (1x1 km). The lynx and rabbit models were also combined through fuzzy logic to forecast the potential for lynx occurrence incorporating information on favourable areas for its main prey. Several areas are proposed as favourable for lynx expansion or re-introduction,

Conservation biogeography of ecologically interacting species: The case of the Iberian lynx and the European rabbit

Aim To relate the recent Iberian lynx decline to changes in the distribution of the European rabbit after the haemorrhagic disease outbreak of 1989. As Iberian rabbits evolved in two geographically separated lineages, being the recent lynx range practically restricted to the southwestern lineage, we also test if differential range dynamics exists for these lineages, with the consequent implications for lynx conservation and reintroduction planning.

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.

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.

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.

Challenges in the rabbit haemorrhagic disease 2 (RHDV2) molecular diagnosis of vaccinated rabbits

Molecular methods are fundamental tools for the diagnosis of viral infections. While interpretation of results is straightforward for unvaccinated animals, where positivity represents ongoing or past infections, the presence of vaccine virus in the tissues of recently vaccinated animals may mislead diagnosis. In this study, we investigated the interference of RHDV2 vaccination in the results of a RT-qPCR for RHDV2 detection, and possible associations between mean Cq values of five animal groups differing in age, vaccination status and origin (domestic/wild). Viral sequences from vaccinated rabbits that died of RHDV2 infection (n = 14) were compared with the sequences from the commercial vaccines used in those animals. Group Cq means were compared through Independent t-test and One-way ANOVA. We proved that RHDV2 vaccine-RNA is not detected by the RT-qPCR as early as 15 days post- vaccination, an important fact in assisting results interpretation for diagnosis. Cq values of vaccinated and non-vaccinated infected domestic adults showed a statistically significant difference (p < 0.05), demonstrating that vaccination-induced immunity reduces viral loads and delays disease progression. Contrarily, in vaccinated young rabbits higher viral loads were registered compared to non-vaccinated kittens. No significant variation (p = 0.3824) was observed between viral loads of non- vaccinated domestic and wild RHDV2-victimised rabbits. Although the reduced number of vaccinated young animals analysed hampered a robust statistical analysis, this occurrence suggests that passively acquired maternal antibodies may inhibit the active immune response to vaccination, delaying protection and favouring disease progression. Our finding emphasises the importance of adapting kitten RHDV2 vaccination schedules to circumvent this interference phenomenon.