Next Article in Journal
Comparison of High-Performance Liquid Chromatography with Sucrose Density Gradient Ultracentrifugation for the Quantification of Foot-and-Mouth Disease Vaccine Antigens
Previous Article in Journal
Designing of a Novel Multi-Antigenic Epitope-Based Vaccine against E. hormaechei: An Intergraded Reverse Vaccinology and Immunoinformatics Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Immunological Cross-Protection between Different Rabbit Hemorrhagic Disease Viruses—Implications for Rabbit Biocontrol and Vaccine Development

1
Virology Laboratory, Elizabeth Macarthur Agricultural Institute, NSW Department of Primary Industries, Menangle, NSW 2568, Australia
2
Health & Biosecurity, Commonwealth Scientific and Industrial Research Organisation, Acton, ACT 2601, Australia
3
Centre for Invasive Species Solutions, Bruce, ACT 2617, Australia
*
Author to whom correspondence should be addressed.
Vaccines 2022, 10(5), 666; https://doi.org/10.3390/vaccines10050666
Submission received: 29 March 2022 / Revised: 19 April 2022 / Accepted: 20 April 2022 / Published: 22 April 2022

Abstract

:
The use of rabbit hemorrhagic disease virus (RHDV) as a biocontrol agent to control feral rabbit populations in Australia, in combination with circulating endemic strains, provides a unique environment to observe the interactions between different lagoviruses competing for the same host. Following the arrival of RHDV2 (GI.2) in Australia, it became necessary to investigate the potential for immunological cross-protection between different variants, and the implications of this for biocontrol programs and vaccine development. Laboratory rabbits of various immune status—(1) rabbits with no detectable immunity against RHDV; (2) rabbits with experimentally acquired immunity after laboratory challenge; (3) rabbits immunised with a GI.2-specific or a multivalent RHDV inactivated virus prototype vaccine; or (4) rabbits with naturally acquired immunity—were challenged with one of three different RHDV variants (GI.1c, GI.1a or GI.2). The degree of cross-protection observed in immune rabbits was associated with the variant used for challenge, infectious dose of the virus and age, or time since acquisition of the immunity, at challenge. The immune status of feral rabbit populations should be determined prior to intentional RHDV release because of the high survival proportions in rabbits with pre-existing immunity. In addition, to protect domestic rabbits in Australia, a multivalent RHDV vaccine should be considered because of the limited cross-protection observed in rabbits given monovalent vaccines.

1. Introduction

Rabbit hemorrhagic disease virus (RHDV), a member of the genus Lagovirus in the family Caliciviridae, is used effectively as a biocontrol agent in Australia to reduce the significant environmental and economic impact of feral European rabbit (Oryctolagus cuniculus) populations [1]. Specifically, viruses within the GI.1 genotype are deliberately released as part of these control efforts, with a GI.1c (classical RHDV) virus in use since 1995 and a GI.1a (RHDV1a) virus in use since 2017 [2,3,4]. In addition to these intentionally released lagoviruses, a non-pathogenic rabbit calicivirus, GI.4 (RCV-A1), is endemic, while another virus, GI.2 (RHDV2), is currently the dominant pathogenic variant across Australia [5,6]. Since the first detection of GI.2 in Australia, the original GI.2 variant (GI.1bP-GI.2) has been replaced by several other naturally occurring GI.2 recombinants, all of which comprise nearly identical VP60 capsid proteins [5,7,8]. The deliberate release of RHDV, in combination with these naturally circulating lagoviruses in Australia, provides a unique environment to observe the evolving interactions between different lagoviruses. These viruses compete for the same host, and this has implications for both the biological control of feral rabbits as well as appropriate vaccine protection for domestic rabbits in Australia.
Rabbit biocontrol programs are very successful in Australia, with the resulting feral population decline delivering an estimated $1 billion AUD in economic and environmental benefits annually [1,3,9]. Nevertheless, both intentionally released and naturally circulating lagoviruses also have significant impacts on domestic rabbits. In a susceptible rabbit, infection with RHDV causes a fatal, fulminant necrotising hepatitis and disseminated intravascular coagulation [10]. Currently, there is one lagovirus vaccine registered in Australia for domestic rabbits (Cylap® RCD, Zoetis, NJ, USA). This vaccine offers complete protection against both the GI.1c and GI.1a variants used in biocontrol programs [11], but GI.1 vaccines are largely ineffective against GI.2 viruses [12]. In addition, this inactivated vaccine is imported and there can be supply shortages, leaving domestic rabbits in Australia vulnerable. Several GI.2 replication-competent, viral-vectored vaccines are manufactured overseas, but these are not suitable for the Australian context. The theoretical transmission of these vaccine viruses into the feral rabbit population could significantly reduce the economic savings and environmental benefits of rabbit control programs [13]. To date, European inactivated vaccines have not been imported and registered in Australia, possibly due to a small market [14]. Consequently, there is a need for a suitable Australian vaccine to protect domestic rabbits against unintentional infection from all circulating pathogenic lagoviruses (GI.1 and GI.2).
Cross-protection between different RHDV variants was investigated previously in the context of both natural immunity acquired after challenge and vaccinal immunity [11,15]. Australian rabbits previously exposed to GI.4 were partially and transiently protected against GI.1c challenge [16], and lead to the investigation of a GI.1a variant (Korean origin) as an additional biocontrol agent to overcome this cross-protection [4]. There is limited immunological cross-protection between genotypes GI.1a and GI.2 [15,17]. This lack of cross-protection is promising for the selection of a GI.2 variant, GI.1bP-GI.2, as the next biocontrol agent, where this virus may overcome existing immunity in a feral rabbit population. However, limited cross-protection is problematic for vaccine development as an effective vaccine must protect a domestic Australian rabbit against GI.1a, GI.1c and GI.2.
The objectives of this research were to establish whether prior immunity to different RHDV genotypes could offer cross-protection against GI.2. Three types of immunity were investigated: firstly, rabbits with experimentally acquired immunity after primary infection with the current GI.1a biocontrol variant or the non-pathogenic GI.4 virus (Figure 1A); secondly, rabbits vaccinated with either a prototype GI.2-specific vaccine or a multivalent RHDV vaccine (Figure 1B); and, finally, rabbits that serendipitously had naturally acquired immunity to GI.1 or GI.2 (Figure 1C). Rabbits from these three groups were challenged with RHDV variants from either the GI.1a, GI.1c or GI.2 genotypes. The interactions between the different Australian lagoviruses explored here have implications for selection and delivery of virus strains for biocontrol purposes and vaccine composition for optimal protection.

