# Comment on Shah et al. Genetic Characteristics and Phylogeographic Dynamics of Lagoviruses

Joana Abrantes, Stéphan Bertagnoli, Patriz Cavadini, Pedro J Estev, Dolor Gavier-Widén, Robyn N Hall, Anton Lavazza, Ghislai Le Gall-Reculé, Jackie E Mahar, Stéphane Marchandeau, et al.

► To cite this version:

Joana Abrantes, Stéphan Bertagnoli, Patriz Cavadini, Pedro J Estev, Dolor Gavier-Widén, et al. Comment on Shah et al. Genetic Characteristics and Phylogeographic Dynamics of Lagoviruses. Viruses, 2024, 10.3390/v16060927. [https://doi.org/10.3390/v16060927](https://doi.org/10.3390/v16060927)

HAL Id: hal-04604804

https://ofb.hal.science/hal-04604804

Submitted on 7 Jun 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Comment on Shah et al. Genetic Characteristics and Phylogeographic Dynamics of Lagoviruses

Joana Abrantes, Stéphan Bertagnoli, Patriz Cavadini, Pedro J. Estev, Dolor Gavier-Widén, Robyn N. Hall, Anton Lavazza, Ghislaine Le Gall-Reculé, Jackie E. Mahar, Stéphane Marchandeau, et al.

<figure><img src="image_1.png" /></figure>

To cite this version:

Joana Abrantes, Stéphan Bertagnoli, Patriz Cavadini, Pedro J. Estev, Dolor Gavier-Widén, et al. Comment on Shah et al. Genetic Characteristics and Phylogeographic Dynamics of Lagoviruses. Viruses, 2024, 10.3390/v16060927. [https://doi.org/10.3390/v16060927](https://doi.org/10.3390/v16060927)

HAL Id: hal-04604804

https://ofb.hal.science/hal-04604804v1

Submitted on 7 Jun 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Comment

Comment on Shah et al. Genetic Characteristics and Phylogeographic Dynamics of Lagoviruses, 1988–2021. *Viruses* 2023, 15, 815 †

Joana Abrantes 1,2,3 Ⓔ, Stéphane Bertagnoli 4, Patrizia Cavadini 5,6, Pedro J. Esteves 1,2,3,7 Ⓔ, Dolores Gavier-Widén 8,9 Ⓔ, Robyn N. Hall 10,11,12 Ⓔ, Antonio Lavazza 5,6 Ⓔ, Ghislaine Le Gall-Reculé 13 Ⓔ, Jackie E. Mahar 14 Ⓔ, Stéphane Marchandeau 15 Ⓔ and Ana M. Lopes 1,2,16,17,* Ⓔ

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Abrantes, J.; Bertagnoli, S.; Cavadini, P.; Esteves, P.J.; Gavier-Widén, D.; Hall, R.N.; Lavazza, A.; Le Gall-Reculé, G.; Mahar, J.E.; Marchandeau, S.; et al. Comment on Shah et al. Genetic Characteristics and Phylogeographic Dynamics of Lagoviruses, 1988–2021. *Viruses* 2023, 15, 815. *Viruses* **2024**, *16*, 927. https://doi.org/10.3390/v16060927

Received: 16 February 2024

Revised: 22 May 2024

Accepted: 31 May 2024

Published: 7 June 2024

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2024 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

**Keywords:** *Lagovirus europaeus*; phylogeny; recombination; evolutionary analysis

Shah and colleagues [1] published an article in *Viruses* where they performed a phylogenetic and phylogeographic analysis of full-length nucleotide sequences of lagoviruses (N = 240). The analysis included pathogenic and non-pathogenic forms of typically rabbit-associated lagoviruses (rabbit hemorrhagic disease virus, RHDV, and rabbit calicivirus, RCV, both belonging to GI according to the nomenclature proposed by Le Pendu et al. [2]) and hare-associated lagoviruses (European brown hare syndrome virus, EBHSV, and hare calicivirus, HaCV, both assigned to GII). According to the authors, the goal of the study was to provide “an update of the phylogenetic and phylogeographic information of lagoviruses that may be used to map the evolutionary history and provide hints for the genetic basis of their emergence and re-emergence”. In our opinion, and despite the authors’ effort in

<footnote n="1">¹ CIBIO, Centro de Investigação em Biodiversidade e Recursos Genéticos, InBIO Laboratório Associado, Campus de Vairão, Universidade do Porto, 4485-661 Vairão, Portugal; jabrantes@cibio.up.pt (J.A.); pjesteves@cibio.up.pt (P.J.E.)</footnote>

<footnote n="2">² BIOPOLIS Program in Genomics, Biodiversity and Land Planning, CIBIO, Campus de Vairão, Universidade do Porto, 4485-661 Vairão, Portugal</footnote>

<footnote n="3">³ Departamento de Biologia, Faculdade de Ciências, Universidade do Porto, 4099-002 Porto, Portugal</footnote>

<footnote n="4">⁴ Laboratoire Interactions Hôtes-Agents Pathogènes, Université de Toulouse, INRAE, ENVT, CEDEX 3, 31076 Toulouse, France; stephane.bertagnoli@envt.fr</footnote>

<footnote n="5">5 Istituto Zooprofilattico Sperimentale della Lombardia e dell'Emilia Romagna, 25124 Brescia, Italy; patrizia.cavadini@izsler.it (P.C.); antonio.lavazza@izsler.it (A.L.)</footnote>

<footnote n="6">6 WOAH Reference Laboratory for Rabbit Haemorrhagic Disease, Via Bianchi 7/9, 25124 Brescia, Italy</footnote>

<footnote n="7">7 CITS—Center of Investigation in Health Technologies, CESPU, 4585-116 Gandra, Portugal</footnote>

<footnote n="8">8 Swedish Veterinary Agency (SVA), 75189 Uppsala, Sweden; dolores.gavier-widen@sva.se</footnote>

<footnote n="9">⁹ Department of Biomedical Sciences and Veterinary Public Health, Swedish University of Agricultural Sciences (SLU), Box 7028, 75007 Uppsala, Sweden</footnote>

<footnote n="10">¹⁰ Commonwealth Scientific and Industrial Research Organisation, Health and Biosecurity, Canberra, ACT 2601, Australia; robyn.hall@ausvet.com.au</footnote>

<footnote n="11">¹¹ Centre for Invasive Species Solutions, Bruce, ACT 2617, Australia</footnote>

<footnote n="12">12 Ausvet Pty Ltd., Canberra, ACT 2617, Australia</footnote>

<footnote n="13">13 Ploufragan-Plouzané-Niort Laboratory, Avian & Rabbit Virology, Immunology & Parasitology Unit, French Agency for Food, Environmental and Occupational Health and Safety (Anses), 22440 Ploufragan, France; ghlislaine.legall-recule@anses.fr</footnote>

<footnote n="14">14 Commonwealth Scientific and Industrial Research Organisation, Australian Animal Health Laboratory and Health and Biosecurity, Geelong, VIC 3220, Australia; Jackie.mahar@csiro.au</footnote>

<footnote n="15">15 French Biodiversity Agency (OFB), 94300 Nantes, France; stephane.marchandeau@ofb.gouv.fr</footnote>

<footnote n="16">16 UMIB—Unit for Multidisciplinary Research in Biomedicine, ICBAS—School of Medicine and Biomedical Sciences, University of Porto, 4050-313 Porto, Portugal</footnote>

<footnote n="17">17 ITR—Laboratory for Integrative and Translational Research in Population Health, 4050-600 Porto, Portugal</footnote>

<footnote>* Correspondence: analopes@cibio.up.pt</footnote>

<footnote>† Authors are listed in alphabetic order except for the first and senior authors.</footnote>
Comment

# Comment on Shah et al. Genetic Characteristics and Phylogeographic Dynamics of Lagoviruses, 1988–2021. *Viruses* 2023, 15, 815 †

Joana Abrantes 1,2,3 Ⓒ, Stéphane Bertagnoli 4, Patrizia Cavadini 5,6, Pedro J. Esteves 1,2,3,7 Ⓒ, Dolores Gavier-Widén 8,9 Ⓒ, Robyn N. Hall 10,11,12 Ⓒ, Antonio Lavazza 5,6 Ⓒ, Ghislaine Le Gall-Reculé 13 Ⓒ, Jackie E. Mahar 14 Ⓒ, Stéphane Marchandeau 15 Ⓒ and Ana M. Lopes 1,2,16,17,* Ⓒ

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Abrantes, J.; Bertagnoli, S.; Cavadini, P.; Esteves, P.J.; Gavier-Widén, D.; Hall, R.N.; Lavazza, A.; Le Gall-Reculé, G.; Mahar, J.E.; Marchandeau, S.; et al. Comment on Shah et al. Genetic Characteristics and Phylogeographic Dynamics of Lagoviruses, 1988–2021. *Viruses* 2023, 15, 815. *Viruses* **2024**, *16*, 927. https://doi.org/10.3390/v16060927

Received: 16 February 2024

Revised: 22 May 2024

Accepted: 31 May 2024

Published: 7 June 2024

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2024 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

**Keywords:** *Lagovirus europaeus*; phylogeny; recombination; evolutionary analysis

Shah and colleagues [1] published an article in *Viruses* where they performed a phylogenetic and phylogeographic analysis of full-length nucleotide sequences of lagoviruses (N = 240). The analysis included pathogenic and non-pathogenic forms of typically rabbit-associated lagoviruses (rabbit hemorrhagic disease virus, RHDV, and rabbit calicivirus, RCV, both belonging to GI according to the nomenclature proposed by Le Pendu et al. [2]) and hare-associated lagoviruses (European brown hare syndrome virus, EBHSV, and hare calicivirus, HaCV, both assigned to GII). According to the authors, the goal of the study was to provide “an update of the phylogenetic and phylogeographic information of lagoviruses that may be used to map the evolutionary history and provide hints for the genetic basis of their emergence and re-emergence”. In our opinion, and despite the authors’ effort in

<footnote n="1">¹ CIBIO, Centro de Investigação em Biodiversidade e Recursos Genéticos, InBIO Laboratório Associado, Campus de Vairão, Universidade do Porto, 4485-661 Vairão, Portugal; jabrantes@cibio.up.pt (J.A.); pjesteves@cibio.up.pt (P.J.E.)</footnote>

<footnote n="2">² BIOPOLIS Program in Genomics, Biodiversity and Land Planning, CIBIO, Campus de Vairão, Universidade do Porto, 4485-661 Vairão, Portugal</footnote>

<footnote n="3">³ Departamento de Biologia, Faculdade de Ciências, Universidade do Porto, 4099-002 Porto, Portugal</footnote>

<footnote n="4">⁴ Laboratoire Interactions Hôtes-Agents Pathogènes, Université de Toulouse, INRAE, ENVT, CEDEX 3, 31076 Toulouse, France; stephane.bertagnoli@envt.fr</footnote>

<footnote n="5">5 Istituto Zooprofilattico Sperimentale della Lombardia e dell'Emilia Romagna, 25124 Brescia, Italy; patrizia.cavadini@izsler.it (P.C.); antonio.lavazza@izsler.it (A.L.)</footnote>

<footnote n="6">6 WOAH Reference Laboratory for Rabbit Haemorrhagic Disease, Via Bianchi 7/9, 25124 Brescia, Italy</footnote>

<footnote n="7">7 CITS—Center of Investigation in Health Technologies, CESPU, 4585-116 Gandra, Portugal</footnote>

<footnote n="8">8 Swedish Veterinary Agency (SVA), 75189 Uppsala, Sweden; dolores.gavier-widen@sva.se</footnote>

<footnote n="9">⁹ Department of Biomedical Sciences and Veterinary Public Health, Swedish University of Agricultural Sciences (SLU), Box 7028, 75007 Uppsala, Sweden</footnote>

<footnote n="10">¹⁰ Commonwealth Scientific and Industrial Research Organisation, Health and Biosecurity, Canberra, ACT 2601, Australia; robyn.hall@ausvet.com.au</footnote>

<footnote n="11">¹¹ Centre for Invasive Species Solutions, Bruce, ACT 2617, Australia</footnote>

<footnote n="12">¹² Ausvet Pty Ltd., Canberra, ACT 2617, Australia</footnote>

<footnote n="13">¹³ Ploufragan-Plouzané-Niort Laboratory, Avian & Rabbit Virology, Immunology & Parasitology Unit, French Agency for Food, Environmental and Occupational Health and Safety (Anses), 22440 Ploufragan, France; ghlislaine.legall-recule@anses.fr</footnote>

<footnote n="14">¹⁴ Commonwealth Scientific and Industrial Research Organisation, Australian Animal Health Laboratory and Health and Biosecurity, Geelong, VIC 3220, Australia; Jackie.mahar@csiro.au</footnote>

<footnote n="15">¹⁵ French Biodiversity Agency (OFB), 94300 Nantes, France; stephane.marchandeau@ofb.gouv.fr</footnote>

<footnote n="16">¹⁶ UMIB—Unit for Multidisciplinary Research in Biomedicine, ICBAS—School of Medicine and Biomedical Sciences, University of Porto, 4050-313 Porto, Portugal</footnote>

<footnote n="17">¹⁷ ITR—Laboratory for Integrative and Translational Research in Population Health, 4050-600 Porto, Portugal</footnote>

<footnote>* Correspondence: analopes@cibio.up.pt</footnote>

<footnote>† Authors are listed in alphabetic order except for the first and senior authors.</footnote>
contributing to the current knowledge of lagovirus evolution, some of the methodologies used in the analysis were inappropriate, leading to incorrect and misleading inferences.

Recombination is a major evolutionary force, well documented to occur at a high frequency in some RNA viruses [3–5]. It often results from the viral polymerase switching to a new template during recombination, frequently close to the original position due to sequence similarity/homology [6]. This reshuffling of viral genomes results in a mix of genes with different evolutionary histories, and the new genomic combinations might confer novel features to the resulting virus, such as higher pathogenic potential, evasion from the host immune system, resistance to antiviral therapies and vaccines, the expansion of host range (species jump), or the modification of cell or host tropism. As in other RNA viruses, recombination has been described for lagoviruses, with a major recombination hotspot located at the boundary between the non-structural and structural gene-encoding regions and other not-so-frequent breakpoints scattered throughout the genome e.g., [7–13]. Notably, recombination seems to have contributed to the emergence of GI.2 in 2010 by combining the non-structural proteins of a non-pathogenic GI.3 [2] with the structural proteins of a novel virus, GI.2 [14]. This novel virus was also shown to have recombined with pathogenic GI.1a and GI.1b, as well as non-pathogenic GI.4 [10,15], all contributing to the pool of circulating lagoviruses. Hares have also been a source of new viruses, with both pathogenic and non-pathogenic GII recombining with GI.2 [7,13]. Nonetheless, recombination is not restricted to this genotype, and other key examples of the importance of this mechanism in lagoviruses include the combination of non-pathogenic with classical pathogenic strains, such as the GI.4eP-GI.1a and the GI.3P-GI.1d [12,16] that were able to spread and cause outbreaks in Australia and Europe, respectively.

Evolutionary analyses on recombinant lagoviruses conducted so far considered the impact of recombination in phylogenetic inference. Indeed, in recombinant strains, as the genomic regions on either side of a recombination breakpoint have different evolutionary origins, their phylogenetic history needs to be inferred independently [17,18]. A single phylogenetic tree inferred using the full genome sequence will not reflect the true evolutionary history of a recombinant strain, and no robust conclusions on branch length or topology can be drawn from such analysis. The impact of recombination on tree topology is robustly demonstrated in a recent paper by Zaman and colleagues [19] (see Figure 2 of that paper), which concludes that a reliable phylogeny can only be inferred from a non-recombinant or minimally recombinant region. Moreover, this finding has been thoroughly demonstrated in previous studies [17,20,21], including for other caliciviruses [22], and has motivated, for instance, the currently accepted classification for noroviruses [23,24]. This classification accounts for recombination by using a dual typing system to reflect the origin of both the polymerase and major capsid genes (which flank the recombination hotspot), clearly distinguishing wild-type strains from recombinants, which are “distinct epidemiological entities”. This system was further used as the basis of the nomenclature for lagoviruses proposed by Le Pendu et al. [2] that has been broadly adopted by those studying lagoviruses (e.g., [3,14]).

While it has been demonstrated that uncommon recombination events might be “accommodated” in phylogenetic inference [25], the high frequency of recombination events throughout lagoviruses evolution, particularly in GI.2, is likely to distort conclusions inferred from a phylogeny based on full-length genomes. In addition, by performing phylogenetic reconstruction with full-length genomes, recombination events, which evidently shape lagoviruses evolution, might be missed, and the complete epidemiological history obscured. Failure to consider recombination when inferring phylogenies also leads to overestimating the number of point mutations and the substitution rate, and the loss of the molecular clock [17], as has likely occurred in Shah et al. Their study presented a scale bar of 30.0 in the maximum likelihood phylogenies presented in Figure 1 (tree based on the full-length genome sequences) and in Supplementary Figures S1–S3 (trees based on full-length genome sequences, VP60 and VP10 coding sequences, respectively). Although the authors fail to indicate what the scale represents, the default assumption
contributing to the current knowledge of lagovirus evolution, some of the methodologies used in the analysis were inappropriate, leading to incorrect and misleading inferences.

Recombination is a major evolutionary force, well documented to occur at a high frequency in some RNA viruses [3–5]. It often results from the viral polymerase switching to a new template during recombination, frequently close to the original position due to sequence similarity/homology [6]. This reshuffling of viral genomes results in a mix of genes with different evolutionary histories, and the new genomic combinations might confer novel features to the resulting virus, such as higher pathogenic potential, evasion from the host immune system, resistance to antiviral therapies and vaccines, the expansion of host range (species jump), or the modification of cell or host tropism. As in other RNA viruses, recombination has been described for lagoviruses, with a major recombination hotspot located at the boundary between the non-structural and structural gene-encoding regions and other not-so-frequent breakpoints scattered throughout the genome e.g., [7–13]. Notably, recombination seems to have contributed to the emergence of GI.2 in 2010 by combining the non-structural proteins of a non-pathogenic GI.3 [2] with the structural proteins of a novel virus, GI.2 [14]. This novel virus was also shown to have recombined with pathogenic GI.1a and GI.1b, as well as non-pathogenic GI.4 [10,15], all contributing to the pool of circulating lagoviruses. Hares have also been a source of new viruses, with both pathogenic and non-pathogenic GII recombining with GI.2 [7,13]. Nonetheless, recombination is not restricted to this genotype, and other key examples of the importance of this mechanism in lagoviruses include the combination of non-pathogenic with classical pathogenic strains, such as the GI.4eP-GI.1a and the GI.3P-GI.1d [12,16] that were able to spread and cause outbreaks in Australia and Europe, respectively.

Evolutionary analyses on recombinant lagoviruses conducted so far considered the impact of recombination in phylogenetic inference. Indeed, in recombinant strains, as the genomic regions on either side of a recombination breakpoint have different evolutionary origins, their phylogenetic history needs to be inferred independently [17,18]. A single phylogenetic tree inferred using the full genome sequence will not reflect the true evolutionary history of a recombinant strain, and no robust conclusions on branch length or topology can be drawn from such analysis. The impact of recombination on tree topology is robustly demonstrated in a recent paper by Zaman and colleagues [19] (see Figure 2 of that paper), which concludes that a reliable phylogeny can only be inferred from a non-recombinant or minimally recombinant region. Moreover, this finding has been thoroughly demonstrated in previous studies [17,20,21], including for other caliciviruses [22], and has motivated, for instance, the currently accepted classification for noroviruses [23,24]. This classification accounts for recombination by using a dual typing system to reflect the origin of both the polymerase and major capsid genes (which flank the recombination hotspot), clearly distinguishing wild-type strains from recombinants, which are “distinct epidemiological entities”. This system was further used as the basis of the nomenclature for lagoviruses proposed by Le Pendu et al. [2] that has been broadly adopted by those studying lagoviruses (e.g., [3,14]).

While it has been demonstrated that uncommon recombination events might be “accommodated” in phylogenetic inference [25], the high frequency of recombination events throughout lagoviruses evolution, particularly in GI.2, is likely to distort conclusions inferred from a phylogeny based on full-length genomes. In addition, by performing phylogenetic reconstruction with full-length genomes, recombination events, which evidently shape lagoviruses evolution, might be missed, and the complete epidemiological history obscured. Failure to consider recombination when inferring phylogenies also leads to overestimating the number of point mutations and the substitution rate, and the loss of the molecular clock [17], as has likely occurred in Shah et al. Their study presented a scale bar of 30.0 in the maximum likelihood phylogenies presented in Figure 1 (tree based on the full-length genome sequences) and in Supplementary Figures S1–S3 (trees based on full-length genome sequences, VP60 and VP10 coding sequences, respectively). Although the authors fail to indicate what the scale represents, the default assumption
is that the scale bar represents 30.0 nucleotide substitutions/site, an unrealistic level of diversity, particularly given the length of the scale bar. This level of phylogenetic distance is incompatible with the current knowledge on RNA viruses' evolutionary rates of around $10^{-3}$ to $10^{-4}$ substitutions/site/year [5,26].

In their work, Shah and colleagues identified six recombination events: one in VP60, one in the non-structural protein RNA-dependent RNA polymerase (RdRp), and four near the boundary between RdRp and the capsid protein [1], the latter four being consistent with the hotspot of recombination previously described [3,7,10,13]. Despite acknowledging the existence of recombination in the dataset and having included additional VP60- and VP10-based phylogenies, their main conclusions were still inferred from the full-length phylogeny, which led to erroneous conclusions. For example, the authors classified the pathogenic isolate CBMad17-1 (GenBank accession number MF407655) as belonging to non-pathogenic RCV, but this virus has been demonstrated to have resulted from a recombination event between the non-structural region of a GI.4 (non-pathogenic) and the structural part of a GI.2 (pathogenic), i.e., GI.4P-GI.2, using the dual typing nomenclature [27]. Likewise, two GII.1P-GI.2 recombinant German variants (GenBank accession numbers LR899142 and LR899187) [7], critically important as the first examples of intergenogroup recombination in lagoviruses, were identified in Shah et al. as belonging to the HaCV/EBHSV cluster (i.e., GII hare-associated cluster). The authors further highlighted that they observed significant differences between the full-length genome and VP60-based phylogenies, in particular, for GI.1d/Ocun/FR/2009/09-03, P95 and GI.2/Ocun/FR/2013/13-165 (GenBank accession numbers MT628290, KJ943791, and MN737112, respectively), which is congruent with their recombinant origin previously described [12,14,28] and disregarded in this paper. On the other hand, several (previously reported) recombinant viruses remained apparently undetected by their methods, such as those associated with the first reported GI.2 outbreaks in 2010 (GenBank accession numbers KM878681, MN737113, MN737114, MN738377 and MN786321) [14] or the Hartmannsdorf isolate (GenBank accession number EF558586) [11,29]. Thus, it is our opinion that it is not correct to conclude (as they did) that "the full-length genomic sequence-based phylogeny provides a more reliable and robust classification of all the lagoviruses", and, indeed, we vehemently disagree with this statement.

The topology/clustering demonstrated in both the full-genome and VP60 phylogenies and the phylogeographic network analysis in Shah et al. is also a reason for concern. As stated in their phylogeographic analysis results section, "all HaCV/EBHSV strains clustered together and were sharing their ancestor with the GI.1 (RHDV) strains, (...) [and] all the RCV strains clustered with the GI.2 (RHDV2)" strains. Moreover, their maximum likelihood phylogeny grouped GI.4 (RCV) with GII (HaCV/EBHSV). This contradicts previous phylogenetic and phylogeographic inferences that consistently grouped lagoviruses into two main clusters, one comprising GI variants (RHDV and RCV) and the other composed of GII variants (EBHSV and HaCV) [2,7,30]. Indeed, there appears to be an overlap between phylogenetic clustering and biological host, with GII viruses (EBHSV and HaCV) being restricted to hare species and GI viruses (RHDV and RCV) being mostly found in rabbits. GI.2 is the exception as it has a wider host tropism, fatally infecting other European and American leporid genera (Lepus, Sylvilagus, and Brachylagus) [31–36], although it first emerged in rabbits (Oryctolagus cuniculus) [37]. Nonetheless, hare species seem to remain as dead-end hosts, and rabbits appear to have a higher susceptibility to lethal disease than other leporids, notwithstanding the limited availability of comparative studies [38]. Phylogenetic analysis in countless previous publications from different labs globally clearly confirm the separation between "rabbit lagoviruses" and "hare lagoviruses", e.g., [2,7,13,30,33,39]. Furthermore, this partition is confirmed by genetic distances between sequences from these groups (Figure 1).

In another example, their GI.2 clade includes viruses recovered from leporids found dead in 1994–1996 (GenBank accession numbers KY765609–11 and KJ943791) [28], long before the emergence of GI.2 in 2010 [40]. The evolutionary history of these four viruses
is that the scale bar represents 30.0 nucleotide substitutions/site, an unrealistic level of diversity, particularly given the length of the scale bar. This level of phylogenetic distance is incompatible with the current knowledge on RNA viruses' evolutionary rates of around $10^{-3}$ to $10^{-4}$ substitutions/site/year [5,26].

In their work, Shah and colleagues identified six recombination events: one in VP60, one in the non-structural protein RNA-dependent RNA polymerase (RdRp), and four near the boundary between RdRp and the capsid protein [1], the latter four being consistent with the hotspot of recombination previously described [3,7,10,13]. Despite acknowledging the existence of recombination in the dataset and having included additional VP60- and VP10-based phylogenies, their main conclusions were still inferred from the full-length phylogeny, which led to erroneous conclusions. For example, the authors classified the pathogenic isolate CBMad17-1 (GenBank accession number MF407655) as belonging to non-pathogenic RCV, but this virus has been demonstrated to have resulted from a recombination event between the non-structural region of a GI.4 (non-pathogenic) and the structural part of a GI.2 (pathogenic), i.e., GI.4P-GI.2, using the dual typing nomenclature [27]. Likewise, two GII.1P-GI.2 recombinant German variants (GenBank accession numbers LR899142 and LR899187) [7], critically important as the first examples of intergenogroup recombination in lagoviruses, were identified in Shah et al. as belonging to the HaCV/EBHSV cluster (i.e., GII hare-associated cluster). The authors further highlighted that they observed significant differences between the full-length genome and VP60-based phylogenies, in particular, for GI.1d/Ocun/FR/2009/09-03, P95 and GI.2/Ocun/FR/2013/13-165 (Gen-Bank accession numbers MT628290, KJ943791, and MN737112, respectively), which is congruent with their recombinant origin previously described [12,14,28] and disregarded in this paper. On the other hand, several (previously reported) recombinant viruses remained apparently undetected by their methods, such as those associated with the first reported GI.2 outbreaks in 2010 (GenBank accession numbers KM878681, MN737113, MN737114, MN738377 and MN786321) [14] or the Hartmannsdorf isolate (GenBank accession number EF558586) [11,29]. Thus, it is our opinion that it is not correct to conclude (as they did) that "the full-length genomic sequence-based phylogeny provides a more reliable and robust classification of all the lagoviruses", and, indeed, we vehemently disagree with this statement.

The topology/clustering demonstrated in both the full-genome and VP60 phylogenies and the phylogeographic network analysis in Shah et al. is also a reason for concern. As stated in their phylogeographic analysis results section, "all HaCV/EBHSV strains clustered together and were sharing their ancestor with the GI.1 (RHDV) strains, (...) [and] all the RCV strains clustered with the GI.2 (RHDV2)" strains. Moreover, their maximum likelihood phylogeny grouped GI.4 (RCV) with GII (HaCV/EBHSV). This contradicts previous phylogenetic and phylogeographic inferences that consistently grouped lagoviruses into two main clusters, one comprising GI variants (RHDV and RCV) and the other composed of GII variants (EBHSV and HaCV) [2,7,30]. Indeed, there appears to be an overlap between phylogenetic clustering and biological host, with GII viruses (EBHSV and HaCV) being restricted to hare species and GI viruses (RHDV and RCV) being mostly found in rabbits. GI.2 is the exception as it has a wider host tropism, fatally infecting other European and American leporid genera (Lepus, Sylvilagus, and Brachylagus) [31–36], although it first emerged in rabbits (Oryctolagus cuniculus) [37]. Nonetheless, hare species seem to remain as dead-end hosts, and rabbits appear to have a higher susceptibility to lethal disease than other leporids, notwithstanding the limited availability of comparative studies [38]. Phylogenetic analysis in countless previous publications from different labs globally clearly confirm the separation between "rabbit lagoviruses" and "hare lagoviruses", e.g., [2,7,13,30,33,39]. Furthermore, this partition is confirmed by genetic distances between sequences from these groups (Figure 1).

In another example, their GI.2 clade includes viruses recovered from leporids found dead in 1994–1996 (GenBank accession numbers KY765609–11 and KJ943791) [28], long before the emergence of GI.2 in 2010 [40]. The evolutionary history of these four viruses
has been robustly examined previously, demonstrating them to be recombinants with GI structural genes and non-structural genes from a distinct, unclassified genetic group [28]. By ignoring previous analyses and their own VP10 phylogeny, the authors accept that GI.2 was circulating in Portugal in the 1990s whilst causing disease outbreaks only about 15 years later, a much less parsimonious explanation than the existence of a GI.xP-GI.1b recombinant with lower epidemiological fitness, causing local infections and being outcompeted to extinction [28,41]. It is noteworthy that, at the time those viruses were collected, the strains circulating in the Iberian Peninsula were mostly GI.1b [42–44]. Hence, it is not surprising that recombination with GI.1b structural regions emerged at this time. Other errors that lead to misinterpretations include the wrong assignment of dates or countries (e.g., RHDV-SD was collected in 1989 rather than in 1993, which corresponds to the submission date in the databases, and RHDV pjG-RHDV-DD06 was collected in the UK rather than in Germany, where the laboratory that submitted the sequence was located). Importantly, the authors claim that “the FRG-USA strain (GenBank ID: NC_001543.1) genome isolated in 2000 is haplotype of the FRG-Germany strain (GenBank ID: M67473.1) isolated in 1991”, while, in fact, they are the same isolate, as clearly stated in Meyers et al. (2000) [45]: “the reference sequence was derived from” Meyers et al. (1991) [46]. As reported by Hicks and Duffy, a single misdated taxon can impact the inferences of evolutionary rates [47].

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Pairwise genetic distances of the dataset used by Shah et al. [1] (n = 240 sequences, 7337 nucleotides), calculated in MEGAX [48]. Distances were calculated as the number of base substitutions per site between sequences. Ambiguous positions were removed from each sequence pair. Darker color corresponds to the highest distances in the analysis, with a range of values between 0 (white) and 40% (dark orange).</figcaption></figure>

Finally, the criteria proposed by Shah et al. to assign viruses to groups and define variants in the phylogeny based on the full-length genome sequences is unintelligible, and these are critical points that may cause misinterpretation of lagovirus classification to the reader. For example, in Shah et al. [1], GI.1 sequences do not form a monophyletic group
has been robustly examined previously, demonstrating them to be recombinants with GI structural genes and non-structural genes from a distinct, unclassified genetic group [28]. By ignoring previous analyses and their own VP10 phylogeny, the authors accept that GI.2 was circulating in Portugal in the 1990s whilst causing disease outbreaks only about 15 years later, a much less parsimonious explanation than the existence of a GI.xP-GI.1b recombinant with lower epidemiological fitness, causing local infections and being outcompeted to extinction [28,41]. It is noteworthy that, at the time those viruses were collected, the strains circulating in the Iberian Peninsula were mostly GI.1b [42–44]. Hence, it is not surprising that recombination with GI.1b structural regions emerged at this time. Other errors that lead to misinterpretations include the wrong assignment of dates or countries (e.g., RHDV-SD was collected in 1989 rather than in 1993, which corresponds to the submission date in the databases, and RHDV pjG-RHDV-DD06 was collected in the UK rather than in Germany, where the laboratory that submitted the sequence was located). Importantly, the authors claim that “the FRG-USA strain (GenBank ID: NC_001543.1) genome isolated in 2000 is haplotype of the FRG-Germany strain (GenBank ID: M67473.1) isolated in 1991”, while, in fact, they are the same isolate, as clearly stated in Meyers et al. (2000) [45]: “the reference sequence was derived from” Meyers et al. (1991) [46]. As reported by Hicks and Duffy, a single misdated taxon can impact the inferences of evolutionary rates [47].

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Pairwise genetic distances of the dataset used by Shah et al. [1] (n = 240 sequences, 7337 nucleotides), calculated in MEGAX [48]. Distances were calculated as the number of base substitutions per site between sequences. Ambiguous positions were removed from each sequence pair. Darker color corresponds to the highest distances in the analysis, with a range of values between 0 (white) and 40% (dark orange).</figcaption></figure>

Finally, the criteria proposed by Shah et al. to assign viruses to groups and define variants in the phylogeny based on the full-length genome sequences is unintelligible, and these are critical points that may cause misinterpretation of lagovirus classification to the reader. For example, in Shah et al. [1], GI.1 sequences do not form a monophyletic group
in either the full-genome or the VP60 tree. Additionally, in the full genome tree, the GI.2 sequences are not monophyletic. Therefore, the classification system used by the authors is based on random/non-supported clusters that are not congruent with the epidemiology and known evolution of lagoviruses. Further, the “updated classification” system they provide lacks a substantial number of recently published near-complete genome sequences, as they exclude genome sequences that were published without the 5' and 3' end primer binding regions. This results in a less comprehensive and incomplete classification system even if the phylogenetic analysis had been conducted correctly.

In our view, at the time of publication, no evidence exists to support that the approach of using genome-wide alignments for maximum likelihood phylogenetic analysis will produce more reliable “hints for the genetic basis of [lagoviruses] emergence and re-emergence”, given the strong history of recombination in lagoviruses. Indeed, we suggest the opposite: that disregarding recombination when constructing phylogenies leads to missed and/or erroneous inferences by implying a shared evolutionary history of genomic regions that have, in fact, evolved separately. Tracing the evolutionary history of viruses is complex and requires accurate detection and consideration of recombination. In short, while full-genome information is highly useful for complete classification and evolutionary analysis, for genus-wide phylogenetic analysis, it is vital that both the structural and the non-structural encoding regions of lagovirus genomes are inferred independently. Ultimately, the work of Shah et al. does not provide “an update of the phylogenetic and phylogeographic information of lagoviruses”; rather, it proposes a new classification system designed on flawed analyses that will only generate confusion around an already complex classification system for these viruses.

**Funding:** Fundação para a Ciência e Tecnologia (FCT) supported the Assistant Researcher grant awarded to Joana Abrantes (CEECIND/00078/2017) and the Junior Researcher grant awarded to Ana M. Lopes (CEECIND/01388/2017).

**Conflicts of Interest:** Author Robyn N. Hall was employed by the company Ausvet Pty Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## References

1. Shah, P.T.; Bahoussi, A.N.; Yang, C.; Yao, G.; Dong, L.; Wu, C.; Xing, L. Genetic Characteristics and Phylogeographic Dynamics of Lagoviruses, 1988–2021. *Viruses* **2023**, *15*, 815. [CrossRef] [PubMed]

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. [CrossRef] [PubMed]

3. Mahar, J.; Jenckel, M.; Huang, N.; Smertina, E.; Holmes, E.; 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. [CrossRef] [PubMed]

4. Parra, G.I. Emergence of norovirus strains: A tale of two genes. *Virus Evol.* **2019**, *5*, vez048. [CrossRef] [PubMed]

5. Parra, G.I.; Squires, R.B.; Karangwa, C.K.; Johnson, J.A.; Lepore, C.J.; Sosnovtsev, S.V.; Green, K.Y. Static and Evolving Norovirus Genotypes: Implications for Epidemiology and Immunity. *PLoS Pathog.* **2017**, *13*, e1006136. [CrossRef] [PubMed]

6. Bull, R.A.; Tanaka, M.M.; White, P.A. Norovirus recombination. *J. Gen. Virol.* **2007**, *88*, 3347–3359. [CrossRef]

7. Szillat, K.P.; Höper, D.; Beer, M.; König, P. Full-genome sequencing of German rabbit haemorrhagic disease virus uncovers recombination between RHDV (GI.2) and EBHSV (GII.1). *Virus Evol.* **2020**, *6*, veaa080. [CrossRef] [PubMed]

8. Mahar, J.E.; Hall, R.N.; Shi, M.; Mourant, R.; Huang, N.; Strive, T.; Holmes, E.C. The discovery of three new hare lagoviruses reveals unexplored viral diversity in this genus. *Virus Evol.* **2019**, *5*, vez005. [CrossRef] [PubMed]

9. Silvério, D.; Lopes, A.M.; Melo-Ferreira, J.; Monterroso, P.; Serronha, A.; Alves, P.C.; Magalhaes, M.; Esteves, P.J.; Abrantes, J. Insights into the evolution of the new variant rabbit hemorrhagic disease virus (GI.2) and the identification of novel recombinant strains. *Transbound. Emerg. Dis.* **2017**, *65*, 983–992. [CrossRef]

10. Lopes, A.M.; Dalton, K.P.; Magalhães, M.J.; Parra, F.; Esteves, P.J.; Holmes, E.C.; Abrantes, J. Full genomic analysis of new variant rabbit hemorrhagic disease virus revealed multiple recombination events. *J. Gen. Virol.* **2015**, *96*, 1309–1319. [CrossRef]

11. Forrester, N.L.; Moss, S.R.; Turner, S.L.; Scirrmeier, H.; Gould, E.A. Recombination in rabbit haemorrhagic disease virus: Possible impact on evolution and epidemiology. *Virology* **2008**, *376*, 390–396. [CrossRef]
in either the full-genome or the VP60 tree. Additionally, in the full genome tree, the GI.2 sequences are not monophyletic. Therefore, the classification system used by the authors is based on random/non-supported clusters that are not congruent with the epidemiology and known evolution of lagoviruses. Further, the “updated classification” system they provide lacks a substantial number of recently published near-complete genome sequences, as they exclude genome sequences that were published without the 5' and 3' end primer binding regions. This results in a less comprehensive and incomplete classification system even if the phylogenetic analysis had been conducted correctly.

In our view, at the time of publication, no evidence exists to support that the approach of using genome-wide alignments for maximum likelihood phylogenetic analysis will produce more reliable “hints for the genetic basis of [lagoviruses] emergence and re-emergence”, given the strong history of recombination in lagoviruses. Indeed, we suggest the opposite: that disregarding recombination when constructing phylogenies leads to missed and/or erroneous inferences by implying a shared evolutionary history of genomic regions that have, in fact, evolved separately. Tracing the evolutionary history of viruses is complex and requires accurate detection and consideration of recombination. In short, while full-genome information is highly useful for complete classification and evolutionary analysis, for genus-wide phylogenetic analysis, it is vital that both the structural and the non-structural encoding regions of lagovirus genomes are inferred independently. Ultimately, the work of Shah et al. does not provide “an update of the phylogenetic and phylogeographic information of lagoviruses”; rather, it proposes a new classification system designed on flawed analyses that will only generate confusion around an already complex classification system for these viruses.

**Funding:** Fundação para a Ciência e Tecnologia (FCT) supported the Assistant Researcher grant awarded to Joana Abrantes (CEECIND/00078/2017) and the Junior Researcher grant awarded to Ana M. Lopes (CEECIND/01388/2017).

**Conflicts of Interest:** Author Robyn N. Hall was employed by the company Ausvet Pty Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## References

1. Shah, P.T.; Bahoussi, A.N.; Yang, C.; Yao, G.; Dong, L.; Wu, C.; Xing, L. Genetic Characteristics and Phylogeographic Dynamics of Lagoviruses, 1988-2021. *Viruses* **2023**, *15*, 815. [CrossRef] [PubMed]

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. [CrossRef] [PubMed]

3. Mahar, J.; Jenckel, M.; Huang, N.; Smertina, E.; Holmes, E.; 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. [CrossRef] [PubMed]

4. Parra, G.I. Emergence of norovirus strains: A tale of two genes. *Virus Evol.* **2019**, *5*, vez048. [CrossRef] [PubMed]

5. Parra, G.I.; Squires, R.B.; Karangwa, C.K.; Johnson, J.A.; Lepore, C.J.; Sosnovtsev, S.V.; Green, K.Y. Static and Evolving Norovirus Genotypes: Implications for Epidemiology and Immunity. *PLoS Pathog.* **2017**, *13*, e1006136. [CrossRef] [PubMed]

6. Bull, R.A.; Tanaka, M.M.; White, P.A. Norovirus recombination. *J. Gen. Virol.* **2007**, *88*, 3347-3359. [CrossRef]

7. Szillat, K.P.; Höper, D.; Beer, M.; König, P. Full-genome sequencing of German rabbit haemorrhagic disease virus uncovers recombination between RHDV (GI.2) and EBHSV (GII.1). *Virus Evol.* **2020**, *6*, veaa080. [CrossRef] [PubMed]

8. Mahar, J.E.; Hall, R.N.; Shi, M.; Mourant, R.; Huang, N.; Strive, T.; Holmes, E.C. The discovery of three new hare lagoviruses reveals unexplored viral diversity in this genus. *Virus Evol.* **2019**, *5*, vez005. [CrossRef] [PubMed]

9. Silvério, D.; Lopes, A.M.; Melo-Ferreira, J.; Monterroso, P.; Serronha, A.; Alves, P.C.; Magalhaes, M.; Esteves, P.J.; Abrantes, J. Insights into the evolution of the new variant rabbit hemorrhagic disease virus (GI.2) and the identification of novel recombinant strains. *Transbound. Emerg. Dis.* **2017**, *65*, 983-992. [CrossRef]

10. Lopes, A.M.; Dalton, K.P.; Magalhães, M.J.; Parra, F.; Esteves, P.J.; Holmes, E.C.; Abrantes, J. Full genomic analysis of new variant rabbit hemorrhagic disease virus revealed multiple recombination events. *J. Gen. Virol.* **2015**, *96*, 1309-1319. [CrossRef]

11. Forrester, N.L.; Moss, S.R.; Turner, S.L.; Scirrmeier, H.; Gould, E.A. Recombination in rabbit haemorrhagic disease virus: Possible impact on evolution and epidemiology. *Virology* **2008**, *376*, 390-396. [CrossRef]
12. Abrantes, J.; Lopes, A.M.; Lemaitre, E.; Ahola, H.; Banihashem, F.; Droillard, C.; Marchandeau, S.; Esteves, P.J.; Neimanis, A.; Le Gall-Reculé, G. Retrospective Analysis Shows That Most RHDV GI.1 Strains Circulating Since the Late 1990s in France and Sweden Were Recombinant GI.3P-GI.1d Strains. *Genes* **2020**, *11*, 910. [CrossRef] [PubMed]

13. Cavadini, P.; Trogu, T.; Velarde, R.; Lavazza, A.; Capucci, L. Recombination between non-structural and structural genes as a mechanism of selection in lagoviruses: The evolutionary dead-end of an RHDV2 isolated from European hare. *Virus Res.* **2024**, *339*, 199257. [CrossRef] [PubMed]

14. Abrantes, J.; Droillard, C.; Lopes, A.M.; Lemaitre, E.; Lucas, P.; Blanchard, Y.; Marchandeau, S.; Esteves, P.J.; Le Gall-Reculé, G. Recombination at the emergence of the pathogenic rabbit haemorrhagic disease virus *Lagovirus europaeus* /GI.2. *Sci. Rep.* **2020**, *10*, 14502. [CrossRef] [PubMed]

15. Li, Y.; Du, D.; Zhou, L.; Su, L.; You, C.; Zhang, H.; Yu, J.; Xiao, L.; Huang, J. First report of GI.1aP-GI.2 recombinants of rabbit hemorrhagic disease virus in domestic rabbits in China. *Front. Microbiol.* **2023**, *14*, 1188380. [CrossRef] [PubMed]

16. Mahar, J.E.; Read, A.J.; Gu, X.; Urakova, N.; Mourant, R.; Piper, M.; Haboury, S.; Holmes, E.C.; Strive, T.; Hall, R.N. Detection and circulation of a novel rabbit hemorrhagic disease virus in Australia. *Emerg. Infect. Dis.* **2018**, *24*, 22–31. [CrossRef] [PubMed]

17. Schierup, M.H.; Hein, J. Consequences of recombination on traditional phylogenetic analysis. *Genetics* **2000**, *156*, 879–891. [CrossRef] [PubMed]

18. Posada, D. How does recombination affect phylogeny estimation? *Trends Ecol. Evol.* **2000**, *15*, 489–490. [CrossRef] [PubMed]

19. Zaman, S.; Sledzieski, S.; Berger, B.; Wu, Y.-C.; Bansal, M.S. virDTL: Viral Recombination Analysis through Phylogenetic Reconciliation and Its Application to Sarbecoviruses and SARS-CoV-2. *J. Comput. Biol.* **2023**, *30*, 3–20. [CrossRef]

20. Posada, D.; Crandall, K.A. The effect of recombination on the accuracy of phylogeny estimation. *J. Mol. Evol.* **2002**, *54*, 396–402. [CrossRef]

21. Sallard, E.; Halloy, J.; Casane, D.; Decroly, E.; van Helden, J. Tracing the origins of SARS-COV-2 in coronavirus phylogenies: A review. *Environ. Chem. Lett.* **2021**, *19*, 769–785. [CrossRef] [PubMed]

22. Tohma, K.; Lepore, C.J.; Martinez, M.; Degiuseppe, J.I.; Khamrin, P.; Saito, M.; Mayta, H.; Nwaba, A.U.A.; Ford-Siltz, L.A.; Green, K.Y.; et al. Genome-wide analyses of human noroviruses provide insights on evolutionary dynamics and evidence of coexisting viral populations evolving under recombination constraints. *PLoS Pathog.* **2021**, *17*, e1009744. [CrossRef] [PubMed]

23. Kroneman, A.; Vega, E.; Vennema, H.; Vinjé, J.; White, P.A.; Hansman, G.; Green, K.; Martella, V.; Katayama, K.; Koopmans, M. Proposal for a unified norovirus nomenclature and genotyping. *Arch. Virol.* **2013**, *158*, 2059–2068. [CrossRef]

24. Chhabra, P.; de Graaf, M.; Parra, G.I.; Chan, M.C.-W.; Green, K.; Martella, V.; Wang, Q.; White, P.A.; Katayama, K.; Vennema, H.; et al. Updated classification of norovirus genogroups and genotypes. *J. Gen. Virol.* **2019**, *100*, 1393–1406. [CrossRef]

25. Stott, C.M.; Bobay, L.-M. Impact of homologous recombination on core genome phylogenies. *BMC Genom.* **2020**, *21*, 829. [CrossRef] [PubMed]

26. Holmes, E.C. *The Evolution and Emergence of RNA Viruses*; Oxford Series in Ecology and Evolution (OSEE); Oxford University Press: Oxford, UK, 2009.

27. Lopes, A.M.; Blanco-Aguiar, J.; Martín-Alonso, A.; Leitão, M.; Foronda, P.; Mendes, M.; Gonçalves, D.; Abrantes, J.; Esteves, P.J. Full genome sequences are key to disclose RHDV2 emergence in the Macaronesian islands. *Virus Genes* **2018**, *54*, 1–4. [CrossRef]

28. Lopes, A.M.; Silvério, D.; Magalhães, M.J.; Areal, H.; Alves, P.C.; Esteves, P.J.; Abrantes, J. Characterization of old RHDV strains by complete genome sequencing identifies a novel genetic group. *Sci. Rep.* **2017**, *7*, 13599. [CrossRef]

29. Abrantes, J.; Esteves, P.J.; van der Loo, W. Evidence for recombination in the major capsid gene VP60 of the rabbit haemorrhagic disease virus (RHDV). *Arch. Virol.* **2008**, *153*, 329–335. [CrossRef]

30. Droillard, C.; Lemaitre, E.; Chatel, M.; Quéméner, A.; Briand, F.-X.; Guitton, J.-S.; Marchandeau, S.; Le Pendu, J.; Eterradossi, N.; Le Gall-Reculé, G. Genetic diversity and evolution of Hare Calicivirus (HaCV), a recently identified lagovirus from *Lepus europaeus*. *Infect. Genet. Evol.* **2020**, *82*, 104310. [CrossRef]

31. Puggioni, G.; Cavadini, P.; Maestrale, C.; Scivoli, R.; Botti, G.; Ligios, C.; Le Gall-Reculé, G.; Lavazza, A.; Capucci, L. The new French 2010 *Rabbit Hemorrhagic Disease Virus* causes an RHD-like disease in the Sardinian Cape Hare (*Lepus capensis mediterraneus*). *Vet. Res.* **2013**, *44*, 96. [CrossRef]

32. Hall, R.N.; Peacock, D.E.; Kovaliski, J.; Mahar, J.E.; Mourant, R.; Piper, M.; Strive, T. Detection of RHDV2 in European brown hares (*Lepus europaeus*) in Australia. *Vet. Rec.* **2017**, *180*, 121. [CrossRef] [PubMed]

33. Neimanis, A.S.; Ahola, H.; Larsson Pettersson, U.; Lopes, A.M.; Abrantes, J.; Zohari, S.; Esteves, P.J.; Gavier-Widén, D. Overcoming species barriers: An outbreak of Lagovirus europaeus GI.2/RHDV2 in an isolated population of mountain hares (Lepus timidus). *BMC Vet. Res.* **2018**, *14*, 367. [CrossRef] [PubMed]

34. Asin, J.; Nyaoke, A.C.; Moore, J.D.; Gonzalez-Astudillo, V.; Clifford, D.L.; Lantz, E.L.; Mikolon, A.B.; Dodd, K.A.; Crossley, B.; Uzal, F.A. Outbreak of rabbit hemorrhagic disease virus 2 in the southwestern United States: First detections in southern California. *J. Vet. Diagn. Investig.* **2021**, *33*, 728–731. [CrossRef] [PubMed]

35. Lankton, J.; Knowles, S.; Keller, S.; Shearn-Bochsler, V.I.; Ip, H.S. Pathology of *Lagovirus europaeus* GI.2/RHDV2/b (Rabbit Hemorrhagic Disease Virus 2) in Native North American Lagomorphs. *J. Wildl. Dis.* **2021**, *57*, 694–700. [CrossRef] [PubMed]

36. Crowell, M.; LaHue, N.; Heath, E.; Shoemaker, K.; Matocq, M. Detection of *Rabbit hemorrhagic disease virus 2* in the Pygmy Rabbit (*Brachylagus idahoensis*) in Nevada, USA. *J. Wildl. Dis.* **2023**, *59*, 342–346. [CrossRef] [PubMed]
12. Abrantes, J.; Lopes, A.M.; Lemaitre, E.; Ahola, H.; Banihashem, F.; Droillard, C.; Marchandeau, S.; Esteves, P.J.; Neimanis, A.; Le Gall-Reculé, G. Retrospective Analysis Shows That Most RHDV GI.1 Strains Circulating Since the Late 1990s in France and Sweden Were Recombinant GI.3P-GI.1d Strains. *Genes* **2020**, *11*, 910. [CrossRef] [PubMed]

13. Cavadini, P.; Trogu, T.; Velarde, R.; Lavazza, A.; Capucci, L. Recombination between non-structural and structural genes as a mechanism of selection in lagoviruses: The evolutionary dead-end of an RHDV2 isolated from European hare. *Virus Res.* **2024**, *339*, 199257. [CrossRef] [PubMed]

14. Abrantes, J.; Droillard, C.; Lopes, A.M.; Lemaitre, E.; Lucas, P.; Blanchard, Y.; Marchandeau, S.; Esteves, P.J.; Le Gall-Reculé, G. Recombination at the emergence of the pathogenic rabbit haemorrhagic disease virus *Lagovirus europaeus* /GI.2. *Sci. Rep.* **2020**, *10*, 14502. [CrossRef] [PubMed]

15. Li, Y.; Du, D.; Zhou, L.; Su, L.; You, C.; Zhang, H.; Yu, J.; Xiao, L.; Huang, J. First report of GI.1aP-GI.2 recombinants of rabbit hemorrhagic disease virus in domestic rabbits in China. *Front. Microbiol.* **2023**, *14*, 1188380. [CrossRef] [PubMed]

16. Mahar, J.E.; Read, A.J.; Gu, X.; Urakova, N.; Mourant, R.; Piper, M.; Haboury, S.; Holmes, E.C.; Strive, T.; Hall, R.N. Detection and circulation of a novel rabbit hemorrhagic disease virus in Australia. *Emerg. Infect. Dis.* **2018**, *24*, 22–31. [CrossRef] [PubMed]

17. Schierup, M.H.; Hein, J. Consequences of recombination on traditional phylogenetic analysis. *Genetics* **2000**, *156*, 879–891. [CrossRef] [PubMed]

18. Posada, D. How does recombination affect phylogeny estimation? *Trends Ecol. Evol.* **2000**, *15*, 489–490. [CrossRef] [PubMed]

19. Zaman, S.; Sledzieski, S.; Berger, B.; Wu, Y.-C.; Bansal, M.S. virDTL: Viral Recombination Analysis through Phylogenetic Reconciliation and Its Application to Sarbecoviruses and SARS-CoV-2. *J. Comput. Biol.* **2023**, *30*, 3–20. [CrossRef]

20. Posada, D.; Crandall, K.A. The effect of recombination on the accuracy of phylogeny estimation. *J. Mol. Evol.* **2002**, *54*, 396–402. [CrossRef]

21. Sallard, E.; Halloy, J.; Casane, D.; Decroly, E.; van Helden, J. Tracing the origins of SARS-COV-2 in coronavirus phylogenies: A review. *Environ. Chem. Lett.* **2021**, *19*, 769–785. [CrossRef] [PubMed]

22. Tohma, K.; Lepore, C.J.; Martinez, M.; Degiuseppe, J.I.; Khamrin, P.; Saito, M.; Mayta, H.; Nwaba, A.U.A.; Ford-Siltz, L.A.; Green, K.Y.; et al. Genome-wide analyses of human noroviruses provide insights on evolutionary dynamics and evidence of coexisting viral populations evolving under recombination constraints. *PLoS Pathog.* **2021**, *17*, e1009744. [CrossRef] [PubMed]

23. Kroneman, A.; Vega, E.; Vennema, H.; Vinjé, J.; White, P.A.; Hansman, G.; Green, K.; Martella, V.; Katayama, K.; Koopmans, M. Proposal for a unified norovirus nomenclature and genotyping. *Arch. Virol.* **2013**, *158*, 2059–2068. [CrossRef]

24. Chhabra, P.; de Graaf, M.; Parra, G.I.; Chan, M.C.-W.; Green, K.; Martella, V.; Wang, Q.; White, P.A.; Katayama, K.; Vennema, H.; et al. Updated classification of norovirus genogroups and genotypes. *J. Gen. Virol.* **2019**, *100*, 1393–1406. [CrossRef]

25. Stott, C.M.; Bobay, L.-M. Impact of homologous recombination on core genome phylogenies. *BMC Genom.* **2020**, *21*, 829. [CrossRef] [PubMed]

26. Holmes, E.C. *The Evolution and Emergence of RNA Viruses*; Oxford Series in Ecology and Evolution (OSEE); Oxford University Press: Oxford, UK, 2009.

27. Lopes, A.M.; Blanco-Aguiar, J.; Martín-Alonso, A.; Leitão, M.; Foronda, P.; Mendes, M.; Gonçalves, D.; Abrantes, J.; Esteves, P.J. Full genome sequences are key to disclose RHDV2 emergence in the Macaronesian islands. *Virus Genes* **2018**, *54*, 1–4. [CrossRef]

28. Lopes, A.M.; Silvério, D.; Magalhães, M.J.; Areal, H.; Alves, P.C.; Esteves, P.J.; Abrantes, J. Characterization of old RHDV strains by complete genome sequencing identifies a novel genetic group. *Sci. Rep.* **2017**, *7*, 13599. [CrossRef]

29. Abrantes, J.; Esteves, P.J.; van der Loo, W. Evidence for recombination in the major capsid gene VP60 of the rabbit haemorrhagic disease virus (RHDV). *Arch. Virol.* **2008**, *153*, 329–335. [CrossRef]

30. Droillard, C.; Lemaitre, E.; Chatel, M.; Quéméner, A.; Briand, F.-X.; Guitton, J.-S.; Marchandeau, S.; Le Pendu, J.; Eterradossi, N.; Le Gall-Reculé, G. Genetic diversity and evolution of Hare Calicivirus (HaCV), a recently identified lagovirus from *Lepus europaeus*. *Infect. Genet. Evol.* **2020**, *82*, 104310. [CrossRef]

31. Puggioni, G.; Cavadini, P.; Maestrale, C.; Scivoli, R.; Botti, G.; Ligios, C.; Le Gall-Reculé, G.; Lavazza, A.; Capucci, L. The new French 2010 *Rabbit Hemorrhagic Disease Virus* causes an RHD-like disease in the Sardinian Cape Hare (*Lepus capensis mediterraneus*). *Vet. Res.* **2013**, *44*, 96. [CrossRef]

32. Hall, R.N.; Peacock, D.E.; Kovaliski, J.; Mahar, J.E.; Mourant, R.; Piper, M.; Strive, T. Detection of RHDV2 in European brown hares (*Lepus europaeus*) in Australia. *Vet. Rec.* **2017**, *180*, 121. [CrossRef] [PubMed]

33. Neimanis, A.S.; Ahola, H.; Larsson Pettersson, U.; Lopes, A.M.; Abrantes, J.; Zohari, S.; Esteves, P.J.; Gavier-Widén, D. Overcoming species barriers: An outbreak of Lagovirus europaeus GI.2/RHDV2 in an isolated population of mountain hares (Lepus timidus). *BMC Vet. Res.* **2018**, *14*, 367. [CrossRef] [PubMed]

34. Asin, J.; Nyaoke, A.C.; Moore, J.D.; Gonzalez-Astudillo, V.; Clifford, D.L.; Lantz, E.L.; Mikolon, A.B.; Dodd, K.A.; Crossley, B.; Uzal, F.A. Outbreak of rabbit hemorrhagic disease virus 2 in the southwestern United States: First detections in southern California. *J. Vet. Diagn. Investig.* **2021**, *33*, 728–731. [CrossRef] [PubMed]

35. Lankton, J.; Knowles, S.; Keller, S.; Shearn-Bochsler, V.I.; Ip, H.S. Pathology of *Lagovirus europaeus* GI.2/RHDV2/b (Rabbit Hemorrhagic Disease Virus 2) in Native North American Lagomorphs. *J. Wildl. Dis.* **2021**, *57*, 694–700. [CrossRef] [PubMed]

36. Crowell, M.; LaHue, N.; Heath, E.; Shoemaker, K.; Matocq, M. Detection of *Rabbit hemorrhagic disease virus 2* in the Pygmy Rabbit (*Brachylagus idahoensis*) in Nevada, USA. *J. Wildl. Dis.* **2023**, *59*, 342–346. [CrossRef] [PubMed]
37. Le Gall-Reculé, G.; Zwingelstein, F.; Boucher, S.; Le Normand, B.; Plassiart, G.; Portejoie, Y.; Decors, A.; Bertagnoli, S.; Guérin, J.L.; Marchandeau, S. Detection of a new variant of rabbit haemorrhagic disease virus in France. *Vet. Rec.* **2011**, *168*, 137–138.
[CrossRef] [PubMed]

38. Mohamed, F.; Gidlewski, T.; Berninger, M.L.; Petrowski, H.M.; Bracht, A.J.; de Rueda, C.B.; Barrette, R.W.; Grady, M.; O'Hearn, E.S.; Lewis, C.E.; et al. Comparative susceptibility of eastern cottontails and New Zealand white rabbits to classical rabbit haemorrhagic disease virus (RHDV) and RHDV2. *Transbound. Emerg. Dis.* **2022**, *69*, e968–e978. [CrossRef] [PubMed]

39. Koh, E.Y.; Ong, J.; Wang, Y.; Toh, X.; Fernandez, C.J.; Huangfu, T.; Hall, R.N.; Toh, S.; Lim, K.; Sng, W.; et al. Rabbit haemorrhagic disease virus 2 from Singapore 2020 outbreak revealed an Australian recombinant variant. *Virus Evol.* **2023**, *9*, vead029. [CrossRef] [PubMed]

40. Le Gall-Reculé, 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. [CrossRef]

41. Lopes, A.M.; Marques, S.; Silva, E.; Magalhaes, M.J.; Pinheiro, A.; Alves, P.C.; Le Pendu, J.; Esteves, P.J.; Thompson, G.; Abrantes, J. Detection of RHDV strains in the Iberian hare (*Lepus granatensis*): Earliest evidence of rabbit lagovirus cross-species infection. *Vet. Res.* **2014**, *45*, 94.

42. Müller, A.; Freitas, J.; Silva, E.; Le Gall-Reculé, G.; Zwingelstein, F.; Abrantes, J.; Esteves, P.J.; Alves, P.C.; van der Loo, W.; Kolodziejek, J.; et al. Evolution of rabbit haemorrhagic disease virus (RHDV) in the European rabbit (*Oryctolagus cuniculus*) from the Iberian Peninsula. *Vet. Microbiol.* **2009**, *135*, 368–373. [CrossRef] [PubMed]

43. Lopes, A.M.; Magalhães, M.J.; Alves, P.C.; Esteves, P.J.; Abrantes, J. An update on the rabbit hemorrhagic disease virus (RHDV) strains circulating in Portugal in the 1990s: Earliest detection of G3-G5 and G6. *Arch. Virol.* **2017**, *162*, 2061–2065. [CrossRef] [PubMed]

44. Alda, F.; Gaitero, T.; Suárez, M.; Merchán, T.; Rocha, G.; Doadrio, I. Evolutionary history and molecular epidemiology of rabbit haemorrhagic disease virus in the Iberian Peninsula and Western Europe. *BMC Evol. Biol.* **2010**, *10*, 347. [CrossRef] [PubMed]

45. Meyers, G.; Wirblich, C.; Thiel, H.-J.; Thumfart, J.O. Rabbit hemorrhagic disease virus: Genome organization and polyprotein processing of a calicivirus studied after transient expression of cDNA constructs. *Virology* **2000**, *276*, 349–363. [CrossRef] [PubMed]

46. Meyers, G.; Wirblich, C.; Thiel, H.J. Rabbit hemorrhagic disease virus—Molecular cloning and nucleotide sequencing of a calicivirus genome. *Virology* **1991**, *184*, 664–676. [CrossRef] [PubMed]

47. Hicks, A.L.; Duffy, S. One misdated sequence of rabbit hemorrhagic disease virus prevents accurate estimation of its nucleotide substitution rate. *BMC Evol. Biol.* **2012**, *12*, 74. [CrossRef]

48. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. *Mol. Biol. Evol.* **2018**, *35*, 1547–1549. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
37. Le Gall-Reculé, G.; Zwingelstein, F.; Boucher, S.; Le Normand, B.; Plassiart, G.; Portejoie, Y.; Decors, A.; Bertagnoli, S.; Guérin, J.L.; Marchandeau, S. Detection of a new variant of rabbit haemorrhagic disease virus in France. *Vet. Rec.* **2011**, *168*, 137–138.
[CrossRef] [PubMed]

38. Mohamed, F.; Gidlewski, T.; Berninger, M.L.; Petrowski, H.M.; Bracht, A.J.; de Rueda, C.B.; Barrette, R.W.; Grady, M.; O'Hearn, E.S.; Lewis, C.E.; et al. Comparative susceptibility of eastern cottontails and New Zealand white rabbits to classical rabbit haemorrhagic disease virus (RHDV) and RHDV2. *Transbound. Emerg. Dis.* **2022**, *69*, e968–e978. [CrossRef] [PubMed]

39. Koh, E.Y.; Ong, J.; Wang, Y.; Toh, X.; Fernandez, C.J.; Huangfu, T.; Hall, R.N.; Toh, S.; Lim, K.; Sng, W.; et al. Rabbit haemorrhagic disease virus 2 from Singapore 2020 outbreak revealed an Australian recombinant variant. *Virus Evol.* **2023**, *9*, vead029. [CrossRef] [PubMed]

40. Le Gall-Reculé, 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. [CrossRef]

41. Lopes, A.M.; Marques, S.; Silva, E.; Magalhaes, M.J.; Pinheiro, A.; Alves, P.C.; Le Pendu, J.; Esteves, P.J.; Thompson, G.; Abrantes, J. Detection of RHDV strains in the Iberian hare (*Lepus granatensis*): Earliest evidence of rabbit lagovirus cross-species infection. *Vet. Res.* **2014**, *45*, 94.

42. Müller, A.; Freitas, J.; Silva, E.; Le Gall-Reculé, G.; Zwingelstein, F.; Abrantes, J.; Esteves, P.J.; Alves, P.C.; van der Loo, W.; Kolodziejek, J.; et al. Evolution of rabbit haemorrhagic disease virus (RHDV) in the European rabbit (*Oryctolagus cuniculus*) from the Iberian Peninsula. *Vet. Microbiol.* **2009**, *135*, 368–373. [CrossRef] [PubMed]

43. Lopes, A.M.; Magalhães, M.J.; Alves, P.C.; Esteves, P.J.; Abrantes, J. An update on the rabbit hemorrhagic disease virus (RHDV) strains circulating in Portugal in the 1990s: Earliest detection of G3-G5 and G6. *Arch. Virol.* **2017**, *162*, 2061–2065. [CrossRef] [PubMed]

44. Alda, F.; Gaitero, T.; Suárez, M.; Merchán, T.; Rocha, G.; Doadrio, I. Evolutionary history and molecular epidemiology of rabbit haemorrhagic disease virus in the Iberian Peninsula and Western Europe. *BMC Evol. Biol.* **2010**, *10*, 347. [CrossRef] [PubMed]

45. Meyers, G.; Wirblich, C.; Thiel, H.-J.; Thumfart, J.O. Rabbit hemorrhagic disease virus: Genome organization and polyprotein processing of a calicivirus studied after transient expression of cDNA constructs. *Virology* **2000**, *276*, 349–363. [CrossRef] [PubMed]

46. Meyers, G.; Wirblich, C.; Thiel, H.J. Rabbit hemorrhagic disease virus—Molecular cloning and nucleotide sequencing of a calicivirus genome. *Virology* **1991**, *184*, 664–676. [CrossRef] [PubMed]

47. Hicks, A.L.; Duffy, S. One misdated sequence of rabbit hemorrhagic disease virus prevents accurate estimation of its nucleotide substitution rate. *BMC Evol. Biol.* **2012**, *12*, 74. [CrossRef]

48. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. *Mol. Biol. Evol.* **2018**, *35*, 1547–1549. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.