2. Materials and Methods

2.1. Animals and Experimental Design

Young domestic rabbits with various immunity and infection statuses were obtained from three geographically distinct regions in Australia. Both the cross-protection experiment on experimentally acquired immunity (Figure 1A, Table 1) and vaccinal immunity (Figure 1B) were planned. For these experiments, sample sizes were estimated based on demonstrating a difference in survival between groups. The investigation of cross-protection with rabbits with naturally acquired immunity (Figure 1C, Table 1) was serendipitous and sample sizes were determined by the number of available rabbits. All rabbits were randomly allocated into their respective treatment groups. All procedures involving animals were carried out in accordance with the “Australian Code of Practice for the Care and Use of Animals for Scientific Purposes” and approved by the NSW Department of Primary Industries, Elizabeth Macarthur Agricultural Institute (EMAI, Menangle, Australial), Animal Ethics Committee (M20-02). Challenge studies were performed at EMAI.

2.2. Cross-Protection in Rabbits with Experimentally Acquired Immunity

New Zealand White rabbits of both sexes were obtained from a colony known to be free of benign and pathogenic lagoviruses. Prior to both experimental challenges, all rabbits were confirmed to be seronegative to Australian lagoviruses by ELISA [18].
To obtain rabbits with GI.1a immunity, five-week-old rabbits were orally inoculated with 100 ‘50% Rabbit Infectious Doses’ (RID50) of a commercially available GI.1a infectious virus preparation (GenBank acc. Number MF598301.1) [11]. The RID50 of the GI.1a challenge was titrated as part of the quality control tests for commercial preparations of RHDV. To obtain rabbits with GI.4c immunity, five-week-old rabbits were given a 1 mL oral dose (with an estimated 2300 genome copies) of a clarified homogenate of duodenal tissue from a GI.4c-infected rabbit (GenBank acc. Number KX357705.1). Titration to determine RID50 of GI.4c was not performed. Therefore, ‘genome copies’ is used as an in vitro approximate to quantify the GI.4c challenge.
After initial challenge, all rabbits seroconverted. To simulate infection against peak immunity, rabbits were then challenged with GI.2 within five weeks of initial infection [16]. Thus, 39 days after infection with GI.1a or GI.4c viruses (i.e., at 12 weeks of age), rabbits were challenged per os (PO) with 50 RID50 of GI.1bP-GI.2 (GenBank acc. Number MW467791) (Figure 1A).

2.3. Cross-Protection in Rabbits Vaccinated with a Prototype Vaccine

Rabbits from a commercial breeding facility that were seronegative to GI.1 and GI.2 were vaccinated with a single dose of either a prototype, inactivated GI.2-specific vaccine or a prototype inactivated multivalent RHDV vaccine (with GI.1a, GI.1c and GI.2 components). Rabbits were vaccinated by subcutaneous injection at 10–12 weeks of age and then challenged PO with a heterologous lagovirus 28 days after vaccination (Figure 1B).

2.4. Cross-Protection in Rabbits with Naturally Acquired Immunity

Rabbits with naturally acquired immunity to lagoviruses were obtained from three commercial rabbit farms in Australia. These rabbits were initially obtained for other purposes. Following the detection of anti-GI.1 or anti-GI.2 antibodies based on serology, or GI.2 disease, these rabbits were then serendipitously selected for inclusion in this study. Although rabbits were not tested for anti-GI.4 antibodies, given the widespread distribution of this virus in Australian rabbit farms, these rabbits could possibly also have had anti-GI.4 antibodies [16]. New Zealand White rabbits from one farm (farm 1) were unvaccinated but were seropositive to GI.1. This was likely to be due to either previous non-lethal exposure to a GI.1c or GI.1a virus or persisting maternal antibodies. Most of these rabbits with GI.1 immunity were challenged PO at 12 weeks of age with GI.2 virus, while some were retained and challenged PO at 33 weeks of age to assess if this immunity waned over time (Figure 1C). New Zealand White rabbits from another farm (farm 2) were unvaccinated but were seropositive to GI.2. The origin of this antibody was likely to be either previous non-lethal exposure to a GI.2 virus or persisting maternal antibodies. These rabbits with GI.2 immunity were challenged PO at 12 weeks of age (Figure 1C). From a third farm (farm 3), New Zealand White cross Flemish Giant rabbits were unvaccinated and within 48 h of arrival at EMAI, sudden deaths occurred in this group. Molecular testing revealed infection with a GI.4cP-GI.2 (4c-recombinant) virus, a variant that was (at the time of this study) endemic to the location of this third farm and not previously detected in the region where EMAI is located (GenBank acc. Number MW460156) [19]. One month after the last signs of disease, surviving rabbits from this GI.4cP-GI.2 infection were challenged PO with a GI.1a or a GI.2 variant (GI.1bP-GI.2) (Figure 1C).

2.5. Virus Challenge

To mimic a natural route of infection, all rabbits were challenged by the oral administration of a 1 mL volume of virus solution. To assess cross-protection between circulating Australian lagoviruses, rabbits were challenged with various doses (Table 1) of (1) the GI.1c biocontrol preparation (GenBank acc. Number KT344772.1) (EMAI, Menangle, Australia); (2) the GI.1a biocontrol preparation (EMAI, Menangle, Australia); or (3) a GI.2 (GI.1bP-GI.2) preparation [20].
Initially, a dose of 50 RID50 was selected as the challenge dose for rabbits with experimentally or naturally acquired immunity. This was based on an estimate of the amount of RHDV a rabbit may be exposed to during a control program [21]. However, to assess whether responses were dose dependent, in other experiments, this was increased to 150 RID50 (low dose in Figure 1C) or 1500 RID50 (high dose in Figure 1C). Additionally, an extreme dose of 150,000 RID50 was used to determine if this could overcome natural immunity to a different GI.2 variant (that is, a GI.1bP-GI.2 challenge of rabbits with natural GI.4cP-GI.2 immunity). All vaccinated rabbits were given a high 1500 RID50 PO dose (Figure 1B) as a routine and substantial demonstration of protection.

2.6. Vaccine

Prototype, inactivated, whole-virus vaccines were produced by inoculating naive rabbits with a virus of either the GI.1a, GI.1c or GI.2 genotype. In short, a semi-purified preparation of each virus was then prepared from liver homogenates and inactivated using binary ethylenimine [22]. The antigen preparation, measured in hemagglutination (HA) units per mL [23], was diluted to the desired concentration in phosphate buffered saline and adjuvanted with an equal volume of Montanide ISA-201 (Seppic, Courbevoie, France). Rabbits were injected subcutaneously with either 1 mL of the G1.2-specific vaccine (representing a 100 HA dose), or 0.5 mL of the multivalent RHDV vaccine (consisting of a 32 HA dose of each genotype—GI.1a, GI.1c and GI.2).

2.7. Monitoring, Sample Collection

All rabbits were group housed pre-challenge and then held individually in separated cages in climate-controlled, biosafety level 2 rooms during their challenge period to avoid cross-infection between animals. Additionally, experiments were staggered to allow adequate decontamination and fumigation between challenges to avoid accidental cross-infection.
To establish the immune status to RHDV, serum samples were collected at several time points from all rabbits: (1) on arrival; (2) prior to vaccination or intentional exposure to RHDV to obtain natural immunity; (3) prior to challenge; and (4) after euthanasia.
Rabbits were monitored twice daily during the seven-day or fourteen-day challenge period for the development of clinical signs of RHDV (such as sudden death, inappetence, reduced water intake, and reduced faecal output). While it was anticipated that moribund rabbits would be euthanised with an overdose of intravenous barbiturate, euthanasia was not required during the challenge period (i.e., all rabbits that succumbed to RHDV infection in this study died peracutely). All surviving rabbits were humanely killed after this challenge period with an overdose of intravenous barbiturate; however, the rabbits with immunity to GI.2 that survived heterologous challenge were repurposed into other research projects.

2.8. Data Analyses

Survival analyses was conducted with R Version 4.0.4 [24]. Packages used for these experiments include survminer [25] and survival [26]. A statistical difference between groups was considered to be significant if the log-rank test gave a p-value less than 0.05.

2.9. Serological Testing

Three specific ELISAs were used to detect antibodies against GI.1, GI.4 and GI.2 to determine the immune status of rabbits prior to challenge. Details of GI.1 and GI.4 antibody ELISAs are as published [27,28]. A blocking ELISA for GI.2 antibodies, similar to that reported in a previous publication [18], was developed with polyclonal RHDV rabbit anti-sera as the capture antibody, semi-purified GI.2 antigen, and a conjugated GI.2-specific monoclonal antibody (4H12, supplied from IZSLER, Brescia, Italy). The assay was validated on serum samples collected from known seronegative and seropositive rabbits. The cut-off point for a positive antibody result was set at a percentage inhibition of greater than 60%.
To help clarify the immune status of individual rabbits, where cross-reactivity between the different antigens was suspected, samples were also tested in three separate hemagglutination inhibition (HI) tests against GI.1a, GI.1c and GI.2 antigens [23].

2.10. Molecular Testing

The presence of viral RNA in swabs of the cut surface of livers was determined using reverse-transcription quantitative PCR (RT-qPCR) assays specific to GI.1, GI.2 and GI.4c [6,29,30]. Details of the nucleic acid extraction and PCR amplification methods and equipment used follow those described previously [31].

3. Results

3.1. Survival of Seronegative Rabbits Following Various RHDV Challenges

Following RHDV challenge, all seronegative rabbits challenged with GI.1a or GI.1c died, and 7/8 challenged with GI.2 died within 14 days of infection (Table 1, Figure 2). The surviving GI.2-infected rabbit (challenged with a low (50 RID50) dose of GI.2) remained seronegative after challenge and GI.2 was not detected by RT-qPCR testing on collected liver swab, indicating that it did not become infected. No seronegative rabbits challenged with a 1500 RID50 dose of GI.2 (0/12), GI.1a (0/12) or GI.1c (0/12) survived, with all succumbing to infection within 4 days post challenge (dpc), 6 dpc and 8 dpc, respectively. Median survival for a 50 RID50 dose of GI.2 was 4 dpc, compared to 3 dpc for a 1500 RID50 dose of GI.2. Median survival times for a 1500 RID50 dose of GI.1a and GI.1c were 3.5 and 3 dpc, respectively. These differences were not statistically significant (p = 0.3).

3.2. Survival of Rabbits with Experimentally Acquired Immunity Following GI.2 Challenge

Following a GI.2 challenge of 50 RID50, all (7/7) rabbits with GI.1a acquired immunity survived, as did 7/9 rabbits with GI.4c acquired immunity (Table 1, Figure 3). The survival probability was not statistically different between the two variants (p = 0.2). Notably, there was strong evidence for cross-protection against GI.2 challenge for both variants, with a statistically significant survival probability compared to seronegative controls (p < 0.0001) given the same 50 RID50 GI.2 challenge.

3.3. Survival of Vaccinated Rabbits Following Various RHDV Challenges

The GI.2-specific vaccine afforded complete protection against GI.2 challenge (Table 2, Figure 4C), with all nine rabbits (9/9) surviving a homologous challenge of 1500 RID50. No GI.2-vaccinated rabbits survived infection with either 1500 RID50 of GI.1a (0/6) (Figure 4A) or 1500 RID50 of GI.1c (0/6) (Figure 4B). All rabbits that had received the multivalent vaccine survived infection against GI.1a (6/6), GI.1 (6/6) and GI.2 (6/6) challenge of 1500 RID50 (Table 2, Figure 4D–F).

3.4. Survival of Rabbits with Naturally Acquired Immunity Following Various RHDV Challenges

For rabbits from farm 1, with naturally acquired GI.1 immunity, cross-protection was significantly associated with the variant used for challenge (p = 0.002). That is, at 33 weeks of age and when given 1500 RID50, all GI.1-immune rabbits survived challenge with GI.1a (7/7) but only one survived challenge with GI.2 (1/7) (Table 1, Figure 5A).
There was a trend towards a dose-dependent difference in survival following GI.2 challenge in rabbits with GI.1 immunity. At 12 weeks of age, when given GI.2, nine rabbits (9/10) given 150 RID50 survived challenge but only five rabbits (5/9) given 1500 RID50 survived challenge (Table 1, Figure 5A). However, this was not statistically significant (p = 0.08).
The cross-protection against GI.2 challenge in rabbits with GI.1 immunity was found to reduce with age and/or time since acquisition of GI.1 immunity. When given the same 1500 RID50 dose of GI.2, at 12 weeks of age, five rabbits (5/9) survived challenge; but at 33 weeks of age, only one rabbit (1/7) survived (Table 1, Figure 5A). This difference between age groups was statistically difference (p = 0.03).
Rabbits from farm 2, with naturally acquired GI.2 immunity, showed varying degrees of cross-protection against infection with 1500 RID50 of either GI.1a, GI.1c or GI.2. One rabbit survived GI.1a challenge (1/4), two rabbits survived GI.1c challenge (2/3) and three rabbits survived GI.2 challenge (3/4) (Table 1, Figure 5B). There was no statistical difference in the observed survival between these different challenges (p = 0.3); however, sample sizes were small in this experiment due to the serendipitous nature of acquiring these rabbits. Future robust studies to assess the relationship between GI.2 immunity and heterologous RHDV challenges should be undertaken.
Rabbits from farm 3, that had survived an RHDV (GI.4cP-GI.2) outbreak, had complete cross-protection against a heterologous GI.2 (GI.1bP-GI.2) challenge, with all (13/13) rabbits surviving a 150,000 RID50 dose (Table 1, Figure 5C). Moderate cross-protection was observed against GI.1a, with nine rabbits (9/12) surviving a 150 RID50 dose challenge and nine rabbits (9/10) surviving a 1500 RID50 challenge. There was no statistical difference in the observed survival between these different challenge doses (p = 0.7).

4. Discussion

Immunological cross-protection between GI.1a, GI.1c and GI.2 lagoviruses has implications for both the use of these viruses as biocontrol agents in feral rabbit populations and the composition of vaccines to provide optimal protection for domestic rabbits in Australia. Here, we present a series of experimental studies in rabbits with various immune statuses and their subsequent response to challenge to Australian lagoviruses.
As expected, all seronegative rabbits succumbed to disease following infection with GI.1c, GI.1a, or GI.2 viruses, demonstrating a very high case fatality rate for all viruses investigated. This is consistent with previous virulence studies with these variants [11,20]. In contrast, we observed significant immunological cross-protection in rabbits with either experimentally or naturally acquired immunity. The degree of cross-protection observed in these immune rabbits was associated with the variant given, the virus dose and the time between the challenge and the previous infection. This is reminiscent of previous reports on the partial cross-protection to GI.1 challenge conveyed by GI.4 immunity, which rapidly declined within weeks following the GI.4 infection [16]. Overall, our findings contrast with previous reports in Europe where there appears to be more limited cross-protection between variants [15,17], as well as reports of GI.1-immune rabbits succumbing to GI.2 in Australia [32], although the time between infections with the two viruses was unknown in that study.
The origin of the naturally acquired immunity for the rabbits from farm 1, with detectable GI.1 antibodies, or for the rabbits from farm 2, with detectable GI.2 antibodies, remains unknown. The GI.1 ELISA is unable to distinguish between GI.1a or GI.1c antibodies, and between antibodies derived from natural infection or GI.1 vaccination. The known presence of GI.4 virus circulating at farm 2 also complicates the serological picture in these rabbits. Rabbits from farm 3 arrived at EMAI unvaccinated and infected with GI.4cP-GI.2. Therefore, the origins of their immunity may be assumed to be a consequence of surviving GI.4cP-GI.2 infection.
The source of the naturally acquired immunity may be due to maternal antibodies. A breeding doe may confer maternal antibodies to her kittens for up to one year after RHDV vaccination [33]. The rabbits from farm 1 were acquired from does vaccinated against GI.1 within one month of their birth. For GI.1, the duration of the protective effects from maternal immunity after infection varies according to the antibody titre of the doe, but was not found to persist beyond approximately 14 weeks of age [34]. The rabbits from farm 2 were acquired from does that survived GI.2 within the previous six months before their birth. For GI.2, this maternally derived antibody may persist until 28 days of life [33].
Alternatively, the naturally acquired immunity could be from previous sub-lethal exposure to these viruses. Both suppliers had no recent reports of disease; but given that GI.1 does not generally cause disease in rabbits less than nine weeks of age, natural infection with these viruses cannot be excluded. GI.1a is known to circulate naturally in the geographic region from which these rabbits were obtained [19]. A cohort of these rabbits with GI.1 immunity were kept until 33 weeks of age. At this older age, antibodies were still detectable, confirming that immunity was not due to maternal antibodies. In addition, these rabbits had complete cross-protection against a substantial PO challenge (1500 RID50) with GI.1a at 33 weeks of age, suggesting that this immunity was not waning. Therefore, it is more likely that these rabbits acquired immunity due to natural exposure to GI.1 at farm 1.
Likewise, this same scenario could exist for the rabbits with GI.2 immunity from farm 2. GI.2 can cause disease in younger rabbits from four weeks of age, and the pathogenicity of Australian GI.2 variants is close to 100% in susceptible animals [20]. However, young rabbits with GI.2 maternally derived immunity are resistant to disease with GI.2 [21,35]. Experimental work using passively immunised kittens, to mimic the presence of maternal antibodies, demonstrated a robust antibody response after virus challenge [21].
Despite an extremely high 150,000 RID50 dose, all rabbits that had survived a GI.4cP-GI.2 infection did not succumb to a heterologous GI.1bP-GI.2 challenge. As both viruses share the same GI.2 capsid protein, it is expected that previous infection with GI.4cP-GI.2 should protect against a lethal challenge of GI.1bP-GI.2 [36]. This immune protection has implications for the use of the GI.2 virus (GI.1bP-GI.2) as a future biocontrol agent. GI.2 is currently the dominant genotype across Australia [19]. It is expected that GI.1bp-GI.2 would not be an effective biocontrol agent in a population of feral rabbits with pre-existing GI.2 immunity. Therefore, the immunity of this population should be established prior to intentional virus release when using RHDV as a biocontrol agent in feral rabbit populations, to achieve optimal population control.
The prototype inactivated GI.2-specific vaccine did not offer any cross-protection against GI.1a or GI.1c (Table 1 and Table 2) in contrast to those rabbits with either experimentally or naturally acquired immunity resulting from infection. This is consistent with previous studies where GI.1 vaccinal protection does not offer an acceptable level of cross-protection to GI.2 [11,15]. The lack of cross-protection from inactivated vaccines is likely due to a difference in immune response as observed with other viral infections [37]. Inactivated vaccines primarily stimulate humoral immunity, but virus infection generates a broader response including strong cell mediated and mucosal immunity [17]. In Australia, to protect domestic rabbits, a multivalent RHDV vaccine that offers complete protection against infection with GI.1a, GI.1c and GI.2 should be considered.

5. Conclusions

The results presented here demonstrate that, in laboratory rabbits, there is cross-protective immunity, even between heterologous strains, that protect against lethal disease. The extent of cross-protection conferred depends on the source of immunity (experimentally acquired, vaccinal protection, or naturally acquired), the virus genotype that has induced the immunity (GI.1a, GI.1c, GI.2 or GI.4c), the timing between acquisition of immunity and challenge, and in some cases the infectious dose of the challenge virus. Consequently, for biocontrol programs, the immunity of feral rabbit population should be established prior to intentional RHDV release. For optimal protection of domestic rabbits from unintentional infection with lagoviruses, a multivalent vaccine providing broad immunity to the recognized viruses should be considered.

Author Contributions

Conceptualization, T.W.O., A.J.R., R.N.H., T.S. and P.D.K.; methodology, T.W.O., A.J.R., R.N.H., T.S. and P.D.K.; formal analysis, T.W.O. and R.N.H.; investigation, T.W.O. and A.J.R.; resources, T.S. and P.D.K.; data curation, T.W.O.; writing—original draft preparation, T.W.O.; writing—review and editing, T.W.O., A.J.R., R.N.H., T.S. and P.D.K.; supervision, A.J.R., R.N.H., T.S. and P.D.K.; project administration, R.N.H., T.S. and P.D.K.; funding acquisition, R.N.H., T.S. and P.D.K. All authors have read and agreed to the published version of the manuscript.

Funding

Funding for this work was provided through the Centre for Invasive Species Solutions to investigate RHDV2 as a potential additional biocontrol agent to manage invasive wild rabbits in Australia (Project# P01-B001).

Institutional Review Board Statement

The animal study protocol was approved by the Animal Ethics Committee for the NSW Department of Primary Industries, Elizabeth Macarthur Agricultural Institute, Menangle, Australia (M20-02, approved 9 February 2020).

Informed Consent Statement

Not applicable.

Data Availability Statement

All supporting data are included in this manuscript.

Acknowledgments

We wish to acknowledge the assistance of staff members at Elizabeth Macarthur Agricultural Institute and at CSIRO, Acton, ACT for their help in this project: Sylvia Vrankovic, Nina Huang, Tegan King, Sammi Chong, Ros Mourant, and Megan Pavy for undertaking the serological testing, and Kim Koeford, Nataly Boskoski and Sarah Moore for their assistance in the care and sampling of the rabbits. Bernie Dominiak for his helpful comments on this manuscript.

Conflicts of Interest

The authors declare no conflict of interest. The funding body had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Cooke, B.; Chudleigh, P.; Simpson, S.; Saunders, G. The economic benefits of the biological control of rabbits in Australia, 1950–2011. Aust. Econ. Hist. Rev. 2013, 53, 91–107. [Google Scholar] [CrossRef]
  2. Le Pendu, J.; Abrantes, J.; Bertagnoli, S.; Guitton, J.S.; Le Gall-Recule, G.; Lopes, A.M.; Marchandeau, S.; Alda, F.; Almeida, T.; Celio, A.P.; et al. Proposal for a unified classification system and nomenclature of lagoviruses. J. Gen. Virol. 2017, 98, 1658–1666. [Google Scholar] [CrossRef] [PubMed]
  3. Cox, T.E.; Ramsey, D.S.; Sawyers, E.; Campbell, S.; Matthews, J.; Elsworth, P. The impact of RHDV-K5 on rabbit populations in Australia: An evaluation of citizen science surveys to monitor rabbit abundance. Sci. Rep. 2019, 9, 15229. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Strive, T.; Cox, T.E. Lethal biological control of rabbits—The most powerful tools for landscape-scale mitigation of rabbit impacts in Australia. Aust. Zool. 2019, 40, 118–128. [Google Scholar] [CrossRef]
  5. Mahar, J.E.; Jenckel, M.; Huang, N.; Smertina, E.; Holmes, E.C.; Strive, T.; Hall, R.N. Frequent intergenotypic recombination between the non-structural and structural genes is a major driver of epidemiological fitness in caliciviruses. Virus Evol. 2021, 7, veab080. [Google Scholar] [CrossRef] [PubMed]
  6. Strive, T.; Wright, J.D.; Robinson, A.J. Identification and partial characterisation of a new lagovirus in Australian wild rabbits. Virology 2009, 384, 97–105. [Google Scholar] [CrossRef] [Green Version]
  7. Hall, R.N.; Mahar, J.E.; Haboury, S.; Stevens, V.; Holmes, E.C.; Strive, T. Emerging Rabbit Hemorrhagic Disease Virus 2 (RHDVb), Australia. Emerg. Infect. Dis. 2015, 21, 2276–2278. [Google Scholar] [CrossRef]
  8. Taggart, P.L.; Hall, R.N.; Cox, T.E.; Kovaliski, J.; McLeod, S.R.; Strive, T. Changes in virus transmission dynamics following the emergence of RHDV2 shed light on its competitive advantage over previously circulating variants. Transbound. Emerg. Dis. 2021. online ahead of print. [Google Scholar] [CrossRef]
  9. Mutze, G.; Cooke, B.; Alexander, P. The initial impact of rabbit hemorrhagic disease on European rabbit populations in South Australia. J. Wildl. Dis. 1998, 34, 221–227. [Google Scholar] [CrossRef] [Green Version]
  10. Cooke, B.; Fenner, F. Rabbit haemorrhagic disease and the biological control of wild rabbits, Oryctolagus cuniculus, in Australia and New Zealand. Wildl. Res. 2002, 29, 689–706. [Google Scholar] [CrossRef]
  11. Read, A.; Kirkland, P. Efficacy of a commercial vaccine against different strains of rabbit haemorrhagic disease virus. Aust. Vet. J. 2017, 95, 223–226. [Google Scholar] [CrossRef] [PubMed]
  12. Le Gall-Recule, G.; Lavazza, A.; Marchandeau, S.; Bertagnoli, S.; Zwingelstein, F.; Cavadini, P.; Martinelli, N.; Lombardi, G.; Guerin, J.L.; Lemaitre, E.; et al. Emergence of a new lagovirus related to Rabbit Haemorrhagic Disease Virus. Vet. Res. 2013, 44, 81. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Torres, J.M.; Ramírez, M.A.; Morales, M.; Bárcena, J.; Vázquez, B.; Espuña, E.; Pagès-Manté, A.; Sánchez-Vizcaíno, J.M. Safety evaluation of a recombinant myxoma-RHDV virus inducing horizontal transmissible protection against myxomatosis and rabbit haemorrhagic disease. Vaccine 2000, 19, 174–182. [Google Scholar] [CrossRef]
  14. Animal Medicines Australia. Pets in Australia: A National Survey of Pets and People. Available online: https://animalmedicinesaustralia.org.au/wp-content/uploads/2019/10/ANIM001-Pet-Survey-Report19_v1.7_WEB_high-res.pdf (accessed on 25 August 2021).
  15. Calvete, C.; Mendoza, M.; Alcaraz, A.; Sarto, M.P.; Jiménez-de-Bagüéss, M.P.; Calvo, A.J.; Monroy, F.; Calvo, J.H. Rabbit haemorrhagic disease: Cross-protection and comparative pathogenicity of GI.2/RHDV2/b and GI.1b/RHDV lagoviruses in a challenge trial. Vet. Microbiol. 2018, 219, 87–95. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Strive, T.; Elsworth, P.; Liu, J.; Wright, J.D.; Kovaliski, J.; Capucci, L. The non-pathogenic Australian rabbit calicivirus RCV-A1 provides temporal and partial cross protection to lethal Rabbit Haemorrhagic Disease Virus infection which is not dependent on antibody titres. Vet. Res. 2013, 44, 51. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Müller, C.; Hrynkiewicz, R.; Bębnowska, D.; Maldonado, J.; Baratelli, M.; Köllner, B.; Niedźwiedzka-Rystwej, P. Immunity against Lagovirus europaeus and the Impact of the Immunological Studies on Vaccination. Vaccines 2021, 9, 255. [Google Scholar] [CrossRef]
  18. Strive, T.; Piper, M.; Huang, N.; Mourant, R.; Kovaliski, J.; Capucci, L.; Cox, T.E.; Smith, I. Retrospective serological analysis reveals presence of the emerging lagovirus RHDV2 in Australia in wild rabbits at least five months prior to its first detection. Transbound. Emerg. Dis. 2020, 67, 822–833. [Google Scholar] [CrossRef]
  19. Mahar, J.E.; Hall, R.N.; Peacock, D.; Kovaliski, J.; Piper, M.; Mourant, R.; Huang, N.; Campbell, S.; Gu, X.; Read, A. Rabbit hemorrhagic disease virus 2 (RHDV2; GI. 2) is replacing endemic strains of RHDV in the Australian landscape within 18 months of its arrival. J. Virol. 2018, 92, e01374-17. [Google Scholar] [CrossRef] [Green Version]
  20. Hall, R.N.; King, T.; Connor, T.; Read, A.J.; Arrow, J.; Trought, K.; Duckworth, J.; Piper, M.; Strive, T. Age and Infectious Dose Significantly Affect Disease Progression after RHDV2 Infection in Naïve Domestic Rabbits. Viruses 2021, 13, 1184. [Google Scholar] [CrossRef]
  21. Hall, R.N.; King, T.; O’Connor, T.W.; Read, A.J.; Vrankovic, S.; Piper, M.; Strive, T. Passive Immunisation against RHDV2 Induces Protection against Disease but Not Infection. Vaccines 2021, 9, 1197. [Google Scholar] [CrossRef]
  22. Bahnemann, H.G. Binary ethylenimine as an inactivant for foot-and-mouth disease virus and its application for vaccine production. Arch. Virol. 1975, 47, 47–56. [Google Scholar] [CrossRef] [PubMed]
  23. Abrantes, J.; Lopes, A.M. A Review on the Methods Used for the Detection and Diagnosis of Rabbit Hemorrhagic Disease Virus (RHDV). Microorganisms 2021, 9, 972. [Google Scholar] [CrossRef] [PubMed]
  24. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020; Available online: https://www.R-project.org/ (accessed on 1 September 2021).
  25. Kassambara, A.; Kosinski, M.; Biecek, P. Survminer: Drawing Survival Curves Using ‘Ggplot2’. 2020. Available online: https://CRAN.R-project.org/package=survminer (accessed on 1 September 2021).
  26. Therneau, T. A Package for Survival Analysis in R; R Package Version 3.2-13; 2021; Available online: https://CRAN.R-project.org/package=survival (accessed on 1 September 2021).
  27. Liu, J.; Kerr, P.J.; Strive, T. A sensitive and specific blocking ELISA for the detection of rabbit calicivirus RCV-A1 antibodies. Virol. J. 2012, 9, 182. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Collins, B.J.; White, J.R.; Lenghaus, C.; Morrissy, C.J.; Westbury, H.A. Presence of rabbit haemorrhagic disease virus antigen in rabbit tissues as revealed by a monoclonal antibody dependent capture ELISA. J. Virol. Methods 1996, 58, 145–154. [Google Scholar] [CrossRef]
  29. Hall, R.N.; Mahar, J.E.; Read, A.J.; Mourant, R.; Piper, M.; Huang, N.; Strive, T. A strain-specific multiplex RT-PCR for Australian rabbit haemorrhagic disease viruses uncovers a new recombinant virus variant in rabbits and hares. Transbound. Emerg. Dis. 2018, 65, e444–e456. [Google Scholar] [CrossRef]
  30. Gall, A.; Hoffmann, B.; Teifke, J.P.; Lange, B.; Schirrmeier, H. Persistence of viral RNA in rabbits which overcome an experimental RHDV infection detected by a highly sensitive multiplex real-time RT-PCR. Vet. Microbiol. 2007, 120, 17–32. [Google Scholar] [CrossRef]
  31. Qin, S.; Underwood, D.; Driver, L.; Kistler, C.; Diallo, I.; Kirkland, P.D. Evaluation of a duplex reverse-transcription real-time PCR assay for the detection of encephalomyocarditis virus. J. Vet. Diagn. Investig. 2018, 30, 554–559. [Google Scholar] [CrossRef] [Green Version]
  32. Peacock, D.; Kovaliski, J.; Sinclair, R.; Mutze, G.; Iannella, A.; Capucci, L. RHDV2 overcoming RHDV immunity in wild rabbits (Oryctolagus cuniculus) in Australia. Vet. Rec. 2017, 180, 280. [Google Scholar] [CrossRef]
  33. Baratelli, M.; Molist-Badiola, J.; Puigredon-Fontanet, A.; Pascual, M.; Boix, O.; Mora-Igual, F.X.; Woodward, M.; Lavazza, A.; Capucci, L. Characterization of the Maternally Derived Antibody Immunity against Rhdv-2 after Administration in Breeding Does of an Inactivated Vaccine. Vaccines 2020, 8, 484. [Google Scholar] [CrossRef]
  34. Robinson, A.J.; So, P.T.M.; Müller, W.J.; Cooke, B.D.; Capucci, L. Statistical models for the effect of age and maternal antibodies on the development of rabbit haemorrhagic disease in Australian wild rabbits. Wildl. Res. 2002, 29, 663–671. [Google Scholar] [CrossRef]
  35. Calvete, C.; Capucci, L.; Lavazza, A.; Sarto, M.P.; Calvo, A.J.; Monroy, F.; Calvo, J.H. Changes in European wild rabbit population dynamics and the epidemiology of rabbit haemorrhagic disease in response to artificially increased viral transmission. Transbound. Emerg. Dis. 2021. online ahead of print. [Google Scholar] [CrossRef] [PubMed]
  36. Bárcena, J.; Guerra, B.; Angulo, I.; González, J.; Valcárcel, F.; Mata, C.P.; Castón, J.R.; Blanco, E.; Alejo, A. Comparative analysis of rabbit hemorrhagic disease virus (RHDV) and new RHDV2 virus antigenicity, using specific virus-like particles. Vet. Res. 2015, 46, 106. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Cox, R.J.; Brokstad, K.A.; Ogra, P. Influenza Virus: Immunity and Vaccination Strategies. Comparison of the Immune Response to Inactivated and Live, Attenuated Influenza Vaccines. Scand. J. Immunol. 2004, 59, 1–15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. (A) Schematic diagram for the production and challenge of rabbits immune to GI.1a or GI.4c. Rabbits with antibodies against GI.1a or GI.4c were produced by initially infecting by the oral route of infection at five weeks of age. Rabbits were subsequently challenged orally with GI.1bP-GI.2 at 12 weeks of age. (B). Schematic diagram for the production and challenge of rabbits vaccinated with a prototype GI.2-specific vaccine or a multivalent RHDV vaccine. Rabbits were vaccinated via subcutaneous injection at 10–12 weeks of age, and subsequently challenged per os with heterologous RHDV variants 28 days after vaccination. (C) Schematic diagram for the challenge of rabbits with naturally acquired immunity from three commercial farms. Rabbits were challenged with various heterologous RHDV variants (at various ages and with various infectious doses).
Figure 1. (A) Schematic diagram for the production and challenge of rabbits immune to GI.1a or GI.4c. Rabbits with antibodies against GI.1a or GI.4c were produced by initially infecting by the oral route of infection at five weeks of age. Rabbits were subsequently challenged orally with GI.1bP-GI.2 at 12 weeks of age. (B). Schematic diagram for the production and challenge of rabbits vaccinated with a prototype GI.2-specific vaccine or a multivalent RHDV vaccine. Rabbits were vaccinated via subcutaneous injection at 10–12 weeks of age, and subsequently challenged per os with heterologous RHDV variants 28 days after vaccination. (C) Schematic diagram for the challenge of rabbits with naturally acquired immunity from three commercial farms. Rabbits were challenged with various heterologous RHDV variants (at various ages and with various infectious doses).
Vaccines 10 00666 g001
Figure 2. Survival curve of seronegative rabbits after challenge with GI.1c (red), GI.1a (yellow) or two different doses of GI.2 (purple-dashed, dark purple). For seronegative rabbits, there was no statistical difference in survival between the variant used nor the infectious dose.
Figure 2. Survival curve of seronegative rabbits after challenge with GI.1c (red), GI.1a (yellow) or two different doses of GI.2 (purple-dashed, dark purple). For seronegative rabbits, there was no statistical difference in survival between the variant used nor the infectious dose.
Vaccines 10 00666 g002
Figure 3. Survival curve of rabbits with experimentally acquired GI.1a (yellow) or GI.4c (green) immunity compared to seronegative (grey) rabbits after a GI.2 challenge (50 RID50 PO). Immunity with either GI.1a or GI.4c variants provided significant cross-protection against GI.2 challenge.
Figure 3. Survival curve of rabbits with experimentally acquired GI.1a (yellow) or GI.4c (green) immunity compared to seronegative (grey) rabbits after a GI.2 challenge (50 RID50 PO). Immunity with either GI.1a or GI.4c variants provided significant cross-protection against GI.2 challenge.
Vaccines 10 00666 g003
Figure 4. Survival curve of rabbits given a GI.2-specific vaccine (AC, magenta), multivalent RHDV vaccine (DF, blue), or no vaccine (AF, grey) and challenged with 1500 RID50 of GI.1a, GI.1c or GI.2. While the GI.2-specific vaccine provided complete protection against GI.2 challenge, only the rabbits given the multivalent RHDV vaccine were protected against GI.1a, GI.1c and GI.2 challenge.
Figure 4. Survival curve of rabbits given a GI.2-specific vaccine (AC, magenta), multivalent RHDV vaccine (DF, blue), or no vaccine (AF, grey) and challenged with 1500 RID50 of GI.1a, GI.1c or GI.2. While the GI.2-specific vaccine provided complete protection against GI.2 challenge, only the rabbits given the multivalent RHDV vaccine were protected against GI.1a, GI.1c and GI.2 challenge.
Vaccines 10 00666 g004
Figure 5. (A) Survival curve of rabbits with naturally acquired GI.1 immunity challenged with GI.1a (yellow), or with various doses of GI.2 (150 RID50—dotted-purple, 1500 RID50—purple) at different ages (12 weeks—purple dashed line, 33 weeks—purple solid line). All rabbits survived GI.1a challenge. A dose-dependent and age-related difference in survival was observed against GI.2 challenge, with rabbits given a higher dose or challenged at an older age less likely to survive. (B) Survival curve of rabbits with naturally acquired GI.2 immunity after GI.1c (red), GI.1a (yellow) or GI.2 (purple) challenge. Fewer rabbits survived a GI.1a challenge compared to the other challenge viruses. (C) Survival curve of rabbits with RHDV (GI.4cP-GI.2) immunity after heterologous GI.2 (purple) or GI.1a challenge at various doses (150 RID50—yellow-dashed, 1500 RID50—yellow-solid). All rabbits survived GI.2 challenge, whereas fewer rabbits survived GI.1a challenge.
Figure 5. (A) Survival curve of rabbits with naturally acquired GI.1 immunity challenged with GI.1a (yellow), or with various doses of GI.2 (150 RID50—dotted-purple, 1500 RID50—purple) at different ages (12 weeks—purple dashed line, 33 weeks—purple solid line). All rabbits survived GI.1a challenge. A dose-dependent and age-related difference in survival was observed against GI.2 challenge, with rabbits given a higher dose or challenged at an older age less likely to survive. (B) Survival curve of rabbits with naturally acquired GI.2 immunity after GI.1c (red), GI.1a (yellow) or GI.2 (purple) challenge. Fewer rabbits survived a GI.1a challenge compared to the other challenge viruses. (C) Survival curve of rabbits with RHDV (GI.4cP-GI.2) immunity after heterologous GI.2 (purple) or GI.1a challenge at various doses (150 RID50—yellow-dashed, 1500 RID50—yellow-solid). All rabbits survived GI.2 challenge, whereas fewer rabbits survived GI.1a challenge.
Vaccines 10 00666 g005
Table 1. Immunity status, age at challenge, challenge virus, infectious dose and the subsequent survival proportion of rabbits with experimentally or naturally acquired immunity.
Table 1. Immunity status, age at challenge, challenge virus, infectious dose and the subsequent survival proportion of rabbits with experimentally or naturally acquired immunity.
Immunity StatusAge at Challenge (Weeks)Challenge VirusInfectious Dose (RID50 1)Survived/Total
Seronegative11GI.2 2501/8 3
Seronegative16GI.1a15000/12
Seronegative16GI.1c15000/12
Seronegative16GI.2215000/12
GI.4c11GI.22507/9
GI.1a11GI.22507/7
GI.1 412GI.221509/10
GI.1 412GI.2215005/9
GI.1 433GI.2215001/7
GI.1 433GI.1a15007/7
GI.210–12GI.1a15001/4
GI.210–12GI.1c15002/3
GI.210–12GI.2215003/4
GI.2 512GI.1a1509/12
GI.2 512GI.1a15009/13
GI.2 512GI.2 2150,00012/12
1 The 50% rabbit infectious dose (RID50). 2 This GI.2 refers to a naturally occurring RHDV2 recombinant virus designated GI.1bP-GI.2 (GenBank acc. Number MW467791). 3 This one surviving rabbit did not seroconvert after challenge and was not infected with GI.2. 4 These rabbits could have antibodies from either GI.1a or GI.1c exposure. 5 This GI.2 refers to a naturally occurring RHDV2 recombinant virus designated GI.4cP-GI.2 (GenBank acc. Number MW460156).
Table 2. Proportion of vaccinated rabbits that succumbed to rabbit hemorrhagic disease following oral inoculation with Australian lagoviruses.
Table 2. Proportion of vaccinated rabbits that succumbed to rabbit hemorrhagic disease following oral inoculation with Australian lagoviruses.
Vaccination StatusVaccine Dose (HAU 1)Challenge Virus (Dose 1500 RID50)Survived/Total
SeronegativeNo vaccineGI.1c0/12
SeronegativeNo vaccineGI.1a0/12
SeronegativeNo vaccineGI.2 20/12
GI.2 2128GI.1c0/6
GI.2 2128GI.1a0/6
GI.2 2100GI.2 29/9
Multivalent96GI.1c6/6
Multivalent96GI.1a6/6
Multivalent96GI.2 26/6
1 Total hemagglutination units. 2 This GI.2 refers to a RHDV2 recombinant designated GI.1bP-GI.2.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

O’Connor, T.W.; Read, A.J.; Hall, R.N.; Strive, T.; Kirkland, P.D. Immunological Cross-Protection between Different Rabbit Hemorrhagic Disease Viruses—Implications for Rabbit Biocontrol and Vaccine Development. Vaccines 2022, 10, 666. https://doi.org/10.3390/vaccines10050666

AMA Style

O’Connor TW, Read AJ, Hall RN, Strive T, Kirkland PD. Immunological Cross-Protection between Different Rabbit Hemorrhagic Disease Viruses—Implications for Rabbit Biocontrol and Vaccine Development. Vaccines. 2022; 10(5):666. https://doi.org/10.3390/vaccines10050666

Chicago/Turabian Style

O’Connor, Tiffany W., Andrew J. Read, Robyn N. Hall, Tanja Strive, and Peter D. Kirkland. 2022. "Immunological Cross-Protection between Different Rabbit Hemorrhagic Disease Viruses—Implications for Rabbit Biocontrol and Vaccine Development" Vaccines 10, no. 5: 666. https://doi.org/10.3390/vaccines10050666

APA Style

O’Connor, T. W., Read, A. J., Hall, R. N., Strive, T., & Kirkland, P. D. (2022). Immunological Cross-Protection between Different Rabbit Hemorrhagic Disease Viruses—Implications for Rabbit Biocontrol and Vaccine Development. Vaccines, 10(5), 666. https://doi.org/10.3390/vaccines10050666

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop